Early Cosmology Constrained
Abstract
We investigate our knowledge of early universe cosmology by exploring how much additional energy density can be placed in different components beyond those in the CDM model. To do this we use a method to separate early and lateuniverse information enclosed in observational data, thus markedly reducing the modeldependency of the conclusions. We find that the 95% credibility regions for extra energy components of the early universe at recombination are: nonaccelerating additional fluid density parameter and extra radiation parameterised as extra effective neutrino species when imposing flatness. Our constraints thus show that even when analyzing the data in this largely modelindependent way, the possibility of hiding extra energy components beyond CDM in the early universe is seriously constrained by current observations. We also find that the standard ruler, the sound horizon at radiation drag, can be well determined in a way that does not depend on latetime Universe assumptions, but depends strongly on earlytime physics and in particular on additional components that behave like radiation. We find that the standard ruler length determined in this way is Mpc if the radiation and neutrino components are standard, but the uncertainty increases by an order of magnitude when nonstandard dark radiation components are allowed, to Mpc.
a,b,c,d,e]Licia Verde, f,a]Emilio Bellini, g]Cassio Pigozzo, h]Alan F. Heavens, a,b,c,d]Raul Jimenez
Prepared for submission to JCAP
Early Cosmology Constrained

ICC, University of Barcelona, IEECUB, Martí Franquès, 1, E08028 Barcelona, Spain

ICREA, Pg. Lluis Companys 23, 08010 Barcelona, Spain

Radcliffe Institute for Advanced Study, Harvard University, MA 02138, USA

Institute for Theory and Computation, HarvardSmithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA

Institute of Theoretical Astrophysics, University of Oslo, 0315 Oslo, Norway

University of Oxford, Denys Wilkinson Building, Keble Road, Oxford, OX1 3RH, UK

Instituto de Física, Universidade Federal da Bahia, Salvador, BA, Brasil

