Symmetry, Integrability and Geometry: Methods and Applications Novel Possibility of Nonstandard Statistics in the Inflationary Spectrum of Primordial Inhomogeneities ⋆

Inflation is considered one of the cornerstones of modern cosmology. However, the account of the origin of cosmic structure, as provided by the standard inflationary paradigm, is not fully satisfactory. The fundamental issue is the inability of the usual account to point out the physical mechanism responsible for generating the inhomogeneity and anisotropy of our Universe, starting from the exactly homogeneous and isotropic vacuum state associated with the early inflationary regime. We briefly review this issue here together with the proposal to address this shortcoming in terms of a dynamical collapse of the vacuum state of the inflaton field, which has been considered in previous works. The main goal of this manuscript is to discuss certain statistical aspects associated with the collapse and its implications in the primordial spectrum, particularly those connected with the possible appearance of a novel type of unusual correlations.


Introduction
The inflationary paradigm represents one of the pillars of the ΛCDM cosmological model. It was initially proposed as a solution to the classical naturalness problems of the the big bang models [1,13,15,18,28,40], but its impact became even more significant when it came to be regarded as the natural mechanism to account for the seeds of cosmic structure [14,16,17,31,32,41]. However, as originally presented in [37], the standard picture suffers from a conceptual drawback: we know that our Universe is quantum mechanical and, thus, the classical descriptions must be regarded as nothing but short-hand and imprecise characterizations of complicated quantum mechanical states. The Universe that we inhabit today is certainly very well described at the classical level by an inhomogeneous and anisotropic classical state. That is, such description must be nothing but a concise and imperfect description of an equally inhomogeneous and anisotropic quantum state. Consequently, if we want to seriously consider a theory, in which the early quantum state of the Universe was perfectly symmetric (the symmetry being the homogeneity and isotropy), then we need to provide and explanation for why the quantum state that describes our actual Universe lacks such symmetries. Since there is nothing in the dynamical evolution (as given by the standard inflationary approach) of the quantum state that can break those This paper is a contribution to the Special Issue "Loop Quantum Gravity and Cosmology". The full collection is available at http://www.emis.de/journals/SIGMA/LQGC.html symmetries, then we are left with an incomplete theory. Therefore a natural inquiry arises: How does the actual inhomogeneities in the Universe we inhabit emerged?. The standard and phenomenologically successful accounts on this matter rely on the identification of the quantum fluctuations of certain observables associated to a homogeneous and isotropic Universe, with the averages over an ensemble of inhomogeneous Universes of their analogue classical quantities, i.e. the identification of quantum uncertainties and classical stochastic perturbations. Nevertheless, as discussed in [37,44], one can not avoid concluding that, using standard physics, there is no clear justification for such kind of identification. In fact, these shortcomings have been recognized by others in the literature [29,30,33,47].
The idea that has been considered in previous works [37,42,43,44,45], as a possibility to deal with that problem, involves adding to the standard inflationary model the hypothesis that the collapse of the wave function is an actual physical process that occurs independently of external observers. The collapse hypothesis relies on the idea that something that effectively might be described in terms of such a collapse of the wave function could have its origins in the passage from the atemporal regime of quantum gravity to the classical space-time description underlining the general theory of relativity. That is, in going from one description to the other, we might be forced to characterize some aspects of the underlying physics in terms of sudden jumps which are not compatible with the unitary Schrödinger evolution, and which modify the state of the Universe in an stochastic way. Therefore, being capable of transforming a condition that was initially homogeneous and isotropic into another one that is not. It is worth mentioning that, in fact, conceptually similar ideas have been considered by many physicists concerned with the measurement problem in quantum mechanics [3,12,34], and the proposal that some quantum aspects associated to the gravitational interaction might lie at the root of this effect only serves to make the idea even more attractive [8,9,35,36].
The detailed discussion of the conceptual problems and those associated with the standard inflationary paradigm have been presented in previous works by some of us and by others in [25,34,37,44]. We will not reproduce those arguments here. On the other hand, a detailed discussion of the statistical aspects concerning both our modified proposal involving the quantum collapse of the wave function and as well as various aspects of the standard picture are presented in [25]; and again we will not reproduce most of those arguments here. In this paper our aim is to present one novel possibility regarding the statistical aspects which might be present in the primordial spectrum. It would correspond to the existence of non-trivial correlations between the random variables characterizing the collapse of certain modes and its statistical properties. In particular, these features would give raise to a very particular deviation from the flat shape of the spectrum of the primordial fluctuations, and thus to a characteristic signature temperature anisotropies in the CMB. This is a crucial difference with the standard accounts, as the distribution of the different quantum field modes has previously been assumed to be Gaussian, i.e. that all the fundamental variables are uncorrelated.
As we will show next, the consideration of this possibility is enabled by the fact that the statistical aspects, in the collapse picture, are quite transparent and clearly identifiable. In fact, the motivation for the work presented in this article, is the formal and rigorous analysis presented in [7]. There, the collapse hypothesis is formally considered within the full framework of Quantum Field Theory in curved space-times, as this is the most reliable approach we can take in the absence of a complete theory of quantum gravity.
The paper is organized as follows. In Section 2, we will review some of the considerations contrasting the collapse proposal and the standard approach. We give a brief account of way we think that these ideas would fit into the more standard views of the quantum/gravity interface. We also introduce the Semiclassical Self-consistent Configuration program which is the main motivation behind the results presented in this article. In Section 3, we review the collapse models description for the inflationary origin of the seeds of the cosmic structure. In Section 4, we present how the stochastic nature of the collapse proposal might induce new statistical features in the observable spectrum of anisotropies in the CMB. Finally, in Section 5, we discuss our findings. The conventions we will be using include a (−, +, +, +) signature for the spacetime metric. We will use units where c = 1, but will keep the gravitational constant G and explicit throughout the paper.

