Issue |
A&A
Volume 586, February 2016
|
|
---|---|---|
Article Number | A133 | |
Number of page(s) | 25 | |
Section | Cosmology (including clusters of galaxies) | |
DOI | https://doi.org/10.1051/0004-6361/201425034 | |
Published online | 09 February 2016 |
Planck intermediate results
XXX. The angular power spectrum of polarized dust emission at intermediate and high Galactic latitudes
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
African Institute for Mathematical Sciences,
6–8 Melrose Road, Muizenberg,
Cape Town, South
Africa
3
Agenzia Spaziale Italiana Science Data Center,
via del Politecnico snc,
00133
Roma,
Italy
4
Agenzia Spaziale Italiana, viale Liegi 26, 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 KwaZulu-Natal,
Westville Campus, Private Bag
X54001, 4000
Durban, South
Africa
8
Atacama Large Millimeter/submillimeter Array, ALMA Santiago Central
Offices, Alonso de Cordova 3107,
Vitacura, Casilla
763 0355
Santiago,
Chile
9
CITA, University of Toronto, 60 St. George St., Toronto, ON
M5S 3H8,
Canada
10
CNRS, IRAP, , 9
Av. colonel Roche, BP
44346, 31028
Toulouse Cedex 4,
France
11
CRANN, Trinity College, Dublin, Ireland
12
California Institute of Technology, Pasadena, California, USA
13
Centre for Theoretical Cosmology, DAMTP, University of
Cambridge, Wilberforce
Road, Cambridge
CB3 0WA,
UK
14
Centro de Estudios de Física del Cosmos de Aragón
(CEFCA), Plaza San Juan, 1, planta
2, 44001
Teruel,
Spain
15
Computational Cosmology Center, Lawrence Berkeley National
Laboratory, Berkeley,
California,
USA
16
Consejo Superior de Investigaciones Científicas
(CSIC), Madrid,
Spain
17
DSM/Irfu/SPP, CEA-Saclay, 91191
Gif-sur-Yvette Cedex,
France
18
DTU Space, National Space Institute, Technical University of
Denmark, Elektrovej
327, 2800
Kgs. Lyngby,
Denmark
19
Département de Physique Théorique, Université de
Genève, 24 Quai E.
Ansermet, 1211
Genève 4,
Switzerland
20
Departamento de Física, Universidad de Oviedo,
Avda. Calvo Sotelo s/n,
Oviedo,
Spain
21
Department of Astronomy and Astrophysics, University of
Toronto, 50 Saint George Street,
Toronto, Ontario,
Canada
22
Department of Astrophysics/IMAPP, Radboud University
Nijmegen, PO Box
9010, 6500 GL
Nijmegen, The
Netherlands
23
Department of Physics & Astronomy, University of British
Columbia, 6224 Agricultural Road,
Vancouver, British
Columbia, Canada
24
Department of Physics and Astronomy, Dana and David Dornsife College
of Letter, Arts and Sciences, University of Southern California,
Los Angeles, CA
90089,
USA
25
Department of Physics and Astronomy, University College
London, London
WC1E 6BT,
UK
26
Department of Physics, Florida State University,
Keen Physics Building, 77 Chieftan
Way, Tallahassee,
Florida,
USA
27
Department of Physics, Gustaf Hällströmin katu 2a, University of
Helsinki, Helsinki,
Finland
28
Department of Physics, Princeton University,
Princeton, New Jersey, USA
29
Department of Physics, University of California,
Berkeley, California, USA
30
Department of Physics, University of California,
One Shields Avenue, Davis, California, USA
31
Department of Physics, University of California,
Santa Barbara, California, USA
32
Department of Physics, University of Illinois at
Urbana-Champaign, 1110 West Green
Street, Urbana,
Illinois,
USA
33
Dipartimento di Fisica e Astronomia G. Galilei, Università degli
Studi di Padova, via Marzolo
8, 35131
Padova,
Italy
34
Dipartimento di Fisica e Scienze della Terra, Università di
Ferrara, via Saragat
1, 44122
Ferrara,
Italy
35
Dipartimento di Fisica, Università La Sapienza,
P.le A. Moro 2, Roma, Italy
36
Dipartimento di Fisica, Università degli Studi di
Milano, via Celoria,
16
Milano,
Italy
37
Dipartimento di Fisica, Università degli Studi di
Trieste, via A. Valerio
2, Trieste,
Italy
38
Dipartimento di Fisica, Università di Roma Tor
Vergata, via della Ricerca Scientifica
1, Roma,
Italy
39
Dipartimento di Matematica, Università di Roma Tor
Vergata, via della Ricerca
Scientifica 1, Roma,
Italy
40
Discovery Center, Niels Bohr Institute, Blegdamsvej 17, Copenhagen, Denmark
41
Dpto. Astrofísica, Universidad de La Laguna (ULL),
38206, La Laguna, Tenerife,
Spain
42
European Southern Observatory, ESO Vitacura, Alonso de Cordova 3107, Vitacura, Casilla
19001
Santiago,
Chile
43
European Space Agency, ESAC, Planck Science Office, Camino bajo del
Castillo, s/n, Urbanización
Villafranca del Castillo, Villanueva de la Cañada, Madrid, Spain
44
European Space Agency, ESTEC, Keplerlaan 1, 2201 AZ
Noordwijk, The
Netherlands
45 Facoltà di Ingegneria, Università
degli Studi e-Campus, via Isimbardi
10, 22060
Novedrate ( CO), Italy
46
Gran Sasso Science Institute, INFN,viale F. Crispi 7, 67100
L’Aquila,
Italy
47
HGSFP and University of Heidelberg, Theoretical Physics
Department, Philosophenweg
16, 69120
Heidelberg,
Germany
48
Haverford College Astronomy Department, 370 Lancaster Avenue, Haverford, Pennsylvania, USA
49
Helsinki Institute of Physics, Gustaf Hällströmin katu 2, University
of Helsinki, Helsinki,
Finland
50
INAF–Osservatorio Astronomico di Padova,
Vicolo dell’Osservatorio 5,
Padova,
Italy
51
INAF–Osservatorio Astronomico di Roma, via di Frascati 33, Monte Porzio Catone,
Italy
52
INAF–Osservatorio Astronomico di Trieste,
via G.B. Tiepolo 11, Trieste, Italy
53
INAF/IASF Bologna, via Gobetti
101, Bologna, Italy
54
INAF/IASF Milano, via E. Bassini 15, Milano, Italy
55
INFN, Sezione di Bologna, via Irnerio 46, 40126
Bologna,
Italy
56
INFN, Sezione di Roma 1, Università di Roma Sapienza,
Piazzale Aldo Moro 2,
00185
Roma,
Italy
57
INFN, Sezione di Roma 2, Università di Roma Tor
Vergata, via della Ricerca
Scientifica 1, Roma,
Italy
58
INFN/National Institute for Nuclear Physics,
via Valerio 2, 34127
Trieste,
Italy
59
IPAG: Institut de Planétologie et d’Astrophysique de Grenoble,
Université Grenoble Alpes, CNRS, IPAG, 38000
Grenoble,
France
60
IUCAA, Post Bag 4, Ganeshkhind, Pune University
Campus, 411 007
Pune,
India
61
Imperial College London, Astrophysics group, Blackett
Laboratory, Prince Consort
Road, London,
SW7 2AZ,
UK
62
Infrared Processing and Analysis Center, California Institute of
Technology, Pasadena,
CA
91125,
USA
63
Institut Universitaire de France, 103 Bd Saint-Michel, 75005
Paris,
France
64
Institut d’Astrophysique Spatiale, CNRS (UMR 8617) Université
Paris-Sud 11, Bâtiment
121, Orsay,
France
65
Institut d’Astrophysique de Paris, CNRS (UMR 7095),
98bis Boulevard Arago,
75014
Paris,
France
66
Institute for Space Sciences, Bucharest-Magurale,
Romania
67
Institute of Astronomy, University of Cambridge,
Madingley Road, Cambridge
CB3 0HA,
UK
68 Institute of Theoretical
Astrophysics, University of Oslo, Blindern, Oslo,
Norway
69
Instituto de Astrofísica de Canarias, C/Vía Láctea s/n, La Laguna, Tenerife, Spain
70
Instituto de Física de Cantabria (CSIC-Universidad de
Cantabria), Avda. de los Castros
s/n, Santander,
Spain
71
Istituto Nazionale di Fisica Nucleare, Sezione di
Padova, via Marzolo
8, 35131
Padova,
Italy
72
Jet Propulsion Laboratory, California Institute of
Technology, 4800 Oak Grove
Drive, Pasadena,
California,
USA
73
Jodrell Bank Centre for Astrophysics, Alan Turing Building, School
of Physics and Astronomy, The University of Manchester, Oxford Road, Manchester, M13
9PL, UK
74
Kavli Institute for Cosmology Cambridge,
Madingley Road, Cambridge, CB3 0HA, UK
75
LAL, Université Paris-Sud, CNRS/IN2P3, Orsay, France
76
LAPTh, Univ. de Savoie, CNRS, BP110, 74941
Annecy-le-Vieux,
France
77
LERMA, CNRS, Observatoire de Paris, 61 Avenue de l’Observatoire, Paris, France
78
Laboratoire AIM, IRFU/Service d’Astrophysique – CEA/DSM – CNRS –
Université Paris Diderot, Bât. 709,
CEA-Saclay, 91191
Gif-sur-Yvette Cedex,
France
79
Laboratoire Traitement et Communication de l’Information, CNRS (UMR
5141) and Télécom ParisTech, 46 rue
Barrault, 75634
Paris Cedex 13,
France
80
Laboratoire de Physique Subatomique et de Cosmologie, Université
Joseph Fourier Grenoble I, CNRS/IN2P3, Institut National Polytechnique de
Grenoble, 53 rue des
Martyrs, 38026
Grenoble Cedex,
France
81
Laboratoire de Physique Théorique, Université Paris-Sud 11 &
CNRS, Bâtiment 210,
91405
Orsay,
France
82
Lawrence Berkeley National Laboratory, Berkeley, California, USA
83
Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Str. 1, 85741
Garching,
Germany
84
McGill Physics, Ernest Rutherford Physics Building, McGill
University, 3600 rue
University, Montréal,
QC, H3A 2T8, Canada
85
National University of Ireland, Department of Experimental
Physics, Maynooth,
Co. Kildare,
Ireland
86
Niels Bohr Institute, Blegdamsvej 17, Copenhagen, Denmark
87
Observational Cosmology, Mail Stop 367-17, California Institute of
Technology, Pasadena,
CA, 91125, USA
88
Optical Science Laboratory, University College London,
Gower Street, London, UK
89
SB-ITP-LPPC, EPFL, 1015
Lausanne,
Switzerland
90
SISSA, Astrophysics Sector, via Bonomea 265, 34136
Trieste,
Italy
91
School of Physics and Astronomy, Cardiff University,
Queens Buildings, The Parade,
Cardiff, CF24 3AA, UK
92
School of Physics and Astronomy, University of
Nottingham, Nottingham
NG7 2RD,
UK
93
Space Research Institute (IKI), Russian Academy of
Sciences, Profsoyuznaya Str,
84/32, 117997
Moscow,
Russia
94
Space Sciences Laboratory, University of California,
Berkeley, California, USA
95
Special Astrophysical Observatory, Russian Academy of
Sciences, Nizhnij Arkhyz,
Zelenchukskiy region, 369167
Karachai-Cherkessian Republic,
Russia
96
Sub-Department of Astrophysics, University of Oxford,
Keble Road, Oxford
OX1 3RH,
UK
97 Theory Division, PH-TH,
CERN, 1211,
Geneva 23,
Switzerland
98
UPMC Univ. Paris 6, UMR7095, 98bis Boulevard Arago, 75014
Paris,
France
99
Université de Toulouse, UPS-OMP, IRAP, 31028
Toulouse Cedex 4,
France
100
Universities Space Research Association, Stratospheric Observatory
for Infrared Astronomy, MS
232-11, Moffett
Field, CA
94035,
USA
101
University Observatory, Ludwig Maximilian University of
Munich, Scheinerstrasse
1, 81679
Munich,
Germany
102
University of Granada, Departamento de Física Teórica y del Cosmos,
Facultad de Ciencias, Granada, Spain
103
University of Granada, Instituto Carlos I de Física Teórica y
Computacional, Granada,
Spain
104
Warsaw University Observatory, Aleje Ujazdowskie 4, 00-478
Warszawa,
Poland
Received: 19 September 2014
Accepted: 1 December 2014
The polarized thermal emission from diffuse Galactic dust is the main foreground present in measurements of the polarization of the cosmic microwave background (CMB) at frequencies above 100 GHz. In this paper we exploit the uniqueness of the Planck HFI polarization data from 100 to 353 GHz to measure the polarized dust angular power spectra CℓEE and CℓBB over the multipole range 40 <ℓ< 600 well away from the Galactic plane. These measurements will bring new insights into interstellar dust physics and allow a precise determination of the level of contamination for CMB polarization experiments. Despite the non-Gaussian and anisotropic nature of Galactic dust, we show that general statistical properties of the emission can be characterized accurately over large fractions of the sky using angular power spectra. The polarization power spectra of the dust are well described by power laws in multipole, Cℓ ∝ ℓα, with exponents αEE,BB = −2.42 ± 0.02. The amplitudes of the polarization power spectra vary with the average brightness in a way similar to the intensity power spectra. The frequency dependence of the dust polarization spectra is consistent with modified blackbody emission with βd = 1.59 and Td = 19.6 K down to the lowest Planck HFI frequencies. We find a systematic difference between the amplitudes of the Galactic B- and E-modes, CℓBB/CℓEE = 0.5. We verify that these general properties are preserved towards high Galactic latitudes with low dust column densities. We show that even in the faintest dust-emitting regions there are no “clean” windows in the sky where primordial CMB B-mode polarization measurements could be made without subtraction of foreground emission. Finally, we investigate the level of dust polarization in the specific field recently targeted by the BICEP2 experiment. Extrapolation of the Planck 353 GHz data to 150 GHz gives a dust power 𝒟ℓBB ≡ ℓ(ℓ+1)CℓBB/(2π) of 1.32 × 10-2 μKCMB2 over the multipole range of the primordial recombination bump (40 <ℓ< 120); the statistical uncertainty is ± 0.29 × 10-2 μKCMB2 and there is an additional uncertainty (+0.28, −0.24) × 10-2 μKCMB2 from the extrapolation. This level is the same magnitude as reported by BICEP2 over this ℓ range, which highlights the need for assessment of the polarized dust signal even in the cleanest windows of the sky.
Key words: cosmic background radiation / cosmology: observations / ISM: structure / ISM: magnetic fields / polarization
© ESO, 2016
1. Introduction
The sky at high Galactic latitude and frequencies above about 100 GHz is dominated by thermal emission from the Galactic interstellar medium, specifically arising from dust grains of about 0.1 μm in size. Asymmetrical dust grains align with the Galactic magnetic field to produce polarized emission. This polarized submillimetre emission has been measured from ground-based and balloon-borne telescopes (e.g., Hildebrand et al. 1999; Benoît et al. 2004; Ponthieu et al. 2005; Vaillancourt 2007; Matthews et al. 2014). The observed polarization relates to the nature, size, and shape of dust grains and the mechanisms of alignment, discussed for example by Draine (2004) and Martin (2007). It also probes the structure of the Galactic magnetic field, which is an essential component of models of Galactic dust polarization (Baccigalupi 2003; Fauvet et al. 2011, 2012; O’Dea et al. 2012; Jaffe et al. 2013; Delabrouille et al. 2013).
The polarized emission from dust is also of interest in the context of foregrounds (Tucci et al. 2005; Dunkley et al. 2009a; Gold et al. 2011) to the cosmic microwave background (CMB). On angular scales between 10′ and a few tens of degrees, cosmological B-mode polarization signals may be present that were imprinted during the epoch of inflation. The discovery of a primordial B-mode polarization signature is a major scientific goal of many CMB experiments. These include ground-based experiments (ACTPol, Niemack et al. 2010; BICEP2, BICEP2 Collaboration 2014a; Keck-array, Staniszewski et al. 2012; POLARBEAR, Arnold et al. 2010; QUBIC, Ghribi et al. 2014; QUIJOTE, Rubiño-Martín et al. 2010; and SPTpol, Austermann et al. 2012), stratospheric balloon missions (EBEX, Grainger et al. 2008; and SPIDER, Fraisse et al. 2013), and the ESA Planck1 satellite (Tauber et al. 2010). Accurate assessment and, if necessary, subtraction of foreground contamination is critical to the measurement of CMB E- and B-mode polarization because the expected signals from inflation and late-time reionization are expected to be small.
Planck has measured the all-sky dust polarization at 353 GHz, where the
dust emission dominates over other polarized signals. These data have been presented in a
first set of publications in which the focus was on the structure of the Galactic magnetic
field and the characterization of dust polarization properties (Planck Collaboration Int. XIX 2015; Planck Collaboration Int. XX 2015; Planck
Collaboration Int. XXI 2015; Planck Collaboration
Int. XXII 2015). Here, we use the Planck polarized data to compute
the and
power spectra of dust polarization over the
multipole range 40 <ℓ<
600, on large fractions of the sky away from the Galactic plane. We also
investigate dust polarization in sky patches at high Galactic latitude with sizes comparable
to those surveyed by ground-based CMB experiments. We derive statistical properties of dust
polarization from these spectra, characterizing the shape of the spectra and their amplitude
with respect to both the observing frequency and the mean dust intensity of the sky region
over which they are computed. We verify that these properties hold in low-column-density
patches at high Galactic latitude and we explore statistically the potential existence of
“clean” patches on the sky that might be suitable for cosmology.
Our analysis of dust polarization is relevant to the present generation of CMB polarization observations, and to the design of future experiments. It gives a statistical description of Galactic dust polarization, providing input for the modelling of Galactic dust as part of component separation methods and for CMB polarization likelihood analysis parameterization.
The BICEP2 collaboration has recently reported a significant detection of the B-mode power spectrum around the expected angular scale of the recombination bump, namely a few degrees (BICEP2 Collaboration 2014a,b). Their analysis was based on dust polarization models that predicted subdominant contamination of their B-mode signal by dust polarization. We use information from our detailed analysis of Planck polarization data at 353 GHz to assess the potential dust contamination.
The paper is organized as follows. In Sect. 2, we
present the Planck HFI polarization data used in this work and describe the
general properties of the polarization maps in terms of emission components and systematic
effects. In Sect. 3, we describe our method for
computing the dust and
angular power spectra, including the
selected science regions of interest on the sky. We assess and compare the two methods we
use to compute the power spectra in Appendix A. In
Sect. 4, we present power spectra of dust polarization
for multipoles ℓ>
40, computed with high signal-to-noise ratio (S/N) on large fractions of
the sky, and characterize their shape and amplitude. Complementary angular power spectra
involving temperature and polarization,
and
, and cross-polarization,
are given in Appendix B. We extend this analysis to smaller sky patches at high Galactic
latitude in Sect. 5. In Appendix C we discuss some complementary aspects of this analysis of patches.
These results are used specifically in Sect. 5.3 to
build a map of the expected dust contamination of the
power spectrum at 150 GHz and
ℓ = 80. In
Sect. 6 we present a study of the polarized dust
emission in the vicinity of the BICEP2 region. Systematic effects relating to the
Planck angular power spectrum estimates are assessed in Appendix D and we discuss the decorrelation of the dust signal
between frequencies in Appendix E. Section 7 summarizes the main conclusions and discusses the
implications of this work for future CMB experiments2.
2. Planck polarization maps
2.1. Planck data
The Planck Collaboration recently released the Planck satellite nominal mission temperature data and published a set of papers describing these data and their cosmological interpretation (e.g., Planck Collaboration I 2014; Planck Collaboration XVI 2014). These results are based on the data from the two instruments on board the satellite (Low Frequency Instrument, LFI, Mennella et al. 2011, and High Frequency Instrument, HFI, Planck HFI Core Team 2011). The data processing of the nominal mission data (Surveys 1 and 2, 14 months) was summarized in Planck Collaboration II (2014) and Planck Collaboration VI (2014).
Planck HFI measures the linear polarization at 100, 143, 217, and 353 GHz (Rosset et al. 2010). The properties of the detectors (sensitivity, spectral response, noise properties, beams, etc.) are described in detail in Lamarre et al. (2010) and their in-flight performance is reported in Planck HFI Core Team (2011), Planck Collaboration VII (2014), Planck Collaboration VIII (2014), Planck Collaboration IX (2014), and Planck Collaboration X (2014), while Planck Collaboration VI (2014) describes the general processing applied to the data to measure polarization. In this paper, we make use of full-mission (Surveys 1 to 5, 30 months, Planck Collaboration I 2014) polarization maps of Planck HFI (internal data release “DX11d”), projected into the HEALPix pixelization scheme (Górski et al. 2005). This is one of the first publications to use these maps, which will be described in the Planck cosmology 2015 release.
To compute polarization angular power spectra, we use Q and U maps at 100, 143, 217, and 353 GHz. Specifically, we calculate power spectra using the so-called “Detector-Set” maps (hereafter “DetSets”), constructed using two subsets of polarization sensitive bolometers (PSBs) at a given frequency (see Table 3 in Planck Collaboration VI 2014). Each DetSet polarization map is constructed using data from two pairs of PSBs, with the angle between the two PSBs in a pair being 90°, and the angle between pairs being 45°. In this paper we concentrate on the Q and U maps at 353 GHz. The Stokes Q and U maps at lower frequencies (100, 143, and 217 GHz) are only used to determine the spectral energy distribution (SED) of the dust emission in polarization.
To quantify systematic effects, we also use maps made from other data subsets (Planck Collaboration VI 2014). We use the ring halves (hereafter “HalfRing”), where the approximately 60 circles performed for each Planck telescope ring (also called a stable pointing period) are divided into two independent subsets of 30 circles. We use observational years (hereafter “Years”), consisting of Surveys 1 and 2 on the one hand and Surveys 3 and 4 on the other, to build two additional maps with independent noise.
The Planck maps we use are in thermodynamic units (KCMB). To characterize the SED of the dust emission in polarization we express the data as the specific intensity (such as that for Stokes I dust emission, Id(ν)) at the Planck reference frequencies, using the conversion factors and colour corrections from Planck Collaboration IX (2014)3. For the average dust SED at intermediate Galactic latitudes, the colour correction factor is 1.12 at 353 GHz (see Table 3 in Planck Collaboration Int. XXII 2015).
In addition to these basic products, a Planck CO map from Planck Collaboration XIII (2014), the so-called “type 3” map, and the Planck 857 GHz map, are also used in the selection of the large intermediate latitude analysis regions (see Sect. 3.3.1).
2.2. Emission contributions to the Planck HFI polarization maps
2.2.1. Polarized thermal dust emission
Thermal dust emission is partially linearly polarized (e.g., Hildebrand et al. 1999; Benoît et al. 2004; Ponthieu et al. 2005; Vaillancourt 2007). It is the dominant polarized foreground signal in the high frequency Planck bands (Tucci et al. 2005; Dunkley et al. 2009b; Fraisse et al. 2009; Fauvet et al. 2011; Planck Collaboration Int. XXII 2015).
Dust polarization arises from alignment of non-spherical grains with the interstellar magnetic field (e.g., Hildebrand 1988; Draine 2004; Martin 2007). The structure of the dust polarization sky has already been described using maps of the polarization fraction (p) and angle (ψ) derived from the Planck HFI 353 GHz data (Planck Collaboration Int. XIX 2015; Planck Collaboration Int. XX 2015). The map of p shows structure on all scales, with polarization fractions ranging from low (less than 1%) to high values (greater than 18%). Planck Collaboration Int. XIX (2015) and Planck Collaboration Int. XX (2015) report an anti-correlation between p and the local dispersion of ψ, which indicates that variations in p arise mainly from depolarization associated with changes in the magnetic field orientation within the beam, rather than from changes in the efficiency of grain alignment.
Planck Collaboration Int. XXII (2015) showed
that the SED of polarized dust emission over the four Planck HFI
frequencies from 100 to 353 GHz is consistent with a modified blackbody emission law of
the type Id(ν) ∝
νβdBν(Td),
with spectral index βd = 1.59 for Td = 19.6
K 4, and where Bν is the Planck
function. About 39% of the sky at intermediate Galactic latitudes was analysed5. Among 400 circular patches with 10° radius (equivalent to a sky fraction
) the 1σ dispersion of
βd was 0.17 for constant Td = 19.6 K.
We scale this uncertainty on βd to larger sky areas by using the
factor
. This is a conservative choice because
this uncertainty includes the effects of noise in the data and so is an upper limit to
the true regional variations of βd on this scale. This polarization
spectral index can be compared to variations in the spectral index
for the intensity SED. For that quantity
the S/N of the data is higher than for polarization and Planck Collaboration Int. XXII (2015) report a dispersion of 0.07
(1σ) over
the same sized circular patches. Planck Collaboration
Int. XVII (2014) extend this analysis for intensity to high Galactic latitudes
in the southern Galactic cap, using the dust-H i correlation to separate the
faint emission of dust from the anisotropies of the cosmic infrared background, and find
a dispersion of about 0.10 in
. We expect spectral variations to be
correlated in polarization and intensity, unless the dust emission has a significant
component that is unpolarized.
2.2.2. CMB
The CMB temperature anisotropies have been measured with unprecedented accuracy by the Planck Collaboration (Planck Collaboration I 2014; Planck Collaboration XV 2014), and preliminary Planck polarization results have been demonstrated to be in very good agreement with the cosmology inferred from temperature measurements (Planck Collaboration I 2014; Planck Collaboration XVI 2014).
For , the ΛCDM concordance model has been shown to
be a very good fit to all the available data (including preliminary
Planck results at ℓ ≳ 50; see Barkats et al. 2014, for a recent compendium). For 353 GHz data at small
angular scales (ℓ ≳
400), the E-mode CMB polarization is comparable to the
power of dust polarization at high Galactic latitudes.
The CMB B-mode power, even for the highest primordial tensor
perturbation models, is negligible with respect to the dust polarization at 353 GHz.
Since no reliable published CMB polarization maps are available, we have chosen not to
remove the CMB polarization from the Planck HFI Q and U maps. Nevertheless,
because the CMB E-mode polarization is significant with respect to
the dust at 353 GHz at high multipoles (and even at lower multipoles for the lower
frequencies), when studying the Planck HFI bands we subtract the
Planck best-fit ΛCDM model (Col. 2 of Table 2 in Planck Collaboration XVII 2014) from the dust power
spectra, paying the price of an increased error due to sample variance. No CMB is
removed in this work when computing the dust
spectra.
2.2.3. Synchrotron emission
Synchrotron emission is known to be significantly polarized (up to 75% for typical relativistic electron spectra, Rybicki & Lightman 1979). Since its specific intensity scaling with frequency follows a power law with a spectral index close to −3 (Gold et al. 2011; Macellari et al. 2011; Fuskeland et al. 2014), synchrotron polarized emission is expected to be subdominant in the Planck HFI channels in general and negligible at 353 GHz (Tucci et al. 2005; Dunkley et al. 2009a; Gold et al. 2011; Fauvet et al. 2011; Fuskeland et al. 2014; Planck Collaboration Int. XXII 2015). Hence, we neither subtract nor mask any synchrotron contribution before estimating the angular power spectra of dust polarization. The justification of this assumption will be demonstrated below by studying the frequency dependence of the polarized dust power between 100 and 353 GHz (but see also Appendix D.4).
2.2.4. Polarized point sources
Radio sources have been shown to have a fractional polarization of a few percent (e.g., Battye et al. 2011; Massardi et al. 2013).Their contribution to the polarization angular power spectra in the Planck HFI bands is expected to be negligible at low and intermediate multipoles (Battye et al. 2011; Tucci & Toffolatti 2012). Upper limits have been set on the polarization of infrared galaxies, and their contribution to the polarization power spectra is also expected to be negligible (e.g., Seiffert et al. 2007). However, the brightest of the polarized point sources can be responsible for ringing in the angular power spectra estimation, and therefore need to be masked (see Sect. 3.3.1).
2.2.5. CO emission
The first three carbon monoxide (CO) Galactic emission lines at 115 GHz (J = 1 → 0), 230 GHz (J = 2 → 1), and 345 GHz (J = 3 → 2), contribute significantly to the power in the Planck HFI bands at 100, 217, and 353 GHz, respectively (Planck Collaboration IX 2014). The Planck data were used to produce the first all-sky maps of Galactic CO emission (Planck Collaboration XIII 2014). It is known that CO emission can be intrinsically polarized (Goldreich & Kylafis 1982; Li & Henning 2011). Furthermore, CO emission can induce spurious polarization because of the differences in spectral transmission at the CO frequencies between Planck HFI detectors (Planck Collaboration IX 2014). For these reasons, we mask CO-emitting regions (Sect. 3.3.1). Outside this mask, as has been shown in Planck Collaboration XIII (2014) by comparing the Planck CO maps to high Galactic latitude ground-based CO observations (Hartmann et al. 1998; Magnani et al. 2000), the CO emission is negligible in the Planck channels (lower than one fourth of each channel noise rms at 95% CL). We will check that our polarization analysis is not contaminated by CO by examining the frequency dependence of the polarized angular power spectra of the dust (see Sect. 4.5).
2.3. Systematics of the Planck HFI polarization maps
The first CMB polarization results from Planck were presented in Planck Collaboration I (2014) and Planck Collaboration XVI (2014). The EE power spectrum at ℓ> 200 was found to be consistent with the cosmological model derived from temperature anisotropies. Systematic effects in the data have so far limited the use of Planck HFI polarization data on large angular scales. The polarization systematics in the 2013 data were discussed and estimated in Planck Collaboration VI (2014; see the power spectra shown in their Fig. 27). The same data were used in the Planck Galactic polarization papers (Planck Collaboration Int. XIX 2015; Planck Collaboration Int. XX 2015; Planck Collaboration Int. XXI 2015; Planck Collaboration Int. XXII 2015) and there the low brightness regions of the 353 GHz sky were masked.
In this paper, we use a new set of Planck polarization maps for which the systematic effects have been significantly reduced. Corrections will be fully described and applied in the Planck 2015 cosmology release, as well as the remaining systematic effects that we describe briefly here.
Two main effects have been corrected in the time-ordered data prior to mapmaking. The
correction for the nonlinearity of the analogue-to-digital converters was improved, and we
have also corrected the data for very long time constants that were not previously
identified (see Planck Collaboration VI 2014; Planck Collaboration X 2014). After these corrections,
the main systematic effects over the multipole range 40 <ℓ< 600 relevant
to this analysis result from leakage of intensity into polarization maps.Effects arising
from polarization angle and polarization efficiency uncertainties have been shown to be
second order (Planck Collaboration VI 2014). The
leakage effect can be expressed as (1)where ν is the frequency, and
, the leakage pattern for bolometer
b, is
multiplied by the different s leakage source maps
and their associated scaling coefficients
. The leakage patterns are fully determined
by the scanning strategy. They represent the cumulative result of all the systematic
effects that lead to a leakage of intensity to polarization. In the
Planck HFI bands, there are three main sources of intensity to
polarization leakage: (i) monopole differences between detectors not corrected by data
destriping (the intensity source term is constant over the sky); (ii) bolometer
inter-calibration mismatch (the source term is the full intensity map, including the CMB
dipole and the Galactic emission); and (iii) a dust spectral mismatch term. The spectral
bandpass varies from one bolometer to another for a given band. As the bolometer gain is
calibrated on the CMB dipole, the differential gains on the dust emission produce the
bandpass mismatch term (for which the source term is the dust intensity map).
Results from a global fit of the Q and U maps with these three leakage terms
are used to quantify the leakage. This fit
yields estimates of the scaling coefficients at each frequency ν in Eq. (1), which allows us to compute angular power
spectra of the leakage terms in Sect. 4. We point
out that the fit captures any emission in the maps that has a pattern on the sky similar
to one of the leakage patterns. These global fit maps are used to assess the level of
systematics, but are not removed from the data.
An independent and complementary estimate of systematic effects in the data is provided by the null tests that we can build from the DetSets, HalfRings, and Years data subsets (see Sect. 2.1). These null tests are a good way of determining the level of any systematics other than intensity to polarization leakage. These are pursued in Sect. 4.1 and Appendices C.1 and D.3.
3. Computation of angular power spectra of polarized dust emission
3.1. Methods
We can use the Planck data to compute the polarization angular power
spectra ( and
) of the polarized dust emission within
selected sky regions. Even if the statistical properties of the dust emission on the sky
might not be entirely captured by a two-point function estimator, because the scope of
this paper is to assess the level of dust in the framework of CMB data analysis, we follow
the approximation that is generally made when processing such data, i.e., that a large
fraction of the information is contained in power spectra.
On an incomplete sky, a polarization field can be divided into three classes of modes: pure E-modes; pure B-modes; and “ambiguous” modes, which are a mixture of the true E- and B-modes (Bunn et al. 2003).
The ambiguous modes represent a cross-talk between E and B, which is often referred
to as “E-to-B leakage” for the CMB (because for the CMB
). Methods used to estimate the CMB angular
power spectrum for polarization account for and correct analytically for the incomplete
sky coverage. However, the presence of the ambiguous modes yields a biased estimate of the
variance of the spectra, unless so-called “pure” power spectrum estimators are used (Smith 2006). For dust polarization, the power in
E- and
B-modes are
comparable, and we do not expect a significant variance bias.
The two specific approaches we use are Xpol, our main method, which we describe as a “classical” pseudo-Cℓ estimator6 and, for comparison, Xpure, a pure pseudo-Cℓ estimator. Since they give similar results for the present study, as we demonstrate below, the former is chosen as our main method because it is computationally less expensive.
3.1.1. Xpol
Xpol is an extension of the Xspect method (Tristram et al. 2005) to polarization. Xspect computes the pseudo power spectra and corrects them for incomplete sky coverage, filtering effects, and pixel and beam window functions. Correction for incomplete sky coverage is performed using a Master-like algorithm (Hivon et al. 2002), consisting of the inversion of the mode-mode coupling matrix Mℓℓ′ that describes the effect of the partial sky coverage; Mℓℓ′ is computed directly from the power spectrum of the mask that selects the data in the analysis region of interest (Sect. 3.3). The HEALPix pixel window functions (Górski et al. 2005) are used to correct for pixelization effects. For the beams we use the Planck HFI individual detector beam transfer functions described in Planck Collaboration VII (2014).
Xpol estimates the error bars analytically, without requiring Monte Carlo simulations.
Using simulated data comprising inhomogeneous white noise and a Gaussian map with a dust
power spectrum, we have checked that this analytical estimate is not biased for
multipoles ℓ>
40. The analytical error bars combine the contributions from
instrumental noise and sample variance. The Gaussian approximation of the sample
variance is (2)where X = { E,B
}, fsky is the retained sky fraction
(which can be
if the sky field is apodized), and
Δℓbin is the size of the multipole bin
ℓbin. To estimate the error relevant to
the polarized dust signal measured within a given region, we subtract quadratically this
estimate of the contribution from sample variance. The performance of Xpol on Gaussian
simulations that have EE and BB dust-like angular power spectra is presented
in Appendix A.
3.1.2. Xpure
Xpure is a numerical implementation of the pure pseudo-spectral approach described and validated in Grain et al. (2009). The method is optimized for computing CMB B-mode power spectra over small sky patches. It uses a suitably chosen sky apodization that vanishes (along with its first derivative) at the edges of the patch in order to minimize the effects of E-to-B leakage. For the estimation of the angular power spectra of the Planck data we compute the cross-correlation of two different DetSets. Uncertainties are obtained by performing Monte Carlo inhomogeneous white noise simulations, considering the diagonal terms of the pixel-pixel covariance for the two data sets.
We have applied this algorithm to Planck 353 GHz maps and simulations to estimate the B-mode power spectrum of the Galactic thermal dust emission in regions ranging from 1% to 30% of the sky. We have used Xpure here as a cross-check for the robustness of the Xpol method presented in the previous section. Validation and comparison of the performance of the two algorithms is presented in Appendix A on simulations and in Appendix D.1 on application to the data.
3.2. Computing cross-spectra
To avoid a bias arising from the noise, we compute all the Planck power
spectra from cross-correlations of DetSets maps (see Sect. 2). The noise independence of the two DetSet maps at a given frequency was
quantified in Planck Collaboration XV (2014) and
the resulting level of noise bias in the cross-power spectra between them has been shown
to be negligible. The cross-power spectrum at a given frequency ν is
(3)where
and
are the two independent DetSet 1 and
DetSet 2 maps at the frequency ν. The Planck cross-band spectrum
between the frequencies ν and ν′ is
(4)or equivalently the cross-spectrum between
ν and
ν′ of the
averaged frequency maps
.
As we noted in Sect. 2.2.2, from each computed
spectrum we subtract the
Planck best-fit ΛCDM
model (Planck Collaboration XVI 2014), i.e., the theoretical
model obtained from the temperature data
fit, while
spectra are kept unaltered.
![]() |
Fig. 1 Masks and complementary selected large regions that retain fractional coverage of the sky fsky from 0.8 to 0.3 (see details in Sect. 3.3.1). The grey is the CO mask, whose complement is a selected region with fsky = 0.8. In increments of fsky = 0.1, the retained regions can be identified by the colours yellow (0.3) to black (0.8), inclusively. Also shown is the (unapodized) point source mask used. |
3.3. Selection of regions
To measure the dust polarization power spectra with high S/N, we select six large regions, the analysis regions of interest at intermediate Galactic latitude, which have effective coverage of the sky from 24 to 72% (see Sect. 3.3.1)7. For statistical studies at high Galactic latitude, we compute spectra on a complete set of smaller regions or patches (Sect. 3.3.2), similar in size to the patches observed in typical CMB experiments.
3.3.1. Large regions
For selection of all of the large regions, we used the Planck CO map from Planck Collaboration XIII (2014), smoothed to a 5° resolution, to mask the sky wherever the CO line brightness ICO ≥ 0.4 K km s-1 8. This mask is shown in Fig. 1. The complement to this mask by itself defines a preliminary region that retains a sky fraction fsky = 0.8.
We then mask the sky above successively lower thresholds of I857 in the Planck 857 GHz intensity map, smoothed to a 5° resolution, chosen such that together with the CO mask we select five more preliminary regions that retain fsky from 0.7 to 0.3 in steps of 0.1. These six regions are shown in Fig. 1.
To avoid power leakage, these six masks are then apodized by convolving with a
5° FWHM Gaussian which
alters the window function by gradually reducing the signal towards the edges of the
retained regions and thus lowers the effective retained sky coverage. The
value is simply defined as the mean sky
coverage of the window function map.
Finally, we mask data within a radius 2σbeam of point sources selected from the Planck Catalogue of Compact Sources (PCCS, Planck Collaboration XXVIII 2014) at 353 GHz.The selected sources have S/N> 7 and a flux density above 400 mJy9. Selection of spurious infrared sources in bright dust-emitting regions is avoided by using contamination indicators of infrared cirrus listed in the PCCS description (Planck Collaboration XXVIII 2014). This point-source masking is done to prevent the brightest polarized sources from producing ringing in the power spectrum estimation, while avoiding the removal of dust emitting regions and their statistical contribution to the angular power spectra. The details of this source selection will be presented in the Planck 2014 release papers. The edges of the masks around point sources were apodized with a 30′ FWHM Gaussian, further reducing the retained net effective sky coverage.
In combination these masking and apodization procedures result in six large retained
(LR) regions, which we label using the percentage of the sky retained (the net effective
fractional sky coverages, , are listed in Table 1), e.g., LR72 for the largest region and LR24 for
the smallest.
Table 1 also lists other properties of the regions, including ⟨I353⟩, the mean specific intensity at 353 GHz within the region in MJy sr-1, and NH i, the mean H i column density in units of 1020 cm-2 computed on the LAB H i survey data cube (Kalberla et al. 2005).
3.3.2. Small patches at high Galactic latitude
To examine the statistics of the angular power spectra at high Galactic latitude, we
also analyse the Planck polarization maps within patches with a size
similar to those of typical ground-based and balloon-borne CMB experiments (e.g., QUIET Collaboration et al. 2012; Hanson et al. 2013; BICEP2 Collaboration 2014b; Ade et al.
2014; Naess et al. 2014). Specifically,
we consider 400 deg2 circular areas (radius 11°̣3) centred on the central pixel
positions of the HEALPix Nside = 8 grid that have Galactic
latitude | b | >
35°. This results in 352 such patches. These are apodized with a
2° FWHM Gaussian, which
reduces the retained sky fraction to for each patch. For these small patches,
we do not mask point sources as was done in selecting the large regions, since we want
to preserve the same
for each mask, nor is the CO mask
needed, since we use these masks only on the 353 GHz high Galactic latitude data. We
note that for Nside = 8, pixels have an area of
about 54 deg2 and
characteristic centre-to-centre spacing of about 7°̣4. Therefore, this HEALPix grid
oversamples the sky relative to this patch size, so that the patches overlap and are not
independent.
4. Dust polarized angular power spectra at intermediate Galactic latitude
In this section, we quantify the 353 GHz dust polarization in the power spectrum domain,
achieving a high S/N through the use of LR regions. For convenience we present results for
and
, where
.
4.1. Description of the spectra
Using Xpol we have computed and
from the two DetSets at 353 GHz, as a
function of the multipole ℓ in the range 40–60010. These represent the first measurements of the thermal dust
and
power spectra on large fractions of the
sky for ℓ>
40 (see Ponthieu et al. 2005;
Gold et al. 2011, for earlier related studies).
In Fig. 2 we present the results for fsky = { 0.3,0.5,0.7 }. The amplitudes of the spectra increase with increasing fsky because the polarized emission is brighter on average when more sky is retained. We point out that the LR regions presented in Sect. 3.3.1 overlap because they form a nested set. However, the power spectra derived from them are almost independent of each other because each spectrum is dominated by the brightest areas in the corresponding region, namely the parts closest to the Galactic plane and hence the areas in which each LR region differs from those nested inside it.
The and
spectra are characterized by a power-law
dependence on multipole ℓ over the range ℓ = 40 to 600; furthermore, the slope is similar for
different regions that retain from 24 to 63% of the sky.
Figure 2 also shows the
power spectrum computed from the
Planck 2013 best-fit ΛCDM model of the CMB temperature data (Planck Collaboration XVI 2014), and the much lower expectation for
from the CMB model with primordial
gravitational waves with amplitude r = 0.2. At 353 GHz, the
angular power spectra of the dust are
about 3−4 orders of magnitude
larger than the CMB model at ℓ
= 30, 1−2
orders of magnitude larger at ℓ = 100, and about the same order of magnitude as
the CMB at ℓ>
300. At 353
GHz, the
angular power spectra for dust are much
greater than the CMB model power spectrum for all ℓ values in Fig. 2. The dust power spectra are larger than the
r = 0.2 CMB
spectrum by 4−5 orders of
magnitude at ℓ =
30, and by 3−4 orders of magnitude at ℓ = 100. At ℓ = 500, where the lensing
of CMB anisotropies is the dominant contribution to the CMB model spectrum, the dust is
still 2−3 orders of magnitude
higher.
![]() |
Fig. 2 Planck HFI 353 GHz |
As discussed in Appendix A, we do not expect any
significant bias, or E-to-B leakage, from the computation of the dust angular
power spectra using Xpol. Figure 2 also includes the
and
spectra at 353 GHz, computed from our
estimate of the leakage terms from intensity to polarization (discussed in Sect. 2.3), for the same three LR regions. The dust spectra
are much higher than the corresponding spectra for the leakage, which represent the main
systematic effects over the ℓ range of interest for this work. The largest
contamination of the dust signal by leakage is about 3.5% at ℓ = 50, for both the
and
spectra. Because we consider that our
estimate of the leakage maps is conservative, we conclude that contamination of the dust
and
power spectra by systematic effects
amounts to a maximum of 4% at ℓ = 50, and is less at higher multipoles.
Therefore, we have not corrected the Planck data for
intensity-to-polarization leakage in this work.
Finally, in Fig. 2 we present for the three LR
regions the absolute value of the null-test spectra anticipated in Sect. 2.3, here computed from the cross-spectra of the
HalfRing/DetSet differences, i.e., (353DS1,HR1−353DS1,HR2) /2 × (353DS2,HR1 −
353DS2,HR2)/2. These
and
spectra show a behaviour that is close to
what is expected from a white-noise dominated (thus ℓ2) spectrum.
The amplitudes of these error estimates are consistent with the noise expectations; there
is no evidence for any effects of systematics.
For completeness, in Appendix B we present a
further quantification of the power spectrum of thermal dust emission at 353 GHz via the
spectra involving temperature and polarization, and
, and the polarization cross-spectrum,
.
![]() |
Fig. 3 Best-fit power-law exponents αEE (red squares)
and αBB (blue
circles) fitted to the 353 GHz dust |
4.2. Power-law fit
To assess the apparent power-law dependence quantitatively, we made a χ2 fit to the
spectra at 353 GHz using the form
, where X ∈ { E,B
}. For
and
we fit the 22 band-powers in the range
60 <ℓ<
500. For both power spectra we restricted the fit to ℓ> 60 to avoid a
possible bias from systematic effects in the data, and also because the angular power
spectra of the dust polarization exhibit more spatial variation than is expected for a
Gaussian random field, particularly on large angular scales.
The exponents of the power-law fits, αXX, are plotted in
Fig. 3 for each of the six LR regions identified with
fsky. All exponents are consistent with
constant values of αEE = −2.41 ± 0.02
and of αBB = −2.45 ± 0.03.
While there is a slight indication of a steeper slope for
than for
, hereafter we adopt the mean exponent
−2.42 ± 0.02. This exponent
is consistent with the value αTT fitted to the
Planck 353 GHz dust intensity power spectra in this range of
ℓ (Planck Collaboration XV 2014; Planck Collaboration Int. XXII 2015) on similar-sized regions at
intermediate Galactic latitude, but slightly flatter than the αTT fitted at higher
ℓ and
higher frequency (Miville-Deschênes et al. 2007,
2010; Planck
Collaboration XV 2014).
For the fits with fixed exponent αEE,BB = −2.42, the
values of the χ2 (with number of degrees of freedom,
Nd.o.f.
= 21) are listed in Table 1 for
the six LR regions. For the spectra, the χ2 values range
from 26.3 (probability to exceed, PTE =
0.2) to 44.8 (PTE =
0.002), with a trend for a quality of fit that degrades with increasing
fsky. For the
spectra, the χ2 values range
from 14.0 (PTE = 0.87) to 22.1
(PTE = 0.39), with a trend
related to fsky.
Possible explanations for this difference between and
spectrum shape descriptions, are the
chance correlation between dust and CMB polarization, which is not taken into account in
the subtraction of the CMB
spectrum, and the increasing S/N degrading
the overall quality of the fit when going from
to
and for
from fsky = 0.3 to
fsky =
0.8 as the amplitude of the dust polarized signal increases. We stress
that while the power law in ℓ is a good general description of the shape of the
dust polarized power spectra, a full characterization would have to consider the detailed
features that can be seen in Fig. 2.
![]() |
Fig. 4 Amplitude of the dust AEE (red squares) and ABB (blue circles) power spectra, normalized with respect to the largest amplitude for each mode. These are plotted versus the mean dust intensity ⟨I353⟩ for the six LR regions (top panel). A power-law fit of the form AXX(⟨I353⟩)= KXX⟨I353⟩1.9, X ∈ { E,B }, is overplotted as a dashed line of the corresponding colour (these almost overlap). The bottom panel presents the ratio of the data and the fitted ⟨I353⟩1.9 power law; the range associated with the ± 1σ uncertainty in the power-law exponent of 1.9 is shown in grey. For details see Sect. 4.3. |
4.3. Amplitude dependence on ⟨I353⟩
Similarly to what has been done for intensity power spectra in, e.g., Gautier et al. (1992) and Miville-Deschênes et al. (2007), we investigate how the amplitude of
the dust polarization power spectrum scales with the dust intensity. To quantify the dust
emission at 353 GHz, we use
the model map derived from a modified blackbody fit to the Planck data at
ν ≥ 353 GHz
and IRAS at λ = 100
μm, presented in Planck Collaboration XI (2014). This map is corrected for zodiacal emission and
the brightest extragalactic point sources are subtracted. The Galactic reference offsets
of the underlying IRAS and Planck data were obtained through a method
based on correlation with 21 cm data from the LAB H i survey (Kalberla et al. 2005) integrated in velocity,
effectively removing the CIB monopole in the model map. The mean dust intensity,
⟨I353⟩, listed for each LR region in Table
1, ranges from 0.068 to 0.227 MJy sr-1 for increasing
. The mean column density calculated from
the LAB survey data is also listed in Table 1, with
values ranging from 1.65 × 1020
cm-2 to 6.05 ×
1020 cm-2.
For all of the LR regions, we fit the and
spectra with a power law in
ℓ, using
the fixed exponent αEE,BB = − 2.42, over
the ℓ-ranges
defined in Sect. 4.2. The amplitudes AEE derived from these
fits are listed in Table 1; the ABB amplitudes can be
retrieved from the ABB/AEE
ratio. These amplitudes are plotted as a function of ⟨I353⟩ in Fig.
4 after normalization by the maximum value found
for the largest region (LR72). We fit the empirical dependence of these amplitudes on
⟨I353⟩ as a power law of the form
AXX(⟨I353⟩)=
KXX⟨I353⟩
ϵXX
where X ∈
{E,B}. The two fitted exponents are quite similar,
ϵEE = 1.88 ± 0.02 and
ϵBB = 1.90 ± 0.02. The
exponent that we find for polarization is close to the one observed in the diffuse
interstellar medium for the dust intensity, consistent with
, where ⟨Iν⟩ is the mean value of
the dust specific intensity (Miville-Deschênes et al.
2007). Values close to 2 are expected, because we compute angular power spectra,
which deal with squared quantities.
Although the data points roughly follow this ⟨I353⟩1.9 dependence, the empirical law
fails to fully describe individual dust amplitudes (e.g., the estimate is off by about 20%
for on LR33). The scaling can help to asses
the order of magnitude of the dust polarization level on a specific region, but is not a
substitute for actually characterizing the polarized angular power spectra.
![]() |
Fig. 5 Ratio of the amplitudes of the |
4.4. Amplitude of
relative to
We examine the ratio of the amplitudes of the fitted power laws found in Sect. 4.3. The ABB/AEE
ratios are listed in Table 1, and are plotted for
different values of fsky in Fig. 5. For all of the LR regions, we observe more power in the
dust spectrum than in
. All ratios are consistent with a value of
ABB/AEE
= 0.52 ± 0.03, significantly different from unity, over various large
fractions of the intermediate latitude sky.
This result is not taken into account in existing models of polarized microwave dust
emission that have been developed to test component separation methods. For example, we
have computed the and
spectra over the LR regions for the
Planck Sky Model (Delabrouille et al.
2013) and the model of O’Dea et al.
(2012); for both models and all LR regions we find a ratio ABB/AEE
close to 1. However, these two models are based on a very simplified picture of the
Galactic magnetic field geometry and assumptions on how the polarized emission depends on
it. Further insight into the structure of the dust polarization sky is required to account
for the observed ratio.
4.5. Amplitude dependence on frequency
Finally, we explore the frequency dependence of the amplitude of the angular power
spectra. We compute the and
angular power spectra from the
Q and
U DetSet
maps at 100, 143, 217, and 353 GHz (see Sect. 3.2).
From these four sets of polarization maps, we compute ten power spectra: 100 × 100; 100 × 143; 100 × 217; 100 × 353; 143 × 143; 143 × 217; 143 × 353; 217 × 217; 217 × 353; and 353 × 353.
The ten angular cross-power spectra are consistent with a power law in ℓ, with the exponent
αEE,BB = −2.42 measured
at 353 GHz (Sect. 4.2). Therefore, to each of these spectra we fit the
amplitudes of a power-law function that has a fixed exponent αEE,BB = −
2.42, in the range 40
<ℓ< 500, for
and
. As an illustration of the quality of the
fit, for the smallest region (LR24, fsky = 0.3) the averages and dispersions
of the χ2 (21 degrees of freedom) of the fits are
and
for the ten cross-frequency spectra.
![]() |
Fig. 6 Frequency dependence of the amplitudes AEE,BB of the
angular power spectra, relative to 353 GHz (see details in Sect. 4.5). Results for
|
To compare the frequency dependence of the results of the fits to that expected from the
SED for dust polarization from Planck Collaboration Int.
XXII (2015), we converted the fitted amplitudes AEE,BB from
μK to units of (MJy sr-1)2, taking
into account the Planck colour corrections11. For all regions, we examined the frequency dependence by plotting the
amplitudes normalized to unity at 353
GHz versus the effective frequency12. A representative example is shown in Fig. 6 for the smallest region, LR24 (fsky = 0.3).
For all of the LR regions the frequency dependence found is in good agreement with the square of the adopted dust SED, which is a modified blackbody spectrum having βd = 1.59 and Td = 19.6 K (Planck Collaboration Int. XXII 2015). We note that no fit was performed at this stage and that the square of the adopted dust SED goes through the 353 GHz data point. However, if we fit the amplitude of the dust frequency dependence, with fixed βd and Td, the χ2 (Nd.o.f. = 9) is 13.1 for EE (PTE = 0.16) and 3.2 for BB (PTE = 0.96). This good agreement supports the assumption made in the present work about the faintness of the synchrotron and CO emission at high Galactic latitude (see Sect. 2.2). The residuals to the fit do not show any evidence for excess power at 100 GHz, which might arise from polarized synchrotron emission. This result is consistent with the study of synchrotron polarization at high Galactic latitudes by Fuskeland et al. (2014), and confirmed in Appendix D.4. Furthermore, we do not see any excess at either 100 or 217 GHz, such as could arise from leakage and/or polarization associated with CO line emission in these bands (Sect. 2.2.5).
5. Statistical study of the dust power spectra at high Galactic latitudes
In Sect. 4, statistical properties of dust
polarization were derived from angular power spectra computed with the
Planck353 GHz
data on large fractions of the sky ( from 24% to 72%). Most of the CMB
experiments target fields at high Galactic latitude with a smaller size than this. So now we
evaluate the dust B-mode power in such patches.
We perform a statistical analysis computing the and
power spectra of the 353 GHz Planck polarization
data for high Galactic latitude circular patches of size 400 deg2, similar to the size of the
fields observed in ground-based and balloon-borne CMB polarization experiments (see Sect.
3.3.2). We verify that the empirical scaling law of
the power spectra amplitudes as a function of the mean intensity, derived in Sect. 4.3, holds on small independent fractions of the sky that
are not nested.
5.1. Data processing
We compute the dust and
angular power spectra at 353 GHz from the
cross-spectra of the two DetSets with independent noise in this channel, using Xpol on all
of the 352 high Galactic latitude (|
b | > 35°) patches defined in Sect. 3.3.2. As shown in Appendix A, we do not expect any significant cut-sky leakage from
E
polarization to B polarization for the dust in such patches. The
and
spectra are computed in the range
40 <ℓ<
370, using top-hat binning in the intervals defined between multipoles
40, 70, 110, 160, 220, 290, and 370.
To each of the 352 and
spectra we fit the power law in
ℓ presented
in Sect. 4.2 over the range 40 <ℓ< 370 with a
slope fixed to αEE,BB = −2.42 in
order to compute the amplitudes AEE,BB and their
associated errors. Since we are interested in quantifying the possible contamination by
dust polarization for ground-based CMB experiments, here we express the amplitude of the
fitted power laws at the position of the recombination bump at ℓ = 80
(
). In Appendix C.1, using a null test, we show that the potential systematic effects
in the data do not affect the computation of the
and
spectrum amplitudes.
![]() |
Fig. 7 Fitted dust |
5.2. Results
For patches covering 1% or more of the sky, we do not observe any significant departure
in the shapes of the and
spectra from the results presented in
Sect. 4.2, even if a large dispersion of the
power-law exponent (attributable to the low S/N for some patches) is observed. The
/
ratios are on average consistent with what
has been derived on the larger LR regions in Sects. 4.4. These ratios are presented in Appendix C.2.
The fitted dust and
amplitudes (AEE and ABB) at ℓ = 80 are presented in
Fig. 7 in units of μK
at 353 GHz and as a function of the mean
dust intensity ⟨I353⟩ of each patch.
We see a clear correlation of the dust and
amplitudes with ⟨I353⟩. The data
points are consistent with the same ⟨I353⟩ dependence of the scaling law
found in Sect. 4.3, AEE,BB
∝⟨I353⟩1.9, indicating that this
empirical law also applies in the faintest regions of the sky and holds reasonably well
for the description of the amplitudes for any approximately 1% patch of the sky at
latitudes above 35°.
The amplitudes of the two empirical scaling laws are adjusted to the data points by computing the median over all the patches of AEE,BB/⟨I353⟩1.9. The amplitudes that we derive from Fig. 7 are in agreement with the amplitudes reported for the LR regions in Sect. 4.3 and Table 1.
The dust brightness ⟨I353⟩ in the cleanest 400 deg2 patches (⟨I353⟩ down to 0.038 MJy sr-1) is about a factor of 2 lower than the value listed in Table 1 for the LR24 region (0.068 MJy sr-1). Applying the empirical scaling derived in Sect. 4.3, the expected level of the dust polarization angular power spectrum is a factor of (0.068/0.038)1.9 = 3 higher on 24% of the sky than on the faintest patches.
A dispersion around this empirical power law for the patches is observed. The shaded
regions in Fig. 7 show the expected scatter in this
relation coming from the noise, together with Gaussian sample variance. The former is
computed from Planck inhomogeneous noise Monte Carlo simulations
presented in Appendix C.1 and the latter is computed
given the and the binning used to determine the
spectra (see Sect. 5.1). Thus, the general trend of
the dust
and
amplitudes among these patches follows the
empirical scaling in ⟨I353⟩, but the inhomogeneous nature of
the polarized dust emission, including variations in magnetic field orientation and grain
alignment, is responsible for the large dispersion around it, larger than expected for a
stationary Gaussian process plus instrumental noise.
Even at low ⟨I353⟩, the statistical error on the
amplitude at ℓ = 80 for such small
patches is at least 7.5 μK
(3σ).
![]() |
Fig. 8 Top: map in orthographic projection of the 150 GHz
|
5.3. Optimizing the search for a primordial B-mode signal at the recombination peak
In Sect. 5.1 the 353 GHz ABB amplitudes were
obtained for the 352 patches of 400 deg2 at Galactic latitudes | b | > 35°. These
amplitudes are expressed at ℓ
= 80, the approximate position of the maximum of the CMB tensor
B-mode
recombination peak. In order to extrapolate these amplitudes to 150 GHz, an observing
frequency typical of CMB experiments searching for B-modes, we use the dust
SED of Planck Collaboration Int. XXII (2015)
presented in Sect. 2.2.1, which is a modified
blackbody spectrum with βd = 1.59 ± 0.17 and Td = 19.6 K. The
conversion factor from the Planck 353 GHz band to the 150 GHz central
frequency for this extrapolation, taking into account the Planck colour
corrections and the statistical uncertainty of βd, is
in the SED, and this squared in
.
These extrapolated estimates are divided by the value of the r = 1 primordial tensor CMB
spectrum at ℓ = 80, 6.71 × 10-2 μK
, to express the estimated power in units
that we denote rd. Because the CMB primordial tensor
B-mode
power scales linearly with r 13, a value
rd =
0.1 would mean that the expected contamination from dust at
ℓ = 80 is
equal to the amplitude of the primordial tensor CMB
for r = 0.1. For each of these
estimates we also compute σ(rd), the quadratic
sum of the fit errors on ABB and the above
uncertainty from the extrapolation to 150 GHz. We note that the fitted amplitudes
ABB for five of these
patches are negative14, but are consistent with
rd =
0 at 1σ. The rd values vary from −0.17 to more than 10 and σ(rd) ranges from as
small as 5.6 × 10-2
to larger than 10. Taking the smallest value of σ(rd) we see that
Planck measurements of the dust
spectra for such small patches have, at
best, a statistical uncertainty of 0.17 (3 σ) in rd units.
To reveal the spatial dependence over the high latitude sky, a map of the absolute value of rd is shown in Fig. 8. Each computed value of rd is given in a pixel at the position of the centre of each of the 352 patches defined in Sect. 3.3.2 on an Nside = 8HEALPix map. The accompanying map of σ(rd) is also presented.
We can see in Fig. 8 that there is a high latitude region in the southern Galactic hemisphere for which rd is quite low. This region is also associated with a small estimated uncertainty. For example, the six HEALPix Nside = 8 pixels numbered 741, 742, 754, 755, 762, and 763 (in the HEALPix “ring” ordering scheme) have rd = 0.053 ± 0.096, 0.027 ± 0.098, −0.062 ± 0.052, −0.020 ± 0.127, 0.057 ± 0.122, and −0.031 ± 0.121, respectively. These pixels are located around Galactic coordinates l = −30°, b = −70°.
We stress that the expectation for this low level of dust contamination is valid only for
these particular patches, including their positions, sizes, shapes, and apodizations. In
addition, because we found the amplitudes of the dust
spectra associated with these patches
based on a power-law fit, our estimate does not take into account possible features in the
power spectra that might alter the precise value of dust contamination. Nevertheless,
there are clearly some patches that appear to be “optimal," i.e., cleaner than the
others. It needs to be emphasized, however, that finding the cleanest areas of the
polarized sky for primordial B-mode searches cannot be accomplished accurately
using the Planck total intensity maps alone.
6. BB angular power spectrum of dust in the BICEP2 field
In this section, we use the Planck data and the results presented above to
assess the dust polarization in the field observed by the BICEP2 experiment (BICEP2 Collaboration 2014a,b). As above, finding at 150 GHz involves two steps, measuring the
power spectrum at 353 GHz and then extrapolating the amplitude to 150 GHz.
6.1. An approximation to the observed BICEP2 field
To define a field similar to and representative of the actual BICEP2 field, we carried
out the following steps (at the HEALPix Nside = 2048 pixelization): (i) we
constructed a mask M by filling the inside of the BICEP2 deep-field
outline15 with 1 and the outside with 0; (ii) we
took the complement of this mask, M′ = 1−M; (iii) inside
M′ we
computed the distance to the border using the HEALPix “iprocess_mask” procedure; (iv) we
smoothed this “distance map” with a 7°̣5 FWHM Gaussian; (v) we apodized, computing the 7°̣5 FWHM Gaussian weight from the distance
map to obtain M′′; and (vi) we took the complement of
M′′ to be
our definition of the BICEP2 field (i.e., MB2 = 1−M′′). The
resulting field, which in the rest of this paper we will refer to as MB2, has
(689 deg2) and ⟨I353⟩= 0.060 MJy sr-1. Its
is larger than the BICEP2 373 deg2 deep-field region, but
MB2 is similar to the BICEP2 inverse noise
variance map presented in BICEP2 Collaboration
(2014a), even if it extends further in declination but less in right ascension.
In Appendix D.2 we show that the main results in
this section do not depend significantly on the definition of this field. However, we note
that the various applications of filtering in the BICEP2 pipeline mean that their observed
spatial modes are a subset of those present within their nominal observed region.
![]() |
Fig. 9 Planck 353 GHz |
6.2. Statistical estimate of the dust B-mode level in the BICEP2 field
We have seen in Sect. 5.2 that the dust
spectra amplitudes statistically follow an
empirical scaling law of the mean dust intensity ⟨I353⟩ for the
patch on which they are computed. We can use this scaling law to assess the most probable
dust
value for MB2. For this
field, which has a mean dust intensity of ⟨I353⟩ = 0.060 MJy sr-1, the expected value at
353 GHz is (13.4 ±
0.26) μK
at ℓ = 80, taking into account
the uncertainty on the fitted ⟨I353⟩1.9 amplitude from Sect. 4.3.
In order to extrapolate this value to the BICEP2 observing frequency, 150 GHz, we proceed
as in Sect. 5.3, using the same typical dust SED
from Planck Collaboration Int. XXII (2015). For
patches the size of MB2 the expected dispersion of
βd (Sect. 2.2.1) is 0.11, which introduces an uncertainty in the extrapolation. The
extrapolation, unlike in that section, is made to 150 GHz taking into account the BICEP2
bandpass16. The conversion factor from the
Planck 353 GHz band to the 150 GHz BICEP2 band, taking into account
both Planck and BICEP2 colour corrections, as well as the statistical
error on βd, is thus
in the SED, and therefore this value
squared in
.
The resulting expected dust amplitude at 150 GHz is
μK
at ℓ = 80. The uncertainty
quoted here applies to the most probable value for patches with a similar ⟨I353⟩, but it
ignores the intrinsic dispersion around this value (see Sect. 5.2). Therefore, a direct measurement of the dust
spectrum in MB2 is required
to complete the assessment of its polarized dust level.
In Sect. 5.3, we have localized the
amplitudes from the analysis of patches on
a map. The pixels falling inside the approximate BICEP2 deep-field region shown in Fig.
8 give a mean value of rd = 0.207,
i.e., an expected dust power of
= 1.39 × 10-2 μK
at 150 GHz and ℓ = 80; we note that in
this region the Planck estimate of the dust contamination is
significantly higher, by a factor of about 2, than it is in some more optimal patches.
Nevertheless, a more accurate estimate of the dust
polarization amplitude in the BICEP2
observed field requires the direct computation of the angular power spectrum on a similar
region.
6.3. BB angular power spectrum of dust in the BICEP2 field
We now compute the spectrum of the dust at 150 GHz on
MB2 presented in Sect. 6.1, using the Planck 353 GHz data. For
this purpose, we compute the 353 GHz
spectrum by cross-correlating the two
DetSets that have independent noise at this frequency in the three ℓ-bins defined by the
intervals between multipoles 40, 120, 250, and 400. The ℓ-bin sizes were increased
with respect to the previous sections in order to increase the S/N, especially around the
range of the CMB recombination bump. Appendix D.1
confirms that the result does not depend on the method of computing the power spectrum.
This power spectrum is extrapolated to 150 GHz as in Sect. 6.2, with an extrapolation uncertainty estimated from the inferred dispersion of
βd. Our final estimate of the
spectrum is presented in Fig. 9, together with its 1σ error budget. For the
first bin, ℓ =
40−120, the expected level of dust polarized
, as extrapolated to 150 GHz, is
1.32 ×
10-2 μK
(Fig. 9). The statistical error, estimated from Monte Carlo simulations of
inhomogeneous Planck noise (presented in Appendix A for this particular binning), is ± 0.29 × 10-2 μK
, so that the dust
spectrum is statistically detected at
4.5σ in
this broad ℓ
bin.
In order to assess the potential contribution from systematics, we have computed the dust
spectrum on MB2 on different
subsets of the data and performed null tests, which are presented in Appendix D.3. In this lowest bin of ℓ, we do not observe any
departure from what is allowed by noise. Nevertheless, we stress that below the noise
level our cross-spectra could be subject to a positive or negative bias due to systematic
effects. For example, if instead of taking the DetSets cross-spectra (as we have done
throughout this paper) we take the mean value computed from the DetSets, HalfRings, and
Years cross-spectra (presented in Appendix D.3), the
statistical significance of our measurement is decreased from 4.5σ to 3.6σ.
The uncertainty coming from the MB2 definition (presented in Appendix
D.2) is 0.04 × 10-2 μK for this bin, thus much less than the
statistical error. For this reason, it is not added to the error budget. However, the
spectral extrapolation to 150 GHz adds an additional uncertainty (+ 0.28,−0.24) ×
10-2 μK
to the estimated power in MB2, added in
quadrature in Fig. 9.
![]() |
Fig. 10 Frequency dependence of the amplitude ABB of the
angular power spectrum |
The expected value in this lowest-ℓ bin from direct computation of the
power spectrum on MB2, as shown in
Fig. 9, is lower than (but consistent with) the
statistical expectation from the analysis of the 352 high Galactic latitude patches
presented in Sects. 5.2 and 6.2. This indicates that MB2 is not one of the outliers of Fig.
7 and therefore its dust B-mode power is well
represented by its mean dust intensity through the empirical scaling law
⟨I353⟩1.9.
These values of the amplitude in the ℓ range of the primordial
recombination bump are of the same magnitude as those reported by BICEP2 Collaboration (2014b). Our results emphasize the need for a
dedicated joint analysis of the B-mode polarization in this region incorporating
all pertinent observational details of the Planck and BICEP2 data sets,
which is in progress.
6.4. Frequency dependence
We complement the power spectrum analysis of the 353 GHz map with Planck data at lower frequencies. As in the analysis in Sect. 4.5, we compute the frequency dependence of the BB power measured by Planck at HFI frequencies in the BICEP2 field, using the patch MB2 as defined in Sect. 6.1.
We compute on MB2 the
Planck auto- and cross-power spectra from the
three Planck HFI bands at 100, 143, 217, and 353 GHz, using the two
DetSets with independent noise at each frequency, resulting in ten angular power spectra
(100 × 100, 100 × 143, 100 × 217, 100 × 353, 143 × 143, 143 × 217, 143 × 353, 217 × 217, 217 × 353, and 353 × 353), constructed by combining the
cross-spectra as presented in Sect. 3.2. We use the
same multipole binning as in Sect. 6.3. To each of
these
spectra, we fit the amplitude of a power
law in ℓ with
a fixed exponent αBB = −0.42 (see Sect. 4.2). In Fig. 10
we plot these amplitudes as a function of the effective frequency from 143 to 353 GHz, in
units of sky brightness squared, as in Sect. 4.5.
Data points at effective frequencies below 143 GHz are not presented, because the dust
polarization is not detected at these frequencies. An upper limit on the synchrotron
contribution at 150 GHz from the Planck LFI data is given in Appendix
D.4.
We can see that the frequency dependence of the amplitudes of the Planck
HFI spectra is in very good agreement with a
squared dust modified blackbody spectrum having βd = 1.59 and Td = 19.6 K
(Planck Collaboration Int. XXII 2015). We note
that this emission model was normalized only to the 353 GHz point and that no global fit
has been performed. Nevertheless, the χ2 value from the amplitudes relative to
this model is 4.56 (Nd.o.f.
= 7). This shows that dust dominates in the specific MB2 region
defined where these cross-spectra have been computed. This result emphasizes the need for
a dedicated joint Planck-BICEP2 analysis.
7. Conclusions
We have presented the first nearly all-sky statistical analysis of the polarized emission
from interstellar dust, focussing mostly on the characterization of this emission as a
foreground contaminant at frequencies above 100 GHz. Our quantitative analysis of the
angular dependence of the dust polarization relies on measurements at 353 GHz of the
and
(alternatively
and
) angular power spectra for multipoles
40 <ℓ<
500. At this frequency only two polarized components are present: dust
emission; and the CMB, which is subdominant in this multipole range. We have found that the
statistical, spatial, and spectral distribution properties can be represented accurately by
a simple model over most of the sky, and for all frequencies at which
Planck HFI measures polarization.
-
The angular power spectra
and
at 353 GHz are well fit by power laws in ℓ with exponents consistent with αEE,BB = −2.42 ± 0.02, for sky fractions ranging from 24% to 72% for the LR regions used.
-
The amplitudes of
and
in the LR regions vary with mean dust intensity at 353 GHz, ⟨I353⟩, roughly as ⟨I353⟩1.9.
-
The frequency dependence of the dust
and
from 353 GHz down to 100 GHz, obtained after removal of the
prediction from the Planck best-fit CMB model (Planck Collaboration XVI 2014), is accurately described by the modified blackbody dust emission law derived in Planck Collaboration Int. XXII (2015), with βd = 1.59 and Td=19.6 K.
-
The ratio between the amplitudes of the two polarization power spectra is
, which is not consistent with the simplest theoretical models.
-
Dust
and
spectra computed for 352 high Galactic latitude 400 deg2 patches satisfy the above general properties at 353 GHz and have the same frequency dependence.
We have shown that the Planck determination of the 353 GHz dust polarization properties is unaffected by systematic errors for ℓ> 40. This enables us to draw the following conclusions relevant for CMB polarization experiments aimed at the detection of primordial CMB tensor B-modes.
-
Extrapolating the Planck 353 GHz
spectra computed on the 400 deg2 circular patches at high Galactic latitude to 150 GHz shows that we expect significant contamination by dust over most of the high Galactic latitude sky in the ℓ range of interest for detecting a primordial
spectrum.
-
Even for the cleanest of these regions, the Planck statistical error on the estimate of
amplitude at ℓ = 80 for such small regions is at best 0.17 (3σ) in units of rd.
-
Our results show that subtraction of polarized dust emission will be essential for detecting primordial B-modes at a level of r = 0.1 or below.
-
There is a significant dispersion of the polarization
amplitude for a given dust total intensity. Choices of the cleanest areas of the polarized sky cannot be made accurately using the Planck total intensity maps alone.
-
At present, component separation, or template cleaning, can best be done with the Planck HFI 353 GHz data, but the accuracy of such cleaning is limited by Planck noise in small fields. Ground-based or balloon-borne experiments should include dust channels at high frequency. Alternatively, if they intend to rely on the Planck data to remove the dust emission, they should optimize the integration time and area so as to have a similar S/N level for the CMB and dust power spectra.
-
Over the multipole range 40 <ℓ< 120, the Planck 353 GHz
power spectrum extrapolated to 150 GHz yields a value 1.32 × 10-2 μK
, with statistical error ± 0.29 × 10-2 μK
and a further uncertainty (+0.28, −0.24) × 10-2 μK
from the extrapolation. This value is comparable in magnitude to the BICEP2 measurements at these multipoles that correspond to the recombination bump.
-
The frequency dependence of
across the Planck bands is consistent with the typical SED of dust polarization (Planck Collaboration Int. XXII 2015).
-
Assessing the dust contribution to the B-mode power measured by the BICEP2 experiment requires a dedicated joint analysis with Planck, incorporating all pertinent observational details of the two data sets, such as masking, filtering, and colour corrections.
-
We have identified regions in which the dust polarization
amplitude may be significantly lower, by about a factor of 2, than in the BICEP2 observing region.
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 (in particular the lead countries France and Italy), with contributions from NASA (USA) and telescope reflectors provided by a collaboration between ESA and a scientific consortium led and funded by Denmark.
While this paper was in preparation three papers have used publicly available polarization information from Planck to infer potentially high levels of dust contamination in the BICEP2 field (Flauger et al. 2014; Mortonson & Seljak 2014; Colley & Gott 2015).
The conversion factor from KCMB to MJy sr-1 is computed for a standard specific intensity Iν ∝ ν−1. The colour correction depends on the actual dust SED; it is the scaling factor used to transform from the specific intensity of the dust emission, at the reference frequency, to the Planck brightness in MJy sr-1 (see Eq. (19) in Planck Collaboration Int. XXII 2015). The conversion factors and the colour corrections are computed via Eq. (32) in Planck Collaboration IX (2014) using the Planck HFI filters and the PlanckUcCC software available through the Planck Explanatory Supplement (http://www.sciops.esa.int/wikiSI/planckpla/index.php?title=Unit_conversion_and_Color_correction&instance=Planck_Public_PLA); we use the band-average values.
More specifically, for the latitude range 10°< | b | < 60°, with patches contained within the region in Fig. 1 (below) defined by including that with fsky = 0.8 and then removing that with fsky = 0.4.
Although the selection process is similar to that in Planck Collaboration XV (2014), there are differences in detail.
Conversion factors were computed as described in Sect. 2.1, here using colour corrections corresponding to a dust modified blackbody spectrum with βd = 1.59 and Td = 19.6 K.
http://bicepkeck.org This outline encloses the effective 373 deg2 of deep integration (87 nK-degrees, BICEP2 Collaboration 2014b).
Acknowledgments
The development of Planck has been supported by: ESA; CNES and CNRS/INSU-IN2P3-INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MICINN, 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); and PRACE (EU). A description of the Planck Collaboration and a list of its members, including the technical or scientific activities in which they have been involved, can be found at http://www.rssd.esa.int/index.php?project=PLANCK&page=Planck_Collaboration. Some of the results in this paper have been derived using the HEALPix package. 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.
References
- Ade, P. A. R., Akiba, Y., Anthony, A. E., et al. 2014, Phys. Rev. Lett., 113, 021301 [NASA ADS] [CrossRef] [Google Scholar]
- Arnold, K., Ade, P. A. R., Anthony, A. E., et al. 2010, in SPIE Conf. Ser., 7741 [Google Scholar]
- Austermann, J. E., Aird, K. A., Beall, J. A., et al. 2012, in SPIE Conf. Ser., 8452 [Google Scholar]
- Baccigalupi, C. 2003, New Astron. Rev., 47, 1127 [Google Scholar]
- Barkats, D., Aikin, R., Bischoff, C., et al. 2014, ApJ, 783, 67 [NASA ADS] [CrossRef] [Google Scholar]
- Battye, R. A., Browne, I. W. A., Peel, M. W., Jackson, N. J., & Dickinson, C. 2011, MNRAS, 413, 132 [NASA ADS] [CrossRef] [Google Scholar]
- Benoît, A., Ade, P., Amblard, A., et al. 2004, A&A, 424, 571 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- BICEP2 Collaboration 2014a, ApJ, 792, 62 [NASA ADS] [CrossRef] [Google Scholar]
- BICEP2 Collaboration 2014b, Phys. Rev. Lett., 112, 241101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Bunn, E. F., Zaldarriaga, M., Tegmark, M., & de Oliveira-Costa, A. 2003, Phys. Rev. D, 67, 023501 [NASA ADS] [CrossRef] [Google Scholar]
- Colley, W. N., & Gott, III, J. R. 2015, MNRAS, 447, 2034 [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]
- Draine, B. T. 2004, in The Cold Universe, Saas-Fee Advanced Course 32 (Springer), eds. A. W. Blain, F. Combes, B. T. Draine, D. Pfenniger, & Y. Revaz, 213 [Google Scholar]
- Dunkley, J., Amblard, A., Baccigalupi, C., et al. 2009a, in AIP Conf. Ser. 1141, eds. S. Dodelson, D. Baumann, A. Cooray, et al., 222 [Google Scholar]
- Dunkley, J., Komatsu, E., Nolta, M. R., et al. 2009b, ApJS, 180, 306 [NASA ADS] [CrossRef] [Google Scholar]
- Fauvet, L., Macías-Pérez, J. F., Aumont, J., et al. 2011, A&A, 526, A145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fauvet, L., Macías-Pérez, J. F., & Désert, F. X. 2012, Astropart. Phys., 36, 57 [NASA ADS] [CrossRef] [Google Scholar]
- Flauger, R., Hill, J. C., & Spergel, D. N. 2014, J. Cosmol. Astropart. Phys., 8, 39 [Google Scholar]
- Fraisse, A. A., Brown, J.-A. C., Dobler, G., et al. 2009, in AIP Conf. Ser. 1141, eds. S. Dodelson, D. Baumann, A. Cooray, et al., 265 [Google Scholar]
- Fraisse, A. A., Ade, P. A. R., Amiri, M., et al. 2013, J. Cosmol. Astropart. Phys., 4, 47 [Google Scholar]
- Fuskeland, U., Wehus, I. K., Eriksen, H. K., & Næss, S. K. 2014, ApJ, 790, 104 [NASA ADS] [CrossRef] [Google Scholar]
- Gautier, III, T. N., Boulanger, F., Perault, M., & Puget, J. L. 1992, AJ, 103, 1313 [NASA ADS] [CrossRef] [Google Scholar]
- Ghribi, A., Aumont, J., Battistelli, E. S., et al. 2014, J. Low Temp. Phys., 176, 698 [NASA ADS] [Google Scholar]
- Gold, B., Odegard, N., Weiland, J. L., et al. 2011, ApJS, 192, 15 [Google Scholar]
- Goldreich, P., & Kylafis, N. D. 1982, ApJ, 253, 606 [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]
- Grain, J., Tristram, M., & Stompor, R. 2009, Phys. Rev. D, 79, 123515 [Google Scholar]
- Grainger, W., Aboobaker, A. M., Ade, P., et al. 2008, ASP Conf. Ser., 7020 [Google Scholar]
- Hanson, D., Hoover, S., Crites, A., et al. 2013, Phys. Rev. Lett., 111, 141301 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Hartmann, D., Magnani, L., & Thaddeus, P. 1998, ApJ, 492, 205 [NASA ADS] [CrossRef] [Google Scholar]
- Hildebrand, R. H. 1988, Roy. Astron. Quant. J., 29, 327 [Google Scholar]
- Hildebrand, R. H., Dotson, J. L., Dowell, C. D., Schleuning, D. A., & Vaillancourt, J. E. 1999, ApJ, 516, 834 [NASA ADS] [CrossRef] [Google Scholar]
- Hivon, E., Górski, K. M., Netterfield, C. B., et al. 2002, ApJ, 567, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Jaffe, T. R., Ferrière, K. M., Banday, A. J., et al. 2013, MNRAS, 431, 683 [NASA ADS] [CrossRef] [Google Scholar]
- Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lamarre, J., Puget, J., Ade, P. A. R., et al. 2010, A&A, 520, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Li, H.-B., & Henning, T. 2011, Nature, 479, 499 [NASA ADS] [CrossRef] [Google Scholar]
- Macellari, N., Pierpaoli, E., Dickinson, C., & Vaillancourt, J. E. 2011, MNRAS, 418, 888 [NASA ADS] [CrossRef] [Google Scholar]
- Magnani, L., Hartmann, D., Holcomb, S. L., Smith, L. E., & Thaddeus, P. 2000, ApJ, 535, 167 [NASA ADS] [CrossRef] [Google Scholar]
- Martin, P. G. 2007, in EAS Pub. Ser. 23, eds. M.-A. Miville-Deschênes, & F. Boulanger, 165 [Google Scholar]
- Massardi, M., Burke-Spolaor, S. G., Murphy, T., et al. 2013, MNRAS, 436, 2915 [NASA ADS] [CrossRef] [Google Scholar]
- Matthews, T. G., Ade, P. A. R., Angilè, F. E., et al. 2014, ApJ, 784, 116 [NASA ADS] [CrossRef] [Google Scholar]
- Mennella, A., Butler, R. C., Curto, A., et al. 2011, A&A, 536, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Miville-Deschênes, M.-A., Lagache, G., Boulanger, F., & Puget, J.-L. 2007, A&A, 469, 595 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Miville-Deschênes, M.-A., Martin, P. G., Abergel, A., et al. 2010, A&A, 518, L104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mortonson, M. J., & Seljak, U. 2014, JCAP, 10, 035 [NASA ADS] [CrossRef] [Google Scholar]
- Naess, S., Hasselfield, M., McMahon, J., et al. 2014, J. Cosmol. Astropart. Phys., 10, 7 [Google Scholar]
- Niemack, M. D., Ade, P. A. R., Aguirre, J., et al. 2010, in SPIE Conf. Ser., 7741 [Google Scholar]
- O’Dea, D. T., Clark, C. N., Contaldi, C. R., & MacTavish, C. J. 2012, MNRAS, 419, 1795 [NASA ADS] [CrossRef] [Google Scholar]
- Planck Collaboration I. 2014, A&A, 571, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration II. 2014, A&A, 571, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration VI. 2014, A&A, 571, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration VII. 2014, A&A, 571, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration VIII. 2014, A&A, 571, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration IX. 2014, A&A, 571, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration X. 2014, A&A, 571, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration XI. 2014, A&A, 571, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration XIII. 2014, A&A, 571, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration XV. 2014, A&A, 571, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration XVI. 2014, A&A, 571, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration XVII. 2014, A&A, 571, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration XXVIII. 2014, A&A, 571, A28 [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. XXI. 2015, A&A, 576, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration Int. XXII. 2015, A&A, 576, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck HFI Core Team 2011, A&A, 536, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ponthieu, N., Macías-Pérez, J. F., Tristram, M., et al. 2005, A&A, 444, 327 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- QUIET Collaboration,Araujo, D., Bischoff, C., et al. 2012, ApJ, 760, 145 [NASA ADS] [CrossRef] [Google Scholar]
- Rosset, C., Tristram, M., Ponthieu, N., et al. 2010, A&A, 520, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rubiño-Martín, J. A., Rebolo, R., Tucci, M., et al. 2010, in Highlights of Spanish Astrophysics V, ed. J. M. Diego, L. J. Goicoechea, J. I. González-Serrano, & J. Gorgas, 127 [Google Scholar]
- Rybicki, G. B., & Lightman, A. P. 1979, Radiative processes in Astrophysics (Wiley-Interscience) [Google Scholar]
- Seiffert, M., Borys, C., Scott, D., & Halpern, M. 2007, MNRAS, 374, 409 [NASA ADS] [CrossRef] [Google Scholar]
- Smith, K. M. 2006, New Astron Rev., 50, 1025 [NASA ADS] [CrossRef] [Google Scholar]
- Staniszewski, Z., Aikin, R. W., Amiri, M., et al. 2012, J. Low Temp. Phys., 167, 827 [NASA ADS] [CrossRef] [Google Scholar]
- Tassis, K., & Pavlidou, V. 2014, MNRAS, submitted [arXiv:1410.8136] [Google Scholar]
- Tauber, J. A., Mandolesi, N., Puget, J.-L., et al. 2010, A&A, 520, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Tristram, M., Macías-Pérez, J. F., Renault, C., & Santos, D. 2005, MNRAS, 358, 833 [NASA ADS] [CrossRef] [Google Scholar]
- Tucci, M., & Toffolatti, L. 2012, Adv. Astron., 2012, 624987 [Google Scholar]
- Tucci, M., Martínez-González, E., Vielva, P., & Delabrouille, J. 2005, MNRAS, 360, 935 [NASA ADS] [CrossRef] [Google Scholar]
- Vaillancourt, J. E. 2007, in EAS Pub. Ser., 23, eds. M.-A. Miville-Deschênes, & F. Boulanger, 147 [Google Scholar]
Appendix A: Power spectrum estimator performance
Using Xpol (see Sect. 3.1.1) we have estimated the
dust polarization and
angular power spectrum, not on the full
sky, but on particular cuts of the sky, such as the LR24 region defined in Sect. 3.3.1 and the much smaller MB2 patch
defined in Sect. 6.1. To validate the performance of
Xpol on such cuts we compare the results of this algorithm to those from Xpure (see Sect.
3.1.2), using simulated data in both cases.
![]() |
Fig. A.1 Upper parts of both panels: performance of the
Xpol (blue points, see Sect. 3.1.1) and Xpure
(red squares, see Sect. 3.1.2) algorithms on
the |
Simulating the dust polarization for assessment of the performance of polarization angular power spectra estimators is an important issue. One either relies on statistically isotropic, Gaussian simulations with polarization power spectra similar to those of the dust, or alternatively one uses non-Gaussian, anisotropic simulations. The problem in the first case is that the hypotheses of Gaussianity and statistical isotropy might apply only partially to a dust polarization map. The problem in the second case is that one has to produce a non-Gaussian and anisotropic simulation, which has a defined input angular power spectrum, in order to be able to characterize what is retrieved in the output. In this section, we will use statistically isotropic, Gaussian simulations to explore the performance of our algorithms for recovering the dust polarization angular power spectra on a given masked sky. This choice is motivated by the fact that cosmologists have tended to make this assumption when analysing the CMB anisotropies from a sky contaminated by non-Gaussian and anisotropic processes.
We generated 1000 full sky Q and U Gaussian map simulations from
and
power spectra having
, the typical dust power spectrum shape
that we measured in Sect. 4.2. Additionally, we gave
to these spectra the B/E amplitude hierarchy that we have measured in
Sect. 4.4, i.e.,
. For each of these simulated dust
polarized signal Q and U maps, we computed the
and
angular power spectra on the cut sky LR24
region and the much smaller MB2 patch. We note that the structure in
the simulations is driven by the assumed power spectra and the random phases drawn for
each mode and is not correlated with the geometry of the mask.
The results are shown in Fig. A.1. For each region,
we give the angular power spectra recovered by Xpol
and Xpure from the simulations, together with the input power spectrum and the sample
variance associated with each region.
We can see that for both methods there is no overall bias in the recovered
angular power spectra, demonstrating that
there is no substantial leakage from E to B. At multipoles above 40, the recovered values
agree with the input values and the maximum excursions (5%) seen in some bins are
consistent with the expected accuracy for this number of simulations. Moreover, Xpol and
Xpure give very similar results in terms of the mean value and uncertainty.In the
multipole bin ℓ ≃
40 on MB2 the small offsets for both Xpol and
Xpure (−4% and 8%, respectively) might be hinting at a
limitation of the estimation procedure in this multipole bin because of the size of the
patch.
A related comparative assessment of the two methods on actual Planck data is presented in Appendix D.1.
Appendix B: Complementary dust spectra at 353 GHz:
,
, and
Using the same procedure as described in Sects. 3
and 4 for the dust polarized angular power spectra
for the LR regions at intermediate Galactic latitude, we have computed the 353 GHz dust
spectra involving temperature and polarization, and
, and the cross-spectrum for polarization,
. This completes the entire set of
polarization-related power spectra. Spectra such as
and
necessarily vanish for cosmic fields that
satisfy statistical isotropy and parity invariance. However, the emission from our Galaxy
satisfies neither of these and so the
and
spectra are generally expected to be
non-zero.
![]() |
Fig. B.1 Planck HFI 353 GHz |
In Fig. B.1 we present the
results for fsky = { 0.3,0.5,0.7
}; we note that the Planck best-fit CMB
power spectrum was removed, as in Sect.
3.2. Again, the amplitudes of the spectra increase
with increasing fsky because the polarized emission is
brighter on average when more sky is retained. The amplitudes of the spectra are about a
factor of 3 higher than those reported for
in Fig. 2. For the largest region there is evidence of a power-law dependence on
multipole ℓ
over the range ℓ =
40 to 600,
with a slope compatible with that found for
and
in Sect. 4.2 (αTE = −2.37, −2.43,
and −2.47 for fsky = 0.3, 0.5,
and 0.7, respectively). The fact that the slopes of the spectra are similar to those found
in Sect. 4.2 for the dust polarization and in Planck Collaboration Int. XXII (2015) for dust
intensity, indicates that the TE cross-correlation is dominated by the dust
correlations in this multipole range, and not by the other components in the intensity
maps (e.g., CIB or point sources). We also show the
power spectrum computed from the
Planck 2013 best-fit ΛCDM model of the CMB temperature data (Planck Collaboration XVI 2014). As in Fig. 2 for
, the CMB model is below the 353 GHz
angular power spectra at low ℓ and about the same order of magnitude at
ℓ> 300.
![]() |
Fig. B.2 Planck HFI 353 GHz |
Similarly, the results for are shown in Fig. B.2, but now on a linear scale. For the two largest regions
(fsky = 0.5 and 0.7), we detect a
significant signal in the range ℓ = 60−130 that increases with the sky area. In the
same multipole range, the dust is only marginally detected for the fsky = 0.3
region. At lower and higher ℓ, the individual binned measurements of
on all the regions are compatible with
zero.
![]() |
Fig. B.3 Planck HFI 353 GHz |
Finally, we show the results for in Fig. B.3. There is a marginal detection of the
spectrum in the range ℓ = 60−130 for the
fsky = 0.7 mask. For the other masks and
at other multipoles, the
spectra are compatible with zero.
Appendix C: Notes on the analysis of the high Galactic latitude patches
Appendix C.1: Assessment of Planck polarization systematic uncertainties
![]() |
Fig. C.1 Top: same as Fig. 7, but
instead of computing the |
![]() |
Fig. C.2 Ratio of the dust ABB and AEE angular power spectra amplitudes computed on the 352 patches at high Galactic latitudes (described in Sect. 5, circles), as a function of the mean dust intensity ⟨I353⟩ of the patch. The mean value, ⟨ABB/AEE⟩ = 0.51 ± 0.18 here, is represented as a dashed-blue line and the associated ± 1σ dispersion as a light blue shaded area. The ABB/AEE ratio for the MB2 spectra presented in Sect. 6.3 is shown as an orange diamond. |
To assess the level of systematic effects in the results presented in Sect. 5, we use a null-test analysis (see also Appendix
D.3). We repeated the statistical analysis for
the 353 GHz and
cross-spectra on the high-Galactic
latitude patches. The steps of the analysis are exactly the same as in Sect. 5.1, except that instead of using the 353 GHz
353DS1 ×
353DS2 cross-power spectra, we use the (353DS1,HR1−353DS1,HR2)/2 ×
(353DS2,HR1−353DS2,HR2)/2 half-difference
cross-spectra. These cross-spectra are expected to be consistent with zero for
uncorrelated noise and subdominant systematics.
The computed dust and
amplitudes AEE,BB (at
ℓ = 80
computed on each of these patches) for the null-test spectra, in μK
at 353 GHz, are shown in Fig. C.1 as a function of ⟨I353⟩. We
also report in this figure the expected noise level computed from Monte Carlo
simulations of Planck inhomogeneous white noise.
We can see that the fitted amplitudes for the null-test spectra on the 352 high Galactic latitude patches are scattered around zero, without any obvious bias, for the entire range of ⟨I353⟩; this is the expectation for half-difference cross-spectra. Taking into account the individual error bars, there is no patch that has a null-test spectrum amplitude inconsistent with the ± 3σ range of the instrumental noise dispersion.
We thus conclude from this analysis that the results presented in Sect. 5 are not significantly affected by instrumental systematics.
Appendix C.2: Dust B-mode amplitudes compared to E-mode
In Sect. 5, we computed the 353 GHz dust
and
amplitudes on the 352 high Galactic
latitude patches defined in Sect. 3.3.2. We
present in Fig. C.2 the ratio of the amplitudes of
these spectra, ABB/AEE
as a function of the mean dust intensity ⟨I353⟩ of the patch on which it was
computed.
These ratios are scattered around the inverse noise variance weighted ⟨ABB/AEE⟩ = 0.52 ± 0.19 value. This value is consistent with the one computed in Sect. 4.4 on the LR regions. Nevertheless, it can be seen from Fig. C.2 that there are patches with a significantly different ABB/AEE ratio, e.g., some smaller values for ⟨I353⟩ in the range 0.08−0.1 MJy sr-1.
We additionally report in Fig. C.2 the
ABB/AEE
value determined for the MB2 spectra computed in Sect. 6.3. The AEE and
ABB amplitudes are
computed as they are in Sect. 5, fitting a power
law in ℓ to
the and
spectra with a fixed exponent
αEE,BB = −2.42.
Thus, for MB2, we find that under these
assumptions ABB/AEE
= 0.70 ± 0.19. This value is higher than the average for such fields,
but consistent with the mean expected value.
Appendix D: Systematic effects relating to the power spectrum estimate in the BICEP2 field
Appendix D.1: Dust
spectrum measurements in the BICEP2
field using Xpol and Xpure
We present in Fig. D.1 the difference of the
power spectra computed by Xpol and Xpure
in MB2
(
). The processing is the same as
described in Sect. 6.3.
The Xpol and Xpure angular power spectra are consistent with each other in all three ℓ bins. The maximum difference observed is in the ℓ = 40–120 bin, where the Xpure spectrum is 10% higher than the Xpol spectrum presented in Sect. 6.3. Moreover, the Xpol recovered error bars and the Xpure error bars from Monte Carlo simulations of Planck inhomogeneous noise, also presented in Fig. D.1, are consistent.
We conclude from this comparison that the results presented in Sect. 6.3 do not depend significantly on the method used for
estimating the dust spectrum.
![]() |
Fig. D.1 Xpol minus XpurePlanck 353 GHz
|
![]() |
Fig. D.2 Planck 353 GHz |
Appendix D.2: Impact of uncertainties in the specification of the BICEP2 field
Our estimation of the dust angular power spectrum in
MB2 can depend on the detailed
definition of the field. To explore the dependence of our results on the extent and
specification of our approximation of the BICEP2 field, we repeat our computation of the
dust
power spectrum on the 353 GHz data
extrapolated to 150 GHz, as in Sect. 6.3, for
slight variations in the definition of this field.
For this purpose, we start from the mask M defined in Sect. 6.1, which consists of zeros and ones. We increase its size (with a smoothing
and a cut) in order to have an outline that extends exactly 1° further in every direction. This
defines (“b” for binary, since it contains zeros
and ones like M). We then decrease the size of
(with smoothing and cut) by exactly
1° in every direction to
obtain
. After that
is shrunk in the same way to obtain
and we proceed iteratively down to
. Then to each
mask we apply exactly the same procedure
for apodizing that we applied to M to produce MB2 in Sect.
6.1 to yield the final Mi field. The effective
area (
) is 805 deg2 for M5,
677 deg2 for
M4, 564 deg2 for M3,
460 deg2 for
M2, and 365 deg2 for M1. We note
that because of the forward and backward smoothing, MB2 defined in
Sect. 6.1 is similar, but not exactly equal to
M4.
We compute the 353 GHz power spectrum on all five fields using
Xpol and extrapolate the results to 150 GHz with the frequency dependence from Planck Collaboration Int. XXII (2015), as in Sect.
6.3. We show in Fig. D.2 the difference between these spectra and the
spectrum computed on MB2 presented
in Sect. 6.3. We also show the error budget for
MB2, centred on zero because we are
looking at differences between power spectra.
The recovered scatter is very low in the bin with ℓ = 40–120 (the standard
deviation of the mean value among the fields is 0.04 × 10-2
μK), much lower than the statistical
uncertainty in this particular bin. In the two other bins, the scatter is more
important, but still consistent with the statistical uncertainty.
We conclude from this analysis that the main result presented in Sect. 6.3 does not depend on our precise definition of MB2.
![]() |
Fig. D.3 Planck 353 GHz |
Appendix D.3: Assessment of the Planck HFI systematics in the BICEP2 field
To assess the potential contribution of systematic effects to the dust 353 GHz
spectrum computed on MB2, we repeat
the analysis presented in Sect. 6.3 on null tests
coming from subtracting different Planck 353 GHz data subsets.
The cross-DetSets, cross-HalfRings, and cross-Years spectra are computed for the
353 GHz data extrapolated to 150 GHz. The differences of these cross-spectra (DetSets
minus Years, Years minus HalfRings, and HalfRings minus DetSets) are presented in Fig.
D.3. In the first ℓ bin (ℓ = 40−120), the three
differences of the dust power spectra are consistent with zero,
given the statistical uncertainty. They exhibit a dispersion consistent with what is
expected from statistical noise. In the other bins the dispersion is larger, but still
consistent with the noise. We conclude that the systematic errors are subdominant.
In addition to these data cuts, we create null-test
spectra on MB2 from the
four data subsets at 353 GHz involving detector sets and ring halves. We compute the
cross-DetSet spectra between their two HalfRing half-differences (353DS1, HR1−353DS1, HR2)/2 ×
(353DS2, HR1−353DS2, HR2)/2 and the
cross-HalfRings between the two DetSet half-differences (353DS1, HR1−353DS2, HR1)/2 ×
(353DS1, HR2−353DS2, HR2)/2. These null-test
spectra are also shown in Fig. D.3. We can see
that these spectra are consistent with zero in all the ℓ bins, given the
statistical and systematic uncertainties.
We list in Table D.1 the χ2 values for the null hypothesis, as well as the PTE for these null tests (plus additional values not presented in Fig. D.3). All these tests are consistent with the null hypothesis.
Appendix D.4: Synchrotron contribution
An estimate can be made of the level of a synchrotron contribution in MB2 by a
simple extrapolation from the lowest frequency Planck LFI channels
using a reasonable range of spectral indices. Starting with the Planck
LFI 28.4 GHz Q and U maps masked to retain MB2, a pure
pseudo-Cℓ power spectrum
estimator was used to calculate the B-mode band-powers
. The noise bias and significance was
estimated from 100 FFP8 noise simulations17, which
show the 28.4 GHz power spectrum to be completely within the expected ± 2σ noise envelope. We
therefore use the 3σ upper limit on the scatter of the simulations as
the maximum allowable synchrotron signal at 28.4 GHz.
A spectral index of −3.12 ± 0.04 for synchrotron emission at high Galactic latitudes was found by Fuskeland et al. (2014) using TT plots of WMAP 9 yr data. Adopting a spectral index of −3.1 for a direct extrapolation to 150 GHz, we find 3σ upper limits on the possible contribution of synchrotron B-mode power at ℓ = 45 and ℓ = 74 to be 2.1%, expressed relative to the primordial spectrum with r = 0.2. This estimate is consistent with the values from Fuskeland et al. (2014). A more detailed and complete analysis of diffuse polarized synchrotron emission will appear in a forthcoming Planck paper.
Appendix E: Dust polarization decorrelation
![]() |
Fig. E.1 Ratio of the 217 ×
353 |
Extrapolation of the polarized dust signal observed at 353 GHz to lower frequencies (relevant to the CMB) hinges on the assumption that the pattern of dust polarization on the sky is independent of frequency, with only the amplitude changing, scaling as the SED. However, because the dust emission is the sum of many contributions from grains of different composition, size, shape, and alignment, which will change along the line of sight, there can in principle be some decorrelation of the signal from frequency to frequency (e.g. Tassis & Pavlidou 2014). We have seen in Sect. 4.5 that the dust polarization power spectra amplitudes are compatible with the typical dust SED, indicating that if there is such a decorrelation it has to be relatively weak. Nevertheless, in this Appendix, we will assess this possible effect quantitatively.
If there is a decorrelation of the dust polarization with frequency, one expects the
cross-spectrum between two frequencies to be smaller than the geometric mean of the two
auto-spectra. Hence, we compute the ratio between the cross- and auto-spectra for 217 and
353 GHz, where we have the sensitivity to quantify the effect. We compute
at these two frequencies and for the
cross-spectrum and fit a power-law amplitude ABB in each case, with
a slope fixed to αBB = −2.42, as in
Sect. 4.2.
The ratios are shown in Fig. E.1 as a function of ⟨I353⟩ for the
352 patches used in Sect. 5. Although there is a
large dispersion at low ⟨I353⟩ one can see that the ratios are
centred on unity, indicating that on average there is no decorrelation; among the 353
patches, we find a mean decorrelation ratio of
.
In Fig. E.1 we also show this ratio for the six LR
regions of Sect. 4 and the MB2 region of
Sect. 6, finding no evidence for a significant
decorrelation; for the LR regions, we specifically find
.
All Tables
All Figures
![]() |
Fig. 1 Masks and complementary selected large regions that retain fractional coverage of the sky fsky from 0.8 to 0.3 (see details in Sect. 3.3.1). The grey is the CO mask, whose complement is a selected region with fsky = 0.8. In increments of fsky = 0.1, the retained regions can be identified by the colours yellow (0.3) to black (0.8), inclusively. Also shown is the (unapodized) point source mask used. |
In the text |
![]() |
Fig. 2 Planck HFI 353 GHz |
In the text |
![]() |
Fig. 3 Best-fit power-law exponents αEE (red squares)
and αBB (blue
circles) fitted to the 353 GHz dust |
In the text |
![]() |
Fig. 4 Amplitude of the dust AEE (red squares) and ABB (blue circles) power spectra, normalized with respect to the largest amplitude for each mode. These are plotted versus the mean dust intensity ⟨I353⟩ for the six LR regions (top panel). A power-law fit of the form AXX(⟨I353⟩)= KXX⟨I353⟩1.9, X ∈ { E,B }, is overplotted as a dashed line of the corresponding colour (these almost overlap). The bottom panel presents the ratio of the data and the fitted ⟨I353⟩1.9 power law; the range associated with the ± 1σ uncertainty in the power-law exponent of 1.9 is shown in grey. For details see Sect. 4.3. |
In the text |
![]() |
Fig. 5 Ratio of the amplitudes of the |
In the text |
![]() |
Fig. 6 Frequency dependence of the amplitudes AEE,BB of the
angular power spectra, relative to 353 GHz (see details in Sect. 4.5). Results for
|
In the text |
![]() |
Fig. 7 Fitted dust |
In the text |
![]() |
Fig. 8 Top: map in orthographic projection of the 150 GHz
|
In the text |
![]() |
Fig. 9 Planck 353 GHz |
In the text |
![]() |
Fig. 10 Frequency dependence of the amplitude ABB of the
angular power spectrum |
In the text |
![]() |
Fig. A.1 Upper parts of both panels: performance of the
Xpol (blue points, see Sect. 3.1.1) and Xpure
(red squares, see Sect. 3.1.2) algorithms on
the |
In the text |
![]() |
Fig. B.1 Planck HFI 353 GHz |
In the text |
![]() |
Fig. B.2 Planck HFI 353 GHz |
In the text |
![]() |
Fig. B.3 Planck HFI 353 GHz |
In the text |
![]() |
Fig. C.1 Top: same as Fig. 7, but
instead of computing the |
In the text |
![]() |
Fig. C.2 Ratio of the dust ABB and AEE angular power spectra amplitudes computed on the 352 patches at high Galactic latitudes (described in Sect. 5, circles), as a function of the mean dust intensity ⟨I353⟩ of the patch. The mean value, ⟨ABB/AEE⟩ = 0.51 ± 0.18 here, is represented as a dashed-blue line and the associated ± 1σ dispersion as a light blue shaded area. The ABB/AEE ratio for the MB2 spectra presented in Sect. 6.3 is shown as an orange diamond. |
In the text |
![]() |
Fig. D.1 Xpol minus XpurePlanck 353 GHz
|
In the text |
![]() |
Fig. D.2 Planck 353 GHz |
In the text |
![]() |
Fig. D.3 Planck 353 GHz |
In the text |
![]() |
Fig. E.1 Ratio of the 217 ×
353 |
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.