Imperial Centre for Inference and Cosmology (ICIC), Imperial College, Blackett Laboratory, Prince Consort Road, London SW7 2AZ, U.K.
Contents
1 Introduction
The cosmic microwave background (CMB) is one of the most important cosmological probes and is the key observable in establishing the standard (CDM) cosmological model. Lower redshift observations such as clustering of largescale structure or probes of the expansion history (supernovae, baryon acoustic oscillations, cosmic chronometers), have been essential to confirm this picture and to constrain the properties of dark energy and other energy components of the Universe.
When fitting the CMB power spectrum (temperature and polarization) one usually has to make simultaneous assumptions about the early and late cosmology, with the implication that the physics of both epochs are entwined in the resulting constraints. Models of the late stages of the cosmological evolution rely on less solid physical grounds than the early stages: the physics at decoupling is well understood (atomic physics, general relativity and linearly perturbed FriedmannRobertsonWalker metric) but the late time cosmic acceleration is not, relying, as it does, on a new ingredient such as dark energy.
Most constraints on cosmological parameters are formally modeldependent, with CDM or extensions being routinely assumed. It is therefore natural to ask: “how much do we know about early cosmology?”, “how much of what we know about the early time depends on assumptions about the late cosmology?” And conversely: “how much of what we know about late cosmology depends on assumptions about early cosmology?” These questions are key in the epoch of precision cosmology: once the cosmological parameters of a model are known at the % level, going beyond parameter fitting becomes of interest and it is important to directly test the model itself. One way to achieve this is to perform modelindependent analyses as much as possible.
Ref. [1] proposed that it is possible to analyse the CMB in a manner which is independent of the details of latetime cosmology. Essentially, varying the latetime physics introduces simple changes to the amplitude and angular scaling of the microwave background fluctuations, in the high angular wavenumber range where the integrated SachsWolfe effect is unimportant. This approach, further developed by [2, 3], can be readily exploited, and Refs. [2, 3] propose that it avoids making assumptions about the most relevant latecosmology effects, yielding “pure” or “disentangled” information on early cosmology alone.
Earlytime observables are not the only window into early cosmology. In fact Ref. [4] (see also [5]) provides a measurement of the Baryon Acoustic Oscillation (BAO) scale – the sound horizon at radiation drag, a quantity fixed by early cosmology – from latetime observations. The measurement itself is independent of assumptions of latetime physics, beyond the applicability of the FriedmannRobertsonWalker metric, and makes no assumptions about earlytime physics either.
Here we demonstrate how much we can know about early cosmology in a way that is independent on assumptions about the latetime physics. We use stateoftheart CMB observations to constrain the composition of the early Universe and the properties of the major energy components. Then we explore how constraints on the sound horizon at radiation drag depend on early time physics assumptions. We finally combine these results with constraints on the same sound horizon obtained in a modelindependent way from the latest probe of the Universe expansion history.
This paper is organized as follows. After presenting the methodology and the data sets in §2 we study, analyzing CMB data in the traditional way, an earlydark energy model in §3 to report the latest constraints on the model and to explicitly show the limitations of a modeldependent approach. Pure information about early cosmology is obtained in §4. Finally we present our results and conclusions in §5.
2 Data and Methods
We use the newest Planck Collaboration data release from 2015 [6, 7, 8] which we refer to as “Planck 2015”. We consider low (), temperature and polarization data (referred to as lowTEB) and high () temperature (TT) and polarization (TEEE) data. In the analysis for section 3, we also include the effects on parameter constraints of Planck lensing power spectrum reconstruction (lensing).
Unless otherwise stated, we use the parameter inference code Monte Python [2] interfaced with the Boltzmann code CLASS [9, 10] to generate samples from the posterior via Monte Carlo Markov Chains (MCMC).
We also use the modelindependent measurement of the sound horizon at the end of radiation drag pioneered by Ref. [4] and updated to the latest data by [11]. The sound horizon at the redshift of radiation drag is a standard ruler defined by early Universe physics:
(2.1) 
where is the sound speed and the Hubble parameter. When estimated from CMB observations, is a socalled derived quantity: and are given in terms of other cosmological parameters within a specified cosmological model. Here, for simplicity, we adopt the notation . Ref. [4, 5] showed that the sound horizon, although a property set in the early Universe, could be measured directly from latetime data.
As discussed in [12], can be used as an anchor for the cosmic distance ladder (where the different rungs are supernovae type Ia, BAO etc.), the anchor being set by earlyUniverse observations. Its determination is indirect and therefore dependent on the adopted assumptions about early Universe physics. On the other hand, the cosmic distance ladder can be anchored at using the local determination of . This “direct” ladder, when including BAO and measurements yields an estimate of which does not depend on assumptions about early Universe physics (but relies on the observation of standard candles and a standard ruler).
The results of Ref. [4] have been updated by [11], obtaining Mpc (using JLA type IA SNe data [13], BAO measurements [14, 15, 12] and the Hubble constant determination of Ref.[16]) and Mpc when including, instead of , cosmic chronometer measurements from passive elliptical galaxy ages [17], yielding . These measurements are used to importance sample the Planck 2015 chains and obtain combined constraints.
When performing parameter inference from CMB data in the standard way, simultaneous assumptions about early and latetime physics must be made. As a consequence, in the resulting constraints, early and latetime physics, which should be independent, are inextricably entwined. In §3, we demonstrate this point using a popular model for early dark energy. This motivates us to analyze CMB data in a way that is independent of –or at least robust to detailed assumptions about– latetime cosmology.
3 A practical example: Early Dark Energy
We begin with a modeldependent analysis for Early Dark Energy (EDE), which we make modelindependent in §4. While in standard CDM the contribution of dark energy at early times is fully negligible, dynamical dark energy models might yield a nonnegligible contribution at high redshifts. In particular such models may be similar to CDM at late times (thus being compatible with latetime geometric probes such as SNe) but have significant dark energy density early on and as such they go under the name of early dark energy (EDE) models [18]. EDE slows the growth of structures at early times and changes the heights and positions of the acoustic peaks in the CMB power spectrum. A widely used model is the one proposed by [19], which assumes a constant fraction of early dark energy until a transition at recent times. This model has the advantage of requiring only two extra parameters compared to CDM. It has been used extensively and in particular it was explored in the official analysis of the Planck 2013 [20] and Planck 2015 data [21].
Instead of parameterising the dark energy equation of state, , the authors of [19] proposed a parametrization of the dark energy density as a function of scale factor . Thus the Friedmann equation for a flat universe can be rewritten as
(3.1) 
where is the density parameter of relativistic species (i.e., neutrinos and photons) extrapolated to the present day, and is the presentday matter density parameter (baryonic plus dark matter).
The parametrization is given in terms of the presentday value of the dark energy parameter , the dark energy content in the early Universe and the presenttime dark energy equation of state parameter ,
(3.2) 
The time evolution of the equation of state can be derived from the conservation equation of the dark energy fluid [18]
(3.3) 
were is the scale factor at the equality between radiation and matter. The evolution of both and are represented in Figure 1 for selected sets of parameters. It should be noted that, even if this model is classified as EDE, it has a tracking behaviour. Indeed, during radiation domination the equation of state approaches ; during matter domination it approaches and only at late times it behaves as a cosmological constant if . It is possible to verify that the dark energy density described by (3.2) varies smoothly from to , and the redshift of transition is determined by the model’s parameters , and .
According to [19] this parameterization gives a monotonic function for , as long as
(3.4) 
Contrary to the standard CDM, this model has an evolving dark energy density. This means that it has to be interpreted as a parametrization for a dynamical degree of freedom (e.g. a scalar field) whose scope is to modify the expansion history of the Universe. However, any scalar field has fluctuations, and their contribution in principle should be considered when evolving the perturbations. This has been done in e.g. [21] using the Boltzmann code CAMB [22] and considering a parametrized quintessence fluid at the perturbative level, i.e. a fluid with sound speed and no anisotropic stress. Then it is clear that the differences between this model and a cosmological constant can come from two different sectors, i.e. i) a different expansion history and ii) a pure perturbative effect characterized by the additional degree of freedom. The choice of a quintessence field is arbitrary, but a relativistic sound speed ensures that DE perturbations are suppressed on subhorizon scales and the only effect that potentially modifies the evolution of the perturbations is a different expansion history.
The presence of dark energy perturbations couples earlytime and latetime constraints: perturbations in the dark energy component are most relevant at latetime but their effect is seen in the CMB. To illustrate this point we also consider another approach where we take into account only the contribution given by the different expansion history and not by the perturbed quintessence field using a modification of the public code “CLASS”. The effect can be appreciated in Fig. 2, where we show the temperature and lensing power spectra for different values of . “CAMB” lines are obtained considering perturbations of a quintessence field, while for the “CLASS” lines we modified just the background evolution. Since the effects of a quintessence field could be relevant on very large scales (close to the cosmological horizon), it is not surprising to see that the temperature power spectrum is affected mostly at low , while the inclusion of perturbations modifies the gravitational potential power spectrum uniformly.
The main motivation to consider just a modified expansion history is that the standard equations for the evolution of a minimally coupled scalar field (as quintessence) are singular when the equation of state approaches (and thus a hard prior in has to be imposed). This is a well known property of simple single scalar field models (e.g. [23]), but it is not true for more complicated theories. Then, the “background only” approach (hereafter, CLASS implementation) ensures that we can cross the phantom divide and explore models with .
So while the Planck Collaboration, when dealing with Eq. 3.2, could only obtain upper limits for the relevant parameters, adopting hard priors and [21], we explore, for the first time, the region where and .^{1}^{1}1Of course for models where the density parameter goes negative this component cannot be interpreted as a fluid.
In Fig. 3 we show the marginalized posterior constraints for this early dark energy model accounting only for contributions to the background (CLASS implementation). Different combinations of Planck likelihoods (with or without high polarization and with or without lensing) are considered. In Fig. 4 the results of the analysis including the effects of perturbations (CAMB implementation) is shown. This can be compared directly with the analysis performed by the Planck team.
We note that there are differences in the constraints on this model between the CLASS (background only) implementation and the CAMB (with perturbations) implementation. Besides the fact that the CLASS implementation crosses the phantom divide and constrains values of , , note that, even in the common region (, ) there are differences. In particular in Fig. 4 (with perturbations) there is a correlation and degeneracy between and which is not present in Fig. 3 (backgroundonly implementation). Larger values of early dark energy are allowed when perturbations are included.
Given the fact that so little is known about the nature of dark energy, each of the two approaches should be considered as a phenomenological description of dark energy, each of them capturing different physics. For example, in the nonperturbative case, the effect of dark energy on the background is captured but it is assumed that perturbations in the new degree of freedom have no observable effects. In the perturbative case a very particular perfect fluid, with sound speed equal to 1, is considered.
As shown in Fig. 2, besides the effect on the gravitational potential power spectrum – which affects the CMB lensing signal – differences between the two implementations are evident at low , and arise through the Integrated SachsWolfe (ISW) effect. The ISW effect is a latetime effect, which depends on physics at redshifts much below the last scattering surface, yet it is affected by an earlytime quantity: the amount of early dark energy. The same can be said about the CMB lensing signal, which arises from latetime physical effects (see [24]). The constraints on early dark energy so obtained are therefore highly modeldependent: they depend on assumptions of physics in the latetime Universe.
Parameter  TT  TT_lens  TTTEEE  TTTEEE_lens  Code 