General framework of the collapse model
The collapse model for the emergence of the seeds of structure during inflation was originally proposed in [37], and since then has been further developed [5,6,24,26]. In this section we will present a review of the ideas that we have explored and which are connected to the signature of novel statistical features explored in this paper.
Within the standard approach [2,20,21,27,48], the possibility of deviations from simple statistics, is often referred as the study of non-Gaussianity, and is a feature thought to be characterized only in the higher order statistical n-point functions. We will discuss a novel type of deviation from the standard statistics that is suggested by our approach and that could lead to modifications directly observable in CMB the spectrum.
In the standard inflationary paradigm, the perturbations of the field δφ and the perturbations of the curvature Ψ are treated as standard quantum fields 1 evolving in a classical quasi-De Sitter background space-time. The statistical aspects involve a series of identification between quantum n-point functions, statistical n-point functions for ensembles of universes and averages over our universe. Those issues are often left unspecified and the clear characterization of these aspects in the theoretical analysis is often missing. On the other hand, within the collapse proposal, it becomes quite evident where the stochastic nature of the problem resides. That is, the dynamical collapse of the wave function is taken to be a physical mechanism governed essentially by the stochastic rules inherent to the Quantum Theory. Although, we of course do not know what is the physics behind the collapse, we can parameterize it (in the next section we will detail this parametrization) and hope to learn something about its basic characteristics. The important point is that, if the initial state (say the Bunch-Davies (BD) vacuum state) collapses (or is reduced) to a different state, which lacks the symmetries of the original one, and in a manner in which nontrivial correlations might be generated, then the energy momentum tensor of this resulting state would lead, not only to the emergence of anisotropies and inhomogeneities, but might result in a particular type of signature appearing in the primordial spectrum of fluctuations as observable features. We will explain in the next sections how such primordial non-flat spectrum could arise, with the departures from flatness representing the novel signatures we have been discussing.

