Cosmic tango between the very small and the very large: Addressing CMB anomalies through Loop Quantum Cosmology

While the standard, six-parameter, spatially flat $\Lambda$CDM model has been highly successful, certain anomalies in the cosmic microwave background bring out a tension between this model and observations. The statistical significance of any one anomaly is small. However, taken together, the presence of two or more of them imply that according to standard inflationary theories we live in quite an exceptional universe. We revisit the analysis of the PLANCK collaboration using loop quantum cosmology, where an unforeseen interplay between the ultraviolet and the infrared makes the \emph{primordial} power spectrum scale dependent at very small $k$. Consequently, we are led to a somewhat different $\Lambda$CDM universe in which anomalies associated with large scale power suppression and the lensing amplitude are both alleviated. The analysis also leads to new predictions for future observations. This article is addressed both to cosmology and LQG communities, and we have attempted to make it self-contained.


I. INTRODUCTION
The quantum geometry effects underlying loop quantum gravity (LQG) lead to a natural resolution of the big bang singularity (see, e.g., [1,2] for reviews). Therefore, one can hope to meaningfully extend the standard inflationary paradigm to the Planck regime. Over the past decade, several closely related approaches have been used to carry out this task, leading to a striking interplay between theory and observations (see, in particular, [3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21]). In this article we will focus on the recent results that shed new light on the anomalous features seen in the cosmic microwave background (CMB). Specifically, we will show that in our approach two of the anomalies seen in the CMB can be accounted for using the pre-inflationary dynamics of loop quantum cosmology (LQC). This phase of dynamics alters the quantum state of cosmological perturbations at the onset of the (relevant part of the) slow roll, leading to revised values of the six parameters that characterize the ΛCDM universe. The revision alleviates the tension due to two anomalies that have received considerable attention, while leaving the successes of standard inflation intact. Main results were reported in [22]. The purpose of this paper is to provide details and also present some supplementary material to put the results in a broader context. These results illustrate that LQC has matured sufficiently to lead to testable predictions.
The paper is addressed both to the LQG community and cosmologists. For the benefit of the LQG community, that primarily focuses on mathematical physics, we have included a discussion of the interplay between theory and observations that leads to the six parameter ΛCDM cosmological model. We will summarize the underlying procedure and point out certain subtleties in data analysis. For cosmologists, we will summarize the key features of LQC that lead to new observable predictions. Specifically we will explain how the quantum geometry effects in the ultraviolet, that lead to the singularity resolution, have unforeseen and interesting consequences on the dynamics of cosmological perturbations in the infrared. It is this 'cosmic tango' between the very small and the very large that alleviates anomalies. Overall, in terms of conceptual flow, we have attempted to make this paper self-contained. In particular, within the page limits of this special issue, we clarify apparently conflicting statements in the LQC literature. In order to make the material accessible to both communities, we will have to briefly review ideas and results that are likely to be well-known in one community but not the other.
The two anomalies we focus on arise as follows. Motivated in large part by inflationary scenarios, the CMB analysis generally begins by assuming that the primordial scalar power spectrum has a nearly scale invariant form, characterized by just two parameters, the scalar perturbation amplitude A s and the scalar spectral index n s . We will refer to this form as the standard ansatz (SA). A s , n s and 4 other parameters (discussed in section II A) characterize a specific ΛCDM universe. Given these six parameters one can evolve the primordial perturbations using known astrophysics and predict the observable power spectra. By varying the values of the 6 parameters, and confronting the theoretical prediction with observations, one finds the posterior probability distributions of the six cosmological parameters. By and large the CMB observations can be well explained using the ΛCDM universe determined by the marginalized mean values of these parameters. However, one also finds some anomalous features. The first is power suppression at large angular scales: the observed power in the temperature-temperature (TT) spectrum is suppressed for 30 in the spherical harmonic decomposition, relative to what the theory predicts. The second anomaly we focus on is associated with the so-called lensing amplitude, A L , associated with gravitational lensing that the CMB photons experience as they propagate from the surface of last scattering to us. The ΛCDM cosmology based on the SA assumes A L = 1 while, when it is allowed to vary, A L prefers a value larger than unity. This tension hints at an internal inconsistency. To alleviate it, one can introduce a positive spatial curvature but then there are inconsistencies with the low z measurements, prompting a recent suggestion [23] that this anomaly gives rise to a "possible crisis in cosmology." As we will see, both the anomalies are simultaneously alleviated in our approach. The key new element is the following: Pre-inflationary dynamics of LQC leads to a primordial power spectrum that differs from SA, but only at large angular scales. While it continues to be nearly scale invariant -and essentially indistinguishable from the one given by the SAfor 30, there is a specific power suppression for 30. As a result, the best-fit values of the six cosmological parameters change. Interestingly, the change in 5 of the 6 parameters is extremely small, 0.4%. But the value of the 6th parameter -the optical depth τ -is increased by ∼ 9.8%! We will see that this change then leads to the alleviation of the tension between observations and theoretical predictions based on the SA. Note that in spite of this significant change in the value of τ , LQC leaves the highly successful predictions of standard inflation at small angular scales unaffected. In particular, all the finer features of various power spectra predicted by standard inflation for > 30 -where the observational error bars are small-are present also in the LQC prediction. Thus the LQC analysis provides an explicit example supporting a conclusion of [24] that trans-Planckian effects are not a "threat to inflation".
The paper is organized as follows. Section II summarizes the procedure used in observational cosmology to arrive at the 6 parameter ΛCDM model and explains the two anomalies and their significance in greater detail. Section III summarizes the basic results from LQC that are used in the subsequent analysis. In particular, we explain the origin of the surprising interplay between the ultraviolet and the infrared that is a rather robust feature of the LQC approaches. The main results are presented in Section IV. They include a discussion of: (i) the LQC corrected primordial power spectrum for scalar perturbations; (ii) the TT, temperature-electric polarization (TE), the electric polarization (EE), and the lensing potential (φφ) power spectra we predict, and comparisons with those obtained using the SA as well as with the observed power spectra reported by the PLANCK team in their final analysis [25]. As usual these power spectra are presented in terms of the spherical harmonic components C of the respective correlation functions; (iii) the TT correlation function C(θ) predicted by LQC and its comparison with the prediction of the SA as well as PLANCK observations; (iv) the A L versus τ plots that show that the observed values fall in the 1σ contour in LQC, but not if one uses the SA; and, (v) the power spectrum for BB polarization predicted by LQC and its comparison with that predicted by the SA. The detailed calculations underlying these plots were performed using the Starobinsky and quadratic potentials. The first is preferred phenomenologically while the second has been used often because of its simplicity. Close agreement between the two sets of results is an indication of robustness of the LQC results. In Section V we summarize the main results and put them in a broader context. So far the discrepancy between the results of the SHOES team [26] and CMB measurements [27] associated with the value of the Hubble parameter has not been systematically addressed in LQC. This is in large part because it is not yet clear whether there is a definitive tension, or if the discrepancy is primarily due to systematic calibration offsets [28]. Observations may decide on this issue in the near future.

