Planck 2013 results. XXII. Constraints on inflation
^{1}
APC, AstroParticule et Cosmologie, Université Paris Diderot,
CNRS/IN2P3, CEA/lrfu, Observatoire de Paris, Sorbonne Paris Cité, 10 rue Alice Domon et Léonie
Duquet, 75205
Paris Cedex 13,
France
^{2}
Aalto University Metsähovi Radio Observatory,
Metsähovintie 114, 02540
Kylmälä,
Finland
^{3}
African Institute for Mathematical Sciences,
68 Melrose Road, Muizenberg,
7701
Cape Town, South
Africa
^{4}
Agenzia Spaziale Italiana Science Data Center, via del Politecnico
snc, 00133
Roma,
Italy
^{5}
Agenzia Spaziale Italiana, Viale Liegi 26,
00198
Roma,
Italy
^{6}
Astrophysics Group, Cavendish Laboratory, University of
Cambridge, J J Thomson
Avenue, Cambridge
CB3 0HE,
UK
^{7}
Astrophysics & Cosmology Research Unit, School of
Mathematics, Statistics & Computer Science, University of
KwaZuluNatal, Westville Campus,
Private Bag X54001, 4000
Durban, South
Africa
^{8}
Atacama Large Millimeter/submillimeter Array, ALMA Santiago
Central Offices, Alonso de Cordova 3107, Vitacura, Casilla
763 0355
Santiago,
Chile
^{9}
CITA, University of Toronto, 60 St. George St., Toronto, ON
M5S 3H8,
Canada
^{10}
CNRS, IRAP, 9
Av. colonel Roche, BP
44346, 31028
Toulouse Cedex 4,
France
^{11}
California Institute of Technology, Pasadena, California, USA
^{12}
Centre for Theoretical Cosmology, DAMTP, University of
Cambridge, Wilberforce
Road, Cambridge
CB3 0WA,
UK
^{13}
Centro de Estudios de Física del Cosmos de Aragón (CEFCA), Plaza
San Juan, 1, planta 2, 44001
Teruel,
Spain
^{14}
Computational Cosmology Center, Lawrence Berkeley National
Laboratory, Berkeley,
California,
USA
^{15}
Consejo Superior de Investigaciones Científicas
(CSIC), 28006
Madrid,
Spain
^{16}
DSM/Irfu/SPP, CEASaclay, 91191
GifsurYvette Cedex,
France
^{17}
DTU Space, National Space Institute, Technical University of
Denmark, Elektrovej
327, 2800
Kgs. Lyngby,
Denmark
^{18}
Département de Physique Théorique, Université de
Genève, 24 quai E.
Ansermet, 1211
Genève 4,
Switzerland
^{19}
Departamento de Física Fundamental, Facultad de Ciencias,
Universidad de Salamanca, 37008
Salamanca,
Spain
^{20}
Departamento de Física, Universidad de Oviedo,
Avda. Calvo Sotelo s/n,
33007
Oviedo,
Spain
^{21}
Department of Astronomy and Astrophysics, University of
Toronto, 50 Saint George Street,
Toronto, Ontario,
Canada
^{22}
Department of Astrophysics/IMAPP, Radboud University
Nijmegen, PO Box
9010, 6500 GL
Nijmegen, The
Netherlands
^{23}
Department of Electrical Engineering and Computer Sciences,
University of California, Berkeley, California, USA
^{24}
Department of Physics & Astronomy, University of British
Columbia, 6224 Agricultural Road,
Vancouver, British
Columbia, Canada
^{25}
Department of Physics and Astronomy, Dana and David Dornsife
College of Letter, Arts and Sciences, University of Southern California,
Los Angeles, CA
90089,
USA
^{26}
Department of Physics and Astronomy, University College
London, London
WC1E 6BT,
UK
^{27}
Department of Physics and Astronomy, University of
Sussex, Brighton
BN1 9QH,
UK
^{28}
Department of Physics, Florida State University,
Keen Physics Building, 77 Chieftan
Way, Tallahassee,
Florida,
USA
^{29}
Department of Physics, Gustaf Hällströmin katu 2a, University of
Helsinki, Helsinki,
Finland
^{30}
Department of Physics, Princeton University,
Princeton, New Jersey, USA
^{31}
Department of Physics, University of California,
Berkeley, California, USA
^{32}
Department of Physics, University of California,
One Shields Avenue,
Davis, California, USA
^{33}
Department of Physics, University of California,
Santa Barbara, California, USA
^{34}
Department of Physics, University of Illinois at
UrbanaChampaign, 1110 West Green
Street, Urbana,
Illinois,
USA
^{35}
Dipartimento di Fisica e Astronomia G. Galilei, Università degli
Studi di Padova, via Marzolo
8, 35131
Padova,
Italy
^{36}
Dipartimento di Fisica e Scienze della Terra, Università di
Ferrara, via Saragat
1, 44122
Ferrara,
Italy
^{37}
Dipartimento di Fisica, Università La Sapienza,
P.le A. Moro 2, 00185
Roma,
Italy
^{38}
Dipartimento di Fisica, Università degli Studi di
Milano, via Celoria,
16, 20133
Milano,
Italy
^{39}
Dipartimento di Fisica, Università degli Studi di
Trieste, via A. Valerio
2, 34127
Trieste,
Italy
^{40}
Dipartimento di Fisica, Università di Roma Tor
Vergata, via della Ricerca
Scientifica, 1, 00133
Roma,
Italy
^{41}
Discovery Center, Niels Bohr Institute,
Blegdamsvej 17, 2100
Copenhagen,
Denmark
^{42}
Dpto. Astrofísica, Universidad de La Laguna (ULL),
38206
La Laguna, Tenerife,
Spain
^{43}
European Southern Observatory, ESO Vitacura, Alonso de Cordova 3107, Vitacura, Casilla
19001
Santiago,
Chile
^{44}
European Space Agency, ESAC, Planck Science Office, Camino bajo
del Castillo, s/n, Urbanización Villafranca del Castillo, Villanueva de la
Cañada, 28692
Madrid,
Spain
^{45}
European Space Agency, ESTEC, Keplerlaan 1,
2201 AZ
Noordwijk, The
Netherlands
^{46}
Haverford College Astronomy Department,
370 Lancaster Avenue,
Haverford, Pennsylvania, USA
^{47}
Helsinki Institute of Physics, Gustaf Hällströmin katu 2,
University of Helsinki, 00014
Helsinki,
Finland
^{48}
INAF – Osservatorio Astronomico di Padova, Vicolo
dell’Osservatorio 5, 35122
Padova,
Italy
^{49}
INAF – Osservatorio Astronomico di Roma, via di Frascati
33, 00040
Monte Porzio Catone,
Italy
^{50}
INAF – Osservatorio Astronomico di Trieste, via G.B. Tiepolo
11, 34131
Trieste,
Italy
^{51}
INAF Istituto di Radioastronomia, via P. Gobetti
101, 40129
Bologna,
Italy
^{52}
INAF/IASF Bologna, via Gobetti 101, 40129
Bologna,
Italy
^{53}
INAF/IASF Milano, via E. Bassini 15, 20133
Milano,
Italy
^{54}
INFN, Sezione di Bologna, via Irnerio 46,
40126
Bologna,
Italy
^{55}
INFN, Sezione di Roma 1, Università di Roma
Sapienza, P.le Aldo Moro
2, 00185
Roma,
Italy
^{56}
IPAG: Institut de Planétologie et d’Astrophysique de Grenoble,
Université Joseph Fourier, Grenoble 1/CNRSINSU, UMR 5274,
38041
Grenoble,
France
^{57}
IUCAA, Post Bag 4, Ganeshkhind, Pune University
Campus, 411 007
Pune,
India
^{58}
Imperial College London, Astrophysics group, Blackett
Laboratory, Prince Consort
Road, London,
SW7 2AZ,
UK
^{59}
Infrared Processing and Analysis Center, California Institute of
Technology, Pasadena,
CA
91125,
USA
^{60}
Institut Néel, CNRS, Université Joseph Fourier Grenoble
I, 25 rue des
Martyrs, 38042
Grenoble,
France
^{61}
Institut Universitaire de France, 103 bd SaintMichel, 75005
Paris,
France
^{62}
Institut d’Astrophysique Spatiale, CNRS (UMR 8617) Université
ParisSud 11, Bâtiment
121, 91405
Orsay,
France
^{63}
Institut d’Astrophysique de Paris, CNRS (UMR 7095),
98bis boulevard Arago,
75014
Paris,
France
^{64}
Institute for Space Sciences, BucharestMagurale,
Romania
^{65}
Institute of Astronomy and Astrophysics, Academia
Sinica, 10617
Taipei,
Taiwan
^{66}
Institute of Astronomy, University of Cambridge,
Madingley Road, Cambridge
CB3 0HA,
UK
^{67}
Institute of Theoretical Astrophysics, University of
Oslo, Blindern,
0315
Oslo,
Norway
^{68}
Instituto de Astrofísica de Canarias, C/Vía Láctea s/n, 38200 La Laguna,
Tenerife,
Spain
^{69}
Instituto de Física de Cantabria (CSICUniversidad de
Cantabria), Avda. de los Castros
s/n, Santander,
Spain
^{70}
Jet Propulsion Laboratory, California Institute of
Technology, 4800 Oak Grove
Drive, Pasadena,
California,
USA
^{71}
Jodrell Bank Centre for Astrophysics, Alan Turing Building, School
of Physics and Astronomy, The University of Manchester, Oxford Road, Manchester, M13
9PL, UK
^{72}
Kavli Institute for Cosmology Cambridge,
Madingley Road, Cambridge, CB3 0HA, UK
^{73}
LAL, Université ParisSud, CNRS/IN2P3,
Orsay,
France
^{74}
LERMA, CNRS, Observatoire de Paris, 61 avenue de
l’Observatoire, 75014
Paris,
France
^{75}
Laboratoire AIM, IRFU/Service d’Astrophysique – CEA/DSM – CNRS –
Université Paris Diderot, Bât. 709, CEASaclay, 91191
GifsurYvette Cedex,
France
^{76}
Laboratoire Traitement et Communication de l’Information, CNRS
(UMR 5141) and Télécom ParisTech, 46 rue Barrault, 75634
Paris Cedex 13,
France
^{77}
Laboratoire de Physique Subatomique et de Cosmologie, Université
Joseph Fourier Grenoble I, CNRS/IN2P3, Institut National Polytechnique de
Grenoble, 53 rue des
Martyrs, 38026
Grenoble Cedex,
France
^{78}
Laboratoire de Physique Théorique, Université ParisSud 11
& CNRS, Bâtiment
210, 91405
Orsay,
France
^{79}
Lawrence Berkeley National Laboratory,
Berkeley, California, USA
^{80}
Leung Center for Cosmology and Particle Astrophysics, National
Taiwan University, 10617
Taipei,
Taiwan
^{81}
MaxPlanckInstitut für Astrophysik, KarlSchwarzschildStr. 1, 85741
Garching,
Germany
^{82}
McGill Physics, Ernest Rutherford Physics Building, McGill
University, 3600 rue University, Montréal, QC,
H3A 2T8,
Canada
^{83}
MilliLab, VTT Technical Research Centre of Finland, Tietotie
3, 02044
Espoo,
Finland
^{84}
National University of Ireland, Department of Experimental
Physics, Maynooth,Co.
Kildare, Ireland
^{85}
Niels Bohr Institute, Blegdamsvej 17, 2100
Copenhagen,
Denmark
^{86}
Observational Cosmology, Mail Stop 36717, California Institute of
Technology, Pasadena
CA
91125,
USA
^{87}
Optical Science Laboratory, University College
London, Gower
Street, London,
UK
^{88}
SBITPLPPC, EPFL, 1015
Lausanne,
Switzerland
^{89}
SISSA, Astrophysics Sector, via Bonomea 265,
34136
Trieste,
Italy
^{90}
School of Physics and Astronomy, Cardiff University,
Queens Buildings, The Parade,
Cardiff, CF24 3AA, UK
^{91}
School of Physics and Astronomy, University of
Nottingham, Nottingham
NG7 2RD,
UK
^{92}
Space Research Institute (IKI), Russian Academy of
Sciences, Profsoyuznaya Str,
84/32, 117997
Moscow,
Russia
^{93}
Space Sciences Laboratory, University of California,
Berkeley, California, USA
^{94}
Special Astrophysical Observatory, Russian Academy of
Sciences, Nizhnij Arkhyz,
Zelenchukskiy region, 369167
KarachaiCherkessian Republic,
Russia
^{95}
Stanford University, Dept of Physics, Varian Physics Bldg, 382 via Pueblo
Mall, Stanford,
California,
USA
^{96}
SubDepartment of Astrophysics, University of
Oxford, Keble Road,
Oxford
OX1 3RH,
UK
^{97}
Theory Division, PHTH, CERN, 1211
Geneva 23,
Switzerland
^{98}
UPMC Univ Paris 06, UMR7095, 98bis boulevard Arago, 75014
Paris,
France
^{99}
Université de Toulouse, UPSOMP, IRAP, 31028
Toulouse Cedex 4,
France
^{100}
University of Granada, Departamento de Física Teórica y del
Cosmos, Facultad de Ciencias, 18071
Granada,
Spain
^{101}
Warsaw University Observatory, Aleje Ujazdowskie 4, 00478
Warszawa,
Poland
Received:
25
March
2013
Accepted:
28
January
2014
We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined with the WMAP largeangle polarization, constrain the scalar spectral index to be n_{s} = 0.9603 ± 0.0073, ruling out exact scale invariance at over 5σ.Planck establishes an upper bound on the tensortoscalar ratio of r< 0.11 (95% CL). The Planck data thus shrink the space of allowed standard inflationary models, preferring potentials with V′′< 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential models of degree n ≥ 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, obtaining dn_{s}/ dlnk = − 0.0134 ± 0.0090. We verify these conclusions through a numerical analysis, which makes no slowroll approximation, and carry out a Bayesian parameter estimation and modelselection analysis for a number of inflationary models including monomial, natural, and hilltop potentials. For each model, we present the Planck constraints on the parameters of the potential and explore several possibilities for the postinflationary entropy generation epoch, thus obtaining nontrivial datadriven constraints. We also present a direct reconstruction of the observable range of the inflaton potential. Unless a quartic term is allowed in the potential, we find results consistent with secondorder slowroll predictions. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit by Δχ^{2}_{eff} ≈ 10; however, Bayesian evidence does not prefer these models. We constrain several singlefield inflation models with generalized Lagrangians by combining power spectrum data with Planck bounds on f_{NL}. Planck constrains with unprecedented accuracy the amplitude and possible correlation (with the adiabatic mode) of nondecaying isocurvature fluctuations. The fractional primordial contributions of cold dark matter (CDM) isocurvature modes of the types expected in the curvaton and axion scenarios have upper bounds of 0.25% and 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelated isocurvature component can improve the χ^{2}_{eff} by approximately 4 as a result of slightly lowering the theoretical prediction for the ℓ ≲ 40 multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions.
Key words: cosmic background radiation / inflation / early Universe
© ESO, 2014
1. Introduction
This paper, one of a set associated with the 2013 release of data from the Planck^{1} mission (Planck Collaboration I–XXXI 2014), describes the implications of the Planck measurement of cosmic microwave background (CMB) anisotropies for cosmic inflation. In this first release only the Planck temperature data resulting from the nominal mission are used, which includes 2.6 full surveys of the sky. The interpretation of the CMB polarization as seen by Planck will be presented in a later series of publications. This paper exploits the data presented in Planck Collaboration II (2014), Planck Collaboration XII (2014), Planck Collaboration XV (2014), and Planck Collaboration XVII (2014). Other closely related papers discuss the estimates of cosmological parameters in Planck Collaboration XVI (2014) and investigations of nonGaussianity in Planck Collaboration XXIV (2014).
In the early 1980s inflationary cosmology, which postulates an epoch of nearly exponential expansion, was proposed in order to resolve a number of puzzles of standard big bang cosmology such as the entropy, flatness, horizon, smoothness, and monopole problems (Brout et al. 1978; Starobinsky 1980; Kazanas 1980; Sato 1981; Guth 1981; Linde 1982; Albrecht & Steinhardt 1982; Linde 1983). During inflation, cosmological fluctuations resulting from quantum fluctuations are generated and can be calculated using the semiclassical theory of quantum fields in curved spacetime (Mukhanov & Chibisov 1981, 1982; Hawking 1982; Guth & Pi 1982; Starobinsky 1982; Bardeen et al. 1983; Mukhanov 1985).
Cosmological observations prior to Planck are consistent with the simplest models of inflation within the slowroll paradigm. Recent observations of the CMB anisotropies (Story et al. 2013; Bennett et al. 2013; Hinshaw et al. 2013; Hou et al. 2014; Das et al. 2014) and of largescale structure (Beutler et al. 2011; Padmanabhan et al. 2012; Anderson et al. 2012) indicate that our Universe is very close to spatially flat and has primordial density fluctuations that are nearly Gaussian and adiabatic and are described by a nearly scaleinvariant power spectrum. PrePlanck CMB observations also established that the amplitude of primordial gravitational waves, with a nearly scaleinvariant spectrum (Starobinsky 1979; Rubakov et al. 1982; Fabbri & Pollock 1983), is at most small.
Most of the results in this paper are based on the twopoint statistics of the CMB as measured by Planck, exploiting the data presented in Planck Collaboration XV (2014), Planck Collaboration XVI (2014), and Planck Collaboration XVII (2014). The Planck results testing the Gaussianity of the primordial CMB component are described in the companion papers Planck Collaboration XXIII (2014), Planck Collaboration XXIV (2014), and Planck Collaboration XXV (2014). Planck finds values for the nonGaussian f_{NL} parameter of the CMB bispectrum consistent with the Gaussian hypothesis (Planck Collaboration XXIV 2014). This result has important implications for inflation. The simplest slowroll inflationary models predict a level of f_{NL} of the same order as the slowroll parameters and therefore too small to be detected by Planck.
The paper is organized as follows. Section 2 reviews inflationary theory, emphasizing in particular those aspects used later in the paper. In Sect. 3 the statistical methodology and the Planck likelihood as well as the likelihoods from the other astrophysical data sets used here are described. Section 4 presents constraints on slowroll inflation and studies their robustness under generalizations of the minimal assumptions of our baseline cosmological model. In Sect. 5 Bayesian model comparison of several inflationary models is carried out taking into account the uncertainty from the end of inflation to the beginning of the radiation dominated era. Section 6 reconstructs the inflationary potential over the range corresponding to the scales observable in the CMB. In Sect. 7 a penalized likelihood reconstruction of the primordial perturbation spectrum is performed. Section 8 reports on a parametric search for oscillations and features in the primordial scalar power spectrum. Section 9 examines constraints on noncanonical singlefield models of inflation including the f_{NL} measurements from Planck Collaboration XXIV (2014). In Sect. 10 constraints on isocurvature modes are established, thus testing the hypothesis that initial conditions were solely adiabatic. We summarize our conclusions in Sect. 11. Appendix A is dedicated to the constraints on slowroll inflation derived by sampling the Hubble flow functions (HFF) in the analytic expressions for the scalar and tensor power spectra. Definitions of the most relevant symbols used in this paper can be found in Tables 1 and 2.
2. Lightning review of inflation
Cosmological parameter definitions.
Conventions and definitions for inflation physics.
Before describing cosmic inflation, which was developed in the early 1980s, it is useful to review the state of theory prior to its introduction. Lifshitz (1946; see also Lifshitz & Khalatnikov 1963) first wrote down and solved the equations for the evolution of linearized perturbations about a homogeneous and isotropic FriedmannLemaîtreRobertsonWalker spacetime within the framework of general relativity. The general framework adopted was based on two assumptions:

(i)
The cosmological perturbations can be described by asinglecomponent fluid, at very early times.