The semiclassical self-consistent conf iguration
As mentioned in the Introduction, the inflationary scenario is a very peculiar setting where one can explore some of the fundamental aspects of modern physics. In particular, it provides us with a situation in which the effects of gravity and quantum mechanics are strongly coupled. However, we do not have a complete theory of Quantum Gravity (QG). In fact, one of the problems that arises when trying to unify General Relativity with Quantum Mechanics is the so called problem of time in quantum gravity [19].
It is a known fact that the proposals for the canonical Quantization of Gravity leads to timeless theories. This is simply because the lapse and shift vectors, naturally associated with the notion of time in General Relativity, disappear from the theory (the technical reason is that the physical states of the theory need to satisfy the Hamiltonian and momentum constraints). The Loop Quantum Gravity approach (LQG), as the modern version of such ideas, does not escape from such problem. Let us analyze in brief detail our view on how to deal with it and its relation with our proposal (see [7,43]).
In LQG, one starts with the characterization of gravity in terms of the canonical variables given by he triad E a i and the connection variables A i a . These variables are serve to describe the geometry of a 3-spatial hypersurface, and characterize its embedding in a 4-dimensional (spacetime) manifold. At the quantum level, however, one passes to a description where the basic operators are associated with holonomies and fluxes (U γ , F s ) which correspond to integrals of the connection and triads over closed paths γ and 2-surfaces s respectively. The quantum states of the theory are thus characterized in terms of these variables, and the recovery of space-time notions from this setting is rather non trivial. In fact it seems likely one would need special kinds of semiclassical states which nevertheless should be solutions of the constraints of the theory including the Hamiltonian constraint. It is still a challenge to construct these states and the full recovery of space-time from them seems to be further away. However we can try to envision how this process would look when such technical difficulties are overcome.
Consider now, that we have incorporated into the theory, along with the geometrical variables, some matter fields, each described by a pairs of cannonical variables (ϕ i , π i ). In most attempts to recover time, the crucial step is the construction a quantum variableT (ϕ i , π i , U γ , F s ), which might depend on both the matter and gravitational degrees of freedom, and which shouls play the role of as a "physical clock". On the other hand, we assume we are given a wave function for the configuration variables Φ(ϕ 1 , . . . , ϕ n , U γ ). Following the canonical quantization procedure, the wave function Φ must satisfy the Hamiltonian constrain, i.e.Ĥ α Φ(ϕ 1 , . . . , ϕ n , U γ ) = 0; α = 0, 1, 2, 3, . . .. One can then obtain an effective wave function Ψ for the remaining variables by using a projection operatorP T, [t,t+δt] . The projector operator acts on the wave function Φ projecting it onto the subspace where the spectrum ofT (ϕ i , π i , U γ , F s ), corresponding to the region between t and t+δt, acquires some particular values. After an effective wave function Ψ(t) has been obtained by this procedure, one should be able, in principle, compute the expectations values of geometrical operators like areas and volumes (as well as others) to characterize the spatial metric and extrinsic curvature associated with the "hypersurface" corresponding to the selected value of t. We would do this for a large enough range of such values and obtain data characterizing the space-time in terms of data for a slicing. In other words, the procedure described allows us to obtain an "average" 3+1 decomposition for the space-time, where the slicing would be identified with the hypersurfaces on which the geometrical quantities are given by the expectation values (calculated using the wave function Ψ(t)). Therefore, one would be able to use the standard lapse and shift functions to specify the space-time and the slicing. The aspect of this approach that we want to focus on, is that, as shown in [10,11], the evolution of Ψ(t) would be given by a modified version of Schrödinger's equation including departures from the standard unitary evolution. Thus it is not unnatural that an effective sort of Schrödinger's equation, involving something that looks, from the effective point of view, as a collapse of the wave-function, could find its fundamental explanation in the setting we have proposed. In fact, these types of modifications have been previously proposed, with different motivations, and in different contexts, by the Quantum Foundations community leading to the so called collapse models [3,12,34,35,36].
As we do not have at this stage a complete workable theory of QG, the framework under which our proposal has been studied is directly that of semiclassical general relativity, which we have taken to be just an effective description of limited range of applicability, and which might well be the only setting where standard space-time notions can be recovered from an underlying quantum gravity theory. However, we will assume that its applicability includes the inflationary cosmological regime which is the situation of interest to us here. As we have shown in [7], the nature of the formal description one needs to consider, even though we do not take it to be in any sense fundamental, is not as trivial as one might initially imagine.
In that approach, one considers that the Universe can be described, by what we call a Semiclassical Self-consistent Configuration (SSC). That is, a space-time geometry characterized by a classical space-time metric and a standard quantum field theory constructed on that fixed space-time background, together with a particular state in that construction such that the semiclassical Einstein equations hold. In other words, we will say that the set g µν (x),φ(x),π(x), H | ξ ∈ H represents a SSC if and only ifφ(x),π(x) and H correspond to a quantum field theory constructed over a space-time with metric g µν (x) (as described in, say [46]), and the state |ξ in H is such that for all the points in the space-time manifold.
Such description is thought to be appropriate in the regime of interests except in those times when a collapse occurs. The resulting description is such that when the initial situation is homogeneous and isotropic, the evolution can not generate any departure from that symmetry. The emergence of inhomogeneities and anisotropies requires therefore something akin to a collapse process, something which in turn clearly necessitates a modified description. In [7] we developed what we consider to be a natural prescription for the description of such collapses, and we analyzed in detail the emergence through such process of an excitation of a single mode.
The proposal is as follows: One considers, first within the Hilbert space associated to the given SSC-i, that a transition |ξ (i) → |ζ (i) target "is about to happen", with both |ξ (i) and |ζ (i) In order to have a sensible picture, our recipe is, as discussed in [7], to relate this state |ζ (i) target with another one |ζ (ii) "living" in a new Hilbert space We denote the new SSC by SSC-ii. Thus, first we need to determine the "target" (non-physical) state in H (i) to which the initial state is in a sense "tempted" to jump, and after that, one relates such target state with a corresponding state in the Hilbert space of a new SSC, the SSC-ii. Following our previous treatments on the subject (see for instance [37]), one then considers that the target state is chosen stochastically, guided by the quantum uncertainties of designated field operators, evaluated on the initial state |ξ (i) , at the collapsing time.
It is natural to assume that after the collapse, the expectation values of the field and momentum operators in each mode will be related to the uncertainties of the pre-collapse state (recall that the expectation values in the vacuum state are zero). In the vacuum state, the operatorsŷ k andπ k are individually characterized by Gaussian wave functions centered at 0 with widths (uncertainties) (∆ŷ k ) 2 0 and (∆π k ) 2 0 , respectively. The collapse proposals that we have been considering are based on the hypothesis that each mode would jump to a new state where the expectation value would be determined by both, the scale of the uncertainties and some random variable. The idea is that the collapse would be a process taking place at a certain time, and resembling an imprecise measurement, which, of course, would be brought upon by novel physics requiring no observer or measuring apparatus.
A more precise characterization of the collapse process maintains the essence of such prescription, but now takes into account that the change of the state would need to be accompanied by a change in the space-time metric and that would require a new Hilbert space, and a new quantum state in that Hilbert space. As shown in [7], the scheme serves to characterize a complete post-collapse SSC (i.e. the SSC-ii). In fact, regarding the identification between the two different SSC's involved in the collapse, we took in [7] what turned out to be a very promising prescription: Consider that the collapse takes place along a Cauchy hypersurface Σ. A transition from the physical state |ξ (i) in H (i) to the physical state |ζ (ii) in H (ii) (associated to a certain target non-physical state |ζ (i) target in H (i) ) will occur in a way that i.e. in such a way that the expectation value of the energy momentum tensor, associated to the states |ζ (i) target and |ζ (ii) evaluated on the Cauchy hypersurface Σ, coincides. Note that the left hand side in the expression above is meant to be constructed from the elements of the SSC-i (although |ζ (i) target is not really the state of the SSC-i), while the right hand side correspond to quantities evaluated using the SSC-ii.
In the situation of interest studied in [7], the SSC-i corresponded to the homogenous and isotropic space-time characterized by Ψ = 0 with the state of the quantum field corresponding to the BD vacuum. Meanwhile, the SSC-ii corresponded to an excitation of a single mode (with comoving wave number k 0 ), characterized by a Newtonian potential given by Ψ = F (η) cos( k 0 · x), and a specific quantum state for the inflaton field. The energy momentum tensor for that state is compatible with this space-time metric according to the SSC recipe.
The point is that, as shown in detail in [7], when the SSC-ii corresponds to a potential including spatial dependences with wavenumber k 0 the normal modes of the corresponding Hilbert space H (ii) , which would otherwise be characterized by the typical spatial dependence e i k· x , and would be called the k modes, would now contain also corrective contributions of the form e i( k± k 0 )· x .
In turn, this implies that if we want to ensure that the energy momentum tensor does not lead to the excitation of other modes besides the ± k 0 in the Newtonian potential, then the state of the quantum field needs to be excited, not only in the modes ± k 0 , but also in the modes ±n k 0 , with n integer (with the degree of excitation decreasing with increasing n). In those circumstances, one can not longer argue that at the time when these higher modes collapse, there should be complete symmetry in the a priori statistical distribution of the post-collapse parameters of the state.
That is, there seems to be ample justification to assume that the value of the random parameter x n k 0 (n = 1) would not be completely independent from the value taken by the random parameters x k 0 ,y and x k 0 ,π . Focussing for simplicity on the modes, which in view of the preceding discussion one expects to be more closely correlated, we limit hereafter consideration to the case n = 2, and characterize the level of correlation by an unknown quantity we label ε. We will show in detail how these ideas are implemented in Section 4.
We should note here we will not be using at this point the full fledged formal treatment developed in [7]. This is because, as can be seen there, the problem becomes extremely cumbersome even in the treatment of a single mode. Thus, even though it is in principle possible to use such detailed formalism to treat the complete set of relevant modes, when studying the CMB spectrum the task quickly becomes a practical impossibility. We will instead rely on the less formal treatments we had employed in previous works, and retain from the results of [7] only the unusual correlations we discussed above. Before doing so, we will briefly review the basic details, of the simple version of our formalism. The reader should be warned that this is not (and can not expected to be) a fully self contained work, as doing that would involve repeating extensive discussions and analysis already published in the literature.

