Quasi-matter bounce and inflation in the light of the CSL model

The Continuous Spontaneous Localization (CSL) model has been proposed as a possible solution to the quantum measurement problem by modifying the Schr\"{o}dinger equation. In this work, we apply the CSL model to two cosmological models of the early Universe: the matter bounce scenario and slow roll inflation. In particular, we focus on the generation of the classical primordial inhomogeneities and anisotropies that arise from the dynamical evolution, provided by the CSL mechanism, of the quantum state associated to the quantum fields. In each case, we obtained a prediction for the shape and the parameters characterizing the primordial spectra (scalar and tensor), i.e. the amplitude, the spectral index and the tensor-to-scalar ratio. We found that there exist CSL parameter values, allowed by other non-cosmological experiments, for which our predictions for the angular power spectrum of the CMB temperature anisotropy are consistent with the best fit canonical model to the latest data released by the Planck Collaboration.


I. INTRODUCTION
After approximately three decades since the cosmological inflationary paradigm was conceived [1][2][3][4], all of its generic predictions have withstood the confrontation with observational data, in particular, those coming from the Cosmic Microwave Background (CMB) radiation [5][6][7]. That has led a large group of cosmologists to consider inflation as a well established theory of the early Universe. Inflation was originally proposed to provide a solution to the puzzles of the hot Big Bang theory (e.g. the horizon and flatness problems). However, the modern success of inflation is that, allegedly, it can offer us an explanation about the origin of the primordial inhomogeneities [8][9][10][11]. The standard argument is also rather pictorial: the quantum fluctuations of the vacuum associated to the inflaton are stretched out to cosmological scales due to the accelerated expansion of the spacetime; those fluctuations are considered the seeds of all large scale structures observed in the Universe. Furthermore, in Ref. [12] is investigated the detectability of possible traces of the quantum nature regarding the primordial perturbations.
On the other hand, proponents of alternative scenarios to inflation argue that even if it is the most fashionable model of the early Universe, that does not mean it is necessarily true [13,14]. Furthermore, another feature that would make alternative models worthwhile of study is that they might avoid some long standing puzzles of the inflationary paradigm. Among those issues, we can mention: the subject of eternal inflation, a feature that is present in almost every model of inflation [15] and which also leads to the controversial topic of the multiverse; the initial singularity problem and the trans-Planckian problem for primordial perturbations [16], which are related by the fact that one is interpolating the solutions provided by General Relativity in regimes where it may no longer be valid; and finally, it has been argued that the potentials associated to the inflaton, that best fit the latest observed data, need to be fine-tuned [17,18]. Although the aforementioned problems are not considered real problems by some scientists [19,20], others seem to disagree [18,21]. However, we think that if other alternative models can reproduce the main results linked to inflation, we should make use of the observational data available to test them.
One of the alternative models to inflation that seems to be consistent with the latest data is the so called matter bounce scenario (MBS) [21][22][23][24][25][26][27][28]. In this cosmological model, the initial singularity of the standard model is replaced by a non-singular bounce. That is, instead of an ever-expanding Universe, it assumes an early contracting matterdominated Universe, which continues to evolve towards a bouncing phase and, later, enters into the expanding-phase of standard cosmology. The Universe described by the MBS relies on a single scalar field satisfying an equation of state that mimics that of a dust-like fluid. Additionally, in order to describe successfully a bouncing phase with a single scalar field, one needs to use cosmologies beyond the realm of General Relativity, such as, loop quantum cosmologies, teleparallel F (T ) gravity or F (R) gravity. Proponents of the MBS claim that the potential associated to the scalar field is less fine-tuned than that of inflation, and also solves the historical problem of requiring very special initial conditions for the Big Bang model [22,23], which originally motivated the development of the inflationary framework. However, the MBS is not exactly problem free. A complete assessment of the present conceptual issues is provided in Ref. [23]. In spite of not being completely finished from a theoretical point of view, the MBS is quite simple in its treatment of the primordial perturbations. That makes it an interesting case of study for the purpose of this article. In particular, the generation of the primordial perturbations is depicted during the contracting phase, i.e. in a regime where gravity is well described by General Relativity, and the perturbations correspond to inhomogeneities of a single scalar field.
In addition to the prior puzzles and successes of inflation and the MBS, there remains an important question: what is the precise physical mechanism that converts quantum fluctuations of the vacuum into classical perturbations of the spacetime? This question has been the subject of numerous works in the past and the consensus seems to favor the decoherence framework [29][30][31][32][33]. 1 Nevertheless, decoherence cannot address that question by itself. 2 In other words, even if one would choose (or not) to embrace the decoherence program, a particular interpretation of Quantum Mechanics must be selected (implicitly or explicitly). The Copenhagen-orthodox-interpretation requires to identify a notion of observer that performs a measurement on the system; which, in the decoherence framework, is equivalent to identify the unobservables or external degrees of freedom of the system. It is not clear how to do such identifications if the system is the early Universe. Other interpretations such as many-worlds, consistent histories and hidden variables formulations, might be adopted with varying degrees of success (see for instance [38][39][40]).
In the present article, in order to address the quantum-to-classical transition of the primordial perturbations, we will choose to work with the Continuous Spontaneous Localization (CSL) model. The CSL model belongs to a large class of models known as objective dynamical reduction models or simply called collapse models. Collapse models attempt to provide a solution to the measurement problem of Quantum Mechanics [41][42][43][44][45]. The proponents of these models state that the measurement problem originates from the linear character of the quantum dynamics encoded in the Schrödinger equation. The common idea shared in these collapse models is to introduce some nonlinear stochastic corrections to the Schrödinger equation that breaks its linearity. According to the collapse models, a noise field couples nonlinearly with the system (usually with the spatial degree of freedom of the system), inducing a spontaneous random localization of the wave function in a sufficiently small region of the space. Suitably chosen collapse parameters make sure that micro-systems evolve essentially (but not exactly) following the dynamics provided by the Schrödinger equation, while macro-systems are extremely sensible to the nonlinear effects resulting in a perfectly localization of the wave function. Furthermore, there is no need to mention or to introduce a notion of an observer/measurement device as in the Copenhagen interpretation, which is a desired feature in the context of the early Universe and cosmology in general.
The CSL model has been applied before to the inflationary Universe in an attempt to explain the quantumto-classical transition of the primordial perturbations [46][47][48][49][50]. Also, recently a new effective collapse mechanism, independent of the CSL model, has been proposed to deal with the measurement problem during the inflationary era [51]. However, among those works, the key role played by the collapse mechanism varies and also yields different predictions for the primordial power spectrum, some of which might be consistent with the observational data. In the present article, we will subscribe to the conceptual point of view first presented in [47,52], which was developed within the semiclassical gravity framework, and latter in [50,53] was extended to the standard quantization procedure of the primordial perturbations using the so called Mukhanov-Sasaki variable [8,54]. The main role that we advocate for the dynamical reduction mechanism of the state vector, modeled in this paper by the CSL model, is to directly generate the primordial curvature perturbations. Specifically, the initial state of the quantum field-the vacuum state-evolves dynamically according to the modified Schrödinger equation provided by the CSL model. This evolution leads to a final state that does not share the initial symmetry of the vacuum, i.e. it is not homogeneous and isotropic. 3 In this way, the collapse mechanism generates the inhomogeneities and anisotropies of the matter fields. These asymmetries are codified in the evolved quantum state and, thus, are responsible for generating the perturbations of the spacetime. 4 Note that the previous prescription, regarding our approach to address the birth of the primordial perturbations, does not require the inclusion of an exponential expansion phase in the Universe that "stretches out" the quantum fluctuations of the vacuum (or the squeezing of the field variables as usually argued). Therefore, in principle, it should be possible to extend our picture to alternative scenarios dealing with the origin of the cosmological perturbations. Moreover, since the cosmic observations are well constrained, it should also be feasible to test the predictions that result from applying our framework in those alternative cosmological models. In the present article, we focus on the implementation of the CSL model within the framework of the MBS and, in parallel, we present the same appliance of the CSL model to the slow roll inflationary model of the early Universe. In this way, we can appreciate more clearly where the CSL model enters into the picture; particularly at the moment when computing the theoretical predictions. The main motivation behind the present work is that if the CSL model can be truly considered as a physical model of the quantum world, which also avoids the standard quantum measurement problem, then it should also be possible to use it in different contexts from the traditional laboratory settings. The cosmological context provides a rich avenue to explore such foundational issues and, more important, there exist sufficient precise data to test the initial hypotheses. As a consequence, we will analyze the predictions resulting from implementing the CSL model in the MBS and in the inflationary model of the early Universe, and we will compare the corresponding results with the one provided by the best fit standard cosmological model. Additionally, we will focus on the range of values allowed for the parameters of the CSL model, experimentally tested [55,56] in non-cosmological frameworks.
The paper is organized as follows: in Sect. II, we provide a very brief synopsis of the main features of the CSL model, with particular emphasis on those that will be useful for the next sections. In Sect. III, we present the characterization of the primordial perturbations within the two cosmological models that we are considering, i.e. the MBS and standard slow roll inflation. In Sect. IV, we show the connection between the observational quantities and the theoretical predictions that result from adopting our conceptual point of view concerning the CSL model. In Sect. V, we explicitly show the implementation of the CSL model to the MBS and inflation, and we also present the predictions for the primordial power spectra (scalar and tensor) in each case. In Sect. VI, we discuss the implications of the results obtained; additionally, we compare the predicted scalar power spectra with the standard one. In Sect. VII, we analyze the viability of the CSL model using the data extracted from the CMB when considering the best fit cosmological model. Finally, in Sect. VIII, we end with our conclusions. We include an Appendix containing the computational details that led to the results presented in Sect. V.