CAMB  
–  CLASS  

CAMB  
–  CLASS  

CAMB  
–  CLASS  
[Mpc] 
CAMB  
–  CLASS  
[km s Mpc] 
CAMB  
–  CLASS 
4 Separating Early Cosmology from Late Cosmology
In the previous section we provided constraints by fitting an EDE model against recent CMB observations. Clearly, the results we presented are model dependent. But, assuming that the model we choose is sufficiently generic, one could naïvely expect that our results can be legitimately interpreted as “early times constraints”. This is not correct, since our EDE model fixes the entire evolution of the universe and the physics of the CMB as we observe it depends on both early and late time cosmology. In other words, the model we chose gives predictions also for the late time expansion history, and the CMB contains this information. It is then clear that, in order to decouple early from late time physics, two steps are needed. The first one is to assume a model to describe cosmology at early times that has to be as general as possible (this point will be described in greater detail in the next section). The second step is to decouple early from late time physics in the CMB spectra.
The key to progress is that latetime physics affects the CMB power spectrum in simple ways at high . [1, 2] show that the effects are limited to: (i) projection effects from real space to harmonic space, controlled by the angular diameter distance at the recombination epoch, i.e. ; (ii) the late ISW effect, affecting only small ; (iii) reionization, suppressing equally all multipoles at . It is possible to show [1, 2] that all these effects produce just a rescaled amplitude () and position () in the CMB high multipoles. Then, to remove the dependence on the latetime cosmology we must not only ignore low multipoles and take into account (and marginalize over) the degeneracy between a direct rescaling of the amplitude and position of the CMB multipoles and late time cosmology. In addition, one has to remove the effects of late times appearing in the lensing signal. It has been argued [1, 2] that this can be done by marginalizing over an overall amplitude and tilt of the Newtonian potential. The resulting approach that we follow is:

we neglect low multipoles for temperature and polarization spectra. We will then consider only “high” Planck data, i.e. ; this effect is illustrated in the top panel of Fig. 5, where it is shown how removing these low multiples suppresses drastically the influence of latetime effects;

we fix the optical depth to a typical value, i.e. , since any value of can be compensated by changing ;

we do not interpret the parameters and as the scalar amplitude of fluctuations and the presentday Hubble parameter. Here effectively sets the lastscattering distance and is only connected with the current expansion rate if we assume a model connecting early to late times. For we interpret only the combination , and for interpret only , which represent vertical and horizontal scale factors of the CMB spectra. In Fig. 6 we show that it is possible to compensate the effects of a variation of and with suitable rescaling and shift of the temperature ;

we marginalize over lensing by rescaling the Newtonian potential as
where (, ) are two new free parameters and we fixed Mpc.
In Fig. 5 we show how the CMB power spectra are modified by choosing different dark energy/modified gravity models with the same earlytime evolution but with different latetime physics. We employ the hi_class [25] code, which implements Horndeski’s theory for dark energy/modified gravity models; Horndeski is the most general scalartensor theory described by secondorder equations of motion, and contains many wellknown dark energy/modified gravity models as special cases (see [27] for more details). It can be seen (upper panels) that the temperature power spectrum is affected mostly by a change in the ISW effect. In the bottom panels we show the polarization power spectrum (EE) and the behaviour of the lensing potential (). As expected, the polarization power spectrum is affected mostly at low multipoles, while the lensing changes at all (even if the effect is larger at low). As it will be clear below, in the simple implementation described in [1, 2] and adopted here, this approach removes the latetime information for models whose latetime evolution is not too far away from a vanilla CDM model. It is not a problem for models with small or vanishing dark energy. However, for models with very large values of the dark energy density parameter (roughly equivalent to at ), the adopted priors matter and the procedure is not guaranteed to remove all the latetime signal. In fact, as we discuss below, in these cases we find residual latetime effects in the form of ISW effects up to . For such models, little can be said about dark energy in a modelindependent way. Constraints on other quantities that are subdominant at late time are indeed modelindependent.
Finally, as explained in [5, 4] low redshift observations can be analyzed so that they yield a cosmologyindependent estimate of , an earlytime quantity. Below we will use this constraint in combination with the above analyses to help constrain earlytime physics.
4.1 Testing Early Cosmology
Here we concentrate on developing a set of models generic enough to test the components of the early Universe and their properties. The idea is to give as much freedom as possible to the earlytime evolution of the Universe with simple and wellmotivated extensions of the standard cosmological model. Then, the models under consideration are

