arxiv: v1 [astro-ph.co] 28 Jun 2017

Transcription

1 YITP-SB Infationary Features and Shifts in Cosmoogica Parameters from Panck 015 Data arxiv: v1 [astro-ph.co] 8 Jun 017 Georges Obied, 1 Cora Dvorkin, 1 Chen Heinrich, Wayne Hu, and Vinicius Miranda 3 1 Harvard University, Department of Physics, Cambridge, MA 0138, USA Kavi Institute for Cosmoogica Physics, Department of Astronomy & Astrophysics, Enrico Fermi Institute, University of Chicago, Chicago, IL Center for Partice Cosmoogy, Department of Physics and Astronomy, University of Pennsyvania, Phiadephia, Pennsyvania 19104, USA We expore the reationship between features in the Panck 015 temperature and poarization data, shifts in the cosmoogica parameters, and features from infation. Residuas in the temperature data at ow mutipoe, which are responsibe for the high H 0 70 km s 1 Mpc 1 and ow σ 8Ω 1/ m vaues from < 1000 in power-aw ΛCDM modes, are better fit to infationary features with a 1.9σ preference for running of the running of the tit or a stronger 99% CL oca significance preference for a sharp drop in power around k = 04 Mpc 1 in generaized sow ro and a ower H 0 67 km s 1 Mpc 1. The same in-phase acoustic residuas at > 1000 that drive the goba H 0 constraints and appear as a ensing anomay aso favor running parameters which aow even ower H 0, but not once ensing reconstruction is considered. Poarization spectra are intrinsicay highy sensitive to these parameter shifts, and even more so in the Panck 015 TE data due to an outier at 165, which disfavors the best fit H 0 ΛCDM soution by more than σ, and high H 0 vaue at amost 3σ. Current poarization data aso sighty enhance the significance of a sharp suppression of arge-scae power but eave room for arge improvements in the future with cosmic variance imited E-mode measurements. I. INTRODUCTION The anisotropies of the cosmic microwave background (CMB) continue to be one of the most powerfu probes we have to study physica conditions in the eary universe. With the reease of Panck 015 data, we now have access to precise a-sky measurements of the poarization as we as the temperature fuctuations of the CMB. Observations from supernovae, baryon acoustic osciations, Big Bang Nuceosynthesis, among other datasets support, in addition to the CMB, a broady consistent cosmoogica mode termed ΛCDM. Nonetheess, evidence for features beyond ΛCDM in the CMB temperature data have been caimed (e.g. [1 9]), as we as disagreements between CMB predictions for oca observabes under ΛCDM and their measured vaues (e.g. the oca expansion rate H 0 [7, 10]). The two are reated because these oca cosmoogica inferences depend on the assumptions in the ΛCDM mode, in particuar the form of the infationary power spectrum. It is, thus, essentia to consider the impact of one type of deviation from ΛCDM on the other. Moreover poarization data pay an essentia roe in breaking degeneracies between the two. It is the aim of this paper to examine these reationships between features in the temperature and poarization data, shifts in cosmoogica parameters and infationary features. Previous studies have tested the consistency of Panck data at different stages in the chain of inference from the raw data to cosmoogica parameters (e.g. [11 13]). In particuar, Ref. [13] spit the temperature power spectrum into two disjoint mutipoe ranges, the ower being simiar to the range of WMAP, and anayzed the two ranges separatey. They then find a 3σ discrepancy between parameters derived from the two parts of the same dataset. The Panck coaboration [14] carried out a meticuous investigation of the cause of this tension and discovered that it was mainy due to a deficit in power at ow- and osciatory residuas in the mutipoe range More recenty, Ref. [15] examined the consistency of the poarization and temperature datasets, finding no strong evidence to disfavor the ΛCDM mode. In addition, severa works [5, 16 6] have expored the effects of infationary features on CMB observabes, incuding the centra roe poarization data pay in confirming or rejecting such features. For instance, stepike features in the infationary potentia coud cause a dip and bump in the temperature data at ow. Such modes coud arise from symmetry-breaking phase transitions in the eary universe, among other reasons [7 30]. These features invariaby vioate the sow-ro approximation and, thus, require more sophisticated modeing of the resutant infationary power spectrum. In this paper, we adopt both the traditiona running of the tit type parameters and the generaized sow ro (GSR) formaism [3, 31 33], which aows the incusion of extra spine basis parameters that accommodate order unity power spectrum variations from sow ro predictions. We study the shifts in cosmoogica parameters, such as H 0 and the ampitude of oca structure, when these extra parameters are incuded and expore the aspects of the temperature and poarization data that drive them. This paper is organized as foows. Sec. II describes the datasets and parameter inference techniques. The resuts of our anayses are presented in Sec. III, which focuses on the shifts in cosmoogica parameters, and in Sec. IV, which focuses on the impications for infationary