II. A CONCISE SYNOPSIS OF THE CSL MODEL
In this section, we provide a brief summary of the relevant features of the CSL model; for a detailed review, we refer the reader to Refs. [43,44].
In the CSL model, the modification of the Schrödinger equation induces a collapse of the wave function towards one of the possible eigenstates of an operatorΘ, called the collapse operator, with certain rate λ. The self-induced collapse is due to the interaction of the system with a background noise W(t) that can be considered as a continuous-time stochastic process of the Wiener kind. The modified Schrödinger equation drives the time evolution of an initial state as withT the time-ordering operator. The probability associated with a particular realization of W(t) is, The norm of the state |Ψ, t evolves dynamically, and Eq. (2) implies that the most probable state will be the one with the largest norm. From Eqs. (1) and (2), it can be derived the evolution equation of the density matrix operator ρ. That is, The density matrix operator can be used to obtain the ensemble average of the expectation value of an operator Ô = Tr [Ôρ]. Henceforth, from Eq. (3) it follows that The average is over possible realizations of the noise W(t), each realization corresponding to a single outcome of the final state |Ψ, t . One of the most important features of collapse models is the so-called amplification mechanism. That is, assuming that the reduction (collapse) rates for the M constituents of a macroscopic object are equal (λ i = λ), it can be proved that the reduction rate for the center of mass of an M -particle system is amplified by a factor of M with respect to that of a single constituent [41,57]; in other words, λ macro = M λ.
The parameter λ sets the strength of the collapse process. In the original model, proposed by Ghirardi-Rimmini-Webber (GRW), the authors suggested a value of λ GRW ≃ 10 −16 s −1 for r C ≃ 100 nm. However, Adler suggested a greater value λ Adler ≃ 10 −8 s −1 for r C ≃ 100 nm [58] (the parameter r C is called the correlation length of the noise and provides a measure for the spatial resolution of the collapse [41,43,57]). Recent experiments have been devised to set bounds on the parameter λ [59,60]. Furthermore, it is claimed that matter-wave interferometry provides the most generic way to experimentally test the collapse models [55,56]. Those results suggest that the range between λ GRW and λ Adler is still viable for some variations of the original CSL model (e.g. by considering non-white noise).
Henceforth, the main characteristics of the CSL model are: (1) The modification to the Schrödinger equation is nonlinear and leads to a breakdown of the superposition principle for macroscopic objects; (2) The random nature of Quantum Mechanics is concealed in the noise W(t) and is consistent with Born's rule; (3) An amplification mechanism exists, through the parameter λ which is related to the strength of the collapse. This strength is weak for microscopic objects and strong for macroscopic bodies.
Another main aspect of the collapse models is that the collapse mechanism injects energy into the system. In fact, previous works have performed a preliminary analysis using cosmological data to set bounds on the value of λ [61]. The energy increase is minimal, e.g. for a particle of mass m = 10 −23 g, one obtains δE/t ≃ 10 −25 eV s −1 [43]. In other words, an increase of 10 −8 eV will take 10 10 years. However, even if the energy increase can be ignored at the phenomenological level, a more realistic model should remove this issue.
Moreover, the increase of energy in the collapse models leads to difficulties when trying to formulate relativistic collapse models. Additionally, the collapse mechanism occurs in such a way that is nonlocal. This implies that the collapse of the wave function must be instantaneous or superluminal (but the nonlocal features cannot be exploited to send signals at superluminal speed). Also, the nonlocality is necessary to ensure that the models are consistent with the violation of Bell's inequalities. Several relativistic models have been proposed so far [62][63][64], none of which can be considered completely finished. In spite of the lack of a relativistic collapse model, we will apply the CSL model to the primordial Universe, i.e. to inflation and the MBS, but in order to provide a more detailed picture, we need first to establish the mathematical framework of the primordial Universe in the two approaches considered in this work.

III. TWO APPROACHES: ACCELERATED EXPANSION OR QUASI-MATTER CONTRACTION
This section presents the details of the two cosmological approaches, describing the dynamics of the Universe, that we will be considering in the rest of the manuscript. In particular, we are going to work with the following two scenarios: 1. An accelerated expansion of the early Universe given by the simplest inflationary model, that is, a single scalar field in the slow roll approximation with canonical kinetic term. Since such a model is probably very well known for most readers, we will not dwell into much detail here.
2. The MBS [22][23][24][25][26], a cosmological model in which the Universe undertakes a quasi-matter contracting phase, then experiences a non-singular bounce and finally enters into the standard cosmological expansion. Since in this model the primordial perturbations are born during the contracting stage of the Universe, we will focus exclusively on that cosmic stage. We will refer to such a stage as the quasi-matter contracting Universe (QMCU).

