arxiv: v1 [astro-ph.co] 1 Aug 2018

Transcription

1 Received XXXX; Revised XXXX; Accepted XXXX DOI: xxx/xxxx ARTICLE TYPE Confirmation of the detection of B-modes in the Panck poarization maps H. U. Nørgaard - Niesen arxiv: v1 [astro-ph.co] 1 Aug Nationa Space Institute (DTU Space), Technica University of Denmark Correspondence H. U. Nørgaard - Niesen, Emai: hunn@space.dtu.dk Present Address Eektrovej, DK-2800 Kgs. Lyngby, Denmark 1 INTRODUCTION Since the discovery of the Cosmic Microwave Background radiation (CMB) by Penzias and Wison(1965), investigations of this radiation has been a centra feature in order to get new information about the eariest evoution of the Universe. Due to technica imitations, unti recenty most of the observationa effort has been concentrated to obtain very accurate temperature measurements. Rees (1968) showed that poarization observations of CMB on arger scaes coud be a very important way to study the very eary history of the Universe. Symmetries in the production and growth of the poarization signa are constraining the configurations of the CMB poarization. Density (scaar) perturbations produce temperature (T) One of the main probems for extracting the Cosmic Microwave Background (CMB) from submm/mm observations is to correct for the Gaactic components, mainy synchrotron, free - free and therma dust emission with the required accuracy. Through a series of papers, it has been demonstrated that this task can be fufied by means of simpe neura networks with high confidence. The main purpose of this paper is to demonstrate that the CMB BB power spectrum detected in the Panck 2015 poarization maps is present in the improved Panck 2017 maps with higher signato-noise ratio. Two features have been detected in the EB power spectrum in the new data set, both with S/N 4. The origin of these features is most ikey eakage from E to B with a eve of about 1 per cent. This eakage gives no significant contribution to the detected BB power spectrum. The TB power spectrum is consistent with a zero signa. Atogether, the BB power spectrum is not consistent with the canonica tensor-to-scaar modes combined with gravitationa ensing spectra. These resuts wi give additiona strong arguments for support to the proposed poarization sateite projects to foow up on the Panck mission. KEYWORDS: Cosmic Microwave Background, Component Separation, method:neura network fuctuations and E (cur - free) poarization modes, whie tensor (gravitationa waves) produce both T, E and B (divergence free) modes. It is generay assumed that the initia inhomogeneities in the Universe were Gaussian distributed. Linear theory predicts that the CMB fuctuations are aso Gaussian, and that the CMB spectrum can be fuy described by 4 power spectra TT, EE, BB and TE whie the TB and EB power spectra shoud zero due to parity constraints (e.g. Kamionkowski et. (1997). CMB poarization measurements were provided by the WMAP sateite, the fina resuts given in Bennett et a. (2013). For r, the tensor-to-scaar ratio, Bennett et a. found, for WMAP-ony data, an upper imit of Ade et a. (2014) report a detection of B-modes in a 400 sq. degree area on the sky with the BICEP2 instruments, observing from the South Poe. These observations were done at ony one frequency, 150 GHz, impying that their estimate of

2 2 AUTHOR ONE poarized Gaactic emission in this sma area on the sky is uncertain. Ade et a. (2015) attacked this probem by combining BICEP2 + Keck 150 GHz data and Panck 30 GHz-353 GHz observations in this area. This investigation determined the ampitude of the ensing spectrum 1.12 ± 0.18, reative to the standard λcdm mode, and an upper imit on the tensor-toscaar ratio r < 0.13, with 95 percent confidence. The Panck Coaboration XV (2016) determined the CMB ensing potentia at a eve of 40σ, whie the Panck Coaboration XIII (2016) found an upper imit of 0.11 on the tensor-to-scaar ratio. For the Panck Coaboration, a key issue has aways been to carefuy investigate the data fows from the detectors and find ways to correct for systematic errors. Since the reease of the Panck data in 2015, this effort has been continued and significant improvement in the fina Panck frequency maps has been obtained. Since no new observations have been obtained between 2015 and 2017, the improvements are mainy due to remova of remaining systematic errors. In a series of papers, the capabiities of neura networks for deaing with mm/submm observations have been investigated by Nørgaard - Niesen & Jørgensen (2008), Nørgaard - Niesen & Hebert (2009), Nørgaard - Niesen (2010), and Nørgaard - Niesen (2012). The ast two papers showed, by using data from the WMAP sateite, that simpe neura networks can extract the CMB temperature and poarization signas with exceent accuracy. In Nørgaard - Niesen (2016, hereafter NN1) the Panck poarization frequency maps, reeased in 2015, the Stoke Q and U parameters were extracted by means of simpe neura networks. The BB power spectrum was detected in with S/N = 4.5. It was demonstrated that the contribution from Gaactic emission and remaining systematic errors in the maps were very sma, but they coud not be competey rued out. Due to the improvement in the Panck 2017 poarizations maps, it is feasibe to repeat the 2016 anaysis. As in NN1, the this paper is concentrated on the reiabiity of the detected poarization power spectra. The structure of the paper is the foowing: Sect. 2 and 3 give a short overview of the different anaysis toos appied in NN1, the power spectra from the extracted Q and U maps are presented in Sect. 4, whie the fuy caibrated power spectra are shown in Sect. 5, the contamination of the power spectra is discussed in Sect. 6. Concusions are given in Sect.7. 2 THE FINAL PLANCK POLARIZATION MAPS In the present anaysis, the fina Panck frequency maps (named DX12) have been expoited: ony the fina Q and U maps in the frequency range 30 GHz GHz has been expoited no attempt has been performed to correct the Panck frequency maps to the same point spread function the anaysis has been performed excusivey on the Panck frequency maps, impying that no auxiiary data or information has been expoited From the fina set of Panck frequency maps, 4 sets have been expoited in this paper: the 2 haf mission maps (HM1 and HM2) and the maps combined from the ODD years and the EVEN years of the mission. The precise definition of this data sets can be found in Panck Coaboration I (2017) 3 THE ROAD FROM FREQUENCY MAPS TO CMB POWER SPECTRA In this paper, the procedure for extracting the Stoke parameters Q and U from the fina Panck frequency maps foows cosey the procedure used in NN1. In this section, the key features of this procedure are outined. 3.1 The appied foreground mode As for the previous investigations, the main purpose is ony to extract the CMB signa. Therefore, meaningfu physica modes of the Gaactic foreground components are not required, ony a coherent mathematica mode of the non - CMB signa (a signas that are not reated to the CMB) is needed. The Independent Component Anaysis (Hyvarinen et a. 2001, ICA) has turned out to be we suited to provide a simpe mode of the non - CMB signa. Briefy, ICA is a mathematica method to extract a set of independent components, in this case from the 7 Panck frequency maps. Like in NN1, each of the HM1, HM2, ODD and EVEN data sets are spit into independent components by ICA. By anaysing these component maps, it is evident that 2 components contain neary a the extended Gaactic emission. For each of these components, ICA provides a normaized spectrum and an ampitude per sky pixe, which are used in setting up the neura networks, see Sect The neura network method In the series of papers referred above, it has been demonstrated that neura networks are a very effective too for extracted

3 AUTHOR ONE 3 CMB from mm/submm observations. The setup of the neura networks is the same as previousy: simpe muti ayer perceptron neura networks with 2 hidden ayers have been expoited the network is set up to extract signas with the same spectra behaviour as CMB the avaiabe frequency range is expoited the input for the neura networks are seven inputs channes (30 GHz GHz) and one output channe (Q or U), a in units of μk the standard non - inear activation function tansig has been expoited The scheme for setting up the data for training the neura network is: 1. defining the required intervas of the ampitudes of the 2 ICA components and for Q or U 2. draw 3 random numbers within these intervas 3. cacuate the combined spectrum 4. for each input channe, add random Gaussian noise within the range corresponding to the Panck hit maps in the different data sets 5. run these steps N(train) times, to get reasonabe number of spectra to train the network 6. train the network to estabish the transformation between the input channes and the output channe 7. test the accuracy of the network by running an independent TEST data set, with N(test) spectra obtained in the same way as the TRAIN data set, through the network. 8. if the accuracy is satisfactory (meaning the systematic errors as function of the input parameters are very sma and the residuas are Gaussian distributed), the Panck data set is run through the network The main advantages of using neura networks are: the number of weights needed to set up the fu transformation between the 7 input channes and 1 output channe (Q or U) is in our case 80, sma compared to the number of pixes in each of the Panck frequency maps ( ) the foreground mode is ony using the input frequency maps themseves (neither auxiiary data nor assumptions about e.g. the spectra behaviour of the different foreground components are appied) a sma number of parameters are needed in the foreground mode, since a physica meaningfu mode is not required FIGURE 1 The CC HM and OE EE power spectra together with the mean. It is seen that these spectra fit nicey together up to = Y - axis: ( + 1) (2π)C EE 3.3 Appied sky mask In order to avoid residua contamination by emission from the MIky Way in the extracted power spectra, the apodized CAM- SPEC mask used in NN1, covering about 63 per cent of the sky, has aso been appied here. Since the contribution from point sources is expected to be insignificant, it is aso negected in this work. 3.4 The Panck simuations In the process of estabishing the Panck 2015 data, the whoe Panck data chain has been extensivey simuated (FFP8). In NN1, this data set was mainy used for caibrating the derived power spectra to an absoute scae. The basic simuations behind FFP8 has been improved, but from a caibration point of view, the improvements are not significant. Therefore, the NN1 caibration scheme has aso been foowed here. 3.5 The high pass fiter It turned out that specia treatments of the Panck 2015 maps at mutipoes ess than 20 were needed, mainy due to the uncertainty in the contribution of arge scae systematic errors sti in the frequency maps. Athough there have been significant improvements in the frequency maps since 2015, the discussion of this range of mutipoe moments is outside the goa of

4 4 AUTHOR ONE TABLE 1 Definition of the standard Panck high pass fiter, defined with 1 = 20 and 2 = 40 (Panck 2015 IX) range vaue < < < 2 0.5[1 cos(π 1 )] 2 1 > FIGURE 3 The CC EB power spectrum. Y - axis: ( + 1) (2π)C EB 4.1 The EE power spectra The CC HM and EO EE power spectra together with the mean are shown in Fig. 1. It is seen that the HM and OE power spectra are very simiar up to FIGURE 2 The CC HM and OE BB power spectra, compared with the 2015 BB power spectrum (NN1). Taking the uncertainty into account, a reasonabe agreement is seen between the <2015> and <2017> spectra. Y - axis: (+1) (2π)C BB this investigation, and the standard Panck high pass fiter has been introduced (Tabe 1) as in NN1. 4 THE POWER SPECTRA OF THE OBSERVED PLANCK POLARIZATION MAPS The accuracy of cross correation (CC) power spectra is superior compared to auto correation (AC) power spectra. Therefore, the discussions in this paper are concentrated on the HM1 <-> HM2 (HM) and the ODD <-> EVEN (OE) CC power spectra. A power spectra anaysed in this paper have been extracted with the HEALPix anafast IDL routine (Gorski et a. (2005)). The - interva of the power spectra are the same as for the Panck 2015 EE power spectrum. 4.2 The BB power spectra The derived CC HM and OE BB power spectra and the mean are shown in Fig. 2, together with the CC BB power spectra from NN1. Taking the accuracy of these power spectra into account a reasonabe agreement is seen. It is seen that the main features in the 2015 spectrum (NN1) is sti present. 4.3 The EB power spectrum The CC EB power spectrum shows 2 features around = 675 and = 375, each detected with a S/N 4 (Fig.3 ). The positions of these feaures correpond to the postions of the first 2 peaks in the EE power spectrum (Fig.1 ). 4.4 The TB power spectrum The CC HM TB spectrum is presented in Fig. 4. It is cear that no part of this spectrum has been significanty detected. For the average power is ± 0.13, 5 THE CALIBRATION OF THE POWER SPECTRA The basic caibration scheme used in NN1 (mainy corrections for the sky mask, the sensitivity of the Panck detectors and for the effective point spread function of the NN method) has been appied. Briefy, these corrections were estimated by combining a set of simuated FFP8 frequency maps with reaistic noise

5 AUTHOR ONE 5 FIGURE 4 The CC TB power spectrum, it is evident that this spectrum is consistent with a zero signa. Y - axis: ( + T B 1) (2π)C FIGURE 6 The fuy caibrated CC <HM, OE> BB power spectrum, compared with PLANCK 2015 estimates of the BB power spectrum. CAMB BB tensor spectra with tensor to scaar ratios of 0.5 and 0.3 are aso shown. Y - axis: ( + 1) (2π)C BB FIGURE 5 The fuy caibrated CC <HM, OE> EE power spectrum, compared with the Panck 2015 EE power spectrum (Panck 2015 XI) incuding error bars. An exceent agreement is seen for < Y - axis: ( + 1) (2π)C EE maps and theoretica tensor-to-scaar maps, cacuated by the CAMB software package (see 6 THE FULLY CALIBRATED EE, BB AND EB POWER SPECTRA The mean, fuy caibrated, CC <HM, OE> EE spectrum are given in Fig. 5 together with the Panck 2015 EE power spectrum, a good agreement is seen. As in NN1, the error bars given in the figures are extracted from the scatter around the spectra in each - interva. In Fig. 6, the mean, fuy caibrated CC <HM, OE> BB power spectrum is shown. The excess from = 100 to = 275, discussed in NN1, are confirmed with improved S/N = 6.6 FIGURE 7 The Signa to Noise of the 2017 BB power spectrum, compared with previous obtain spectrum. A considerabe improvement is evident (compared to 4.5 in NN1 ), the feature has a S/N of 4.1. The theoretica CAMB spectra with tensor to scaar ratios of 1.0, 0.5, 0.3 are aso shown together with the CAMB ensing spectrum. It is cear that a combination of these CAMB spectra wi have difficuties in fitting the observed spectrum. In Fig. 7, the difference in S/N between the BB power spectrum found in NN1 and the spectrum given here is highighted. The improvement is evident. The fuy caibrated CC EB power spectrum is shown in Fig. 8. The 2 features around = 675 and = 375 have S N = 4.4, and 4.2, respectivey. These features coud indicate a vioation of the parity assumption, but the coincidence with the

6 6 AUTHOR ONE positions with the first 2 peaks of the EE power spectrum coud aso indicate a eakage from E to B. To expain the strength of the EB features, a eakage of 1-2 per cent is sufficient. At this stage, it is difficut to excude a eakage of this eve. Such a eakage woud give a insignificant contribution to the BB power spectrum compared to the noise eve at the positions of the features. In Fig. 1, it is seen that the EE power is much weaker than the first peaks in the - range reevant for the BB spectrum discussed in this paper, impying that the E to B eakage gives no significant contribution to the detected BB power spectrum. 7 REMAINING SYSTEMATICS EFFECTS FROM THE REDUCTION OF THE PLANCK FREQUENCY MAPS As emphasized above, important improvements in removing residua systematics have been obtained by the LFI and HFI instrument teams since As emphasized in the two Panck 2015 detector papers (Panck Coaboration III and VIII (2016)) the systematic errors on arger scaes are the most difficut to remove, but the high-pass-fiter (Tabe 1 ) is designed to remove this kind of systematics. Panck Coaboration III (2017) gives a detaied summary of a the systematic errors found during the reduction phase of the HFI data and imits on their ampitudes. From their Tabe 7 it is seen that main contributor of remaining systematic errors in the frequency maps is the uncertainties in the ADC inear corrections. For the 100 GHz, 143 GHz and 217 GHz channes (the most sensitive channes to the CMB signa on Panck) the residua eve of systematic errors is μk 2, peaking at sma s, impying that they give ony minima contributions to the cross power spectra presented here. For the other important contribution to systematics e.g. caibration and bandpass mismatch are neary an order of magnitude ower. Further, the neura networks are design to extract components with a fat frequency spectrum. Atogether, it is reasonabe to assume that the observed BB and EB cross spectra are not pouted significanty by residua systematic errors. The E/B eakage due the uncertainties in the beams of the Panck optica system is estimated to be ess than 3x10 3 μk 2, which is sma compared to the detected BB power spectrum (see Fig. 6). 8 THE E/B LEAKAGE Due to the fact that most, if not a, CMB experiments extract the power spectra from ony a fraction of the sky (due to FIGURE 8 The fuy caibrated CC <HM, OE> EB power spectrum. Y - axis: ( + 1) (2π)C EB the strong signa from the MIky Way in the reevant frequency range) the E modes wi to some extend eak into the B modes. This probem has been discussed in a number of papers e.g. Kendrick (2006) and Grain et a. (2009). They show the probem can be significanty reduced by a proper apodization of the mask. Most of these investigations have anaysed the probem if sma areas on the sky, say 1 per cent, are observed, reevant for e.g. baoon fights. Since the Panck Mission is covering the whoe sky, it is possibe to anayse a arge part of the sky, it is ony necessary to excude the brightest area of the Miky Way. As emphasized in Section 3.3, data from about 60 percent of the sky has been expoited in this paper. To investigate the eve of eakage in the detected BB power spectrum, the BB power spectra has been cacuated from the CAMB T/S = 1.0 modes with different masks appied. The eakage (defined as the difference of the power spectrum incuded a mask minus the fu sky power spectrum). The masks appied are simiar to the mask described in Section 3.3, except for the size of the sky coverage. In Fig. 9, it is seen that the eve of the E/B eakage for the mask used in cacuated the BB power spectrum (the purpe ine) is insignificant compared to the detected spectrum(fig. 6). 9 CONTAMINATION BY GALACTIC EMISSION As emphasized above, the neura networks are setup to extract a signa corresponding to a fat spectrum, but at ow fux eves the noise of the observations presents a chaenge for the neura networks. A simpe test to estimate the contribution of the Gaactic components is to investigate how the cross power spectra are depending of the Gaactic atitude. In Fig. 10 the BB power spectra obtained within different CAMSPEC masks with different sky coverage ( ) are shown. It is seen that

7 AUTHOR ONE 7 FIGURE 9 The BB(mask) corrected - BB(fu sky), for the CAMB T/S = 1 mode. It is seen that the E/B eakage due to the appied mask is much smaer than the detected BB power spectrum. Y - axis: ( + 1) (2π)C EB FIGURE 11 The CC synchrotron HM 1 and NN HM 2 EE power spectrum. The synchrotron HM 1 map is provided by the Panck Commander team. Due to the observationa errors in the LFI channes, they have appied a ow pass fiter removing mutipoemoments above 250. It is evident that there is no evidence for a contamination of the NN maps by synchrotron emission from the Miky Way. X - axis: δ = 1, Y - axis: ( + 1) (2π)C EE FIGURE 10 The <HM, OE> CC BB power spectra obtained with a set of CAMSPEC masks, with sky covering factors of 0.26, 0.40, 0.52, 0.61, The spectra have been normaized to the covering factor of the mask used in obtained the HM and ODd/EVEN spectra. If the accuracy of the spectra is taken into account thee are no trend as function of the sky coverage. Y - axis: ( + 1) (2π)C BB within the accuracy there is no variation of the spectra as function of the sky coverage, indicating a non - significant Gaactic contribution. 9.1 The synchrotron emission It is generay accepted that ony the Gaactic synchrotron emission and the therma dust emission are poarized. The best avaiabe maps of these components have been derived by the FIGURE 12 The CC synchrotron HM 1 and NN HM 2 BB power spectrum. The synchrotron HM 1 map is provided by the Panck Commander team. Due to the observationa errors in the LFI channes, they have appied a ow pass fiter removing mutipoemoments above 250. It is evident that there is no evidence for a contamination of the NN maps by synchrotron emission from the Miky Way. X - axis: δ = 1,Y - axis: ( + 1) (2π)C BB Commander team from the Panck 2017 data set. Due to the uncertainties in the LFI channes they have a appied a ow pass fiter. In Figs. 11 and 12, the CC EE and BB power spectra between the Commander synchrotron maps and the corresponding NN maps are shown. It is evident that there is

8 8 AUTHOR ONE negigibe synchrotron contribution to the derived NN power spectra. is aso negigibe contribution from dust to the derived NN power spectra. 10 CONCLUSIONS FIGURE 13 The CC dust HM 1 and NN HM 2 EE power spectrum. The dust HM 1 map is provided by the Panck Commander team. It is evident that there is no sign of a contamination of the NN maps by therma dust emission from the Miky Way. X - axis: δ = 1,Y - axis: ( + 1) (2π)C EE It has been demonstrated that with the improved accuracy of the fina Panck poarization maps compared the 2015 reease, the detection of the BB power spectrum in NN1 is confirmed, with a higher confidence. Possibe contamination from Gaactic emission and remaining systematic errors in the Panck frequency maps has been rued out. Two features in the EB power spectrum have been found each detected with a S/N 4. At this stage, the most ikey origin of these features is eakage from E modes to B modes of the order of 1 per cent. This eve of eakage gives no significant contribution to the detected BB power spectrum. The TB power spectrum is found to be consistent with a zero spectrum. The confirmation of the BB power spectrum wi, no doubt, give new strong arguments for the proposed poarization missions to foow up on Panck. 11 ACKNOWLEDGEMENTS The author acknowedges that this work woud not have been possibe without the massive efforts by a ot of strongy committed scientists and engineers within the Panck Coaboration. A anonymous referee is acknowedged for vauabe comments. This work has taken advantage of the HEAPix software package. REFERENCES FIGURE 14 The CC dust HM 1 and NN HM 2 BB power spectrum. The dust HM 1 map is provided by the Panck Commander team. It is evident that there is no evidence for a contamination of the NN maps by therma dust emission from the Miky Way. X - axis: δ = 1,Y - axis: ( + 1) (2π)C BB 9.2 The dust therma emission Simiary, in Figs. 13 and 14, the CC EE and BB power spectra between the Commander therma dust emission maps and the corresponding NN maps are shown. It is evident that there Ade et a.: 2014, Phys. Rev. Lett.112, Ade et a.: 2015, arxiv Bennett et a., the WMAP Team 2013, ApJS, 208, 20 Gorski et a., 2005, ApJ, 622, p.759 Grain et a., 2006, arxiv Hyvarinen et., 2001, New York Wiey Kamionkowski et a.: 1997, Phys. Rev. Lett.78, 2058 Nørgaard-Niesen, H. U., Jørgensen, H. E., 2008, Ap&SS, 318, 195 Nørgaard-Niesen, H. U., Hebert, K.: 2009, AN 330, 863 Nørgaard-Niesen, H. U.: 2010, A&A, 520, A2 Nørgaard-Niesen, H. U.: 2012, Ap&SS, 340, 161 Nørgaard-Niesen, H. U.: 2016, AN, 337, 662 Penzias, A. A. & Wison, R. W., 1965, ApJ, 142, 419 Panck 2015 Resuts III: 2016, A&A, 594, A3 Panck 2015 Resuts VIII: 2016, A&A, 594, A8 Panck 2015 Resuts XI: 2016, A&A, 594, A11 Panck 2015 Resuts XIII: 2016, A&A, 594, A13 Panck 2015 Resuts XV: 2016, A&A, 594, A15 Panck 2017 Resuts I: 2017 Panck 2017 Resuts III: 2017 Rees, M: 1968, ApJ153, 11

9 AUTHOR ONE 9 Smith, KK. M.: arxiv How cite this artice: Nørgaard - Niesen, H. U. 2017, Confirmation of the detection of B-modes in the Panck poarization maps, AN, XXX. How cite this artice: Nørgaard - Niesen, H. U. 2017, Confirmation of the detection of B-modes in the Panck poarization maps, AN, XXX.