2 Dataset Likeihood EE φφ binned + ow TEB binned TEEE + ow TEB ens reconstruction TABLE I. Combinations of Panck 015 data sets and ikeihoods. a In addition, we spit the binned and TEEE ikeihoods mutipoe ranges, the atter with a modified ikeihood code. a Specificay: ow TEB = ow SMW 70 dx11d v5c Ap binned = pik dx11dr HM v18 binned TEEE = pik dx11dr HM v18 TEEE ens reconstruction = smica.g30 ft fu pp, Mode Parameters {τ, θ MC, Ω b h, Ω ch } ΛCDM +{A s, n s} rλcdm +{A s, n s, n run, n run,run} SB +{A s, n s, p i} rsb +{A s, n s, n run, n run,run, p i} TABLE II. Modes and parameterizations of the infationary power spectrum. The parameters p i are spine basis coefficients that generaize the tit as defined in Appendix A. In each case, the 4 ate time cosmoogica parameters are prepended to the parameter ist. features. We concude in Sec. V. II. DATA AND MODELS For the anayses described in this paper, we use the pubicy avaiabe Panck 015 data, which incude the power spectra of cosmic microwave background (CMB) temperature and poarization fuctuations. For the ow mutipoe ( 9) range, we use the standard pik owteb ikeihood code suppied by the Panck coaboration. For the high mutipoe ( 30) range, we simiary use the Panck pik binned *1 ikeihood for the baseine and the TEEE ikeihood to assess the impact of the poarization data. Note that the ow- poarization data are incuded in a anayses, even though we refer to these cases as and EE respectivey. We aso ater the range of the different ikeihoods to probe the impact of specific regions of the data, in the case of the TEEE ikeihood using a custom modification of the code. In some cases, we use information from the Panck ensing power spectrum C φφ in the mutipoe range We aso incude the standard foreground parameters in each anaysis. * These datasets and ikeihoods are summarized in Tab. I. *1 We have separatey tested that the unbinned ikeihoods give statisticay indistinguishabe resuts. * Premarginaizing foregrounds assuming ΛCDM with pik-ite as in Ref. [14] can change inferences on the infationary running parameters. Conversey, not using the whoe range in to con- We empoy a Bayesian approach to infer parameter posterior distributions and to derive confidence imits, namey the Markov Chain Monte Caro (MCMC) technique impemented via the pubicy avaiabe CosmoMC code [34] inked to a modified version of the Botzmann code CAMB [35]. For each combination of datasets and modes, we run 4 MCMC chains unti convergence, determined by the Geman-Rubin criterion R 1 < 1. For the parameters, we take in a cases the standard 4 ate-time ΛCDM cosmoogica parameters: the optica depth to reionization τ, the (approximate) anguar scae of the sound horizon θ MC, the baryon density Ω b h and the cod dark matter density Ω c h. Tensions in the ΛCDM cosmoogy are mainy associated with Ω c h, which determines the caibration of the physica scae of the sound horizon and, hence, the CMB inference of the Hubbe constant H 0 as we as the growth of structure and ensing observabes represented by σ 8 Ω 1/ m. These two parameters expose tensions in the ΛCDM mode when compared to oca, non-cmb measurements. Here and throughout the paper, σ 8 is the rms of the inear density fied smoothed with a tophat at 8h 1 Mpc and Ω m = Ω c + Ω b + Ω ν. We assume one massive neutrino with m ν = 6 ev and the usua N eff = We summarize our mode and parameter choices in Tab. II. In Fig. 1 we show the parameter derivatives or responses of the, TE and EE power spectra around the best fit ΛCDM mode to the Panck dataset given in Tab. III. Here we normaize the derivatives to σ, the cosmic variance per mutipoe, where, T T ; σ = +1 CT T 1 +1 C T T +1 CEE C EE + (C T E ), T E;, EE. Here and beow σ is fixed by the fiducia mode when varying parameters. For and EE this rescaing produces fractiona derivatives scaed by 1/ at 1, which has the benefit of removing features from parameters ike normaization that simpy echo the acoustic peaks in the spectrum. We have aso rescaed the parameters to either the error in the parameter from the ΛCDM mode, e.g. for Ω c h we show the parameter p = Ω c h /σ(ω c h ), or by a fixed number, e.g. for n run, p = 50n run for carity. For the cases where we scae to σ, this has the added benefit of showing where current imits can be improved with better measurements. Increasing Ω c h decreases the overa power in the first few acoustic peaks in due to the reduction of radiation driving and the eary ISW effects from decaying potentias inside the sound horizon (see [36] Fig. 4). Since the former effect changes the ampitude of acoustic osciations, it aso reduces the temperature peaks, which (1) strain foregrounds aows < 1000 modes that do not provide a good goba fit to foregrounds (see Tab. III and [14]).

3 3 ( + 1)C/π TE EE σ 1 dc/dp Ω b h /σ Ω c h /σ θ MC /σ 500τ n A s e τ /σ n s /σ 50n run 70n run,run p 1 p 5 p 10 p FIG. 1. Response of the, TE, EE power spectra to various parameters around the ΛCDM mode of Tab. III. Derivatives are scaed by σ, the cosmic variance error per of the mode, and parameters are scaed by either their Panck 015 ΛCDM error σ or by a fixed number for visibiity. Vertica dotted ines are drawn at the acoustic peak ocations and horizonta dashed ines at zero eves. are extrema of the osciations, reative to troughs, which are nodes. Before the damping tai, these changes are neary in phase with the acoustic peaks in contrast to θ MC, which makes out-of-phase changes that shift peak positions. Damping tends to shift the observed peak ocations to arger anges, the phasing of the two drifts at high, but we wi nonetheess refer to them as in phase and out of phase osciatory changes since our focus is on the peaks that are we measured by Panck. Note that with higher mutipoe information from ground based ex-