A. The background
The inflationary Universe and the QMCU are both described by Einstein equations G ab = 8πGT ab (c = 1), while the matter fields are characterized by a single scalar field. In the case of inflation the scalar field is the inflaton φ, and in the QMCU the scalar field will be denoted by ϕ.
As mentioned earlier, for inflation, we will consider standard slow roll inflation. In that case, the background spacetime is described by a quasi-de Sitter Universe, characterized by H ≃ −1/[η(1 − ǫ)], with H ≡ a ′ /a the conformal expansion rate, a being the scale factor and the slow roll parameter is defined as ǫ ≡ 1−H ′ /H 2 ; a prime denotes partial derivative with respect to conformal time η. The energy density of the Universe is dominated by the potential of the inflaton V , and during slow roll inflation the condition ǫ ≃ M 2 P /2(∂ φ V /V ) 2 ≪ 1 is satisfied, with M 2 P ≡ (8πG) −1 the reduced Planck mass. Since we will work in a full quasi-de Sitter expansion, another useful parameter to characterize slow roll inflation is the second slow roll parameter, i.e. δ ≡ ǫ − ǫ ′ /2Hǫ ≪ 1.
In the case of the QMCU, the starting point is also a flat FLRW geometry that leads to the Friedmann and conservation equations. The field ϕ is separated into an homogeneous part ϕ 0 (η) plus small inhomogeneities δϕ(x, η). The homogeneous part satisfies where W is the potential associated to the field ϕ.
In the QMCU, it is assumed that the equation of state associated to the scalar field almost mimics that of ordinary matter, i.e. P = ωρ such that |ω| ≪ 1; the latter implies ϕ ′2 0 ≃ 2a 2 W . Consequently, the scale factor (in conformal time) evolves as a(η) ≃ η 2 /9.
The quasi-matter contraction is characterized with a small parameter |ǫ| ≪ 1, which plays the same role as the slow roll parameter in inflation. The parameter ǫ is defined as (see e.g. [22]) The case ǫ = 0 corresponds to an exact matter-dominated contracting phase (note that ǫ = ω). Furthermore, for sake of completeness we introduce another parameter such that |δ 2 | ≪ |ǫ|. The parameter |δ 2 | is analog to the δ parameter of slow roll inflation and it is related to the running of the spectral index in the QMCU model (see e.g. [22]). As is well known, it is not straightforward to accomplish a non-singular bounce within the framework of General Relativity by considering a single canonical scalar field, since the null energy condition (NEC) is violated (see for instance [14,23]). As a consequence, one possible option is to work with cosmologies within the context of modified gravity theories. In the case of the QMCU presented in [22,23], the authors worked within the framework of holonomy corrected loop quantum cosmology and teleparallel F (T ) gravity.
It is also important to note that even if a non-singular bounce cannot be achieved within General Relativity, the origin of the primordial perturbations is assumed to take place during the contracting (pre-bounce) phase of the Universe, where the curvature and energy scales are low enough to be described by General Relativity. On the other hand, one must present the conditions that need to be fulfilled such that the shape of the primordial spectrum, associated to the perturbations, remains practically unchanged when passing through the bounce. We will discuss this subject in more detail in the next section.

B. Perturbations
In the inflationary Universe and in the QMCU, one can separate the scalar field into an homogeneous part plus small inhomogeneous perturbations. Moreover, the metric associated to the spacetime, in both cases, is described by a FLRW background metric plus perturbations; which are classified as scalar, vector and tensor types (in this paper we will not consider vector perturbations). One useful quantity to describe the scalar (an also the tensor) perturbations is the so called Mukhanov-Sasaki (MS) variable. During inflation, the MS variable is defined by with Φ the gauge invariant quantity known as the Bardeen potential [65], which, in the longitudinal gauge, corresponds to the curvature perturbation. A similar expression to Eq. (8) can be used in the QMCU by replacing the fields φ ′ 0 and δφ with ϕ ′ 0 and δϕ, respectively. The advantage of relying on the MS variable is that, when expanding the action of a scalar field minimally coupled to gravity into second order scalar perturbations, one obtains δ (2) with v k the Fourier modes associated to the MS variable, z = aφ ′ 0 /H during inflation, and z = aϕ ′ 0 /H when considering the QMCU. However, it is important to note that during the bouncing phase, the action given by the Lagrangian in Eq. (9) remains the same but the expression for z changes (see Ref. [25] for an explicit calculation within F (T ) theories). On the other hand, during the contraction phase, the quantity z ′ /z can be written explicitly in terms of the QMCU parameterǭ, as in a similar fashion using the slow roll inflation parameters, i.e.
Note that, since |δ 2 | ≪ |ǫ|, the δ 2 parameter does not enter into the expression z ′ /z at first order for the QMCU.
The CSL model is based on a nonlinear modification to the Schrödinger equation; consequently, it will be advantageous to perform the quantization of the perturbations in the Schrödinger picture, where the relevant physical objects are the Hamiltonian and the wave functional. The Hamiltonian associated to L in Eq. (9) where the indexes R, I denote the real and imaginary parts of v k and p k . The canonical conjugated momentum Since In the Schrödinger picture, the wave functional Ψ[v(x, η)] characterizes the state of the system. Furthermore, in Fourier space, the wave functional can be factorized into modes components Ψ . From now on, we will deal with each mode separately. Henceforth, each mode of the wave functional, associated to the real and imaginary parts of the canonical variables, satisfies the Schrödinger equationĤ R,I k Ψ R,I k = i∂Ψ R,I k /∂η, with the Hamiltonian provided by (11). Note that one can also choose to work with the wave functional in the momentum . The standard assumption is that, at an early conformal time τ → −∞, the modes are in their adiabatic ground state, which is a Gaussian centered at zero with certain spread. This applies to both, the inflationary Universe and the QMCU. In addition, this ground state is commonly referred to as the Bunch-Davies vacuum. Thus, the conformal time η is in the range [τ, 0 − ).
Given that the initial quantum state is Gaussian, its shape will be preserved during its evolution. The explicit expression of the Gaussian state, in the field representation, is: and, equivalently, in the momentum representation Therefore, the wave functional evolves according to Schrödinger equation, with initial conditions given by = 0 corresponding to the Bunch-Davies vacuum, which is perfectly homogeneous and isotropic in the sense of a vacuum state in quantum field theory. The fact that we are introducing the wave functional in the field and momentum representations is related to the choice of the collapse operator in the CSL model, i.e., since there is no physical reason to choose one over the other, both choices are equally acceptable (at least from the phenomenological point of view). In the next section, we will show how to extract the physical quantities from the theory to be compared with the observations.