The collapse model for the inf lationary f ield
The starting point of the specific analysis is the same as the standard picture, i.e. the action of a scalar field coupled to gravity: where φ stands for the inflaton and V for the inflaton's potential. One then splits both, metric and scalar field into a spatially homogeneous (background) part and an in-homogeneous part (fluctuation), i.e. g = g 0 + δg, φ = φ 0 + δφ. Note that in the full fledged formalism introduced in [7], the quantum description is applied, not only to the inflaton mode perturbation, but also to the "background field" which is characterized in terms of the zero mode of a quantum field theory construction on the inflating space-time. At the same time, the metric sector is treated exclusively at the classical level as corresponding to the semiclassical description we have relied on (for more on the justification for this see [7]). Here however, we will return to the simpler but more convenient formalism of previous treatments, retaining as we have said only the interesting novel correlations we discussed. Thus, here the zero mode of the inflaton field will be treated as a classical background field.
The background is taken to be the spatially flat Friedmann-Robertson Universe with line element ds 2 = a(η) 2 −dη 2 + δ ij dx i dx j , and the homogeneous scalar field φ 0 (η). The evolution equations for this background are scalar field equations, The scale factor solution corresponding to the inflationary era of standard inflationary cosmology, written using a conformal time, is: is the slow-roll parameter which during the inflationary stage is considered to be very small 1, thus H I const, and with the scalar field φ 0 in the slow roll regime, i.e.φ 0 = −(a 3 /3ȧ)V . According to the standard inflationary scenario, this era is followed by a reheating period in which the Universe is repopulated with ordinary matter fields, a regime that then evolves towards a standard hot big bang cosmology regime leading up to the present cosmological time. The functional form of a(η) during these latter periods changes, but we will ignore those details because most of the change in the value of a occurs during the inflationary regime. We will set a = 1 at the "present cosmological time", and assume that inflationary regime ends at a value of η = η R , negative and very small in absolute terms (η R −10 −22 Mpc).
Next, we consider the perturbations. We will ignore vectorial perturbations, and focus on the scalar and tensorial ones. First, we will focus on the scalar perturbations. Working in the so called longitudinal gauge, the perturbed metric is written as: where Ψ stands for the scalar perturbation usually known as the Newtonian potential.
The perturbation of the scalar field is related to a perturbation of the energy momentum tensor, and reflected into Einstein's equations which at lowest order lead to the following constraint equation for the Newtonian potential: where we introduced the abbreviation s ≡ 4πGφ 0 . Regarding the field δφ, it is convenient to work with the rescaled field variable y = aδφ and its conjugate momentum π =ẏ −ȧ a y = aδφ. We consider the quantum theory of this field setting the problem, for simplicity, in a finite box of side L, which can be taken to ∞ at the end of all calculations. Thus we write the expansions for field and conjugate momentum operators aŝ where the sum is over the wave vectors k satisfying k i L = 2πn i for i = 1, 2, 3 with n i integer and whereŷ k (η) ≡ y k (η)â k + y * with the usual choice of modes: which leads to what is known as the Bunch-Davies vacuum.
Note that according to the point of view we discussed in Section 2.1, and having at this point the quantum theory for the relevant matter fields, the effects of the quantum fields on the geometrical variables are codified in the semiclassical Einstein's equations. Thus equation (1) must be replaced by At this point, one can clearly observe that if the state of the quantum field is in the vacuum state, the metric perturbations vanish and thus the space-time is homogeneous and isotropic. The proposal is based on the consideration of a self induced collapse, which we take to operate in close analogy with a "measurement" (but evidently, with no external measuring apparatus or observer involved). This leads us to want to work with Hermitian operators, which in ordinary quantum mechanics are the ones susceptible of direct measurement. Therefore, we must separate bothŷ k (η) andπ k (η) into their real and imaginary partsŷ k (η) =ŷ k So far, the treatment is similar to the standard one, except in that only the scalar field is being treated at the quantum level, while the metric fluctuation is treated classically. At this point, it is worthwhile to emphasize that the vacuum state defined byâ k R,I |0 = 0 is 100% translational and rotationally invariant, as can be seen from the fact that it is annihilated by the generators of translations and rotations (i.e. by the linear and the angular momentum operators, P andˆ L, respectively).
Next we need to specify in more detail the modeling of the collapse. Then, take into account that after the collapse has taken place, we must consider the continuing evolution of the expectation values of the field variables until the end of inflation, and eventually up to the last scattering hypersurface (in fact, if we want to actually compare our analysis with the detailed observations, we must evolve also through the reheating period and through the decoupling era up to today's Universe. This, however, is normally taken into account through the use of appropriate transfer functions, and we will assume that the same procedure can be implemented after the present analysis).
We will further assume that the collapse is somehow analogous to an imprecise measurement 3 of the operatorsŷ R,I k (η) andπ R,I k (η). Now, we will specify the rules according to which the 2 Note however that he commutation relations of the real and imaginary creation and annihilation operators are however non-standard: , with all other commutators vanishing. This is known to indicate that the operators corresponding to k and − k are identical in the real case and identical up to a sign in the imaginary case. 3 An imprecise measurement of an observable, is one in which one does not end with an exact eigenstate of that observable but rather with a state which is only peaked around the eigenvalue. Thus, we could consider measuring a certain particle's position and momentum, ending up with a state characterized by wave packet with both position and momentum defined to a limited extent, and which, of course, does not entail a conflict with Heisenberg's uncertainty bound. collapse happens. Again, at this point our criteria will be simplicity and naturalness. What we have to describe is the state |Θ after the collapse. It turns out that, for the goals at hand, all we need to specify is the expectation values immediately after the collapse as this determines the expectation value of the field and momentum operator for the mode k at all times after the collapse. That is using the result in the evolution equations for the expectation values (i.e. using Ehrenfest's theorem), one obtains ŷ k R,I (η) and π k R,I (η) for the state that resulted from the collapse, for all later times.
As we indicated, at this point we assume that after the collapse, the expectation values of the field and momentum operators in each mode will be related to the uncertainties of the precollapse state (recall that the expectation values in the vacuum state are zero). In the vacuum state,ŷ k andπ k individually are distributed according to Gaussian wave functions centered at 0 with spread (∆ŷ k ) 2 0 and (∆π k ) 2 0 , respectively. We might consider various possibilities for the detailed form of this collapse. Thus, for their generic form, associated with the ideas above, we write: where η c k represents the time of collapse for each mode. The random variables x R,I k,y , x R,I k,π are selected randomly from a distribution which, in principle, might be non-Gaussian, i.e. different modes could present particular correlations (we will elaborate on this in the next section). At this point, we must emphasize that our Universe corresponds to a single realization of these random variables, and thus each of these quantities x R,I k,y , x R,I k,π have a single specific value. The two parameters λ y and λ π are real numbers we use to characterize the different collapse schemes, e.g.: i) λ y = 0, λ π = 1, ii) λ y = λ π = 1 (which we call the "symmetric scheme"). It is clear that one can devise many other models of collapse, a good fraction of them can be described within the scheme above, while others require a slightly modified treatment [5]. Still, there are surely many other possibilities which we have not even thought about and which might require drastically modified formalisms.
The explicit expressions for the ŷ R,I k (η) Θ and π R,I k (η) Θ are then where D k (η) ≡ kη − z k and z k ≡ kη c k . With this information at hand we can now compute the perturbations of the metric after the collapse of all the modes 4 .