II. THEORETICAL PREDICTIONS AND PLANCK OBSERVATIONS
This section is addressed primarily to the LQG community. In II A we summarize the procedure used by the PLANCK team to arrive at the six parameter ΛCDM model and in II B we explain the power suppression and the lensing amplitude anomalies in a bit more detail.
A. The 6 parameter ΛCDM model The six parameters that characterize the ΛCDM universe can be neatly divided in three groups. The first two parameters -the amplitude A s and the spectral index n s for scalar modes-feature in the SA for the primordial power spectrum: Here k is the wave number in the Fourier decomposition and k is a pivot scale (set to k = 0.002Mpc −1 in the WMAP analysis and k = 0.05Mpc −1 in the PLANCK analysis). If we had n s = 1, the primordial spectrum would be scale invariant, i.e., it would be independent of the wave number k of the cosmological perturbation. If n s is less than 1 (as observations imply) then there is more power at small k, i.e., the power spectrum has a red tilt. (One can also consider the possibility of a running n s , where it has a k dependence but we will not need this generality.) The second set of parameters, the baryonic and cold matter densities Ω b h 2 and Ω c h 2 , are important for the propagation of cosmological perturbations starting from the end of inflation. The last group of parameters are 100θ M C , that characterizes the angular scale of acoustic oscillations, and the optical depth τ that characterizes the reionization epoch. These two parameters govern the propagation of perturbations from the last scattering surface to now. Thus, given these six parameters one can use the known astrophysics to propagate the cosmological perturbations starting from the end of inflation, providing us with the theoretical predictions for power spectra we should observe now. More precisely, for each choice of the six parameters, one can calculate 4 correlation functions C TT , C TE , C EE , C φφ that are the spherical harmonic decompositions of the corresponding correlation function C XY (θ) (where θ characterizes the angular separation of two points in the sky). These correlation functions can be measured and compared with the theoretical predictions, expressed as functions of the six parameters. The statistical analysis, usually done by employing Markov-Chain Monte Carlo method, then leads to the posterior probability distribution for the six parameters. In particular, the maximum likelihood value of the marginalized probability distributions yield the values of six parameters that determine a ΛCDM universe. Observations of the PLANCK collaboration provide these values (together with the corresponding 1-sigma spreads); this is the "universe according to PLANCK" (within the 68% confidence level, characterized by the 1σ contours).
Once these parameters are determined, one can calculate additional observable quantities assuming that model and, by carrying out measurements, one can subject the model to consistency tests. For example, the lensing amplitude A L is set to unity in this construction. One can let this parameter vary and test if this value is consistent with observations. Another type of test is provided by the (odd parity) B-modes. In any one model, one can calculate the correlation function C BB . As we discuss in Section V, several observational missions will soon measure this correlation function with accuracy that may be sufficient to distinguish one model from another [29,30,32]. Similarly, the reionization depth τ will be measured by missions that are unrelated to the CMB [33]. They will constrain τ , providing us with independent checks on the current ΛCDM model.

B. Anomalies
As Fig. 1 shows, the TT power spectrum is in excellent agreement with the theoretical predictions using the SA at small angular scales ( > 50). This is especially noteworthy because the instrumental errors are truly minuscule in this range.
However, for 30 the observed power is lower than the theoretical prediction. This power suppression was evident already in the WMAP data, and is reenforced by the PLANCK findings. Over the years it has been argued [34][35][36][37] that this anomaly is brought to forefront if one carries out the comparison using the quantity S 1/2 := harmonic coefficients C . Qualitatively, large angular scales correspond to small in the spherical harmonic decomposition. However, for any given θ o the value C(θ o ) receives contributions from all . Therefore S 1/2 is a more direct measure of the cumulative power for θ ≥ 60 • than the C 's for low s. Indeed, as the left panel in Fig 2 vividly shows, the observed C(θ) (the black, dotted curve) is very close to zero for θ > 60 • , in contrast to the theoretical (solid, red) curve. The observed value of S 1/2 is ∼ 1209, while the theoretical prediction from the SA is 42, 496.5, some 35 times larger. Since the extent of this discrepancy is not immediate from Fig. 1, one may wonder why the power spectrum are not reported using C TT (θ) in place of C TT . The reason is that the C TT for distinct s are (almost) uncorrelated and can therefore be treated as independent observables, while C TT (θ) for distinct θ have massive cross-covariance, whence the statistical significance of power suppression is only 2−3σ in spite of the large deviation seen in the left panel of Fig. 2. Also, because of these correlations, to obtain the 1 and 2σ contours in this plot one has to take into account a large covariance matrix which in turn requires a detailed understanding of the instruments and the masking procedure near θ = 180 • used in the data analysis to remove the contamination coming from the galactic plane. The second anomaly is associated with the lensing amplitude A L depicted in the right panel of Fig. 2. As it propagates from the last scattering surface at z ≈ 1100 to us, the CMB is lensed due to inhomogeneities. The lensing potential is nearly Gaussian because there are many lenses along the line of sight. As explained in Section II A, the six parameter ΛCDM universe is determined using best fits to all four power spectra. Once this is done, one can compare each observed power spectrum, one by one, with the theoretically predicted power spectrum for that specific ΛCDM universe. Just as this comparison revealed an anomalous suppression of power in C TT for < 30, one finds an anomaly also in the lensing potential power spectrum C φφ : Relative to the prediction of the best-fit ΛCDM model, there is power enhancement in the range 8 ≤ ≤ 400 used by the PLANCK collaboration to report the baseline cosmological results [27]. (In this range, the reconstruction procedure is robust and the impact of systematics is reduced). As a consistency check on the 6-parameter ΛCDM model, one introduces a 7th phenomenological parameter A L -the lensing amplitude, normalized so that A L = 1 in the 6 parameter ΛCDM model-and allows it to vary. Varying A L can be considered as a conservative way of marginalizing over the systematics of the PLANCK data. Departure of A L from unity signals a tension with predictions based on the standard 6-parameter ΛCDM . One finds that A L is higher than 1 at ∼ 1.9σ level [27]. This is the lensing amplitude anomaly and its occurrence has been interpreted as a hint of new physics [38]. The right panel of Fig. 2 illustrates this tension. Here A L is plotted against the optical depth τ because this plot will be useful when we compare the results from LQC with those from the SA in section IV: Of the six parameters, only τ receives significant corrections from LQC. τ is singled out even within the SA by the fact that the relative error (as measured by the ratio of the standard deviation to the mean value) in τ is ∼ 13% while that in the other five ΛCDM parameters are less than 1%. In the plot, the tension is manifested in the fact that the line A L = 1 lies outside the 1σ contour. Attempts to alleviate this tension within the standard paradigm based on general relativity (GR) -e.g. changing the background geometry by introducing spatial curvature-are not supported by lensing reconstruction or Baryonic oscillations (BAO) data (since the joint constraint with BAO is consistent with flat universe, with Ω K = 0.001 ± 0.002).
As noted in the introduction, while the statistical significance of either of these anomalies is low, together the two imply that the observed universe will emerge only once in ∼ 10 6 realizations of the posterior probability distributions. Therefore, as the PLANCK collaboration suggested both in its 2015 and 2018 data releases, alleviation of this tension is of considerable interest especially if the mechanism is rooted in physics beyond GR [39,40].