“LCDM”: here we use the standard 6parameter cosmological model with a cosmological constant. Even if this model cannot be considered very general, it is instructive because here cannot be interpreted as the usual cosmological constant that causes the accelerated expansion of the Universe at late times. In principle it should be considered as an early cosmological constant decoupled from the late one;

Dark energy fluid ,“DE fld”: in this case we replace with a standard perfect fluid. We fix the sound speed of the scalar field to 1 and we impose a hard bound on the equation of state parameter, i.e. in order to avoid degeneracy with matter. The freedom here is represented by the fact that the additional scalar field is not constrained to have . While is a free parameter (hereafter ) it is considered constant in time;

“k”: here we add to the cosmological constant an arbitrary spatial curvature with density parameter . The idea of having a spatial curvature is often used in the literature, and it is very much constrained to be close to 0 by latetime observations (e.g. [8]. However, here we focus just on the earlytime evolution and thus expect poor constraints on this parameter. It is useful to introduce it since the additional freedom given by a parameter that scales differently from a cosmological constant is able to modify substantially the earlytime expansion history of the Universe;

“nu”: This is effectively a standard CDM model but with the addition of nonstandard neutrinos. The idea [28] is to modify the properties of neutrinos by varying two parameters, an effective sound speed and a viscosity parameter . This case is reduced to the standard one by fixing both values to 1/3. In this model we fix one massive neutrino with mass but we let vary the number of massless neutrinos, i.e. for the standard case. The massive and massless neutrinos share the same phenomenological properties, i.e. and are in common. It is possible to show that nonstandard neutrinos behave as a standard scalar field that scales as radiation at the background level by setting . Then, in this model the contribution of an additional scalar field at the perturbation level can come directly from the neutrino sector at the price of renouncing the standard neutrino behaviour;

“matteriation”: This is similar to the “DE fld” model, but with an equation of state parameter that can take values intermediate between that of pressureless nonrelativistic matter and that of radiation . For simplicity here we assume that this fluid contributes to a modification of the expansion history only, while the perturbations are the standard matter/metric perturbations. This case complements the “nu” case above.
Figures 7, 8, 9, 10 and 11 show the 68% and 95% joint posterior credibility regions for the relevant parameters for the “LCDM”, “DE fld”, “k”,“nu” and “matteriation” cases respectively using Planck 2015 data (TT and TEEE). Since we wish to obtain constraints on earlytime cosmology only, the density parameters in the plots are reported at recombination (). In all cases the densities of nonstandard species are constrained to be very small compared to the matter density (). Special care must be made when interpreting the constraints on especially for the curvature case. and are degenerate but the degeneracy is not perfect. Positive values of are bounded from above due to the prior that (as well as all other densities) must be positive and a minimal amount of is needed to produce the acoustic peaks. However the degeneracy does not appear to extend to positive values beyond and does not extend beyond , which essentially restricts the presentday dark energy density parameter to be no more than . Larger values may introduce ISW fluctuations on scales where our analysis assumes there are none, so we caution against overinterpreting this limit. Constraints on this parameter are expected to be of the order of or below of the matter density at recombination [29]. The unphysical constraint on is due to residual latetime constraints. We note that in the nonflat case the constraint on is much weaker () than in the other cases. This does not invalidate the derived constraints on the other parameters, as and are in any case subdominant at recombination. Our findings confirm that the Universe at early times is well described by the standard cosmological model, and there is no much room for additional components that may be relevant at early times.
What the observations constrain remarkably well about the early cosmology is the expansion rate at recombination (especially for standard radiation and neutrinos contributions). This can be appreciated in Tab. 2 (third line, entry). As described in [30] the faster the Universe is expanding at recombination, the more difficult it is for the hydrogen to recombine, increasing the ionization fraction. Larger ionization fraction yields a broader visibility function. This in turn has two effects on observable quantities. On the CMB temperature it leads to a larger (Silk) damping of small scales anisotropies: the first acoustic peak is unaffected but higher peaks are damped. On the polarization, because of the increased photon mean path, largescale polarization anisotropies () are enhanced. Planck was forecasted to constrain at the 0.9% level [30] and Tab. 2 confirms it. Adjustments in the composition of the Universe around recombination time cannot therefore drive changes larger than that.
In summary, the expansion rate from matterradiation equality through recombination and especially at recombination is very well constrained. This tight constraint on leads in turn to constraints, at the same level, on the physical densities. Moreover, the amount of radiation, neutrinos and (physical) cold dark matter densities are constrained also by the perturbations. As a result, our analysis shows the residual (small) freedom in the early Universe composition.
The results for the “LCDM” case are reported in more detail in Appendix A. These results may be of interest beyond the scope of this paper. In fact, because they are robust to detailed assumptions about dark energy properties they can be used as a darkenergy “hardened” CMB prior for dark energy analyses.
As discussed above, one caveat in interpreting these results is that we have excluded multipoles to remove the latetime ISW effect. In principle, for nonstandard models, there could still be some latetime contribution at which, in this approach, could then mistakenly be interpreted as early cosmology signal. This affects the constraints on dark energy but not on the other parameters.
In Fig. 12 we show how latetime effects that become important when at is cannot be fully removed by the approach of [1, 2]. For example, the late ISW effect extends up to , modifying the shape of the first peak and making the rescaling imperfect, when are included. Excluding all would yield greatly degraded constraints.
In all the plots we show also the posterior distribution of the soundhorizon at radiation drag. This is an earlytime quantity which with this approach is now measured in a way that is independent of assumptions about latetime physics. The sound horizon at radiation drag, , remains well constrained for all the models under consideration, with a remarkable exception. The only one that allows for a large variation of is the “nu” model. A possible explanation is given by the fact that this is the only model that modifies directly the properties of species (and thus the expansion rate) that are “naturally” important at early times (neutrinos). In the other examples, we added species that have a density that scales more rapidly than radiation in the past, and thus become progressively unimportant at early times (e.g., ).
The precision in this measurement is better than 0.5%, which is just a factor two larger than for the standard analysis [8] being the central value fully consistent. It is interesting to note that the modelindependent measurement of the “local standard ruler” [11] is in good agreement with the value of found from the earlytime only analysis and has errorbars comparable to those on found in the “nu” model. As shown in Fig. 10, including this additional constraint slightly alters the bounds on this model, in particular it slightly reduces the error bars on bringing the 95% credible interval to . Forecasted improvements on Baryon Acoustic Oscillation data and on measurements of are expected to further reduce the error on the ruler (see e.g., [11] for a discussion), promising therefore more stringent tests on the physics in the early Universe.
LCDM  DE fld  k  nu  matteriation  
[Mpc]  
        
      
      
        
        
        
5 Conclusions
We have explored in a modelindependent way, how well the early universe conditions (i.e. universe components and their properties) are constrained by current CMB observations. To do so, we have used a method that attempts to decouple the lateuniverse from the earlyuniverse. We have compared the results with modelindependent determinations of the standard ruler corresponding to the sound horizon at radiation drag from the literature, finding good agreement.
In our method, we allow for the possibility that energydensity components different from those in the CDM model are present in the early Universe. In particular, we allowed for extra dark energy (“DE fld”, with equation of state ), extra fluid matter (“MR”, ), extra radiation or other relativistic species (, parameterized by ) and extra curvature (“k”), covering many possibilities for the expansion history of the universe at early times. Our methodology, as presented and implemented here works very well for models that at late time are not too drastically different from the standard CDM model. In particular, there are residual late time effects when at . For dark energy constraints, the methodology as currently implemented, should be interpreted as “robust” to detailed assumptions about dark energy rather than fully modelindependent.
Stateoftheart CMB observations including temperature and polarization information, effectively constrain the expansion rate at recombination with precision. The expansion rate governs the ionization fraction and the width of the visibility function, which in turn affects the Silk damping and the amplitude of the polarization signal [30]. Adjustments in the composition of the Universe around recombination time are constrained by this effect. Moreover, components that cluster are constrained also through their effects on the perturbations.
We found that current observations constrain surprisingly tightly these extra components: and when imposing spatial flatness. These energy densities are all reported as 95% credible intervals, at . For the latter case of extra radiation, when we use the local measurement of the standard ruler to further limit the amount of dark radiation, we obtain that using SNe, BAOs and a Hubble prior (SBH), and replacing the Hubble prior with cosmic chronometers (CSB).
Our conclusion is that current CMB (temperature and polarization) observations, alone or in conjunction with lowredshift, modelindependent measurements of the standard ruler, tell us that the earlyuniverse (up to recombination) is very well known: the data require the presence of baryons, radiation, a fluid that behaves like neutrinos and dark matter; dark energy and curvature are negligible. There is no evidence –and not much room– for exotic matterenergy components besides those in the CDM; the standard CDM model describes extremely well the early cosmology.
We envision that modelindependent approaches like the one presented here will be of value as cosmology moves beyond parameter fitting and towards model testing.
Acknowledgements
Funding for this work was partially provided by the Spanish MINECO under projects AYA201458747P and MDM2014 0369 of ICCUB (Unidad de Excelencia Maria de Maeztu), Royal Society grant IE140357 and by CNPq (Brazil). We acknowledge the support of Imperial College through the CosmoCLASSIC collaboration.
Appendix A
Largescale structure, supernovae or other low redshift observations analyses on dark energy properties are often slowed down significantly by the fact that, to include CMB information, a Boltzmann code must be suitably modified and ran for each specific model. Not including CMB information simplifies the analysis greatly, but also reduces dramatically the constraints.
One approach adopted so far has been to use the socalled CMB distance priors ([31, 32] and [33] and references therein).
Here we propose a slightly more sophisticated approach of using the results provided by the (latetime model independent) “LCDM” model analysis and reported in Tab. 2. These results are by construction independent of latetime assumptions about perturbations and are “robust” to detailed assumptions about the expansion history and in particular about the detailed dark energy properties, yet they fully capture the CMB information assuming standard earlytime physics. The values for the density parameters are reported at the redshift of recombination. They can then be easily extrapolated to lower redshifts for any given model for the Universe expansion history. The recombination redshift is well constrained and does not show significant degeneracy with any of the parameters, so for this purpose it can be assumed to be fixed at .
References
 [1] Marc Vonlanthen, Syksy Rasanen, and Ruth Durrer, Modelindependent cosmological constraints from the CMB, JCAP, 1008:023, 2010.
 [2] Benjamin Audren, Julien Lesgourgues, Karim Benabed, and Simon Prunet, Conservative Constraints on Early Cosmology: an illustration of the Monte Python cosmological parameter inference code, JCAP, 1302:001, 2013.
 [3] Benjamin Audren. Separate Constraints on Early and Late Cosmology, Mon. Not. Roy. Astron. Soc., 444(1):827–832, 2014,
 [4] Alan Heavens, Raul Jimenez, and Licia Verde, Standard rulers, candles, and clocks from the lowredshift Universe, Phys. Rev. Lett., 113(24):241302, 2014.
 [5] Will Sutherland, On measuring the absolute scale of baryon acoustic oscillations, Mon. Not. Roy. Astron. Soc., 426:1280, 2012.
 [6] R. Adam et al. [Planck Collaboration], Planck 2015 results. I. Overview of products and scientific results, Astron. Astrophys. 594 (2016) A1 doi:10.1051/00046361/201527101 [arXiv:1502.01582 [astroph.CO]].
 [7] N. Aghanim et al. [Planck Collaboration], Planck 2015 results. XI. CMB power spectra, likelihoods, and robustness of parameters, Astron. Astrophys. 594 (2016) A11 doi:10.1051/00046361/201526926 [arXiv:1507.02704 [astroph.CO]].
 [8] P. A. R. Ade et al. [Planck Collaboration], Planck 2015 results. XIII. Cosmological parameters, Astron. Astrophys. 594 (2016) A13 doi:10.1051/00046361/201525830 [arXiv:1502.01589 [astroph.CO]].
 [9] J. Lesgourgues, “The Cosmic Linear Anisotropy Solving System (CLASS) I: Overview,” arXiv:1104.2932 [astroph.IM].
 [10] Diego Blas, Julien Lesgourgues, and Thomas Tram, The Cosmic Linear Anisotropy Solving System (CLASS) II: Approximation schemes, JCAP, 1107:034, 2011.
 [11] L. Verde, J. L. Bernal, A. F. Heavens and R. Jimenez, The length of the lowredshift standard ruler, arXiv:1607.05297 [astroph.CO].
 [12] Antonio J. Cuesta et al, The clustering of galaxies in the SDSSIII Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the correlation function of LOWZ and CMASS galaxies in Data Release 12, Mon. Not. Roy. Astron. Soc., 457:1770, 2016.
 [13] M. Betoule et al, Improved cosmological constraints from a joint analysis of the SDSSII and SNLS supernova samples, Astron. Astrophys., 568:A22, 2014.
 [14] Eyal A. Kazin et al, The WiggleZ Dark Energy Survey: improved distance measurements to z = 1 with reconstruction of the baryonic acoustic feature, Mon. Not. Roy. Astron. Soc., 441(4):3524–3542, 2014.
 [15] Florian Beutler, Chris Blake, Matthew Colless, D. Heath Jones, Lister StaveleySmith, Lachlan Campbell, Quentin Parker, Will Saunders, and Fred Watson. The 6dF Galaxy Survey: Baryon Acoustic Oscillations and the Local Hubble Constant. Mon. Not. Roy. Astron. Soc., 416:3017–3032, 2011.
 [16] A. G. Riess et al., A 2.4% Determination of the Local Value of the Hubble Constant, Astrophys. J. 826 (2016) no.1, 56 doi:10.3847/0004637X/826/1/56 [arXiv:1604.01424 [astroph.CO]].
 [17] Michele Moresco, Lucia Pozzetti, Andrea Cimatti, Raul Jimenez, Claudia Maraston, Licia Verde, Daniel Thomas, Annalisa Citro, Rita Tojeiro, and David Wilkinson, A 6% measurement of the Hubble parameter at : direct evidence of the epoch of cosmic reacceleration, 2016.
 [18] Christof Wetterich, Phenomenological parameterization of quintessence, Phys. Lett., B594:17–22, 2004.
 [19] Michael Doran and Georg Robbers, Early dark energy cosmologies, JCAP, 0606:026, 2006.
 [20] P. A. R. Ade et al. [Planck Collaboration], Planck 2013 results. I. Overview of products and scientific results, Astron. Astrophys. 571 (2014) A1 doi:10.1051/00046361/201321529 [arXiv:1303.5062 [astroph.CO]].
 [21] P. A. R. Ade et al. [Planck Collaboration], Planck 2015 results. XIV. Dark energy and modified gravity, Astron. Astrophys. 594 (2016) A14 doi:10.1051/00046361/201525814 [arXiv:1502.01590 [astroph.CO]].
 [22] Antony Lewis, Anthony Challinor, and Anthony Lasenby, Efficient computation of CMB anisotropies in closed FRW models, Astrophys. J., 538:473–476, 2000.
 [23] Guillermo Ballesteros and Julien Lesgourgues. Dark energy with nonadiabatic sound speed: initial conditions and detectability, JCAP, 1010:014, 2010.
 [24] A. Kosowsky, M. Milosavljevic and R. Jimenez, Efficient cosmological parameter estimation from microwave background anisotropies, Phys. Rev. D 66 (2002) 063007 doi:10.1103/PhysRevD.66.063007 [astroph/0206014].
 [25] M. ZumalacÂrregui, E. Bellini, I. Sawicki and J. Lesgourgues, hiclass: Horndeski in the Cosmic Linear Anisotropy Solving System, arXiv:1605.06102 [astroph.CO].
 [26] Emilio Bellini and Ignacy Sawicki, Maximal freedom at minimum cost: linear largescale structure in general modifications of gravity, JCAP, 1407:050, 2014.
 [27] E. Bellini and I. Sawicki, Maximal freedom at minimum cost: linear largescale structure in general modifications of gravity, JCAP, 7:050, July 2014.
 [28] Wayne Hu, Daniel J. Eisenstein, Max Tegmark, and Martin J. White, Observationally determining the properties of dark matter, Phys. Rev., D59:023512, 1999.
 [29] T. Karwal and M. Kamionkowski, Early dark energy, the Hubbleparameter tension, and the string axiverse, arXiv:1608.01309 [astroph.CO].
 [30] O. Zahn and M. Zaldarriaga, Probing the Friedmann equation during recombination with future cosmic microwave background experiments, Phys. Rev. Lett, 67(6):063002, 2003.
 [31] Yun Wang and Pia Mukherjee, Observational Constraints on Dark Energy and Cosmic Curvature, Phys. Rev., D76:103533, 2007.
 [32] Hong Li, JunQing Xia, GongBo Zhao, ZuHui Fan, and Xinmin Zhang, On using the WMAP distance priors in constraining the time evolving equation of state of dark energy, Astrophys. J., 683:L1–L4, 2008.
 [33] QingGuo Huang, Ke Wang, and Sai Wang, Distance priors from Planck 2015 data, JCAP, 1512(12):022, 2015.