The observable quantities
Now, we must put together our semi-classical description of the gravitational degrees of freedom and the quantum mechanics description of the inflaton field. We recall that this entails the semiclassical version of the perturbed Einstein's equation that, in our case, leads to equation (2).
The Fourier components at the conformal time η are given by: where we used that during inflation s ≡ 4πGφ 0 = /2aH I /M P , with M P the reduced Planck's mass M 2 P ≡ 2 /(8πG). The expectation value depends on the state of the quantum field, therefore, as we have already noted, prior to the collapse, we have Ψ k (η) = 0, and the space-time is still homogeneous and isotropic at the corresponding scale. However, after the collapse takes place, the state of the field is a different state with new expectation values which generically will not vanish, indicating that after this time the Universe becomes anisotropic and in-homogeneous at the corresponding scale. We now can reconstruct the space-time value of the Newtonian potential using Ψ(η, x) = 1 L 3 k e i k· x Ψ k (η), to extract the quantities of observational interest. In order to make contact with the observations, we shall relate the expression (7) for the evolution of the Newtonian potential during the early phase of accelerated expansion, to the small anisotropies observed in the temperature of the cosmic microwave background radiation, δT (θ, ϕ)/T 0 with T 0 ≈ 2.725K the temperature average. They are considered as the fingerprints of the small perturbations pervading the Universe at the time of decoupling, and undoubtedly any model for the origin of the seeds of cosmic structure should account for them.
The observational data is described in terms of the coefficients α lm of the multipolar series expansion Here θ and ϕ are the coordinates on the celestial two-sphere, with Y lm (θ, ϕ) the spherical harmonics. The different multipole numbers l correspond to different angular scales; low l to large scales and high l to small scales. At large angular scales (l 20) the anisotropies in the CMB arise due to the Sachs-Wolfe effect. That effect relates the anisotropies in the temperature observed today on the celestial sphere to the inhomogeneities in the Newtonian potential on the last scattering surface: δT T 0 (θ, ϕ) = 1 3 Ψ(η D , x D ), where η D is the conformal time of decoupling which lies in the matter-dominated epoch, and x D = R D (sin θ sin ϕ, sin θ cos ϕ, cos θ), with R D the radius of the last scattering surface. Furthermore, using Ψ(η, x) = 1 L 3 k e i k· x Ψ k (η) and e i k· x D = 4π lm i l j l (kR D )Y lm (θ, ϕ)Y * lm (k), the expression (8) for α lm can be rewritten in the form with j l (kR D ) the spherical Bessel function of order l. The transfer function ∆(k) represents the evolution of the Newtonian potential from the end of inflation η R to the last scattering surface at the time of decoupling η D , i.e. Ψ k (η D ) = ∆(k)Ψ k (η R ). We will be ignoring this aspect from this point onward, despite the fact that this transfer functions are behind the famous acoustic peaks, the most noteworthy feature of the CMB spectrum, as they relate to aspects of plasma physics that are well understood and thus uninteresting from our point of view. This will mean that, from our point of view, the observational spectrum would have such features removed before comparing with our results (this is in the same spirit that one removes the uninteresting imprint of our galaxy on the observations of the CMB, or the dipole associated with our peculiar motion). This would be equivalent to assume that the observed spectrum fits well with the flat Harrison-Zel'dovich spectrum and setting ∆(k) = 1, and thus setting for our purposes here Substituting (6) in (7) and using (3), (4) gives where z k ≡ kη c k , we note that the quantity D k is being evaluated at η R , i.e. D k (η R ) = kη R − z k . Finally using, (10) in (9) yields: We will further simplify the treatment by noting that the value of −kη R is exceedingly small for the modes of interest. Thus, this will be ignored in comparison to z k , i.e. we will be assuming that we can make the approximation D k → −z k . It is worthwhile to mention that the relation of α lm with the Newtonian potential, as obtained in (11) within the collapse framework has no analogue in the usual treatments of the subject. It provides us with a clear identification of the aspects of the analysis where the "randomness" is located. In this case, it resides in the randomly selected values x R,I k,y , x R,I k,π that appear in the expressions of the collapses associated with each of the modes. Here, we also find a clarification of how is it that, in spite of the intrinsic randomness, we can make any prediction at all. The individual complex quantities α lm correspond to large sums of complex contributions, each one having a certain randomness, but leading in combination to a characteristic value, in just the same way as a random walk that is made up of multiple steps. In other words, the justification for the use of statistics in our approach is: The quantity α lm is the sum of contributions from the collection of modes, each contribution being a random number leading to what is in effect a sort of "two-dimensional random walk", which total displacement corresponds to the observational quantity. Nothing like this can be found in the most popular accounts, in which the issues we have been focussing on are hidden in a maze of often unspecified assumptions and unjustified identifications [44].
Finally, let us discuss the implications of our proposal regarding the tensorial perturbations. The analysis we present is for the sake of completeness, as we have already presented it in previous works [24,37], and is not the main subject of this article.
The view we take, regarding the relation between the metric (classical) description of gravity and the quantum matter degrees of freedom, has a direct implication on the estimates for the production of gravitational waves (from inflation). The idea is that the mechanism that generates the actual the fluctuations is tied to on the quantum uncertainties for the scalar field. Due to some unknown quantum gravitational effect, the state of each mode of the scalar field collapses. After the collapse, the density inhomogeneities and anisotropies feed into the gravitational degrees of freedom resulting in perturbations (scalar, vector and tensor) of the metric. Nevertheless, the metric itself does not induce a quantum gravitational collapse. Therefore, as the scalar field does not act as a source for the gravitational tensor modes -at least not at the order considered here -the tensor modes cannot be excited. Thus, the collapse scheme naturally predicts 5 a zero -or at least a strongly suppressed -amplitude of gravitational waves into the CMB.

