The Knotted Sky II: Does BICEP2 require a nontrivial primordial power spectrum?
Abstract
An inflationary gravitational wave background consistent with BICEP2 is difficult to reconcile with a simple powerlaw spectrum of primordial scalar perturbations. Tensor modes contribute to the temperature anisotropies at multipoles with , and this effect — together with a prior on the form of the scalar perturbations — was the source of previous bounds on the tensortoscalar ratio. We compute Bayesian evidence for combined fits to BICEP2 and Planck for three nontrivial primordial spectra: a) a running spectral index, b) a cutoff at fixed wavenumber, and c) a spectrum described by a linear spline with a single internal knot. We find no evidence for a cutoff, weak evidence for a running index, and significant evidence for a “broken” spectrum. Taken at facevalue, the BICEP2 results require two new inflationary parameters in order to describe both the broken scale invariance in the perturbation spectrum and the observed tensortoscalar ratio. Alternatively, this tension may be resolved by additional data and more detailed analyses.
a]Kevork N. Abazajian, b]Grigor Aslanyan, b]Richard Easther, b]and Layne C. Price
Prepared for submission to JCAP
The Knotted Sky II: Does BICEP2 require a nontrivial primordial power spectrum?

Department of Physics, University of California at Irvine, Irvine, CA 92697