4 4 Parameter ΛCDM ΛCDM ( < 1000) rλcdm rλcdm + φφ SB rsb 100Ω b h Ω ch θ MC τ n(10 10 A s) n s H 0 [km s 1 Mpc 1 ] σ 8Ωm 1/ TABLE III. Cosmoogica parameters of best fit modes for various datasets. In the ΛCDM ( < 1000) mode the best fit assumes the foreground mode of the best fit ΛCDM mode, unike the MCMC resuts where foregrounds are jointy fit. This assures a good goba fit when foregrounds are i-constrained by subsets of the data. periments this effect can be distinguished from entirey in phase effects. Indeed increasing Ω c h aso increases the gravitationa ensing, which aso reduces the peaks reative to the troughs and remains entirey in phase through the damping tai. For the Panck data, the acoustic effect of Ω c h is arger than the ensing effect but they both contribute [14]. The impact of Ω c h on the EE spectrum is much sharper near the sound horizon ( 00, first poarization trough, see Fig. 3 in Ref. [37]) since potentia decay does not change it after recombination as it does temperature through the eary ISW effect. This enhanced sensitivity carries over to the TE spectrum, which dominates the Panck poarization constraints at high. To infer the 4 ate-time parameters from CMB data we must assume a parameterized form for the infationary curvature spectrum R, and vice versa. We choose to mode R in 4 ways to highight its impact. In a cases, A s is the ampitude R (k 0) = A s where k 0 = 8 Mpc 1. Fig. 1 shows that this choice paces the pivot point at 10 3, which is near the best constrained scae of the Panck data. It aso shows the effect of raising ensing by raising A s in the osciatory response beyond the third peak. For the sope d n R /d n k, we consider a constant tit n s 1 as the baseine case, which we ca ΛCDM. Unike the ate-time cosmoogica parameters, A s and n s coherenty change the spectrum across a mutipoes rather than introducing ocaized features in C /σ (see Fig. 1). On the other hand, the dominant source of tension in the ΛCDM mode comes from anomaies in the spectrum that are more ocaized. Fig. 1 shows that by aowing running and running of the running of the tit: d n ( ) R k d n k = n s 1 + n run n + n ( ) run,run k n, k 0 k 0 () one can change the ow-mutipoe part of the spectrum reative to the high-mutipoe part, but ony if these additiona parameters are of order the n s 1 itsef. In standard sow ro infation, each order of running is suppressed by an additiona factor of n s 1, which shoud make them safey unobservabe in the Panck data. We ca the case where these parameters are not constrained by the sow ro hierarchy rλcdm. In these modes, R wi have a prominent but broad feature on CMB scaes. Whie these parameters can mode arge-ampitude deviations from sow ro, they cannot mode rapid variations that occur on a time scae of an efod or shorter. Such rapid variations woud be required to produce sharp features in mutipoe space in the observed power spectra. To distinguish between rapid variations that are confined to ow mutipoes and a smooth change in the overa amount of arge to sma-scae power, we empoy the generaized sow-ro (GSR) approximation. Here, d n R d n k G, (3) where G is a function of efods of the sound horizon during infation, as detaied in Appendix A. We consider arbitrary variations δg around both the standard n s 1 ( SB ) and running forms ( rsb ) for scaes between the horizon and the first acoustic peak. Specificay, we take a spine basis with 0 knots p i across decades in scae. In Fig. 1 we show the response of the observabe power spectrum to these parameters. Note that individuay these parameters produce step-ike features in the power spectrum and their superposition can produce fine scae features to fit the ow- anomaies in the data. In Fig. 1 these steps appear tapered due to rescaing by cosmic variance errors σ 1 1/. They produce sharper features in EE at a sighty higher mutipoe due to projection effects. Since the GSR approximation is sti perturbative, we impose an integra constraint on the second order impact of the p i parameters, max(i 1 ) < 1/ (see Eq. A7). Finay, we assume no tensors, r = 0 throughout. Cosmoogica parameter constraints derived from the MCMC ikeihood anayses are summarized in Fig., and best fit modes are isted in Tab. III. We use these modes in the next two sections to highight the interpay of features in the data, shifts in cosmoogica parameters,

5 5 ΛCDM rλcdm SB rsb ( < 1000) EE ( < 1000) EE Ω b h Ω c h θ MC ( < 1000) EE ( < 1000) EE τ og (10 10 A s ) n s ( < 1000) EE ( < 1000) EE H 0 [km s 1 Mpc 1 ] A s e τ FIG.. Cosmoogica parameter constraints with the 4 combinations of datasets in Tab. I on the 4 modes for the infationary power spectrum in Tab. II σ 8 Ω 0.5 m and features in infation. III. H 0 AND COSMOLOGICAL PARAMETERS Of the 4 ate-time parameters, the one most dependent on the dataset and mode assumptions of the previous section is the cod dark matter density Ω c h (see Fig. ), which directy contros the inferences on H 0 and the amount of structure today σ 8 Ω 1/ m. Raising Ω c h increases the scae of the sound horizon reative to the current horizon because of the reduced effect of radiation on the former. Compensating this effect requires a ower Hubbe constant to increase the distance to recombination. In ΛCDM the net effect is to keep Ω m h 3 neary constant [38]. As discussed in Ref. [14], this sensitivity is argey driven by the Panck residuas from the ΛCDM best fit mode in two regions: power deficits and gitches at ow 40 and osciatory residuas at As shown in Fig. 1, Ω c h contros the overa ampitude of the acoustic peaks reative to ow mutipoes, and aso impacts the height of the peaks reative to the troughs. Given the abiity of the infationary power spectrum to change the

6 6 C T T /σ n L ΛCDM ( < 1000) SB FIG. 3. Top pane: residuas in the data (points) and modes reative to the best fit ΛCDM mode of Tab. III (dashed ine). Bottom pane: change in the ikeihood n L of modes reative to ΛCDM starting from the owest bin at < 30 (dotted ines) and adding each individua data point thereafter. The ΛCDM ( < 1000) mode (red ine) fits the residuas at ow 40 and better due to changes in H 0 and other parameters, but it is a poor fit to the in-phase residuas at > Conversey, the infationary feature SB mode (bue) fits the ow- residuas yieding a arge n L 1 improvement at ow, but ony sma changes at high. Note that the ΛCDM ( < 1000) mode has the same foreground parameters as the ΛCDM mode, and the data points assume foregrounds are fixed to the best fit ΛCDM vaues in Ref. [9]. ΛCDM ΛCDM ( < 1000) SB SB ( < 1000) rλcdm rλcdm +φφ rsb 0.5 Posterior C T T /σ H 0 [km s 1 Mpc 1 ] FIG. 4. H 0 from the fu vs ( < 1000) datasets. Low- residuas in Fig. 3 drive H 0 higher in the ΛCDM ( < 1000) anaysis (red vs back curves). In the SB mode, these residuas are fit by infationary rather than cosmoogica parameter changes, so ow vaues of H 0 that are favored in the fu anaysis (gray shaded) are aso compatibe in the ( < 1000, red shaded) anaysis FIG. 5. residuas at > 1000 reative to the best fit ΛCDM mode of Tab. III (dashed ine). With the running parameters, rλcdm and rsb seek to fit the osciatory residuas and make the peaks smoother. Adding the ens reconstruction data φφ argey removes this preference. reative power and amount of smoothing of the peaks due to ensing, the interpretation of these residuas depends on the infationary mode assumptions. In addition, as shown in Ref. [9], the TE data have a surprisingy arge impact on Ω c h and H 0, favoring a ow vaue for the atter. We further expore these issues and their reationship to infationary features in turn.