Modif ied statistical features in the primordial spectrum
In this section, we will show how the statistical aspects of the collapse model might induce non-trivial features in the primordial spectrum. This is a novel aspect and contrasts with what has been done in the previous approaches, where the different modes that contribute to the primordial spectrum are considered to be uncorrelated.
Let us start by noting that according to (11) all the modes contribute to α lm with a complex number. If we had the outcomes characterizing each of the individual collapses we could, of course, predict the value of these quantities. However, we have at this point no other access to such information than the observational quantities α lm themselves.
We hope to be able to say something about these, but doing so requires the consideration of further hypothesis regarding the statistical aspects of the physics behind the collapse as well as the conditions previous to them.
As is generally the case with random walks, one can not hope to estimate the direction of the final displacement. However, one might say something about its estimated magnitude. It is for that reason we will be focussing on estimating the most likely value of the magnitude |α lm | 2 . We can make some progress, for instance, by making some assumption allowing us to regard the specific outcomes characterizing our Universe as a typical member of some hypothetical ensemble of Universes.
For example, we are interested in estimating the most likely value of the magnitude of |α lm | 2 above, and in such hypothetical ensemble we might hope that it comes very close to our single sample. It is worthwhile emphasizing that, for each l and m, we have one single complex number characterizing the actual observations (and thus the real Universe we inhabit). For a given l for instance, we should avoid confusing ensemble averages with averages of such quantities over the 2l + 1 values of m. The other universes in the ensemble are just constructs of our imagination and there is nothing in our theories that would indicate that they are real.
We can simplify things even further by taking the ensemble average |α lm | 2 (recall that the bar represents an ensemble average) and identifying it with the most likely value of the quantity, and needless is to say that these two notions are different in many types of ensembles. However, let us, for the moment, ignore this issue and assume the identity of those two values. We thus have are completely uncorrelated (leading thus to the cancelation of the cross terms in (12)). As mentioned in Section 2.1, our refined approach suggests that, within each set of random variables, there would be a correlation between the modes k and 2 k (and also of course between k quantum interface faced in the current problem. It is of course conceivable, although seems harder to understand in a wider context (see the discussion in Section 8 of [43]), that a collapse might be incorporated into a setting where both the gravitation and scalar field perturbations are simultaneously treated at the quantum level. If the latter happened to be the correct approach, something that would be possible to ascertain when we have a fully satisfactory theory of quantum gravity, our conclusion about the tensor modes would be modified. and k/2). This correlation would induce a "small" departure from the normal probability distribution function of the random variables, where the "smallness" is characterized by the introduction of a parameter ε (do not confuse ε with the slow-roll parameter ). That is, we would assume that the ensemble average of the product of two random variables is given by x I k,i with i = y, π. Note that this suggests also many more complex and higher order correlations between various sets of variables, but we will neglect those higher order effects here, in part because we do not have a specific self-content mold for those at this point. We recall that the random variables corresponding to k and − k are not independent. Using (13) and (14) in (12) we obtain As we have found in previous works, the data indicates that, to a good degree of accuracy [5], the assumption that z k is independent of k, i.e. the time of collapse goes as η c k ∝ k −1 , is a good approximation. Therefore, assuming z k = z = const, (15) takes the form where C(z) ≡ λ 2 π M 2 (z) + λ 2 y N 2 (z). Now we write the sum as an integral by noting that the allowed values of the components of k are separated by ∆k i = 2π/L, that is where in the previous expression we performed the integral over the angular part of k.
Rewriting the remaining integrals, we find where we use a change of variable x ≡ kR D . Evaluating the integrals of the products of spherical Bessel's functions gives the final result where Γ(l) is the Gamma function, and 2 F 1 (α, β; γ, z) is the hypergeometric function 6 . 6 The hypergeometric function is defined as 2F1(α, β; γ, z) ≡ Now, as (16) does not depend on m it is clear that the expectation of the quantity of observational interest, i.e. the angular spectrum l(l + 1)C l ≡ l(l + 1)(2l + 1) −1 m |α lm | 2 , is just |α lm | 2 and thus this quantity must be compared with: where we have defined are uncorrelated, random, and Gaussian (see (13) and (14)), then we recover the standard flat spectrum, i.e. l(l + 1)C l is independent of l. The departure of the flat spectrum has a very specific signature that, in principle, can be searched for observationally. There are of course many other possibilities of unexpected correlations and each one will have a particular signature.
The behavior of G(l) for 2 ≤ l ≤ 20 is shown in Fig. 1. As we can see, the effect is stronger for low l and it decreases in a nearly exponential fashion for large l (we have made a direct confirmation of this fact). The next step, of a more rigorous analysis, would be to consider the acoustic oscillations, due to plasma physics, which play an important role after the end of the inflationary regime. Then, we could use the recent observational data available (e.g. WMAP7), and thus put some constraints on the parameter ε.
Recall, however, that the result (17) we obtained depends on the behavior of another parameter we have introduced, namely the time of collapse η c k . Let us recall, that in order to recover the detailed form of the spectrum, we have had to make the additional assumption that η c k ∝ 1/k (i.e. z k is independent of k). At this point we only have a very heuristic proposal that seems to account for this assumption (see [37], in terms of an argument that attempts to connect that pattern with certain ideas of Penrose's and Diosi's linking the collapse of the wave function with aspects of quantum gravity). In fact, the assumption that η c k ∝ 1/k seems necessary to recover something close to the observed spectral shape, as has already been noted in previous works [5,37]. The point is that, in order to constrain the values of ε and z k , we should perform an analysis such as the one presented in [22]. That work, consisted of a detailed statistical analysis testing some specific "collapse schemes", and establishing constraints on the free parameters of the modes considered. The analysis was carried out comparing the details predictions obtained by a modified CMB fast code incorporating out models, and recent data from the CMB, including the 7-yr data release of the WMAP [23] and the matter power spectrum measured employing the Luminous Red Galaxies by the Sloan Digital Sky Survey [38]. A similar analysis for the results of the present paper lies outside the scope of this work. However, our main result of this work, i.e. the pattern noted in (17), indicates that the appearance of the correlations between the modes k and 2 k; affecting the statistics of the random variables , and thus our estimates for the observational quantities. It seems that looking for such specific pattern of deviation in the form of CMB spectrum, would be a straight forward (but highly resource consuming) computational endeavor.
Finally, let us end this section by noting a heuristic relation between our result (17) and the so called "parametric resonance" effect in quantum optics. The set up involves a parametric nonlinear crystal excited by an intense laser beam at angular frequency ω p = 2ω. The interesting fact is that a weak signal beam at angular frequency ω is spontaneously generated. In a particular setting, the nonlinear crystal can mix a second input signal ω s = ω with that resulting from the pump, and this produces an idler beam at angular frequency ω. The nonlinear process either amplifies or de-amplifies the excitation of the signal modes depending on the phases. Thus, the output field is the result of delicate quantum interference resulting from a highly correlated state. In short, the nonlinear nature of the crystal's response leads to photons of one frequency being converted into photons with a different one, leading to the appearance in the spectrum of the beam of important connections between the components with a given frequency and its first harmonics. The connection with our result (17) would lead us to view the gravitational response to the collapse of matter fields, as described by the SSC scheme, playing the role of the "non-linear crystal". That is, the collapse seems to correlate the mode k the modes n k, (n an integer), in the same sense as the crystal correlates the angular frequency of the signal and the pump. The term proportional to ε in (17) would correspond to the modification on the amplitude in the previous example.