III. LOOP QUANTUM COSMOLOGY
This section is addressed primarily to the cosmology community. In III A we briefly recall how quantum geometry effects underlying LQG lead to a resolution of the big bang singularity, replacing it with a big bounce. Since physical quantities do not diverge anywhere, one can extend the standard inflationary scenario all the way to the bounce. In III B, we discuss the pre-inflationary dynamics of cosmological perturbations, specifically the propagation of quantum fields representing these perturbations on the quantum background geometry provided by LQC. In III C, we explain why, contrary to one's initial expectations, this pre-inflationary dynamics can leave observable signatures at large scales in the CMB. We will use this framework in section IV to extract the LQC corrections to the primordial power spectrum. As mentioned in section I, together with observations, these corrections imply that we live in a somewhat different ΛCDM universe in which the two anomalies are naturally alleviated.
A. The big bounce of LQC Investigations of the early universe are often carried out assuming that spacetime is well approximated by a spatially flat FLRW background metric of GR, together with first order cosmological perturbations that are described by quantum fields. Consider the inflationary paradigm and, for brevity, let us refer to the time when the pivot mode k = k exits the Hubble horizon simply as 'the onset of inflation'. At this onset, while spacetime curvature is huge by astrophysical standards -some 10 65 times that at the horizon of a solar mass black hole-it is only ∼ 10 −12 times the Planck scale. Therefore, at the level of accuracy of current interest, it is safe to ignore the quantum gravity effects even at the onset of inflation. Since spacetime geometry is well approximated by a (perturbed) de Sitter metric at this time, one assumes that the quantum fields representing cosmological perturbations are in the Bunch-Davies (BD) vacuum that is selected by the isometries of the de Sitter metric and evolves the perturbations to the future (as curvature decreases further).
However, conceptually it is rather ad-hoc to begin, so to say, 'in the middle' of evolution. If we go further back in the past, curvature attains the Planck scale, and then diverges at the big bang. During this pre-inflationary epoch, spacetime geometry is not at all wellapproximated by the de Sitter geometry. Why, then, can we assume the state to be the BD at the onset of inflation? Should we not start in the deep Planck regime and check whether the state is in fact in the BD vacuum at this onset? This would require quantum cosmology, where the Friedmann, Lemaître, Robertson, Walker (FLRW) solution of Einstein's equations, characterized by the scale factor a(t) and a matter field φ(t), is replaced by a quantum state Ψ(a, φ) subject to an appropriate quantum version of Einstein-matter field equations. Note that reference to the proper time t has disappeared -quantum dynamics is relational,à la Leibnitz: for example, one can use the matter field φ as an internal clock, and describe how the scale factor evolves with respect to it. Quantum fields representing cosmological perturbations are now to propagate on a quantum FLRW geometry Ψ(a, φ) which assigns probability amplitudes to various metrics, rather than on a single FLRW spacetime.
While this general viewpoint is common to all quantum cosmologies, LQC has two key features that distinguish it from the older Wheeler-DeWitt (WDW) theory, often called quantum geometrodynamics. First, as explained below, the mathematical framework of LQC descends from the well-developed kinematics of LQG, using a symmetry reduction tailored to homogeneity and isotropy. In the WDW theory one is yet to develop rigorous kinematics for full quantum geometrodynamics; because issues related to the presence of an infinite number of degrees of freedom are generally ignored, the underlying mathematical framework has remained formal. In quantum cosmology, then, one introduces structures like the WDW equation without guidance from a more complete framework. This leads to the second key difference. The LQC quantum Einstein's equation is qualitatively different from the WDW equation, in that it mirrors features of the quantum Riemannian geometry of full LQG. As a direct result, strong cosmological singularities -and in particular the big bangare naturally resolved in LQG [2,41,42].
We will now explain these differences in some detail. As is common in quantum field theories, in full LQG one begins with the Heisenberg algebra A of basic ('canonically conjugate') observables, called the holonomy-flux algebra [43][44][45]. We then have a highly non-trivial result that ensures that A admits a unique representation by operators of a Hilbert space H that respects the 'background independence' of 'diffeomorphism covariance' of the theory [46,47]. This representation underlies the rigorous kinematical framework of LQG. In particular, one finds that geometrical observables are well-defined self-adjoint operators with discrete eigenvalues. Of particular interest is the area gap -the first non-zero eigenvalue ∆ of the area operator. It is a fundamental microscopic parameter of the theory that then governs important macroscopic phenomena in LQC that lead, e.g., to finite upper bounds for curvature. 1 In LQC, one first reduces the holonomy-flux algebra A used in full LQG to a smaller symmetry reduced algebra A red . Again there is a uniqueness theorem that guarantees that A red admits a unique representation on a Hilbert space H LQC that respects the action of the (residual) diffeomorphism group on A red [48,49]. This representation is qualitatively different (i.e. unitarily inequivalent) from the Schrödinger representation used in the WDW theory. In particular, the differential operator representing the gravitational part in the WDW equation is not even defined on H LQC ; it is naturally replaced by a certain difference operator that explicitly involves the area gap ∆ [42,50,51]. One can now start with a quantum state Ψ(a, φ) that is peaked on the classical dynamical trajectory at a suitably late time when curvature is low, and evolve it back in time towards the big bang using either the WDW equation or the LQC evolution equation. Interestingly the wave function continues to remain sharply peaked in both cases. In the WDW theory it follows the classical trajectory all the way into the singularity, while in LQC it ceases to follow the classical trajectory once the curvature is about ∼ 10 −4 times the Planck curvature. Then the quantum geometry corrections dominate and the wave function Ψ(φ, a) bounces when the curvature and matter density attain their upper bounds. In this backward evolution, curvature starts decreasing after the bounce and the universe expands. Once the curvature falls below ∼ 10 −4 times the Planck curvature, the wave function again follows a classical trajectory which is now expanding in the past direction. (For details, see [2,42,51,52].) Thus, key differences between LQC and the WDW theory arise from the fact that the WDW theory has no knowledge of the quantum nature of Riemannian geometry that LQC inherits from LQG. Indeed, there is a precise sense in which the LQC evolution equation reduces to the WDW differential equation in the limit in which the area gap goes to zero [52]. The upper bound of curvature in LQC is given by curv any LQC solution Ψ lqc (a, φ), the curvature attains its maximum value at the bounce and this value is extremely well approximated by curv max if the state is sharply peaked. Note that the upper bound diverges as ∆ → 0, in line with the finding that curvature grows unboundedly as one evolves the WDW state Ψ wdw (a, φ) to the past. By contrast, in LQC, while the quantum geometry effects are negligible away from the Planck regime, they become dominant in the Planck regime, creating an effective repulsive force of quantum origin that causes the universe to bounce.
It is interesting that this force rises and falls extremely rapidly, making the agreement with GR excellent outside the Planck regime. However, it has a very non-trivial global effect, in that physics does not stop at the big bang as in GR. Rather, there is an expanding FLRW universe to the future of the bounce and a contracting FLRW universe to the past, joined by a 'quantum bridge'. These qualitatively new features arise without having to introduce matter that violates any of the standard energy conditions, and without having to introduce new boundary conditions, such as the Hartle-Hawking 'no-boundary proposal'; they are consequences just of the quantum corrected Einstein's equations. Thus, the existence of the bounce and the upper bound on curvature and matter density can be directly traced back to quintessential features of quantum geometry. These considerations have been extended beyond the spatially flat FLRW models to include spatial curvature, non-zero cosmological constant, anisotropies (see, e.g., [2,42] and references therein) as well as the simplest inhomogeneities captured by the Gowdy models in GR, and also to the Brans-Dicke theory (see, e.g., [53,54]). Taken together, these results bring out the robustness of the LQC bounce.
Since the area gap plays an important role in the LQC dynamics, before concluding this subsection, we will make a small detour to explain how its numerical value∆ 5.17 2 Pl is arrived at. Recall, first, that in QCD there is a quantization ambiguity -parametrized by an angle θ-because of the freedom in adding a topological term to the action. One encounters a similar quantization ambiguity in LQG (again associated with the freedom to add a term to the action that does not affect equations of motion), encoded in the so-called Barbero-Immirzi parameter, γ > 0, which trickles down to the expressions of observables on H, such as the area operatorÂ S . The eigenvalues ofÂ S are discrete in all γ-sectors. But their numerical values are proportional to γ and vary from one γ sector to another. Observables also have a θ dependence in QCD and the value of θ that Nature has selected is determined experimentally. In LQG, a direct measurement of eigenvalues of geometric operators would determine γ. But of course such a measurement is far beyond the current technological limits. However one can use thought experiments. Specifically, in LQG the number of microstates of a black hole horizon grows exponentially with the area, whence one knows that the entropy is proportional to the horizon area [55,56]. But the proportionality factor depends on the value of γ. Therefore if one requires that the leading term in the statistical mechanical entropy of a spherical black hole should be given by the Bekenstein-Hawking formula S = A/4 2 Pl , one determines γ and thus the LQG sector Nature prefers. In this sector the explicit value of the the area gap yields∆ 5.17 2 Pl [57][58][59][60] (and the leading term in the entropy of more general black holes -not necessarily spherical-agrees with the Bekenstein-Hawking formula). This is the value used in LQC calculations. This concludes our broad-brush overview of how quantum geometry considerations lead to a natural resolution of the big bang singularity in LQC. The resolution has been analyzed in detail in a large number of LQC papers, using Hamiltonian, cosmological-spinfoam and 'consistent histories' frameworks (see, e.g., [2,42,51,61,62]).
Remark: Recently some concerns have been expressed about the simplicity of the LQC description of the early universe, and on whether "general physics principles of effective field theory and covariance" have been appropriately incorporated [63]. Many of the specific technical points were already addressed, e.g., in [42,64,65] and in the Appendix of [66]. In addition, we would like to clarify possible confusion on the following points. First, although 'effective equations' are often used in LQC, conceptually they are on a very different footing from those used in effective field theories: One does not integrate out the UV modes of cosmological perturbations. The term 'effective' is used in a different sense in LQC: these equations carry some of the leading-order information contained in sharply peaked quantum FLRW geometries Ψ(a, φ). As we will see in Section III B, equations satisfied by the cosmological perturbations are indeed covariant. On the issue of simplicity of the LQC description, we note that in the 1980s it was often assumed that the early universe is irregular at all scales and therefore quite far from being as simple as is currently assumed at the onset of inflation. Yet now observations support the premise that the early universe is exceedingly simple in that it is well modeled by a FLRW spacetime with first order cosmological perturbations [24]. Therefore, although a priori one can envisage very complicated quantum geometries, it is far from being clear that they are in fact realized in the Planck regime. Nonetheless, one should keep in mind that, as in other approaches to quantum cosmology, in LQC the starting point is the symmetry reduced, cosmological sector of GR. Difference from the Wheeler-DeWitt theory is that one follows the same systematic procedure in this sector as one does in full LQG. But the much more difficult and fundamental issue of systematically deriving LQC from full LQG is still open mainly because dynamics of full LQG itself is still a subject of active investigation. See, e.g., [67,68] as illustrations of the current status.