(ii)
The initial cosmological perturbations were statistically homogeneous and isotropic, and Gaussian.
These are the simplest–but by no means unique–assumptions for defining a stochastic process for the initial conditions. Assumption (i), where only a single adiabatic mode is excited, is just the simplest possibility. In Sect. 10 we shall describe isocurvature perturbations, where other available modes are excited, and report on the constraints established by Planck. Assumption (ii) is a priori more questionable given the understanding at the time. An appeal can be made to the fact that any physics at weak coupling could explain (ii), but at the time these assumptions were somewhat ad hoc.
Even with the strong assumptions (i) and (ii), comparisons with observations cannot be made without further restrictions on the functional form of the primordial power spectrum of largescale spatial curvature inhomogeneities ℛ, , where n_{s} is the (scalar) spectral index. The notion of a scaleinvariant (i.e., n_{s} = 1) primordial power spectrum was introduced by Harrison (1970), Zeldovich (1972), and Peebles & Yu (1970) to address this problem. These authors showed that a scaleinvariant power law was consistent with the crude constraints on large and smallscale perturbations available at the time. However, other than its mathematical simplicity, no compelling theoretical explanation for this Ansatz was put forth. An important current question, addressed in Sect. 4, is whether n_{s} = 1 (i.e., exact scale invariance) is consistent with the data, or whether there is convincing evidence for small deviations from exact scale invariance. Although the inflationary potential can be tuned to obtain n_{s} = 1, inflationary models generically predict deviations from n_{s} = 1, usually on the red side (i.e., n_{s}< 1).
2.1. Cosmic inflation
Inflation was developed in a series of papers by Brout et al. (1978), Starobinsky (1980), Kazanas (1980), Sato (1981), Guth (1981), Linde (1982, 1983), and Albrecht & Steinhardt (1982). By generating an equation of state with a significant negative pressure (i.e., w = p/ρ ≈ − 1) before the radiation epoch, inflation solves a number of cosmological conundrums (the monopole, horizon, smoothness, and entropy problems), which had plagued all cosmological models extrapolating a matterradiation equation of state all the way back to the singularity. Such an equation of state (p ≈ − ρ) and the resulting nearly exponential expansion are obtained from a scalar field, the inflaton, with a canonical kinetic term (i.e., ), slowly rolling in the framework of Einstein gravity.
The homogeneous evolution of the inflaton field φ is governed by the equation of motion (1)and the Friedmann equation (2)Here H = ȧ/a is the Hubble parameter, the subscript φ denotes the derivative with respect to φ,M_{pl} = (8πG)^{− 1 / 2} is the reduced Planck mass, and V is the potential. (We use units where c = ħ = 1.) The evolution during the stage of quasiexponential expansion, when the scalar field rolls slowly down the potential, can be approximated by neglecting the second time derivative in Eq. (1) and the kinetic energy term in Eq. (2), so that
Necessary conditions for the slowroll described above are ϵ_{V} ≪ 1 and  η_{V}  ≪ 1, where the slowroll parameters ϵ_{V} and η_{V} are defined as The analogous hierarchy of HFF slowroll parameters measures instead the deviation from an exact exponential expansion. This hierarchy is defined as ϵ_{1} = − Ḣ/H^{2}, , with i ≥ 1. By using Eqs. (3) and (4), we have that ϵ_{1} ≈ ϵ_{V}, ϵ_{2} ≈ − 2η_{V} + 4ϵ_{V}.
2.2. Quantum generation of fluctuations
Without quantum fluctuations, inflationary theory would fail. Classically, any initial spatial curvature or gradients in the scalar field, as well as any inhomogeneities in other fields, would rapidly decay away during the quasiexponential expansion. The resulting universe would be too homogeneous and isotropic compared with observations. Quantum fluctuations must exist in order to satisfy the uncertainty relations that follow from the canonical commutation relations of quantum field theory. The quantum fluctuations in the inflaton and in the transverse and traceless parts of the metric are amplified by the nearly exponential expansion yielding the scalar and tensor primordial power spectra, respectively.
Many essentially equivalent approaches to quantizing the linearized cosmological fluctuations can be found in the original literature (see, e.g., Mukhanov & Chibisov 1981; Hawking 1982; Guth & Pi 1982; Starobinsky 1982; Bardeen et al. 1983). A simple formalism, which we shall follow here, was introduced by Mukhanov (1988), Mukhanov et al. (1992), and Sasaki (1986). In this approach a gaugeinvariant inflaton fluctuation Q is constructed and canonically quantized. This gaugeinvariant variable Q is the inflaton fluctuation δφ(t,x) in the uniform curvature gauge. The mode function of the inflaton fluctuations δφ(t,x) obeys the evolution equation (7)where . The gaugeinvariant field fluctuation is directly related to the comoving curvature perturbation^{2}(8)Analogously, gravitational waves are described by the two polarization states (+ , ×) of the transverse traceless parts of the metric fluctuations and are amplified by the expansion of the universe as well (Grishchuk 1975). The evolution equation for their mode function is (9)Early discussions of the generation of gravitational waves during inflation include Starobinsky (1979), Rubakov et al. (1982), Fabbri & Pollock (1983), Abbott & Wise (1984), and Starobinsky (1985a).
Because the primordial perturbations are small, of order 10^{5}, the linearized Eqs. (7) and (9) provide an accurate description for the generation and subsequent evolution of the cosmological perturbations during inflation. In this paper we use two approaches for solving for the cosmological perturbations. Firstly, we use an approximate treatment based on the slowroll approximation described below. Secondly, we use an almost exact approach based on numerical integration of the ordinary differential Eqs. (7) and (9) for each value of the comoving wavenumber k. For fixed k the evolution may be divided into three epochs: (i) subHubble evolution; (ii) Hubble crossing evolution; and (iii) superHubble evolution. During (i) the wavelength is much smaller than the Hubble length, and the mode oscillates as it would in a nonexpanding universe (i.e., Minkowski space). Therefore we can proceed with quantization as we would in Minkowski space. We quantize by singling out the positive frequency solution, as in the BunchDavies vacuum (Bunch & Davies 1978). This epoch is the oscillating regime in the WKB approximation. In epoch (iii), by contrast, there are two solutions, a growing and a decaying mode, and the evolution becomes independent of k. We care only about the growing mode. On scales much larger than the Hubble radius (i.e., k ≪ aH), both curvature and tensor fluctuations admit solutions constant in time^{3}. All the interesting, or nontrivial, evolution takes place between epochs (i) and (iii) – that is, during (ii), a few efolds before and after Hubble crossing, and this is the interval where the numerical integration is most useful since the asymptotic expansions are not valid in this transition region. Two numerical codes are used in this paper, ModeCode (Adams et al. 2001; Peiris et al. 2003; Mortonson et al. 2009; Easther & Peiris 2012), and the inflation module of Lesgourgues & Valkenburg (2007) as implemented in CLASS (Lesgourgues 2011; Blas et al. 2011)^{4}.
It is convenient to expand the power spectra of curvature and tensor perturbations on superHubble scales as
where A_{s} (A_{t}) is the scalar (tensor) amplitude and n_{s} (n_{t}), dn_{s}/ dlnk (dn_{t}/ dlnk) and d^{2}n_{s}/ dlnk^{2} are the scalar (tensor) spectral index, the running of the scalar (tensor) spectral index, and the running of the running of the scalar spectral index, respectively.
The parameters of the scalar and tensor power spectra may be calculated approximately in the framework of the slowroll approximation by evaluating the following equations at the value of the inflation field φ_{∗} where the mode k_{∗} = a_{∗}H_{∗} crosses the Hubble radius for the first time. (For a nice review of the slowroll approximation, see for example Liddle & Lyth 1993.) The number of efolds before the end of inflation, N_{∗}, at which the pivot scale k_{∗} exits from the Hubble radius, is (12)where the equality holds in the slowroll approximation, and subscript e denotes the end of inflation.
The coefficients of Eqs. (10) and (11) at their respective leading orders in the slowroll parameters are given by
where the slowroll parameters ϵ_{V} and η_{V} are defined in Eqs. (5) and (6), and the higher order parameters are defined as (20)and (21)In singlefield inflation with a standard kinetic term, as discussed here, the tensor spectrum shape is not independent from the other parameters. The slowroll paradigm implies a tensortoscalar ratio at the pivot scale of (22)referred to as the consistency relation. This consistency relation is also useful to help understand how r is connected to the evolution of the inflaton: (23)The above relation, called the Lyth bound (Lyth 1997), implies that an inflaton variation of the order of the Planck mass is needed to produce r ≳ 0.01. Such a threshold is useful to classify large and smallfield inflationary models with respect to the Lyth bound.
2.3. Ending inflation and the epoch of entropy generation
The greatest uncertainty in calculating the perturbation spectrum predicted from a particular inflationary potential arises in establishing the correspondence between the comoving wavenumber today and the inflaton energy density when the mode of that wavenumber crossed the Hubble radius during inflation (Kinney & Riotto 2006). This correspondence depends both on the inflationary model and on the cosmological evolution from the end of inflation to the present.
After the slowroll stage, becomes as important as the cosmological damping term . Inflation ends gradually as the inflaton picks up kinetic energy so that w is no longer slightly above −1, but rather far from that value. We may arbitrarily deem that inflation ends when w = − 1 / 3 (the value dividing the cases of an expanding and a contracting comoving Hubble radius), or, equivalently, at ϵ_{V} ≈ 1, after which the epoch of entropy generation starts. Because of couplings to other fields, the energy initially in the form of scalar field vacuum energy is transferred to the other fields by perturbative decay (reheating), possibly preceded by a nonperturbative stage (preheating). There is considerable uncertainty about the mechanisms of entropy generation, or thermalization, which subsequently lead to a standard w = 1 / 3 equation of state for radiation.
On the other hand, if we want to identify some k_{∗} today with the value of the inflaton field at the time this scale left the Hubble radius, Eq. (12) needs to be matched to an expression that quantifies how much k_{∗} has shrunk relative to the size of the Hubble radius between the end of inflation and the time when that mode reenters the Hubble radius. This quantity depends both on the inflationary potential and the details of the entropy generation process and is given by (24)where ρ_{end} is the energy density at the end of inflation, ρ_{th} is an energy scale by which the universe has thermalized, a_{0}H_{0} is the present Hubble radius, V_{∗} is the potential energy when k_{∗} left the Hubble radius during inflation, w_{int} characterizes the effective equation of state between the end of inflation and the energy scale specified by ρ_{th}, and g_{th} is the number of effective bosonic degrees of freedom at the energy scale ρ_{th}. In predicting the primordial power spectra at observable scales for a specific inflaton potential, this uncertainty in the reheating history of the universe becomes relevant and can be taken into account by allowing N_{∗} to vary over a range of values. Note that w_{int} is not intended to provide a detailed model for entropy generation, but rather to parameterize the uncertainty regarding the expansion rate of the universe during this intermediate era. Nevertheless, constraints on w_{int} provide observational limits on the uncertain physics during this period.
The first two terms of Eq. (24) are model independent, with the second term being roughly 5 for k_{∗} = 0.05 Mpc^{1}. If thermalization occurs rapidly, or if the reheating stage is close to radiationlike, the magnitude of the second to last term in Eq. (24) is less than roughly unity. The magnitude of the ln(g_{th}) / 12 term is negligible, giving a shift of only 0.58 for the extreme value g_{th} = 10^{3}. For most reasonable inflation models, the fourth term is and the third term is approximately −10, motivating the commonly assumed range 50 <N_{∗}< 60. Nonetheless, more extreme values at both ends are in principle possible (Liddle & Leach 2003). In the figures of Sect. 4 we will mark the range 50 <N_{∗}< 60 as a general guide.
2.4. Perturbations from cosmic inflation at higher order
To calculate the quantum fluctuations generated during cosmic inflation, a linearized quantum field theory in a timedependent background can be used. The leading order is the twopoint correlation function (25)but the inflaton selfinteractions and the nonlinearity of Einstein gravity give small higherorder corrections, of which the nexttoleading order is the three point function (26)which is in general nonzero.
For singlefield inflation with a standard kinetic term in a smooth potential (with initial fluctuations in the BunchDavies vacuum), the nonGaussian contribution to the curvature perturbation during inflation is (Acquaviva et al. 2003; Maldacena 2003), i.e., at an undetectable level smaller than other general relativistic contributions, such as the crosscorrelation between the integrated SachsWolfe effect and weak gravitational lensing of the CMB. For a general scalar field Lagrangian, the nonGaussian contribution can be large enough to be accessible to Planck with f_{NL} of order (Chen et al. 2007), where c_{s} is the sound speed of inflaton fluctuations (see Sect. 9). Other higher order kinetic and spatial derivative terms contribute to larger nonGaussianities. For a review of nonGaussianity generated during inflation, see, for example, Bartolo et al. (2004a) and Chen (2010) as well as the companion paper Planck Collaboration XXIV (2014).
2.5. Multifield models of cosmic inflation
Inflation as described so far assumes a single scalar field that drives and terminates the quasiexponential expansion and also generates the largescale curvature perturbations. When there is more than one field with an effective mass smaller than H, isocurvature perturbations are also generated during inflation by the same mechanism of amplification due to the stretching of the spacetime geometry (Axenides et al. 1983; Linde 1985). Cosmological perturbations in models with an Mcomponent inflaton φ_{i} can be analysed by considering perturbations parallel and perpendicular to the classical trajectory, as treated for example in Gordon et al. (2001). The definition of curvature perturbation generalizing Eq. (8) to the multifield case is (27)where Q_{i} is the gaugeinvariant field fluctuation associated with φ_{i} and . The above formula for the curvature perturbation can also be obtained through the δN formalism, i.e., , where the number of efolds to the end of inflation N is generalized to the multifield case (Starobinsky 1985b; Sasaki & Stewart 1996). The M − 1 normal directions are connected to M − 1 isocurvature perturbations δs_{ij} according to (28)If the trajectory of the average field is curved in field space, then during inflation both curvature and isocurvature fluctuations are generated with nonvanishing correlations (Langlois 1999).
Isocurvature perturbations can be converted into curvature perturbations on large scales, but the opposite does not hold (Mollerach 1990). If such isocurvature perturbations are not totally converted into curvature perturbations, they can have observable effects on CMB anisotropies and on structure formation. In Sect. 10, we present the Planck constraints on a combination of curvature and isocurvature initial conditions and the implications for important twofield scenarios, such as the curvaton (Lyth & Wands 2002) and axion (Lyth 1990) models.
Isocurvature perturbations may lead to a higher level of nonGaussianity compared to a single inflaton with a standard kinetic term (Groot Nibbelink & van Tent 2000). There is no reason to expect the inflaton to be a singlecomponent field. The scalar sector of the Standard Model, as well as its extensions, contains more than one scalar field.
3. Methodology
3.1. Cosmological model and parameters
The parameters of the models to be estimated in this paper fall into three categories: (i) parameters describing the initial perturbations, i.e., characterizing the particular inflationary scenario in question; (ii) parameters determining cosmological evolution at late times (z ≲ 10^{4}); and (iii) parameters that quantify our uncertainty about the instrument and foreground contributions to the angular power spectrum. These will be described in Sect. 3.2.1.
Unless specified otherwise, we assume that the late time cosmology is the standard flat sixparameter ΛCDM model whose energy content consists of photons, baryons, cold dark matter, neutrinos (assuming N_{eff} = 3.046 effective species, one of which is taken to be massive, with a mass of m_{ν} = 0.06 eV), and a cosmological constant. The primordial helium fraction, Y_{P}, is set as a function of Ω_{b}h^{2} and N_{eff} according to the big bang nucleosynthesis consistency condition (Ichikawa & Takahashi 2006; Hamann et al. 2008b), and we fix the CMB mean temperature to T_{0} = 2.7255 K (Fixsen 2009). Reionization is modelled to occur instantaneously at a redshift z_{re}, and the optical depth τ is calculated as a function of z_{re}. This model can be characterized by four free cosmological parameters: Ω_{b}h^{2},Ω_{c}h^{2},θ_{MC}, and τ, defined in Table 1, in addition to the parameters describing the initial perturbations.
3.2. Data
The primary CMB data used for this paper consist of the Planck CMB temperature likelihood supplemented by the Wilkinson Microwave Anisotropy Probe (WMAP) largescale polarization likelihood (henceforth Planck+WP), as described in Sect. 3.2.1. The largeangle Emode polarization spectrum is important for constraining reionization because it breaks the degeneracy in the temperature data between the primordial power spectrum amplitude and the optical depth to reionization. In the analysis constraining cosmic inflation, we restrict ourselves to combining the Planck temperature data with various combinations of the following additional data sets: the Planck lensing power spectrum, other CMB data extending the Planck data to higher ℓ, and BAO data. For the higherresolution CMB data we use measurements from the Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT). These complementary data sets are among the most useful to break degeneracies in parameters. The consequences of including other data sets such as Supernovae Type Ia (SN Ia) or the local measurement of the Hubble constant H_{0} on some of the cosmological models discussed here can be found in the compilation of cosmological parameters for numerous models included in the online Planck Legacy archive^{5}. Combining Planck+WP with various SN Ia data compilations (Conley et al. 2011; Suzuki et al. 2012) or with a direct measurement of H_{0} (Riess et al. 2011) does not significantly alter the conclusions for the simplest slowroll inflationary models presented below. The approach adopted here is the same as in the parameters paper Planck Collaboration XVI (2014).
3.2.1. Planck CMB temperature data
The Planck CMB likelihood is based on a hybrid approach, which combines a Gaussian likelihood approximation derived from temperature pseudo crossspectra at high multipoles (Hamimeche & Lewis 2008), with a pixelbased temperature and polarization likelihood at low multipoles. We summarize the likelihood here. For a detailed description the reader is referred to Planck Collaboration XV (2014).
The smallscale Planck temperature likelihood is based on pseudo crossspectra between pairs of maps at 100, 143, and 217 GHz, masked to retain 49%, 31%, and 31% of the sky, respectively. This results in angular auto and crosscorrelation power spectra covering multipole ranges of 50 ≤ ℓ ≤ 1200 at 100 GHz, 50 ≤ ℓ ≤ 2000 at 143 GHz, and 500 ≤ ℓ ≤ 2500 at 217 GHz as well as for the 143 × 217 GHz crossspectrum. In addition to instrumental uncertainties, mitigated here by using only crossspectra among different detectors, smallscale foreground and CMB secondary anisotropies need to be accounted for. The foreground model used in the Planck highℓ likelihood is described in detail in Planck Collaboration XV (2014) and Planck Collaboration XVI (2014), and includes contributions to the crossfrequency power spectra from unresolved radio point sources, the cosmic infrared background (CIB), and the thermal and kinetic SunyaevZeldovich effects. There are eleven adjustable nuisance parameters: (, , ). In addition, the calibration parameters for the 100 and 217 GHz channels, c_{100} and c_{217}, relative to the 143 GHz channel, and the dominant beam uncertainty eigenmode amplitude are left free in the analysis, with other beam uncertainties marginalized analytically. The Planck highℓ likelihood therefore includes 14 nuisance parameters^{6}.
The lowℓPlanck likelihood combines the Planck temperature data with the large scale 9year WMAP polarization data for this release. The procedure introduced in Page et al. (2007) separates the temperature and polarization likelihood under the assumption of negligible noise in the temperature map. The temperature likelihood uses Gibbs sampling (Eriksen et al. 2007), mapping out the distribution of the ℓ < 50 CMB temperature multipoles from a foregroundcleaned combination of the 30 − 353 GHz maps (Planck Collaboration XII 2014). The polarization likelihood is pixelbased using the WMAP 9year polarization maps at 33, 41, and 61 GHz and includes the temperaturepolarization crosscorrelation (Page et al. 2007). Its angular range is ℓ ≤ 23 for TE, EE, and BB.
3.2.2. Planck lensing data
The primary CMB anisotropies are distorted by the gravitational potential induced by intervening matter. Such lensing, which broadens and smooths out the acoustic oscillations, is taken into account as a correction to the observed temperature power spectrum. The lensing power spectrum can also be recovered by measuring higherorder correlation functions.
Some of our analysis includes the Planck lensing likelihood, derived in Planck Collaboration XVII (2014), which measures the nonGaussian trispectrum of the CMB and is proportional to the power spectrum of the lensing potential. As described in Planck Collaboration XVII (2014), this potential is reconstructed using quadratic estimators (Okamoto & Hu 2003), and its power spectrum is used to estimate the lensing deflection power spectrum. The spectrum is estimated from the 143 and 217 GHz maps, using multipoles in the range 40 <ℓ < 400. The theoretical predictions for the lensing potential power spectrum are calculated at linear order.
3.2.3. ACT and SPT temperature data
We include data from ACT and SPT to extend the multipole range of our CMB likelihood. ACT measures the power spectra and cross spectrum of the 148 and 218 GHz channels (Das et al. 2014), and covers angular scales 500 <ℓ < 10 000 at 148 GHz and 1500 <ℓ < 10 000 at 218 GHz. We use these data in the range ℓ > 1000 in combination with Planck. SPT measures the power spectrum for angular scales 2000 <ℓ < 10 000 at 95, 150, and 220 GHz (Reichardt et al. 2012). The spectrum at larger scales is also measured at 150 GHz (Story et al. 2013), but we do not include this data in our analysis. To model the foregrounds for ACT and SPT we follow a similar approach to the likelihood described in Dunkley et al. (2013), extending the model used for the Planck highℓ likelihood. Additional nuisance parameters are included to model the Poisson source amplitude, the residual Galactic dust contribution, and the interfrequency calibration parameters. More details are provided in Planck Collaboration XV (2014) and Planck Collaboration XVI (2014).
3.2.4. BAO data
The BAO (Baryon Acoustic Oscillation) angular scale serves as a standard ruler and allows us to map out the expansion history of the Universe after last scattering. The BAO scale, extracted from galaxy redshift surveys, provides a constraint on the latetime geometry and breaks degeneracies with other cosmological parameters. Galaxy surveys constrain the ratio , where is the spherically averaged distance scale to the effective survey redshift and r_{s} is the sound horizon (Mehta et al. 2012).
In this analysis we consider a combination of the measurements by the 6dFGRS (Beutler et al. 2011, ), SDSSII (Padmanabhan et al. 2012, ), and BOSS CMASS (Anderson et al. 2012, ) surveys, assuming no correlation between the three data points. This likelihood is described further in Planck Collaboration XVI (2014).
3.3. Parameter estimation
Given a model ℳ with free parameters x ≡ { x_{1},··· ,x_{k} } and a likelihood function of the data ℒ(data  x), the (posterior) probability density as a function of the parameters can be expressed as (29)where P(x  ℳ) represents the dataindependent prior probability density. Unless specified otherwise, we choose wide tophat prior distributions for all cosmological parameters.
We construct the posterior parameter probabilities using the Markov Chain Monte Carlo (MCMC) sampler as implemented in the CosmoMC(Lewis & Bridle 2002) or MontePython (Audren et al. 2012) packages. In some cases, when the calculation of the Bayesian evidence (see below) is desired or when the likelihood function deviates strongly from a multivariate Gaussian, we use the nested sampling algorithm provided by the MultiNest addon module (Feroz & Hobson 2008; Feroz et al. 2009) instead of the MetropolisHastings algorithm.
Joint twodimensional and onedimensional posterior distributions are obtained by marginalization. Numerical values and constraints on parameters are quoted in terms of the mean and 68% central Bayesian interval of the respective onedimensional marginalized posterior distribution.
3.4. Model selection
Two approaches to model selection are commonly used in statistics. The first approach examines the logarithm of the likelihood ratio, or effective χ^{2},(30)between models ℳ_{1} and ℳ_{2}, corrected for the fact that models with more parameters provide a better fit due to fitting away noise, even when the more complicated model is not correct. Various information criteria have been proposed based on this idea (Akaike 1974; Schwarz 1978); see also Liddle (2007). These quantities have the advantage of being independent of prior choice and fairly easy to calculate. The second approach is Bayesian (Cox 1946; Jeffreys 1998; Jaynes & Bretthorst 2003), and is based on evaluating ratios of the model averaged likelihood, or Bayesian evidence, defined by (31)Evidence ratios, also known as Bayes factors, B_{12} ≡ ℰ_{1}/ ℰ_{2}, are naturally interpreted as betting odds between models^{7}. Nested sampling algorithms allow rapid numerical evaluation of ℰ. In this paper we will consider both the effective χ^{2} and the Bayesian evidence^{8}.
4. Constraints on slowroll inflationary models
In this section we describe constraints on slowroll inflation using Planck+WP data in combination with the likelihoods described in Sects. 3.2.2–3.2.4. First we concentrate on characterizing the primordial power spectrum using Planck and other data. We start by showing that the empirical preinflationary HarrisonZeldovich (HZ) spectrum with n_{s} = 1 does not fit the Planck measurements. We further examine whether generalizing the cosmological model, for example by allowing the number of neutrino species to vary, allowing the helium fraction to vary, or admitting a nonstandard reionization scenario could reconcile the data with n_{s} = 1. We conclude that n_{s} ≠ 1 is robust.
We then investigate the Planck constraints on slowroll inflation, allowing a tilt for the spectral index and the presence of tensor modes, and discuss the implications for the simplest standard inflationary models. In this section the question is studied using the slowroll approximation, but later sections move beyond the slowroll approximation. We show that compared to previous experiments, Planck significantly narrows the space of allowed inflationary models. Next we consider evidence for a running of n_{s} and constrain it to be small, although we find a preference for negative running at modest statistical significance. Finally, we comment on the implications for inflation of the Planck constraints on possible deviations from spatial flatness.
4.1. Ruling out exact scale invariance
The simplest Ansatz for characterizing the statistical properties of the primordial cosmological perturbations is the socalled HZ model proposed by Harrison (1970), Zeldovich (1972), and Peebles & Yu (1970). These authors pointed out that a power spectrum with exact scale invariance for the Newtonian gravitation potential fitted the data available at the time, but without giving any theoretical justification for this form of the spectrum. Under exact scale invariance, which would constitute an unexplained new symmetry, the primordial perturbations in the Newtonian gravitational potential look statistically the same whether they are magnified or demagnified. In this simple model, vector and tensor perturbations are absent and the spectrum of curvature perturbations is characterized by a single parameter, the amplitude A_{s}. Inflation, on the other hand, generically breaks this rescaling symmetry. Although under inflation scale invariance still holds approximately, inflation must end. Therefore as different scales are imprinted, the physical conditions must evolve.
Although a detection of a violation of scale invariance would not definitively prove that inflation is responsible for the generation of the primordial perturbations, ruling out the HZ model would confirm the expectation of small deviations from scale invariance, almost always on the red side, which are generic to all inflationary models without fine tuning. We examine in detail the viability of the HZ model using statistics to compare to the more general model where the spectral index is allowed to vary, as motivated by slowroll inflation.
Constraints on cosmological parameters and best fit −2Δln(ℒ) with respect to the standard ΛCDM model, using Planck+WP data, testing the significance of the deviation from the HZ model.
When the cosmological model with n_{s} = 1 is compared with a model in which n_{s} is allowed to vary, we find that allowing n_{s} to deviate from one decreases the best fit effective χ^{2} by 27.9 with respect to the HZ model. Thus the significance of the finding that n_{s} ≠ 1 is in excess of 5σ. The parameters and maximum likelihood of this comparison are reported in Table 3.
One might wonder whether n_{s} = 1 could be reconciled with the data by relaxing some of the assumptions of the underlying cosmological model. Of particular interest is exploring those parameters almost degenerate with the spectral index such as the effective number of neutrino species N_{eff} and the primordial helium fraction Y_{P}, which both alter the damping tail of the temperature spectrum (Trotta & Hansen 2004; Hou et al. 2013), somewhat mimicking a spectral tilt. Assuming a HarrisonZeldovich spectrum and allowing N_{eff} or Y_{P} to float, and thus deviate from their standard values, gives almost as good a fit to Planck+WP data as the ΛCDM model with a varying spectral index, with and 2.2, respectively. However, as shown in Table 3, the HZ, HZ+N_{eff}, and HZ+Y_{P} models require significantly higher baryon densities and reionization optical depths compared to ΛCDM. In the HZ+Y_{P} model, one obtains a helium fraction of Y_{P} = 0.3194 ± 0.013. This value is incompatible both with direct measurements of the primordial helium abundance (Aver et al. 2012) and with standard big bang nucleosynthesis (Hamann et al. 2008b). (For comparison, we note that the value Y_{P} = 0.2477 was obtained as best fit for the ΛCDM model.) The HZ+N_{eff} model, on the other hand, would imply the presence of ΔN_{eff} ≈ 1 new effective neutrino species beyond the three known species. When BAO measurements are included in the likelihood, increases to 39.2 (HZ), 4.6 (HZ+Y_{P}), and 8.0 (HZ+N_{eff}), respectively, for the three models. The significance of this detection is also discussed in Planck Collaboration XVI (2014).
4.2. Constraining inflationary models using the slowroll approximation
We now consider all inflationary models that can be described by the primordial power spectrum parameters consisting of the scalar amplitude, A_{s}, the spectral index, n_{s}, and the tensortoscalar ratio r, all defined at the pivot scale k_{∗}. We assume that the spectral index is independent of the wavenumber k. Negligible running of the spectral index is expected if the slowroll condition is satisfied and higher order corrections in the slowroll approximations can be neglected. In the next subsection we relax this assumption.
Sampling the power spectrum parameters A_{s}, n_{s}, and r is not the only method for constraining slowroll inflation. Another possibility is to sample the Hubble flow functions in the analytic expressions for the scalar and tensor power spectra (Stewart & Lyth 1993; Gong & Stewart 2001; Leach et al. 2002). In the Appendix, we compare the slowroll inflationary predictions by sampling the HFF with Planck data and show that the results obtained in this way agree with those derived by sampling the power spectrum parameters. This confirms similar studies based on previous data (Hamann et al. 2008c; Finelli et al. 2010).
The spectral index estimated from Planck+WP data is (32)This tight bound on n_{s} is crucial for constraining inflation. The Planck constraint on r depends slightly on the pivot scales; we adopt k_{∗} = 0.002 Mpc^{1} to quote our results, with r_{0.002}< 0.12 at 95% CL. This bound improves on the most recent results, including the WMAP 9year constraint of r< 0.38 (Hinshaw et al. 2013), the WMAP 7year + ACT limit of r< 0.28 (Sievers et al. 2013), and the WMAP 7year + SPT limit of r< 0.18 (Story et al. 2013). The new bound from Planck is consistent with the theoretical limit from temperature anisotropies alone (Knox & Turner 1994). When a possible tensor component is included, the spectral index from Planck+WP does not significantly change, with n_{s} = 0.9624 ± 0.0075.
The Planck constraint on r corresponds to an upper bound on the energy scale of inflation (33)at 95% CL. This is equivalent to an upper bound on the Hubble parameter during inflation of H_{∗}/M_{pl}< 3.7 × 10^{5}. In terms of slowroll parameters, Planck+WP constraints imply ϵ_{V}< 0.008 at 95% CL, and .
Constraints on the primordial perturbation parameters in the ΛCDM+tensor model from Planck combined with other data sets.
Fig. 1 Marginalized joint 68% and 95% CL regions for n_{s} and r_{0.002} from Planck in combination with other data sets compared to the theoretical predictions of selected inflationary models. 

Open with DEXTER 
The Planck results on n_{s} and r are robust to the addition of external data sets (see Table 4). When the highℓ CMB ACT + SPT data are added, we obtain n_{s} = 0.9600 ± 0.0071 and r_{0.002}< 0.11 at 95% CL. Including the Planck lensing likelihood we obtain n_{s} = 0.9653 ± 0.0069 and r_{0.002}< 0.13, and adding BAO data gives n_{s} = 0.9643 ± 0.0059 and r_{0.002}< 0.12.
The above bounds are robust to small changes in the polarization likelihood at low multipoles. To test this robustness, instead of using the WMAP polarization likelihood, we impose a Gaussian prior τ = 0.07 ± 0.013 to take into account small shifts due to uncertainties in residual foreground contamination or instrument systematic effects in the evaluation of τ, as performed in Appendix B of Planck Collaboration XVI (2014). We find at most a reduction of 8% for the upper bound on r.
It is useful to plot the inflationary potentials in the n_{s}  r plane using the first two slowroll parameters evaluated at the pivot scale k_{∗} = 0.002 Mpc^{1} (Dodelson et al. 1997). Given our ignorance of the details of the epoch of entropy generation, we assume that the number of efolds N_{∗} to the end of inflation lies in the interval [50,60]. This uncertainty is plotted for those potentials predicting an exit from inflation without changing the potential.
Figure 1 shows the Planck constraints in the n_{s}  r plane and indicates the predictions of a number of representative inflationary potentials (see Lyth & Riotto 1999, for a review of particle physics models of inflation). The sensitivity of Planck data to high multipoles removes the degeneracy between n_{s} and r found using the WMAP data. Planck data favour models with a concave potential. As shown in Fig. 1, most of the joint 95% allowed region lies below the convex potential limit, and concave models with a red tilt in the range [0.945–0.98] are allowed by Planck at 95% CL. In the following we consider the status of several illustrative and commonly discussed inflationary potentials in light of the Planck observations.
Power law potential and chaotic inflation
The simplest class of inflationary models is characterized by a single monomial potential of the form (34)This class of potentials includes the simplest chaotic models, in which inflation starts from large values for the inflaton, φ>M_{pl}. Inflation ends when slowroll is no longer valid, and we assume this to occur at ϵ_{V} = 1. According to Eqs. (5), (6), and (15), this class of potentials predicts to lowest order in slowroll parameters , , . The λφ^{4} model lies well outside the joint 99.7% CL region in the n_{s}  r plane. This result confirms previous findings from, for example, Hinshaw et al. (2013), in which this model lies outside the 95% CL for the WMAP 9year data and is further excluded by CMB data at smaller scales.
The model with a quadratic potential, n = 2 (Linde 1983), often considered the simplest example for inflation, now lies outside the joint 95% CL for the Planck+WP+highℓ data for N_{∗} ≲ 60efolds, as shown in Fig. 1.
A linear potential with n = 1 (McAllister et al. 2010), motivated by axion monodromy, has η_{V} = 0 and lies within the 95% CL region. Inflation with n = 2 / 3 (Silverstein & Westphal 2008), however, also motivated by axion monodromy, now lies on the boundary of the joint 95% CL region. More permissive entropy generation priors allowing N_{∗}< 50 could reconcile this model with the Planck data.
Exponential potential and power law inflation
Inflation with an exponential potential (35)is called power law inflation (Lucchin & Matarrese 1985), because the exact solution for the scale factor is given by a(t) ∝ t^{2 /λ2}. This model is incomplete since inflation would not end without an additional mechanism to stop it. Under the assumption that such a mechanism exists and leaves predictions for cosmological perturbations unmodified, this class of models predicts r = − 8(n_{s} − 1) and now lies outside the joint 99.7% CL contour.
Inverse power law potential
Intermediate inflationary models (Barrow 1990; Muslimov 1990) with inverse power law potentials (36)lead to inflation with a(t) ∝ exp(At^{f}), with A> 0 and 0 <f< 1, where f = 4 / (4 + β) and β> 0. In intermediate inflation there is no natural end to inflation, but if the exit mechanism leaves the inflationary predictions for the cosmological perturbations unmodified, this class of models predicts and r ≈ − 8β(n_{s} − 1) / (β − 2) at lowest order in the slowroll approximation (Barrow & Liddle 1993)^{9}. Intermediate inflationary models lie outside the joint 95% CL contour for any β.
Hilltop models
In another interesting class of potentials, the inflaton rolls away from an unstable equilibrium as in the first new inflationary models (Albrecht & Steinhardt 1982; Linde 1982). We consider (37)where the ellipsis indicates higher order terms that are negligible during inflation but ensure positiveness of the potential later on. An exponent of p = 2 is allowed only as a large field inflationary model, predicting and . This potential leads to predictions in agreement with Planck+WP+BAO joint 95% CL contours for super Planckian values of μ, i.e., μ ≳ 9 M_{pl}.
Models with p ≥ 3 predict n_{s} − 1 ≈ − (2 /N)(p − 1) / (p − 2) when r ≪ 1. The hilltop potential with p = 3 lies outside the joint 95% CL region for Planck+WP+BAO data. The case with p = 4 is also in tension with Planck+WP+BAO, but allowed within the joint 95% CL region for N_{∗} ≳ 50 when r ≪ 1. For larger values of r these models provide a better fit to the Planck+WP+BAO data. The p = 4 hilltop model – without extra terms denoted by the ellipsis in Eq. (37) – is displayed in Fig. 1 in the standard range 50 <N_{∗}< 60 at different values of μ (this model approximates the linear potential for large μ/M_{pl}).
A simple symmetry breaking potential
The symmetry breaking potential (Olive 1990) (38)can be considered as a selfconsistent completion of the hilltop model with p = 2 (although it has a different limiting largefield branch for nonzero r). This potential leads to predictions in agreement with Planck + WP + BAO joint 95% CL contours for super Planckian values of μ (i.e. μ ≳ 13 M_{pl}).
Natural inflation
Another interesting class of potentials is natural inflation (Freese et al. 1990; Adams et al. 1993), initially motivated by its origin in symmetry breaking in an attempt to naturally give rise to the extremely flat potentials required for inflationary cosmology. In natural inflation the effective onedimensional potential takes the form (39)where f is a scale which determines the slope of the potential (see also Binétruy & Gaillard 1986, for an earlier motivation of a cosine potential for the inflaton in the context of superstring theory). Depending on the value of f, the model falls into the large field (f ≳ 1.5 M_{pl}) or small field (f ≲ 1.5 M_{pl}) categories. Therefore, holds for small f, while n_{s} ≈ 1−2 /N, r ≈ 8 /N holds for large f, approximating the m^{2}φ^{2} potential in the latter case (with ). This model agrees with Planck+WP data for f ≳ 5 M_{pl}.
Hybrid inflation
In hybrid inflationary models a second field, χ, coupled to the inflaton, undergoes symmetry breaking. The simplest example of this class is (40)Over most of their parameter space, these models behave effectively as singlefield models for the inflaton φ. The second field χ is close to the origin during the slowroll regime for φ, and inflation ends either by breakdown of slow roll for the inflaton at or by the waterfall transition of χ. The simplest models with (41)are disfavoured for most of the parameter space (Cortês & Liddle 2009). Models with m^{2}φ^{2}/ 2 ~ Λ^{4} are disfavoured due to a high tensortoscalar ratio, and models with U(φ) ≪ Λ^{4} predict a spectral index n_{s}> 1, also disfavoured by the Planck data.
We discuss hybrid inflationary models predicting n_{s}< 1 separately. As an example, the spontaneously broken SUSY model (Dvali et al. 1994) (42)predicts n_{s} − 1 ≈ − (1 + 3α_{h}/ 2) /N_{∗} and r ≈ 8α_{h}/N_{∗}. For α_{h} ≪ 1 and N_{∗} ≈ 50, n_{s} ≈ 0.98 is disfavoured by Planck+WP+BAO data at more than 95% CL. However, more permissive entropy generation priors allowing N_{∗}< 50 or a nonnegligible α_{h} give models consistent with the Planck data.
R^{2} inflation
Inflationary models can also be accommodated within extended theories of gravity. These theories can be analysed either in the original (Jordan) frame or in the conformallyrelated Einstein frame with a KleinGordon scalar field. Due to the invariance of curvature and tensor perturbation power spectra with respect to this conformal transformation, we can use the same methodology described earlier.
The first inflationary model proposed was of this type and was based on higher order gravitational terms in the action (Starobinsky 1980) (43)with the motivation to include semiclassical quantum effects. The predictions for R^{2} inflation were first studied in Mukhanov & Chibisov (1981) and Starobinsky (1983), and can be summarized as n_{s} − 1 ≈ − 8(4N_{∗} + 9) / (4N_{∗} + 3)^{2} and r ≈ 192 / (4N_{∗} + 3)^{2}. Since r is suppressed by another 1 /N_{∗} with respect to the scalar tilt, this model predicts a tiny amount of gravitational waves. This model predicts n_{s} = 0.963 for N_{∗} = 55 and is fully consistent with the Planck constraints.
Nonminimally coupled inflaton
A nonminimal coupling of the inflaton to gravity with the action (44)leads to several interesting consequences, such as a lowering of the tensortoscalar ratio.
The case of a massless selfinteracting inflaton (φ_{0} = 0) agrees with the Planck+WP data for ξ ≠ 0. Within the range 50 <N_{∗}< 60, this model is within the Planck+WP joint 95% CL region for ξ> 0.0019, improving on previous bounds (Tsujikawa & Gumjudpai 2004; Okada et al. 2010).
The amplitude of scalar perturbations is proportional to λ/ξ^{2} for ξ ≫ 1, and therefore the problem of tiny values for the inflaton selfcoupling λ can be alleviated (Spokoiny 1984; Lucchin et al. 1986; Salopek et al. 1989; Fakir & Unruh 1990). The regime φ_{0} ≪ M_{pl} is allowed and φ could be the Standard Model Higgs as proposed in Bezrukov & Shaposhnikov (2008) at tree level (see Barvinsky et al. 2008; Bezrukov & Shaposhnikov 2009, for the inclusion of loop corrections). The Higgs case with ξ ≫ 1 has the same predictions as the R^{2} model in terms of n_{s} and r as a function of N_{∗}. The entropy generation mechanism in the Higgs case can be more efficient than in the R^{2} case and therefore predicts a slightly larger n_{s} (Bezrukov & Gorbunov 2012). This model is fully consistent with the Planck constraints.
The case with ξ< 0 and was also recently emphasized in Linde et al. (2011). With the symmetry breaking potential in Eq. (44), the large field case with φ>φ_{0} is disfavoured by Planck data, whereas the small field case φ<φ_{0} is in agreement with the data.
Fig. 2 Marginalized joint 68% and 95% CL for (dn_{s}/ dlnk,n_{s}) using Planck+WP+BAO, either marginalizing over r or fixing r = 0 at k_{∗} = 0.038 Mpc^{1}. The purple strip shows the prediction for single monomial chaotic inflationary models with 50 <N_{∗}< 60 for comparison. 