IV. THEORETICAL PREDICTIONS AND OBSERVATIONAL QUANTITIES
We begin this section by making some key remarks about the conceptual aspects of our approach and, then, we proceed to identify the relevant physical quantities that will be related with the observed data. We encourage the reader to consult Refs. [34,35,52] for a complete discussion regarding our full picture of the role played by the dynamical reduction of wave function in the cosmological setting. As a matter of fact, the relation between the observables and the predictions from the theory, using the Mukhanov-Sasaki variable during inflation and the CSL model, has been previously exposed in [50]; however, in this section we reproduce the key arguments of such a reference to make the present paper as self-contained as possible. Thus, there is no original work in the following of this section.
The main role for invoking the collapse of the wave function is to find a physical mechanism for breaking the initial homogeneity and isotropy associated to both, the quantum state and the spacetime. More specifically, we assume that a nonlinear modification to the Schrödinger equation, which in the present work is provided by the CSL model, can break the homogeneity and isotropy associated to the vacuum state and, in turn, it can generate the metric perturbations, which correspond to the primordial curvature perturbation.
Note that in the literature one can found statements suggesting that the vacuum fluctuations somehow become classical when the proper wavelength associated to the perturbations becomes larger than the Hubble radius [29,66]. Nevertheless, there is nothing in the dynamics governed by the traditional Schrödinger equation that can change the symmetry of the vacuum state, the symmetry being the homogeneity and isotropy. As a consequence, if the quantum state is perfectly symmetric and the Quantum Theory teaches us that the symmetries of a physical system must be encoded in the quantum state, then there is no clear way to describe the inhomogeneities and anisotropies of the spacetime in the quantum sense. If the quantum state of the system is perfectly symmetric, then its classical description must also be exactly symmetric. Thus, there is a lack of a proper explanation concerning the emergence of the primordial inhomogeneities and anisotropies in the Universe. That is why some non-standard interpretations of Quantum Mechanics, that make use of the Schrödinger equation (e.g. many-worlds, consistent histories, etc.), cannot provide a satisfactory answer to the problem at hand. It is important to note that the previous discussion applies to both cosmological models, the QMCU and inflation.
The modified Schrödinger equation given by the CSL model can successfully change the symmetries of the vacuum state and, at the same time, be responsible for the birth of the primordial curvature perturbation.
Specifically, in the comoving gauge, the curvature perturbation R(x, η) and the MS variable v(x, η) are related by Thus, a quantization of v(x, η) implies a quantization of R(x, η). The question that arises now is: how to relate the quantum objectsv(x, η) andR(x, η)? Furthermore, one may wonder how to relate the physical observables, such as the temperature anisotropies of the CMB, with the quantum objects that emerge from the quantum theory? The traditional answer relies on the quantum correlation functions, in particular, the two-point quantum correlation function 0|R(x, η)R(x ′ , η)|0 and its relation with the two-point angular correlation function of the temperature anisotropies δT /T 0 (n 1 )δT /T 0 (n 2 ), where the bar denotes an average over different directions in the celestial sky andn 1 andn 2 are two unitary vectors denoting some particular directions. We do not find the previous answer to be completely satisfactory, and for a detailed explanation we invite the reader to consult Refs. [34,35].
In order to illustrate our approach, we begin by focusing on the temperature anisotropies of the CMB observed today and its relation to the comoving classical curvature perturbation encoded in the quantity R. Such a relation is approximately given by (i.e. for large angular scales) On the other hand, the observational data are described in terms of the coefficients a lm of the multipolar series here θ and ϕ are the coordinates on the celestial two-sphere, with Y lm (θ, ϕ) the spherical harmonics. Given Eq. (15), the coefficients a lm can be further re-expressed in terms of the Fourier modes associated to R, i.e.
where R D is the comoving radius of the last scattering surface and j l (kR D ) the spherical Bessel function of order l of the first kind. Finally, we can include the effects of late time physics that give rise to so called acoustic peaks. These effects are encoded in the transfer functions ∆ l (k), and thus the coefficients a lm are given by where R k is the primordial comoving curvature perturbation. Also note that for large angular scales ∆ l (k) → j l (kR D ). The next step is to relate R k with the quantum operatorR k . Clearly, if one computes the vacuum expectation value 0|R k |0 and makes it exactly equal to R k , then one obtains precisely zero; while it is clear that for any given l, m, the measured value of the quantity a lm is not zero. As matter of fact, the standard argument is that it is not the quantity a lm that is zero but the average a lm . However, the notion of average is subtle, since in the CMB one has an average over different directions in the sky, while the average that one normally associates to the quantum expectation value of an operator is related to an average over possible outcomes of repeatedly measurements of an observable associated to an operator in the Hilbert space of the system (it is evident that concepts such as measurements, observers, etc. are not well defined in the early Universe).
On the other hand, we will assume that the quantity R k , i.e. the classical value associated to the Fourier mode of the comoving curvature perturbation R(x, η), is an adequate description if the quantum state associated to each mode is sharply peaked around some particular value. In consequence, the classical value corresponds to the expectation value ofR in that particular "peaked" state [53]. In other words, our assumption is that the CSL mechanism will lead to a final state such that the relation holds. Therefore, in our approach, the coefficients a lm in Eq. (18), will be given by where |Ψ corresponds to the evolved state according to the non-unitary modification of the Schrödinger equation provided by the CSL mechanism (see Refs. [46,48,49] for other ways to relate R k andR k using the CSL model, and [50] for a discussion on those approaches). Note also that |Ψ does not share the same symmetries as the vacuum state, i.e. the inhomogeneity and isotropy of the system is encoded in the quantum state |Ψ . Furthermore, Eq. (20) shows how the expectation value of the quantum fieldR k in the state |Ψ acts as a source for the coefficients a lm . A well known observational quantity is the angular power spectrum defined by We will assume that we can identify the observed value |a lm | 2 with the most likely value of |a lm | 2 ML obtained from the theory and, in turn, assume that the most likely value coincides approximately with the average |a lm | 2 . This average is over possible realizations or outcomes of the state |Ψ that results from the CSL evolution. Thus, the observed C obs. l approximately coincides with the theoretical prediction C l given in terms of the average |a lm | 2 , i.e.
Using Eq. (20), the theoretical prediction for the angular power spectrum is Moreover, if the CSL evolution is such that there is no correlation between modes (which can be justified by the fact that we are working at linear order in cosmological perturbation theory), then whereR R,I k denotes the real and imaginary part of the fieldR k (also, we assume that there is no correlation between R R k andR I k ). Therefore, Performing the integral over the angular part of k and summing over m, we obtain On the other hand, the standard relation between the primordial power spectrum and the C l is given by where P s (k) is the dimensionless scalar power spectrum defined as Henceforth, Eqs. (26) and (27) imply that the power spectrum in our approach is given by Note that the definition of the power spectrum, Eq. (28), is the canonical definition when dealing with classical random fields, where the average is over possible realizations of the random fields. In cosmology, the usual identification of the two-point quantum correlation function 0|R kRk ′ |0 with R k R k is subtle and concepts such as ergodicity, decoherence and squeezing of the vacuum state are normally invoked.
Thus, in terms of the MS variable, the scalar power spectrum in our approach is: Equation (30) is the key result from this section. It shows explicitly how to relate the quantities obtained from the quantum theory with the observed temperature anisotropies of the CMB. It also exhibits the difference between our approach and the traditional one.