B. Cosmological perturbations in the pre-inflationary era of LQC
In inflationary paradigms the Mukhanov-Sasaki scalar modes 3 of cosmological perturbations are represented by quantum fieldsQ that propagate on a background FLRW metric g ab . The use of a classical background geometry is justified since, as explained above, spacetime curvature is twelve orders of magnitude below the Planck scale even at the onset of inflation. However, to extend the paradigm all the way to the LQC bounce, one has to replace the metric g ab of GR with an LQC wave function Ψ(a, φ) because assumptions underlying quantum field theory (QFT) on curved spacetimes fail in the Planck regime. At first the task seems daunting: How do you evolve quantum fields when you have only a probability distribution Ψ(a, φ) for various spacetime geometries rather than a single metric g ab ? Fortunately, there is an unexpected simplification [3,5,69]: So long as Ψ(a, φ) is sharply peaked, and the back reaction of the perturbationsQ on the background quantum geometry Ψ remains negligible, dynamics of quantum fieldsQ on Ψ is extremely well-approximated by that of quantum fieldsQ propagating on a smooth, quantum corrected FLRW metricg ab which is constructed in a precise manner from Ψ. As one would expect, coefficients ofg ab depend on . In the literature,g ab is often called the dressed metric. It is 'dressed' by certain quantum fluctuations in Ψ(a, φ) specified below; it carries the information in the quantum geometry Ψ(a, φ) that the propagation of cosmological perturbations is sensitive to.
The construction of the dressed metricg ab can be summarized as follows. Recall first that in the standard inflationary scenario, the Mukhanov-Sasaki quantum fieldQ satisfies a wave equation ( + U/a 2 )Q = 0 where is the d'Alembertian w.r.t. to the background FLRW metric g ab (satisfying the unperturbed, zeroth order Einstein's equations) and U is constructed from the inflationary potential and the background FLRW solution (see, e.g., [5]). At the classical level, this evolution equation can be derived starting with the full Hamiltonian constraint of GR coupled with the scalar field, and then appropriately truncating it to second order in perturbations [4]. In LQC, the background quantum geometry Ψ(a, φ) satisfies the zeroth-order LQC Hamiltonian constraint. The scalar modeQ propagates on this Ψ(a, φ) and its dynamics is governed by the appropriate second order truncation of the full Hamiltonian constraint. If the state Ψ(a, φ) is sharply peaked and the back reaction of the perturbationQ is negligible, then one has the following result [69]: Propagation of Q on the quantum geometry Ψ(a, φ) is very well approximated by that of a quantum field Q satisfying (˜ +Ũ/ã 2 )Q = 0. Here˜ is the d'Alembertian with respect to the dressed metricg andŨ(φ) is the dressed effective potential All operators and their expectation values refer to the Hilbert space of the background FLRW quantum geometry: the expectation values are taken in the state Ψ,Ĥ is the 'free' Hamiltonian in absence of the inflaton potential, andâ(φ) is the (Heisenberg) scale factor operator [5,69]. At first, the result seems surprising. But physically it can be understood using a simple analogy with propagation of light in a medium such as water. In the full quantum description, individual photons interact with the molecules of the material. However, the key features of propagation can be extracted simply by computing a few macroscopic parameters such as the refractive index and birefringence that can be extracted from the microstructure of the material. Other details of the quantum state of the medium are not important to study propagation. In this analogy, the cosmological perturbation plays the role of light and quantum geometry, the role of the medium. To determine the propagation ofQ, one needs to extract onlyã,η andŨ from the quantum state Ψ(a, φ). The rest of the very rich information of quantum geometry it contains is not directly relevant. (Incidentally, the tensor modes satisfy the wave equation for the same dressed metricg ab ; as in standard inflation, there is no dependence on the potential.) It is clear from the form of equations (3.2) and (3.3) that the expressions of the dressed metric and the dressed potential could not have been guessed a priori. They resulted from explicit, detailed calculations [5,69]. The observable predictions reported in section IV are obtained by first calculating the dressed metric and the dressed potential starting from the given quantum geometry Ψ(a, φ) and then evolving the scalar mode using (˜ +Ũ/ã 2 )Q = 0.