Open with DEXTER 
4.3. Running spectral index
We have shown that the single parameter HarrisonZeldovich spectrum does not fit the data and that at least the first two terms A_{s} and n_{s} in the expansion of the primordial power spectrum in powers of ln(k) given in Eq. (10) are needed. Here we consider whether the data require the next term known as the running of the spectral index (Kosowsky & Turner 1995), defined as the derivative of the spectral index with respect to lnk, dn_{s ,t}/ dlnk for scalar or tensor fluctuations. If the slowroll approximation holds and the inflaton has reached its attractor solution, dn_{s}/ dlnk and dn_{t}/ dlnk are related to the potential slowroll parameters, as in Eqs. (17) and (18). In slowroll singlefield inflation, the running is second order in the Hubble slowroll parameters, for scalar and for tensor perturbations (Kosowsky & Turner 1995; Leach et al. 2002), and thus is typically suppressed with respect to n_{s} − 1 and n_{t}, which are first order. Given the tight constraints on the first two slowroll parameters ϵ_{V} and η_{V} (ϵ_{1} and ϵ_{2}) from the present data, typical values of the running to which Planck is sensitive (Pahud et al. 2007) would generically be dominated by the contribution from the third derivative of the potential, encoded in (or ϵ_{3}).
While it is easy to see that the running is invariant under a change in pivot scale, the same does not hold for the spectral index and the amplitude of the primordial power spectrum. It is convenient to choose k_{∗} such that dn_{s}/ dlnk and n_{s} are uncorrelated (Cortês et al. 2007). This approach minimizes the inferred variance of n_{s} and facilitates comparison with constraints on n_{s} in the power law models. Note, however, that the decorrelation pivot scale depends on both the model and the data set used.
Constraints on the primordial perturbation parameters for ΛCDM + dn_{s}/ dlnk, ΛCDM + dn_{s}/ dlnk + r, and ΛCDM + dn_{s}/ dlnk + d^{2}n_{s}/ dlnk^{2} models from Planck combined with other data sets.
We consider a model parameterizing the power spectrum using A_{s}(k_{∗}) ,n_{s}(k_{∗}), and dn_{s}/ dlnk, where k_{∗} = 0.05 Mpc^{1}. The joint constraints on n_{s} and dn_{s}/ dlnk at the decorrelation scale of Mpc^{1} are shown in Fig. 2. The Planck+WP constraints on the running do not change significantly when complementary data sets such as Planck lensing, CMB highℓ, and BAO data are included. We find (45)which is negative at the 1.5σ level. This reduces the uncertainty compared to previous CMB results. Error bars are reduced by 60% compared to the WMAP 9year results (Hinshaw et al. 2013), and by 20–30% compared to WMAP supplemented by SPT and ACT data (Hou et al. 2014; Sievers et al. 2013). Planck finds a smaller scalar running than SPT + WMAP7 (Hou et al. 2014), and larger than ACT + WMAP7 (Sievers et al. 2013). The best fit likelihood improves by only (3 when highℓ data are included) with respect to the minimal case in which n_{s} is scale independent, indicating that the deviation from scale independence is not very significant. The constraint for the spectral index in this case is 0.9630 ± 0.0065 at 68% CL at the decorrelation pivot scale k_{∗} = 0.038 Mpc^{1}. This result implies that the third derivative of the potential is small, i.e., , but compatible with zero at 95% CL, for inflation at low energy (i.e., with ϵ_{V} ≈ 0).
We also test the possibility that the running depends on the wavelength so that d^{2}n_{s}/ dlnk^{2} is nonzero. With Planck+WP data, we find . This result is stable with respect to the addition of complementary data sets, as can be seen from Table 5 and Fig. 3. When d^{2}n_{s}/ dlnk^{2} is allowed in the fit, we find a value for the running dn_{s}/ dlnk consistent with zero.
Fig. 3 Marginalized joint 68% and 95% CL regions for (d^{2}n_{s}/ dlnk^{2}, dn_{s}/ dlnk) using Planck+WP+BAO. 

Open with DEXTER 
Fig. 4 Marginalized joint 68% and 95% CL regions for (r,n_{s}), using Planck+WP+BAO with and without a running spectral index. 

Open with DEXTER 
Finally we allow a nonzero primordial gravitational wave spectrum together with the running. The tensor spectral index and its running are set by the slowroll consistency relations to second order, with n_{t} = − r(2 − r/ 8 − n_{s}) / 8 and dn_{t}/ dlnk = r(r/ 8 + n_{s} − 1) / 8. Planck measures the running to be dn_{s}/ dlnk = − 0.016 ± 0.010 when tensors are included (see Table 5 and Fig. 4). The constraints on the tensortoscalar ratio are relaxed compared to the case with no running, due to an anticorrelation between r and dn_{s}/ dlnk, as shown in Fig. 4 for Planck+WP+BAO.
Varying both tensors and running, Planck+WP implications for slowroll parameters are ϵ_{V}< 0.015 at 95% CL, , and .
In summary, the Planck data prefer a negative running for the scalar spectral index of order dn_{s}/ dlnk ≈ − 0.015, but at only the 1.5σ significance level. This is for Planck alone and in combination with other astrophysical data sets. Weak statistical evidence for negative values of dn_{s}/ dlnk has been claimed in several previous investigations with the WMAP data and smaller scale CMB data (e.g., Spergel et al. 2003; Peiris et al. 2003; Dunkley et al. 2011; Hinshaw et al. 2013; Hou et al. 2014).
If primordial, negative values for dn_{s}/ dlnk of order 10^{2} would be interesting for the physics of inflation. The running of the scalar spectral index is a key prediction for inflationary models. It is strictly zero for power law inflation, whose fit to Planck was shown to be quite poor in the previous section. Chaotic monomial models with V(φ) ∝ φ^{n} predict dn_{s}/ dlnk ≈ − 8(n + 2) / (4N + n)^{2} ≈ (n_{s} − 1)^{2}, and the same order of magnitude (10^{3}) is quite typical for many slowroll inflationary models, such as natural inflation (Adams et al. 1993) or hilltop inflation (Boubekeur & Lyth 2005). It was pointed out that a large negative running of dn_{s}/ dlnk ≲ − 10^{2} would make it difficult to support the N_{∗} ≈ 50efoldings required from inflation (Easther & Peiris 2006), but this holds only without nonzero derivatives higher than the third order in the inflationary potential. Designing inflationary models that predict a negative running of with an acceptable n_{s} and number of efolds is not impossible, as the case with modulated oscillations in the inflationary potential demonstrates (Kobayashi & Takahashi 2011). This occurs, for instance, in the axion monodromy model when the instanton contribution is taken into account (McAllister et al. 2010), giving the potential (46)
4.4. Open inflation
Most models of inflation predict a nearly flat spatial geometry with small deviations from perfect spatial flatness of  Ω_{K}  ~ 10^{5}. Curvature fluctuations may be regarded as local fluctuations in the spatial curvature, and even in models of inflation where the perturbations are calculated about a spatially flat background, the spatial curvature on the largest scales accessible to observation now are subject to fluctuations from perfect spatial flatness (i.e., Ω_{K} = 0). This prediction for this fluctuation is calculated by simply extrapolating the power law spectrum to the largest scale accessible today, so that Ω_{K} as probed by the CMB roughly represents the local curvature fluctuation averaged over our (causal) horizon volume. Although it has sometimes been claimed that spatial flatness is a firm prediction of inflation, it was realized early on that spatial flatness is not an inexorable consequence of inflation and large amounts of spatial curvature (i.e., large compared to the above prediction) can be introduced in a precise way while retaining all the advantages of inflation (Gott 1982; Gott & Statler 1984) through bubble nucleation by false vacuum decay (Coleman & De Luccia 1980). This proposal gained credence when it was shown how to calculate the perturbations in this model around and beyond the curvature scale (Bucher et al. 1995; Bucher & Turok 1995; Yamamoto et al. 1995; Tanaka & Sasaki 1994). See also Ratra & Peebles (1995, 1994) and Lyth & Stewart (1990). For more refined later calculations see for example Garriga et al. (1998, 1999), Gratton & Turok (1999), and references therein. For predictions of the tensor perturbations see for example Bucher & Cohn (1997), Sasaki et al. (1997), and Hertog & Turok (2000).
An interesting proposal using singular instantons and not requiring a false vacuum may be found in Hawking & Turok (1998), and for calculations of the resulting perturbation spectra see Hertog & Turok (2000) and Gratton et al. (2000). Models of this sort have been studied more recently in the context of the string landscape. (See, for example, Vilenkin 2007, for a nice review.) Although some proposals for universes with positive curvature within the framework of inflation have been put forth (Gratton et al. 2002), it is much harder to obtain a closed universe with a spatial geometry of positive spatial curvature (i.e., Ω_{K}< 0) (Linde 2003).
Theoretically, it is of interest to measure Ω_{K} to an accuracy of approximately 10^{4} or slightly better to test the prediction of simple flat inflation for this observable. A statistically significant positive value would suggest that open inflation, perhaps in the context of the landscape, was at play. A statistically significant negative value could pose difficulties for the inflationary paradigm. For a recent discussion of these questions, see for example Freivogel et al. (2006), Kleban & Schillo (2012), and Guth & Nomura (2012).
Fig. 5 Joint posterior for tensors and running of the scalar spectral index marginalizing over other parameters. The contours are set at 68% and 95%. The red contours apply for Planck+WP+highℓ data. The colour of the scattered points indicates the distribution of Ω_{K}. The blue contours apply when BAO data is also included. (Ω_{K} is then found to be well constrained close to zero.) The dashed vertical line shows the norunning solution. 

Open with DEXTER 
In order to see how much spatial curvature is allowed, we consider a rather general model including the parameters r,n_{s}, and dn_{s}/ dlnk as well as Ω_{K}. We find that with Planck+WP, and Ω_{K} = − 0.004 ± 0.0036 with Planck+WP+BAO. More details can be found by consulting the parameter tables available online^{10}. Figure 5 shows r and dn_{s}/ dlnk for this family of models. We conclude that any possible spatial curvature is small in magnitude even within this general model and that the spatial curvature scale is constrained to lie far beyond the horizon today. Open models predict a tensor spectrum enhanced at small wavenumber k^{<}_{~}1, where k = 1 corresponds to the curvature scale, but our constraint on Ω_{K} and cosmic variance imply that this aspect is likely unobservable.
4.5. Relaxing the assumption of the latetime cosmological concordance model
The joint constraints on n_{s} and r shown in Fig. 1 are one of the central results of this paper. However, they are derived assuming the standard ΛCDM cosmology at late times (i.e., z ≲ 10^{4}). It is therefore natural to ask how robust our conclusions are to changes of the late time cosmological model. We discuss two classes of models: firstly, changes to the ΛCDM energy content; and secondly, a more general reionization model. These extensions can lead to degeneracies of the additional parameters with n_{s} or r^{11}.
Fig. 6 Marginalized joint 68% and 95% CL regions for Planck+WP+BAO data for ΛCDM + r + N_{eff} and ΛCDM + r + Y_{P} (top); and ΛCDM + r + ∑ m_{ν} and ΛCDM + r + w (bottom). Shown for comparison are the corresponding contours for the ΛCDM + r model. 

Open with DEXTER 
4.5.1. Extensions to the energy content
We consider the ΛCDM + r + N_{eff}, ΛCDM + r + Y_{P}, ΛCDM + r + ∑ m_{ν}, and ΛCDM + r + w extensions of the standard model. This selection is motivated by the impact on the CMB damping tail of the first two and the effect on the SachsWolfe plateau at low multipoles for the latter two. The resulting contours are shown in Fig. 6. While the lower limit on n_{s} is stable under all extensions considered here, the models that alter the highℓ part of the spectrum permit significantly bluer spectral tilts, and accordingly also lead to a weaker bound on the tensortoscalar ratio. By allowing N_{eff} to float, we obtain n_{s} = 0.9764 ± 0.0106 and r_{0.002}< 0.15 at 95% CL^{12}. For ΛCDM + r + Y_{P} we obtain n_{s} = 0.9810 ± 0.0111 and r_{0.002}< 0.18 at 95% CL. The models modifying the largescale part of the power spectrum, on the other hand, do not lead to a notable degradation of constraints on either n_{s} or r (n_{s} = 0.9648 ± 0.0061 and r_{0.002}< 0.13 at 95% CL for ΛCDM + r + ∑ m_{ν}, and n_{s} = 0.9601±0.0070 and r_{0.002}< 0.11 at 95% CL for ΛCDM+r+w).
4.5.2. General reionization scenario
In the standard rapid reionization scenario typically used in CMB analysis, the Universe is assumed to be completely transparent after recombination, but the ionization fraction increases from zero to one over a duration Δz ≈ 1 at a certain redshift z_{reion}, which is the only unknown parameter of the reionization model. This model is obviously simplistic, but for CMB analysis it works quite well because the CMB has little sensitivity to the details of how the ionization fraction changes from 0 to 1. In this section we study to what extent allowing more general reionization scenarios may alter some of the conclusions concerning the constraints on n_{s} and r as well as on τ. As discussed in Mortonson & Hu (2008a,b), CMB anisotropies constrain reionization almost entirely by using the shape of the largescale EE power spectrum, and the power is redistributed from larger to smaller scales for reionization processes which take place during a nonnegligible redshift interval, since they start at an earlier epoch.
We use the method developed by Mortonson & Hu (2008a) to describe and constrain the reionization history. A complete principal component basis serves to describe the effect of reionization on the largescale Emode polarization power spectrum. Following Mortonson & Hu (2008a) we bin the ionization history x_{e}(z_{i}) using 95 equal width bins with Δz = 0.25 ranging from z_{min} = 6 to z_{max} = 30. For the redshifts z<z_{min} we assume values for x_{e} which take into account first (and possibly second) helium ionization and complete hydrogen ionization (x_{e} = 1.16 for z< 3 and x_{e} = 1.08 for 3 <z< 6). For z> 30 we fix x_{e} = 2 × 10^{4} as the value of x_{e} expected before reionization (and after primordial recombination). Any reionization history can be parameterized as a free function of redshift by decomposing the ionization fraction as , where the principal components, S_{μ}(z), are the eigenfunctions of the Fisher matrix computed by taking the derivatives of the EE polarization power spectrum with respect to x_{e}(z) and is a fiducial ionization fraction. Following Mortonson & Hu (2008a) we consider here the first five eigenfunctions, S_{μ}(z) with μ = 1,...,5, which will be varied with the other cosmological parameters.
In Fig. 7 we plot the 68% and 95% CL regions for n_{s} and r. The constraint on the tensortoscalar ratio is not significantly affected by this additional marginalization, increasing to r_{0.002}< 0.13 at 95% CL. The scalar spectral index is increased to n_{s} = 0.9650 ± 0.0080, compared with n_{s} = 0.9603 ± 0.0073 obtained with the rapid reionization scheme. This is the same trend as noted in Pandolfi et al. (2010) using WMAP data, but the effect is less significant due to the improved measurement of the temperature spectrum by Planck. The larger freedom in the reionization history increases the width of the posterior on the derived optical depth, which is still partially degenerate with the scalar spectral index of primordial perturbations. The n_{s} = 1 model is still excluded at high significance; however, we find compared to the ΛCDM model.
Fig. 7 Marginalized joint 68% and 95% CL regions for Planck+ WP data for ΛCDM + r for instantaneous and general reionization. 

Open with DEXTER 
5. Inflationary model comparison
In Sect. 4.2 several representative families of parameterized models for the inflationary potential were analysed within the slowroll approximation in the neighbourhood of the pivot scale k_{∗}. Approximate constraints were applied to reject models for which there is no plausible scenario for entropy generation. In this section we revisit some of the parametric models defined in Sect. 4.2. Here, however, the modes for the first order perturbations, as described in Sect. 2.2, are integrated numerically. Thus there is no slowroll approximation, and the issue of the existence of a plausible scenario for entropy generation is examined more carefully. We perform a statistical model comparison between the competing parameterized potentials, both within the framework of Bayesian model comparison and in terms of the relative likelihoods of the best fit models from each parameterized family.
As noted in Sect. 2.3, considerable uncertainty surrounds what occurred during the epoch of entropy generation, particularly with respect to the energy scale at which entropy generation ends and the average equation of state between that epoch and the end of inflation. For this reason, we explore a number of scenarios for this intermediate era during which entropy generation takes place. The models compared in this section include inflation with power law potentials, defined in Eq. (34), with several representative values for n; natural inflation, in Eq. (39); and hilltop inflation, with p = 4 and λ = 4Λ^{4}/μ^{4} in Eq. (37).
The free parameters in these inflationary potentials may vary over several orders of magnitude corresponding to unknown scales in high energy particle physics. Consequently a logarithmic prior is a sensible choice for these parameters. However, there is no theoretical guidance on how to truncate these model priors. We therefore adopt broad priors initially and then truncate them as follows in order to compare all the models on an equal footing.
The strongest constraint on the inflationary parameter space comes from the amplitude of the primordial power spectrum. This is a free parameter in most models, and successful structure formation in a universe dominated by cold dark matter has long been known to require primordial fluctuations with ℛ ≈ 10^{5}, or A_{s} ≈ 10^{10} (see e.g., Zeldovich 1972; Linde 1990). We can therefore immediately reject models for which A_{s} is far from this value, so regions of parameter space which do not yield 10^{11} ≤ A_{s} ≤ 10^{7} are a priori excluded. This range is generous relative to estimates of A_{s} prior to Planck (e.g., Komatsu et al. 2011), but the results do not depend strongly on the range chosen. This effectively truncates the logarithmic priors on the model parameters, leaving a parameter subspace compatible with basic structure formation requirements.
For the single parameter models this requirement defines the range of λ in Eq. (34). However, for generic multiparameter models, an A_{s}based cut may select a nontrival region of parameter space, as happens for the two cases considered here (see e.g., Fig. 1 in Easther & Peiris 2012). Without the A_{s}based cut in the prior, the parameter volume for both natural and hilltop inflation would be rectangular, and the corresponding Bayesian evidence values computed for these models would be lowered.
As discussed in Sect. 2.3, specifying an inflationary potential does not enable us to predict the late time CMB angular power spectra. The subsequent expansion history and details of the epoch of entropy generation are required to relate the value of the inflaton field at Hubble radius crossing to comoving wavenumbers in today’s Universe, through Eq. (24). Physically, the fundamental parameter that sets the observable perturbation spectrum is the value of the field φ_{∗} at which the pivot mode leaves the Hubble radius. It can be rescaled by a shift φ_{∗} → φ_{∗} + φ_{0}, and the range over which φ changes during inflation varies greatly between models. Consequently we treat the remaining number of efolds, N_{∗}, after the pivot scale leaves the Hubble radius as a free parameter with a wide uniform prior, since this quantity has a consistent interpretation across models. The pivot scale used to compute N_{∗} is k_{∗} = 0.05 Mpc. However, given our ignorance concerning the epoch of entropy generation, a multitude of entropy generation scenarios for each inflationary potential can occur. Some possibilities are as follows, with parameters referring to Eq. (24).

1.
Instantaneous entropy generation scenario. l Atthe end of inflation, all the energy in the inflaton field isinstantaneously converted into radiation.

2.
Restrictive entropy generation scenario (narrow range for w_{int}). l GeV, and w_{int} ∈ [ − 1 / 3,1 / 3] .

3.
Permissive entropy generation scenario (wide range for w_{int}). l GeV, and w_{int} ∈ [ − 1 / 3,1] .
The equations of state with w_{int} in the range [1 / 3,1] appear less plausible, but models with these values have been put forward (Pallis 2006), so this possibility cannot be completely excluded. Moreover the w_{int} parameterization captures a variety of scenarios in which the postinflationary Universe is thermalized, but not radiation dominated, including phases of coherent oscillations (Martin & Ringeval 2010; Easther et al. 2011), resonance (Traschen & Brandenberger 1990; Kofman et al. 1994, 1997; Allahverdi et al. 2010), kination (Spokoiny 1993; Chung et al. 2007), secondary or thermal inflation (Lyth & Stewart 1996), moduli domination (Banks et al. 1994; de Carlos et al. 1993), primordial black hole domination (Anantua et al. 2009), or a frustrated cosmic string network (Burgess et al. 2005), all of which lead to an expansion rate differing from that of a radiation dominated universe.
Model priors.
At the other extreme, the decision to exclude w_{int}< − 1 / 3, as done here, is not completely justifiable. We cannot, for example, rule out a first order phase transition at a lower energy scale that would drive w_{int} below −1 / 3, but here we neglect this possibility. Our analysis does not preclude a secondary period of inflation, but does require that the average expansion during the postinflationary regime parameterized by w_{int} should not be inflationary. This caveat should be kept in mind.
For some of the parameterized models, tighter constraints can, in principle, be placed on w_{int}. It has been argued (see e.g., Liddle & Leach 2003) that for the λφ^{4} potential, the uncertainties concerning entropy generation contribute almost no uncertainty in the determination of φ_{∗}. This is because according to the virial theorem, a field sloshing about φ = 0 in a quartic potential has the same average w, namely w = 1 / 3, as the radiation equation of state. More generally, for a potential of the form φ^{n} around the minimum, w_{vir} = ⟨ w ⟩ = (n − 2) / (n + 2) (Turner 1983); therefore, one may argue that w_{int} should be restricted to the interval whose endpoints are 1 / 3 and w_{vir}. This approach was taken by Martin & Ringeval (2010) in obtaining Bayesian constraints on the reheating temperature for monomial potentials from the CMB. However this scenario requires a carefully tuned potential that has approximately a φ^{n} shape, both at large field values and near the origin far below the inflationary scale. Typically, potentials for which V(φ) ~ φ^{n} at large field values can have very different shapes near the origin. Thus, following Easther & Peiris (2012), in this paper we explore a broader range of w_{int} for these models (including the cases above as subsets) in order to obtain datadriven constraints on w_{int}.
In this paper we focus on the three representative scenarios itemized above, referred to hereafter as scenarios (1), (2), and (3). Our algorithm draws a value of N_{∗} and then given the value of ρ_{th}, computes w_{int}, which is a derived parameter. Models for which w_{int} lies outside the specified range of each scenario under consideration are excluded.
Inflationary model comparison results.
The full set of priors for the inflationary physics is given in Table 6. Dimensionful quantities are expressed in units with reduced Planck mass M_{pl} set to unity.
Due to the nontrivial likelihood surfaces and the large dimensionality of the parameter spaces explored in this section, we use ModeCode coupled to MultiNest v3.0^{13} to map out the parameter space. In addition to the standard nested sampling (NS) algorithm, MultiNest v3.0 enables nested importance sampling (NIS), resulting in substantial speed gains^{14} and significant enhancements in the accuracy of the Bayesian evidence computation compared to NS alone for the same computational setup.
5.1. Results
Table 7 presents model comparison results for the ensemble of parameterized potential families described above. We report the Bayesian evidence (model averaged likelihood) ratio, which provides a selfconsistent framework for calculating the betting odds between models (see Sect. 3.4). The uncertainty in these logarithmic evidence values is approximately 0.2. We also report the values computed from the 2lnℒ_{max} values found by the sampler.
The monomial models have a single parameter potential, and the natural and hilltop inflation models have two parameters each. All entropy generation scenarios except case (1) contribute one additional parameter to the inflationary sector. The evidence ratios and values are presented with respect to the ΛCDM cosmological model.
None of the inflationary models tested here fit the data as well as the ΛCDM model. This mostly reflects that there is no evidence in the data for r different from zero. Furthermore, the priors listed in Table 6 for the ΛCDM primordial sector are purely phenomenological, roughly corresponding to ranges somewhat broader than WMAP constraints. Narrowing them around the best fit model arbitrarily increases the evidence in its favour. Instead it is instructive to compare the relative evidence for the inflationary models presented.
Table 7 shows that the λφ^{4} model is decisively ruled out by Planck, confirming previous analyses by the WMAP team (Peiris et al. 2003; Spergel et al. 2007; Dunkley et al. 2009; Komatsu et al. 2011) based on the model track plotted on the n_{s}  r plane. Recent model selection analyses (Martin et al. 2011; Easther & Peiris 2012) with WMAP 7year data found that the model was already disfavoured by odds of about 400:1 against. With Planck, the odds against this model are at least 500 000:1 compared to ΛCDM for a broad range of entropy generation scenarios. The same conclusion is confirmed by the extremely poor values for the model. Given the strength of our results, in the flexible setting of the permissive entropy generation scenario, it is possible not just to rule out models where the potential is of the quartic form in the full range from the origin to the inflationary scales, but also a general class where the potential is of the n = 4 form in the φrange where the cosmological perturbations are generated, but exhibits a different shape near the origin.
Two other large field models, the quadratic potential and natural inflation, are somewhat disfavoured by the Planck data, especially when broader entropy generation scenarios are considered. Compared with the ΛCDM model, these models are disfavoured by –6 depending on the entropy generation scenario. This reflects the analysis of Sect. 4, where the overlap of the model predictions and the data constraints on the n_{s}  r plane is seen to be mostly outside the joint 68% CL contour. However, from the Bayesian evidence point of view, it is too early to declare these models incompatible with the data. To make this judgement, it is more prudent to compare these models to the n = 1 case, which has the best evidence with respect to ΛCDM, rather than to ΛCDM itself, which provides our reference point for the evidence calculation, but has arbitrary prior ranges. In their simplest forms – instantaneous entropy generation – the n = 2 and natural inflation models are only disfavoured by odds of about 1–12:1 against, which does not rise to a high level of significance^{15}.
Fig. 8 Constraints on n_{s} vs. r at k_{∗} = 0.002 Mpc^{1} for the inflationary models considered (i.e., power law potentials with n = 2 / 3,1,2, and 4, natural inflation, and hilltop inflation), showing joint 68% and 95% CL. Blue and grey distributions correspond to the restrictive and permissive entropy generation scenarios, respectively. The instantaneous entropy generation case corresponds to the thin (red) contours in the natural and hilltop panels; for the single parameter models, this case corresponds to the lowestr extremity of the restrictive case. The difference between the natural inflation region in Fig. 1 and the natural inflation constraints shown here is due to the strong projection effect described in the text. 