Department of Physics, University of Auckland, Private Bag 92019, Auckland, New Zealand
Contents
1 Introduction
The BICEP2 experiment [1, 2] has reported a detection of primordial Bmodes in the cosmic microwave background (CMB).^{1}^{1}1Bmode polarization from CMB lensing was detected earlier by the POLARBEAR experiment [3]. The most natural explanation for a Bmode signal is a stochastic background of longwavelength gravitational waves, or tensor perturbations [4]. This constitutes strong prima facie evidence for an inflationary phase in the early universe, the most widelystudied source for a stochastic background of gravitational waves. If the observational data and theoretical explanation are confirmed, the Bmode signal will provide unprecedented insight into the mechanism responsible for inflation.
The measured tensortoscalar ratio has a confidenceinterval (CI) of and differs from zero with a statistical significance of . However, the temperature data from Planck [5, 6, 7], SPT [8], and ACT [9], combined with WMAP [10] polarization, yields at the CI, in significant tension with the BICEP2 result.
There are several potential explanations for this discrepancy. The first is that the BICEP2 analysis overestimates the amplitude of the Bmode itself [11]. The second possibility is that the primordial Bmode is accurately measured, but sourced by a mechanism unrelated to the standard assumptions for the primordial inflationary phase [12, 13, 14, 15, 16, 17, 18, 19, 20, 21]. Conversely, existing CMB data may have been misanalysed, although this appears unlikely given the agreement of Planck with WMAP at large and intermediate scales and with ACT and SPT at small scales.
Another suggestion is that the largescale scalar power spectrum is suppressed relative to that predicted by the bestfit CDM scenario. PreBICEP2 constraints on the inflationary gravitational wave background were driven primarily by the contribution of tensor modes to the temperaturetemperature (TT) anisotropies. This is illustrated in Fig. 1, which shows the contribution to at low from a tensor background with . The tensor contribution to individual is small, but systematically increases the TT multipoles for all . However, to constrain the primordial tensor background using this signal we must have an independent estimate of the contributions from scalar perturbations alone. Moreover, the measured lowmultipole typically lie below the best fit values for simple powerlaw spectra, so the inclusion of a tensor background is likely to reduce the likelihood relative to . Fig. 1 also shows a sample angular power spectrum derived from a primordial scalar spectrum with a sharp cutoff in power at a comoving wavenumber . This particular scenario has a very low likelihood relative to the Planck and WMAP datasets, but provides an extreme illustration of how a scalar spectrum with a cutoff could compensate for a tensor contribution to the for .
In this paper we focus on the implications of the BICEP2 result for the scalar power spectrum , performing joint analyses of the BICEP2 and Planck datasets. We consider three possibilities: (i) a running spectral index, (ii) a sharp cutoff in power at scale , and (iii) a discontinuity in the spectral index at scale . The latter two scenarios are implemented via the algorithm described by us in Ref. [22].
The BICEP2 analysis [1, 2] presents joint constraints from BICEP2 and Planck with a running index, but focusses primarily on the polarization and Bmode amplitude and does not discuss the issue in detail. We reproduce the BICEP2 constraints on a running index, and compute Bayesian evidence (relative to CDM) of for the running case. The cutoff spectrum does not give a significant improvement, since it suppresses the scalar power by a factor far larger than the corresponding increase in power due to tensor contributions. Finally, a break in the spectral index — implemented as a splined with a single “interior” point at an arbitrary amplitude and location — provides the best fit to the data. With a broad, uninformative prior we find the change in the logarithmic evidence ratio is . However, using a prior that includes information from our investigation of the Planck data alone, which disfavors knots at or dramatic changes in amplitude, we infer an evidence ratio of , a “significant” to “strong” detection according to conventional model selection criteria.
Consequently, this analysis would suggest that BICEP2 has actually made three significant discoveries about inflation: the first is to confirm its existence, and the second is to show that it took place at a relatively high energy scale. Thirdly, these results also suggest that the primordial scalar spectrum has a nontrivial structure which is inconsistent with simple models of inflation. Alternatively, the tensor spectrum may differ from the “standard” inflationary form, due to either a variant model of inflation [23] or a mechanism that is independent of inflation [12, 13, 14, 17, 20]. Of course, the conservative explanation of our findings is that they point to tension between the Planck and BICEP2 datasets which will be resolved by more complete analyses and/or additional data.
2 Method
2.1 Likelihoods, priors, and Bayesian evidence
We combine the mode results from BICEP2 [1, 2] with temperature and lensing data from Planck [5]. We use the Planck likelihood code [24] with Commander, CamSpec, and lensing likelihood files for the datalikelihood evaluation. We use the COSMO++ library [25] to combine the modified form of the primordial scalar power spectrum (Fig. 2) with the Planck likelihood code and calculate the CMB angular power spectra with CLASS [26, 27]. We employ multimodal nested sampling for parameter estimation and the computation of evidence, using the publicly available code MultiNest [28, 29, 30].
We use model posterior probabilities and Bayesian evidences to compare the statistical significance of two competing reconstruction models. This approach penalizes models with parameters for which the likelihood is large only in small regions of parameter space and protects against overfitting the scalar spectrum with too many knots or bins. This approach has previously been employed for power spectrum reconstruction [31, 32, 33, 34, 35, 36] and gives conservative and robust assessments of possible physical features in the data.
We use the posterior probability to assess the statistical significance of a model . If two models and have the same prior probability, Bayes’ theorem yields the relative betting odds between the models via the Bayes factor
(2.1) 
where the Bayesian evidence (marginalized likelihood) is
(2.2) 
for the model parameters . Here, is the datalikelihood and is the parameter prior probability. Equation (2.2) is a function of our choice of prior, so to be conservative we allow the parameters that define the power spectrum features to vary over a wide range of values in order to thoroughly search the parameter space [37]. Because there is little likelihood for features at , a model with more parameters must produce a very large improvement in the likelihood relative to a featureless scenario for evidence to yield “betting odds” that strongly support the more complex model. Consequently, we also examine the improvement in fit () at the maximum likelihood point for each scenario.
(Posterior Odds)  Jeffreys Scale  Cosmology Scale 

0.0 to 1.0  Not worth more than a bare mention  
1.0 to 2.5  Substantial  Weak 
2.5 to 5.0  Strong  Significant 
Decisive  Strong 
We use uniform priors for all of the standard cosmological parameters , , , and , as well as the “nuisance” parameters in the CamSpec likelihood. For convenience, we set the prior probability distribution for to the posterior distribution obtained by BICEP2 [1, 2], which is equivalent to the direct evaluation of the underlying likelihood with a uniform prior on . Consistent with BICEP2, we assume a flat tensor power spectrum and a pivot scale of [1, 40]. We check that we recover the joint constraints on and and the marginalised posterior for the running index reported by BICEP.
Broad Priors  