7 7 Posterior ΛCDM rλcdm rλcdm +φφ rsb H 0 [km s 1 Mpc 1 ] FIG. 6. H 0 and infationary parameters. In phase > 1000 residuas in Fig. 3 drive H 0 ower than the ΛCDM mode (back) but require running parameters (red curve, rλcdm; red shaded, rsb) in order to produce a goba fit. These fits aso predict more gravitationa ensing and are not favored by ensing reconstruction data (gray shaded, rλcdm +φφ.) A. Low 40 residuas We begin by repeating the test discussed in Ref. [13] of comparing the inferences from the and ( < 1000) datasets on H 0 in the ΛCDM mode context. In Fig. 4, we show that constraints shift from H 0 = [67.3 ± 0.98] km s 1 Mpc 1 to [69.7 ± 1.8] km s 1 Mpc 1, as expected. In Ref. [14], this shift is argey attributed to the ow- residuas between the data and the best fit ΛCDM mode to the fu- data, which we dispay in Fig. 3. *3 These residuas favor ower ow- power compared with the acoustic peaks. As seen in Fig. 3 (red curve), they can be fit better by the ΛCDM ( < 1000) mode by owering Ω c h, which increases the radiation driving to raise the peaks over ower s with adjustments of other parameters, such as n s (see Tab. III). In addition, the outof-phase residuas peaking between the third and fourth acoustic peaks drive a sma change in θ MC. With the fu range in, these same changes cannot fit the data, as we discuss further in the next section. To further test this interpretation, we show the impact of generaizing the mode cass to SB which aows the ow- residuas to be fit by infationary features instead. The maximum ikeihood SB mode of Tab. III is shown in Fig. 3 (bue curve). By design, *3 These data points have a fixed best fit foreground mode under ΛCDM subtracted. In order to compare the ΛCDM ( < 1000) mode to the residuas, we take the best fit with the foreground mode fixed by the fu range in Tab. III. C/C [( + 1)] C φφ /π [ 10 7 ] ΛCDM rλcdm rλcdm + φφ FIG. 7. Lensing power spectrum in the ΛCDM mode vs. running modes rλcdm and rλcdm +φφ. Lowering H 0 in the rλcdm mode simutaneousy raises ensing and is not favored once the ensing reconstruction data (points) are taken into account. the SB parameters do not significanty ater the acoustic peaks (see Fig. 1) eaving the high- residuas and cosmoogica parameters neary unchanged from ΛCDM. Correspondingy, in Fig. 4 we show that for the fu- data set, H 0 constraints aso remain argey unchanged, whereas for < 1000 they broaden, making these ower H 0 vaues compatibe. From Fig., there is a simiar broadening in the running casses of modes rλcdm and rsb. We return to the distinctions between these various means of fitting residuas at ow in Section IV (see aso Fig. 1). B. High 1000 residuas In the fu dataset, the ower Ω c h and higher H 0 that is preferred by the ow- residuas causes even more significant in-phase high- deviations from the data. Higher radiation driving beyond the 4th peak impies sharper peaks than those seen in the data (see Fig. 5). Indeed, there remains in-phase osciatory residuas of opposite sign at high > 1000 even with respect to the fu ΛCDM maximum ikeihood mode. The data favor smoother peaks than those that can be achieved in the ΛCDM context, producing the so-caed ensing anomay, since additiona ensing acts to smooth the peaks (see Fig. 1). Specificay, when aowing the ampitude of ensing C φφ A L C φφ to foat independenty, A L = 1. ± 0.10

