arxiv: v1 [astro-ph.co] 21 Mar 2013

Transcription

1 Nonparametric Forecasts of the CMB Anguar Power Spectra for the Panck Mission Amir Aghamousa Centre for Modeing and Simuation, University of Pune, Pune India arxiv: v [astro-ph.co] Mar 03 Department of Physics, University of Pune, Pune India amir@cms.unipune.ac.in Mihir Arjunwadkar Nationa Centre for Radio Astrophysics, University of Pune Campus, Pune India Centre for Modeing and Simuation, University of Pune, Pune India, mihir@ncra.tifr.res.in and Tarun Souradeep Inter-University Centre for Astronomy and Astrophysics, University of Pune Campus, Pune India tarun@iucaa.ernet.in ABSTRACT The Panck mission, designed for making measurements of the cosmic microwave background (CMB) radiation with unprecedented accuracy and anguar resoution, is expected to reease its first data set in the near future. For the first time in the CMB history, extensive measurements of the CMB poarization wi be made avaiabe for the entire sky. Such precise and rich data are expected to contain a great weath of information about the Universe. The information in the CMB data is convenienty represented in terms of anguar power spectra for temperature and poarization. A proper estimation of these CMB power spectra from data is the first step in making inferences about the Universe and, in particuar, cosmoogica parameters that govern the Universe. In this paper, we provide forecasts for the T T, EE, and T E anguar power spectra for the Panck mission. Our forecasts are made using synthetic data based on the best-fit ΛCDM mode whie conforming to the characteristics and parameters of the Panck mission. We have anaysed such synthetic data sets using a nonparametric function estimation formaism that is otherwise asymptoticay modeindependent. Indeed, the power spectra estimated from such synthetic data turn out to

2 be sufficienty cose to the corresponding ΛCDM spectra. Our anaysis further indicates that the three anguar power spectra are determined sufficienty we for max, where max = 500 (T T ), 377 (EE), and 77 (T E) respectivey. A signature of reionization is evident in the EE and T E power spectra in the form of a bump at ow s. Nonparametric confidence bands in the phase shift (φ m ) vs. acoustic scae ( A ) pane, corresponding to the first eight peaks in the T T power spectrum, show a confuence region for 300 A 305, which is in good agreement with the estimate A = 300 based on the best-fit ΛCDM mode. From our resuts, we expect that the rea Panck data, when reeased, shoud ead to accurate mode-independent estimates of CMB anguar power spectra using our nonparametric regression formaism. Subject headings: cosmic background radiation cosmoogica parameters Methods: data anaysis Methods: statistica. Introduction Observation of the microwave sky reveas that the temperature of the cosmic microwave background (CMB) radiation is not exacty the same in a directions. These sma fuctuations in the temperature are imprinted on the entire sky, impying that the CMB is anisotropic. These primordia anisotropies were first discovered in 99 by the COsmic Background Exporer (COBE). This was foowed up by a remarkabe series of ground-based and baoon-borne experiments, and most recenty by the Wikinson Microwave Anisotropy Probe (WMAP). These fuctuations are beieved to have been generated within 0 35 seconds of the Big Bang. CMB anisotropies are therefore a rich source of information about the eary Universe, and have revoutionized the way we understand our Universe. A study of CMB anisotropies aso heps in probing the fundamenta physics at energy scaes much higher in magnitude compared to those accessibe to partice acceerators. CMB anisotropies are sensitive to the cassica cosmoogy parameters such as expansion rate, curvature, cosmoogica constant, matter content, radiation content, and baryon fraction, and provide insights for modeing structure formation in the Universe (Hu 996). For exampe, measurements of the CMB anisotropies with ever-increasing precision have made it possibe to estabish a standard cosmoogica mode that asserts that the Universe is neary spatiay fat (Smoot 007). The CMB contains a far greater weath of information about the Universe through its (inear) poarization. The CMB has acquired inear poarization through Thomson scattering during either decouping or reionization, sourced by the quadrupoe anisotropy in the radiation distribution at that time (Rees 968; Hu & White 997). This quadrupoe anisotropy coud originate in different sources (Hu & White 997): the scaar or density fuctuations, vector or vortica fuctuations, and tensor or gravitationa wave perturbations. Scaar mode represents the fuctuations in energy density of the cosmoogica pasma at the ast scattering that causes a veocity distribution which eads to bue-shifted photons. Vector mode represents the vorticity on the pasma which causes