Linear Spline  Cutoff  Running 
Informative Priors  

Linear Spline  Cutoff 
2.2 Non–powerlaw scalar spectrum
We perform a complete marginalization of the Planck temperature and lensing likelihood, varying the remaining CDM variables and the full set of Planck “nuisance parameters.” We calculate the Bayesian evidence ratios and combine this with the parameter posterior probabilities to give a conservative modelselection guideline. We use a scaleinvariant primordial tensor spectrum and draw from a prior defined by the BICEP2 posterior, as noted previously.
We summarize all power spectrum priors in Table 2. We reconstruct the primordial scalar power spectrum in the range , using a generalization of the “knotspline” procedure, developed in Refs [31, 32, 41, 33, 42, 43, 34, 35, 44, 22]. This process is illustrated in Fig. 2. A complete discussion is given in Ref. [22], but it can be summarised as follows:

Fix the endpoints at and , but allow their amplitudes and to vary, with logarithmic prior, in the ranges , where .

Add a “knot” with logarithmic prior in , between and allow its amplitude to vary in the same ranges as the endpoints in Step 1.

Interpolate between the endpoints and the knot with linearspline interpolation.
Varying the knotlocation compensates for the “lookelsewhere” effect, as the knot can move over the whole range of . This permits the reconstruction of features in the scalar primordial power spectrum that could appear at any scale in . The broad ranges on the priors for and indicate that possible features are not restricted to large scales. We call this the “broad” prior on the knot’s position and amplitude and it gives conservative values for the Bayesian evidence.
We can also use the information gained during our Planckonly analysis [22] to update the prior on the knot position and amplitude. The Planck data indicates that features should only appear at large scales with . This reduces the (logarithmic) range of the knot position to twothirds of the original volume. Furthermore, we know the spectral index is red at large scales, and a low value of will generate a blue spectrum at , so we can further stipulate that . This gives the Informative Prior for the broken spectrum, which yields less conservative and more significant Bayesian evidences. The Planck likelihood is almost zero in the excluded regions, so the evidence for the informative prior relative to the broad prior is scaled by the ratio of the relative parameter space volume, giving an increase of . In what follows we report evidence values with both priors.
We also analyze the standard powerlaw spectra with both (a) a sharp cutoff at and (b) a running spectral index, defined by
(2.3) 
where is the running. A running spectral index was considered in the BICEP2 analysis [1]; we repeat this analysis both as a check on our inclusion of the BICEP2 results in our likelihood and in order to calculate Bayesian evidence for the running index.
These two forms of the power spectrum are also illustrated in Fig. 2. The pivot scale for the power law is . The cutoff prior is from a logarithmic prior, and we have a uniform prior on with the range . Note that when the running becomes large the Taylor expansion in Eq. (2.3) can be an inaccurate parametrization for the inflationary power spectrum [45]. As above, we also use an Informative Prior for the cutoff spectrum based on the analysis in Ref. [22] with . We report the evidences for both priors, although this makes little difference to the conclusions.
3 Results
Model  

No Knots  —  —  — 
Knot  
Model  
—  —  —  
Cutoff  
Running  — 
CDM  No knots  knot  Cutoff  Running  