8 8 Posterior ΛCDM ΛCDM ( < 1000) ΛCDM EE ΛCDM EE ( < 1000) H 0 [km s 1 Mpc 1 ] FIG. 8. H 0 and poarization data. TE residuas shown in Fig. 9, especiay near 00, drive H 0 ower, even in the < 1000 datasets (red curves to red shaded), whereas the fu data range provide vaues consistent with ΛCDM, with smaer errors. [9]. Aternatey, if ony > 1000 data are considered in the ΛCDM context, H 0 = [64.1 ± 1.7] km s 1 Mpc 1 as pointed out by Ref. [13]. However, given an infationary power spectrum described by A s and n s as in ΛCDM, fitting the osciatory residua in this way does not ead to a consistent goba soution [14]. Reaxing the assumptions on the infationary power spectrum aows for greater fexibiity. In Fig. 5, we see that moving to the rλcdm and rsb casses does indeed fit the osciatory residuas better. Correspondingy, in Fig. 6, these casses aow even ower vaues of H 0 = 66.6±1.1 and 66.±1. km s 1 Mpc 1 respectivey, corresponding to the higher vaues of Ω c h, as shown in Fig.. Note aso that this enhanced freedom is aso associated with a higher vaue of τ and, hence, A s. Given the suppressed curvature power spectrum on arge scaes, it takes a arger τ to provide the same ampitude of poarization from reionization. In Fig. 7, we show that the combination of A s and Ω c h changes in the rλcdm mode from Tab. III aso enhances the ensing power spectrum. However, ensing reconstruction data do not favor such an enhancement. In Fig. 6 we see that, correspondingy, the preference for shifting H 0 even ower in the rλcdm mode mainy goes away with the +φφ dataset, where H 0 = 67. ± 1.1 km s 1 Mpc 1. C. Intermediate 00 TE residuas Ref. [9] showed that the TE data aone constrain H 0 to comparabe precision as the fu data, and that it favors the ow H 0 of the atter as opposed to the high H 0 of ( < 1000). We can attribute some of this abiity to the enhanced, degeneracy breaking, sensitivity to Ω c h in the TE and EE spectra around the first poarization trough at 00 (see Fig. 1). Given Panck noise, TE is more constraining than EE, and so we focus on its impact. In Fig. 8, we show the EE data set inferences on H 0 from the fu range and compare it to the < 1000 range. Most of the extra information from poarization comes from < 1000, as expected given Panck noise. Whie discarding the > 1000 range in the anaysis doubes the H 0 errors, the impact on the EE anaysis is notaby smaer: H 0 = 68.1 ± 0.83 km s 1 Mpc 1 compared with 67. ± 0.66 km s 1 Mpc 1 Moreover, the TE data at < 1000 pus H 0 coser to the vaue preferred by the dataset of the fu rather than < 1000 range. On the face of it, this provides independent support for ow H 0 that does not invove the mutipoe ranges known to have anomaies in and provides some proof against instrumenta systematics. However, even given the arge and sharp poarization responses to Ω c h in Fig. 1, the noisy TE Panck data seem unusuay sensitive to Ω c h and H 0. To better understand the sensitivity, in Fig. 9 we show the residuas of the TE data reative to the best fit ΛCDM mode and aso compare them to the best fit ΛCDM ( < 1000) mode. Note that neither cosmoogica mode is optimized to fit the poarization data. Poarized and unpoarized foregrounds are jointy reoptimized in the ΛCDM and SB cases, whereas for ΛCDM ( < 1000) they are hed fixed to their reoptimized ΛCDM vaues. The difference between the two modes is reativey sma compared with the measured errors. However, the data show a significant outier in the 165 bin reative to both modes at σ for the fu ΛCDM mode (ow H 0 ) and 3σ for the ΛCDM ( < 1000; high H 0 ) mode. We therefore attribute the enhanced sensitivity to H 0 in arge part to the presence of this one outier. To test this attribution, we aso show in Fig. 9 (bottom pane) the change in the EE ikeihood between the ΛCDM and ΛCDM ( < 1000) modes as a function of the maximum. The ΛCDM ( < 1000) mode jumps from a better fit to a worse fit across the = 165 bin, unike with the ikeihood (cf. Fig. 3). We aso show the best fit SB mode in Fig. 9. Intriguingy, even though this mode is not optimized to fit the outier, it does marginay reduce tension with the = 165 bin. It is aso interesting to note that even in the rλcdm where H 0 = 66.9 ± 0.69 km s 1 Mpc 1 and rsb casses where H 0 = 66.8 ± 0.7 km s 1 Mpc 1, the EE data set no onger aows the even ower H 0 vaues found from fitting the > 1000 osciatory residuas even with no ensing reconstruction data (see Fig. ). However, the Panck coaboration considers the TE and EE data sets as preiminary and so the anomaous sensitivity to H 0 in the ΛCDM context, whie intriguing, shoud not be overinterpreted.

9 9 1 C T E /σ 0 1 n L EE ΛCDM ( < 1000) SB FIG. 9. Top pane: residuas in the TE data (points) and the same modes as in Fig. 3 reative to the best fit ΛCDM mode of Tab. III (dashed ine). Bottom pane: n L EE of modes reative to ΛCDM in the same form as Fig. 3. Cosmoogica parameters are fixed to their best fit vaues and the ( < 1000) mode aso has the same foreground parameters as the mode incuding poarized foregrounds. The ΛCDM ( < 1000) mode becomes a worse fit to the EE dataset when passing the outier point at 165 (arrow), whereas the SB mode sighty improves the fit. ΛCDM ΛCDM ( < 1000) rλcdm rλcdm +φφ Parameter rλcdm rλcdm + φφ SB rsb n run n run,run m m m Posterior TABLE IV. Infationary parameters of best fit modes for various datasets. The running parameters are defined in Eq. and m 1, m, m 3 are are the first 3 principa components of the SB parameters rank ordered by increasing variance σ 8 Ω 1/ m FIG. 10. Ampitude of structure σ 8 Ωm 1/. Shifts in Ω ch from ΛCDM ( < 1000; red curve) to the fu range (back curve) drive σ 8 Ω 1/ m higher. Fitting the in-phase high- residuas with the enhanced freedom of the rλcdm mode drives it even higher (gray shaded), but in a manner inconsistent with the ensing reconstruction data (red shaded). Panck (fu) cosmoogy causes an increase in the amount of oca structure σ 8 Ω 1/ m (see Fig. 10). This increase is in moderate tension with measurements of the oca custer abundance, depending on the mass caibration empoyed [9, 39], and of weak gravitationa ensing [40, 41]. D. Loca ampitude σ 8Ω 1/ m The owering of H 0 and raising of Ω m between the WMAP9 or Panck ( < 1000) cosmoogy and the In Fig. 10, we show the effect of generaizing the mode cass to rλcdm. The aowed further raising of Ω c h and owering of H 0 with just the data set shifts σ 8 Ω 1/ m to higher vaues. However, once the φφ ensing reconstruction data are incuded, this preference disappears.

10 10 rλcdm rsb ( < 1000) EE ( < 1000) EE n s n run n run,run FIG. 11. Constraints on infationary parameters tit, its running, and the running of the running from Eq. for the 4 combinations of datasets in Tab. I rλcdm SB rsb 1.0 s [Mpc] C T T /σ pi /σi FIG. 1. Low- residuas (points) and modes reative to the best fit ΛCDM mode of Tab. III (back dashed). Whie the running mode rλcdm (green soid) can suppress power, the SB mode (bue dashed) accommodates the sharp dip favored by the data. Aowing both running and SB in the rsb mode (red dotted) eads to a simiar fit as SB SB i FIG. 13. SB parameters of the SB anaysis in units of their errors. No singe SB parameter p i is significanty detected due to covariances hidden in the marginaized errors. IV. INFLATIONARY FEATURES The same residuas in the Panck data that drive shifts in H 0 and σ 8 Ω 1/ m provide hints that infation may be more compex than singe-fied sow-ro infation woud impy. No singe set of ΛCDM parameters with a power-aw curvature power spectrum simutaneousy gives a very good fit to the data at ow and high, as shown in Fig. 3. Athough singe-fied infation aows for running of the tit and running of the running, the sow-ro approximation predicts that each order shoud be successivey suppressed by O(n s 1), i.e. n run,run How- ever, in the rλcdm anaysis n run,run = 9 ± 15 [9], a 1.9σ preference for a vioation of sow-ro infation (see Fig. 11). In Tab. IV, we give the best fit vaues for running parameters in the various mode casses, and in Fig. 1 we show how we these modes fit the data residuas at ow. Whie the rλcdm mode does indeed ower power on arge scaes, it fais to fit the sharp decine in power at As shown in Fig. 5, there is aso a sight preference for nonzero running parameters due to their abiity to accommodate cosmoogica parameter changes that fit the high- osciatory residuas. However, this preference does not remain if the ens reconstruction data are considered.