Remarks:
1. For clarity, let us spell out the conceptual elements of the procedure used to extract dynamics of cosmological perturbations since there is occasional confusion on this point. The starting point in LQC is the Hamiltonian formulation of GR coupled to the inflaton (albeit in the connection variables used in LQG). One then extracts the sector of full GR that corresponds to the homogeneous isotropic fields (which serve as the background) together with first order perturbations. It is this classical theory that is then quantized using LQG techniques [4]. Dynamics of quantum perturbations are governed by the Hamiltonian constraint operator of the truncated sector, where both the background geometry and perturbations are treated quantum mechanically. One does not simply assume that perturbations satisfy linearized equations of GR on a bouncing classical metric. That the dynamics of perturbations is well approximated by quantum field satisfying an evolution equation involvingg ab andŨ is a result that holds under conditions spelled out above. Note also that the equation is covariant w.r.t.g ab andg ab rapidly tends to the classical FLRW metric of GR outside the Planck regime.
2. Initially, analysis of [69] suggested that the propagation ofQ on the quantum geometry Ψ(a, φ) would be exactly the same as that on the corresponding dressed metric 3.1 and potential (3.3) for any Ψ(a, φ). However, Kaminski later found [70] that there is a subtle infrared problem (that can be missed in numerical simulations since they have to use an infrared cutoff). Kaminski, Kolanowski and Lewandowski [71] then showed that, as a result, the implications of [69] are not as general; the result would not hold without restrictions on the background quantum geometry Ψ(a, φ). This situation is qualitatively similar to that, e.g., in quantum electrodynamics which also faces infrared issues in rigorous treatments. However, in QED the ensuing difficulties can be avoided by focusing just on those quantities that are 'infrared safe'. One can adopt a similar strategy in LQC by introducing suitable infrared safe observables through regularization. Furthermore, for states Ψ(a, φ) that are sufficiently sharply peaked, the regularization ambiguity is completely negligible. Previous calculations of power spectra in LQC (e.g., [3,5,10,11,19,21]), as well as the current investigation, use states that are sufficiently sharply peaked in this sense, whence the use of dressed metric is justified in spite of the infrared difficulties.

C. Primordial spectrum: Why pre-inflationary dynamics matters
A natural question now is whether the pre-inflationary phase of LQC dynamics described in the last two subsections has any observable consequences. Let us therefore focus on the observable modesQ k . These have co-moving wavenumbers k in the range ∼ (0.1k , 300k ), where k = 0.002Mpc −1 is the WMAP pivot scale. The evolution equation for these modes implies that, they 'experience' curvature in the background metricg ab only if their physical wavelength λ(t) = a(t)/k is comparable or larger than the radius of curvature r curv (t) = (6/R) 1 2 ofg ab corresponding to the scalar curvature R at that time. Let us denote by t = t the time at which the relevant slow roll phase starts; this is our onset of inflation. Therefore, a few e-folds before and after t , the observable modes propagate as though they are in flat spacetime and therefore do not get excited by the background geometry. What happens in the distant past? The left panel of Fig. 3 shows the evolution of r curv (blue solid curve) and of λ of observable modes (the gray shaded band), both in GR. Note that in the pre-inflationary epoch r curv is far from being constant whence the spacetime metric is very different from the de Sitter metric. Since the scalar curvature R diverges at the big bang, r curv goes to zero. Because the scale factor a of the classical FLRW metric goes to zero, physical wavelengths λ also goes to zero at the big bang. However, they do not go to zero as fast as r curv . Therefore, as one approaches the big bang in the past evolution, all observable modes exit the curvature radius, 'experience' curvature at sufficiently early times and get excited. These excitations have to be delicately fine-tuned for the state to be in the BD vacuum later on, at t = t i.e. at the onset of inflation. Put differently, the Heisenberg state representing the BD vacuum at the onset of inflation is an unnatural choice from the perspective of the Planck regime because it carries certain delicately choreographed excitations there. Of course, one can argue that the quantum field theory in curved spacetime cannot be extrapolated to the Planck regime. But by itself this argument does not provide a justification for using the BD state at t = t either. In LQC, the situation is quite different. Because the scalar curvatureR of the dressed metricg has a finite upper boundR max 62 −2 Pl , reached at the bounce, r curv reaches its minimum value r min curv 0.31 Pl , whence it is only those modes which satisfy λ r min curv at the bounce that experience curvature in their evolution from the bounce to the onset of inflation. In our approach, (as discussed below) the background quantum geometry Ψ(a, φ) is such that only the longest wavelength observable modes satisfy this inequality. This feature is shown in the right panel of Fig. 3. Therefore, all but the longest wavelength observable modes propagate from the bounce-time to the onset of inflation as though they are in flat spacetime and hence it is natural that they be in the BD vacuum at t = t . It is only the longest wavelength modes that will be excited and hence not in the BD vacuum. But won't these excitations get just washed away during inflation? The answer is in the negative: because of spontaneous emission, the number density of these excitations remains constant [72,73]. As a result, as we will see in section IV, the primordial spectrum does differ from that based on the SA, but only at largest angular scales.
Note that there is a deep interplay between the UV and IR in LQC. As we saw in section II A, in LQC it is the UV modifications of GR in the Planck regime that tame the big bang singularity and make all physical quantities finite. As a result we have a finitẽ R max and a non-zeror min curv . It is then natural for all but the longest wavelength observable modes to be in the BD vacuum at the onset of inflation. While the LQC corrections to the background geometry are significant in the UV, their effect on cosmological perturbations is non-negligible only in the IR. This is the point that was highlighted in the abstract and section I.
Finally, to obtain specific predictions, one needs a quantum state of geometry Ψ(a, φ) and a quantum state ψ(Q, φ) of scalar modes. At this point, different approaches within LQC make different choices (see, e.g., [5-7, 9-11, 14-21]). In this paper, we use the procedure introduced in [12,13]. Strategy is to select these states by introducing some trial principles that relate properties of quantum geometry in the Planck regime with the late time geometry (which can be taken to be that given by general relativity to an excellent degree of approximation). One can then work out the observable consequences. If any prediction is ruled out by observations, one would return to the drawing board and seek alternate principles that would lead to viable states. If predictions are confirmed by observations, one would build confidence in the general direction and attempt to put the principles on a firmer and more satisfactory footing. The currently used principles are somewhat analogous to the Bohr model of the hydrogen atom in the early days of quantum mechanics. While in retrospect it is naive in some fundamental respects, nonetheless the Bohr model was useful because it captured some essential features of the final, correct description of the hydrogen model. Given any inflationary potential, the first principle constrains Ψ(a, φ), determining the number of e-folds in the pre-inflationary dynamics, i.e., between the LQC bounce and the onset of inflation [13]. This restriction enables one to determine the physical wavelength λ of observable modes at the bounce-time, and hence the angular scale on which the LQC corrections to the power spectrum are significant (which turned out to correspond to 30 for the potentials considered.) The band of these modes is depicted by the (gray) shaded region in the right panel of Fig. 3, together with its relation to the curvature radius r curv , depicted by the solid (blue) line. Finally, the principle that determines the quantum state ψ(Q, φ) of scalar modes [13] involves a quantum generalization of Penrose's Weyl curvature hypothesis [74] in the Planck regime near the bounce, which physically corresponds to requiring that the state should be 'as isotropic and homogeneous in the Planck regime, as the Heisenberg uncertainty principle allows'. 4 While the first principle sets the scale at which the LQC primordial spectrum ceases to be nearly scale invariant, and thus differs from that given by the SA, the second led to the conclusion that there is power suppression with re-spect to the SA rather than power enhancement. Detailed calculations are needed to obtain the precise degree of suppression.

IV. RESULTS
In this section we present the main results of this paper, obtained using the LQC summarized in section III. In terms of more commonly used wavenumbers, the new length scale r min cuv ∼ 0.31 Pl introduced by LQC provides a new physical scale k LQC 3.21 −1 Pl , and the primordial spectrum differs from the SA for k phy k LQC at the bounce. In section IV A we first discuss these LQC corrections and then their effect on the observed C XY correlations (where XY refers to TT, TE, EE and φφ). We also present predictions for the value of the optical depth τ and for the BB power spectrum that could potentially be tested using future observations. In IV B, we show that the LQC predictions for the C XY 's lead to resolution of power suppression and lensing amplitude anomalies. In IV C, we will show that the interplay between LQC and observations is a 2-way bridge, in that the CMB observations can also be used to constrain the value of the area gap ∆, the most important of fundamental microscopic parameters of LQG.