Open with DEXTER 
Fig. 9 Constraints on n_{s} vs. w_{int} at k_{∗} = 0.002 Mpc^{1} for the inflationary models considered, as in Fig. 8. The instantaneous entropy generation case (1) corresponds to w_{int} = 1 / 3. 

Open with DEXTER 
Fig. 10 Constraints on log _{10}(λ) vs. w_{int} for the inflationary models considered, as in Figs. 8 and 9. 

Open with DEXTER 
The models most compatible with the Planck data in the set considered here are the two interesting axion monodromy potentials, n = 1 (McAllister et al. 2010) and n = 2 / 3 (Silverstein & Westphal 2008), which are motivated by inflationary model building in the context of string theory. The p = 4 hilltop model presents an interesting case. This model was previously found to be compatible with WMAP 7year data, performing almost as well as the monodromy potentials (Easther & Peiris 2012). However it exhibits significant tension with the Planck data, both in terms of evidence ratios and the maximum likelihood. The only exception is the entropy generation scenario (3) which has odds of greater than 1000:1 against compared to ΛCDM, and yet the maximum likelihood is not significantly different from the n = 1 case. This indicates that while the extra freedom allowed by the least restrictive entropy generation scenario improves the best fit, this prior is not very predictive of the data. However, the result seems counterintuitive and merits further comment; we will consider this question further at the end of this section when parameter estimation results are discussed.
We now turn to parameter constraints. Figure 8 presents marginalized joint constraints from Planck+WP alone on the derived parameters n_{s,0.002} and r_{0.002}. Figure 9 shows the corresponding joint constraints on n_{s,0.002} and w_{int}, again a derived parameter. Figure 10 shows joint constraints on w_{int} and the potential parameter log _{10}(λ) (log _{10}(Λ) in the case of natural inflation). It is instructive to consider the three sets of figures together. The restrictive and permissive entropy generation scenarios are shown on all panels; the instantaneous entropy generation case is shown for the two parameter models, natural inflation and hilltop, in Fig. 8 – for the monomial potentials this case corresponds to the lowestr extremity of the restrictive case constraints, and to w_{int} = 1 / 3 in the other two figures.
The quartic potential conflicts with the data because it predicts a high tensortoscalar ratio. Hence the model maximizes its likelihood by pushing towards the lowerr, bluern_{s} limits of its parameter space, which corresponds to increasing w_{int} as much as allowed by the entropy generation prior. The contours terminate at the lowestr limit when each entropy generation case hits its w_{int} prior upper limit (i.e., w_{int} = 1 / 3 and w_{int} = 1 for restrictive and permissive entropy generation, respectively). For each case, the lower limit on w_{int}, corresponding to the reddestn_{s} extremity of the confidence contours, is datadriven. The quadratic potential encounters the same difficulty but at a less extreme level.
The two axion monodromy potentials are compatible with a wide range of entropy generation scenarios. The instantaneous entropy generation scenario is compatible with the data for both models. For restrictive entropy generation in the n = 1 model, we obtain a datadriven upper limit on w_{int}, which just touches the w_{int} = 1 / 3 case at the 95% CL. At the lower limit, the w_{int} posterior is truncated by the prior, as for the n = 2 / 3 case. For the latter, there is a datadriven upper limit on w_{int} which is controlled by the upper limit on n_{s}. For permissive entropy generation, the upper and lower limits on w_{int} for both models are datadriven, corresponding to the upper and lower limits on n_{s}, respectively.
The constraints on natural inflation require some interpretation. The relationship between the empirical n_{s} and r parameters and the potential parameters for natural inflation is discussed in detail by Savage et al. (2006) and Mortonson et al. (2011) along with parameter constraints derived from WMAP 3year (Spergel et al. 2007) and WMAP 7year (Larson et al. 2011), respectively. In this model, there is a degeneracy between f and Λ in the limit where these parameters are large so that natural inflation resembles the quadratic model. The priors are chosen to exclude most of this region. The priors on log _{10}(f) and log _{10}(Λ) still allow a region of nearly degenerate models that contribute to a ridge seen in the natural inflation panel of Fig. 8. These models closely match the values of n_{s} and r seen in the quadratic potential constraints. The marginalized constraints on n_{s} and r depend strongly on the prior on log _{10}(f) due to the projection of a large number of degenerate models onto this ridge. Therefore the apparent preference for this region of parameter space over models with lower values of r is largely due to this effect and is not driven by the data. This highly nonlinear mapping between the logarithmic priors on the potential parameters and the power law parameters (which are derived parameters in this analysis) leads to a strong projection effect, which accounts for the difference in visual appearance between these contours and the region labelled natural inflation in Fig. 1.
Generally, for fixed N_{∗}, decreasing Λ and f reduces both n_{s} and r. Thus natural inflation models can have lower values of r than the quadratic potential without increasing n_{s} and N_{∗}. This feature means that the potential parameters for this model are relatively uncorrelated with w_{int}, in contrast with the other models considered here, as illustrated in the Fig. 10. Nevertheless we obtain datadriven bounds on w_{int} in the permissive entropy generation case as well as a lower bound in the restrictive entropy generation case. Both bounds overlap with the instantaneous generation limit.
The upper panel of Fig. 11 shows our lower limit of f ≳ 10.0M_{pl} (95% CL), compared with the WMAP7 limit f ≳ 5.0M_{pl} (95% CL) reported by Mortonson et al. (2011). Indeed, the Planck limit is in agreement with the Planck prediction presented in that work. There is a hint of an upper limit on f as well, driven by the fact that this corresponds to the quadratic inflation limit, which is in tension with the data. However this is only a 1σ effect.
The p = 4 hilltop model has two distinct branches: a smallfield scenario, where r ≲ 0.001 and n_{s}< 0.95; and a largefield limit in which V(φ) ~ φ (Adshead et al. 2011). Physically, the smallfield limit is consistent with the Lyth bound, and we can select it by fixing log _{10}(Λ) < − 2.5 in the prior. We observe that the data select out this small field branch, which requires explanation given that we know that the n = 1 model (the limit of the large field branch) is perfectly compatible with the data. In fact, this can be tested by restricting the log _{10}(Λ) prior by hand to the largefield branch. The hilltop model constitutes a difficult sampling problem, as is apparent from Fig. 4 in Easther & Peiris (2012). An examination of the progress of the nested sampler for the case of the full hilltop prior reveals the reason for the smallfield branch being selected out. In this case, the posterior for the largefield branch is extremely thin compared to the model prior in this regime – much thinner than the posterior for the smallfield branch in comparison to its respective prior. Therefore this region occupies very little probability mass, and is dropped in preference to the more predictive smallfield branch. This high likelihood but extremely thin ridge is also responsible for the counterintuitive result reported in the model selection analysis, where the model was found to have a good with respect to ΛCDM, and yet be highly disfavoured by the Bayesian evidence. In summary, confirming the results of Sect. 4, this model is in agreement with the data in the limit where it overlaps the linear n = 1 model. But since this region of high likelihood occupies a very small fraction of the prior, this model is heavily penalized by the Bayesian evidence for failing to predict the data over most of its prior space.
Fig. 11 Potential parameters for natural inflation and hilltop inflation, as in Figs. 8–10. On the natural inflation panel, instantaneous entropy generation corresponds to a thin diagonal along the top edge of entropy generation case (3). 

Open with DEXTER 
6. Observable window of inflation
Section 4.2 presented an analysis of several representative inflationary potentials within the framework of the slowroll approximation and their compatibility with the Planck data. The results are summarized in Fig. 1. In that case the full potential is considered in order to identify a plausible range for the location of φ_{∗} on the potential V(φ). This requires a complete story. In other words, the potential must be specified starting above the point where the largest observable scales first exited the Hubble radius, and extending to the minimum of the potential.
In this section we explore another approach. We adopt the point of view that we are interested in reconstructing the inflationary potential only over the observable range – that is, the interval of φ corresponding to the scales observable today in the CMB. We constrain the potential over the range where these scales exited the Hubble radius during inflation as well as a few efolds before and after. The cosmological perturbations are not imprinted instantaneously at the moment of Hubble radius crossing, but rather gradually over a few efolds. We expand around φ_{∗}, taking the view that a plausible extension of the potential outside this observable range is always possible, so that one has precisely the number of efolds of inflation needed for φ_{∗} to correspond to k_{∗} today.
The argument is that one can always end inflation abruptly by imposing a sharp waterfall feature where needed, or prolong inflation by inserting a sufficiently long plateau into the potential by hand, for example, for models with a large tilt. A foreseeable objection to this approach is that the extensions of the potential required outside the observable window may render the potential unnatural. This possibility should be kept in mind, although naturalness is an elusive and uncomfortably subjective concept. The analysis in this section does not rely on the slowroll approximation. Instead each k mode is integrated exactly by numerical integration, as described in Sect. 2.2, under the assumption of a canonical kinetic term.
Two complementary approaches to reconstruct the potential have been explored in the literature. The first approach, followed in this paper, expands the potential V(φ) directly in powers of (φ − φ_{∗}). In this case the numerical integration of the slowroll solution must start sufficiently early so that any initial transient has had a chance to decay, and one is in the attractor solution when the dynamics of the largest observable modes in the Universe today start to have an interesting evolution. This is the approach followed, e.g., by Lesgourgues & Valkenburg (2007) and Mortonson et al. (2011) using publicly available codes^{16}.
A second approach expands H(φ) as a Taylor series in (φ − φ_{∗}). As discussed in Sect. 2.2, this has the advantage that H(φ) determines both the potential V(φ) and the solution φ(t), so the issue of having to start sufficiently early in order to allow the initial transient to decay is avoided. This method was used, for example, in Kinney (2002), Kinney et al. (2006); Peiris & Easther (2006a,b, 2008), and Easther & Peiris (2006) using analytic and semianalytic approximations, and in Lesgourgues et al. (2008), Powell & Kinney (2007), Hamann et al. (2008c), and Norena et al. (2012) using a fully numerical approach.
These approaches could lead to results that differ from those in Sect. 4.2. Firstly, if the running of the index is large, the slowroll approximation taken to second order is not necessarily accurate for all models allowed by the data. The relation between the spectral parameters (lnA_{s}, n_{s}, dn_{s}/ dlnk, r) and the underlying inflationary potential V(φ) is therefore uncertain. Secondly, for spectra with a large running, there is no guarantee that an inflationary model giving such a spectrum exists. All allowed models have r ≪ 1, so these models are consistent with ϵ_{V}(k_{∗}) ≪ 1 at the pivot scale. However, towards the edge of the observable range, the potential may become incompatible with ϵ_{V}(k) < 1 (i.e., with the requirement of inflationary expansion). These possible pitfalls are avoided using the methods in this section, since the data have been fit directly by the candidate V(φ) or H(φ), computed numerically without any slow roll assumptions over the entire observable range.
Fig. 12 Posterior distribution for the first three slowroll potential parameters using Planck+WP data. In the n = 2 − 4 cases, the inflaton potential is expanded to nth order, and the spectrum is obtained by fitting the numerically computed power spectrum to the data, with no slowroll approximation and no assumption about the extrapolation of the potential outside the observable window. 