11 11.4 SB. SB EE m R e τ ΛCDM rλcdm 0.3 SB 3PC m k [Mpc 1 ] m m m 3 FIG. 16. Functiona reconstruction of infationary power spectrum ( R : 68%, 95% CL) from the first 3 principa components in the SB (shaded) anaysis compared with the best fit running mode rλcdm. FIG. 14. First 3 principa components of the SB p i covariance matrix in the (thin bue) and EE (thick red) anaysis (68%, 95% CL). Both m and m 3 show non-zero deviations whose significance is sighty enhanced by the addition of poarization information. δg s [Mpc] SB 3PC SB 3PC EE FIG. 15. Functiona reconstruction of generaized tit (δg : 68%, 95% CL) from the first 3 principa components in the SB (shaded) and SB EE (dashed and dotted ines) anayses. Poarization information sighty sharpens the dip near 400 Mpc. Given the arge preferred vaues of n run,run, such modes vioate sow ro and may in their own right be considered infationary features, so there is no theoretica motivation to mode the infationary curvature power spectrum as a truncated Tayor expansion around the pivot point as in Eq. (). Indeed, we see in Fig. 1 that the SB parameters of the generaized tit δg aow a sharper decine at 0 30, which better fits the data. On the other hand, with 0 free parameters, the SB mode tends to overfit the residuas eaving highy correated parameter errors and no singe p i with more than a 1σ deviation from the sow ro prediction of zero (see Fig. 13). From Fig. 3, we see that the net improvement in fitting the residuas with SB at < 30 is n L 1. We aso show that combining both SB and running parameters eads to fits that are very simiar to SB at ow. Correspondingy, once the SB parameters are marginaized, the significance of the preference for finite vaues decreases to n run,run = 35 ± 7 (see Fig. 11). Its best fit vaue 41 actuay increases in rsb, in part due to its abiity to fit the high- in-phase residuas. In order to expose the most significant inear combination of the SB p i parameters and their impications, we examine the m a principa components (PCs) of its covariance matrix C ij = p i p j, rank ordered by increasing eigenvaue or variance, foowing Ref. [5], as detaied in the Appendix. For the SB anaysis, the best constrained PCs aso contain the most significant deviations from sow ro. Note that m = m 3 = 0 is outside the 95% CL region and is in fact excuded at 99% CL in Fig. 14 with m = 0.7 ± 0.11 and m 3 = 0.5 ± Of course one shoud interpret this as a oca significance given the number of we measured PCs that this anomay coud have appeared in. Furthermore, using the SB PC basis, we can assess whether the additiona poarization information in the EE data set supports or weakens this preference. In Fig. 14 we see that poarization enhances the significance of the deviation to σ in m 3 = 0.8 ± 0.14 whie eaving m = 0.7 ± 0.11 essentiay unchanged at.5σ and improving the bounds on m 1. Retaining (or fitering ) the SB to ony the first 3 PCs ( SB 3PC, see Eq. A17), we can better isoate the most significant impications for infation and the curvature power spectrum. In Fig. 15, we show the infationary

12 1 impications for δg, the anaog of the tit for modes with rapid deviations from sow ro. The 3PC preference is for a > 95% CL sharp suppression in δg beginning at 300 Mpc and preceded by a ess significant enhancement. Whie the additiona poarization information in EE does not further enhance the ampitude or significance of the suppression, it does enhance the sharpness. Finay, in Fig. 16 we show the impications for the curvature power spectrum. Unike the smooth deviations of the rλcdm mode, the SB 3PC mode has a sharp suppression around k 04 Mpc 1. In canonica singe-fied infation, such a suppression woud be associated with an abrupt transition from a faster to a sower roing rate. V. DISCUSSION We have expored the interpay between features from infation, features in the Panck 015 data, and shifts in the cosmoogica parameters. Preference for high H 0 vaues in ΛCDM from the data at < 1000 is driven by ow < 40 residuas in the data, which then prefer the acoustic peaks to be raised by enhancing radiation driving through owering Ω c h, and consequenty raising H 0 to match the anguar position of the peaks. However, these cosmoogica parameter variations do not match the residuas particuary we. If we marginaize their impact on H 0 by fitting, instead, to infationary features of the generaized tit δg form in the SB anaysis, then the ow H 0 favored by the fu data is compatibe with even the < 1000 data. Correspondingy, the SB parameters show principa components that jointy deviate from zero at the 99% CL in oca significance. These deviations in δg prefer a sharp suppression of power around k 04 Mpc 1. The residuas can aso be fit to a running of the running of the tit in the rλcdm mode, but those fits cannot reproduce the sharpness of the suppression. Furthermore, the arge running of the running favored woud itsef represent an infationary feature in the power spectrum that cannot be expained in sow-ro infation. The same > 1000 osciatory residuas which drive H 0 ower and σ 8 Ω 1/ m higher in the fu data set prefer an even ower H 0 if ony > 1000 is considered [13], but this does not provide a good goba fit in ΛCDM with power-aw initia conditions [14]. With running parameters in the tit, H 0 = 66.6 ± 1.1 km s 1 Mpc 1. The combination of ow- and high- anomaies makes running of the running preferred at the 1.9σ eve in the data set. However, on the ow side, once the sharp feature is marginaized with SB this preference drops to the 1.3σ eve. On the high side, once ensing reconstruction is taken into account, preference for fitting osciatory residuas disappears and H 0 = 67.9 ± 0.9 km s 1 Mpc 1 returns to neary the ΛCDM vaue. In poarization, the signature of changes in H 0 through Ω c h is much sharper and somewhat arger. Poarization suffers ess projection effects and acks the eary ISW that broadens the signature in temperature and its measurement at 00 is important for assuring the robustness of H 0 inferences from the CMB in ΛCDM. Indeed, despite being noisier than, Panck TE data aone aready constrain H 0 as we as [9]. They aso are consistent with the ow H 0 constraints from the fu data, even though their impact comes from < On the other hand, the Panck TE data are anomaousy sensitive to H 0 due to an outier at 165, which even more strongy disfavors the high H 0 of the < 1000 data fit. Poarization data are aso beginning to hep constraints on the infationary features favored by the ow- residuas. In Panck 015, ony the 30 HFI data was reeased, and so the main effect is to strengthen the confinement of the features to arge scaes, thus favoring a sharper feature in the joint anaysis. This sharpening in turn enhances the significance of the deviations in one of the principa components to σ eaving the other at.5σ. At < 30, ony LFI poarization data was reeased, and the impact of infationary features is to change inferences on the optica depth τ, since with smaer infationary power a higher τ is required to produce the same poarization power. In the fina Panck reease, the incusion of HFI poarization data at < 30 shoud aso provide insight on the ow- residuas and their impications for infationary features. It wi be important then to distinguish infationary features from reionization features beyond the simpe near instantaneous reionization modes used here [4]. We eave these considerations to a future work. ACKNOWLEDGMENTS WH and CH were supported by grants NSF PHY , NSF PHY , U.S. Dept. of Energy contract DE-FG0-13ER41958, and NASA ATP NNX15AKG. Computing resources were provided by the University of Chicago Research Computing Center. VM was supported in part by the Chares E. Kaufman Foundation, a supporting organization of the Pittsburgh Foundation. Appendix A: GSR Spine Basis The GSR approach [31] is a generaization of the sowro approximation that aows for features in the power spectrum of n k 1 in a genera singe-fied mode that incudes canonica [3], P (X, φ) [43], and Horndeski (Gaieon) infation [44]. In this approximation, the curvature power spectrum R is sourced by a singe function f = πz c s s, where z is the Mukhanov variabe in the Mukhanov-Sasaki equation, c s is the sound speed and s