A. Power spectra
In our LQC model, the physical principles used to select the background quantum geometry imply that the corresponding ΛCDM universe has undergone approximately 141 e-folds of expansion since the quantum bounce until today [13]. Therefore, the characteristic LQC scale k LQC corresponds to the co-moving wavenumber k o 3.6 × 10 −4 Mpc −1 which sets the scale below which LQC corrections to the primordial scalar power spectrum become important. In particular, the calculations show that for scales k 10k o the power is suppressed whereas for scales k k o the power spectrum is essentially scale invariant as in the SA. This behavior is captured in following modification to the SA for the primordial power spectrum: where f (k) is the correction factor (which is equivalent to the ratio of the power spectra in LQC to that in the SA).
The left panel of Fig. 4 shows the correction factor f (k) plotted with respect to the co-moving wavenumber. It is evident from the plot that for both the Starobinsky and quadratic potentials, relative to the prediction of the SA, there is power suppression for long wavelength modes corresponding to k 10k o . The origin of this difference lies in the fact that, during the pre-inflationary era, the physical wavelength of these modes is sufficiently large to 'experience' the background curvature in the Planck regime near the bounce, leading to excitations over the standard Bunch-Davis (BD) state at the onset of inflation. On the other hand, the physical wavelengths of modes with k k o are much smaller than the curvature radius throughout the pre-inflationary phase, including the deep Planck regime near the bounce. Therefore their state is practically the same as the standard BD state at the onset of inflation. Note that even for modes with k 10k o , LQC corrections are significant only near the bounce. During this epoch the energy density of the scalar field is dominated by the kinetic term. Therefore one would expect the inflationary potentials to have negligible effect on the evolution of the background geometry and perturbations in the deep Planck regime, which is when the LQC corrections are imprinted on the modes of scalar perturbations. This expectation is explicitly borne out in the left panel in Fig 4: the primordial spectra for Starobinsky and quadratic potentials are essentially identical. Analytical considerations of [75] suggest that this feature will persist for a large class of inflationary potentials. Finally, note that the LQC primordial spectrum has a turnaround at very large wavelengths whose origin is related to the sudden spike in r curv near the bounce. Although these modes are not in the observable range for CMB, it is of interest to better understand the origin of this growth in power for very small k because, together with large and strongly scale dependent non-Gaussianity (which is expected from the LQC bounce [15,19]) it could lead to a coupling between the long wavelength modes and the modes observable in CMB, resulting in a modulation of the primordial power spectrum of the observable modes. Such a non-Gaussian modulation could explain the dipolar modulation of CMB and preference for odd parity that has been observed in the CMB [20,21].
Let us now turn to the LQC predictions for observable power spectra. For definiteness, in these plots we work with the Starobinsky potential for the SA and LQC theoretical predictions. On the observational side, 'PLANCK 2018' refers to the TT+TE+EE+lowE+lensing 2018 data released by the PLANCK collaboration. In order to obtain constraints on the parameters in LQC and SA models, we used the publicly available software package COSMOMC [76] which supports the likelihood code used in the original 2018 PLANCK analyses. COSMOMC is based on the Markov-Chain-Monte-Carlo (MCMC) procedure for estimation of parameters based on maximum likelihood analysis. For a given theoretical model with a number of parameters (here the ΛCDM model with 6 parameters), one builds a Markov chain from randomly sampled parameter values from a space predefined by priors. Each Markov chain begins with a random sample from the prior-constrained parameter space. Subsequent steps in the chain are then selected via Metropolis-Hasting algorithm: a Monte Carlo sample from the parameter space is proposed and then accepted or rejected based on the likeli-hood function which quantifies the degree of agreement between the theoretical prediction and the experimental data (in this paper we work with the same likelihood function as used in PLANCK 2018 papers [25].). This procedure is repeated until a convergence criteria is satisfied. The chain of accepted values are further trimmed and thinned in order to remove dependence on the initial data point and correlation between subsequent steps. The converged Markov chain thus obtained approximates the posterior distribution of the parameters to be estimated. In order to obtain one-dimensional constraints on (or twodimensional correlation between) the individual parameters, marginalization procedure is used. For details of the MCMC techniques adapted to cosmological settings see [76]. The right panel of Fig. 4 shows the theoretical predictions as well as the observed power spectrum, with its error bars. The plots clearly show that the primordial power suppression of LQC for small k leads to power suppression in the TT spectrum for multipoles < 30, relative to the predictions from the SA. The same primordial power spectrum can also be used to compute the predicted TE and TT correlation spectra. Fig. 5 shows the TE (left panel) and EE (right panel) spectra as observed by Planck 2018, along with those obtained using the best fit ΛCDM models with LQC and the SA. Note that in the right panel of Fig.  4 and in Fig. 5 the LQC and the SA plots are shown for their corresponding marginalized mean values of the 6 parameters shown in Table I. By contrast, in [13] both the LQC and the SA curves were plotted for the same values of 6 parameters of the ΛCDM model, coming from the SA. As Table I shows, the mean value of optical depth τ is > 9% higher in LQC than with the SA and electric polarization is quite sensitive to τ . Interestingly, this change in the best fit ΛCDM universe has the effect of slightly reducing suppression in the passage from the primordial power spectrum to the observed one. As a result, the LQC power suppression in the TT, TE and EE spectra we now find is somewhat less pronounced than it was in [13].
Finally, the left panel of Fig. 6 shows the lensing correlation spectrum C φφ reported by the PLANCK collaboration with their 2018 data, along with the LQC and the SA predictions. For this figure we have again used the mean marginalized values of the 6 parameters shown in Table I, and the lensing amplitude is fixed to A L = 1, in accordance with the base ΛCDM model. As with the TT, TE and EE spectra, the lensing spectrum also shows suppression at large angular scales corresponding to < 30. However, observations are    quite sparse for low . Interestingly, for 30 < < 100, the LQC prediction for C φφ is slightly larger than that for SA. (This is in fact inevitable; see section V.) This behavior leads the LQC spectrum to fit slightly better with the observed data in the range 30 < < 100 and hints towards resolving the lensing amplitude anomaly in LQC without having to introduce additional modifications to the standard ΛCDM model. In the next subsection we will see that this possibility is in fact realized.
We will conclude this discussion of observable implications of LQC with predictions for future missions. First, as we noted in section II B, currently the optical depth τ is the least accurately measured of the 6 ΛCDM parameters, with a relative error of ∼ 13%. The LQC value is some 9.8% higher than that in 'the universe according to PLANCK'. This prediction will be tested by the future observation of global 21cm evolution at high redshifts that is estimated to reach a percent level accuracy in the measurement of τ [33]. As for the observable power spectra, to date the PLANCK satellite has provided the best full sky measurement of the CMB anisotropies. However, there is still scope for improvement in the measurement of electric polarization, and the odd-parity magnetic polarization is yet to be detected. Therefore, several space-based mission have been planned to further improve the polarization measurements. In addition to the predictions for the TE and EE power spectra discussed above, our LQC model also makes predictions for the BB power spectrum. The right panel of Fig. 6 provides the prediction for the unlensed BB power spectrum both from the SA and in LQC. Recall that the tensor-to-scalar ratio r depends on the potential of the inflationary model. But, being a ratio, it is the same in LQC as from SA [5]. We have set r = 0.0041, its value for the Starobinsky potential. As in the case of other four spectra, we observe a relative suppression of power at low multipoles. However, this is also where the reionization bump occurs. Since LQC predicts a larger value of optical depth, the B-B power suppression is lower than what one might have expected from the primordial power suppression (and using the same value of τ for both SA and LQC). Nonetheless, it may be possible to test this prediction against the data from the future B-mode missions such as LiteBIRD [29], Cosmic Origins Explorer [30], ECHO [31] or Probe Inflation and Cosmic Origins (PICO) [32] (which should observe the BB spectrum if r 0.001). As discussed in Sec. II B, in our approach the pre-inflationary phase of LQC dynamics enables one to address two CMB anomalies. The first, shown in the left panel of Fig. 2, is that the observed angular correlation function C TT (θ) remains close to zero for angles θ > 60 • in contrast to the predictions of the standard ΛCDM model based on the SA. Let us now examine the prediction of C TT (θ) from LQC. This prediction is plotted as the dashed(blue) curve in Fig. 7, along with the prediction from the SA shown as the (red) The tension between observations and the theoretical prediction from the the SA is encapsulated by the fact that S SA 1/2 is almost 35 times larger than S Planck . This is the power suppression anomaly. This discrepancy is appreciably reduced in LQC since S LQC 1/2 is ∼ 1/3 of S SA 1/2 . As indicted in section II B, because C(θ) for different values of θ are correlated, and because one has to use a masking procedure in the data analysis near θ = 180 • , the task of providing the 1σ and 2σ contours around the LQC plots is challenging, requiring manipulations of a large covariance matrix in the data analysis, as well as a detailed understanding of aspects of the instrument. It would be of considerable interest if the CMB experts could provide these plots starting with the LQC TT-power spectrum. Note also that because power is suppressed for 30, if LQC results were used in the PLANCK analysis, error bars for low would also be reduced and we would have a sharp measure of the LQC alleviation of this anomaly.
The second anomaly, shown in the right panel of Fig. 2, is the lensing amplitude anomaly. This anomaly is not directly observed in the C XX plots or in C(θ) plots, but arises when one performs a consistency check of the ΛCDM model. Instead of fixing the lensing amplitude to A L = 1 (as is assumed in standard ΛCDM model) one allows it to vary along with the standard 6 parameters, i.e., one now analyzes a 7 parameter model. The anomaly lies in the finding that A L = 1 is beyond 1σ error bar as shown in the right panel of Fig.  2. This led the authors of [23] to conclude that there is a possible "crisis in cosmology" because, to alleviate this problem, one would need to introduce spatial curvature which creates significant departures from the observed power spectra as small angular scales. As shown in Fig. 8, however, repeating the analysis with the LQC power spectrum restores A L = 1 within 1σ contour thereby resolving the lensing amplitude anomaly and avoiding the hint of a potential "crisis".
To summarize, the observed CMB power spectrum has many non-trivial and interesting features at small angular scales and observational error bars are small in this regime. It is remarkable that these features are correctly predicted by the ΛCDM model based on the SA. The LQC corrections could well have led to discrepancies with these successful predictions of the SA. That does not happen. Rather, there are departures only at large angular scales leading to a suppression of the primordial power for small k. This difference from the nearly scale invariant standard ansatz (2.1) leads to the alleviation of two anomalies in the CMB. Pl that lies to the right of the peak. Thus, there is an unforeseen internal consistency.
As we saw in section III, the area gap ∆ is the key microscopic parameter that determines values of important new, macroscopic observables such as the matter density and the curvature at the bounce. Its specific value,∆ = 5.17 2 Pl , is determined by the statistical mechanical calculation of the black hole entropy in loop quantum gravity (see, e.g., [55,56,59,60]). Results reported in the last two subsections are based on this value. However, with the CMB observations at hand, we can now, so to say, turn the tables and regard ∆ as an additional parameter and use the CMB observations to constrain its value.
Thus, let us consider ∆ to be a free parameter and obtain a posterior probability distribution for its value by letting it vary along with the 6 parameters of the ΛCDM model. Note that this procedure is similar to the self-consistency check that is performed for the lensing amplitude A L . The ΛCDM universe determined by the SA does not pass that test because, when it is allowed to vary, A L prefers a value that is greater than 1 and the discrepancy is significant in that the value 1 lies outside the 1σ contour of the best fit value [27]. Is there perhaps a similar tension here? More precisely, does the value∆ 5.17 2 Pl obtained from black hole entropy calculations lie within the 68% confidence contour of the value preferred by the CMB observations? If it does not, LQC would fail the self-consistency test at the 1σ level. If it does, there would be an unforeseen coherence between detailed conclusions drawn from vary different considerations: the LQG analysis of black hole entropy and the LQC investigation of the very early universe! Fig. 9 shows the one-dimensional posterior distribution of R B . The the best-fit value -the peak in the distribution-is at ∆ = 2.18 2 Pl while the marginalized mean value is 3.86 2 Pl . with the following constraint: Clearly, the value∆ 5.17 2 Pl chosen in Sec. III C and used in this paper, is within 68% (1σ) confidence level of the constraint obtained from Planck 2018. This not only indicates a synergy between the fundamental theoretical considerations and observational data, but also provides internal consistency of the LQC model.