Discussion
The work presented here is a continuation of a long standing program aimed at addressing and proposing a possible cure for a set of shortcomings that afflict the standard accounts for the emergence of the seeds of cosmic structure during inflation. The program has revealed that, although the motivating problems can be characterized as conceptual difficulties, it is intimately connected with what is usually called the measurement problem in quantum mechanics. This problem, which as we have argued, becomes essentially insurmountable within the standard quantum theory in the cosmological context (see [44]). Nevertheless, one is led surprisingly to expect observational consequences associated with some details of the physical unknown mechanisms that lie beyond quantum theory and our current understanding. That is, by parameterizing some aspects of this enigmatic process, we might actually use observations as a guides in learning something about them. In previous works, we have seen, for instance, that the details of the "collapse schemes" can modified the primordial fluctuation spectral shape in very particular ways [37]. We have argued that we can use the observations to place limits on the times of collapse of the various modes [5].
We have discussed that the semiclassical setting should be the appropriate framework to deal with inflationary perturbations. This is simply because the full quantum theory is likely to be one where the space-time notions are not present at all (due in part by the problem of time that afflicts canonical approaches to quantum gravity such as LQG). As has been argued elsewhere in that setting, there would not be tensor modes exited at first order in perturbation theory and thus no primordial gravity waves of a magnitude that could be observed in the forceable future.
Moreover, in the approach we have been using in this and previous works, it becomes quite transparent on how the statistical aspects are related to each one of the observational quantities. That is, the observable α lm is related to the random variables x k 's through a complex "random walk", while the value of |α lm | 2 corresponds to the squared "length" of the random walk. This random walk is associated to a particular realization of a physical quantum process (i.e. associated with one of the multiple possibilities characterized by the collapse of the inflaton's wave-function). However, as we have access to just one such realization -the set of random walks corresponding to our own Universe -the most natural course of action (but certainly not the only one) is to take the ensemble average value of the length of the possible random walks, which corresponds to |α lm | 2 , to be identified with the most likely value, i.e. with |α lm | 2 ML ; and this in turn, is to be compared with the |α lm | 2 of our observable Universe. Of course, our Universe is to be identified with one realization of the random variables, and as such, the resulting picture has an interpretational transparency that seems hard to find in the standard accounts.
The clarification mentioned above has allowed us to take a fresh approach in addressing the possibility of novel statistical features and the manner in which those could be searched for, as discussed in [25]. In the present work we have discussed a novel possibility for nonstandard statistics uncovered by our approach. This novel aspect is tied to the possibility that non trivial statistical features would arise as a result of the of the collapse process itself. Namely that the collapse of one of the modes (say the mode k) might influence the tendency in which the state of a related mode (say the mode 2 k) to collapse in a certain manner. We can heuristically view this possibility as reflecting a modification of the statistical aspects of the state of one mode as the result of the collapse of another.
We noted that the deviation from the standard prediction would be particularly relevant for the first multipoles: It is thus interesting to speculate whether this effect could be linked to the issues raised in [4,39]. In those works, it has been argued that there are specific large-scale anomalies in the temperature map of anisotropies of the CMB for the low modes. According to [39], the observations disagree markedly with the predictions of the standard theory. By that as it might, it seems, that our result (17), which was motivated by purely conceptual issues, shares some of the essential features of those observations, i.e. a departure of the standardflat -spectrum for the lowest l. Evidently, at this stage of the analysis we cannot present more quantitative results, however, we believe that the previous discussion serves to show the potential behind our proposal.
In any event, we have seen that such kind of modifications have particular signatures which could, in principle, be searched for in the data, leading as we have said, either to interesting discoveries, or at least to the setting of bounds on the type of effect we have been considering.