Issue 
A&A
Volume 599, March 2017



Article Number  A51  
Number of page(s)  15  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201629164  
Published online  28 February 2017 
Planck intermediate results
L. Evidence of spatial variation of the polarized thermal dust spectral energy distribution and implications for CMB Bmode analysis
^{1} APC, AstroParticule et Cosmologie, Université Paris Diderot, CNRS/IN2P3, CEA/lrfu, Observatoire de Paris, Sorbonne Paris Cité, 10 rue Alice Domon et Léonie Duquet, 75205 Paris Cedex 13, France
^{2} Aalto University Metsähovi Radio Observatory and Dept of Radio Science and Engineering, PO Box 13000, 00076 Aalto, Finland
^{3} African Institute for Mathematical Sciences, 6–8 Melrose Road, Muizenberg, 7945 Cape Town, South Africa
^{4} Agenzia Spaziale Italiana Science Data Center, via del Politecnico snc, 00133 Roma, Italy
^{5} Aix Marseille Université, CNRS, LAM (Laboratoire d’Astrophysique de Marseille) UMR 7326, 13388 Marseille, France
^{6} Astrophysics Group, Cavendish Laboratory, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, UK
^{7} Astrophysics & Cosmology Research Unit, School of Mathematics, Statistics & Computer Science, University of KwaZuluNatal, Westville Campus, Private Bag X54001, 4000 Durban, South Africa
^{8} CITA, University of Toronto, 60 St. George St., Toronto, ON M5S 3H8, Canada
^{9} CNRS, IRAP, 9 Av. colonel Roche, BP 44346, 31028 Toulouse Cedex 4, France
^{10} California Institute of Technology, Pasadena, 1200 California, USA
^{11} Computational Cosmology Center, Lawrence Berkeley National Laboratory, Berkeley, 94720 California, USA
^{12} DTU Space, National Space Institute, Technical University of Denmark, Elektrovej 327, 2800 Kgs. Lyngby, Denmark
^{13} Département de Physique Théorique, Université de Genève, 24, Quai E. Ansermet, 1211 Genève 4, Switzerland
^{14} Departamento de Astrofísica, Universidad de La Laguna (ULL), 38206 La Laguna, Tenerife, Spain
^{15} Departamento de Física, Universidad de Oviedo, Avda. Calvo Sotelo s/n, 33003 Oviedo, Spain
^{16} Department of Astrophysics/IMAPP, Radboud University Nijmegen, PO Box 9010, 6500 GL Nijmegen, The Netherlands
^{17} Department of Physics & Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, British Columbia, Canada
^{18} Department of Physics and Astronomy, Dana and David Dornsife College of Letter, Arts and Sciences, University of Southern California, Los Angeles, CA 90089, USA
^{19} Department of Physics and Astronomy, University College London, London WC1E 6BT, UK
^{20} Department of Physics, Gustaf Hällströmin katu 2a, University of Helsinki, 00100 Helsinki, Finland
^{21} Department of Physics, Princeton University, 08544 Princeton, New Jersey, USA
^{22} Department of Physics, University of California, 93106 Santa Barbara, California, USA
^{23} Department of Physics, University of Illinois at UrbanaChampaign, 1110 West Green Street, Urbana, Illinois, USA
^{24} Dipartimento di Fisica e Astronomia G. Galilei, Università degli Studi di Padova, via Marzolo 8, 35131 Padova, Italy
^{25} Dipartimento di Fisica e Astronomia, Alma Mater Studiorum, Università degli Studi di Bologna, Viale Berti Pichat 6/2, 40127 Bologna, Italy
^{26} Dipartimento di Fisica e Scienze della Terra, Università di Ferrara, via Saragat 1, 44122 Ferrara, Italy
^{27} Dipartimento di Fisica, Università La Sapienza, P.le A. Moro 2, Roma, Italy
^{28} Dipartimento di Fisica, Università degli Studi di Milano, via Celoria, 16, Milano, Italy
^{29} Dipartimento di Fisica, Università degli Studi di Trieste, via A. Valerio 2, Trieste, Italy
^{30} Dipartimento di Matematica, Università di Roma Tor Vergata, via della Ricerca Scientifica 1, Roma, Italy
^{31} Discovery Center, Niels Bohr Institute, Copenhagen University, Blegdamsvej 17, Copenhagen, Denmark
^{32} European Space Agency, ESAC, Planck Science Office, Camino bajo del Castillo, s/n, Urbanización Villafranca del Castillo, 28692 Villanueva de la Cañada, Madrid, Spain
^{33} European Space Agency, ESTEC, Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands
^{34} Gran Sasso Science Institute, INFN, viale F. Crispi 7, 67100 L’ Aquila, Italy
^{35} HGSFP and University of Heidelberg, Theoretical PhysicsDepartment, Philosophenweg 16, 69120 Heidelberg, Germany
^{36} Haverford College Astronomy Department, 370 Lancaster Avenue, Haverford, Pennsylvania, USA
^{37} Helsinki Institute of Physics, Gustaf Hällströmin katu 2, University of Helsinki, 00100 Helsinki, Finland
^{38} INAF–Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio 5, Padova, Italy
^{39} INAF–Osservatorio Astronomico di Roma, via di Frascati 33, Monte Porzio Catone, Italy
^{40} INAF–Osservatorio Astronomico di Trieste, via G.B. Tiepolo 11, Trieste, Italy
^{41} INAF/IASF Bologna, via Gobetti 101, 40127 Bologna, Italy
^{42} INAF/IASF Milano, via E. Bassini 15, 20133 Milano, Italy
^{43} INFN–CNAF, viale Berti Pichat 6/2, 40127 Bologna, Italy
^{44} INFN, Sezione di Bologna, viale Berti Pichat 6/2, 40127 Bologna, Italy
^{45} INFN, Sezione di Ferrara, via Saragat 1, 44122 Ferrara, Italy
^{46} INFN, Sezione di Roma 1, Università di Roma Sapienza, Piazzale Aldo Moro 2, 00185 Roma, Italy
^{47} INFN, Sezione di Roma 2, Università di Roma Tor Vergata, via della Ricerca Scientifica, 1, 00185 Roma, Italy
^{48} INFN/National Institute for Nuclear Physics, via Valerio 2, 34127 Trieste, Italy
^{49} IUCAA, Post Bag 4, Ganeshkhind, Pune University Campus, 411 007 Pune, India
^{50} Imperial College London, Astrophysics group, Blackett Laboratory, Prince Consort Road, London, SW7 2AZ, UK
^{51} Institut d’Astrophysique Spatiale, CNRS, Univ. ParisSud, Université ParisSaclay, Bât. 121, 91405 Orsay Cedex, France
^{52} Institut d’Astrophysique de Paris, CNRS (UMR7095), 98bis Boulevard Arago, 75014 Paris, France
^{53} Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK
^{54} Institute of Theoretical Astrophysics, University of Oslo, Blindern, 0313 Oslo, Norway
^{55} Instituto de Astrofísica de Canarias, C/Vía Láctea s/n, La Laguna, 38205 Tenerife, Spain
^{56} Instituto de Física de Cantabria (CSICUniversidad de Cantabria), Avda. de los Castros s/n, 39005 Santander, Spain
^{57} Istituto Nazionale di Fisica Nucleare, Sezione di Padova, via Marzolo 8, 35131 Padova, Italy
^{58} Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, California, USA
^{59} Jodrell Bank Centre for Astrophysics, Alan Turing Building, School of Physics and Astronomy, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
^{60} Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA
^{61} Kavli Institute for Cosmology Cambridge, Madingley Road, Cambridge, CB3 0HA, UK
^{62} LAL, Université ParisSud, CNRS/IN2P3, 91400 Orsay, France
^{63} LERMA, CNRS, Observatoire de Paris, 61 Avenue de l’Observatoire, 75014 Paris, France
^{64} Laboratoire AIM, IRFU/Service d’Astrophysique – CEA/DSM – CNRS – Université Paris Diderot, Bât. 709, CEASaclay, 91191 GifsurYvette Cedex, France
^{65} Laboratoire Traitement et Communication de l’Information, CNRS (UMR 5141) and Télécom ParisTech, 46 rue Barrault, 75634 Paris Cedex 13, France
^{66} Laboratoire de Physique Subatomique et Cosmologie, Université GrenobleAlpes, CNRS/IN2P3, 53 rue des Martyrs, 38026 Grenoble Cedex, France
^{67} Laboratoire de Physique Théorique, Université ParisSud 11 & CNRS, Bâtiment 210, 91405 Orsay, France
^{68} Lawrence Berkeley National Laboratory, Berkeley, California, USA
^{69} MaxPlanckInstitut für Astrophysik, KarlSchwarzschildStr. 1, 85741 Garching, Germany
^{70} Mullard Space Science Laboratory, University College London, Surrey RH5 6NT, UK
^{71} Nicolaus Copernicus Astronomical Center, Bartycka 18, 00716 Warsaw, Poland
^{72} Niels Bohr Institute, Copenhagen University, Blegdamsvej 17, 2100 Copenhagen, Denmark
^{73} Nordita (Nordic Institute for Theoretical Physics), Roslagstullsbacken 23, 106 91 Stockholm, Sweden
^{74} SISSA, Astrophysics Sector, via Bonomea 265, 34136, Trieste, Italy
^{75} School of Physics and Astronomy, University of Nottingham, Nottingham NG7 2RD, UK
^{76} Simon Fraser University, Department of Physics, 8888 University Drive, Burnaby BC, Canada
^{77} Sorbonne UniversitéUPMC, UMR7095, Institut d’Astrophysique de Paris, 98bis Boulevard Arago, 75014 Paris, France
^{78} Space Sciences Laboratory, University of California, 94720 Berkeley, California, USA
^{79} SubDepartment of Astrophysics, University of Oxford, Keble Road, Oxford OX1 3RH, UK
^{80} The Oskar Klein Centre for Cosmoparticle Physics, Department of Physics, Stockholm University, AlbaNova, 106 91 Stockholm, Sweden
^{81} UPMC Univ. Paris 06, UMR7095, 98bis Boulevard Arago, 75014 Paris, France
^{82} Université de Toulouse, UPSOMP, IRAP, 31028 Toulouse Cedex 4, France
^{83} University of Granada, Departamento de Física Teórica y del Cosmos, Facultad de Ciencias, 18010 Granada, Spain
^{84} Warsaw University Observatory, Aleje Ujazdowskie 4, 00478 Warszawa, Poland
Received: 22 June 2016
Accepted: 16 November 2016
The characterization of the Galactic foregrounds has been shown to be the main obstacle in thechallenging quest to detect primordial Bmodes in the polarized microwave sky. We make use of the PlanckHFI 2015 data release at high frequencies to place new constraints on the properties of the polarized thermal dust emission at high Galactic latitudes. Here, we specifically study the spatial variability of the dust polarized spectral energy distribution (SED), and its potential impact on the determination of the tensortoscalar ratio, r. We use the correlation ratio of the C^{BB}_{ℓ} angular power spectra between the 217 and 353 GHz channels as a tracer of these potential variations, computed on different high Galactic latitude regions, ranging from 80% to 20% of the sky. The new insight from Planck data is a departure of the correlation ratio from unity that cannot be attributed to a spurious decorrelation due to the cosmic microwave background, instrumental noise, or instrumental systematics. The effect is marginally detected on each region, but the statistical combination of all the regions gives more than 99% confidence for this variation in polarized dust properties. In addition, we show that the decorrelation increases when there is a decrease in the mean column density of the region of the sky being considered, and we propose a simple powerlaw empirical model for this dependence, which matches what is seen in the Planck data. We explore the effect that this measured decorrelation has on simulations of the BICEP2Keck Array/Planck analysis and show that the 2015 constraints from these data still allow a decorrelation between the dust at 150 and 353 GHz that is compatible with our measured value. Finally, using simplified models, we show that either spatial variation of the dust SED or of the dust polarization angle are able to produce decorrelations between 217 and 353 GHz data similar to the values we observe in the data.
Key words: cosmic background radiation / cosmology: observations / submillimeter: ISM / dust, extinction
© ESO, 2017
1. Introduction
The combined BICEP2/Keck Array and Planck^{1} analysis (BICEP2/Keck Array and Planck Collaborations 2015, hereafter BKP15), has confirmed that Galactic foregrounds likely constitute the dominant component of the Bmode polarization signal on a clean patch of the sky and at large scales (ℓ< 200). Characterization of these foregrounds is today the main limitation in the quest for the gravitational wave signature in the Bmode cosmic microwave background (CMB) power spectrum, i.e., for measuring the tensortoscalar ratio, r. This endeavour has been led by the analysis of the Galactic dust foregrounds carried out using Planck data (Planck Collaboration Int. XXX 2016, hereafter PIPXXX), and more recently of the synchrotron foregrounds using Planck and WMAP data (Planck Collaboration XXV 2016; Choi & Page 2015; Krachmalnicoff et al. 2016). It appears that Galactic thermal dust currently represents the major contaminant at high latitude in the spectral bands mainly adopted to search for the CMB signal by groundbased and balloonborne experiments, i.e., between 100 and 220 GHz (Planck Collaboration Int. XXII 2015; Planck Collaboration X 2016).
The efficiency of Galactic dust cleaning for the cosmological Bmode analyses is based on two main factors: the accuracy of the Galactic polarized dust template; and the way it is extrapolated from the submillimetre to the millimetre bands. Recent forecasts of cosmological Bmodes detection for groundbased, balloonborne, and satellite experiments (Creminelli et al. 2015; Errard et al. 2016) appear to allow for very good sensitivity (down to r ≃ 2 × 10^{3}) with only a few observational bands, when assuming a simple modelling of the foreground emission; however, other studies have shown that even in a global component separation framework, accurate modelling of the polarized dust spectral energy distribution (SED) is needed in order to reach the required very low levels of contamination of the cosmological Bmodes by Galactic dust residuals (ArmitageCaplan et al. 2012; Remazeilles et al. 2016). Incorrect modelling of the dust SED would lead to biased estimates of the r parameter.
Several investigations to quantify the spatial variability of the polarized dust SED were initiated using the Planck data. A first estimate of the polarized dust spectral index was discussed in Planck Collaboration Int. XXII (2015) for regions at intermediate Galactic latitudes. Over 39% of the sky the averaged value was shown to be slightly higher than the dust spectral index in intensity, with a value of , compared with . More interestingly, an upper limit of its spatial dispersion was obtained by computing the standard deviation of the mean polarized dust spectral index estimated on 352 patches of 10° diameter, yielding a dispersion of 0.17. However, this estimate has been shown to be dominated by the expected Planck noise, and does not allow us to build a reliable model of these spatial variations.
A second early approach, in PIPXXX, investigated the correlation ratio between the 217 and 353 GHz Planck bands, which is a statistical measurement of the dust SED spatial variation, as we will discuss thoroughly in the following sections. This ratio was computed on large fractions of the sky and led to an upper limit of 7% decorrelation between 217 and 353 GHz. An initial estimate of the impact of a possible dust polarization decorrelation between frequencies, due to a spatial variation of its SED, was performed in the BKP15 analysis. A loss of 10% in the 150 GHz × 353 GHz crossspectrum in the joint analysis, due to the decorrelation, was estimated to produce a positive bias of 0.018 on the determination of r.
In this new study, we present an improved analysis of the correlation ratio between the 353 GHz and 217 GHz Planck bands as a tracer of the spatial variations of the polarized dust SED. We derive new constraints on these variations and look at the impact they would have on the determination of the cosmological Bmode signal.
This paper is organized as follows. We present the data used in this work in Sect. 2. The analysis of the dust polarization correlation ratio is described in Sect. 3. The impact of the polarized dust decorrelation on the determination of the tensortoscalar ratio r is illustrated in Sect. 4. We discuss the possible origin of the spatial variations of the polarized dust SED in Sect. 5, before concluding in Sect. 6.
2. Data and region selection
2.1. Planck data
In this paper, we use the publicly available Planck High Frequency Instrument (HFI) data at 217 GHz and 353 GHz (Planck Collaboration I 2016). The signaltonoise ratio of dust polarization in the PlanckHFI maps at lower frequencies does not allow us to derive significant results in the framework adopted for this study and consequently the other channels are not used. These data consist of a set of maps of the Stokes Q and U parameters at each frequency, projected onto the HEALPix pixelization scheme (Górski et al. 2005). The maps and their properties are described in detail in Planck Collaboration VIII (2016).
We also use subsets of these data in order to exploit the statistical independence of the noise between them. As described in Planck Collaboration VIII (2016), the data were split into either the time or the detector domains. In this paper we use two particular data splits.