—  —  
—  —  
—  —  —  —  
—  —  —  —  
—  —  —  
—  —  —  
—  —  —  —  
—  —  —  — 
Figure 3 shows the reconstructed scalar power spectrum with a standard powerlaw ( knots); a powerlaw with a sharp cutoff; a powerlaw with a running spectral index; and a knot linearspline model. The Bayesian evidences are given in Table 3, along with the improvements in best fit improvements for each case. Although the linear spline with no knots is equivalent to the standard power law case without running, the different parameterisations lead to different prior volumes. Consequently, we report the Bayesian evidence for the cutoff and running power spectra relative to the standard power law case with common priors for and , while we report the Bayesian evidence for the 1knot model compared to the 0knot powerlaw model.
In Figure 4 we show the posterior distributions for obtained with these models, along with the result obtained with the BICEP2 data alone. The posterior distribution for derived for the running power spectrum is shown in Figure 5. Constraints for the standard cosmological parameters and the power spectrum parameters are displayed in Table 4.^{2}^{2}2Although the no knots model is completely equivalent to the standard CDM model in its functional form, the parametrizations and subsequently the priors are different. For this reason the resulting posteriors on all of the parameters could be slightly different, and we present the results for both cases. As noted previously while the spline model with no knots is equivalent to CDM, the different priors associated with the power spectrum parameterisations lead to small differences in ranges for the standard cosmological parameters. Our CI bounds for agree well with those reported in the BICEP2 analysis () [1].
The posteriors for the reconstructed power spectrum in Fig. 3 all recover the standard powerlaw form at small scales . This confirms the Planckonly analysis of Ref. [22], indicating that the BICEP2 detection of does not further imply power spectrum features at intermediate to small scales.
At larger scales () all the nonminimal models indicate a suppression of power in the spectrum of scalar perturbations. The black lines in Fig. 3 are the most likely power spectra, and these all decrease at small . The increase in the likelihood for the bestfit power spectra are reported in Table 3, although we caution that cosmic variance is important at these scales. Also, while a local feature at could also yield posteriors with largescale power suppression (as shown in Section 4 of Ref. [22]), we can be more certain about the posteriors in Fig. 3: since the tensor contribution to is nearly uniform for scales , offsetting this increase in power should require a compensating decrease at all scales and not a local feature at intermediate scales.
The Bayesian evidences in Table 3 show some support for the 1knot linearspline model and a running spectrum with and , respectively. With the Planck temperature and lensing data alone, the evidence for the 1knot model is only [22], indicating that BICEP2 data contributes significantly to the increased evidence. While these models give qualitatively similar spectra, the running powerlaw requires that there is less power at scales in order to achieve the same suppression of largescale power as the 1knot model. The at these scales are wellexplained by scalar contributions only, but the increased likelihood due to largescale suppression is partially offset by the Planck likelihood at . This increases the Bayesian evidence for the 1knot model compared to the running. The cutoff model gives little overall improvement to the data, since a cutoff causes a large decrease in the over all scales larger than the cutoff scale, as seen in Fig. 1. The 1knot model generalizes the cutoff at large scales, and thus gives higher evidence values.
Overall, the Bayesian evidences with the broad, uninformative priors show only a mild increase over the powerlaw prediction. The evidence computed for the 1knot model can be characterized as either “weak” or “substantial,” depending whether one uses the Jeffreys or cosmology scale, as described in Table 1. However, the Informative Prior obtained by incorporating the insight gained from our analysis of the Planck data on its own has a significantly higher evidence due to the reduced parameter volume. This increases by , and the resulting evidence of 3.1 constitutes “strong” or “significant” evidence for the suppression of power in the spectrum of primordial perturbations at large scales.
4 Discussion
Using data from Planck and BICEP2 we perform parameter estimates and calculate Bayesian evidence in order to explore the implications of the BICEP2 result for different parameterisations of the primordial scalar power spectrum. In agreement with the BICEP2 analysis we find that the posterior probability for a running spectral index excludes zero at the 95% confidence interval. Similarly, a spectrum defined by a linear spline with one internal knot shows a distinct preference for a suppression of power in the scalar spectrum at large angular scales, . For an informative prior incorporating the results of the Planckonly analysis [22] the corresponding Bayesian model selection criteria show a pronounced preference for a model with broken scale invariance, with .
This paper extends the analysis of Ref. [22], which reconstructs the primordial scalar power spectrum from Planck temperature data alone. In particular, Ref. [22] shows that the evidence for extra structure in the scalar power spectrum is negligible and that this reconstruction technique successfully recovers artificial signals injected into simulated temperature maps. Consequently, these results quantify the impact of the additional information provided by BICEP2.
If the scalar power spectrum is not well described by the usual powerlaw form, the estimated values of other novel cosmological parameters may also be modified, including sum of the neutrino masses or the number of effective relativistic degrees of freedom, which can be degenerate with features in the scalar power spectrum [46, 47, 48, 49, 50, 44]. Likewise, the BICEP2 data does not put robust constraints on the tensor index . Generically, the expectation from inflation is that to first order in slowroll, so we assume a scaleinvariant tensor spectrum in this analysis. However, while a blue (or sharply peaked) tensor background would also alleviate the tension between Planck and BICEP2, the physical processes that generated this spectrum would at least be as radical as considering a non–powerlaw scalar spectrum.
The recent BICEP2 result provides strong evidence for a primordial tensor background, from which it is inferred that the very early universe underwent an inflationary phase. However, the results presented here imply that BICEP2 also suggests that this inflationary phase yields a nontrivial scalar power spectrum, and that the underlying inflationary mechanism is not welldescribed by a simple, smooth singlefield potential. The inverse problem associated with reconstructing the inflationary potential from data has been widely discussed [51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69] and these methods would have at least three nontrivial input parameters in such a scenario. Likewise, with a large negative running similar to the central value found here, simple inflationary models typically yield an unacceptably small number of efoldings, implying that the running itself must be scale dependent [70].
Needless to say, this analysis takes both the current Planck and BICEP2 dataproducts at facevalue. The most conservative explanation for these result is that future analyses will eliminate much of the apparent tension between BICEP2 and other cosmological datasets. From this perspective our analysis quantifies the extent of that tension.
Acknowledgments
We acknowledge the contribution of the NeSI highperformance computing facilities and the staff at the Centre for eResearch at the University of Auckland, especially Mark Gahegan and Gene Soudlenkov. New Zealand’s national facilities are provided by the New Zealand eScience Infrastructure (NeSI) and funded jointly by NeSI’s collaborator institutions and through the Ministry of Business, Innovation & Employment’s Research Infrastructure programme [http://www.nesi.org.nz]. KNA is supported by NSF CAREER Grant No. PHY1159224.
References
 [1] BICEP2 Collaboration, P. A. R. Ade et. al., BICEP2 I: Detection Of Bmode Polarization at Degree Angular Scales, arXiv:1403.3985 astroph.CO (2014).
 [2] BICEP2 Collaboration, P. A. R. Ade et. al., BICEP2 II: Experiment and ThreeYear Data Set, arXiv:1403.4302 astroph.CO (2014).
 [3] The POLARBEAR Collaboration, P. Ade et. al., A Measurement of the Cosmic Microwave Background BMode Polarization Power Spectrum at SubDegree Scales with POLARBEAR, arXiv:1403.2369.
 [4] U. Seljak and M. Zaldarriaga, Signature of gravity waves in polarization of the microwave background, Phys.Rev.Lett. 78 (1997) 2054–2057, [astroph/9609169].
 [5] Planck Collaboration, P. Ade et. al., Planck 2013 results. I. Overview of products and scientific results, arXiv:1303.5062.
 [6] Planck Collaboration, P. Ade et. al., Planck 2013 results. XVI. Cosmological parameters, arXiv:1303.5076.
 [7] Planck Collaboration, P. Ade et. al., Planck 2013 results. XXII. Constraints on inflation, arXiv:1303.5082.
 [8] SPT Collaboration, Z. Hou et. al., Constraints on Cosmology from the Cosmic Microwave Background Power Spectrum of the 2500square degree SPTSZ Survey, arXiv (Dec., 2012) [arXiv:1212.6267].
 [9] ACT Collaboration, S. Das et. al., The Atacama Cosmology Telescope: Temperature and Gravitational Lensing Power Spectrum Measurements from Three Seasons of Data, arXiv (Jan., 2013) [arXiv:1301.1037].
 [10] WMAP Collaboration, G. Hinshaw et. al., NineYear Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Parameter Results, Astrophys.J.Suppl. 208 (2013) 19, [arXiv:1212.5226].
 [11] M. J. Mortonson and U. Seljak, A joint analysis of Planck and BICEP2 B modes including dust polarization uncertainty, arXiv:1405.5857.
 [12] U. Seljak, U.L. Pen, and N. Turok, Polarization of the microwave background in defect models, Phys.Rev.Lett. 79 (1997) 1615–1618, [astroph/9704231].
 [13] L. Pogosian and M. Wyman, Bmodes from cosmic strings, Phys.Rev. D77 (2008) 083509, [arXiv:0711.0747].
 [14] K. JonesSmith, L. M. Krauss, and H. Mathur, A Nearly Scale Invariant Spectrum of Gravitational Radiation from Global Phase Transitions, Phys.Rev.Lett. 100 (2008) 131302, [arXiv:0712.0778].
 [15] T. Kobayashi and F. Takahashi, Running Spectral Index from Inflation with Modulations, JCAP 1101 (2011) 026, [arXiv:1011.3988].
 [16] J. L. Cook and L. Sorbo, Particle production during inflation and gravitational waves detectable by groundbased interferometers, Phys.Rev. D85 (2012) 023534, [arXiv:1109.0022].
 [17] L. Senatore, E. Silverstein, and M. Zaldarriaga, New Sources of Gravitational Waves during Inflation, arXiv:1109.0542.
 [18] N. Barnaby, J. Moxon, R. Namba, M. Peloso, G. Shiu, et. al., Gravity waves and nonGaussian features from particle production in a sector gravitationally coupled to the inflaton, Phys.Rev. D86 (2012) 103508, [arXiv:1206.6117].
 [19] D. Carney, W. Fischler, E. D. Kovetz, D. Lorshbough, and S. Paban, Rapid field excursions and the inflationary tensor spectrum, JHEP 1211 (2012) 042, [arXiv:1209.3848].
 [20] C. R. Contaldi, M. Peloso, and L. Sorbo, Suppressing the impact of a high tensortoscalar ratio on the temperature anisotropies, arXiv:1403.4596.
 [21] M. Czerny, T. Kobayashi, and F. Takahashi, Running Spectral Index from Largefield Inflation with Modulations Revisited, arXiv:1403.4589.
 [22] G. Aslanyan, L. C. Price, K. N. Abazajian, and R. Easther, The Knotted Sky I: Planck constraints on the primordial power spectrum, arXiv:1403.5849.
 [23] V. Miranda, W. Hu, and P. Adshead, Steps to Reconcile Inflationary Tensor and Scalar Spectra, arXiv:1403.5231.
 [24] P. Collaboration, Planck 2013 results. XV. CMB power spectra and likelihood, arXiv:1303.5075 (2013).
 [25] G. Aslanyan, Cosmo++: An ObjectOriented C++ Library for Cosmology, arXiv:1312.4961 astroph.IM (2013).
 [26] J. Lesgourgues, The Cosmic Linear Anisotropy Solving System (CLASS) I: Overview, arXiv:1104.2932 (2011).
 [27] D. Blas, J. Lesgourgues, and T. Tram, The Cosmic Linear Anisotropy Solving System (CLASS). Part II: Approximation schemes, Journal of Cosmology and Astroparticle Physics 2011 (2011), no. 07 034.
 [28] F. Feroz and M. P. Hobson, Multimodal nested sampling: an efficient and robust alternative to Markov Chain Monte Carlo methods for astronomical data analyses, Monthly Notices of the Royal Astronomical Society 384 (2008), no. 2 449–463.
 [29] F. Feroz, M. P. Hobson, and M. Bridges, MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics, Monthly Notices of the Royal Astronomical Society 398 (2009), no. 4 1601–1614.
 [30] F. Feroz, M. P. Hobson, E. Cameron, and A. N. Pettitt, Importance Nested Sampling and the MultiNest Algorithm, arXiv:1306.2144 (2013).
 [31] M. Bridges, A. Lasenby, and M. Hobson, A bayesian analysis of the primordial power spectrum, Mon.Not.Roy.Astron.Soc. 369 (2006) 1123–1130, [astroph/0511573].
 [32] M. Bridges, A. Lasenby, and M. Hobson, WMAP 3year primordial power spectrum, Mon.Not.Roy.Astron.Soc. (2006) [astroph/0607404].
 [33] M. Bridges, F. Feroz, M. Hobson, and A. Lasenby, Bayesian optimal reconstruction of the primordial power spectrum, arXiv:0812.3541.
 [34] J. A. Vazquez, M. Bridges, M. P. Hobson, and A. N. Lasenby, Model selection applied to reconstruction of the Primordial Power Spectrum, Journal of Cosmology and Astroparticle Physics 06 (2012), no. 0 006.
 [35] J. A. Vazquez, A. Lasenby, M. Bridges, and M. Hobson, A Bayesian study of the primordial power spectrum from a novel closed universe model, Mon.Not.Roy.Astron.Soc. 422 (2012) 1948–1956, [arXiv:1103.4619].
 [36] J. A. Vazquez, M. Bridges, Y.Z. Ma, and M. Hobson, Constraints on the tensortoscalar ratio for nonpowerlaw models, JCAP 1308 (2013) 001, [arXiv:1303.4014].
 [37] R. Easther and H. Peiris, Bayesian Analysis of Inflation II: Model Selection and Constraints on Reheating, arXiv (Dec., 2011) [arXiv:1112.0326].
 [38] H. Jeffreys, Theory of Probability. Oxford University Press, 1961.
 [39] M. P. Hobson, A. H. Jaffe, A. R. Liddle, P. Mukherjee, and D. Parkinson, Bayesian methods in cosmology. Cambridge University Press, 2010.
 [40] B. Audren, D. G. Figueroa, and T. Tram, A note of clarification: BICEP2 and Planck are not in tension, arXiv:1405.1390.
 [41] L. Verde and H. V. Peiris, On minimally parametric primordial power spectrum reconstruction and the evidence for a red tilt, Journal of Cosmology and Astroparticle Physics 07 (2008), no. 0 009.
 [42] H. V. Peiris and L. Verde, The shape of the primordial power spectrum: A last stand before Planck data, Physical Review D 81 (2010), no. 2 21302.
 [43] S. Bird, H. V. Peiris, M. Viel, and L. Verde, Minimally Parametric Power Spectrum Reconstruction from the Lymanalpha Forest, Mon.Not.Roy.Astron.Soc. 413 (2011) 1717–1728, [arXiv:1010.1519].
 [44] R. de Putter, E. V. Linder, and A. Mishra, Inflationary Freedom and Cosmological Neutrino Constraints, arXiv:1401.7022 astroph.CO (2014).
 [45] K. Abazajian, K. Kadota, and E. D. Stewart, Parameterizing the power spectrum: Beyond the truncated Taylor expansion, JCAP 0508 (2005) 008, [astroph/0507224].
 [46] D. K. Hazra, A. Shafieloo, and T. Souradeep, Cosmological parameter estimation with freeform primordial power spectrum, Physical Review D 87 (2013), no. 1 123528.
 [47] A. Shafieloo and T. Souradeep, Assumptions of the primordial spectrum and cosmological parameter estimation, New Journal of Physics 13 (2011), no. 1 3024.
 [48] S. L. Bridle, A. M. Lewis, J. Weller, and G. Efstathiou, Reconstructing the primordial power spectrum, Monthly Notice of the Royal Astronomical Society 342 (2003), no. 4 L72–L78.
 [49] M. Tegmark and M. Zaldarriaga, Separating the early universe from the late universe: Cosmological parameter estimation beyond the black box, Physical Review D 66 (2002), no. 1 103508.
 [50] W. H. Kinney, How to fool cosmic microwave background parameter estimation, Physical Review D 63 (2001), no. 4 43001.
 [51] E. J. Copeland, E. W. Kolb, A. R. Liddle, and J. E. Lidsey, Observing the inflaton potential, Phys.Rev.Lett. 71 (1993) 219–222, [hepph/9304228].
 [52] E. J. Copeland, E. W. Kolb, A. R. Liddle, and J. E. Lidsey, Reconstructing the inflation potential, in principle and in practice, Phys.Rev. D48 (1993) 2529–2547, [hepph/9303288].
 [53] J. E. Lidsey, A. R. Liddle, E. W. Kolb, E. J. Copeland, T. Barreiro, et. al., Reconstructing the inflation potential : An overview, Rev.Mod.Phys. 69 (1997) 373–410, [astroph/9508078].
 [54] S. H. Hansen and M. Kunz, Observational constraints on the inflaton potential, Mon.Not.Roy.Astron.Soc. 336 (2002) 1007–1010, [hepph/0109252].
 [55] S. M. Leach, A. R. Liddle, J. Martin, and D. J. Schwarz, Cosmological parameter estimation and the inflationary cosmology, Phys.Rev. D66 (2002) 023515, [astroph/0202094].
 [56] S. M. Leach and A. R. Liddle, Microwave background constraints on inflationary parameters, Mon.Not.Roy.Astron.Soc. 341 (2003) 1151, [astroph/0207213].
 [57] S. M. Leach and A. R. Liddle, Constraining slow  roll inflation with WMAP and 2dF, Phys.Rev. D68 (2003) 123508, [astroph/0306305].
 [58] W. H. Kinney, E. W. Kolb, A. Melchiorri, and A. Riotto, Inflation model constraints from the Wilkinson Microwave Anisotropy Probe threeyear data, Phys.Rev. D74 (2006) 023502, [astroph/0605338].
 [59] H. Peiris and R. Easther, Recovering the Inflationary Potential and Primordial Power Spectrum With a Slow Roll Prior: Methodology and Application to WMAP 3 Year Data, JCAP 0607 (2006) 002, [astroph/0603587].
 [60] H. Peiris and R. Easther, Slow Roll Reconstruction: Constraints on Inflation from the 3 Year WMAP Dataset, JCAP 0610 (2006) 017, [astroph/0609003].
 [61] J. Martin and C. Ringeval, Inflation after WMAP3: Confronting the SlowRoll and Exact Power Spectra to CMB Data, JCAP 0608 (2006) 009, [astroph/0605367].
 [62] J. Lesgourgues and W. Valkenburg, New constraints on the observable inflaton potential from WMAP and SDSS, Phys.Rev. D75 (2007) 123519, [astroph/0703625].
 [63] J. Lesgourgues, A. A. Starobinsky, and W. Valkenburg, What do WMAP and SDSS really tell about inflation?, JCAP 0801 (2008) 010, [arXiv:0710.1630].
 [64] H. V. Peiris and R. Easther, Primordial Black Holes, Eternal Inflation, and the Inflationary Parameter Space after WMAP5, JCAP 0807 (2008) 024, [arXiv:0805.2154].
 [65] P. Adshead and R. Easther, Constraining Inflation, JCAP 0810 (2008) 047, [arXiv:0802.3898].
 [66] J. Hamann, J. Lesgourgues, and W. Valkenburg, How to constrain inflationary parameter space with minimal priors, JCAP 0804 (2008) 016, [arXiv:0802.0505].
 [67] W. H. Kinney, E. W. Kolb, A. Melchiorri, and A. Riotto, Latest inflation model constraints from cosmic microwave background measurements, Phys.Rev. D78 (2008) 087302, [arXiv:0805.2966].
 [68] J. Martin, C. Ringeval, and R. Trotta, Hunting Down the Best Model of Inflation with Bayesian Evidence, Phys.Rev. D83 (2011) 063524, [arXiv:1009.4157].
 [69] R. Easther and H. V. Peiris, Bayesian Analysis of Inflation II: Model Selection and Constraints on Reheating, Phys.Rev. D85 (2012) 103533, [arXiv:1112.0326].
 [70] R. Easther and H. Peiris, Implications of a Running Spectral Index for Slow Roll Inflation, JCAP 0609 (2006) 010, [astroph/0604214].