V. THE CSL MODEL IN QUASI-MATTER CONTRACTION AND INFLATION
In this section, we will focus on the specific details of implementing the CSL model to the QMCU and the inflationary Universe, and the main goal will be to obtain a prediction for the power spectra.
We begin by noting that, in Eq. (30), the predictions related to the observational data are the objects Ψ|v R,I k |Ψ . Therefore, we will apply the CSL model to each mode of the field and to its real and imaginary parts. As a consequence, we will assume that the evolution of the state vector characterizing each mode of the field, written in conformal time, is given by with H R,I k given in (11). Note that the HamiltonianĤ R,I k depends on the fieldv R,I k which is defined in terms of the inflaton perturbations, but also it can be defined analogously using the perturbations of the scalar field associated to the QMCU [one has to take into account the change in z(η)]. Furthermore, we will consider that Eq. (4) can be extrapolated to a generic quantum modeF k , that is, the real and imaginary parts of the modeF k satisfy: At this point, we have to make a choice regarding the collapse operatorΘ R,I k . At first sight, the natural candidate is the MS variable, namelyΘ R,I k =v R,I k . Nevertheless, we think that in absence of a full relativistic CSL model, there is no a priori choice and, thus, the canonical conjugated momentump R,I k can also be considered as the collapse operator. In fact, in Ref. [50], we have shown that in the framework of the inflationary Universe, the momentum operator can be used as the collapse operator given that, in the longitudinal gauge, the momentum operator is directly related with the curvature perturbation.
Thus, we are going to consider four different cases: (i) The collapse operator isv R,I k during slow roll inflation.
(ii) The collapse operator isv R,I k during the QMCU.
(iii) The collapse operator isp R,I k during slow roll inflation.
(iv) The collapse operator isp R,I k during the QMCU.
We stress that only the third case, that is, the implementation of the CSL model within the inflationary framework using the fieldp R,I k as the collapse operator, was first developed in Ref. [50]. Nevertheless, we are including it in the present work for the sake of completeness. Note however that the analysis in Ref. [50] was done in the longitudinal gauge. In the present paper, we will work in the comoving gauge in all the four cases. The analysis of the three remaining cases, and in particular the implementation of the CSL model during a contracting phase of the early Universe, are presented here for the first time.
For each of these four cases, we will obtain the scalar (and tensor) power spectrum. Furthermore, the calculation of the object Ψ|v R k |Ψ 2 is identical to Ψ|v I k |Ψ 2 . Consequently, we will omit from now on the indexes R, I unless it creates confusion.
Using the Gaussian wave functions in the field representation, Eq. (13), and the probability associated to W(η) in Eq. (2), it can be shown that [47], .
The quantity (4Re[A(η)]) −1 is the standard deviation of the squared field variablev k . It is also the width of every packet in Fourier's space. In a similar manner, using the Gaussian wave function in the momentum representation, Eq. (14), along with Eq. (2), it follows that For cases (i) and (ii), it is convenient to work with Eq. (33); and for cases (iii) and (iv) with Eq. (34). Thus, to calculate v k 2 , we only need to find the two terms on the right hand side of (33) or (34), respectively. The second term on the right hand side of both equations can be found from the CSL evolution equation, Eq. (31), while the first one by using Eq. (32) with the wave function in the corresponding representation. Also, in Eqs. (33) and (34), we consider the regime −kη → 0, which correspond to the range of observational interest, that is, the regime for which the modes are larger than the Hubble radius.
Once we have computed Eqs. (33) and (34), in the corresponding case, we can substitute it into Eq. (30) to give a specific prediction for the scalar power spectrum. The actual calculations are long, so we have included them in Appendix A for the interested reader. In the following, we will show only the main results.
In case (i), our predicted scalar power spectrum during inflation is (at the lowest order in the slow roll parameter): with ν s ≡ 3/2 + 2ǫ − δ and For case (ii), we have with µ s ≡ 3/2 − 6ǭ and In both cases, (i) and (ii), we have also defined: The calculations for obtaining the tensor power spectra are very similar to the one used to obtain the scalar ones (see Appendix A for further details). In case (i), the formula obtained for the tensor power spectrum is where ν t ≡ 3/2 + ǫ. Therefore, the tensor-to-scalar ratio r ≡ P t (k)/P s (k), at the lowest order in the slow roll parameter, is given by which is exactly the same prediction as in the standard inflationary slow roll scenario. Meanwhile, in case (ii), the tensor power spectrum is where µ t ≡ 3/2 − 6ǫ = µ s . Also, for very low energy densities and curvatures, z T = a (see Ref. [22]). The tensor-to-scalar ratio is given by which is also the same as the one presented in Refs. [22,23]. Note that we have evaluated the upper limit of the integrals at η = ∞. The motivation is essentially the same as the one given in Refs. [22,23]. That is, one evaluates the scalar and power spectra at very late times corresponding to when the mode "re-enters the horizon", or more precisely when k = |aH| during the expanding (post-bounce) phase. The previously presented cases (i) and (ii) correspond to selectingv R,I k as the collapse operator. Next, we focus on the results for cases (iii) and (iv), which correspond to choosep R,I k as the collapse operator. For case (iii), we obtain: where we have defined ν s ≡ 1/2 + 2ǫ − δ and In case (iv), the corresponding expression results with the definitions µ s ≡ 3/2 − 6ǫ and The constants c 1 , c 2 and c 3 are shown in Appendix A. In both cases, (iii) and (iv), we have the following definitions The predictions for the tensor-to-scalar ratios are exactly the same as the ones presented in cases (i) and (ii) (see Appendix A).
We end this section by summarizing the main results. We have applied the CSL model to the inflationary Universe and to the QMCU. Moreover, in order to employ the CSL model, we need to choose the collapse operator. We have chosen to work withv R,I k andp R,I k as the collapse operators. Henceforth, we have obtained the scalar power spectra in four different cases Eqs. (35), (37), (44) and (46). On the other hand, introducing the CSL mechanism does not affect the tensor-to-scalar ratio r. Specifically, if one works within the standard inflationary scenario, then the prediction for r is equal to the standard one given by slow roll inflation; meanwhile, if one adopts the QMCU framework, then the predictions are equal to the ones presented in Refs. [22,23].

VI. DISCUSSION ON THE CSL INSPIRED POWER SPECTRA
In this section, we will discuss the implications of the results obtained in the previous section. In particular, we will compare our predicted scalar power spectra with the standard one.
The scalar power spectrum predicted by slow roll inflation is traditionally expressed as [36,67] where k 0 is a pivot scale, and the amplitude A s and the spectral index n s are given by On the other hand, we have four different expressions for the scalar power spectrum, corresponding to the four cases mentioned at the beginning of Sect. V. In the following, we will analyze each one of them, but first we will make a few observations regarding the parameter λ k .
The dependence on k in the parameter λ k encodes the "amplification mechanism", which is characteristic of dynamical reduction models (see Refs. [48,49]). One possible way to determine the exact dependence on k, and perhaps the simplest, is by dimensional analysis. That is, the main evolution equations are given in Eqs. (31) and (32); consequently, in order for those equations to be dimensionally consistent, the fundamental dimensions of λ k change depending on the fundamental units associated to the collapse operatorΘ R,I k . Moreover, we expect that λ k is directly related to λ, i.e. the CSL parameter, which clearly must be the same in all physical situations (cosmological or otherwise). Moreover, taking into account that we are working in units in which = c = 1, the fundamental dimension of λ is [Length] −1 .
Thus, in the case where the collapse operator is chosen to bev R,I k , the most natural expression of λ k , which is consistent with the dimensions of all terms involving the dynamical equations, is And in the case where the selected collapse operator isp R,I k , such an expression is where λ is the CSL parameter, with the same numerical value in all cases. From now on, we will assume that λ k takes the form of Eqs. (51) and (52) depending on the chosen operator acting as the collapse operator.
A. The CSL power spectra during inflation Let us begin the discussion by working within the framework of the inflationary Universe, analyzing cases (i) and (iii).
The scalar power spectrum given in Eq. (35), corresponding to case (i), can be written in a similar form to the one showed in Eq. (49). As usual, the power spectrum can be evaluated at the conformal time where the pivot scale "crosses the horizon"; or more precisely, when −k 0 η = 1 (i.e. k 0 = aH) during the inflationary epoch. Furthermore, the different coefficients that multiply each term of the function F 1 (λ k , ν s ) involve the quantity ν s . For these terms, we can approximate ν s ≃ 3/2 without loss of generality. However, note that such approximation cannot be done to the powers of k involving ν s because these are directly related to the scalar spectral index n s , for which the value n s = 1 is ruled out. Furthermore, in order to provide a suitable normalization for the CSL power spectra, we multiply and divide by the quantity λ|τ |. Thus, the power spectrum in Eq. (35) can be rewritten as: where and C 1 (k) ≡ F 1 (λ k = λk, ν s ≃ 3/2)/λ|τ |; that is, [expressions for ζ k and θ k are given in Eq. (39) with λ k = λk].
Within the inflationary framework, and with the same arguments followed to arrive to Eq. (53), we can write the power spectrum Eq. (44), corresponding to case (iii), in the following form: where A s and n s are the same as in Eq. (54), and C 3 (k) ≡ F 3 (λ k = λ/k, ν s ≃ 1/2)/λ|τ |. Thus, [expressions forζ k andθ k are given in Eq. (48) with λ k = λ/k] Let us make some remarks. Notice that the scalar index predicted by the CSL power spectra is exactly the same as the standard one from slow roll inflation, but the amplitude is slightly different. The difference between the standard amplitude and the one using the CSL model is a factor of λ|τ |/2 [see Eqs. (50) and (54)]. The reason for the factor 1/2 can be traced back to Eq. (30), since in our approach the power spectrum receives an equal contribution from the expectation values v R k 2 and v I k 2 . However, the factor 1/2 will not have any important observational consequences. On the other hand, the factor λ|τ |, which comes from the normalization of C 1 (k) and C 3 (k), does modify the standard predicted amplitude. A quantitative analysis will be done in the next section.
A second remark has to do with the following. It is well known that there is a minimum number of e-foldings for inflation related to the solution of the "horizon problem", and this minimum number depends on the characteristic energy of inflation. A shared characteristic of the functions C 1 (k) and C 3 (k) is that they include the quantity τ , which represents the conformal time at the beginning of inflation. This quantity depends on the energy scale at which inflation ends, which is associated to the inflaton potential V at that time, and the number of e-foldings corresponding to the total duration of inflation.
Third, note that another important feature of the CSL power spectra in inflation is that the function C 1 (k), corresponding to the case in which the collapse operator isv R,I k , depends explicitly on the conformal time η, whilst the function C 3 (k), which corresponds to the case whenp R,I k is the collapse operator, does not exhibit such time dependence. The time dependence on the power spectrum when the collapse operator isv R,I k has been previously noted by other authors [46,48]. Nevertheless, the exact form of their predicted power spectrum is different from the one shown here. As a matter of fact, that difference is illustrated by considering the limiting case λ k = 0. In such mentioned works, for λ = 0 (i.e. standard Schrödigner evolution), their predicted power spectrum is the same as the traditional one. Contrarily, in our approach, if λ k = 0 then P s (k) = 0, which is consistent with our point of view regarding the role played by the CSL model. In any case, even if the pictures used for the role of the CSL model are different between our work and the one in Refs. [46,48], the time dependence on the power spectrum is shared.
In order to continue, we choose to evaluate the power spectrum (or equivalently the function C 1 (k)) at the conformal time when inflation ends, which we denote by η f . We think it is consistent with the previous calculations in which the power spectrum was obtained in the limit −kη → 0, which is satisfied by the value η f . The precise value of η f depends mainly on the characteristic energy scale of inflation and the number of e-foldings assumed for the full inflationary phase N ≡ ln[a(η f )/a(τ )].
Readers familiar with previous works, can check that our expression for the scalar power spectrum Eq. (56), which features the function C 3 (k), is essentially the same as the one obtained in Ref. [50]. The difference is that in the present paper we chose to work in the comoving gauge (where R represents the curvature perturbation), whilst in the aforementioned reference we worked in the longitudinal gauge, where the Bardeen potential Φ corresponds to the curvature perturbation. Therefore, we find reassuring that even having worked in different gauges, the expression for the power spectrum, when the collapse operator is the momentum associated to the Fourier's mode of the MS variable, is the same and it has the attractive feature that does not depends on the conformal time. Figures 1 and 2 show different plots for the functions C 1 (k) and C 3 (k), respectively. In both cases, we have considered the value of the CSL parameter as λ GRW = 1.029 × 10 −2 Mpc −1 , which corresponds to a value favored by experimental data [55,56,59,60]. The various plots in each figure correspond to different values of the characteristic energy of inflation V 1/4 , and the total e-foldings N that inflation is assumed to last, which also set the values of τ and η f . The values of k considered correspond to these of observational interest, i.e. we consider k in the range from 10 −6 to 10 −1 Mpc −1 .
As we can observe, the functions C 1 (k) and C 3 (k) exhibit an oscillatory behavior around the unit. For increasing values of k, the oscillations decrease in amplitude. However, we note that even for decreasing values of k the functions C 1 (k) and C 3 (k) are very close to 1. Consequently, for the chosen values of λ, V 1/4 and N , the functions C 1 (k) ≃ C 3 (k) ≃ 1. That means that the shape of the angular power spectrum C l will not be very different from the standard one, but the amplitude could vary (a complete analysis will be presented in the next section).
Additionally, the fact that C 1 (k) depends on the conformal time does not seem to affect its behavior in a significant manner. In fact, it is closely similar to the one of C 3 (k), which does not depend on the conformal time. That means that the contribution from the time dependent term (i.e. the last term in Eq. (55)), to the total value of the function C 1 (k) is negligible when −kη → 0.