The “detectorset” maps (hereafter DS, also sometimes called “DetSets”). PlanckHFI measures the sky polarization thanks to polarization sensitive bolometers (PSBs), which are each sensitive to one direction of polarization. PSBs are assembled in pairs, the angle between two PSBs in a pair being 90° and the angle between two pairs being 45°, allowing for the reconstruction of Stokes I, Q, and U. Four such pairs are available at 217 and 353 GHz. These were split into two subsets of two pairs (four bolometers) to produce two noiseindependent sets of Q and U maps called “detectorsets”, for each frequency band.

The “halfmission” maps (hereafter HM). PlanckHFI completed five independent fullsky surveys over 30 months. Surveys 1 and 2 constitute the Planck “nominal mission” and this was repeated a second time during Surveys 3 and 4. A fifth survey was performed with a different scanning strategy, but was not included in the released halfmission maps. Thus, PlanckHFI data can be split into two noiseindependent sets of Q and U maps, labelled “HM1” and “HM2”, for each frequency band.
2.2. Region selection
In order to focus on the most diffuse areas of the Galactic dust emission, we have performed our analysis on various fractions of the sky using the set of science regions introduced in PIPXXX. These regions have been constructed to reject regions with CO line brightness I_{CO} ≥ 0.4 K km s^{1}. The complement to this mask by itself defines a preliminary region that retains a sky fraction f_{sky} = 0.8. In combination with thresholding, based on the Planck 857 GHz intensity map, six further preliminary regions are defined over f_{sky} = 0.2 to 0.7 in steps of 0.1. After point source masking and apodization, this procedure leads to the large retained (LR) regions LR16, LR24, LR33, LR42, LR53, LR63, and LR72, where the numbers denote the net effective sky coverages as a percentage, i.e., We note that the highestlatitude and smallest region, LR16, is an addition to those defined in PIPXXX, following the same procedure. Furthermore, LR63 has been split into its north and south Galactic hemisphere portions, yielding LR63N and LR63S, covering and 0.30, respectively. All of these regions are shown in Fig. 1.
As previously stated in PIPXXX, the dust polarization angular power spectra computed on these regions can be considered to be approximately statistically independent because most of the power arises from the brightest 10% of a given region, the same 10% that differentiates one (preliminary) region from another.
More specifically, for the power spectra computed with the 353 GHz Planck data, the difference (nonoverlapping) region LRx_{i}−LRx_{i−1}, i> 1 (about 10% of the sky, by definition), contains more than 75% of the power computed on LRx_{i}, for x ∈ { 24,33,42,53,63,72 }. LR63N and LR63S are independent of one another, but not of LR63. Nevertheless, we prefer to utilize the overlapping regions (which are not fully statistically independent) in order to be able to relate the present work to the previous Planck analyses in general and to PIPXXX in particular. These masks are publicly available through the Planck Legacy Archive^{2}.
2.3. Simulations
The simulated polarization maps presented in this work were built using a simplified twocomponent model consisting of dust plus CMB, both simulated as stochastic realizations with a Gaussian distribution. CMB maps are defined as realizations based on the Planck Collaboration bestfit ΛCDM model (Planck Collaboration XIII 2016), assuming a tensortoscalar ratio r = 0 and an optical depth τ = 0.06 (Planck Collaboration Int. XLVII 2016; Planck Collaboration Int. XLVI 2016). The dust component maps at 353 GHz are defined as Gaussian realizations using a powerlaw model of the EE and BB angular power spectra (with a spectral index equal to −0.42 and amplitudes matching those in Table 1 of PIPXXX), following the prescriptions of PIPXXX, and normalized for each region of the sky introduced in Sect. 2.2. Dust maps at other frequencies are scaled with a constant modified blackbody spectrum with and T_{d} = 19.6 K (Planck Collaboration Int. XVII 2014; Planck Collaboration Int. XXII 2015). The results of this paper have also been confronted with realistic dust simulations (nonGaussian, described in Vansyngel et al. 2016) finding statistical significances and probabilitiestoexceed (PTEs) similar to the ones obtained with the Gaussian realizations of the dust.
Fig. 1
Masks and complementary science regions that retain fractional coverage of the sky, f_{sky}, from 0.8 to 0.2 (see details in Sect. 2.2). The grey is the CO mask, whose complement is a selected region with f_{sky} = 0.8. In increments of f_{sky} = 0.1, the retained regions can be identified by the colours from yellow (0.3) to blue (0.8). The LR63N and LR63S regions are displayed in black and white, respectively, and the LR63 region is the union of the two. Also shown is the (unapodized) point source mask used. 
The instrumental noise component is then introduced in each pixel using the Planck Q and U covariance maps at 217 and 353 GHz, associated with the DS and HM data setups. Using a set of 1000 simulations, we checked that the instrumental noise in the range of multipoles ℓ = 50–700 built from covariances was consistent with the Full Focal Plane Monte Carlo noise simulations, namely FFP8 (Planck Collaboration XII 2016) at 217 and 353 GHz, which are publicly available through the Planck Legacy Archive^{2}. For each of the regions described in Sect. 2.2, we built 1000 independent dust, CMB and noise realizations with these properties.
We note that we did not directly use the FFP8 maps in this analysis because of two issues with the polarized dust component in the Planck Sky Model (PSM, Delabrouille et al. 2013) used to build these simulations. First, the polarized dust maps include Planck instrumental noise because they are computed from the Planck 353 GHz maps smoothed to 30′. This contribution significantly enhances the apparent dust signal at high Galactic latitudes, where the signaltonoise ratio of the data is low at this resolution. Second, the polarized dust component in the PSM already includes spatial variations of the polarized dust spectral index. Thus, it cannot be used to perform a nulltest for no spatial variation of the SED.
3. The dust polarization correlation ratio
3.1. Definition
We statistically estimate the spatial variation of the dust SED by constructing the correlation ratio for the dust between 217 and 353 GHz, as a function of the multipole ℓ. This is defined as the ratio of the crossspectrum between these bands and the geometric mean of the two autospectra in the same bands, i.e., (1)where X ∈ { E,B }. If the maps at 217 and 353 GHz contain only dust and the dust SED is constant over the region for which the power spectra are computed, then . However, if the dust is not the only component to contribute to the sky polarization or if the dust SED varies spatially, then the ratio is expected to deviate from unity. Nevertheless, the spatial variations of the SED, as we see in the next subsection, do not affect the ratio at the largest scales. This is the reason why the study in PIPXXX (where the ratio was computed from the fitted amplitude of power laws dominated by the largest scales) found no significant deviation from unity. In what follows we conduct this analysis for ℓ> 50 to avoid any significant contribution from Planck systematics effects at low ℓ (Planck Collaboration VIII 2016).
In order to avoid any bias issues due to the noise autocorrelation, and in order to minimize the systematic effects, the auto and crossspectra are computed using independent sets of Planck data, namely the halfmission and detectorset maps (Planck Collaboration VIII 2016). Hence the autospectra are computed for frequency ν as (2)where and are the two independent sets of data, e.g., DS1 and DS2, or HM1 and HM2 (see Sect. 2.1). Similarly, the crossspectra between two frequencies, ν_{1} and ν_{2}, are given by (3)where the indices i and j take the values 1 and 2. The results are labelled “HM” or “DS” when obtained using the halfmission or detectorset maps, respectively.
The spectra have been computed using Xpol (Tristram et al. 2005), which is a polarization pseudoC_{ℓ} estimator that corrects for incomplete sky coverage and pixel and beam window functions.
3.2. Planck measurements
Fig. 2
Dust polarization correlation ratios (top panel) and (bottom panel) between 217 and 353 GHz, computed on the LR63 region. The correlation ratios determined from the detectorset data (DS) splits are shown as yellow squares and the ratios computed from the halfmission (HM) splits as red diamonds. The correlation ratio expectations for the model described in Eq. (5) are shown as blue dashed lines (blue circles when binned as the data). Horizontal blue segments in the top panel represent the range of the ℓ bins. The uncertainties have been estimated as the median absolute deviation over a set of 1000 simulations (see Sect. 2.3) of CMB, dust, and Gaussian noise. 
The EE and BB dust correlation ratios obtained in four multipoles ranges (ℓ = [50,160], [160,320], [320,500], and [500,700]) are shown in Fig. 2 for the Planck HM and DS data versions, and a fraction of the sky f_{sky} = 0.7, i.e., the LR63 region.
Since the definition of the correlation ratio, , uses nonlinear operations (such as ratio and square root), the associated uncertainties are not trivial to estimate and so are determined from the median absolute deviation of 1000 Monte Carlo noise realizations, including Planck instrumental noise, as detailed in Sect. 2.3.
The correlation ratios between the 217 and 353 GHz bands that might be expected are computed using the simulation setup introduced in Sect. 2.3, and defined by a twocomponent modelling of dust plus CMB signals, assuming no spatial variations of the dust SED and no noise. If we suppose that the Planck maps at 217 and 353 GHz are a sum of CMB and dust components, then (in thermodynamic units and assuming no instrumental noise) we have
where M_{dust} is the dust map at 353 GHz, M_{CMB} is the CMB map, α is a constant scaling coefficient representing the dust SED, and M can represent the Stokes parameters Q and U. Then, combining Eqs. (1) and (4) and assuming that the dust and CMB components are not spatially correlated, the expected correlation ratio becomes (5)where X ∈ { E,B }. It can clearly be seen that even if α is a constant, the CMB component will make the correlation ratio ℛ ≠ 1. The model power spectra corresponding to Eq. (5) are also shown in Fig. 2.
The PlanckEE data match the expected EE correlation ratio and are strongly dominated by the CMB signal. Two approaches have been considered for removing the CMB component from the correlation ratio in order to see the effect of the dust decorrelation in the EE spectra, namely analysis in pixel space or in multipole space. The noise on the CMB template, subtracted from the Planck maps, would produce an autocorrelation of the noise when building the correlation ratios and would strongly impact our analysis in polarization; this argues against using the first (pixelbased) option. Moreover, the second option, which consists of correcting the 217 and 353 GHz Planck crossspectra by subtracting a model of the CMB power spectrum, is affected by the cosmic variance of the CMB, which is dominant compared to the dust component in the EE correlation ratios. For all these reasons, in the following analysis we focus on the multipolebased BB modes only, where the CMB component (coming from the lensing Bmodes, Planck Collaboration Int. XLI 2016) is subdominant compared to the observed signal.
As can be seen in the lower panel of Fig. 2, the BB correlation ratios, , exhibit a clear deficit compared to the expected model; this is discussed in more detail in the next section. We have also carried out similar analyses of the other regions defined in Sect. 2.2. The same deficit can be seen in Fig. 3, where the results averaged over the lowest multipole bin are shown as a function of the mean column density (see Table 1).
3.3. Significance of the Planck measurements
Fig. 3
Dependence of the Planck correlation ratios with the mean column density of the region on which they are computed in the first multipole bin (ℓ = [50,160]). The HM and DS measurements can be compared to the theoretical expectations (including dust and CMB, as in Eq. (5)) shown as blue segments. The grey bars give the 68% and 95% confidence levels computed over 1000 realizations of dust, CMB, and Gaussian noise. The column density dependence fit by a power law (see Eq. (6)) is shown as a dashed line. 
The significance of the Planck data correlation ratios can be quantified using the distribution of computed with simulations including Planck instrumental noise at 217 and 353 GHz. Since the denominator of can get close to zero, the distribution of the correlation ratios from dust, CMB, and noise can be highly nonGaussian. For this reason, we cannot give symmetrical error bars and significance levels expressed in terms of σ. Instead we use the PTE, which makes no assumption about the shape of the distribution. As we will see, the Planck measurements of the correlation ratios in the DS and HM data sets appear systematically smaller than the most probable correlation ratio in our simulations that include instrumental noise.
The impact of the CMB and the noise on the correlation ratio can be seen in Fig. 3 via the grey vertical histograms. These data are based on the correlation ratios obtained on a set of 1000 Monte Carlo realizations, including Gaussian Planck noise on top of the simulated CMB and dust components (see Sect. 2.3). The distributions of the simulated correlation ratios are available for all multipole ranges and regions in Appendix A (in Figs. A.1–A.9). While noise barely affects the most probable ratios in the first multipole bin (ℓ = 50–160) when compared to the theoretical expectation (blue), it can create an important level of decorrelation in the other bins, particularly in the smallest regions (up to an additional 30% in the fourth multipole bin for LR16).
In order to quantify the significance of the Planck measurements with respect to the simulated decorrelation from CMB and noise, we compute the PTE, defined as the probability of a simulation having more decorrelation (i.e., smaller correlation ratio) than the data. This is computed as the fraction of the 1000 realizations having a correlation ratio smaller than the Planck measurements, for each multipole bin and HM/DS case. These PTEs are listed in Table 1. We also compute a combined PTE over all multipole bins, defined as the probability of obtaining correlation ratios smaller than Planck measurements simultaneously in all four multipole bins (50 < ℓ < 700).
The combined PTEs (for 50 < ℓ < 700) are < 1.5% for the DS case and <0.1% for the HM case. When focusing on individual multipole bins, the detection level is not as strong; the PTE values range between 0.6% and 11.2% in the first multipole bin (50 < ℓ < 160) for both DS and HM cases. In the second multipole bin (160 < ℓ < 320), the DS correlation ratio PTEs range from 11.5 to 47.4%, while for the HM case they range from 1.1 to 39.9% (with significant PTEs on several regions though). The third and fourth multipole bins (320 <ℓ < 700) show no significant evidence of decorrelation.
The significant excess of decorrelation in the Planck data between 217 and 353 GHz, especially in the first multipole bin (50 < ℓ < 160) is consequently very unlikely to be attributable to CMB or to instrumental noise (or to systematic effects, see next section). We therefore conclude that this excess is a statistical measurement of the spatial variation of the polarized dust SED.
3.4. Impact of systematic effects
Since we have restricted our analysis to multipoles ℓ > 50, the 217 and 353 GHz crossspectra are not significantly affected by those systematic effects that are most important at low multipoles, such as the ADC nonlinearity correction or the dipole and calibration uncertainties (Planck Collaboration VII 2016; Planck Collaboration VIII 2016). However, the Planck crossspectra in the multipole range 50 < ℓ < 700 could be affected by beam systematics. Thanks to its definition, the correlation ratio should be approximately independent of the beam uncertainty, because of the presence of the same beam functions, and , in the numerator and denominator. A further systematic contribution could arise from the difference between the beam function of the 353 × 217 crossspectra, , and the product of the independent beam functions, . We have checked that this ratio exhibits a very low departure from unity, at the 10^{5} level, which cannot reproduce the amplitude of the observed Planck correlation ratio. We also checked that the correction of the bandpass mismatch in the 217 and 353 GHz bands does not affect the correlation ratio. The same analysis has been reproduced using two versions of the bandpass mismatch corrections (see Planck Collaboration VII 2016), yielding results consistent down to 0.1%.
We use the two splits, HM and DS, as indicators of the level of residuals due to systematic effects in our analysis. This can be assessed by examining Figs. A.1 to A.9. While the HM and DS correlation ratios are very consistent in the first and last multipole bins (ℓ = 50–160, and ℓ = 500–700), they are not as consistent for the second multipole bin, (ℓ = 160–320). This apparent discrepancy is not explained by the current knowledge of any systematic effects in Planck, and so indicates the need for some caution.
3.5. Dependence on column density
For CMB polarization studies, it is important to characterize the dependence of the observed decorrelation ratio on column density. The Planck correlation ratios in the first multipole bin (ℓ = [50,160]), obtained on the various science regions, are shown in Fig. 3 as a function of the mean column density, N_{H i}, computed for each region as the average over the unmasked pixels of the Planck column density map, assuming a constant opacity τ/N_{H i} (Planck Collaboration XI 2014). The HM and DS measurements of the PlanckBB correlation ratio can be compared to the theoretical expectation (blue segment) and the dispersion due to noise (grey histograms) computed over 1000 Monte Carlo realizations including Gaussian noise (see Sect. 2.3). We recall (see Sect. 2.2) that measurements of the correlation ratio obtained in different regions can be considered statistically independent to a good approximation (except for LR63 with respect to LR63N or LR63S).
In this first multipole bin (50 < ℓ < 160), where the primordial Bmode signal is expected, the BB correlation ratio of the DS and HM cases can be described well by a power law of N_{H i}: (6)with and γ = −1.76 ± 0.69. Hence the more diffuse the Galactic foregrounds, the stronger the decorrelation between 217 and 353 GHz. This is an important issue for CMB analyses, which mainly focus on the most diffuse regions of the sky in order to minimize the contamination by Galactic dust emission.
4. Impact on the CMB Bmodes
Fig. 4
spectra likelihood posteriors on the r and A_{d} parameters derived from the simulations of the BICEP2Keck and Planck 353 GHz data, using: 150 and 353 GHz channels (BKP, red); and 95, 150 and 353 GHz channels (BK(+95), blue). The 1D posteriors of r (marginalized over and A_{d}) and of A_{d} (marginalized over and r) are shown in the left and middle panels (blue and red lines, respectively). The input value in our simulations for the dust amplitude at 353 GHz (4.5 μK^{2} at ℓ = 80) is indicated as a dashed line. The 2D posterior marginalized over is presented in the right panel (68% in darker shading and 95% in lighter shading). 
In Sect. 3 we show that the correlation of the dust polarization Bmodes between 217 and 353 GHz can depart significantly from unity for multipoles ℓ ≳ 50. This decorrelation, as previously noted by, e.g., Tassis & Pavlidou (2015) or BKP15, will have an impact on the search for CMB primordial Bmodes that assume a constant dust polarization SED over the region of sky considered. In order to quantify this effect, we use toy model simulations of the BICEP2/Keck and Planck data, introducing a decorrelation between the 150 and the 353 GHz channels that matches our results in Sect. 3. This analysis is intended to be illustrative only, and not to be an exact reproduction of the work presented in BICEP2/Keck Array and Planck Collaborations (2015) and BICEP2 and Keck Array Collaborations (2016).
We approximate the likelihood analyses presented in BKP15 and BICEP2 and Keck Array Collaborations (2016) using simple simulations of the CMB and dust angular power spectra. We directly simulate the 150 × 150, 150 × 353, and 353 × 353 angular power spectra (where the Planck 353 GHz spectrum comes from two noiseindependent detectorset subsamples). In order to include the noise contribution from these experiments and to perform a Monte Carlo analysis over 2000 simulations, we add the sample variance and the noise contribution to the spectra as a Gaussian realization of (7)for each power spectra C_{b} computed in an ℓcentred bin of size Δℓ_{b}, where is the effective sky fraction, and are the signal plus noise autopower spectra in the frequency bands ν_{1} and ν_{2}, respectively, and is the signalonly crosspower spectrum between these two frequencies (supposing the noise to be uncorrelated between the two bands); although the noise affects the variance of the C_{ℓ}, we note that it does not affect its mean value, since we crosscorrelate noiseindependent maps.
The CMB spectrum is generated from a Planck 2015 bestfit ΛCDM model (Planck Collaboration XIII 2016) with no tensor modes (r = 0). The dust for the 353 × 353 power spectrum is constructed as a power law in ℓ, specifically ℓ^{0.42} following PIPXXX. The amplitude of this spectrum at ℓ = 80 is taken to be A_{d} = 4.5 μK^{2}. This is an ad hoc value chosen to lie between the predicted PIPXXX value in the BICEP2 region (A_{d} = 13.4 ± 0.26 μK^{2}), and the BKP15 value (, marginally compatible with our chosen value, which could be underestimated if some decorrelation exists). A single modified blackbody spectrum is applied to scale the dust 353 GHz spectrum to the other frequencies, with and T_{d} = 19.6 K.
Finally, we introduce a decorrelation factor in our simulated crossspectra, which we chose to be constant in ℓ. If we make the assumption that the SED spatial variations come from spatial variations of the dust spectral index around its mean value (see Sect. 5.1), a firstorder expansion gives a frequency dependence of (1−ℛ^{BB}) that scales as [ln(ν_{1}/ν_{2})] ^{2}. We explored many values for the ratio and we have chosen to present here the results we obtain for a correlation ratio between 150 and 353 GHz of . With the frequency scaling of [ln(ν_{1}/ν_{2})] ^{2}, this ratio becomes , , and .
We construct a threeparameter likelihood function , similar to the one used in BKP15, with a Gaussian prior on . For each simulation, the posteriors on r and A_{d} are marginalized over and we construct the final posterior as the histogram over 2000 simulations of the individual maximum likelihood values for r and A_{d}.
Our results when approximating the BKP15 analysis are presented in Fig. 4. The maximum likelihood values are r = 0.046 ± 0.036, or r< 0.12 at 95% CL, and A_{d} = 3.23 ± 0.85 μK^{2}. The bias on r is higher than the value assessed in section V.A of BKP15 (we tested that our simulations give the same 0.018 bias on r when using ), given that we introduce more decorrelation and that our dust amplitude is higher. The value we find for A_{d} is similar to that of the BKP15 analysis.
Finally, we repeat the same analysis on simulations corresponding to BICEP2 and Keck Array Collaborations (2016), where we add the 95 GHz data and increase the sensitivity in the 150 GHz channel with respect to BKP15. Unlike the analysis used in the BICEP2 and Keck Array Collaborations (2016), our analysis does not parametrize a synchrotron component. The posteriors from these simulations, labelled “BK(+95)”, are also shown in Fig. 4. The maximum likelihood values are r = 0.014 ± 0.027 (or r< 0.07 at 95% CL) and A_{d} = 3.54 ± 0.77 μK^{2}. Even without a synchrotron component in the model of the data, the positions of the peak in the posteriors on r and A_{d} with respect to BKP15 are shifted in the same direction as found in BICEP2 and Keck Array Collaborations (2016).
The region of the sky observed by the BICEP and Keck instruments presented in BKP15 has a mean column density of N_{H i} = 1.6 × 10^{20} cm^{2}, very similar to the value for our region LR24. In the multipole range 50 <ℓ< 160, using the empirical relation derived in Sect. 3.5 (Eq. (6)), we expect to have a decorrelation between 217 and 353 GHz. Introducing the latter decorrelation value in our simple simulations shifts the r posterior towards higher values (up to r = 0.1), making a decorrelation as high as the one we measure on the LR24 mask very unlikely, given the BKP15 data. This shows the limitation of the empirical relation we derived in Sect. 3.5 when dealing with small regions, since the properties of the decorrelation might be very variable over the sky. A specific analysis of these data could quantitatively confirm the amount of decorrelation that is already allowed or excluded by the data.
The results presented in this section stress that a decorrelation between the dust polarization at any two frequencies will result in a positive bias in the r posterior, in the absence of an appropriate modelling in the likelihood parametrization (even if β_{d} is fixed) or in any component separation. The current BICEP2/Keck and Planck limits on r still leave room for a decorrelation of the dust polarization among frequencies that could be enough to lead to spurious detections for future StageIII or StageIV CMB experiments (see, e.g., Abazajian et al. 2016).
5. Discussion
We now quantify how the observed decorrelation of the BB power spectrum between the 217 and 353 GHz bands can be explained by spatial variations of the polarized dust SED using two toy models presented in Sects. 5.1 and 5.2. This simplified description in characterizing spatial variations of the dust SED is a first step, which ignores correlations between dust properties and the structure of the magnetized interstellar medium. Correlations between matter and the Galactic magnetic field have been shown to be essential to account for statistical properties of dust polarization at high Galactic latitudes (Clark et al. 2015; Planck Collaboration Int. XXXVIII 2016). In Sect. 5.3, we use the framework introduced by Planck Collaboration Int. XLIV (2016) to discuss why such correlations are also likely to be an essential element of any physical account of the variations of the dust SED in polarization.
5.1. Spectral index variations
In a first approach, we assume that the variations of the polarized dust SED can be fully explained by spatial variations of the polarized dust spectral index applied simultaneously to the Stokes Q and U components.
We make a simplifying approximation by assuming that the polarized dust spectral index follows a Gaussian distribution centred on the mean value (Planck Collaboration Int. XXII 2015), with a single dispersion, , at all scales over the whole sky. Specifically the dust polarization spectral index is given by (8)where is a Gaussian distribution centred on x_{0} with a standard deviation σ, and is the dispersion of the spectral index map, defined as the standard deviation after smoothing at a resolution of 1°. The dust temperature is kept constant over the sky and equal to T_{d} = 19.6 K (Planck Collaboration Int. XXII 2015).
A main caveat of this approach comes from the power introduced by the spatial variations of the spectral index, which alters the power spectrum of dust polarization. This effect remains small between 217 and 353 GHz for low values of at ℓ< 150, but can lead to dust power spectra that are inconsistent with Planck observations when extrapolated to further bands. This effect is inherent to this modelling approach, which must only be considered as illustrative. However, here the simulated maps at 217 and 353 GHz are built from maps at an intermediate frequency ( GHz) to minimize the addition of power.
The BB correlation ratio model between 217 and 353 GHz is constructed as follows. We start with a set of Q and U dust template maps at 277 GHz, appropriately normalized to Planck data as detailed in Sect. 2.3. The polarization dust maps at 217 and 353 GHz are extrapolated from 277 GHz using a Gaussian realization of the polarized dust spectral index, given a level of the dispersion . These simulated maps do not include noise at this point. The correlation ratio model is finally obtained by averaging 100 realizations of the correlation ratio computed on a pair of simulated maps.
This model is illustrated in Fig. 5 for three nonzero values of and is compared to the Planck HM and DS measurements in the LR42 region. In order to match the Planck data (including uncertainties), an indicative value of around 0.07 is suggested by this simple analysis.
Fig. 5
Correlation ratio modelled with Gaussian spatial variations of the polarized dust spectral index, in LR42 as in Fig. 2. The HM and DS Planck measurements are shown as diamonds and squares, respectively. The model is plotted for four values of , which is defined as the standard deviation of the Gaussian realization of the spectral index after smoothing at 1°. The first model, with , is described in Sect. 3.2. 
We compare our estimate of the spectral index variations for dust polarization to those measured for the total dust intensity. We use the Commander (Eriksen et al. 2006, 2008) dust component maps derived from a modified blackbody fit to Planck data (Planck Collaboration X 2016), providing two separate maps of dust temperature and dust spectral index. From the ratio I_{353}/I_{217} computed from these two maps, we derive an equivalent allsky map of the intensity dust spectral index, , assuming a constant dust temperature of T_{d} = 19.6 K. We use the two halfmission maps of the dust spectral index, computed from 1° resolution maps, instead of the fullsurvey map, in order to reduce the impact of noise and systematic effects when computing the covariance, and we derive an estimate of the standard deviation of the dust spectral index in intensity, in LR42, about half the value measured for polarization in the same region. We note that this value is not corrected for the contribution of the data noise and the anisotropies of the cosmic infrared background (Planck Collaboration Int. XXII 2015). Both of these contribution are much smaller than the empirical value of 0.17 found in Planck Collaboration Int. XXII (2015), which is dominated by noise.
5.2. Polarization angle variations
In a second approach, we assume that the decorrelation of the BB power spectrum between 217 and 353 GHz can be explained by spatial variations of the polarization angle, keeping the dust temperature and polarized spectral index constant over the whole sky. Unlike the first model, this modelling approach conserves the total power in the power spectrum. It is motivated by the nature of the polarized signal, which can be considered as the sum of spin2 quantities over multiple components with varying spectral dependencies along the line of sight. Physical interpretation of polarization angle variations are discussed further in Sect. 5.3.
The spatial variations of the polarization angle are assumed to follow a circular normal distribution (or von Mises distribution) around 0, given by (9)where ℐ_{0}(x) is the modified Bessel function of order 0, and κ is analogous to 1 /σ^{2} for the normal distribution. While the circular normal distribution allows us to define random angles in the range [− π,π], we rescale the realizations of angular variations to match the definition for the range of polarization angles lying between −π/ 2 and π/ 2. Again we start from a set of Q and U dust maps at 277 GHz, which are then extrapolated to 217 and 353 GHz following a modified blackbody spectrum using T_{d} = 19.6 K and . The polarization pseudovector obtained from Q and U is rotated independently at each frequency using two different realizations of the polarization angle variations. The correlation ratio model is then obtained by averaging 100 realizations of the correlation ratio computed from a pair of simulated maps.
This model is illustrated Fig. 6 for three finite values of κ in the LR42 region. The case κ = 2, which matches the Planck data quite well at ℓ ~ 100, represents a 1 σ dispersion of the polarization angle of about 2° between the 217 and 353 GHz polarization maps, after smoothing at 1° resolution.
Fig. 6
As in Fig. 5, but illustrating the correlation ratio modelled with Gaussian spatial variations of the polarization angle. The model is plotted for four values of κ, which sets the level of the circular normal distribution for the polarization angle, analogous to 1 /σ^{2} for the normal distribution. 
5.3. Origin of the spatial variations of the polarized SED
The spatial variations of the Galactic dust SED are found to be greater for dust polarization than for dust intensity. This difference could result from the lack of correlation between spectral variation and emission features in our modelling in Sect. 5.1. It could also reflect the different nature of the two observables. While variations in the dust SED tend to average out in intensity (because it is a scalar quantity), they may not average out as much in polarization (because it is a pseudovector). In other words, dust polarization depends on the magnetic field structure, while dust intensity does not. In this section, we give further details of this interpretation of the decorrelation.
The analysis of Planck data offers several lines of evidence for the imprint of interstellar magnetohydrodynamical (MHD) turbulence on the dust polarization sky. Firstly, Planck Collaboration Int. XIX (2015) has reported an anticorrelation between the polarization fraction and the local dispersion of the polarization angle. Secondly, large filamentary depolarization patterns observed in the Planck 353 GHz maps are associated, in most cases, with large, local, fluctuations of the polarization angles related to the magnetic field structure (Planck Collaboration Int. XIX 2015). Thirdly, the polarization fraction of the Planck 353 GHz map shows a large scatter at high and intermediate latitudes, which can be interpreted as line of sight depolarization associated with interstellar MHD turbulence (Planck Collaboration Int. XX 2015; Planck Collaboration Int. XLIV 2016). Lastly, several studies have reported a clear trend where the magnetic field is locally aligned with the filamentary structure of interstellar matter (Martin et al. 2015; Clark et al. 2015; Planck Collaboration Int. XXXII 2016; Planck Collaboration Int. XXXVIII 2016; Kalberla et al. 2016).
As presented in Planck Collaboration Int. XLIV (2016), the Planck polarized data at high Galactic latitudes can be modelled using a small number of effective polarization layers in each of which the Galactic magnetic field has a turbulent component of the magnetic field with a distinct orientation. Within this model, the integration of dust polarization along the line of sight can be viewed as the result of an oriented random walk in the (Q, U) plane, with a small number of steps tending towards the mean direction of the Galactic magnetic field. The magnetic field orientation of each turbulent layer sets the direction of the step, while the dust polarization properties, including the efficiency of grain alignment, sets their length. Changes in dust properties generate differences in the relative lengths of the steps across frequency, and thereby differences in the polarization fraction and angle, as illustrated graphically in Fig. 7. In this framework, variations of the dust SED along the line of sight impact both the polarization fraction and the polarization angle, because the structure of the magnetic field and of diffuse interstellar matter are correlated (Planck Collaboration Int. XXXII 2016; Planck Collaboration Int. XXXVIII 2016). In the random walk, the variations in length and angle of the polarization pseudovector have the most impact for a small number of steps. Planck Collaboration Int. XLIV (2016) related the number of layers in their model to the density structure of the diffuse interstellar medium and to the correlation length of the turbulent component of the magnetic field along the line of sight.
Fig. 7
Schematic illustration of the decomposition of the line of sight complex polarization vector (in the (Q, U) space) into a random walk process through four turbulent layers, shown for two neighbouring pixels at two frequencies, ν_{1} and ν_{2}. The integrated polarization pseudovectors are affected by the polarization angle fluctuations, leading to a decorrelation between the two frequencies. 
A quantitative modelling of this perspective on dust polarization is required to assess its ability to reproduce the results of our analysis of the Planck data and its impact on the dust/CMB componentseparation task. To our knowledge, none of the studies carried out so far to quantify limits set by polarized dust foregrounds on future CMB experiments in order to search for primordial Bmodes (except for the work of Tassis & Pavlidou 2015) have considered the frequency decorrelation that arises from the interplay between dust properties and interstellar MHD turbulence.
6. Conclusions
We have used the 217 and 353 GHz Planck 2015 data to investigate the spatial variation of the polarized dust SED at high Galactic latitudes. We computed – the BB crossspectrum between the Planck polarization maps in these bands divided by the geometric mean of the two autospectra in the same bands – over large regions of the sky at high latitudes, for multipoles 50 <ℓ< 700. The value of was computed with distinct sets of data to control systematics.
The ratio has been shown to be significantly lower than what is expected purely from the presence of CMB and noise in the data with a confidence, estimated using data simulations, higher than 99%. We interpret this result as evidence of significant spatial variations of the dust polarization SED.
In the multipole bin 50 <ℓ< 160 that encompasses the recombination bump of the primordial Bmode signal, values are consistent for the distinct Planck data splits we used. The measured values exhibit a systematic trend with column density, where decreases with decreasing mean column densities N_{H i} as , with γ = −1.76 ± 0.69 and . This suggests that, statistically speaking, the cleaner a sky area is from Galactic foregrounds, the more challenging it may be to extrapolate dust polarization from submillimetre to CMB bands.
The spatial variations of the dust SED are shown to be stronger in polarization than in intensity. This difference may reflect the interplay, encoded in the data, between the polarization properties of dust grains, including grain alignment, the Galactic magnetic field, and the density structure of interstellar matter.
We have proposed two toy models to quantify the observed decorrelation, one based on spatial Gaussian variations of the polarized dust spectral index and the other on spatial Gaussian variations of the polarization angle. Both models reproduce the general trend observed in Planck data, with reasonable accuracy given the noise level of the Planck measurements. They represent a first step in characterizing variations of the dust SED, as a prelude to further work on a physically motivated model, which will take into account the expected correlation between the spatial variations of the SED and structures in the dust polarization maps.
Spatial variations of the dust SED can lead to biased estimates of the tensortoscalar ratio, r, because of inaccurate extrapolation of dust polarization when cleaning or modelling the Bmode signal at microwave frequencies. As an illustration, we have shown that in a region such as that observed by the BICEP2Keck Array, a decorrelation of 15% between the 353 and 150 GHz bands could result in a likelihood posterior on the tensortoscalar ratio of r = 0.046 ± 0.036, comparable to the joint BICEP2 and Keck Array/Planck Collaboration result.
It appears essential now to place tighter constraints on the spectral dependence of polarized dust emission in the submillimetre in order to properly propagate the information on the Galactic foregrounds into the CMB bands. More specifically, the spatial variations of the polarized dust SED need to be mapped and modelled with improved accuracy in order to be able to confidently reach a level of Galactic dust residual lower than r ~ 10^{2}.
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/INSUIN2P3INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MINECO, J.A., 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/planckcollaboration. The research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP7/20072013)/ERC grant agreement No. 267934.
References
 Abazajian, K. N., Adshead, P., Ahmed, Z., et al. 2016, ArXiv eprints [arXiv:1610.02743] [Google Scholar]
 ArmitageCaplan, C., Dunkley, J., Eriksen, H. K., & Dickinson, C. 2012, MNRAS, 424, 1914 [NASA ADS] [CrossRef] [Google Scholar]
 BICEP2 and Keck Array Collaborations 2016, Phys. Rev. Lett., 116, 031302 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 BICEP2/Keck Array and Planck Collaborations 2015, Phys. Rev. Lett., 114, 101301 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Choi, S. K., & Page, L. A. 2015, J. Cosmol. Astropart. Phys., 12, 020 [NASA ADS] [CrossRef] [Google Scholar]
 Clark, S. E., Hill, J. C., Peek, J. E. G., Putman, M. E., & Babler, B. L. 2015, Phys. Rev. Lett., 115, 241302 [NASA ADS] [CrossRef] [Google Scholar]
 Creminelli, P., López Nacir, D., Simonović, M., Trevisan, G., & Zaldarriaga, M. 2015, J. Cosmol. Astropart. Phys., 11, 31 [NASA ADS] [CrossRef] [Google Scholar]
 Delabrouille, J., Betoule, M., Melin, J.B., et al. 2013, A&A, 553, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eriksen, H. K., Dickinson, C., Lawrence, C. R., et al. 2006, ApJ, 641, 665 [NASA ADS] [CrossRef] [Google Scholar]
 Eriksen, H. K., Jewell, J. B., Dickinson, C., et al. 2008, ApJ, 676, 10 [NASA ADS] [CrossRef] [Google Scholar]
 Errard, J., Feeney, S. M., Peiris, H. V., & Jaffe, A. H. 2016, JCAP, 03, 052 [NASA ADS] [CrossRef] [Google Scholar]
 Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 [NASA ADS] [CrossRef] [Google Scholar]
 Kalberla, P. M. W., Kerp, J., Haud, U., et al. 2016, ApJ, 821, 117 [NASA ADS] [CrossRef] [Google Scholar]
 Krachmalnicoff, N., Baccigalupi, C., Aumont, J., Bersanelli, M., & Mennella, A. 2016, A&A, 588, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Martin, P. G., Blagrave, K. P. M., Lockman, F. J., et al. 2015, ApJ, 809, 153 [NASA ADS] [CrossRef] [Google Scholar]
 Planck Collaboration XI. 2014, A&A, 571, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration I. 2016, A&A, 594, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VII. 2016, A&A, 594, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VIII. 2016, A&A, 594, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration X. 2016, A&A, 594, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XII. 2016, A&A, 594, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIII. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXV. 2016, A&A, 594, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XVII. 2014, A&A, 566, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XIX. 2015, A&A, 576, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XX. 2015, A&A, 576, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XXII. 2015, A&A, 576, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XXX. 2016, A&A, 586, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XXXII. 2016, A&A, 586, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XXXVIII. 2016, A&A, 586, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XLI. 2016, A&A, 596, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XLIV. 2016, A&A, 596, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XLVI. 2016, A&A, 596, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XLVII. 2016, A&A, 596, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Remazeilles, M., Dickinson, C., Eriksen, H. K. K., & Wehus, I. K. 2016, MNRAS, 458, 2032 [NASA ADS] [CrossRef] [Google Scholar]
 Tassis, K., & Pavlidou, V. 2015, MNRAS, 451, L90 [NASA ADS] [CrossRef] [Google Scholar]
 Tristram, M., MacíasPérez, J. F., Renault, C., & Santos, D. 2005, MNRAS, 358, 833 [NASA ADS] [CrossRef] [Google Scholar]
 Vansyngel, F., Boulanger, F., Ghosh, T., et al., 2016, A&A, submitted [arXiv:1611.02577] [Google Scholar]