3 3 Dopper shifts that resut into quadrupoar obes. However, such vorticity wi be damped by expansion (as are a motions that are not enhanced by gravity), and its effect are expected to be negigibe. Tensor perturbation are the effect of gravity waves which stretch and squeeze space in orthogona directions. This aso stretches the photon waveength, and hence produces a quadrupoar variation in temperature. The dependence of CMB poarization on cosmoogica parameters differs from that of temperature anisotropies. As such, it provides additiona constraints on cosmoogica parameters that hep break degeneracies. Accurate measurements of CMB poarization wi therefore enabe us to affirm the vaidity/consistency of different cosmoogica modes (Zadarriaga 004; Zadarriaga et a. 997; Eisenstein et a. 999). CMB poarization information is further expected to hep determine initia conditions and evoution of structure in the Universe, origin of primordia fuctuations, existence of any topoogica defects, composition of the Universe, etc. The Panck mission (The Panck Coaboration 006) is a space-based fu-sky probe for thirdgeneration CMB experiments designed for this purpose. Indeed, one of the main objectives of this mission is to measure the primordia fuctuations of the CMB with an accuracy prescribed by the fundamenta astrophysica imits, through improvements in sensitivity and anguar resoution, and through better contro over noise and confounding foregrounds. The higher anguar resoution of Panck impies that higher-order peaks in the CMB anguar spectra can be determined with better precision, which in turn transates to determination of cosmoogica parameters (such as baryon and dark matter densities) with improvement in statistica precision by an order of magnitude. Furthermore, the Panck mission is designed to make extensive measurements of the E-mode poarization spectrum over mutipoes up to 500 with unprecedented precision, together with good contro over poarized foreground noise. These measurements are expected to provide insights into the physics of eary Universe, epoch of recombination, structure formation, aowabe modes of primordia fuctuations (adiabatic vs. isocurvature modes), reionization history of the Universe, and hep in estabishing constraints on the primordia power spectrum. Panck wi aso hep constrain the fundamenta physics at high energies that are not possibe to probe through terrestria experiments (The Panck Coaboration 006). In our previous work (Aghamousa et a. 0), we estimated the CMB T T power spectrum from four WMAP data reaizations using a nonparametric function estimation methodoogy (Beran 000; Genovese et a. 004). This methodoogy does not impose any specific form or mode for the power spectrum, and determines the fit by optimizing a measure of smoothness that depends ony on characteristics of the data. This ensures that the fit and the subsequent anaysis is approximatey mode-independent for sufficienty arge data sizes. Further, this methodoogy quantifies the uncertainties in the fit in the form of a high-dimensiona eipsoida confidence set that is centered at this fit and captures the true but unknown power spectrum with a pre-specified probabiity. This confidence set is the prime inferentia object of this methodoogy which aows addressing compex inferentia questions about the data meaningfuy in a unified framework.

4 4 The Panck mission is expected to reease high-accuracy CMB data sets in the near future. In this paper, we therefore attempt to forecast the three CMB power spectra for the CMB to see what can be expected (and inferred) from rea Panck data when it gets reeased, when anaysed using this nonparametric regression methodoogy. For this purpose, we use synthetic Panck-ike data conforming to the specifications and parameters of the Panck mission. This synthetic data is based on the assumption that the best-fit ΛCDM mode is the true mode of the Universe. In what foows, Sec. briefy describes the three CMB anguar power spectra, foowed by a description (Sec. 3) of the synthetic data used in this work. This is foowed by the resuts (Sec. 4) and a concusion (Sec. 5).. Temperature and poarization power spectra Poarized radiation is typicay characterized in terms of the Stokes parameters I, Q, U and V (Jackson 998). The parameter V describing circuar poarization is not generated by Thomson scattering for the CMB. The two Stokes parameters Q and U describing inear poarization form the components of a rank- symmetric trace-free tensor P ab. Since any two-dimensiona symmetric tensor can be represented using two scaar fieds, P ab can be represented using to the eectric (i.e., gradient) P E and the magnetic (i.e., cur) P B modes of poarization. This decomposition is unique, and is simiar to the decomposition of a vector fied into a gradient and a divergence-free vector fied (Chainor 004). The P E and P B scaars are defined on a sphere, and can therefore be expanded in a spherica harmonics basis as ( )! ( + )! ae m Y m(θ, φ) () P E (θ, φ) = P B (θ, φ) = m m ( )! ( + )! ab m Y m(θ, φ), () which define the E- and B-mode mutipoes a E m and ab m respectivey (Chainor 004; The Panck Coaboration 006; Hu & White 997). Any CMB measurement can be decomposed into three maps (T, E and B respectivey). A tota of six anguar power spectra (T T, EE, BB, T E, T B, and EB) can be obtained from these three components (de Oiveira-Costa 005). These six power spectra are defined by expanding the T, E and B maps in terms of spherica harmonics, resuting into the foowing correation structure: < a Y m ay Y m >= CY δ δ mm, (3) where Y, Y are E, B or T. In the absence of parity vioation, and under the assumption of Gaussian fuctuations, the temperature and poarization anisotropies of the CMB are described by the C T T, C EE, C BB, C T E power spectra competey (Chainor 004; Hu & White 997; Zadarriaga 004). Since the B-mode poarization is not expected to be detected by the Panck mission