B. The CSL power spectra in the QMCU
We switch now the discussion to the framework of the QMCU, i.e. cases (ii) and (iv), which correspond to selectinĝ v R,I k orp R,I k as the collapse operator, respectively. The scalar power spectra given in both cases, i.e. Eqs. (37) and (46), can also be written in a manner similar to the standard spectrum Eq. (49). Once again, following Refs. [22,23], we choose to evaluate the spectrum at the conformal time where the pivot scale "reenters the horizon" k 0 = |aH|, which happens at late times during the expansion phase of the Universe (consequently the upper limit of the integral is evaluated at η → ∞). We approximate (for the same arguments as in the previous subsection) µ s ≃ 3/2 in the coefficients of the terms in expressions F 2 (λ k , µ s ) and F 4 (λ k , µ s ) (but not in the powers of k as these powers are directly related to the spectral index n s ). Additionally, the parameter λ k is assumed to be λ k = λk for case (ii) and λ k = λ/k in case (iv). Moreover, we multiply and divide by a factor of λ|τ | in order to properly normalize the expressions F 2 and F 4 .
Henceforth, the scalar power spectrum for case (ii), Eq. (37), will be written as where and C 2 (k) ≡ F 2 (λ k = λk, µ s ≃ 3/2)/λ|τ |; thus, On the other hand, the power spectrum for case (iv), Eq. (46), will be written as with A s and n s the same as in Eq. (59), and C 4 (k) ≡ F 4 (λ k = λ/k, µ s ≃ 3/2)/λ|τ |. Hence, As in the case of the inflationary Universe, the predicted value of the scalar spectral index n s is not affected by the CSL model. In fact, it has the same expression as that of the QMCU original models presented in Refs. [22,23].
Nevertheless, as can be seen in Eq. (59), the amplitude of the spectrum is modified by an extra factor of λ|τ |/2 with respect to the original QMCU model, that is, In this case, τ corresponds to the beginning of the quasi-matter dominated period. Regarding the amplitude of the spectrum in the QMCU model, when the background evolution is driven by a matter dominated Universe, it can be obtained analytically working within F (T ) gravity or LQC. In the teleparallel F (T ) case, the original amplitude is while in the LQC case, the original amplitude is where ρ P is the Planck energy density, C ≃ 0.9159 is Catalan's constant and ρ c is called the critical density, which corresponds to the energy density at which the Universe bounces, both expressions for the amplitude can be consulted in Refs. [22,23]. Thus, in order to obtain an amplitude in both cases (the teleparallel gravity case and the LQC case) that is consistent with that obtained from the CMB data (i.e. A s ≃ 10 −9 ), and taking into account the extra factor of λ|τ |/2 coming from the CSL model, the value of the energy density at the bouncing point must satisfy Generically λ|τ | ≫ 1; hence, the CSL model introduces an extra constriction to the QMCU, that is, ρ c ≪ ρ P . In the next section we will perform a more quantitative analysis.
The functions C 2 (k) and C 4 (k) share a characteristic feature, namely they depend explicitly on −kη, which comes from a series expansion around −kη → 0. Consequently, we choose to evaluate η = η f corresponding to the end of the quasi-matter domination stage or the onset of the bouncing phase. One can also define the total number of e-foldings N ≡ ln[a(τ )/a(η f )] for the duration of the quasi-matter dominated phase; however, notice that in this case, since there is no horizon problem, there is no minimum value of N . Another important aspect is that if λ = 0 then C 2 (k) = C 4 (k) = 0. In other words, if the evolution of the state vector is completely unitary, then there are no perturbations of the spacetime at all and the state vector continues being perfectly symmetric, which is consistent with our conceptual framework.
In Figs. 3 and 4, we show different plots for the functions C 2 (k) and C 4 (k), respectively. In both cases, we have considered the value λ GRW = 1.029 × 10 −2 Mpc −1 . The various plots in each figure correspond to different values of τ and η f . The values of k considered correspond to the values of observational interest, hence, we consider k in the range from 10 −6 to 10 −1 Mpc −1 . As we can see, the functions C 2 (k) and C 4 (k) exhibit the same oscillatory behavior around the unity as its counterparts during inflation. Also, the amplitude of each oscillation decreases for increasing values of k.
We end this section with a few comments regarding the dependence of the power spectrum on η in cases (i), (ii) and (iv). In case (i), which corresponds to selecting the MS variable as the collapse operator during inflation, the term containing the η dependence is the last one of C 1 (k) in Eq. (55). On the other hand, since the amplitude associated to the modes R k is "frozen" on super-Hubble scales, the behavior of C 1 (k) will not change for super-horizon modes. As a matter of fact, the plots in Fig. 1 show that C 1 (k) is essentially a constant in the limit −kη → 0, which means that the term containing the η dependence is sub-dominant in such a limit. In cases (ii) and (iv), corresponding to the framework of the QMCU, the behavior of the functions C 2 (k) and C 4 (k) are very similar to that of C 1 (k) (see Figs. 3 and 4). That is, they are practically a constant in the limit −kη → 0, which means that the terms involving η, i.e. the last terms of Eqs. (60) and (62), are sub-dominant in the super-Hubble limit. Nevertheless, since in cases (ii) and (iv) the Universe approaches a non-singular bounce, it might be the case that, when the mode "reenters the horizon" (k ≫ |aH|) during the bouncing phase, a modification of the dynamical evolution of the functions C 2 (k) and C 4 (k) would occur. However, if the duration of the bouncing phase is short enough then one could intuitively consider that the spectrum is left unchanged (although counterexamples exist in the literature [68]). Therefore, one could perform a full analysis regarding the CSL model during the bounce within the QMCU. Nonetheless, we will take a pragmatical approach and assume that the shape of the spectrum, provided by the functions C 2 (k) and C 4 (k), survives the bouncing phase and, then, we will use the observational data to further constraint or completely discard the predicted spectra. In case that the predicted spectra are consistent with observations, one can proceed to perform the full-fledged analysis of implementing the CSL model to the QMCU during the bouncing phase and study the possible corrections that may arise from passing the perturbations through the bounce. This subject, however, will not be explored in the present paper.
In the next section, we will explore the implications of the predicted spectra using the observational data.