Appendix A: Distribution of the from Planck simulations
We present in this appendix the distribution of the ratio for our simulations of the Planck 217 and 353 GHz maps, including Gaussian CMB, dust, and Planck noise, which are described in Sect. 3.3. Since the quantity is a ratio, we expect distributions that are strongly nonGaussian, even for Gaussian simulations, as soon as the denominator gets close to zero.
In Figs. A.1–A.9, we show the histogram of the correlation ratio from 1000 realizations of detectorset and halfmission simulations. The histogram “occurrences” are the fractions of simulations that fall in each bin. These are compared to the values obtained from the Planck data. Since we compute only noiseindependent crossspectra that can take negative values, when has an imaginary denominator, the data value is absent. The PTE for each data set, computed as the percentage of simulations that have a smaller than the data, is given in Table 1.
Fig. A.1
Distribution of the correlation ratios on the LR16 region from Planck simulations for detectorset (DS, red histograms) and halfmission (HM, orange histograms) data splits, for the four multipole bins we use in the analysis: top left, ℓ ∈ [50,160]; top right, ℓ ∈ [160,320]; bottom left, ℓ ∈ [320,500]; and bottom right, ℓ ∈ [500,700]. The values derived from the data are shown as vertical lines of the corresponding colour. The PTEs corresponding to these data values with respect to the simulations are reported in Table 1. 
All Tables
All Figures
Fig. 1
Masks and complementary science regions that retain fractional coverage of the sky, f_{sky}, from 0.8 to 0.2 (see details in Sect. 2.2). The grey is the CO mask, whose complement is a selected region with f_{sky} = 0.8. In increments of f_{sky} = 0.1, the retained regions can be identified by the colours from yellow (0.3) to blue (0.8). The LR63N and LR63S regions are displayed in black and white, respectively, and the LR63 region is the union of the two. Also shown is the (unapodized) point source mask used. 