13 13 is the sound horizon: s = aend a da a c s ah. (A1) Here a end is the scae factor at the end of infation and H is the Hubbe parameter during infation. For a canonica scaar fied in Panck units φ, z = a dφ/da. In the ordinary sow-ro approximation, R f, whereas in GSR: where n R(k) G(n s ) + s + n [ 1 + I 1 (k) ], ds s W (ks)g (n s) G(n s) = n f + 3 (n f), (A) (A3) with = d/d n s and s is an arbitrary epoch during infation such that a reevant k-modes are we outside the sound horizon, ks 1. The window functions, W (x) = 3 sin(x) x 3 at eading order and at second order through 3 cos(x) x X(x) = 3 x 3 (sin x x cos x), I 1 (k) = sin(x), (A4) x ds s G (n s)x(ks), (A5) (A6) characterize the freezeout of the source function around sound horizon crossing. This form for the power spectrum remains a good approximation if the second order term [3] I 1 < 1, (A7) and hence aows for up to order unity features in the curvature power spectrum. The ordinary sow-ro approximation corresponds to G (n s) = 1 n s, and resuts in a power-aw curvature power spectrum. In the main text, we use GSR to fit the ow- anomaies in the power spectrum, and so we choose to restrict our parameterization of fuctuations, δg, around this constant to 00 < s Mpc < (A8) Next, we foow Ref. [3] in defining a band imit for the frequency of deviations by samping p i = δg (n s i ) at a rate of 10 per decade for a tota of 0 parameters. This rate is sufficient to capture the ow- anomaies. We then construct the smooth function using the natura spine of these points. Specificay, we expoit the inearity of spines to define the spine basis (SB) functions B i (n s) from the natura spine of the set of unit ampitude perturbations p i = 1, p j i = 0. In the SB cass of modes we then describe the source function as: G (n s) = (1 n s ) + i p i B i (n s). (A9) The advantage of the GSR form in Eq. (A) is that the integras are inear in G and, hence, the impact of the individua components can be precomputed separatey as: W i (k) = X i (k) = s 0 ds s W (ks)b i(n s), ds s X(ks)B i(n s), (A10) so that the power spectrum becomes a sum over the basis ( ) ns 1 k n R(k) = n A s + [ p i Wi (k) W i (k 0 ) ] k 0 i [ ] 1 + I + n 1 (k) 1 + I1 (k, (A11) 0) where I 1 (k) = π (1 n s) + 1 p i X i (k). i (A1) Note that we have absorbed the normaization constant G(n s ) into the ampitude of the power spectrum at the scae k 0 : A s = R(k 0 ). (A13) In the rsb cass of modes we repace the first term in Eq. (A11) with the running form defined by Eq. (). Since our choice of parameters oversampes δg reative to what the data can constrain, individua measurements of p i from the MCMC mainy fit noise with any true signa buried in the covariance between parameters. We therefore construct the principa components derived from an eigenvaue decomposition of the MCMC covariance matrix estimate: C ij = p i p j p i p j = a S ia σ as ja, (A14) where S ia is an orthonorma matrix of eigenvectors. Specificay, we define the PC parameters m a = i S ia p i, such that their covariance matrix satisfies: (A15) m a m b m a m b = δ ab σ a. (A16)