5 5 accuratey (The Panck Coaboration 006), we focus ony on the C T T, C EE spectra. and C T E power The standard deviation of the C T T, C EE, C T E Panck Coaboration 006; Knox 995) ( C T T ) ( C EE ) ( C T E ) power spectra is approximatey given by (The (C T ( + )f T + ω T W ) sky (4) (C E E + ω P ( + )f ) sky (5) ( (C E ( + )f E + ω P W )(C T T + ω T W ) + (C T E)), sky (6) where ω T = (σ p,t θ F W HM ), and ω P = (σ p,p θ F W HM ) are the weights per soid ange for temperature and poarization, and f sky is the fraction of observed sky in the experiment. The σ p,t and σ p,p are noise standard deviations per resoution eement (θ F W HM θ F W HM ). The window function for a Gaussian beam is where beam = 8 n (θ F W HM ). W = exp ( ( + )/( beam )), (7) The importance of decomposing the P ab tensor in terms of the E and B modes comes from the fact that inear scaar perturbations do not generate any B-mode poarization (Hu & White 997). The tensor mode contributes to both E as we as B modes, whereas the vector mode contributes ony to the B-mode poarization. Therefore the E part of the decomposition stems from the scaar and tensor modes, and the B part originates in the vector and tensor modes (The Panck Coaboration 006; Chainor 004; Hu & White 997). Cosmoogica reevance of these poarization modes is as foows. The E-mode mainy foows the veocity of the cosmoogica pasma at decouping. Compared to the temperature anisotropies, which originate in photon density fuctuations at the ast scattering, the E mode therefore contains more information about some of the cosmoogica parameters (Zadarriaga et a. 997; Chainor 004), and can ead to better estimates of parameters such as the baryon and cod dark matter densities. Poarization power spectra have an osciatory structure that is anaogous to that of the T T power spectrum. For exampe, peaks in EE power spectrum are out of phase with those in the T T power spectrum due to anisotropy generated at the ast scattering. The T E power spectrum, which has a higher ampitude compared to the EE power spectrum, is a measure of the correations (positive or negative) between density and veocity fuctuations (Scott & Smoot 0; Chainor 004). The phase difference between acoustic peaks in the T T, EE and T E power spectra can be used as a mode-independent check for the physics of acoustic osciations (Chainor 004). Adiabatic and isocurvature perturbations aso have different effects on the phase of the poarization spectra: Predicted poarization power spectra for isocurvature perturbations show out-ofphase peaks and dips compared to those for adiabatic perturbations such that the power spectra

6 6 from isocurvature perturbations appear to be -shifted versions of those for adiabatic perturbations (Sievers et a. 007). A meaningfu estimation of the poarization power spectra can therefore be used to determine which of the two scenarios is coser to truth. Another cosmoogica phenomenon that affects the poarization spectra is reionization: Reionization of Universe started when the first generation of stars started producing a fux photons. The resuting free eectrons started re-scattering the CMB radiation. Athough ony a sma fraction of CMB photons got scattered this way during the reionization era, the imprints of reionization are expected to be seen as distortions in the poarization power spectra at arge anguar scaes of the order of 0 degrees. The height and ocation of the reionization bump (Zadarriaga 997; Kapinghat et a. 003; Hoder et a. 003) expected at ow mutipoes ( 0) has information reated to tota optica depth and the reionization epoch redshift (Zadarriaga et a. 997; Kapinghat et a. 003). Athough the precision of reionization bump detection is imited by cosmic variance at ow (Hu & Hoder 003), constraining it wi hep understand the reionization history better, and break degeneracies between severa cosmoogica parameters by constraining the optica depth better (Zadarriaga et a. 997; Eisenstein et a. 999). 3. Synthetic data for the Panck mission For forecasting the CMB anguar power spectra for the Panck mission, we generate the synthetic Panck-ike data using the FuturCMB code (Perotto et a. 006). FuturCMB generates a simuated anguar power spectrum using a user-provided theoretica power spectrum C true (which is assumed to be the true spectrum) for frequency channes representing Panck measurements, and generates the corresponding noise power spectrum N conforming to the Panck characteristics. This is done by generating a random reaization of the spherica harmonics coefficients a m, assumed to be Gaussian random variabes with mean zero and variance where C true Var(a m ) = C true + N, (8) is a proxy for the true but otherwise unknown anguar power spectrum. For C true, we use spectra generated using CAMB (Lewis et a. 000) for the best-fit ΛCDM cosmoogica parameters obtained from the WMAP 7-year data. We aso imit FuturCMB to 500, a range that corresponds to the three Panck frequency channes (00, 4 and 7 GHz). N is the noise power spectrum given by N = ω W, (9) where ω is ω T and ω P for temperature and poarization respectivey, and W is the window function for a Gaussian beam (Eq. 7). The noise in the T E power spectrum is taken to be zero because noise contributions from different maps are uncorreated (Perotto et a. 006). FuturCMB then

7 7 ( + )C π Fig.. A reaization of simuated T T power spectrum data for the Panck mission, generated using FuturCMB (Perotto et a. 006). Back points: data incuding noise; red points: simuated data after subtracting noise. ( + )C π Fig.. A reaization of simuated EE power spectrum data for the Panck mission, generated using FuturCMB (Perotto et a. 006). Back points: data incuding noise; red points: simuated data after subtracting noise.

8 8 ( + )C π Fig. 3. A reaization of simuated T E power spectrum data for the Panck mission, generated using FuturCMB (Perotto et a. 006). FuturCMB assumes the noise to be zero for the T E data. cacuates the power spectra data C map as C map = ( + ) + m= a m. (0) Eq. 0 is an unbiased estimator of Eq. 8; its expected vaue is therefore equa to C true + N. The noise spectra are expected to dominate over the true spectrum for sufficienty high vaues of. This is seen in figures 3, where the back points represent the FuturCMB output C map of the T T, EE and T E spectra: the upward trends in the tai of T T and EE spectra are the resut of noise dominating the data at high s. The T E power spectrum, on the other hand, does not show any such upward trend because the noise in the T E spectrum is assumed zero. To obtain synthetic T T and EE data, we therefore subtract the corresponding noise spectra from the C map output of FutureCMB (Figures, red points). The covariance matrix of the simuated anguar power spectra is taken to be a diagona matrix with diagona eements defined in Eq. 4, 5 and 6 for T T, EE, and T E respectivey. 4. Resuts and discussion Nonparametric fits to synthetic Panck data. For estimating the power spectra from synthetic Panck data, we use the nonparametric regression method described in (Aghamousa et a.

9 9 ( + )C π Fig. 4. T T nonparametric fits. Bue, fu-freedom fit (EDoF 7); red, restricted-freedom fit (EDoF=7); back: best-fit ΛCDM spectrum; grey: simuated data reaization. ( + )C π Fig. 5. EE nonparametric fits. Bue, fu-freedom fit (EDoF 90); red, restricted-freedom fit (EDoF=4); back: best-fit ΛCDM spectrum; grey: simuated data reaization.

10 0 ( + )C π Fig. 6. T E nonparametric fits. Bue, fu-freedom fit (EDoF 95); red, restricted-freedom fit (EDoF=40); back: best-fit ΛCDM spectrum; grey: simuated data reaization. 0). Whie our synthetic data is generated under the assumption that the ΛCDM mode as estimated from the WMAP 7-year data is the true mode of the Universe, this formaism for nonparametric regression and inference itsef does not make any assumptions about the shape of the true regression function underying the data; it is asymptoticay mode-independent. A nonparametric fit, under this methodoogy, can be characterized by its effective degrees of freedom (EDoF), which can be thought of as the equivaent of the number of parameters in a parametric regression probem. Using this methodoogy, we obtain nonparametric fits to the synthetic T T, EE and T E data by appropriatey constraining the EDoF of the fits. (Additiona detais about our nonparametric fits can be found in Sec. A and B.) Figures 4, 5 and 6 show nonparametric fits (red curves) to the T T, EE, and T E data respectivey, which are in good agreement with the underying ΛCDM spectra (C true, back curves) used to generate the synthetic data. This shows that this nonparametric regression methodoogy, which does not assume any specific form of the true (but generay unknown) regression function, can recover the underying true spectrum with high accuracy especiay where noise eves are not too high. How we wi the Panck fits be determined by data aone? To see how noise in the data affects oca uncertainties in a fitted spectrum, we compute approximate 95% confidence intervas for each fitted C using 5000 randomy samped spectra from the corresponding confidence set. The ratio of this confidence interva to the absoute vaue of the fitted C (assumed to be nonzero) is a reative measure of how we each fitted C is determined (Genovese et a. 004; Aghamousa et a.

11 EE =377 TE =77 TT =46 Tota Vertica Variation in the Fit/Fit Fig. 7. The resuts of a probe of the confidence sets for the T T (red), EE (bue), and T E (green) nonparametric restricted-freedom fits to the synthetic Panck data, to determine how we the fits are expected to be determined by the data aone. The quantity potted for each data reaization is the tota vertica variation at each within the respective 95% (σ) confidence set, divided by the absoute vaue of the fit (assumed nonzero): Vaues indicate that the fit is tighty determined by the data, whereas vaues indicate that the data contain very itte information about the height of the anguar power spectrum at that. Disregarding the ow- region for the EE fit, and spikes for the T E fit (which arise from neary zero fitted C vaues), the marked vertica ines indicate the approximate -vaue at which each curve rises above. 0): A vaue impies that the fit is we determined by the data, and a vaue impies that the data contain very itte information about height of the power spectrum at that. In Figure 7, we pot this reative error (95%) for a three spectra as a function of the mutipoe index. We see that, by this criterion, the Panck power spectra are expected to be we-determined up to 46(T T ), 377(EE) and 77(T E). Since the T E fit osciates around zero (Figure 6), this quantity takes very high vaues at s where the T E spectrum has a neary zero vaue. This resuts in mutipe spikes in Figure 7 (green dash curve), but this does not impy that the fit is i-determined at these s. Ignoring these spikes, we see that the reative error in the T E fit is beow unity up to 77, which indicates the range over which this fit is expected to be we-determined by data aone. Uncertainties on the ocations and heights of peaks and dips. Locations and heights of peaks and dips in the CMB power spectra contain information about cosmoogica modes and

12 parameters. Uncertainties in the ocation and height of a peak or a dip in a fitted spectrum can thus hep assess uncertainties in the vaues of reated parameter. Foowing the procedure outined in (Aghamousa et a. 0), we samped the 95% confidence set of each fitted spectrum uniformy to generate spectra variations whie ensuring that at east 5000 of these are acceptabe (see Sec. B for detais). Figures 8, 9, and 0 show the resuts of this exercise, together with tabuated vaues in Tabes,, and 3 respectivey. The box around a peak or a dip represents the argest horizonta and vertica variations in the scatter; these represent the 95% confidence intervas on the ocation and height of a peak or a dip. Fig % confidence boxes the ocations and heights of peaks and dips in the T T fit. Back curve is the restricted-freedom monotone fit to the synthetic Panck T T data (grey points). The number of acceptabe spectra variations samped from the 95% confidence set is These uncertainties are tabuated in Tabe. For the T T fit, these boxes around peaks and dips are tiny (Figure 8), which is a refection of the accuracy of the Panck T T data. Such precise determination of peaks and dips shoud ead to more robust estimates of reated cosmoogica parameters than what is currenty avaiabe. The theoretica T T power spectrum (Figure 4, back curve) shows a sma upturn at ow. This upturn is primariy the resut of the integrated Sachs-Wof (ISW) effect. In Figure 8, this upturn corresponds to the first dip in the spectra variations samped from the 95% confidence set. Peaks and dips in the EE power spectrum show reasonaby ow uncertainties up to 00 (Figure 9). Beyond this, the data contain high eves of noise, and therefore a uncertainties become much arger. This is in agreement with the behavior of the EE curve in Figure 7 (bue curve). We aso expect a sma bump in the EE power spectrum which is reated to the epoch of reionization.

13 3 Fig % confidence boxes the ocations and heights of peaks and dips in the EE fit. Back curve is the restricted-freedom monotone fit to the synthetic Panck T T data (grey points). The number of acceptabe spectra variations samped from the 95% confidence set is These uncertainties are tabuated in Tabe. This bump is indeed seen in both nonparametric EE fits in Figure 5. Figure 9 shows a (tiny) 95% uncertainty box for this bump. A precise determination of the height and ocation of this peak (Tabe ) shoud revea usefu information about the epoch of reionization. The uncertainty boxes on peaks and dips in the T E fit (Figure 0) are reasonaby sma ti 800, again in agreement with the resut depicted in Figure 7 (green dashed curve). The somewhat pecuiar uncertainty boxes on the ast three peaks in the T E fit are due to the fact that there are two spurious peaks at = 90 and = 070 in the restricted-freedom fit (see Sec. A) which are a resut of the high noise eves at high s. At ow mutipoes, the T E fit aso shows a bump (Figure 6) which is reated to the epoch of reionization. Athough a simiar bump in the EE fit is known to be more informative (The Panck Coaboration 006), a determination of this peak in the T E spectrum shoud aso ead to usefu information about reionization. Are the acoustic peaks in the T T and EE spectra out of phase with respected to each other? From the fundamenta physics of the CMB anisotropies, we expect the acoustic peaks in the T T and EE power spectra to be out phase with respect to each other. One way of estabishing this is by considering the ratio of peak ocations in the EE fit to the corresponding ones in the T T fit, which shoud be (m + 0.5)/m for the mth peak (Sievers et a. 007). We depict the 95%

14 4 Fig % confidence boxes the ocations and heights of peaks and dips in the T E fit. Back curve is the restricted-freedom monotone fit to the synthetic Panck T T data (grey points). The number of acceptabe spectra variations samped from the 95% confidence set is These uncertainties are tabuated in Tabe 3. peaks ocations of TT power spectra versus EE one in Figure shows the 95% confidence intervas on peak ocations in the EE fit (red) potted against the corresponding uncertainties in the T T fit (bue). Aso potted are the peak ocation pairs corresponding to the best-fit ΛCDM mode (back dots), and points (green) based on the theoretica expectation (m + 0.5)/m. A the potted quantities are by and arge consistent with each other, indicating that the expected behavior of out-of-phase peak ocations is indeed vindicated by the data. An estimate of the acoustic scae parameter A. Tabe ists the 95% confidence intervas on peak and dip ocations and heights for the T T power spectrum fit. As a way of iustrating the roe of these uncertainties in the estimation of cosmoogica parameters, we consider the foowing reationship (Hu et a. 00; Doran & Liey 00) between the ocation m of the mth peak, the acoustic scae A, and the shift parameter φ m for TT power spectrum: m = A (m φ m ). If we substitute the end-points of the 95% confidence interva for the mth peak ocation, then this reationship resuts into hyperboic confidence bands in the A φ m pane (Figure ). The intersection of these bands (for the first 8 peaks in the T T fit) determine an estimated confidence interva for the acoustic scae 300 A 305 which is in agreement with the reported vaue A = 300 by (Page et a. 003). In comparison with our previous estimate based on the WMAP 7-year data (Aghamousa et a. 0), the current estimate has improved remarkaby. Furthermore,

15 5 Fig.. The 95% confidence intervas on peak ocations in the EE fit (red) potted against the corresponding confidence intervas in the T T fit (bue). Aso potted are the peak ocation pairs corresponding to the best-fit ΛCDM mode (back dots), and points (green) based on the theoretica expectation (m + 0.5)/m. A the potted quantities are by and arge consistent with each other, indicating that the expected behavior of out-of-phase peak ocations is indeed vindicated by the data. any additiona information about the phase shifts φ m can ead to an even more refined estimate for acoustic scae. 5. Concusion In this paper, we have addressed the question of what coud be the expected outcome of a nonparametric anaysis of the T T, EE and T E power spectrum data sets from the Panck space mission when they get reeased. For this purpose, we have used synthetic/simuated Panck-ike data sets based on the best-fit ΛCDM mode, and have anaysed them using a nonparametric regression and inference methodoogy. Our resuts show that the T T power spectrum can be estimated with such a high accuracy that a peaks are resoved up to 500. We expect that the EE and T E

16 6 m = 8 φ m m = 7 m = 6 m = 5 m = 4 m = 3 m = m = A Fig.. Confidence bands for the acoustic scae A and the shift φ m for the mth peak, as derived from the 95% confidence intervas on the first eight peak ocations of estimated TT power spectrum. power spectra wi be reasonaby we-determined and can hep, e.g., better understand reionization history, address the issue of adiabatic versus isocurvature perturbations, and estimate of some of the cosmoogica parameters precisey. As a resut, we expect to have a better understanding of the Universe via the Panck data even from an agnostic, nonparametric approach. A. Fu-freedom and restricted-freedom nonparametric fits The fu-freedom fits are obtained by minimizing the risk function subject to monotonicity constraints on the shrinkage parameters (Aghamousa et a. 0). Such fu-freedom fits can be quite osciatory especiay where the noise eves in the data are high. Despite this fit being a reasonabe fit (in the sense that it captures the essentia trend in the data we), a cosmoogica modes expect far smoother anguar power spectra. To account for this, we minimize the risk again,

17 7 Peak Location Peak Height Dip Location Dip Height : (04, 34) h : ( , 595.6) + : (406, 4) h + : (638.5, ) : ( , ) : (55, 549) h : (57.08, 69.86) + 3 : (805, 87) h 3 : ( , ) 3+ 4 : (0, 4) h 4 : (5.64, 67.8) 4+ 5 : (48, 439) h 5 : ( , ) 5+ 6 : (74, 76) h 6 : ( , ) 6+ 7 : (989, 085) h 7 : (9.873, 5.404) 7+ 8 : (8, 498) h 8 : (7.93, ) 8+ : (666, 688) h + : (009, 08) h + : (30, 35) h + : (630, 666) h + : (9, 007) h + : (0, 44) h + : (394, 500) h + : (969.94, 09.34) : ( , ) : ( , ) : (99.838, ) : ( ,.5749) : (5.7434,.0498) Tabe : 95% Confidence Interva on Severa Features of TT Anguar Power Spectrum Peak Location Peak Height Dip Location Dip Height : (7, 6) h : ( 0.005, 0.035) + : (7, 7) h + : ( 0.009, 0.08) : (0.460, ) : (33, 44) h : (0.934,.54) + 3 : (390, 40) h 3 : (0.4654, 3.58) 3+ 4 : (68, 697) h 4 : (35.0, ) 4+ 5 : (978, 009) h 5 : ( , ) 5+ 6 : (58, 360) h 6 : (4.445, ) 6+ 7 : (44, 80) h 7 : (7.339, ) 7+ 8 : (66, 999) h 8 : (.46, ) 8+ : (93, 05) h + : (5, 53) h 3+ : (85, 844) h 4+ : (35, 84) h 5+ : (40, 668) h 6+ : (46, 969) h 7+ : (755, 000) h 8+ : (5.444, ) : (8.669, 4.80) : (4.9459, ) : ( , 0.600) : ( , 8.670) : ( , ) Tabe : 95% Confidence Interva on Severa Features of EE Anguar Power Spectrum restricting the EDoF of the fit to a vaue ess than that for the fu-freedom fit. We continue reducing the EDoF in this way unti we obtain an acceptaby smooth fit. We ca this the restricted-freedom fit. The fu-freedom fit for the T T power spectrum is amost smooth, but with tiny wigges (Figure 4, bue curve), and corresponds to EDoF 7. We achieve a numericay smooth fit for the T T power spectrum at EDoF 7 (Figure 4, red curve). We see that this restricted-freedom fit foows the fu-freedom fit amost exacty, but without the additiona wigges which resuted from noise in the data. For the EE power spectrum, the fu-freedom fit has EDoF 89 (Figure 5, bue curve). This fit is quite wiggy at high mutipoes, as can be expected from high noise eves in data there. The Panck proposa (The Panck Coaboration 006) aso expects the EE power spectrum to have high noise at high mutipoes, and hence unabe to resove peaks beyond 000. Therefore, we have excuded the EE fit beyond = 000 (Figure 5) from further anaysis. The restricted-freedom smoother version of the EE fit corresponds to EDoF= 3.79 (Figure 5, red curve). We see again that the smooth fit foows the fu-freedom fit, but averages over the wigges. For the simuated T E data, we obtain the fu-freedom fit at EDoF 95 (Figure 6, bue curve). The peaks are resoved we up to 000. At higher mutipoes, we see a wiggy fit with fase peaks due to high eves of

18 8 Peak Location Peak Height Dip Location Dip Height : (3, 38) h : (0.609,.58) + : (4, 58) h + : ( , ) : ( , ) : (30, 36) h : (3.89, 36.57) + 3 : (588, 605) h 3 : (6.370, ) 3+ 4 : (903, 93) h 4 : (56.478, ) 4+ 5 : (99, 43) h 5 : (4.335, 3.50) 5+ 6 : (5, 579) h 6 : (.957, 5.787) 6+ 7 : (77, 954) h 7 : ( 4.494, ) 7+ 8 : (884, 08) h 8 : ( 0.883, ) 8+ 9 : (03, 354) h 9 : (.805, ) 9+ 0 : (55, 499) h 0 : ( , 75.59) 0+ : (46, 477) h + : (74, 758) h 3+ : (06, 086) h 4+ : (353, 400) h 5+ : (650, 750) h 6+ : (8, 08) h 7+ : (959, 8) h 8+ : (073, 49) h 9+ : (57, 500) h 0+ : ( , 4.938) : ( , ) : ( , ) : ( , ) : ( 4.00,.888) : ( ,.4636) : ( , 0.379) : ( 6.37, ) Tabe 3: 95% Confidence Interva on Severa Features of TE Anguar Power Spectrum noise. A smoother restricted-freedom fit is obtained at EDoF= 40 (Figure 6, red curve). Despite this additiona smoothing, because of the high noise at high mutipoes, we see two fase bumps at = 90, 070, and a fase peak at = 304. B. Probing the confidence sets A high-dimensiona confidence set (for a prespecified eve of confidence; typicay, 95%) around a fit is the prime inferentia object for this nonparametric methodoogy (Aghamousa et a. 0). By probing the confidence set, we can determine uncertainties on specific feature of a fit, vaidate different cosmoogica modes against the data, etc. For finding uncertainties on the ocation and height of a peak or a dips (Figures 8, 9, and 0) we use the same methodoogy as was used in Aghamousa et a. (0): We uniformy sampe the confidence set of each smooth fit, and record the peaks and dips of acceptabe spectra thus samped. The most extreme variations in the height or ocation define a confidence interva on the respective estimates (represented in the figures as boxes). Compared to Aghamousa et a. (0), we change a few criteria for accepting such samped spectra to account for the greater anguar resoution of the Panck data: For the T T and T E fits, we seect spectra with exacty 8 and 0 peaks respectivey. In the case of the T E fit, this criterion incudes a peak at ow mutipoes ( 50) and a fase peak at = 304. For the EE fit, we sampe spectra with 8 peaks for 000, together with the additiona condition that there must be a peak at ow mutipoes ( 50). The ast peak in the EE restricted-freedom fit consists of two tiny but cose bumps due to noise in the data. Such fine structure is not expected here on theoretica grounds. Therefore, in case we find a samped spectrum with two tiny peaks around the ocation of the eighth peak in the fit, but with separation ess than 0 mutipoes, we consider them as a singe peak and record their average ocation and height.

19 9 REFERENCES Aghamousa, A., Arjunwadkar, M., & Souradeep, T. 0, Astrophys. J., 745, 4 Beran, R. 000, J. Amer. Statist. Assoc., 95, 55 Chainor, A. 004, in Lecture Notes in Physics, Berin Springer Verag, Vo. 653, Lecture Notes in Physics, Berin Springer Verag, ed. Papantonopouos, E., 7 de Oiveira-Costa, A. 005, in Astronomica Society of the Pacific Conference Series, Vo. 343, Astronomica Poarimetry: Current Status and Future Directions, ed. Adamson, A., Aspin, C., Davis, C., & Fujiyoshi, T., 485 Doran, M., & Liey, M. 00, Mon. Not. Roy. Astron. Soc., 330, 965 Eisenstein, D. J., Hu, W., & Tegmark, M. 999, Astrophys. J., 58, Genovese, C. R., Mier, C. J., Nicho, R. C., Arjunwadkar, M., & Wasserman, L. 004, Statist. Sci., 9, 308 Hoder, G. P., Haiman, Z., Kapinghat, M., & Knox, L. 003, The Astrophysica Journa, 595, 3 Hu, W. 996, in Lecture Notes in Physics, Berin Springer Verag, Vo. 470, The Universe at High-z, Large-Scae Structure and the Cosmic Microwave Background, ed. Martinez-Gonzaez, E., & Sanz, J. L., 07 Hu, W., Fukugita, M., Zadarriaga, M., & Tegmark, M. 00, Astrophys. J., 549, 669 Hu, W., & Hoder, G. P. 003, Phys. Rev. D, 68, 0300 Hu, W., & White, M. 997, New Astronomy,, 33 Jackson, J. D. 998, Cassica Eectrodynamics Third Edition, 3rd edn. (Wiey) Kapinghat, M., Chu, M., Haiman, Z., Hoder, G. P., Knox, L., & Skordis, C. 003, Astrophys. J., 583, 4 Knox, L. 995, Phys. Rev. D, 5, 4307 Lewis, A., Chainor, A., & Lasenby, A. 000, Astrophys. J., 538, 473 Page, L., et a. 003, Astrophys. J. Supp. Ser., 48, 33 Perotto, L., Lesgourgues, J., Hannestad, S., Tu, H., & Wong, Y. Y. Y. 006, Journa of Cosmoogy and Astropartice Physics, 006, 03 Rees, M. J. 968, Astrophys. J., Lett., 53, L Scott, D., & Smoot, G. 0, Phys. Rev. D, 86, 0000

20 0 Sievers, J. L., et a. 007, Astrophys. J., 660, 976 Smoot, G. F. 007, Rev. Mod. Phys., 79, 349 The Panck Coaboration. 006, ArXiv Astrophysics e-prints, astro-ph/ Zadarriaga, M. 997, Phys. Rev. D, 55, 8 Zadarriaga, M. 004, Measuring and Modeing the Universe, 309 Zadarriaga, M., Sperge, D. N., & Sejak, U. 997, Astrophys. J., 488, This preprint was prepared with the AAS L A TEX macros v5..