In the text 
Fig. 2
Dust polarization correlation ratios (top panel) and (bottom panel) between 217 and 353 GHz, computed on the LR63 region. The correlation ratios determined from the detectorset data (DS) splits are shown as yellow squares and the ratios computed from the halfmission (HM) splits as red diamonds. The correlation ratio expectations for the model described in Eq. (5) are shown as blue dashed lines (blue circles when binned as the data). Horizontal blue segments in the top panel represent the range of the ℓ bins. The uncertainties have been estimated as the median absolute deviation over a set of 1000 simulations (see Sect. 2.3) of CMB, dust, and Gaussian noise. 

In the text 
Fig. 3
Dependence of the Planck correlation ratios with the mean column density of the region on which they are computed in the first multipole bin (ℓ = [50,160]). The HM and DS measurements can be compared to the theoretical expectations (including dust and CMB, as in Eq. (5)) shown as blue segments. The grey bars give the 68% and 95% confidence levels computed over 1000 realizations of dust, CMB, and Gaussian noise. The column density dependence fit by a power law (see Eq. (6)) is shown as a dashed line. 

In the text 
Fig. 4
spectra likelihood posteriors on the r and A_{d} parameters derived from the simulations of the BICEP2Keck and Planck 353 GHz data, using: 150 and 353 GHz channels (BKP, red); and 95, 150 and 353 GHz channels (BK(+95), blue). The 1D posteriors of r (marginalized over and A_{d}) and of A_{d} (marginalized over and r) are shown in the left and middle panels (blue and red lines, respectively). The input value in our simulations for the dust amplitude at 353 GHz (4.5 μK^{2} at ℓ = 80) is indicated as a dashed line. The 2D posterior marginalized over is presented in the right panel (68% in darker shading and 95% in lighter shading). 