V. DISCUSSION
To determine the six parameter ΛCDM universe we live in, the PLANCK team began with the standard ansatz (2.1) for the primordial spectrum and used known astrophysics to determine the observable power spectra for various values of the six cosmological parameters, A s , n s , Ω b h 2 , Ω c h 2 , 100θ M C . By comparing these theoretical predictions with observations, they determined the marginalized mean values of the six parameters (together with the error bars corresponding to 68% confidence level). This is the ΛCDM universe according to PLANCK. This procedure has had tremendous success, especially with the finer features of the power spectra at small angular scales. However, there are also some anomalies. Since they are only at ∼ 2−3σ level, the statistical significance of any one anomaly is low. However, taken together, two or more anomalies imply that we live in an exceptional realization of the six posterior distributions provided by this procedure.
One can view these anomalies as potential gates to new physics. Indeed, the PLANCK collaboration has emphasized this possibility in its 2015 [39] as well as 2018 data releases [40]. As the second of these papers points out, "...if any anomalies have primordial origin, then their large scale nature would suggest an explanation rooted in fundamental physics. Thus it is worth exploring any models that might explain an anomaly (even better, multiple anomalies) naturally, or with very few parameters." LQC researchers have followed up on this suggestion. In this paper we presented a concrete realization of this idea. Specifically, one begins with the observation that, in the standard procedure summarized above, the theoretical input, beyond known astrophysics, is the SA, motivated by the inflationary scenario. It assumes that the primordial spectrum is nearly scale invariant and can be characterized just by two numbers A s and n s across all wavenumbers k. However, in LQC the resolution of the big bang singularity introduces a new scale k LQC and quantum gravity corrections in the pre-inflationary phase of dynamics appear for k k LQC . For these small wavenumber modes, the LQC primordial spectrum is no longer nearly scale invariant, whence there is departure from the predictions drawn from the SA.
Several closely related approaches have been used in LQC to probe the effects of this pre-inflationary dynamics (see, e.g., [3-16, 19, 21]). Our approach has two main ingredients: (i) the use of sharply peaked quantum states Ψ(a, φ) for the background quantum FLRW geometry, that then enable one to well-approximate the dynamics of cosmological perturbations on the quantum geometry Ψ(a, φ) by that on a quantum corrected 'dressed metric' g ab [5,56]; and, (ii) the use of certain principles to select Ψ(a, φ) and the quantum state ψ(Q, φ) of the scalar mode of cosmological perturbations for any given inflationary potential [12,13]. For any given inflationary potential, the principle used to select Ψ(a, φ) limits the number of e-folds during pre-inflationary dynamics, thereby implying that the modes that receive significant LQC corrections in the primordial spectrum correspond to the large angular scales 30. The principle used to select ψ(Q, φ) implies that in the primordial spectrum there is power reduction (rather than enhancement) in these modes. This then translates to a power suppression for 30 in the observed power spectra. For modes with 30, the LQC power spectra are indistinguishable from those obtained using the SA. Thus, LQC predictions leave the highly successful predictions of standard inflation at small angular scales unaffected, but modify the predictions at large angular scales.
Details of the LQC pre-inflationary dynamics reveal some interesting facts. First, the quantum geometry effects on the background FLRW geometry are dominant only in a short interval around the bounce. Second, it is the modes whose physical wavelength λ phy is longer than the curvature radius r curv during the pre-inflationary evolution that fail to be in the BD vacuum at onset of inflation. In the observable band, only the longest wavelength modes are thus affected and they satisfy λ phy r curv only for 2-3 e-folds after the bounce. Thus, the background quantum geometry as well as the quantum perturbations Planck receive non-negligible LQC corrections during a very short duration. Yet these corrections lead to observable effects in that they alleviate some anomalies. Third, while we did not discuss tensor modes in this paper, their power spectra have the same behavior as that of scalar modes and, given an inflationary potential, the tensor to scalar ratio r does not receive LQC corrections (within accuracies reported here). Finally, there is an unforeseen interplay between the UV and the IR: While it is the UV modifications of GR that lead to the singularity resolution and create the new LQC scale k LQC , the structure of the evolution equations satisfied by cosmological perturbations is such that it is the IR modes with k k LQC that are affected during their pre-inflationary evolution. It is this unforeseen cosmic tango between the very small and the very large that is responsible for the alleviation of the two anomalies discussed in this paper.
Given that LQC simultaneously alleviates the power suppression and the lensing amplitude anomalies, it is worth investigating a more general question: Are the two conceptually related? As reported in [22], the answer is in the affirmative. Since this relation seems not to have been noticed before, we will make a small detour to explain it. Let us begin by assuming that there is some mechanism -not necessarily originating in LQC-that provides a primordial power spectrum of the form with f (k) < 1 for k < k • and f (k) = 1 for k > k • for some k • , and compare and contrast the new ΛCDM universe obtained from this modified ansatz with that given by the SA of Eq. (2.1). In the first step, we can restrict our analysis only to smaller angular scales (k k • ). Then, the primordial spectrum in both schemes would be the same, whence we would obtain the same marginalized mean values of the six cosmological parameters. Denote byÅ s the marginalized mean value of the scalar amplitude A s thus obtained. In the second step, let us consider the full range of observable modes, including k ≤ k • . Now, given that the observations show that the TT power is suppressed at large-scales, i.e., for k k • , if one uses the SA+ΛCDM model the marginalized mean value A SA s using the entire k range will be lower thanÅ s . By contrast if the primordial power spectrum is of the form of Eq. (5.1), the initial power is already suppressed by f (k). Therefore,Å s will not have to be lowered as much to obtain the marginalized mean value A new s . Thus, we have The key point is the last inequality: A new s > A SA s . (We spelled out the argument because at first it seems counter-intuitive that power suppression leads to a larger A new s . But note that power is suppressed only for low k.) Next, we know that for large k, the product A s e −2τ is fixed by observations. Hence, it follows that the best fit values of the optical depth in the two scheme must satisfy τ new > τ SA . Finally, from the very definition of lensing amplitude, the value of A L is anti-correlated to the value of A s . Therefore, it follows that we have the inequality A new L < A SA L . Thus in any theory that has primordial spectrum of the form (5.1), A s , τ and A L will have the same qualitative behavior as in LQC, and hence the tension with observations would be reduced. What LQC provides is a precise form of the suppression factor f (k) from 'first principles', and hence specific quantitative predictions. Also, recall from section IV A that in our analysis the LQC f (k) also came with a specific value k o 3.6 × 10 −4 Mpc −1 for k • . Other mechanisms could well lead to a very different value. If so, in the observed power spectra suppression would arise at a very different value of . Finally, the specific f (k) computed from LQC also leads to other predictions -e.g., for the BB power spectrum discussed in section IV A-that need not be shared by other mechanisms.
In this respect, it would be of interest to compare the LQC predictions with those that result from effective field theoriesà la Ginsburg and Landau, where slow roll inflation is generically preceded by a fast roll phase that leads to a suppression of CMB quadrupole [77]. This will require a calculation of angular power spectrum C(θ), and of the measure S 1/2 of power suppression, using the marginalized posterior probabilities of the 6 cosmological parameters in this effective field theory approach, and a reanalysis of the lensing amplitude along the lines of section IV B. Yet other mechanisms have been proposed to account for power suppression at large angular scales in the context of GR (see, e.g. [78][79][80][81][82][83][84][85]). These are compared and contrasted with LQC in Section 5 of [13]. Finally power suppression at large angular scales has been studied in the context of other bouncing models. In those discussions the bounce is often just assumed (as in [86]), or obtained by adding a scalar field with a negative kinetic term which violates standard energy conditions ( as in [87]). Our approach is different in that: (i) the bounce is a prediction of LQC; (ii) since the mechanism has its roots in quantum geometry effects underlying LQG, additional scalar fields or violations of energy conditions are not involved; (iii) the standard inflationary potentials are used without adjustments to provide a fast roll phase; and, most importantly, (iv) our goal is to investigate whether the CMB observations can inform quantum gravity and vice versa.
In our view, the big bounce and the pre-inflationary dynamics of cosmological perturbations are on a robust footing although the discussion would benefit from a further sharpening of the detailed arguments that led us to the dressed metricg ab . The part of the analysis that is on a less solid footing concerns the specific principles [12,13] that were used to select the wave function Ψ o (a, φ) of the background quantum geometry and the state ψ(a, φ) of perturbations. Note that these choices are necessary to make predictions in any approach that starts in the Planck regime. Indeed, even in standard inflation one has to assume that the perturbations are in the BD vacuum at the start of slow roll, and as discussed in section III C, it is difficult to justify this assumption from first principles. The fact that the principles led us to predictions that not only reproduce the successes of standard inflation, but also alleviate the tension associated with two anomalies, is an indication that they set us on the right track. However, they should be regarded as tentative first steps, to be improved upon and sharpened in future. In particular, the effect of the duration of the reheating epoch is yet to be investigated in detail. Another direction for future work is suggested by a second approach that has been used to alleviate CMB anomalies using pre-inflationary dynamics [19,21]. The point of departure is the same as in our approach: the bounce sets a new scale and modes with the longer wavelength at the bounce get excited during the pre-inflationary evolution and are no longer in the BD vacuum at the onset of inflation. However, the primordial spectrum does not exhibit power suppression (4.1) as in our case. Instead, a key role is played by the superhorizon modes and their non-Gaussian correlations with the longest wavelength modes with those that are observable in the CMB. These correlations enhance the probability of finding certain features in the individual realizations of the primordial probability distribution, thereby alleviating anomalies associated with dipolar asymmetry and power suppression anomalies. It would be of considerable interest to re-examine all CMB anomalies from a perspective that combines strengths of the two approaches. Finally, all approaches to finding observable consequences of LQC in CMB, that we are aware of, assume an inflationary potential. At a certain level of discussion this is justified since a detailed analysis of the criticisms of the inflationary paradigm has concluded that, although important questions remain, the case for inflation has been "strengthened by the PLANCK data" [24]. However, from a fundamental perspective, introduction of an inflaton and a specific potential is ad-hoc. Now, in full LQG, the Einstein-Hilbert Lagrangian receives higher order (in particular R 2 ) corrections because the curvature operator is constructed from Planck scale Wilson loops and hence non-local at the fundamental microscopic scale. Therefore it is important to analyze if these corrections would provide a natural basis for inflation purely from gravitational considerations, as in Starobinsky inflation.