VII. EFFECTS ON THE CMB TEMPERATURE SPECTRUM AND ITS IMPLICATIONS ON THE COSMOLOGICAL PARAMETERS
The aim of this section is to analyze the viability of the CSL model by comparing the corresponding predictions with the ones coming from the best fit canonical model to the CMB data. In particular, we will focus on the power spectra obtained using the CSL model and its effect on the angular power spectrum.
In order to perform our analysis, we start by setting the cosmological parameters of our fiducial model, which will be used as a reference to compare with the CSL inspired spectra. The fiducial cosmology will be the best fitting flat  Table 4 presented by the latest Planck Collaboration [5].
Furthermore, we recall that the primordial power spectrum and the angular power spectrum are related by Eq.
Hence, we will use the CSL predicted power spectra P s (k) = A s (k/k 0 ) ns−1 C i (k) with i = 1, 2, 3, 4, which correspond to the four different cases that we have considered so far. We will focus first on the inflationary model of the early Universe and, then, on the QMCU model. Also, note that the fiducial model corresponds to: P s (k) = A s (k/k 0 ) ns−1 . The precise prediction for the angular power spectrum will be obtained by using the Boltzmann code CAMB [69], with the aforementioned cosmological parameters.
A. The angular power spectrum and the CSL model during inflation During inflation, the CSL power spectra is characterized by the functions C 1 (k) and C 3 (k), Eqs. (55) and (57), with standard spectral index n s − 1 = −4ǫ + 2δ and amplitude The output of the CAMB code, that is, the temperature autocorrelation power spectrum of the fiducial model and the one provided by the CSL model during inflation are indistinguishable; thus, we have decided not to show the plots. Instead, we present the the relative difference, which we define as Figure 5 shows the relative difference between both predictions. On the left, we have chosenv R,I k as the collapse operator while on the right we have chosenp R,I k . In both cases, we have set an energy scale of 10 −5 M P for the energy at which inflation ends, and a total amount of inflation corresponding to 65 e-foldings.
We observe that the relative difference between the fiducial spectrum and the one predicted using the CSL model with, for instance λ GRW , is practically null (the highest difference is around 0.01%). This statement applies to both elections of the collapse operator and for other λ values listed in Table I (not shown in the figure).
We have also checked that the essentially null relative difference between the fiducial model and the CSL model during inflation is also present in the E polarization autocorrelation power spectrum C EE l and the temperature polarization cross correlation power spectrum C TE l .  Table I and text for details. Table I: Estimation of ǫ from Eq. (71). We have used |τ | = 7.8 × 10 7 Mpc, which corresponds to V 1/4 = 10 14 GeV and N = 65, and the four values of λ shown below. Also, we have estimated the tensor-to-scalar ratio using r = 16ǫ. Note that only the λGRW case is compatible with the latest observations from Planck Collaboration [5]. On the other hand, the amplitude of the power spectrum A s consistent with the CMB data is A s ≃ 10 −9 [5]. Henceforth, the amplitude obtained using the CSL model, as shown in Eq. (68), must satisfy Clearly, different values of λ will have an effect on the amplitude of the spectrum. Assuming that the pivot scale k 0 crosses the Hubble radius at an energy scale of V 1/4 0 = 10 −4 M P (i.e. one order of magnitude less than the presumed energy at which inflation ends), an estimate for ǫ can be calculated. Therefore, the above equation leads to Table I shows the different values of ǫ obtained by considering several λ values. Also, in the same table, we provide an estimate for the tensor-to-scalar ratio r (recall that the CSL model predicts the same relation as standard inflation, i.e r = 16ǫ). From Table I, it can be seen that only the value corresponding to λ GRW is consistent with both, the observed shape and amplitude of the spectrum. In particular, assuming a characteristic energy scale of inflation of 10 −4 M P ≃ 10 14 GeV, a total amount of inflation corresponding to N = 65, and the value of λ GRW , we obtain an angular spectrum with a shape and an amplitude that is indistinguishable from the fiducial model, which we know is consistent with the observational data. The amplitude of the spectrum for this particular set of values leads to an estimate for the slow roll parameter and the tensor-to-scalar ratio of ǫ ≃ 10 −4 and r ≃ 10 −2 , respectively. Those values of ǫ and r are consistent with the ones presented by the latest results of the Planck Collaboration [6].
It is also instructive to mention that, if future observations confirm the results of the BICEP2 Collaboration [70], i.e. r ≃ 0.2, then the value of λ GRW would not be compatible with the values of V and N used in Table I. In fact, the  Table II also achieve an excellent fit (not shown), but the analysis done in this work does not allow prefer one value over another.
λ value that might be compatible would be one such that λ ≪ λ GRW . That would open a new range of parameter space to explore in addition to considering other experimental setups different from cosmological ones [55,56].
B. The angular power spectrum and the CSL model during the quasi-matter contracting phase In the QMCU framework, the power spectra are characterized by the functions C 2 (k) and C 4 (k), i.e. Eqs. (60) and (62), respectively. The predicted scalar spectral index is n s − 1 = 12ǫ, and the amplitude is given by Notice we have approximated the integral that appears in the amplitude, corresponding to Eq. (59), by ρ c /ρ P [see Eqs. (64) and (65)]. The output of the CAMB code, that is, the temperature autocorrelation power spectrum of the fiducial model and the one provided by the CSL model during the QMCU are also indistinguishable. Thus, we present again the relative difference defined in Eq. (69) where now C CSL l corresponds to the angular power spectrum during the QMCU. Figure 6 shows the relative difference between both predictions. On the left, we have chosenv R,I k as the collapse operator while on the right we have chosenp R,I k . In both cases, we have assumed a total duration of 50 e-foldings for the quasi-matter contracting phase and a conformal time |τ | ≃ 10 8 Mpc corresponding to the beginning of the contracting stage. We show only the plot of S(l) corresponding to the value of λ Adler merely as an illustrative example; the plots for the values corresponding to λ 1 , λ 2 , λ GRW follow the exactly same behavior as the one shown in Fig. 6.
We found no difference between the fiducial spectrum and the one provided by the CSL model for the four values of λ listed in Table II; the highest relative difference is around 0.1%. This statement applies to both elections of the collapse operator. Finally, we have also checked that the essentially null relative difference between the fiducial model and the CSL model during the QMCU is also present in the E polarization autocorrelation power spectrum C EE l and the temperature polarization cross correlation power spectrum C TE l . On the other hand, the amplitude of the power spectrum A s consistent with the CMB data is A s ≃ 10 −9 [5]. Therefore, the ratio ρ c /ρ P , obtained using the CSL model, as shown in Eq. (72), must satisfy Consequently, an estimate for the energy scale of the critical energy E c is In Table II, we show the different values of κ 1 and κ 2 by using the chosen λ values, with |τ | = 1.15 × 10 8 Mpc. We infer that the four values of λ considered are consistent with a critical energy scale in the range (10 −3 M P , 10 −6 M P ).
It is worthwhile to mention that, in the QMCU, the spectral index n s and the tensor-to-scalar ratio r are not related each other as in the standard inflationary paradigm [see Eqs. (43) and (59)]. However, the spectral index, along with the running of the spectral index α s ≡ dn s /d ln k, are the two main parameters of the QMCU model used to compare with the observational data [22][23][24]. The CSL model applied to the QMCU does not affect those parameters. The only observable affected by the CSL model is the amplitude of the spectrum (which is not related to the parameter r as in the standard spectrum). Consequently, in order to put an upper bound to the energy scale at which the bouncing phase begins, and which would be equivalent to set a constraint on the parameter λ, one should consider a specific theoretical model of the QMCU (i.e. to choose a specific dynamics and a potential of the field ϕ). Therefore, in the QMCU, and with the same degree of accuracy of past works dealing with the same model, all of the four values of the CSL parameter λ considered here yield consistent predictions with observational data; specifically, the predictions regarding the shape of the spectrum, the scalar spectral index and the running of the spectral index.
On the other hand, note that the information that could discriminate among different values of λ is codified in the amplitude of the spectrum. The predicted amplitude of the spectrum depends on the critical energy density ρ c (the value of the energy density at the bouncing time), which is model dependent.