Open with DEXTER 
Numerical reconstruction of potential parameters, compared to results with the slowroll approximation, when tensors and running are included (Planck+WP 95% CL, with k_{∗} = 0.05 Mpc^{1}).
Fig. 13 Posterior distribution for the Taylor expansion coefficients, V_{i}, of the inflaton potential. The potential is expanded to nth order, assuming a flat prior on ϵ_{V}, η_{V}, , and . The coefficients V_{i} are expressed in natural units (where The contours show only half of the allowed regions for potential parameters: the other half is symmetric, with opposite signs for V_{1} and V_{3}. 

Open with DEXTER 
Fig. 14 Observable range of the best fitting inflaton potentials, when V(φ) is Taylor expanded to nth order around the pivot value φ_{∗}, in natural units (where assuming a flat prior on ϵ_{V}, η_{V}, , and , and using Planck+WP data. Potentials obtained under the transformation (φ − φ_{∗}) → (φ_{∗} − φ) leave the same observable signature and are also allowed. The sparsity of potentials with a small V_{0} = V(φ_{∗}) is explained by the flat prior on ϵ_{V} rather than on ln(V_{0}). In fact, V_{0} is unbounded from below. 

Open with DEXTER 
We define a class of models over the observable range based on the expectation that the potential should be smooth. V(φ) is approximated by a Taylor expansion up to order n, and we explore the cases n = 2,3, and 4. For each V(φ), we integrate over inflationary fluctuations using the inflation module implemented in CLASS^{17} (Lesgourgues 2011; Blas et al. 2011) as described in Lesgourgues & Valkenburg (2007). Potentials are rejected for which the attractor solution cannot be reached when the largest observable scales cross the Hubble radius. The parameters sampled are the potential and its derivatives at the pivot scale when k_{∗} crosses the Hubble radius during inflation.
To avoid parameter degeneracies, we impose uniform priors on and at the pivot field value φ_{∗}. The advantage of uniform priors on these parameters is that – to the extent that the slowroll conditions are satisfied – these coefficients relate linearly to observable quantities such as n_{s}, r, dn_{s}/ dlnk, and d^{2}n_{s}/ dlnk^{2}. Figure 12 and Table 8 show the posterior probability for these coefficients, and Fig. 13 shows the posterior probability for the Taylor series coefficients V_{i}. In Fig. 14, we show the observable range of the best fitting inflaton potentials (for a sample extracted randomly from the converged Markov chains). The edges of the observable range correspond to Hubble crossing for the minimum and maximum values of k used in the Boltzmann code. We stress here that the Planck data suggest a flat potential when the lowest order slowroll primordial spectra are considered, as analysed in Sect. 4.2. However, when the restrictions to lowest order in the slowroll approximation are relaxed, the inflaton potential can differ markedly from a plateaulike potential, as the green curves in Fig. 14 show.
When fitting V(φ) for each model in parameter space, we compute (n_{s}, dn_{s}/dlnk, d/dlnk^{2}, and r) at the pivot scale a posteriori directly from the numerical primordial spectra. The results are shown in Table 9 and can be compared to those of Table 5 in Sect. 4. We can also use the results of Sect. 4 for the ΛCDM + r + dn_{s}/ dlnk model to infer the potential parameters (ϵ_{V}, η_{V}, and ) using the secondorder slowroll expressions, and compare different approaches in the space of potential parameters (see Fig. 12 and the last column in Table 8).
Scalar amplitude, tilt, running, running of the running, and tensortoscalar ratio inferred from a numerical reconstruction of the inflaton potential (Planck+WP 95% CL, with k_{∗} = 0.05 Mpc^{1}).
The model with a quadratic potential in the observable window (n = 2) leads to bounds on ϵ_{V}, η_{V}, n_{s}, and r very close to the ΛCDM + r case. This is not a surprise since such potentials cannot give values of n_{s} and r compatible with the data and at the same time a large running. A significant dn_{s}/ dlnk can be generated only in the presence of a large (i.e., with a significant V^{′′′}). Since quadratic potentials produce little running, they are faithfully described by the slowroll approximation.
The model with a cubic term (n = 3) has the freedom to generate a large running, dn_{s}/ dlnk. Indeed one can check that the results for the n = 3 model are close to those of the ΛCDM + r + dn_{s}/ dlnk model presented in Sect. 4. The agreement between these two models remains very good, despite the fact that in the presence of a large running, the slowroll approximation can become inaccurate. The running in a potential with n = 3 is not exactly scale invariant; this is not captured by the ΛCDM + r + dn_{s}/ dlnk parameterization.
The n = 4 model has even more freedom, allowing a considerable running of the running d^{2}n_{s}/ dlnk^{2} (to the extent that inflation holds during the observable efolds). In that case, the spectrum is better fitted when the two parameters r and d^{2}n_{s}/ dlnk^{2} are nonzero. In Fig. 14, we see that most n = 4 potentials have a long and steep tail for φ<φ_{∗}, with a kink around φ_{∗} − 0.4 (in natural units). This shape generates a significant running on the largest observable scales, while preserving a smaller running on smaller scales. With such a feature in the scalar primordial spectrum at large scales combined with a nonzero contribution from tensor fluctuations, the best fit model for n = 4 has a temperature spectrum very close to that of the minimal ΛCDM model for ℓ > 40, but not for smaller multipoles. The amplitude of the SachsWolfe plateau is smaller. This allows the largescale data points from Planck to be fitted slightly better. However, the case r = dn_{s}/ dlnk = 0 still lies at the edge of the 95% CL, and the minimum effective χ^{2} of this model is smaller than in the n = 2 case by only 3.7.
A comparison of the n = 3 and n = 4 results clearly shows that the process of expanding the inflaton potential to various orders and fitting it to the data does not converge (at least not by n = 4). Given the 1–2σ preference of Planck data for a nonzero running and running of the running, we find that a modelindependent reconstruction of the inflaton potential is not possible under the assumptions of this section. In other words, as long as we assume that V(φ) can be described during and after observable inflation by a polynomial of order 2 or 3, we can put strong bounds on ϵ_{V} and η_{V}. But if we introduce more derivatives to describe the observable part of the potential and allow complete freedom to extrapolate V(φ) outside this region, the constraints can be easily evaded.
7. Primordial power spectrum reconstruction
In this section we report on a search for features in the primordial power spectrum. In the basic six parameter model studied in the companion Planck paper Planck Collaboration XVI (2014), the primordial power spectrum , which includes only the adiabatic mode, is modelled using the power law , for which the best fit values are A_{s} = 2.20 × 10^{9} and n_{s} = 0.9603 for a pivot scale k_{∗} = 0.05 Mpc^{1}. An extension of this parameterization is also considered allowing for a running of the spectral index (dn_{s}/ dlnk = − 0.013 ± 0.009). But in all cases considered it was assumed that the power spectrum is smooth and without bumps, sharp features, or wiggles. In this section we investigate whether any statistically significant evidence for features is present in the data when these assumptions are relaxed. Allowing an arbitrary function for the input power spectrum is not an option because in this case the recovered primordial power spectrum is dominated by smallscale noise. Instead we consider here a penalized likelihood approach where a preference for smooth power spectra is imposed. Section 8 of this paper pursues a complementary approach where several parametric models for wiggles and features are explored to see whether a statistically significantly better fit can be obtained.
An extensive literature exists on how to search for features in the power spectrum using a wide range of methods. The following papers and the references therein provide a sampling of the literature on nonparametric reconstruction: RichardsonLucy deconvolution (Lucy 1974; Richardson 1972; Hamann et al. 2010; Shafieloo & Souradeep 2004, 2008), deconvolution (TocchiniValentini et al. 2005, 2006; Ichiki & Nagata 2009; Nagata & Yokoyama 2008, 2009), smoothing splines (Verde & Peiris 2008; Peiris & Verde 2010; Sealfon et al. 2005; Gauthier & Bucher 2012), linear interpolation (Hannestad 2004; Bridle et al. 2003), and Bayesian model selection (Bridges et al. 2009; Vázquez et al. 2012). The approach pursued here follows TocchiniValentini et al. (2006) and Gauthier & Bucher (2012) most closely. More technical details and extensive tests validating the method can be found in the latter reference.
Let be the best fit power spectrum of the six parameter model. We define a general Ansatz for the power spectrum in terms of a fractional variation, f(k), relative to this fiducial model, so that (47)Any features are then described in terms of f(k).
In this analysis we use the Planck+WP likelihood supplemented by the following roughness penalty or prior, which is added to −2lnℒ: (48)where κ = ln(k/ Mpc^{1}) and κ_{min} and κ_{max} delimit the scales probed by the data. The first regularization term penalizes any deviation from a straight line of the function f(κ). The second and third terms drive the f(κ) to zero where there are effectively no constraints from the data. The value of λ controls the smoothness of the reconstruction, but the precise value of α is less important. It must be large enough to force f(k) towards zero when κ<κ_{min} and κ>κ_{max} but not so large as to render the matrices illconditioned. We use α = 10^{4}.
We represent f(k) using a cubic Bspline on a grid of points in kspace uniformly spaced in κ with step size Δκ = 0.025 and extended from κ = − 12.5 to κ = − 0.3 giving us a total of 485 knots so that . The density of grid points is sufficiently large so that artefacts near the scale of the knot spacing are suppressed for the values of λ used here. Given the large number of dimensions, it is not practical to explore the likelihood using MCMC methods. However, for the power spectrum parameters, the predicted C_{ℓ}s are related by a linear transformation given fixed cosmological parameters, allowing us (for fixed cosmological and nuisance parameters) to find the maximum likelihood solution using the NewtonRaphson method^{18}. We define (49)where the vector Θ represents the cosmological parameters unrelated to the power spectrum and the foreground nuisance parameters. We first minimize over f_{i} using NewtonRaphson iteration and then in an outer loop minimize over Θ using the downhill simplex algorithm. To carry out this procedure, the Planck likelihood code was modified to compute the gradient and Hessian of the likelihood with respect to the C_{ℓ}s.
The cosmological Boltzmann solver CAMB was modified to accept the vector of primordial power spectrum knots f. By default CAMB calculates the C_{ℓ}s for a subset of ℓ and interpolates to obtain the full multipole power spectrum. Instead we calculate the C_{ℓ}s at each ℓ explicitly.
The boundaries κ_{min} and κ_{max} defining where f(κ) is allowed to differ from zero are chosen to match the range of ℓ constrained by the highℓ likelihood. The likelihood includes C_{ℓ}s between ℓ = 50 and ℓ = 2500, which roughly corresponds to k ∈ [0.003,0.2] Mpc^{1}. The lowℓ likelihood covers ℓ = 2 to ℓ = 49, which roughly corresponds to k ∈ [10^{4},0.003] Mpc^{1}. In this range of ℓ, cosmic variance is large, making feature detection difficult. Calculating the gradient and Hessian of a pixelbased lowℓ likelihood is computationally time consuming. We therefore use the lowℓ likelihood only to constrain the cosmological parameters. We choose k_{min} = 0.005 Mpc^{1} and k_{max} = 0.3 Mpc^{1}. Within this k range, variations in the C_{ℓ}s due to the f_{i} are too small to affect the overall likelihood through the lowℓ likelihood. We observed that the difference in the lowℓ likelihood between the reconstructed f(k) and f(k) = 0 is small (<1%) compared to the difference in the highℓ likelihood.
The cosmological parameters τ and A_{s} are almost completely degenerate for the temperature anisotropy except at very low ℓ, so we fix τ to its best fit value for the fiducial model. The likelihood contains additional nuisance parameters that model foreground components and beam shapes, as discussed in Planck Collaboration XV (2014). Many of the nuisance parameters, it can be argued, are unlikely to introduce spurious smallscale structure because they represent foreground models with a power law and thus smooth angular power spectrum. However some nuisance parameters, in particular those describing beam uncertainties, could conceivably introduce artefacts into the reconstruction. Unfortunately, converging to the correct maximum likelihood reconstruction with all the beam shape parameters included is prohibitively time consuming. Therefore we fix the nuisance parameters to their fiducial best fit values, leaving a more detailed examination of this issue to future work.
We found that simultaneously allowing extra degrees of freedom for smallscale structure and variations in the cosmological parameters changes the best fit fiducial model – that is, A_{s} and n_{s} – so that the variations with respect to the fiducial model no longer visibly gave the best straight line fit. Therefore we allow A_{s} and n_{s} to vary, so that the fiducial model is indeed the best straight line fit through the plotted data points. Detailed investigation showed that neither the priors, nor lowℓ, nor highℓ data play a significant role in determining the best fit fiducial model. This effect is small and within the error bars for A_{s} and n_{s} established assuming the fiducial model. To summarize, we maximize the likelihood with respect to the control points f_{i} and the three cosmological parameters h, Ω_{c}h^{2}, and Ω_{b}h^{2}. We then update the fiducial model (A_{s} and n_{s}) at each iteration by finding the best fit power law through the current best fit reconstruction.
Once the maximum likelihood solution has been found, the second derivatives about this solution are readily calculated by extracting the relevant matrices for most of the components, and estimating the remaining components using finite differences. The second derivative matrix is used to estimate the error on the reconstructed f_{i} and the three cosmological parameters h, Ω_{c}h^{2}, and Ω_{b}h^{2}. Monte Carlo simulations of a fiducial data set with a simplified CMB likelihood including some of the nonGaussianities^{19} suggest that this approximation of the error is accurate for λ^{>}_{~}10^{3}.
Figure 15 summarizes our results, showing the estimated f(k) in bins and the corresponding 1σ and 2σ errors. Errors in k are also shown to represent the minimum reconstructible width evaluated at the middle of each bin. This is the minimum width that a Gaussian feature must have to be reconstructed with a small enough bias such that the mean reconstruction differs by less than 10% rms. The minimum reconstructible width is closely related to the correlation length, so that the errors between adjacent bins are weakly correlated and the total number of bins represents roughly the effective number of independent degrees of freedom.
While the plots with a significant roughness penalty – that is, with λ = 10^{6} and λ = 10^{5} – do not show any statistically significant evidence for features standing out above the noise of the reconstruction, for a smaller roughness penalty – that is, for λ = 10^{4} and λ = 10^{3} – a nominally statistically significant feature is clearly visible around k ≈ 0.13 Mpc^{1}. We do not understand the origin of this feature, which may be primordial or may arise as a foreground or other systematic error in the highℓ portion of the likelihood. It should be noted that most of the robustness tests described in the likelihood paper assume smooth power spectra. The maximum excursions are locally at 3.2σ and 3.9σ for λ = 10^{4} and λ = 10^{3}, respectively. In each of these two cases we correct for the look elsewhere effect by calculating the probability that one of the plotted error bars deviates by the same number of or more standard deviations. This calculation is carried out using the covariance matrix of the plotted error bars. We obtain p = 1.74% and p = 0.21%, which corresponds to 2.4σ and 3.1σ, respectively. Additional simulations were carried out to validate the method by generating mock data according to the fiducial model and measuring the errors of the reconstruction obtained. These investigations confirm the error model. These tests were carried out both with and without test features. It can be argued that foregrounds are unlikely to explain the observed feature because all the foreground models involve smooth power law templates, whereas this feature is localized in multipole number. It is important to assess by means of a more extensive set of simulations whether the statistical significance assigned to this result is accurate.
We investigate which CMB angular multipoles correspond to this apparent feature. Figure 16 shows the C_{ℓ} residual from the reconstructed power spectrum with the best fit power law power spectrum subtracted together with the data for each of the frequency map correlation combinations used in the CamSpec likelihood. We observe a smooth dip around ℓ ≈ 1800, which is significant compared to the error bars, in particular for the 217 GHz map. To determine whether this dip is in fact responsible for large deviation in the reconstruction, we take the λ = 10^{3} best fit reconstruction and set f(k) = 0 everywhere except for 0.1Mpc^{1}<k< 0.15Mpc^{1} – the region where the large deviation is located – and calculate the corresponding C_{ℓ} spectrum. Figure 17 plots the C_{ℓ} residuals of this test feature, which show a large dip at around ℓ ≈ 1800, thus demonstrating that the dip in the C_{ℓ} residual of the data centred at ℓ ≈ 1800 is responsible for the large excursions in the primordial power spectrum reconstructions.
Fig. 15 Planck primordial power spectrum feature search results. Top four panels: the reconstructed power spectrum at four values for the smoothing parameter λ. The red curves indicate the maximum likelihood configuration for the fractional deviation f(k) of the power spectrum relative to a power law fiducial model (with A_{s} = 2.20 × 10^{9} and n_{s} = 0.9603) for the penalized likelihood. The error bars have a width corresponding to the minimum reconstructible width (the minimum width for a Gaussian feature so that the mean square deviation of the expectation value of the reconstruction differs by less than 10%) and a vertical extent showing the 1σ and 2σ limits for the fractional deviation averaged over the box. The grey hashed regions at the far left and right show where the fixing prior (i.e., α) sets f(k) = 0. The inner grey regions show where the reconstruction bias is so great that the minimum reconstructible width is undefined. Mock features in this region produce reverberations over the entire interval. With λ = 10^{3} and λ = 10^{4}, we find statistically significant fluctuations around k ~ 0.1 Mpc^{1}. Lower panels: the 1σ error bars for three combinations of cosmological parameters at the four values of λ. The maximum likelihood value for the fiducial model is indicated by the dashed line for comparison. 

Open with DEXTER 
Fig. 16 CMB multipole spectrum residuals for best fit primordial power spectrum reconstruction with smoothing parameter λ = 10^{3}. The panels show the C_{ℓ} spectrum residuals (compared to the best fit power law fiducial model represented by the horizontal straight dashed line) for the four auto and crossspectra included in the highℓ likelihood. Here The data points have been binned with Δℓ = 31 and foregrounds subtracted according to the best fit foreground parameters. The solid black line shows the CMB spectrum residual for the maximum likelihood primordial power spectrum reconstruction with λ = 10^{3}. 

Open with DEXTER 
Note added: the broad dip around ℓ = 1800 in the temperature power spectrum, seen in the 217 GHz channel, has been shown to result from residuals that were strongest in the first survey. In work done after submission of this paper, this feature was shown to be associated with imperfectly subtracted electromagnetic interference generated by the drive electronics of the 4 K cooler and picked up by the detector readout electronics. In a recent study, more aggressive measures were applied to remove 4 K contaminated data. When this censured data was propagated all the way to the maps and the power spectrum, the amplitude of the feature was lowered to below the noise for the first survey. A more complete account and analysis will appear in the next round of Planck cosmology papers.
8. Parametric searches for primordial power spectrum features
In this section we continue to investigate deviations of the primordial power spectrum from a smooth, featureless function, in this case by testing a set of theoretically motivated models.
8.1. Models and priors
We consider three models describing features in the primordial power spectrum: adding a global oscillation, a localized oscillation, or a cutoff to the largescale power spectrum.
8.1.1. Wiggles model
Due to the exponential growth of the scale factor during inflation, a periodically recurring event in proper time which affects the amplitude of curvature perturbations would produce features that are periodic in lnk. This occurs, for instance, for non BunchDavies initial conditions (Easther et al. 2001; Danielsson 2002; Martin & Brandenberger 2003; Bozza et al. 2003), or, e.g., in the axion monodromy model (Silverstein & Westphal 2008), as a consequence of instanton induced corrections to the potential (Flauger et al. 2010). In these scenarios the primordial spectrum has an oscillation superimposed on an underlying smooth spectrum.
Fig. 17 CMB multipole spectrum residual for the primordial power spectrum test feature. The test feature (top) is set to the anomalously large deviation of the primordial power spectrum reconstruction for λ = 10^{3} in the interval 0.1Mpc^{1}<k< 0.15Mpc^{1}, and is zero elsewhere. Bottom: the angular spectrum corresponding to this feature. We observe a large dip at ℓ ≈ 1800. 

Open with DEXTER 
Prior ranges imposed for the wiggles, step inflation, and cutoff model parameters.
Here we consider the following parameterization of the primordial spectrum (referred to as the wiggles model): (50)with amplitude α_{w}, frequency ω, and phase ϕ to quantify the superimposed oscillations. The underlying smooth spectrum has the standard power law form (51)The prior ranges for the wiggles model parameters are given in Table 10. The obvious prior for the phase ϕ is uniform over the interval (0,2π). We also choose a uniform prior on α_{w} (a logarithmic prior on α_{w} introduces considerable dependence of the resulting marginalized posteriors on the lower limit and does not contain the smooth spectrum as a limiting case). The sensitivity to primordial wiggles is limited at high frequencies by the width of the transfer function (Hamann et al. 2008a) and at low frequencies by the requirement of at least one full oscillation in the observable part of the power spectrum. Since Planck data are sensitive to wavenumbers over a range of roughly four orders of magnitude, this condition implies ω ≳ 0.5. Here we restrict the analysis to ω< 100 and assume a uniform prior. Larger values of the frequency are theoretically possible, e.g., in axion monodromy models (Flauger et al. 2010), but the amplitude of the oscillations in the C_{ℓ}s will be suppressed with respect to the primordial one. A comprehensive search for higher frequency oscillations is currently underway.
8.1.2. Step inflation model
If the slow roll of inflation is briefly interrupted, for instance by a phase transition (Starobinsky 1992; Hunt & Sarkar 2004), a burst of resonant particle production (Chung et al. 2000), or a step in the inflaton potential (Adams et al. 2001), or if the speed of sound changes suddenly (Achúcarro et al. 2011), a localized oscillatory feature is superimposed on the scalar primordial power spectrum. We adopt the approximate parameterization for such a feature from a step in the potential, introduced by Adshead et al. (2012), with (52)where (53)As in the wiggles model, we choose a uniform prior on the amplitude parameter (see Table 10). The parameter η_{f} determines both the frequency of the feature and its location, which is required to lie in the observable range. The damping envelope of the feature is set by the ratio η_{f}/x_{d}. We impose uniform priors on the logarithms of η_{f} and x_{d}.
8.1.3. Cutoff model
A number of models have been suggested to explain the apparent lack of power in the quadrupole and octupole of the WMAP temperature power spectrum. Typically in these models, the onset of a slowroll phase coincides with the time when the largest observable scales exited the Hubble radius during inflation. This naturally suppresses the primordial power spectrum at large scales (see, e.g., Sinha & Souradeep 2006). We consider a phenomenological parameterization of a cutoff proposed in Contaldi et al. (2003), given by (54)We apply uniform priors on λ_{c}, which determines the steepness of the cutoff, and on the logarithm of the cutoff scale k_{c}.
8.2. Method
To achieve the necessary numerical precision for models with features in the primordial spectra, we modify the standard settings of the CAMB numerical code in order to calculate C_{ℓ} at each ℓ rather than interpolating and refine the grid in wavenumber for the numerical integration. These changes significantly slow down the computation. In the models considered here, the likelihood function has characteristics that make sampling difficult, such as extended plateaus and multiple isolated maxima, which render the MetropolisHastings algorithm inefficient. We therefore use the nested sampling algorithm implemented in the MultiNest addon (Feroz & Hobson 2008; Feroz et al. 2009) to CosmoMC (Lewis & Bridle 2002), which can also calculate the Bayesian evidence and the likelihood profiles.
Improvement in fit and logarithm of the Bayes factor B_{0X} with respect to power law ΛCDM and best fit parameter values for the wiggles, step inflation, and cutoff models (the larger lnB_{0X}, the greater the preference for a featureless power law spectrum).
The signatures of the feature models considered here are unique and cannot be mimicked by varying other parameters, which lead to smooth variations of the power spectrum (with the exception of highly tuned very low frequency oscillations that can change the acoustic peak structure). We thus restrict ourselves to varying only the parameters describing the features and keep all remaining cosmological and nuisance parameters fixed to their ΛCDM best fit values^{20}.
Fig. 18 Top: best fit primordial spectrum of curvature perturbations for the power law (black), wiggles (red), step inflation (green), and cutoff (blue) models. Centre: residuals of the temperature angular power spectrum. Note that the scale of the vertical axis changes at ℓ = 50. Bottom: zoom of region around the first acoustic peak. 

Open with DEXTER 
Fig. 19 Marginalized posterior probability (red) and profile −2Δlnℒ_{max} (black) for selected parameters of the wiggles (top), step inflation (middle), and cutoff model (bottom). 

Open with DEXTER 
8.3. Results
For all three models we find that including these additional features improves the quality of the fit with respect to a pure power law spectrum. For the Planck+WP data, we show the best fit primordial curvature power spectra and temperature angular power spectrum residuals in Fig. 18, and report the best fit parameter values in Table 11. Since in all three cases the likelihood functions do not tend to zero in all directions of the respective parameter spaces, the Bayesian quantities (i.e., posterior distributions and Bayes factors) depend strongly on the choice of prior. For this reason, we also quote two priorindependent quantities, the effective χ^{2} (i.e., and the profile −2Δlnℒ_{max} as a function of selected model parameters plotted alongside the marginalized posteriors in Fig. 19, which illustrates the unconventional shape of the likelihood functions.
For the wiggles model, oscillations around the first acoustic peak and in the 700 <ℓ < 900 range improve the fit to the data, whereas for the best fit step inflation model the spectrum between the SachsWolfe plateau and the first acoustic peak is fit better. Quantitatively, the cutoff model improves the fit only modestly, with , but both the wiggles and step inflation models lead to a larger improvement, with , at the cost of three new parameters. Already for prePlanck data, improvements of have been reported in related analyses (e.g., Peiris et al. 2003; Martin & Ringeval 2004; Elgarøy et al. 2003; Covi et al. 2006; Meerburg et al. 2012; Benetti et al. 2013; Peiris et al. 2013). Note that in the step inflation model, the best fit does not coincide with the maximum of the marginalized posterior probability, indicating that some degree of fine tuning is necessary to reach the maximum of the likelihood. The maximum of the marginalized posterior at ln(η_{f}/ Mpc) ≈ 7.2 actually reproduces the feature at ℓ ≈ 20 − 40 found previously in WMAP data (Peiris et al. 2003). The secondary peak at ln(η_{f}/ Mpc) ≈ 4 corresponds to a feature at multipoles ℓ ≈ 1800, where the analysis of Sect. 7 found a feature. However the model does not account for this feature well, yielding an improvement of only .
Whether or not these findings can be considered statistically significant or arise simply from overfitting noisy data is not a trivial question (see, for instance, the discussion in Bennett et al. 2011). From a frequentist statistics point of view, an answer would require the rather involved procedure of repeating the analysis on a large set of simulations. In designing the test statistic, special care would need to be taken in making sure to take into account the look elsewhere effect (i.e., the fact that a particular observed anomaly may be very unlikely, whereas the probability of observing some anomaly may be much larger). From a Bayesian statistics point of view, it is the Bayesian evidence that can tell us how probable the extended models are, compared to the baseline power law primordial power spectrum.
For the models and the choice of prior probabilities considered here, the Bayesian evidence in fact favours, albeit weakly, the simple power law spectrum over the more complex models. The reason is that the Bayesian evidence punishes a lack of predictivity in these models. Most of the parameter space volume is not compatible with the data. A good match to observations is obtained within only a small subregion. Nonetheless, the observed features remain interesting since if they are real, they will also leave traces in other observabless, most notably, in the Emode polarization spectrum, where the signatures of features in the primordial spectrum are actually less washed out than in the temperature spectrum (Mortonson et al. 2009). The forthcoming Planck polarization data will prove very useful in this regard. Additionally, since strong deviations from power law behaviour typically indicate nonlinear physics, these models generically also predict a nonGaussian signal potentially observable in the bispectrum (Planck Collaboration XXIV 2014). However, the best fit wiggles and step inflation models have oscillations with a frequency too high to be accessible to bispectrum analysis at present.
9. Combined analysis with Planck f_{NL} constraints for single field inflation
In the previous sections we have analysed inflationary models with a canonical kinetic term. This led to the tensortoscalar consistency condition requiring n_{t} = − r/ 8. It is interesting to consider more general classes of inflationary models characterized by a nonstandard kinetic term (Garriga & Mukhanov 1999) or more general higherderivative operators (Kobayashi et al. 2010). An interesting subclass of these models are those in which the Lagrangian is a general function of the scalar inflaton field and its first derivative: L= P(φ,X), where X = − g^{μν}∂_{μ}φ∂_{ν}φ/ 2. A more general extension is provided by the socalled effective field theory of inflation (Cheung et al. 2008), which has a richer phenomenology.
We restrict our analysis to the first class of models (Garriga & Mukhanov 1999; Chen et al. 2007), which includes kinflation models (ArmendárizPicón et al. 1999; Garriga & Mukhanov 1999), and DiracBornInfield (DBI) models introduced in the context of brane inflation (Silverstein & Tong 2004; Alishahiha et al. 2004). In this class of models inflation can take place with a steep potential or can be driven by the kinetic term. One of the main features of inflationary models with a nonstandard kinetic term is that the inflaton fluctuations can propagate at a sound speed c_{s}< 1. As shown in previous analyses (e.g., Peiris et al. 2007; Powell et al. 2009; Lorenz et al. 2008; Agarwal & Bean 2009) there are strong degeneracies between the parameters determining the observable power spectra. Constraints on primordial nonGaussianity can help breaking this degeneracy, and we show how Planck’s combined measurement of the power spectrum and the nonlinearity parameter f_{NL} (Planck Collaboration XXIV 2014) improves constraints on this class of models.
In models with a nonstandard kinetic term the sound speed of the inflaton is given by (Garriga & Mukhanov 1999), so that in the canonical models, where P(φ,X) = V(φ) − X, one finds c_{s} = 1, while in general a nontrivial c_{s}< 1 corresponds to deviations from this standard case. Therefore, in these models, new parameters, such as the sound speed and its running, appear in the expressions for the inflationary observables. For the running of the sound speed it is useful to define an additional slowroll parameter (55)For values of the slowroll parameters much less than unity, the leading order scalar power spectrum is modified (Garriga & Mukhanov 1999) to (56)which is evaluated at kc_{s} = aH. The scalar spectral index gets an additional contribution from the running of the sound speed, (57)The gravitational sector remains unaltered by the nontrivial inflaton sound speed, retaining the same form as for the standard slowroll models. Therefore the usual consistency relation is modified to r ≈ − 8n_{t}c_{s} with n_{t} = − 2ϵ_{1} as usual (Garriga & Mukhanov 1999). The more accurate relation employed in this analysis is (58)This accounts for the difference in freezeout between the scalar and tensor perturbations (Peiris et al. 2007; Powell et al. 2009; Lorenz et al. 2008; Agarwal & Bean 2009) taking place at kc_{s} = aH for the scalar fluctuations, and at k = aH for the tensor modes.
Limiting ourselves to the predictions at lowest order in the slowroll parameters, there are clearly degeneracies between the parameters (A_{s},c_{s},ϵ_{1},ϵ_{2},s), which make the constraints on the inflationary power spectra observables less stringent in terms of these microscopic parameters. However, for models where the inflaton field has a nonstandard kinetic term with c_{s} ≪ 1, a high level of primordial nonGaussianity of the scalar perturbations is generated (see, e.g., Chen et al. 2007). In these models primordial nonGaussianity is produced by the higherderivative interaction terms that arise when expanding the kinetic part of the Lagrangian, P(φ,X). The amplitude of the nonGaussianity, defined by the nonlinearity parameter f_{NL}, receives two dominant contributions, arising from the inflaton interaction terms and . Each of them produces nonGaussianity shapes similar to the socalled equilateral type (Babich et al. 2004), i.e., a signal that peaks for equilateral triangles k_{1} = k_{2} = k_{3}. However, they are sufficiently distinct that the total signal can be very different from the equilateral one (Senatore et al. 2010). The nonlinearity parameter of the second interaction term is , while the other is determined by a second independent amplitude (Chen et al. 2007; Senatore et al. 2010). Constraints on the primordial nonGaussianity, presented in the companion paper Planck Collaboration XXIV (2014), thus allow us to construct a lower limit for the sound speed c_{s}. This helps reduce degeneracies in the parameter space of inflationary models with nonstandard kinetic terms. In particular, without the limits on the sound speed coming from the constraints on primordial nonGaussianity, it is not possible to derive an upper limit on the parameter ϵ_{1}, because the relation between the tensortoscalar ratio and ϵ_{1} also involves the sound speed (see, e.g., Eq. (58)).
In this paper, we consider three cases. One is a general analysis as described above, where we focus on the simplest case of a constant speed of sound with s = 0. From the Planck limits on primordial nonGaussianity in general single field models of inflation (Planck Collaboration XXIV 2014), the most conservative constraint on the sound speed is (59)In this large parameter space, we assume a uniform prior 0.02 ≤ c_{s} ≤ 1 in Eq. (58) within the HFF formalism described in the Appendix. We show the joint constraints on ϵ_{1} and ϵ_{2} in Fig. 20. By including the 95% CL constraint on c_{s} from Eq. (59), Planck+WP constrain ϵ_{1}< 0.053. Such constraints can be compared with the restricted case of c_{s} = 1, also shown in Fig. 20, with ϵ_{1}< 0.008 at 95% CL.
Fig. 20 Marginalized joint 68% and 95% CL for (ϵ_{1},ϵ_{2}) for Planck+WP data comparing the canonical Lagrangian case with c_{s} = 1 to the case of varying c_{s} with a uniform prior 0.02 <c_{s}< 1 derived from the Planck nonGaussianity measurements. 

Open with DEXTER 
The other two cases analysed correspond to specific models where the inflaton has a nonstandard kinetic term. The degeneracy between the different slowroll parameters is broken because these models specifically predict that s = 0, or s ∝ ϵ_{2}. As an example, we first consider the case where the action takes the DiracBornInfield (DBI) form (60)Here V(φ) is the potential and f(φ) is the warp factor determined by the geometry of the extra dimensions. For DBI models a stronger bound on c_{s} is derived (Planck Collaboration XXIV 2014): c_{s}> 0.07 at 95% CL. With the uniform prior 0.07 <c_{s}< 1 and s = 0, Planck + WP constrain ϵ_{1}< 0.042 at 95% CL.
An important case is f(φ) ≈ λ/φ^{4} (for details, see Silverstein & Tong (2004), Alishahiha et al. (2004), Chen et al. (2007), and references therein). There are two possibilities. First, in ultraviolet (UV) DBI models, the inflaton field moves under a quadratic potential V(φ) ≈ m^{2}φ^{2}/ 2 from the UV side of the warped space to the infrared side, with . It is known that this case is already at odds with observations if theoretical internal consistency of the model and constraints on power spectra and primordial nonGaussianity are taken into account (Baumann & McAllister 2007; Lidsey & Huston 2007; Bean et al. 2007, 2008; Peiris et al. 2007). It is therefore interesting to look at the other case, namely infrared DBI models (Chen 2005b,a) where the inflaton field moves from the IR to the UV side, and the inflaton potential is (61)with a wide range of values allowed for β in principle, 0.1 <β< 10^{9} (Bean et al. 2008). Here we focus on a minimal version of the IR DBI models where string effects are neglected, so that the usual field theory computation of the primordial curvature perturbation holds. For IR DBI models accounting for such effects and a more involved treatment of the dynamics, see Chen (2005a,c), and Bean et al. (2008). In this minimal IR DBI model, one finds (Chen 2005c; Chen et al. 2007) c_{s} ≈ (βN_{∗}/ 3)^{1}, n_{s} − 1 = − 4 /N_{∗}, and dn_{s}/d (in this model one can verify that s ≈ 1 /N_{∗} ≈ ϵ_{2}/ 3). Here primordial nonGaussianity of the equilateral type is generated with an amplitude .
If we consider 60 ≤ N_{∗} ≤ 90, then the predicted spectral index lies within the range 0.93 ≤ n_{s} ≤ 0.96, which is consistent with the Planck measurement of the spectral index at the 3σ level, for N_{∗} ≥ 60. The constraints on nonGaussianity give at 68% CL (Planck Collaboration XXIV 2014). Combining these constraints with the power spectrum constraints, marginalizing over 60 ≤ N_{∗} ≤ 90, we obtain (62)This strongly limits the allowed parameter space of these models.
As a final example, we consider a class of powerlaw kinflation models characterized by the Lagrangian (ArmendárizPicón et al. 1999) (63)In this case, for small values of γ one finds: , , n_{s} − 1 = − 3γ. The sound speed is a constant (s = 0), with constant γ. The primordial nonGaussianity in this model has an amplitude . Therefore, all the inflationary observables depend essentially on a single parameter γ. Imposing a prior of 0 <γ< 2 / 3 from the nonGaussianity constraint at 68% CL (Planck Collaboration XXIV 2014), we obtain γ ≥ 0.05 at 95% CL. At the same time, our measurement of the spectral index constrains 0.01 ≤ γ ≤ 0.02 at 95% CL. This class of kinflation models is therefore excluded by the combined constraints on primordial nonGaussianity and the power spectrum.
10. Isocurvature modes
10.1. Theoretical background
In this section we explore the constraints imposed by Planck on scenarios where the primordial cosmological perturbations were not entirely adiabatic. These scenarios also include isocurvature modes, possibly correlated among themselves as well as with the adiabatic mode. The adiabatic mode is characterized by the property that at very early times the universe obeyed a common, spatially uniform equation of state and all components initially shared a common velocity field. For the adiabatic mode the density perturbations in the various components (i.e., baryons, CDM, photons, and neutrinos) are locked together. Here baryons include their accompanying leptons, assumed tightly coupled to maintain charge neutrality.
Isocurvature modes arise from spatial variations in the equation of state or from relative velocities between the components. To analyse how the CMB perturbations were imprinted, it is most convenient to define isocurvature modes at a sufficiently late time, such that the relevant components, according to our present best understanding, consisted of baryons, photons, CDM, and neutrinos. Under this hypothesis, in addition to the adiabatic mode there are four possible nondecaying isocurvature modes: the baryon, CDM, and neutrino density isocurvature modes, and the neutrino velocity isocurvature mode (see, e.g., Bucher et al. 2000, for a discussion and further references).
The impact of isocurvature modes on the CMB was first studied in detail by Peebles & Yu (1970) and Efstathiou & Bond (1986, 1987), who contemplated the possibility that isocurvature perturbations rather than adiabatic perturbations were the sole source of cosmological perturbations. Linde (1985), Polarski & Starobinsky (1994), Linde & Mukhanov (1997), and GarcíaBellido & Wands (1996) pointed out various scenarios in which isocurvature perturbations could be generated within the context of inflation. Bucher et al. (2000) carried out a systematic study of isocurvature modes from a phenomenological perspective, pointing out the relevance of two additional modes: the neutrino density and velocity modes. Lyth & Wands (2002), Moroi & Takahashi (2001), and Bartolo & Liddle (2002) studied an interesting socalled curvaton scenario, in which adiabatic fluctuations from inflation contribute negligibly, but quantum fluctuations in a transverse direction modulate the density of decaying particles, leading to isocurvature perturbations correlated with the adiabatic mode.
Several authors have studied the constraints on isocurvature modes imposed by previous microwave background experiments, including Stompor et al. (1996), Langlois & Riazuelo (2000), Amendola et al. (2002), Peiris et al. (2003), Valiviita & Muhonen (2003), Bucher et al. (2004), Moodley et al. (2004), Beltran et al. (2004), KurkiSuonio et al. (2005), Dunkley et al. (2005), Bean et al. (2006), Trotta (2007), Keskitalo et al. (2007), and Komatsu et al. (2009). A more complete set of references may be found in Valiviita et al. (2012).
Before proceeding we must define precisely how to characterize these isocurvature modes on superHubble scales during the epoch after entropy generation, during which we assume that the stressenergy content of the universe can be modelled as a multicomponent fluid composed of baryons, CDM particles, photons, and neutrinos. If we assume that the evolution of the universe during this epoch was adiabatic (used here in the sense of thermodynamically reversible), then the entropy per unit comoving volume is conserved and serves as a useful reference with respect to which the abundances of the other components can be expressed.
The baryon isocurvature mode may be expressed in terms of fractional fluctuations in the baryontoentropy ratio, which is conserved on superHubble scales during this epoch. The CDM and neutrino density isocurvature (NDI) modes may be defined analogously. The neutrino velocity isocurvature (NVI) mode refers to fluctuations in the neutrino velocity relative to the average bulk velocity of the cosmic fluid. For the CMB, the baryon and CDM isocurvature modes yield almost identical angular spectra because the deficit of one is balanced by an excess of the other, so we do not consider them separately here. In this way the primordial isocurvature modes may be defined as dimensionless stochastic variables ℐ_{CDI}, ℐ_{NDI}, ℐ_{NVI}, like the variable ℛ describing the adiabatic mode^{21}. In this basis, the CDI mode can be seen as an effective isocurvature mode, encoding both CDM and baryon isocurvature fluctuations through (Gordon & Lewis 2003).
Within this framework, Gaussian fluctuations for the most general cosmological perturbation are described by a 4 × 4 positive definite matrixvalued power spectrum of the form (64)Following the conventions used in CAMB (Lewis & Bridle 2002; Lewis 2011) and CLASS (Lesgourgues 2011; Blas et al. 2011), the primordial isocurvature modes are normalized as follows in the synchronous gauge: for the CDI mode, is the primordial power spectrum of the density contrast δρ_{CDM}/ρ_{CDM}; for the NDI mode it is that of δρ_{ν}/ρ_{ν}; and for the NVI mode, that of the neutrino velocity v_{ν} times 4/3^{22}.
If isocurvature modes are present, the most plausible mechanism for exciting them involves inflation with a multicomponent inflaton field. To have an interesting spectrum on the large scales probed by the CMB, isocurvature modes require longrange correlations. Inflation with a multicomponent inflaton provides a well motivated scenario for establishing such correlations. Inflation with a singlecomponent scalar field can excite only the adiabatic mode. In models of inflation with light (compared to the Hubble expansion rate) transverse directions, the scalar field along these transverse directions becomes disordered in a way described by an approximately scaleinvariant spectrum. If the inflaton has M light components, there are (M − 1) potential isocurvature modes during inflation. Whether or not the fluctuations along these transverse directions are subsequently transformed into the latetime isocurvature modes described above depends on the details of what happens after inflation, as described more formally below.
As explained for example in Langlois (1999), Gordon et al. (2001), Groot Nibbelink & van Tent (2000, 2002), and Byrnes & Wands (2006), for inflationary models where the inflaton follows a curved trajectory, correlations are generically established between the isocurvature and curvature degrees of freedom. To lowest order in the slowroll approximation, this leads to a situation where the adiabatic perturbation is the sum of several components each of differing spectral index.
The postinflationary evolution determines how the isocurvature fluctuations generated during inflation transmute into the three specific isocurvature modes studied here. Little is known about the details of what happens during the epoch of entropy generation, but to linear order we may express how the fields ℛ_{inf} (i.e., the curvature perturbation at the end of inflation) and the transverse components of the inflaton field σ_{1},...,σ_{M − 1} (i.e., the components orthogonal to the slowroll direction) transform into curvature perturbations and latetime isocurvature modes at the end of the epoch of entropy generation as the linear transformation (65)where a = BI, CDI, NDI, NVI, while A = 1,...,(M − 1) labels the transverse components of the N component inflaton field. Physically, the fluctuations along the transverse directions modulate particle production during the epoch of entropy generation.
The neutrino density isocurvature can be excited in much the same way as the CDM and baryon isocurvature mode because at least within the standard electroweak model, in which there are no leptonic flavour changing processes, L_{e},L_{μ}, and L_{τ} are separately conserved. Known nonperturbative processes such as the sphaleron can trade lepton and baryon asymmetries with each other and alter flavour asymmetries, but they cannot erase such asymmetries altogether. Plausible models generating the neutrino density mode are therefore possible (Bucher et al. 2000; Gordon & Malik 2004), but as for the neutrino velocity isocurvature mode, to date no plausible generation mechanism has been put forth.
This extension of the adiabatic ΛCDM model to nonadiabatic initial conditions represents an important test of inflation. Singlefield inflation can produce only adiabatic perturbations, since exciting isocurvature perturbations requires additional degrees of freedom during inflation. Therefore a detection of primordial isocurvature perturbations would point to more complicated models of inflation.
10.2. Adiabatic with one isocurvature mode and free spectral indices
In this paper we investigate three of the four possible isocurvature modes of the ΛCDM scenario, since the baryon and CDM isocurvature perturbations are indistinguishable in the CMB angular power spectra. The CDM, neutrino density, and neutrino velocity isocurvature perturbations lead to different power spectra for CMB anisotropies, as shown in Fig. 21. We limit ourselves to studying one isocurvature mode at a time, in the presence of a curvature perturbation. More general combinations with two or three isocurvature modes may be contemplated, but without the Planck high frequency polarization likelihood, it is difficult to constrain such scenarios, so we postpone a discussion of this case to the next release.
Theoretically, one expects the power spectra of the isocurvature modes and their correlations to exhibit near but not necessarily exact scale invariance. As a general test of adiabaticity, it is nevertheless interesting to compare a more general model to the Planck data, assuming that the adiabatic, isocurvature, and crosscorrelation spectra obey power laws with free spectral indices. Blue values of the spectral indices are particularly interesting from the point of view of testing adiabaticity, because the acoustic peaks arising from two of the isocurvature modes are out of phase with the adiabatic peaks by roughly π/ 2 near the first acoustic peak. This is not true for the neutrino velocity isocurvature mode however.
Fig. 21 C_{TT} anisotropy shapes for the three isocurvature modes. Top: the shapes of the CDM isocurvature mode, neutrino density isocurvature mode, and neutrino velocity isocurvature mode are shown together with the adiabatic mode. The modes have the same amplitude parameters ( for the adiabatic mode and for each isocurvature mode). Bottom: the narrower multipole range illustrates the relative phases of the acoustic oscillations for these modes. 

Open with DEXTER 
Isocurvature mode constraints from Planck+WP data.
In the literature, models with one isocurvature as well as the adiabatic mode (possibly correlated) are often parameterized by specifying the 2 × 2 correlation matrix at a certain pivot scale k_{0} wth components along with their respective spectral indices n_{ℛℛ},n_{ℛℐ},n_{ℐℐ} (e.g., Amendola et al. 2002; Beltran et al. 2004). We do not follow this approach because in the absence of a statistically significant detection, the posterior distributions for the spectral indices are difficult to interpret and sensitive to how the prior is chosen. We instead adopt a parameterization where is specified at two scales k = k_{1} and k = k_{2} and interpolated geometrically according to^{23}(66)where a,b = ℐ,ℛ and ℐ = ℐ_{CDI}, ℐ_{NDI}, or ℐ_{NVI}. We set k_{1} = 2 × 10^{3}Mpc^{1} and k_{2} = 0.1Mpc^{1}, so that [k_{1},k_{2}] spans most of the range in k constrained by Planck data. A uniform prior for the components is assumed, where autocorrelation amplitudes are positive, although the crosscorrelation amplitudes may take both signs subject to the constraints (67)to ensure positive definiteness. For the logarithm of the offdiagonal elements in Eq. (66) to be real, we must have . This Ansatz can be trivially modified to admit the case by inserting appropriate minus signs, but this parameterization does not allow the case where the sign of the correlation changes. In practice we deal with this by assuming a uniform prior not on , but on its absolute value, and then we impose . The constraints in Eq. (67) ensure that is positive definite within the interval [k_{1},k_{2}] , but generically positive definiteness is violated sufficiently far outside this interval, either for very small or very large k. Where this happens we reduce the magnitude of so that there is either total correlation or anticorrelation. The kinks thus introduced lie outside the range [k_{1},k_{2}] . Within this range, the spectral indices n_{ℛℛ},n_{ℛℐ},n_{ℐℐ} are scaleindependent. Finally, our sign conventions are such that positive values for correspond to a positive contribution of the crosscorrelation term to the SachsWolfe component of the total temperature spectrum.
When the constraining power of the data is weak, a crucial question is to what extent the posterior distribution results from the data rather than from the prior distribution. The parameterization above is not the only one that could have been adopted, and other possible priors are typically related by a nonconstant Jacobian. For each model, we indicate the loglikelihood for the best fit model, in order to allow model comparison.
Fig. 22 Two dimensional distributions for power in isocurvature modes using Planck+WP data. 

Open with DEXTER 
The Planck+WP results for the three isocurvature modes using this twoscale parameterization are shown in Fig. 22 and included in the summary Table 12. The power spectra and are normalized according to the primordial values of the fields ℛ(x) and ℐ(x) defined above. It is interesting to consider how much isocurvature power is allowed expressed as a fraction of the power in three bands spanning the CMB temperature spectrum observed by Planck. To this end, we define the following derived quantities where (71)The 95% confidence limits from the onedimensional posterior distributions for these fractional contributions in the full range (ℓ_{min},ℓ_{max}) = (2,2500) are shown in Table 12. The range of allowed values for α_{ℛℛ}(2,2500) is a measure of the adiabaticity of fluctuations in the CMB. The posterior distributions of the fractions α_{ℐℐ},α_{ℛℐ} in three multipole ranges are shown in Fig. 23. We also report the primordial isocurvature fraction, defined as (72)at three values of k. Table 12 also shows the effective χ^{2} = − 2lnℒ_{max} for all models, compared to the minimal sixparameter ΛCDM model. In Fig. 24 we show the ratio of temperature spectra for the best fit mixed model to the adiabatic model.
Fig. 23 Fractional contribution of isocurvature modes to the power spectrum. We show the distributions α_{ℐℐ}(2,20), α_{ℛℐ}(2,20), α_{ℐℐ}(21,200), α_{ℛℐ}(21,200), α_{ℐℐ}(201,2500), and α_{ℛℐ}(201,2500), defined in Eq. (69), for the CDI, NDI, and NVI modes as constrained by the Planck+WP data. 

Open with DEXTER 
Fig. 24 Temperature spectrum of best fit models with a mixture of adiabatic and isocurvature modes. Top: spectrum of the best fit mixed models relative to that of the pure adiabatic model. Bottom: zoom of the SachsWolfe plateau of the best fit temperature spectrum for each of the three cases plus the pure adiabatic model, shown together with Planck lowℓ data points. 

Open with DEXTER 
The results for α_{ℛℛ}(2,2500) show that the nonadiabatic contribution to the temperature variance can be as large as 7% (9%, 5%) in the CDI (NDI, NVI) model (at 95% CL). These results are driven by the fact that on large scales, for ℓ ≤ 40, the Planck data points on average have a slightly smaller amplitude than the best fitting ΛCDM model. Hence the data prefer a significant amount of anticorrelated isocurvature modes, leading to a reduction of amplitude of the SachsWolfe plateau and to a decrease of the effective χ^{2} by up to 4.6^{24}. This situation explains the rather loose bounds on the derived parameter α_{ℐℐ}(2,20), as shown in Fig. 23.
A comparison of and shows that best fitting models have an isocurvature spectral index n_{ℐℐ} close to 1.7 for CDI, 1.1 for NDI, and 1.0 for NVI modes.
For CDI and NDI, the amplitude of acoustic peaks quickly decreases with increasing ℓ, so that the constraints are entirely driven by small ℓs. Since the same value of the primordial amplitude leads to different plateau amplitudes for the two isocurvature models (see Fig. 21), the bounds on and are consistently stronger for CDI than for NDI. For NVI, the acoustic peak amplitude is larger than the plateau amplitude. In NVI models, the data cannot allow for a too large amplitude of correlated isocurvature modes at small ℓ, because the total spectrum would be distorted at larger ℓ. This possibility is strongly disfavoured by the data, which are consistent with the peak location predicted by a pure adiabatic model. Hence in the NVI case we obtain slightly stronger bounds and a smaller reduction of the effective χ^{2}.
The fact that the data prefer models with a significant contribution from CDI or NDI modes should be interpreted with care. The detection of a shift in the phase of acoustic oscillations would bring unambigous evidence in favour of isocurvature modes. With Planck data, we are not in this situation. The evidence is driven by a small deficit of amplitude in the SachsWolfe plateau, which could have several different possible explanations (such as a deficit in the largescale primordial power spectrum, as already seen in the previous sections). However, multifield inflationary scenarios can produce the mixture of curvature and isocurvature fluctuations which we have found to provide a good fit to the Planck data.
10.3. Special cases
The sixparameter models of the previous subsection including one isocurvature mode and the adiabatic mode make no assumptions about the spectral indices of each mode or the degree of correlation between the isocurvature mode and the adiabatic mode. This leads to a large number of additional degrees of freedom. There are both theoretical and phenomenological motivations for choosing special values for some of the parameters, leading to special cases with just one more degree of freedom with respect to the adiabatic case. The results are reported in Table 12, for uncorrelated perturbations with n_{ℐℐ} = 1, and fully correlated or anticorrelated perturbations with n_{ℐℐ} = n_{ℛℛ}. In the general case, anticorrelated isocurvature perturbations slightly improve the fit to the Planck data. We consider below the implications of our results for two important cases: the axion and curvaton scenarios.
10.3.1. Constraints on axion isocurvature
The axion field was proposed to solve the strong CP problem and constitutes a wellmotivated dark matter candidate. (See for example Preskill et al. 1983; Turner 1990; Peccei 2008; Sikivie 2008; Raffelt 2008; Kim & Carosi 2010.) The axion is the pseudoGoldstone boson of the broken PecceiQuinn (PQ) symmetry. Under certain assumptions, the axion field may induce significant isocurvature perturbations (Turner et al. 1983; Axenides et al. 1983; Steinhardt & Turner 1983; Linde 1984, 1985, 1991; Seckel & Turner 1985; Kofman 1986; Lyth 1990, 1992; Linde & Lyth 1990; Turner & Wilczek 1991). If inflation takes place after PQ symmetry breaking, the quantum fluctuations of the inflaton are responsible for primordial curvature perturbations, while those of the axion field generate primordial entropy perturbations. After the QCD transition, when one of the vacua becomes preferred giving the axion field a mass, the axions behave as cold dark matter. This way of producing axionic dark matter is called the misalignment angle mechanism. In such a scenario, the CMB anisotropy may include significant power from CDM isocurvature fluctuations. In that case, the fraction of CDM isocurvature modes is related to the energy scale of inflation, H_{inf}, through (Lyth 1990; Beltran et al. 2007; Bae et al. 2008; Hamann et al. 2009) (73)where Ω_{a}h^{2} is the relic axion density, R_{a} the fraction of CDM consisting of axions, and f_{a} the PQ symmetry breaking scale. In this model, CDM isocurvature perturbations should be totally uncorrelated with adiabatic perturbations and have a spectral index n_{ℐℐ} very close to one since in the firstorder slowroll approximation the index reads (1 − 2ϵ_{V}). Since the sensitivity of the data to n_{ℐℐ} is very limited (Beltran et al. 2007), we assume for simplicity that n_{ℐℐ} = 1.
Within the general parametrization presented in Eq. (66), we can select the axion case by imposing , as well as the condition (74)corresponding to n_{ℐℐ} = 1. We therefore have three independent parameters, and , and we sample these parameters with uniform prior distributions. The fraction β_{iso}(k_{∗}) with k_{∗} = 0.05 Mpc^{1} is then a derived parameter. Since the data constrain β_{iso} ≪ 1, the relation between β_{iso} and is nearly linear, so the primordial isocurvature fraction is sampled with a close to uniform prior.
Constraints on this model are shown in Table 12. We find (75)at the scale k_{mid} = 0.05 Mpc^{1}, with a best fit value of zero. Hence there is no evidence for axion generated isocurvature perturbations. This limit significantly improves the previous CMB bounds. At the scale k = 0.002 Mpc^{1}, our result reads β_{iso}< 0.036, to be compared to β_{iso}< 0.15 for WMAP 9year alone, or β_{iso}< 0.061 for WMAP+ACT+SPT at 95% CL (Hinshaw et al. 2013). This bound can be used to exclude regions of parameter space composed of f_{a}, R_{a}, and the energy scale of inflation, but cannot be used to obtain a modelindependent bound on f_{a}. However, if we assume (i) that the PQ symmetry is broken during inflation; (ii) that it is not restored by the quantum fluctuations of the inflaton (which imposes H_{inf}/ (2π) <f_{a}), nor by thermal fluctuations in case of a very efficient reheating stage; and (iii) that all the CDM consists of axions produced by the misalignment angle, then we can derive an upper bound on the energy scale of inflation as (76)
10.3.2. Constraints on the curvaton scenario
In the simplest onefield inflationary models curvature perturbations arise from quantum fluctuations in the inflaton field, but this is not the only way to generate curvature perturbations. Isocurvature perturbations may seed curvature perturbations outside the Hubble radius (Polarski & Starobinsky 1994; Langlois 1999; Gordon et al. 2001), so it is possible that a significant component of the observed adiabatic mode could be strongly correlated with an isocurvature mode. This happens for instance in the curvaton scenario (Mollerach 1990; Enqvist & Sloth 2002; Moroi & Takahashi 2001; Lyth & Wands 2002; Lyth et al. 2003; Gordon & Lewis 2003). The curvaton is an extra light scalar field acquiring a spectrum of fluctuations on cosmological scales during inflation. Depending on its density evolution and decay history, this field could be responsible for part of the observed adiabatic perturbations, or all of them, or for a mixture of correlated adiabatic and isocurvature perturbations.
We focus here on the simplest viable version of this scenario in which the curvaton decays into CDM particles while contributing a nonnegligible fraction (77)to the total energy density of the universe. If the curvaton dominates at decay time (r_{D} = 1), its primordial fluctuations seed curvature perturbations equivalent to a pure adiabatic mode. If r_{D}< 1, curvaton fluctuations are only partially converted into adiabatic perturbations, while CDM particles carry CDI perturbations, which are fully correlated with the adiabatic perturbations since they share a common origin. We recall that with our conventions, “fully correlated” means that the crosscorrelation term contributes constructively to the SachsWolfe component of the total temperature spectrum. Some authors define the correlation with the opposite sign and call this case “fully anticorrelated” (e.g., Komatsu et al. 2011; Hinshaw et al. 2013). In this model, the CDI fraction is related to r_{D} by (Gordon & Lewis 2003) (78)In our notation this is equivalent to (79)Within the general parametrization presented in Eq. (66), we can satisfy this case by imposing (80)together with the condition (81)corresponding to n_{ℐℐ} = n_{ℛℛ}. As in the axion case, this results in three independent parameters and which we sample with uniform priors. The constraints for this model are shown in Table 12. The best fit model is still the pure adiabatic case, and the upper bound (82)is scale independent, since the adiabatic and isocurvature tilts are assumed to be equal. This is a significant improvement over the WMAP 9year bounds, β_{iso}< 0.012 for WMAP alone, or β_{iso}< 0.0076 for WMAP+ACT+SPT at 95% CL (Hinshaw et al. 2013). We conclude that in this scenario, the curvaton should decay when it dominates the energy density of the universe, with r_{D}> 0.983.
The nonlinearity parameter in the curvaton model studied here is (Bartolo et al. 2004c,b) (83)assuming a quadratic potential for the curvaton field (Sasaki et al. 2006). In the pure adiabatic case (r_{D} = 1) this leads to . The constraint 0.98 <r_{D}< 1 then corresponds to . Taking into account the Planck result (Planck Collaboration XXIV 2014), we conclude that the Planck data are consistent with the scenario where the curvaton decays into CDM when it dominates the energy density of the universe, and its fluctuations are almost entirely converted into adiabatic ones.
11. Conclusions
This paper establishes the status of cosmic inflation in the context of the first release of the Planck cosmological results, which includes the temperature data from the first 2.6 sky surveys. CMB polarization as measured by Planck will be the subject of a future release. We find that standard slowroll singlefield inflation is compatible with the Planck data. This result is confirmed by other papers of this series. Planck in combination with WMAP 9year large angular scale polarization (WP) yields Ω_{K} = − 0.006 ± 0.018 at 95% CL by combining temperature and lensing information (Planck Collaboration XVI 2014; Planck Collaboration XVII 2014). The bispectral nonGaussianity parameter f_{NL} measured by Planck is consistent with zero (Planck Collaboration XXIV 2014). These results are compatible with zero spatial curvature and a small value of f_{NL}, as predicted in the simplest slowroll inflationary models.
A key Planck result is the measurement of the scalar perturbation spectral index. Planck+WP data give n_{s} = 0.9603 ± 0.0073 (and n_{s} = 0.9629 ± 0.0057 when combined with BAO). This result disfavours the HarrisonZeldovich (HZ) n_{s} = 1 model at more than 5σ. Even in extended cosmological models, the HZ spectrum cannot be reconciled with the data. Allowing a general reionization scenario yields with respect to ΛCDM for Planck+WP data. When the primordial helium abundance or the effective number of neutrino species are allowed to vary, the best fit of the HZ model to a combination of Planck+WP and BAO data is still worse by and 8.0, respectively.
We find no evidence for Planck data preferring a generalization of a simple power law spectrum to include a running of the spectral index (dn_{s}/ dlnk = − 0.0134 ± 0.0090) or a running of the running ( with Planck+WP). In a model admitting tensor fluctuations, the 95% CL bound on the tensortoscalar ratio is r_{0.002}< 0.12 (< 0.11) using Planck+WP (plus highℓ CMB data). This bound on r implies an upper limit for the inflation energy scale of 1.9 × 10^{16}GeV, or equivalently, for the Hubble parameter H_{∗}< 3.7 × 10^{5}M_{pl}, at 95% CL.
The degeneracy between n_{s} and r, which plagued previous CMB measurements, is now removed by the Planck precision in the determination of the highest acoustic peaks. Inflaton potentials with a concave shape are favoured and occupy most of the 95% confidence region allowed by Planck+WP in the n_{s}  r plane. Models with an exponential potential, a monomial potential with a power larger than two, or hybrid models driven by a quadratic term are disfavoured at more than 95% confidence. The quadratic largefield model, in the past often cited as the simplest inflationary model, now lies at the edge of the 95% CL contours allowed by Planck+WP+highℓ CMB data.
A Bayesian parameter estimation and model comparison analysis of a representative sample of singlefield slowroll models shows that Planck is able to discriminate between these models with results that are robust even when a broad set of entropy generation scenarios are allowed. In addition to confirming the exclusion of the φ^{4} potential, the Bayesian evidence computed from the Planck data provides significant odds (logarithms of the Bayes factor of about −5 or lower relative to ΛCDM) against largefield models compatible with previous cosmological data, such as the φ^{2} potential, and twoparameter potentials such as natural inflation and the hilltop potential. As presented in Sect. 5, Planck establishes strong constraints on the parameter values of specific inflationary scenarios. For example, the scale parameter of the natural inflation potential is constrained to be log (f/M_{pl}) ≳ 1.1 (95% CL), improving upon the WMAP 7year limit on f by a factor of two. The Planck data limit the possibilities for the unexplored physics between the end of inflation and the beginning of the radiation dominated era. Datadriven constraints are obtained on w_{int}, the effective equation of state in the postinflationary era. Particularly for the disfavoured models listed above, their parameters are pushed to unnatural values (w_{int} ≳ 1 / 3) in order to become more compatible with the data.
Using an essentially exact numerical calculation of the predicted primordial spectrum, we reconstruct the observable window of the inflaton potential, expanding the potential as a Taylor series up to a fixed order. For an observable potential described by a polynomial of order three, the reconstruction agrees well with the slowroll predictions. If a quartic term is allowed, the result deviates from the slowroll prediction because the Planck data favour a slightly smaller amplitude for the SachsWolfe plateau relative to the ℓ > 40 part of the power spectrum than the best fitting minimal ΛCDM model with a power law primordial spectrum. A potential with a fourthorder polynomial can fit this feature, thus reducing the effective by approximately four.
A penalized likelihood reconstruction of the primordial power spectrum shows hints of structure at modest statistical significance. However, recent work after submission suggests that this feature can be explained by electromagnetic interference. Parameterized models producing superimposed oscillations (possibly motivated by deviations from the BunchDavies vacuum state, axion monodromy, or a sharp step in the inflaton potential) improve the by roughly 10, where three extra parameters have been added. However, a Bayesian model comparison analysis does not strongly favour the model with oscillations over the standard featureless power spectrum. With Planck polarization data, a more conclusive result on superimposed oscillations is expected.
We combine power spectrum constraints with those on the nonlinearity parameter f_{NL} (Planck Collaboration XXIV 2014) to constrain singlefield inflation with generalized Lagrangians, in which nonGaussianities are larger than those predicted by the simplest slowroll inflationary models. We show how the limits on the inflation sound speed derived in Planck Collaboration XXIV (2014) are crucial to constrain slowroll parameters for generalized Lagrangians. We also show how particular examples of DBI Inflation and kinflation can be constrained by this combination of Planck data.
We test the hypothesis that the primordial cosmological perturbations were exclusively adiabatic. We analyse all nonsingular (i.e., nondecaying) isocurvature modes arbitrarily correlated to the adiabatic mode, using a parameterization where the isocurvature contributions are specified at two scales. The oscillatory pattern in the Planck temperature spectrum is compatible with purely adiabatic perturbations, and therefore constrains any isocurvature contribution to be small at those multipoles. As a consequence, axion and curvaton scenarios, in which the CDM isocurvature mode is uncorrelated or fully correlated with the adiabatic mode, are not favoured by Planck. The upper bounds on the isocurvature fraction at k = 0.05 Mpc^{1} are 0.039 for the axion, and 0.0025 for the curvaton, at 95% CL. However general models with an arbitrarily correlated mixture of adiabatic and (CDM or neutrino) isocurvature modes have the freedom to lower the SachsWolfe plateau relative to the highℓ spectrum, and reduce the effective by more than four.
The simplest inflationary models have passed an exacting test with the Planck data. The full mission data including Planck’s polarization measurements will help answer further fundamental questions, helping to probe nonsmooth power spectra and the energy scale of inflation as well as extensions to more complex models.
Planck (http://www.esa.int/Planck) is a project of the European Space Agency (ESA) with instruments provided by two scientific consortia funded by ESA member states (in particular the lead countries France and Italy), with contributions from NASA (USA) and telescope reflectors provided by a collaboration between ESA and a scientific consortium led and funded by Denmark.
After the Planck March 2013 release, a minor error was found in the ordering of the beam transfer functions applied to the 217 × 217 crossspectra in the Planck highℓ likelihood. An extensive analysis of the corresponding revised Planck highℓ likelihood showed that this error has a negligible impact on cosmological parameters and is absorbed by small shifts in the foreground parameters. See Planck Collaboration XVI (2014) for more details.
Note that since the average is performed over the entire support of the prior probability density, the evidence depends strongly on the probability range for the adjustable parameters. Whereas in parametric inference, the exact extent of the prior ranges often becomes irrelevant as long as they are “wide enough” (i.e., containing the bulk of the highlikelihood region in parameter space), the value of the evidence will generally depend on precisely how wide the prior range was chosen.
After the submission of the first version of this paper, uncertainties arising from the minimization algorithm in the best fit cosmological parameters and the best fit likelihood were studied. The uncertainties found were and therefore do not alter our conclusions. The values for Δχ^{2} reported have not been updated.
See Starobinsky (2005) for the inflationary model producing an exactly scaleinvariant power spectrum with r ≠ 0 beyond the slowroll approximation.
Available at: http://www.sciops.esa.int/index.php?project=planck&page=Planck_Legacy_Archive
We considered a further generalization, which also causes the joint constraints on n_{s} and r to change slightly. We allowed the amplitude of the lensing contribution to the temperature power spectrum A_{L} to vary as a free parameter. In this case we find the following Planck+WP+BAO constraints: n_{s} = 0.972 ± 0.006, r< 0.15 at 95% CL.
Selected nonstandard values for N_{eff} deserve further investigation. For the additional fractional contribution motivated by a Goldstone boson ΔN_{eff} = 0.39 (Weinberg 2013), we obtain n_{s} = 0.9726 ± 0.0057 and r_{0.002}< 0.14 at 95% CL for Planck+WP+BAO+highℓ. For an additional species of neutrinos, we obtain n_{s} = 0.9947 ± 0.0056 and r_{0.002}< 0.17 at 95% CL for Planck+WP+BAO+highℓ.
Made available ahead of public release to the Planck Collaboration by Farhan Feroz and Mike Hobson (Feroz et al. 2013).
In other words, of the neutrino perturbation dipole, F_{ν1} = 4θ_{v}/ (3k) in the notation of Ma & Bertschinger (1995).
For the three general models, the posterior distribution is actually multimodal. Here we are referring to models contributing to the main peak in the posterior, with the highest maximum likelihood. There is another peak with a smaller maximum likelihood, appearing in Fig. 23 as a small bump for positive values of the crosscorrelation amplitude. In this paper, we do not carry out a separate investigation for models contributing to this secondary peak.
Acknowledgments
The development of Planck has been supported by: ESA; CNES and CNRS/INSUIN2P3INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MICINN, JA and RES (Spain); Tekes, AoF and CSC (Finland); DLR and MPG (Germany); CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerland); RCN (Norway); SFI (Ireland); FCT/MCTES (Portugal); and PRACE (EU). A description of the Planck Collaboration and a list of its members, including the technical or scientific activities in which they have been involved, can be found at http://www.sciops.esa.int/index.php?project=planck&page=Planck_Collaboration. We gratefully acknowledge CINECA (http://www.cineca.it/) under the agreement LFI/CINECA and IN2P3 Computer Center (http://cc.in2p3.fr) for providing a significant amount of the computing resources and services needed for this work.
References
 Abbott, L., & Wise, M. B. 1984, Nucl. Phys. B, 244, 541 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Achúcarro, A., Gong, J.O., Hardeman, S., Palma, G. A., & Patil, S. P. 2011, JCAP, 1101, 030 [NASA ADS] [CrossRef] [Google Scholar]
 Acquaviva, V., Bartolo, N., Matarrese, S., & Riotto, A. 2003, Nucl. Phys. B, 667, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Adams, F. C., Bond, J. R., Freese, K., Frieman, J. A., & Olinto, A. V. 1993, Phys. Rev. D, 47, 426 [NASA ADS] [CrossRef] [Google Scholar]
 Adams, J., Cresswell, B., & Easther, R. 2001, Phys. Rev. D, 64, 123514 [NASA ADS] [CrossRef] [Google Scholar]
 Adshead, P., Easther, R., Pritchard, J., & Loeb, A. 2011, JCAP, 1102, 021 [NASA ADS] [CrossRef] [Google Scholar]
 Adshead, P., Dvorkin, C., Hu, W., & Lim, E. A. 2012, Phys. Rev. D, 85, 3531 [NASA ADS] [Google Scholar]
 Agarwal, N., & Bean, R. 2009, Phys. Rev. D, 79, 023503 [NASA ADS] [CrossRef] [Google Scholar]
 Akaike, H. 1974, IEEE Trans. Automat. Cont., 19, 716 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Albrecht, A., & Steinhardt, P. J. 1982, Phys. Rev. Lett., 48, 1220 [NASA ADS] [CrossRef] [Google Scholar]
 Alishahiha, M., Silverstein, E., & Tong, D. 2004, Phys. Rev. D, 70, 123505 [NASA ADS] [CrossRef] [Google Scholar]
 Allahverdi, R., Brandenberger, R., CyrRacine, F.Y., & Mazumdar, A. 2010, Ann. Rev. Nucl. Part. Sci., 60, 27 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Amendola, L., Gordon, C., Wands, D., & Sasaki, M. 2002, Phys. Rev. Lett., 88, 211302 [NASA ADS] [CrossRef] [Google Scholar]
 Anantua, R., Easther, R., & Giblin, J. T. 2009, Phys. Rev. Lett., 103, 111303 [NASA ADS] [CrossRef] [Google Scholar]
 Anderson, L., Aubourg, E., Bailey, S., et al. 2012, MNRAS, 427, 3435 [NASA ADS] [CrossRef] [Google Scholar]
 ArmendárizPicón, C., Damour, T., & Mukhanov, V. 1999, Phys. Lett. B, 458, 209 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Audren, B., Lesgourgues, J., Benabed, K., & Prunet, S. 2012, JCAP, 1302, 001 [NASA ADS] [CrossRef] [Google Scholar]
 Aver, E., Olive, K. A., & Skillman, E. D. 2012, JCAP, 1204, 004 [NASA ADS] [CrossRef] [Google Scholar]
 Axenides, M., Brandenberger, R. H., & Turner, M. S. 1983, Phys. Lett. B, 126, 178 [NASA ADS] [CrossRef] [Google Scholar]
 Babich, D., Creminelli, P., & Zaldarriaga, M. 2004, JCAP, 8, 9 [NASA ADS] [CrossRef] [Google Scholar]
 Bae, K. J., Huh, J.H., & Kim, J. E. 2008, JCAP, 0809, 005 [NASA ADS] [CrossRef] [Google Scholar]
 Banks, T., Kaplan, D. B., & Nelson, A. E. 1994, Phys. Rev. D, 49, 779 [NASA ADS] [CrossRef] [Google Scholar]
 Bardeen, J. M., Steinhardt, P. J., & Turner, M. S. 1983, Phys. Rev. D, 28, 679 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Barrow, J. D. 1990, Phys. Lett. B, 235, 40 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Barrow, J. D., & Liddle, A. R. 1993, Phys. Rev. D, 47, 5219 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Bartolo, N., & Liddle, A. R. 2002, Phys.Rev. D, 65, 121301 [NASA ADS] [CrossRef] [Google Scholar]
 Bartolo, N., Komatsu, E., Matarrese, S., & Riotto, A. 2004a, Phys. Rept., 402, 103 [NASA ADS] [CrossRef] [Google Scholar]
 Bartolo, N., Matarrese, S., & Riotto, A. 2004b, Phys. Rev. Lett., 93, 231301 [NASA ADS] [CrossRef] [Google Scholar]
 Bartolo, N., Matarrese, S., & Riotto, A. 2004c, Phys. Rev. D, 69, 043503 [NASA ADS] [CrossRef] [Google Scholar]
 Barvinsky, A. O., Kamenshchik, A. Y., & Starobinsky, A. A. 2008, JCAP, 0811, 021 [NASA ADS] [CrossRef] [Google Scholar]
 Baumann, D., & McAllister, L. 2007, Phys. Rev. D, 75, 123508 [NASA ADS] [CrossRef] [Google Scholar]
 Bean, R., Dunkley, J., & Pierpaoli, E. 2006, Phys. Rev. D, 74, 063503 [NASA ADS] [CrossRef] [Google Scholar]
 Bean, R., Shandera, S. E., Tye, S.H. H., & Xu, J. 2007, JCAP, 5, 4 [NASA ADS] [CrossRef] [Google Scholar]
 Bean, R., Chen, X., Peiris, H., & Xu, J. 2008, Phys. Rev. D, 77, 023527 [NASA ADS] [CrossRef] [Google Scholar]
 Beltran, M., GarcíaBellido, J., Lesgourgues, J., & Riazuelo, A. 2004, Phys. Rev. D, 70, 103530 [NASA ADS] [CrossRef] [Google Scholar]
 Beltran, M., GarcíaBellido, J., & Lesgourgues, J. 2007, Phys. Rev. D, 75, 103507 [NASA ADS] [CrossRef] [Google Scholar]
 Benetti, M., Pandolfi, S., Lattanzi, M., Martinelli, M., & Melchiorri, A. 2013, Phys. Rev. D, 87, 023519 [NASA ADS] [CrossRef] [Google Scholar]
 Bennett, C., Hill, R., Hinshaw, G., et al. 2011, ApJS, 192, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Bennett, C., Larson, D., Weiland, J., et al. 2013, ApJS, 208, 20 [NASA ADS] [CrossRef] [Google Scholar]
 Beutler, F., Blake, C., Colless, M., et al. 2011, MNRAS, 416, 3017 [NASA ADS] [CrossRef] [Google Scholar]
 Bezrukov, F., & Gorbunov, D. 2012, Phys. Lett. B, 713, 365 [NASA ADS] [CrossRef] [Google Scholar]
 Bezrukov, F., & Shaposhnikov, M. 2008, Phys. Lett. B, 659, 703 [NASA ADS] [CrossRef] [Google Scholar]
 Bezrukov, F., & Shaposhnikov, M. 2009, J. High Energy Phys., 07, 089 [NASA ADS] [CrossRef] [Google Scholar]
 Binétruy, P., & Gaillard, M. K. 1986, Phys. Rev. D, 34, 3069 [NASA ADS] [CrossRef] [Google Scholar]
 Blas, D., Lesgourgues, J., & Tram, T. 2011, JCAP, 1107, 034 [NASA ADS] [CrossRef] [Google Scholar]
 Boubekeur, L., & Lyth, D. 2005, JCAP, 0507, 010 [NASA ADS] [CrossRef] [Google Scholar]
 Bozza, V., Giovannini, M., & Veneziano, G. 2003, JCAP, 0305, 001 [NASA ADS] [CrossRef] [Google Scholar]
 Bridges, M., Feroz, F., Hobson, M. P., & Lasenby, A. N. 2009, MNRAS, 400, 1075 [NASA ADS] [CrossRef] [Google Scholar]
 Bridle, S., Lewis, A., Weller, J., & Efstathiou, G. 2003, MNRAS, 342, L72 [NASA ADS] [CrossRef] [Google Scholar]
 Brout, R., Englert, F., & Gunzig, E. 1978, Ann. Phys., 115, 78 [NASA ADS] [CrossRef] [Google Scholar]
 Bucher, M., & Cohn, J. 1997, Phys. Rev. D, 55, 7461 [NASA ADS] [CrossRef] [Google Scholar]
 Bucher, M., & Turok, N. 1995, Phys. Rev. D, 52, 5538 [NASA ADS] [CrossRef] [Google Scholar]
 Bucher, M., Goldhaber, A. S., & Turok, N. 1995, Phys. Rev. D, 52, 3314 [NASA ADS] [CrossRef] [Google Scholar]
 Bucher, M., Moodley, K., & Turok, N. 2000, Phys. Rev. D, 62, 083508 [NASA ADS] [CrossRef] [Google Scholar]
 Bucher, M., Dunkley, J., Ferreira, P., Moodley, K., & Skordis, C. 2004, Phys. Rev. Lett., 93, 081301 [NASA ADS] [CrossRef] [Google Scholar]
 Bunch, T., & Davies, P. 1978, Proc. Roy. Soc. London A, 360, 117 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Burgess, C., Easther, R., Mazumdar, A., Mota, D. F., & Multamaki, T. 2005, J. High Energy Phys., 0505, 067 [NASA ADS] [CrossRef] [Google Scholar]
 Byrnes, C. T., & Wands, D. 2006, Phys. Rev. D, 74, 043529 [NASA ADS] [CrossRef] [Google Scholar]
 Casadio, R., Finelli, F., Kamenshchik, A., Luzzi, M., & Venturi, G. 2006, JCAP, 0604, 011 [NASA ADS] [CrossRef] [Google Scholar]
 Chen, X. 2005a, J. High Energy Phys., 8, 45 [NASA ADS] [CrossRef] [Google Scholar]
 Chen, X. 2005b, Phys. Rev. D, 71, 063506 [NASA ADS] [CrossRef] [Google Scholar]
 Chen, X. 2005c, Phys. Rev. D, 72, 123518 [NASA ADS] [CrossRef] [Google Scholar]
 Chen, X. 2010, Adv. Astron., id.638979 [Google Scholar]
 Chen, X., Huang, M.X., Kachru, S., & Shiu, G. 2007, JCAP, 1, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Cheung, C., Creminelli, P., Fitzpatrick, A. L., Kaplan, J., & Senatore, L. 2008, J. High Energy Phys., 0803, 014 [NASA ADS] [CrossRef] [Google Scholar]
 Chung, D. J., Kolb, E. W., Riotto, A., & Tkachev, I. I. 2000, Phys. Rev. D, 62, 043508 [NASA ADS] [CrossRef] [Google Scholar]
 Chung, D. J., Everett, L. L., & Matchev, K. T. 2007, Phys. Rev. D, 76, 103530 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Coleman, S. R., & De Luccia, F. 1980, Phys. Rev. D, 21, 3305 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Conley, A., Guy, J., Sullivan, M., et al. 2011, ApJS, 192, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Contaldi, C. R., Peloso, M., Kofman, L., & Linde, A. D. 2003, JCAP, 0307, 002 [NASA ADS] [CrossRef] [Google Scholar]
 Cortês, M., & Liddle, A. R. 2009, Phys. Rev. D, 80, 083524 [NASA ADS] [CrossRef] [Google Scholar]
 Cortês, M., Liddle, A. R., & Mukherjee, P. 2007, Phys. Rev. D, 75, 083520 [NASA ADS] [CrossRef] [Google Scholar]
 Covi, L., Hamann, J., Melchiorri, A., Slosar, A., & Sorbera, I. 2006, Phys. Rev. D, 74, 083509 [NASA ADS] [CrossRef] [Google Scholar]
 Cox, R. T. 1946, Am. J. Phys., 14, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Danielsson, U. H. 2002, Phys. Rev. D, 66, 023511 [NASA ADS] [CrossRef] [Google Scholar]
 Das, S., Louis, T., Nolta, M. R., et al. 2014, JCAP, 04, 014 [NASA ADS] [CrossRef] [Google Scholar]
 de Carlos, B., Casas, J. A., Quevedo, F., & Roulet, E. 1993, Phys. Lett. B, 318, 447 [NASA ADS] [CrossRef] [Google Scholar]
 Dodelson, S., Kinney, W. H., & Kolb, E. W. 1997, Phys. Rev. D, 56, 3207 [NASA ADS] [CrossRef] [Google Scholar]
 Dunkley, J., Bucher, M., Ferreira, P., Moodley, K., & Skordis, C. 2005, Phys. Rev. Lett., 95, 261303 [NASA ADS] [CrossRef] [Google Scholar]
 Dunkley, J., Komatsu, E., Nolta, M. R., et al. 2009, ApJS, 180, 306 [NASA ADS] [CrossRef] [Google Scholar]
 Dunkley, J., Hlozek, R., Sievers, J., et al. 2011, ApJ, 739, 52 [NASA ADS] [CrossRef] [Google Scholar]
 Dunkley, J., Calabrese, E., Sievers, J., et al. 2013, JCAP, 7, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Dvali, G. R., Shafi, Q., & Schaefer, R. K. 1994, Phys. Rev. Lett., 73, 1886 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Easther, R., & Peiris, H. 2006, JCAP, 0609, 010 [NASA ADS] [Google Scholar]
 Easther, R., & Peiris, H. V. 2012, Phys. Rev. D, 85, 103533 [NASA ADS] [CrossRef] [Google Scholar]
 Easther, R., Greene, B. R., Kinney, W. H., & Shiu, G. 2001, Phys. Rev. D, 64, 103502 [NASA ADS] [CrossRef] [Google Scholar]
 Easther, R., Flauger, R., & Gilmore, J. B. 2011, JCAP, 1104, 027 [NASA ADS] [CrossRef] [Google Scholar]
 Efstathiou, G., & Bond, J. R. 1986, MNRAS, 218, 103 [NASA ADS] [CrossRef] [Google Scholar]
 Efstathiou, G., & Bond, J. R. 1987, MNRAS, 227, 33 [NASA ADS] [CrossRef] [Google Scholar]
 Elgarøy, Ø., Hannestad, S., & Haugbølle, T. 2003, JCAP, 0309, 008 [NASA ADS] [CrossRef] [Google Scholar]
 Enqvist, K., & Sloth, M. S. 2002, Nucl. Phys. B, 626, 395 [NASA ADS] [CrossRef] [Google Scholar]
 Eriksen, H. K., Huey, G., Saha, R., et al. 2007, ApJ, 656, 641 [NASA ADS] [CrossRef] [Google Scholar]
 Fabbri, R., & Pollock, M. 1983, Phys. Lett. B, 125, 445 [NASA ADS] [CrossRef] [Google Scholar]
 Fakir, R., & Unruh, W. 1990, Phys. Rev. D, 41, 1783 [NASA ADS] [CrossRef] [Google Scholar]
 Feroz, F., & Hobson, M. 2008, MNRAS, 384, 449 [NASA ADS] [CrossRef] [Google Scholar]
 Feroz, F., Hobson, M., & Bridges, M. 2009, MNRAS, 398, 1601 [NASA ADS] [CrossRef] [Google Scholar]
 Feroz, F., Hobson, M. P., Cameron, E., & Pettitt, A. N. 2013 [arXiv:1306.2144] [Google Scholar]
 Finelli, F., Hamann, J., Leach, S. M., & Lesgourgues, J. 2010, JCAP, 1004, 011 [NASA ADS] [CrossRef] [Google Scholar]
 Fixsen, D. 2009, ApJ, 707, 916 [NASA ADS] [CrossRef] [Google Scholar]
 Flauger, R., McAllister, L., Pajer, E., Westphal, A., & Xu, G. 2010, JCAP, 1006, 009 [NASA ADS] [CrossRef] [Google Scholar]
 Freese, K., Frieman, J. A., & Olinto, A. V. 1990, Phys. Rev. Lett., 65, 3233 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Freivogel, B., Kleban, M., Rodríguez Martínez, M., & Susskind, L. 2006, J. High Energy Phys., 0603, 039 [NASA ADS] [CrossRef] [Google Scholar]
 GarcíaBellido, J., & Wands, D. 1996, Phys. Rev. D, 53, 5437 [NASA ADS] [CrossRef] [Google Scholar]
 Garriga, J., & Mukhanov, V. F. 1999, Phys. Lett. B, 458, 219 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Garriga, J., Montes, X., Sasaki, M., & Tanaka, T. 1998, Nucl. Phys. B, 513, 343 [NASA ADS] [CrossRef] [Google Scholar]
 Garriga, J., Montes, X., Sasaki, M., & Tanaka, T. 1999, Nucl. Phys. B, 551, 317 [NASA ADS] [CrossRef] [Google Scholar]
 Gauthier, C., & Bucher, M. 2012, JCAP, 1210, 050 [NASA ADS] [CrossRef] [Google Scholar]
 Gong, J.O., & Stewart, E. D. 2001, Phys. Lett. B, 510, 1 [NASA ADS] [Google Scholar]
 Gordon, C., & Lewis, A. 2003, Phys. Rev. D, 67, 123513 [NASA ADS] [CrossRef] [Google Scholar]
 Gordon, C., & Malik, K. A. 2004, Phys. Rev. D, 69, 063508 [NASA ADS] [CrossRef] [Google Scholar]
 Gordon, C., Wands, D., Bassett, B. A., & Maartens, R. 2001, Phys. Rev. D, 63, 023506 [NASA ADS] [CrossRef] [Google Scholar]
 Gott, J. 1982, Nature, 295, 304 [NASA ADS] [CrossRef] [Google Scholar]
 Gott, J., & Statler, T. 1984, Phys. Lett. B, 136, 157 [NASA ADS] [CrossRef] [Google Scholar]
 Gratton, S., & Turok, N. 1999, Phys. Rev. D, 60, 123507 [NASA ADS] [CrossRef] [Google Scholar]
 Gratton, S., Hertog, T., & Turok, N. 2000, Phys. Rev. D, 62, 063501 [NASA ADS] [CrossRef] [Google Scholar]
 Gratton, S., Lewis, A., & Turok, N. 2002, Phys. Rev. D, 65, 043513 [NASA ADS] [CrossRef] [Google Scholar]
 Grishchuk, L. 1975, Sov. Phys. JETP, 40, 409 [NASA ADS] [Google Scholar]
 Groot Nibbelink, S., & van Tent, B. 2000, unpublished [arXiv:hepph/0011325] [Google Scholar]
 Groot Nibbelink, S., & van Tent, B. 2002, Class. Quant. Grav., 19, 613 [NASA ADS] [CrossRef] [Google Scholar]
 Guth, A. H. 1981, Phys. Rev. D, 23, 347 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Guth, A. H., & Nomura, Y. 2012, Phys. Rev. D, 86, 023534 [NASA ADS] [CrossRef] [Google Scholar]
 Guth, A. H., & Pi, S. 1982, Phys. Rev. Lett., 49, 1110 [NASA ADS] [CrossRef] [Google Scholar]
 Habib, S., Heitmann, K., Jungman, G., & MolinaParis, C. 2002, Phys. Rev. Lett., 89, 281301 [NASA ADS] [CrossRef] [Google Scholar]
 Hamann, J., Hannestad, S., Sloth, M. S., & Wong, Y. Y. Y. 2008a, JCAP, 0809, 015 [NASA ADS] [CrossRef] [Google Scholar]
 Hamann, J., Lesgourgues, J., & Mangano, G. 2008b, JCAP, 0803, 004 [NASA ADS] [CrossRef] [Google Scholar]
 Hamann, J., Lesgourgues, J., & Valkenburg, W. 2008c, JCAP, 0804, 016 [NASA ADS] [CrossRef] [Google Scholar]
 Hamann, J., Hannestad, S., Raffelt, G. G., & Wong, Y. Y. Y. 2009, JCAP, 0906, 022 [NASA ADS] [CrossRef] [Google Scholar]
 Hamann, J., Shafieloo, A., & Souradeep, T. 2010, JCAP, 1004, 010 [NASA ADS] [CrossRef] [Google Scholar]
 Hamimeche, S., & Lewis, A. 2008, Phys. Rev. D, 77, 103013 [NASA ADS] [CrossRef] [Google Scholar]
 Hannestad, S. 2004, JCAP, 0404, 002 [NASA ADS] [CrossRef] [Google Scholar]
 Harrison, E. R. 1970, Phys. Rev. D, 1, 2726 [NASA ADS] [CrossRef] [Google Scholar]
 Hawking, S. 1982, Phys. Lett. B, 115, 295 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Hawking, S., & Turok, N. 1998, Phys. Lett. B, 425, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Hertog, T., & Turok, N. 2000, Phys. Rev. D, 62, 083514 [NASA ADS] [CrossRef] [Google Scholar]
 Hinshaw, G., Larson, D., Komatsu, E., et al. 2013, ApJS, 208, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Hou, Z., Keisler, R., Knox, L., Millea, M., & Reichardt, C. 2013, Phys. Rev. D, 87, 083008 [NASA ADS] [CrossRef] [Google Scholar]
 Hou, Z., Reichardt, C., Story, K., et al. 2014, ApJ, 782, 74 [NASA ADS] [CrossRef] [Google Scholar]
 Hunt, P., & Sarkar, S. 2004, Phys. Rev. D, 70, 103518 [NASA ADS] [CrossRef] [Google Scholar]
 Ichikawa, K., & Takahashi, T. 2006, Phys. Rev. D, 73, 063528 [NASA ADS] [CrossRef] [Google Scholar]
 Ichiki, K., & Nagata, R. 2009, Phys. Rev. D, 80, 083002 [NASA ADS] [CrossRef] [Google Scholar]
 Jaynes, E. T., & Bretthorst, G. L. 2003, Probability Theory, Cambridge, UK: Probability Theory, eds. E. T. Jaynes & G. L. Bretthorst (Cambridge University Press) [Google Scholar]
 Jeffreys, H. 1998, Theory of Probability, 3rd edn. (Oxford University Press) [Google Scholar]
 Kazanas, D. 1980, ApJ, 241, L59 [NASA ADS] [CrossRef] [Google Scholar]
 Keskitalo, R., KurkiSuonio, H., Muhonen, V., & Valiviita, J. 2007, JCAP, 0709, 008 [NASA ADS] [CrossRef] [Google Scholar]
 Kim, J. E., & Carosi, G. 2010, Rev. Mod. Phys., 82, 557 [NASA ADS] [CrossRef] [Google Scholar]
 Kinney, W. H. 2002, Phys. Rev. D, 66, 083508 [NASA ADS] [CrossRef] [Google Scholar]
 Kinney, W. H., & Riotto, A. 2006, JCAP, 3, 11 [NASA ADS] [CrossRef] [Google Scholar]
 Kinney, W. H., Kolb, E. W., Melchiorri, A., & Riotto, A. 2006, Phys. Rev. D, 74, 023502 [NASA ADS] [CrossRef] [Google Scholar]
 Kleban, M., & Schillo, M. 2012, JCAP, 1206, 029 [NASA ADS] [CrossRef] [Google Scholar]
 Knox, L., & Turner, M. S. 1994, Phys. Rev. Lett., 73, 3347 [NASA ADS] [CrossRef] [Google Scholar]
 Kobayashi, T., & Takahashi, F. 2011, JCAP, 1101, 026 [NASA ADS] [CrossRef] [Google Scholar]
 Kobayashi, T., Yamaguchi, M., & Yokoyama, J. 2010, Phys. Rev. Lett., 105, 231302 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Kofman, L. 1986, Phys. Lett. B, 173, 400 [NASA ADS] [CrossRef] [Google Scholar]
 Kofman, L., Linde, A. D., & Starobinsky, A. A. 1994, Phys. Rev. Lett., 73, 3195 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Kofman, L., Linde, A. D., & Starobinsky, A. A. 1997, Phys. Rev. D, 56, 3258 [NASA ADS] [CrossRef] [Google Scholar]
 Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, ApJS, 180, 330 [NASA ADS] [CrossRef] [Google Scholar]
 Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Kosowsky, A., & Turner, M. S. 1995, Phys. Rev. D, 52, 1739 [NASA ADS] [CrossRef] [Google Scholar]
 KurkiSuonio, H., Muhonen, V., & Valiviita, J. 2005, Phys. Rev. D, 71, 063005 [NASA ADS] [CrossRef] [Google Scholar]
 Langlois, D. 1999, Phys. Rev. D, 59, 123512 [NASA ADS] [CrossRef] [Google Scholar]
 Langlois, D., & Riazuelo, A. 2000, Phys. Rev. D, 62, 043504 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Larson, D., Dunkley, J., Hinshaw, G., et al. 2011, ApJS, 192, 16 [NASA ADS] [CrossRef] [Google Scholar]
 Leach, S. M., Liddle, A. R., Martin, J., & Schwarz, D. J. 2002, Phys. Rev. D, 66, 023515 [NASA ADS] [CrossRef] [Google Scholar]
 Lesgourgues, J. 2011 [arXiv:1104.2932] [Google Scholar]
 Lesgourgues, J., & Valkenburg, W. 2007, Phys. Rev. D, 75, 123519 [NASA ADS] [CrossRef] [Google Scholar]
 Lesgourgues, J., Starobinsky, A. A., & Valkenburg, W. 2008, JCAP, 0801, 010 [NASA ADS] [CrossRef] [Google Scholar]
 Lewis, A. 2011, http://cosmologist.info/notes/CAMB.pdf [Google Scholar]
 Lewis, A., & Bridle, S. 2002, Phys. Rev. D, 66, 103511 [NASA ADS] [CrossRef] [Google Scholar]
 Liddle, A. R. 2007, MNRAS, 377, L74 [NASA ADS] [Google Scholar]
 Liddle, A. R., & Lyth, D. H. 1993, Phys. Rept., 231, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Liddle, A. R., & Leach, S. M. 2003, Phys. Rev. D, 68, 103503 [NASA ADS] [CrossRef] [Google Scholar]
 Lidsey, J. E., & Huston, I. 2007, JCAP, 7, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Lifshitz, E. 1946, J. Phys. (USSR), 10, 116 [Google Scholar]
 Lifshitz, E., & Khalatnikov, I. 1963, Adv. Phys., 12, 185 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Linde, A., Noorbala, M., & Westphal, A. 2011, JCAP, 1103, 013 [NASA ADS] [CrossRef] [Google Scholar]
 Linde, A. D. 1982, Phys. Lett. B, 108, 389 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Linde, A. D. 1983, Phys. Lett. B,129, 177 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Linde, A. D. 1984, JETP Lett., 40, 1333 [NASA ADS] [Google Scholar]
 Linde, A. D. 1985, Phys. Lett. B, 158, 375 [NASA ADS] [CrossRef] [Google Scholar]
 Linde, A. D. 1990, Particle physics and inflationary cosmology (Harwood) [Google Scholar]
 Linde, A. D. 1991, Phys. Lett. B, 259, 38 [NASA ADS] [CrossRef] [Google Scholar]
 Linde, A. D. 2003, JCAP, 0305, 002 [NASA ADS] [CrossRef] [Google Scholar]
 Linde, A. D., & Lyth, D. H. 1990, Phys. Lett. B, 246, 353 [NASA ADS] [CrossRef] [Google Scholar]
 Linde, A. D., & Mukhanov, V. F. 1997, Phys. Rev. D, 56, 535 [NASA ADS] [CrossRef] [Google Scholar]
 Lorenz, L., Martin, J., & Ringeval, C. 2008, Phys. Rev. D, 78, 083513 [NASA ADS] [CrossRef] [Google Scholar]
 Lucchin, F., & Matarrese, S. 1985, Phys. Rev. D, 32, 1316 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Lucchin, F., Matarrese, S., & Pollock, M. 1986, Phys. Lett. B, 167, 163 [NASA ADS] [CrossRef] [Google Scholar]
 Lucy, L. B. 1974, AJ, 79, 745 [NASA ADS] [CrossRef] [Google Scholar]
 Lyth, D. 1992, Phys. Rev. D, 45, 3394 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Lyth, D. H. 1990, Phys. Lett. B, 236, 408 [NASA ADS] [CrossRef] [Google Scholar]
 Lyth, D. H. 1997, Phys. Rev. Lett., 78, 1861 [NASA ADS] [CrossRef] [Google Scholar]
 Lyth, D. H., & Riotto, A. 1999, Phys. Rept., 314, 1 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Lyth, D. H., & Stewart, E. D. 1990, Phys. Lett. B, 252, 336 [NASA ADS] [CrossRef] [Google Scholar]
 Lyth, D. H., & Stewart, E. D. 1996, Phys. Rev. D, 53, 1784 [NASA ADS] [CrossRef] [Google Scholar]
 Lyth, D. H., Ungarelli, C., & Wands, D. 2003, Phys. Rev. D, 67, 023503 [NASA ADS] [CrossRef] [Google Scholar]
 Lyth, D. H., & Wands, D. 2002, Phys. Lett. B, 524, 5 [NASA ADS] [CrossRef] [Google Scholar]
 Ma, C.P., & Bertschinger, E. 1995, ApJ, 455, 7 [NASA ADS] [CrossRef] [Google Scholar]
 Maldacena, J. 2003, J. High Energy Phys., 5, 13 [NASA ADS] [CrossRef] [Google Scholar]
 Martin, J., & Brandenberger, R. 2003, Phys. Rev. D, 68, 063513 [NASA ADS] [CrossRef] [Google Scholar]
 Martin, J., & Ringeval, C. 2004, Phys. Rev. D, 69, 083515 [NASA ADS] [CrossRef] [Google Scholar]
 Martin, J., & Ringeval, C. 2010, Phys. Rev. D, 82, 023511 [NASA ADS] [CrossRef] [Google Scholar]
 Martin, J., & Schwarz, D. J. 2003, Phys. Rev. D, 67, 083512 [NASA ADS] [CrossRef] [Google Scholar]
 Martin, J., Ringeval, C., & Trotta, R. 2011, Phys. Rev. D, 83, 063524 [NASA ADS] [CrossRef] [Google Scholar]
 McAllister, L., Silverstein, E., & Westphal, A. 2010, Phys. Rev. D, 82, 046003 [NASA ADS] [CrossRef] [Google Scholar]
 Meerburg, P. D., Wijers, R. A. M. J., & van der Schaar, J. P. 2012, MNRAS, 421, 369 [NASA ADS] [Google Scholar]
 Mehta, K. T., Cuesta, A. J., Xu, X., Eisenstein, D. J., & Padmanabhan, N. 2012, MNRAS, 427, 2168 [NASA ADS] [CrossRef] [Google Scholar]
 Mollerach, S. 1990, Phys. Rev. D, 42, 313 [NASA ADS] [CrossRef] [Google Scholar]
 Moodley, K., Bucher, M., Dunkley, J., Ferreira, P., & Skordis, C. 2004, Phys. Rev. D, 70, 103520 [NASA ADS] [CrossRef] [Google Scholar]
 Moroi, T., & Takahashi, T. 2001, Phys. Lett. B, 522, 215 [NASA ADS] [CrossRef] [Google Scholar]
 Mortonson, M. J., & Hu, W. 2008a, ApJ, 672, 737 [NASA ADS] [CrossRef] [Google Scholar]
 Mortonson, M. J., & Hu, W. 2008b, ApJ, 686, L53 [NASA ADS] [CrossRef] [Google Scholar]
 Mortonson, M. J., Dvorkin, C., Peiris, H. V., & Hu, W. 2009, Phys. Rev. D, 79, 103519 [NASA ADS] [CrossRef] [Google Scholar]
 Mortonson, M. J., Peiris, H. V., & Easther, R. 2011, Phys. Rev. D, 83, 043505 [NASA ADS] [CrossRef] [Google Scholar]
 Mukhanov, V. F. 1985, JETP Lett., 41, 493 [NASA ADS] [Google Scholar]
 Mukhanov, V. F. 1988, Sov. Phys. JETP, 67, 1297 [Google Scholar]
 Mukhanov, V. F., & Chibisov, G. 1981, JETP Lett., 33, 532 [NASA ADS] [Google Scholar]
 Mukhanov, V. F., & Chibisov, G. 1982, Sov. Phys. JETP, 56, 258 [Google Scholar]
 Mukhanov, V. F., Feldman, H., & Brandenberger, R. H. 1992, Phys. Rept., 215, 203 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Muslimov, A. 1990, Class. Quant. Grav., 7, 231 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Nagata, R., & Yokoyama, J. 2008, Phys. Rev. D, 78, 123002 [NASA ADS] [CrossRef] [Google Scholar]
 Nagata, R., & Yokoyama, J. 2009, Phys. Rev. D, 79, 043010 [NASA ADS] [CrossRef] [Google Scholar]
 Norena, J., Wagner, C., Verde, L., Peiris, H. V., & Easther, R. 2012, Phys. Rev. D, 86, 023505 [NASA ADS] [CrossRef] [Google Scholar]
 Okada, N., Rehman, M. U., & Shafi, Q. 2010, Phys. Rev. D, 82, 043502 [NASA ADS] [CrossRef] [Google Scholar]
 Okamoto, T., & Hu, W. 2003, Phys. Rev. D, 67, 083002 [NASA ADS] [CrossRef] [Google Scholar]
 Olive, K. A. 1990, Phys. Rept., 190, 307 [NASA ADS] [CrossRef] [Google Scholar]
 Padmanabhan, N., Xu, X., Eisenstein, D. J., et al. 2012, MNRAS, 427, 2132 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Page, L., Hinshaw, G., Komatsu, E., et al. 2007, ApJS, 170, 335 [NASA ADS] [CrossRef] [Google Scholar]
 Pahud, C., Liddle, A. R., Mukherjee, P., & Parkinson, D. 2007, MNRAS, 381, 489 [NASA ADS] [CrossRef] [Google Scholar]
 Pallis, C. 2006, Nucl. Phys. B, 751, 129 [NASA ADS] [CrossRef] [Google Scholar]
 Pandolfi, S., Giusarma, E., Kolb, E. W., et al. 2010, Phys. Rev. D, 82, 123527 [NASA ADS] [CrossRef] [Google Scholar]
 Peccei, R. 2008, Lect. Notes Phys., 741, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Peebles, P. J. E., & Yu, J. T. 1970, ApJ, 162, 815 [NASA ADS] [CrossRef] [Google Scholar]
 Peiris, H., & Easther, R. 2006a, JCAP, 0607, 002 [NASA ADS] [CrossRef] [Google Scholar]
 Peiris, H., & Easther, R. 2006b, JCAP, 0610, 017 [NASA ADS] [CrossRef] [Google Scholar]
 Peiris, H., Baumann, D., Friedman, B., & Cooray, A. 2007, Phys. Rev. D, 76, 103517 [NASA ADS] [CrossRef] [Google Scholar]
 Peiris, H., Easther, R., & Flauger, R. 2013, JCAP, 09, 018 [NASA ADS] [CrossRef] [Google Scholar]
 Peiris, H. V., & Easther, R. 2008, JCAP, 0807, 024 [NASA ADS] [CrossRef] [Google Scholar]
 Peiris, H. V., & Verde, L. 2010, Phys. Rev. D, 81, 021302 [NASA ADS] [CrossRef] [Google Scholar]
 Peiris, H. V., Komatsu, E., Verde, L., et al. 2003, ApJS, 148, 213 [NASA ADS] [CrossRef] [Google Scholar]
 Planck Collaboration I. 2014, A&A, 571, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration II. 2014, A&A, 571, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration III. 2014, A&A, 571, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration IV. 2014, A&A, 571, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration V. 2014, A&A, 571, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VI. 2014, A&A, 571, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VII. 2014, A&A, 571, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VIII. 2014, A&A, 571, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration IX. 2014, A&A, 571, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration X. 2014, A&A, 571, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XI. 2014, A&A, 571, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XII. 2014, A&A, 571, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIII. 2014, A&A, 571, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIV. 2014, A&A, 571, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XV. 2014, A&A, 571, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVI. 2014, A&A, 571, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVII. 2014, A&A, 571, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVIII. 2014, A&A, 571, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIX. 2014, A&A, 571, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XX. 2014, A&A, 571, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXI. 2014, A&A, 571, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXII. 2014, A&A, 571, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXIII. 2014, A&A, 571, A23 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Planck Collaboration XXIV. 2014, A&A, 571, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXV. 2014, A&A, 571, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXVI. 2014, A&A, 571, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXVII. 2014, A&A, 571, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXVIII. 2014, A&A, 571, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXIX. 2014, A&A, 571, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXX. 2014, A&A, 571, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXXI. 2014, A&A, 571, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Polarski, D., & Starobinsky, A. A. 1994, Phys. Rev. D, 50, 6123 [NASA ADS] [CrossRef] [Google Scholar]
 Powell, B. A., & Kinney, W. H. 2007, JCAP, 0708, 006 [NASA ADS] [CrossRef] [Google Scholar]
 Powell, B. A., Tzirakis, K., & Kinney, W. H. 2009, JCAP, 4, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Preskill, J., Wise, M. B., & Wilczek, F. 1983, Phys. Lett. B, 120, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Raffelt, G. G. 2008, Lect. Notes Phys., 741, 51 [NASA ADS] [CrossRef] [Google Scholar]
 Ratra, B., & Peebles, P. 1994, ApJ, 432, L5 [NASA ADS] [CrossRef] [Google Scholar]
 Ratra, B., & Peebles, P. 1995, Phys. Rev. D, 52, 1837 [NASA ADS] [CrossRef] [Google Scholar]
 Reichardt, C., Shaw, L., Zahn, O., et al. 2012, ApJ, 755, 70 [NASA ADS] [CrossRef] [Google Scholar]
 Richardson, W. H. 1972, J. Opt. Soc. Am., 62, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Riess, A. G., Macri, L., Casertano, S., et al. 2011, ApJ, 730, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Rubakov, V., Sazhin, M., & Veryaskin, A. 1982, Phys. Lett. B, 115, 189 [NASA ADS] [CrossRef] [Google Scholar]
 Salopek, D., Bond, J., & Bardeen, J. M. 1989, Phys. Rev. D, 40, 1753 [NASA ADS] [CrossRef] [Google Scholar]
 Sasaki, M. 1986, Prog. Theor. Phys., 76, 1036 [NASA ADS] [CrossRef] [Google Scholar]
 Sasaki, M., & Stewart, E. D. 1996, Prog. Theor. Phys., 95, 71 [NASA ADS] [CrossRef] [Google Scholar]
 Sasaki, M., Tanaka, T., & Yakushige, Y. 1997, Phys. Rev. D, 56, 616 [NASA ADS] [CrossRef] [Google Scholar]
 Sasaki, M., Valiviita, J., & Wands, D. 2006, Phys. Rev. D, 74, 103003 [NASA ADS] [CrossRef] [Google Scholar]
 Sato, K. 1981, MNRAS, 195, 467 [NASA ADS] [CrossRef] [Google Scholar]
 Savage, C., Freese, K., & Kinney, W. H. 2006, Phys. Rev. D, 74, 123511 [NASA ADS] [CrossRef] [Google Scholar]
 Schwarz, G. 1978, Ann. Stat., 6, 461 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Sealfon, C., Verde, L., & Jimenez, R. 2005, Phys. Rev. D, 72, 103520 [NASA ADS] [CrossRef] [Google Scholar]
 Seckel, D., & Turner, M. S. 1985, Phys. Rev. D, 32, 3178 [NASA ADS] [CrossRef] [Google Scholar]
 Senatore, L., Smith, K. M., & Zaldarriaga, M. 2010, JCAP, 1, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Shafieloo, A., & Souradeep, T. 2004, Phys. Rev. D, 70, 043523 [NASA ADS] [CrossRef] [Google Scholar]
 Shafieloo, A., & Souradeep, T. 2008, Phys. Rev. D, 78, 023511 [NASA ADS] [CrossRef] [Google Scholar]
 Sievers, J. L., Hlozek, R. A., Nolta, M. R., et al. 2013, JCAP, 10, 060 [NASA ADS] [CrossRef] [Google Scholar]
 Sikivie, P. 2008, Lect. Notes Phys., 741, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Silverstein, E., & Tong, D. 2004, Phys. Rev. D, 70, 103505 [NASA ADS] [CrossRef] [Google Scholar]
 Silverstein, E., & Westphal, A. 2008, Phys. Rev. D, 78, 106003 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Sinha, R., & Souradeep, T. 2006, Phys. Rev. D, 74, 043518 [NASA ADS] [CrossRef] [Google Scholar]
 Spergel, D. N., Verde, L., Peiris, H. V., et al. 2003, ApJS, 148, 175 [NASA ADS] [CrossRef] [Google Scholar]
 Spergel, D. N., Bean, R., Doré, O., et al. 2007, ApJS, 170, 377 [NASA ADS] [CrossRef] [Google Scholar]
 Spokoiny, B. 1984, Phys. Lett. B, 147, 39 [NASA ADS] [CrossRef] [Google Scholar]
 Spokoiny, B. 1993, Phys. Lett. B, 315, 40 [NASA ADS] [CrossRef] [Google Scholar]
 Starobinsky, A. A. 1979, JETP Lett., 30, 682 [NASA ADS] [Google Scholar]
 Starobinsky, A. A. 1980, Phys. Lett. B, 91, 99 [NASA ADS] [CrossRef] [Google Scholar]
 Starobinsky, A. A. 1982, Phys. Lett. B, 117, 175 [NASA ADS] [CrossRef] [Google Scholar]
 Starobinsky, A. A. 1983, Sov. Astron. Lett., 9, 302 [NASA ADS] [Google Scholar]
 Starobinsky, A. A. 1985a, Sov. Astron. Lett., 11, 133 [NASA ADS] [Google Scholar]
 Starobinsky, A. A. 1985b, JETP Lett., 42, 152 [NASA ADS] [Google Scholar]
 Starobinsky, A. A. 1992, JETP Lett., 55, 489 [NASA ADS] [Google Scholar]
 Starobinsky, A. A. 2005, JETP Lett., 82, 169 [NASA ADS] [CrossRef] [Google Scholar]
 Steinhardt, P. J., & Turner, M. S. 1983, Phys. Lett. B, 129, 51 [NASA ADS] [CrossRef] [Google Scholar]
 Stewart, E. D., & Lyth, D. H. 1993, Phys. Lett. B, 302, 171 [NASA ADS] [CrossRef] [Google Scholar]
 Stompor, R., Banday, A. J., & Gorski, K. M. 1996, ApJ, 463, 8 [NASA ADS] [CrossRef] [Google Scholar]
 Story, K., Reichardt, C., Hou, Z., et al. 2013, ApJ, 779, 86 [NASA ADS] [CrossRef] [Google Scholar]
 Suzuki, N., Rubin, D., Lidman, C., et al. 2012, ApJ, 746, 85 [NASA ADS] [CrossRef] [Google Scholar]
 Tanaka, T., & Sasaki, M. 1994, Phys. Rev. D, 50, 6444 [NASA ADS] [CrossRef] [Google Scholar]
 TocchiniValentini, D., Douspis, M., & Silk, J. 2005, MNRAS, 359, 31 [NASA ADS] [CrossRef] [Google Scholar]
 TocchiniValentini, D., Hoffman, Y., & Silk, J. 2006, MNRAS, 367, 1095 [NASA ADS] [CrossRef] [Google Scholar]
 Traschen, J. H., & Brandenberger, R. H. 1990, Phys. Rev. D, 42, 2491 [NASA ADS] [CrossRef] [Google Scholar]
 Trotta, R. 2007, MNRAS, 375, L26 [NASA ADS] [Google Scholar]
 Trotta, R., & Hansen, S. H. 2004, Phys. Rev. D, 69, 023509 [NASA ADS] [CrossRef] [Google Scholar]
 Tsujikawa, S., & Gumjudpai, B. 2004, Phys. Rev. D, 69, 123523 [NASA ADS] [CrossRef] [Google Scholar]
 Turner, M. S. 1983, Phys. Rev. D, 28, 1243 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Turner, M. S. 1990, Phys. Rept., 197, 67 [NASA ADS] [CrossRef] [Google Scholar]
 Turner, M. S., & Wilczek, F. 1991, Phys. Rev. Lett., 66, 5 [NASA ADS] [CrossRef] [Google Scholar]
 Turner, M. S., Wilczek, F., & Zee, A. 1983, Phys. Lett. B, 125, 35 [NASA ADS] [CrossRef] [Google Scholar]
 Valiviita, J., & Muhonen, V. 2003, Phys. Rev. Lett., 91, 131302 [NASA ADS] [CrossRef] [Google Scholar]
 Valiviita, J., Savelainen, M., Talvitie, M., KurkiSuonio, H., & Rusak, S. 2012, ApJ, 753, 151 [NASA ADS] [CrossRef] [Google Scholar]
 Vázquez, J. A., Bridges, M., Hobson, M., & Lasenby, A. 2012, JCAP, 06, 006 [NASA ADS] [CrossRef] [Google Scholar]
 Verde, L., & Peiris, H. V. 2008, JCAP, 0807, 009 [NASA ADS] [CrossRef] [Google Scholar]
 Vilenkin, A. 2007, J. Phys. A, 40, 6777 [NASA ADS] [CrossRef] [Google Scholar]
 Weinberg, S. 2013, Phys. Rev. Lett., 110, 241301 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Yamamoto, K., Sasaki, M., & Tanaka, T. 1995, ApJ, 455, 412 [NASA ADS] [CrossRef] [Google Scholar]
 Zeldovich, Y. 1972, MNRAS, 160, 1 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Sampling the Hubble flow functions
Fig. A.1 Planck constraints on the HFFs (ϵ_{1},ϵ_{2},ϵ_{3}) assuming either ϵ_{3} = 0 and the firstorder slowroll approximation for the computation of the primordial spectra, or ϵ_{3} ≠ 0 and the secondorder slowroll approximation (HFF GFM). 

Open with DEXTER 
Fig. A.2 Comparison of the Planck constraints on the HFFs (ϵ_{1},ϵ_{2},ϵ_{3}) using the GFM and the MCE. 

Open with DEXTER 
Fig. A.3 Planck constraints on the spectral parameters n_{s},dn_{s}/ dlnk, and r. We compare constraints computed with the secondorder slowroll approximation, starting from flat priors on the HFF parameters at the pivot scale (ϵ_{1},ϵ_{2},ϵ_{3}), with those obtained directly from n_{s},dn_{s}/ dlnk, and r. In the latter case we enforce the secondorder consistency conditions for the tensortoscalar ratio and for the running of the tensor spectral index. 

Open with DEXTER 
In this Appendix we briefly review how to constrain slowroll inflation by sampling the Hubble flow functions (HFFs) and discuss how well the results agree with those derived by sampling directly the parameters lnA_{s}, n_{s}, r, and dn_{s}/ dlnk. This method fully exploits an analytic perturbative expansion in terms of the HFFs for the primordial spectra of cosmological fluctuations during slowroll inflation (Stewart & Lyth 1993; Gong & Stewart 2001; Leach et al. 2002), which selfconsistently extends to highest order the first terms presented in Eqs. (13)–(19). Since z′′/z in Eq. (7) and a′′/a in Eq. (9) can be rewritten exactly in terms of the Hubble flow functions, the HFF hierarchy, rather than the potential hierarchy, is best suited for this purpose. The slowroll analytic power spectra have been calculated up to second order using the Green’s function method (Gong & Stewart 2001; Leach et al. 2002). Other approximations are available in the literature, including WKB (Martin & Schwarz 2003), the uniform approximation (Habib et al. 2002), or the method of comparison equations (Casadio et al. 2006).
The dependence of the amplitudes in Eqs. (13) and (14) in terms of HFF is given by (A.1)where X = s,t, and b_{s0},b_{t0} are with C ≡ ln2 + γ_{E} − 2 ≈ − 0.7296 (γ_{E} is the EulerMascheroni constant). At second order the coefficients of the expansion depend on the particular approximation scheme: X = C and Δ_{s0} = Δ_{t0} = 0 apply for the Green’s function method (GFM; Gong & Stewart 2001; Leach et al. 2002), and X = D = 1 / 3 − ln3, Δ_{s0} = (D − C)(D + ln2) − 1 / 18, Δ_{t0} = 2D(D − C) − 1 / 9 apply for the method of comparison equations (MCE; Casadio et al. 2006). As predicted by the consistency relation, A_{h 0} = 16ϵ_{1}A_{ℛ ′}.
The full perturbative expressions up to second order in HFF for the spectral indices and the running of the indices are
We now constrain these parameters using the Planck+WP data. If we restrict ourselves to first order, we obtain ϵ_{1}< 0.0074 at 95% CL and at 68% CL. At second order with GFM, we obtain ϵ_{1}< 0.013 at 95% CL, and at 68% CL. The comparison of results at first and second order is shown in Fig. A.1. Different approximation schemes lead to small differences, as Fig. A.2 shows for the GFM versus the MCE. Figure A.3 shows the agreement between the physical parameters reconstructed from the HFF method and those directly sampled as in Sect. 4.3.
All Tables
Constraints on cosmological parameters and best fit −2Δln(ℒ) with respect to the standard ΛCDM model, using Planck+WP data, testing the significance of the deviation from the HZ model.
Constraints on the primordial perturbation parameters in the ΛCDM+tensor model from Planck combined with other data sets.
Constraints on the primordial perturbation parameters for ΛCDM + dn_{s}/ dlnk, ΛCDM + dn_{s}/ dlnk + r, and ΛCDM + dn_{s}/ dlnk + d^{2}n_{s}/ dlnk^{2} models from Planck combined with other data sets.
Numerical reconstruction of potential parameters, compared to results with the slowroll approximation, when tensors and running are included (Planck+WP 95% CL, with k_{∗} = 0.05 Mpc^{1}).
Scalar amplitude, tilt, running, running of the running, and tensortoscalar ratio inferred from a numerical reconstruction of the inflaton potential (Planck+WP 95% CL, with k_{∗} = 0.05 Mpc^{1}).
Prior ranges imposed for the wiggles, step inflation, and cutoff model parameters.
Improvement in fit and logarithm of the Bayes factor B_{0X} with respect to power law ΛCDM and best fit parameter values for the wiggles, step inflation, and cutoff models (the larger lnB_{0X}, the greater the preference for a featureless power law spectrum).
All Figures
Fig. 1 Marginalized joint 68% and 95% CL regions for n_{s} and r_{0.002} from Planck in combination with other data sets compared to the theoretical predictions of selected inflationary models. 

Open with DEXTER  
In the text 
Fig. 2 Marginalized joint 68% and 95% CL for (dn_{s}/ dlnk,n_{s}) using Planck+WP+BAO, either marginalizing over r or fixing r = 0 at k_{∗} = 0.038 Mpc^{1}. The purple strip shows the prediction for single monomial chaotic inflationary models with 50 <N_{∗}< 60 for comparison. 

Open with DEXTER  
In the text 
Fig. 3 Marginalized joint 68% and 95% CL regions for (d^{2}n_{s}/ dlnk^{2}, dn_{s}/ dlnk) using Planck+WP+BAO. 

Open with DEXTER  
In the text 
Fig. 4 Marginalized joint 68% and 95% CL regions for (r,n_{s}), using Planck+WP+BAO with and without a running spectral index. 

Open with DEXTER  
In the text 
Fig. 5 Joint posterior for tensors and running of the scalar spectral index marginalizing over other parameters. The contours are set at 68% and 95%. The red contours apply for Planck+WP+highℓ data. The colour of the scattered points indicates the distribution of Ω_{K}. The blue contours apply when BAO data is also included. (Ω_{K} is then found to be well constrained close to zero.) The dashed vertical line shows the norunning solution. 

Open with DEXTER  
In the text 
Fig. 6 Marginalized joint 68% and 95% CL regions for Planck+WP+BAO data for ΛCDM + r + N_{eff} and ΛCDM + r + Y_{P} (top); and ΛCDM + r + ∑ m_{ν} and ΛCDM + r + w (bottom). Shown for comparison are the corresponding contours for the ΛCDM + r model. 

Open with DEXTER  
In the text 
Fig. 7 Marginalized joint 68% and 95% CL regions for Planck+ WP data for ΛCDM + r for instantaneous and general reionization. 

Open with DEXTER  
In the text 
Fig. 8 Constraints on n_{s} vs. r at k_{∗} = 0.002 Mpc^{1} for the inflationary models considered (i.e., power law potentials with n = 2 / 3,1,2, and 4, natural inflation, and hilltop inflation), showing joint 68% and 95% CL. Blue and grey distributions correspond to the restrictive and permissive entropy generation scenarios, respectively. The instantaneous entropy generation case corresponds to the thin (red) contours in the natural and hilltop panels; for the single parameter models, this case corresponds to the lowestr extremity of the restrictive case. The difference between the natural inflation region in Fig. 1 and the natural inflation constraints shown here is due to the strong projection effect described in the text. 

Open with DEXTER  
In the text 
Fig. 9 Constraints on n_{s} vs. w_{int} at k_{∗} = 0.002 Mpc^{1} for the inflationary models considered, as in Fig. 8. The instantaneous entropy generation case (1) corresponds to w_{int} = 1 / 3. 

Open with DEXTER  
In the text 
Fig. 10 Constraints on log _{10}(λ) vs. w_{int} for the inflationary models considered, as in Figs. 8 and 9. 

Open with DEXTER  
In the text 
Fig. 11 Potential parameters for natural inflation and hilltop inflation, as in Figs. 8–10. On the natural inflation panel, instantaneous entropy generation corresponds to a thin diagonal along the top edge of entropy generation case (3). 

Open with DEXTER  
In the text 
Fig. 12 Posterior distribution for the first three slowroll potential parameters using Planck+WP data. In the n = 2 − 4 cases, the inflaton potential is expanded to nth order, and the spectrum is obtained by fitting the numerically computed power spectrum to the data, with no slowroll approximation and no assumption about the extrapolation of the potential outside the observable window. 

Open with DEXTER  
In the text 
Fig. 13 Posterior distribution for the Taylor expansion coefficients, V_{i}, of the inflaton potential. The potential is expanded to nth order, assuming a flat prior on ϵ_{V}, η_{V}, , and . The coefficients V_{i} are expressed in natural units (where The contours show only half of the allowed regions for potential parameters: the other half is symmetric, with opposite signs for V_{1} and V_{3}. 

Open with DEXTER  
In the text 
Fig. 14 Observable range of the best fitting inflaton potentials, when V(φ) is Taylor expanded to nth order around the pivot value φ_{∗}, in natural units (where assuming a flat prior on ϵ_{V}, η_{V}, , and , and using Planck+WP data. Potentials obtained under the transformation (φ − φ_{∗}) → (φ_{∗} − φ) leave the same observable signature and are also allowed. The sparsity of potentials with a small V_{0} = V(φ_{∗}) is explained by the flat prior on ϵ_{V} rather than on ln(V_{0}). In fact, V_{0} is unbounded from below. 

Open with DEXTER  
In the text 
Fig. 15 Planck primordial power spectrum feature search results. Top four panels: the reconstructed power spectrum at four values for the smoothing parameter λ. The red curves indicate the maximum likelihood configuration for the fractional deviation f(k) of the power spectrum relative to a power law fiducial model (with A_{s} = 2.20 × 10^{9} and n_{s} = 0.9603) for the penalized likelihood. The error bars have a width corresponding to the minimum reconstructible width (the minimum width for a Gaussian feature so that the mean square deviation of the expectation value of the reconstruction differs by less than 10%) and a vertical extent showing the 1σ and 2σ limits for the fractional deviation averaged over the box. The grey hashed regions at the far left and right show where the fixing prior (i.e., α) sets f(k) = 0. The inner grey regions show where the reconstruction bias is so great that the minimum reconstructible width is undefined. Mock features in this region produce reverberations over the entire interval. With λ = 10^{3} and λ = 10^{4}, we find statistically significant fluctuations around k ~ 0.1 Mpc^{1}. Lower panels: the 1σ error bars for three combinations of cosmological parameters at the four values of λ. The maximum likelihood value for the fiducial model is indicated by the dashed line for comparison. 

Open with DEXTER  
In the text 
Fig. 16 CMB multipole spectrum residuals for best fit primordial power spectrum reconstruction with smoothing parameter λ = 10^{3}. The panels show the C_{ℓ} spectrum residuals (compared to the best fit power law fiducial model represented by the horizontal straight dashed line) for the four auto and crossspectra included in the highℓ likelihood. Here The data points have been binned with Δℓ = 31 and foregrounds subtracted according to the best fit foreground parameters. The solid black line shows the CMB spectrum residual for the maximum likelihood primordial power spectrum reconstruction with λ = 10^{3}. 

Open with DEXTER  
In the text 
Fig. 17 CMB multipole spectrum residual for the primordial power spectrum test feature. The test feature (top) is set to the anomalously large deviation of the primordial power spectrum reconstruction for λ = 10^{3} in the interval 0.1Mpc^{1}<k< 0.15Mpc^{1}, and is zero elsewhere. Bottom: the angular spectrum corresponding to this feature. We observe a large dip at ℓ ≈ 1800. 

Open with DEXTER  
In the text 
Fig. 18 Top: best fit primordial spectrum of curvature perturbations for the power law (black), wiggles (red), step inflation (green), and cutoff (blue) models. Centre: residuals of the temperature angular power spectrum. Note that the scale of the vertical axis changes at ℓ = 50. Bottom: zoom of region around the first acoustic peak. 

Open with DEXTER  
In the text 
Fig. 19 Marginalized posterior probability (red) and profile −2Δlnℒ_{max} (black) for selected parameters of the wiggles (top), step inflation (middle), and cutoff model (bottom). 

Open with DEXTER  
In the text 
Fig. 20 Marginalized joint 68% and 95% CL for (ϵ_{1},ϵ_{2}) for Planck+WP data comparing the canonical Lagrangian case with c_{s} = 1 to the case of varying c_{s} with a uniform prior 0.02 <c_{s}< 1 derived from the Planck nonGaussianity measurements. 

Open with DEXTER  
In the text 
Fig. 21 C_{TT} anisotropy shapes for the three isocurvature modes. Top: the shapes of the CDM isocurvature mode, neutrino density isocurvature mode, and neutrino velocity isocurvature mode are shown together with the adiabatic mode. The modes have the same amplitude parameters ( for the adiabatic mode and for each isocurvature mode). Bottom: the narrower multipole range illustrates the relative phases of the acoustic oscillations for these modes. 

Open with DEXTER  
In the text 
Fig. 22 Two dimensional distributions for power in isocurvature modes using Planck+WP data. 

Open with DEXTER  
In the text 
Fig. 23 Fractional contribution of isocurvature modes to the power spectrum. We show the distributions α_{ℐℐ}(2,20), α_{ℛℐ}(2,20), α_{ℐℐ}(21,200), α_{ℛℐ}(21,200), α_{ℐℐ}(201,2500), and α_{ℛℐ}(201,2500), defined in Eq. (69), for the CDI, NDI, and NVI modes as constrained by the Planck+WP data. 

Open with DEXTER  
In the text 
Fig. 24 Temperature spectrum of best fit models with a mixture of adiabatic and isocurvature modes. Top: spectrum of the best fit mixed models relative to that of the pure adiabatic model. Bottom: zoom of the SachsWolfe plateau of the best fit temperature spectrum for each of the three cases plus the pure adiabatic model, shown together with Planck lowℓ data points. 

Open with DEXTER  
In the text 
Fig. A.1 Planck constraints on the HFFs (ϵ_{1},ϵ_{2},ϵ_{3}) assuming either ϵ_{3} = 0 and the firstorder slowroll approximation for the computation of the primordial spectra, or ϵ_{3} ≠ 0 and the secondorder slowroll approximation (HFF GFM). 

Open with DEXTER  
In the text 
Fig. A.2 Comparison of the Planck constraints on the HFFs (ϵ_{1},ϵ_{2},ϵ_{3}) using the GFM and the MCE. 

Open with DEXTER  
In the text 
Fig. A.3 Planck constraints on the spectral parameters n_{s},dn_{s}/ dlnk, and r. We compare constraints computed with the secondorder slowroll approximation, starting from flat priors on the HFF parameters at the pivot scale (ϵ_{1},ϵ_{2},ϵ_{3}), with those obtained directly from n_{s},dn_{s}/ dlnk, and r. In the latter case we enforce the secondorder consistency conditions for the tensortoscalar ratio and for the running of the tensor spectral index. 

Open with DEXTER  
In the text 