14 Given a rank ordering of the PC modes from smaest to argest variance, we can aso construct a 3 PC fitered reconstruction of δg as in Ref. [3]: 3 δg 3PC(n s i ) = m a S ia, a=1 14 (A17) and simiary use δg 3PC to construct the 3PC fitered used in the main text. R [1] C. L. Bennett, A. Banday, K. M. Gorski, G. Hinshaw, P. Jackson, P. Keegstra, A. Kogut, G. F. Smoot, D. T. Wikinson, and E. L. Wright, Astrophys. J. 464, L1 (1996), arxiv:astro-ph/ [astro-ph]. [] D. N. Sperge et a. (WMAP), Astrophys. J. Supp. 148, 175 (003), arxiv:astro-ph/03009 [astro-ph]. [3] G. Hinshaw et a. (WMAP), Astrophys. J. Supp. 148, 135 (003), arxiv:astro-ph/03017 [astro-ph]. [4] H. V. Peiris et a. (WMAP), Astrophys. J. Supp. 148, 13 (003), arxiv:astro-ph/0305 [astro-ph]. [5] M. J. Mortonson, C. Dvorkin, H. V. Peiris, and W. Hu, Phys. Rev. D79, (009), arxiv: [astroph.co]. [6] P. A. R. Ade et a. (Panck), Astron. Astrophys. 571, A15 (014), arxiv: [astro-ph.co]. [7] P. A. R. Ade et a. (Panck), Astron. Astrophys. 571, A16 (014), arxiv: [astro-ph.co]. [8] Y.-F. Cai, E. G. M. Ferreira, B. Hu, and J. Quintin, Phys. Rev. D9, (015), arxiv: [astroph.co]. [9] P. A. R. Ade et a. (Panck), Astron. Astrophys. 594, A13 (016), arxiv: [astro-ph.co]. [10] A. G. Riess, L. Macri, S. Casertano, H. Lampeit, H. C. Ferguson, A. V. Fiippenko, S. W. Jha, W. Li, and R. Chornock, Astrophys. J. 730, 119 (011), arxiv: [11] D. Larson, J. L. Weiand, G. Hinshaw, and C. L. Bennett, Astrophys. J. 801, 9 (015), arxiv: [astroph.co]. [1] D. N. Sperge, R. Fauger, and R. Hozek, Phys. Rev. D91, (015), arxiv: [astro-ph.co]. [13] G. E. Addison, Y. Huang, D. J. Watts, C. L. Bennett, M. Hapern, G. Hinshaw, and J. L. Weiand, Astrophys. J. 818, 13 (016), arxiv: [astro-ph.co]. [14] N. Aghanim et a. (Panck), (016), arxiv: [astro-ph.co]. [15] A. Shafieoo and D. K. Hazra, JCAP 1704, 01 (017), arxiv: [astro-ph.co]. [16] W. Hu and T. Okamoto, Phys. Rev. D69, (004), arxiv:astro-ph/ [astro-ph]. [17] S. L. Bride, A. M. Lewis, J. Weer, and G. Efstathiou, Mon. Not. Roy. Astron. Soc. 34, L7 (003), arxiv:astro-ph/ [astro-ph]. [18] S. M. Leach, Mon. Not. Roy. Astron. Soc. 37, 646 (006), arxiv:astro-ph/ [astro-ph]. [19] C. Seafon, L. Verde, and R. Jimenez, Phys. Rev. D7, (005), arxiv:astro-ph/ [astro-ph]. [0] C. Dvorkin, H. V. Peiris, and W. Hu, Phys. Rev. D77, (008), arxiv: [astro-ph]. [1] H. V. Peiris and L. Verde, Phys. Rev. D 81, 0130 (010), arxiv: [astro-ph.co]. [] P. Paykari and A. H. Jaffe, apj 711, 1 (010), arxiv: [astro-ph.co]. [3] C. Dvorkin and W. Hu, Phys. Rev. D84, (011), arxiv: [astro-ph.co]. [4] P. A. R. Ade et a. (Panck), Astron. Astrophys. 571, A (014), arxiv: [astro-ph.co]. [5] V. Miranda, W. Hu, and C. Dvorkin, Physica Review D 91, (015). [6] P. A. R. Ade et a. (Panck), Astron. Astrophys. 594, A0 (016), arxiv: [astro-ph.co]. [7] J. Sik and M. S. Turner, Phys. Rev. D35, 419 (1987). [8] D. Poarski and A. A. Starobinsky, Nuc. Phys. B385, 63 (199). [9] J. A. Adams, G. G. Ross, and S. Sarkar, Nuc. Phys. B503, 405 (1997), arxiv:hep-ph/ [hep-ph]. [30] P. Hunt and S. Sarkar, Phys. Rev. D70, (004), arxiv:astro-ph/ [astro-ph]. [31] E. D. Stewart, Phys. Rev. D65, (00), arxiv:astro-ph/01103 [astro-ph]. [3] C. Dvorkin and W. Hu, Phys. Rev. D81, (010), arxiv: [astro-ph.co]. [33] C. Dvorkin and W. Hu, Phys. Rev. D8, (010), arxiv: [astro-ph.co]. [34] A. Lewis and S. Bride, Phys. Rev. D66, (00), arxiv:astro-ph/ [astro-ph]. [35] A. Lewis, A. Chainor, and A. Lasenby, Astrophys. J. 538, 473 (000), arxiv:astro-ph/ [astro-ph]. [36] W. Hu and N. Sugiyama, Astrophys. J. 444, 489 (1995), arxiv:astro-ph/ [astro-ph]. [37] W. Hu and S. Dodeson, Ann. Rev. Astron. Astrophys. 40, 171 (00), arxiv:astro-ph/ [astro-ph]. [38] W. Hu, M. Fukugita, M. Zadarriaga, and M. Tegmark, Astrophys. J. 549, 669 (001), arxiv:astro-ph/ [astro-ph]. [39] P. A. R. Ade et a. (Panck), Astron. Astrophys. 571, A0 (014), arxiv: [astro-ph.co]. [40] S. Joudaki et a., Mon. Not. Roy. Astron. Soc. 465, 033 (017), arxiv: [astro-ph.co]. [41] H. Hidebrandt et a., Mon. Not. Roy. Astron. Soc. 465, 1454 (017), arxiv: [astro-ph.co]. [4] C. H. Heinrich, V. Miranda, and W. Hu, Phys. Rev. D95, (017), arxiv: [astro-ph.co]. [43] W. Hu, Phys. Rev. D84, (011), arxiv: [astro-ph.co]. [44] H. Motohashi and W. Hu, (017), arxiv: [hep-th].