VIII. CONCLUSIONS
The CSL model is a physical mechanism that attempts to provide a solution to the measurement problem of Quantum Mechanics by modifying the Schrödinger equation. The CSL model can be referred as an objective reduction mechanism or "effective collapse" of the wave function, and one of the main elements of this model is the collapse operator, i.e. the operator whose eigenstates correspond to the evolved states by the collapse mechanism. Also, in principle, it is possible to apply such a mechanism to any physical system.
In this work, we have applied the CSL model to the early Universe by considering two cosmological models: the matter bounce scenario (MBS) and standard slow roll inflation. Additionally, we have considered two different collapse schemes, one in which the field variable (given in terms of the Mukhanov-Sasaki variable) serves as the collapse operator, and other scheme where the collapse operator is the conjugated momentum.
In all cases, we have found a prediction for the primordial power spectrum, which is a function of the standard parameters of each cosmological model, and also of the CSL parameter λ. Although the exact expressions for the primordial power spectra are different in each case, there are features that are essentially the same as its standardnon-collapse-counterparts. Specifically, the predictions for the scalar spectral index and the tensor-to-scalar ratio are exactly the same as the ones given in the MBS and slow roll inflation without collapse. On the other hand, in each case, the shape of the spectrum is modified by a function of the wave number k, associated to the modes of the field, and by the inclusion of the λ parameter. However, for a suitable choice of values corresponding to the parameters of the cosmological models, there is no significant change in the prediction for the CMB angular power spectrum (i.e. the C l 's) that can be distinguished from the canonical flat ΛCDM model.
Meanwhile, the prediction for the amplitude of the spectrum is modified directly by the parameter λ. We have empirically explored the range of values of λ, from the originally value suggested by Ghirardi-Rimini-Weber (GRW) λ GRW ≃ 10 −16 s −1 [41], to the one given by Adler λ Adler ≃ 10 −8 s −1 [58]. In the case of slow roll inflation, we have found that for a characteristic energy scale of 10 14 GeV and a total amount of inflation of 65 e-folds, only the value suggested by GRW is compatible with the observational bound of the amplitude; other values of λ greater than λ GRW , e.g. λ Adler , cannot be made compatible with the observed amplitude (because that would require values for the slow roll parameter such that ǫ > 1). In the MBS case, we have found that the modification in the predicted amplitude of the spectrum, given by the λ parameter, causes that the critical energy density ρ c , i.e. the energy density at which the bouncing phase begins, to be several orders of magnitude less than the Planck energy density ρ P . The precise number of orders of magnitude varies according to the value of λ. For instance, by assuming λ GRW and a total amount of ∼ 50 e-folds for the matter dominated contracting phase, we have ρ c ≃ 10 −15 ρ P . The latter relation is obtained by requiring the compatibility between the predicted amplitude of the scalar power spectrum and the one from the Planck CMB data A s ∼ 10 −9 .
In conclusion, it was possible to incorporate the CSL model into the cosmological context again, in particular when dealing with the quantum-to-classical transition of the primordial inhomogeneities. Moreover, it is remarkable that our implementation of the CSL model yields predictions that are also in agreement with experiments in the regimes so far investigated empirically. Those experiments involve values of the CSL parameter λ that have been tested in laboratory settings, quite disengaged from the cosmological framework. We acknowledge that at this stage, the application of CSL model to the early Universe, as done in this manuscript, can be seen as an ad hoc employment. However, the fact that the predictions can be empirically tested make us hopeful that future studies will overcome the perceived shortcomings. and for case (ii) the corresponding expansion is Re A k (η) ≃ λ k 2k(−m − 1) (−kη) + ζ −2m k sin(πm + 2mθ k ) sin(πm) kπ 2 −2m Γ(−m) 2 (−kη) −2m−1 . (A4) In both cases, ζ k and θ k are given in Eq. (39). Next, we focus on the first term of Eq. (33), i.e. v 2 k . It will be useful to define the following quantities: Q ≡ v 2 k , R ≡ p 2 k and S ≡ p kvk +v kpk .
The equations of evolution for Q, R and S are obtained using Eq. (32), withΘ k =v k . That is, Therefore, we have a linear system of coupled differential equations, whose general solution is a particular solution to the system plus a solution to the homogeneous equation (with λ k = 0). After a long series of calculations we find: and the constants C 1 , C 2 and C 3 are found by imposing the initial conditions corresponding to the Bunch-Davies vacuum state: Q(τ ) = 1/(2k), R(τ ) = k/2 and S(τ ) = 0. Equation (A7) is exact; expanding it again around −kη → 0 yields Q(η) ≃ π 2k 2 sin 2 (nπ) where Using the above results, we can compute the quantity v k 2 using Eq. (33). In case (i), we substitute Eqs. (A3) and (A9) into Eq. (33), obtaining v k 2 = Q(η) − 1 4ReA k (η) ≃ π 2k 2 sin 2 (nπ) (A11) With the expression in (A11) at hand (which is valid for v R k 2 and v I k 2 ), and using Eq. (30), our predicted scalar power spectrum during inflation (at the lowest order in the slow roll parameter) is given in Eq. (35). (we have also used that during inflation z 2 (η) ≃ 2ǫM 2 P /(H 2 η 2 ) and m = n = 3/2 + 2ǫ − δ). (A12) Since in this case m = −5/2 + 6ǫ and n = 3/2 − 6ǫ, and considering only the first dominant term in the expansion around −kη → 0, we finally obtain the expression for the power spectrum presented in Eq. (46). The procedure to obtain the tensor power spectra is analogous to the one outlined in the previous subsection, but clearly the difference is thatΘ R,I k =p R,I k . In the following, we will only present the results. For case (iii), the tensor power spectrum is given by and, consequently, the tensor-to-scalar ratio is r = 16ǫ, which is the same as the standard prediction of slow roll inflation.
For case (iv), the formula for the tensor power spectrum is Therefore, the tensor-to-scalar ratio is exactly the same as the one shown in Eq. (43).