In the text 
Fig. 5
Correlation ratio modelled with Gaussian spatial variations of the polarized dust spectral index, in LR42 as in Fig. 2. The HM and DS Planck measurements are shown as diamonds and squares, respectively. The model is plotted for four values of , which is defined as the standard deviation of the Gaussian realization of the spectral index after smoothing at 1°. The first model, with , is described in Sect. 3.2. 

In the text 
Fig. 6
As in Fig. 5, but illustrating the correlation ratio modelled with Gaussian spatial variations of the polarization angle. The model is plotted for four values of κ, which sets the level of the circular normal distribution for the polarization angle, analogous to 1 /σ^{2} for the normal distribution. 

In the text 
Fig. 7
Schematic illustration of the decomposition of the line of sight complex polarization vector (in the (Q, U) space) into a random walk process through four turbulent layers, shown for two neighbouring pixels at two frequencies, ν_{1} and ν_{2}. The integrated polarization pseudovectors are affected by the polarization angle fluctuations, leading to a decorrelation between the two frequencies. 

In the text 
Fig. A.1
Distribution of the correlation ratios on the LR16 region from Planck simulations for detectorset (DS, red histograms) and halfmission (HM, orange histograms) data splits, for the four multipole bins we use in the analysis: top left, ℓ ∈ [50,160]; top right, ℓ ∈ [160,320]; bottom left, ℓ ∈ [320,500]; and bottom right, ℓ ∈ [500,700]. The values derived from the data are shown as vertical lines of the corresponding colour. The PTEs corresponding to these data values with respect to the simulations are reported in Table 1. 

In the text 
Fig. A.2
Same as Fig. A.1, for the LR24 region. 

In the text 
Fig. A.3
Same as Fig. A.1, for the LR33 region. 

In the text 
Fig. A.4
Same as Fig. A.1, for the LR42 region. 

In the text 
Fig. A.5
Same as Fig. A.1, for the LR53 region. 

In the text 
Fig. A.6
Same as Fig. A.1, for the LR63N region. 

In the text 
Fig. A.7
Same as Fig. A.1, for the LR63S region. 

In the text 
Fig. A.8
Same as Fig. A.1, for the LR63 region. 

In the text 
Fig. A.9
Same as Fig. A.1, for the LR72 region. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.