Emergence of inflationary perturbations in the CSL model

The inflationary paradigm is the most successful model that explains the observed spectrum of primordial perturbations. However, the precise emergence of such inhomogeneities and the quantum-to-classical transition of the perturbations has not yet reached a consensus among the community. The Continuous Spontaneous Localization model (CSL), in the cosmological context, might be used to provide a solution to the mentioned issues by considering a dynamical reduction of the wave function. The CSL model has been applied to the inflationary universe before and different conclusions have been obtained. In this letter, we use a different approach to implement the CSL model during inflation. In particular, in addition to accounting for the quantum-to-classical transition, we use the CSL model to generate the primordial perturbations, that is, the dynamical evolution provided by the CSL model is responsible for the transition from a homogeneous and isotropic initial state to a final one lacking such symmetries. Our approach leads to results that can be clearly distinguished from preceding works. Specifically, the scalar and tensor power spectra are not time-dependent, and retains the amplification mechanism of the CSL model. Moreover, our framework depends only on one parameter (the CSL parameter) and its value is consistent with cosmological and laboratory observations.


I. INTRODUCTION
The inflationary paradigm is held among the majority of cosmologists as a successful model for addressing the primordial inhomogeneities that represent the seeds of cosmic structure. In fact, recent observations from the Cosmic Microwave Background (CMB) radiation [1][2][3][4][5][6] are quite consistent with the standard prediction from the simplest inflationary model, namely, the prediction of a nearly scale-invariant power spectrum. On the other hand, the exact physical mechanism responsible for the generation of the primordial curvature perturbations, associated to a highly-Gaussian stochastic classical field is still a matter of debate. In particular, inflation is based on a combination of quantum mechanics and general relativity, two theories that are difficult to merge at both the conceptual and technical level. Therefore, it is expected that when these theories are used in the same footing, some difficulties might arise. Specifically, the quantumto-classical transition, a subject which has been present since the conception of the quantum theory, is an issue that has not been fully resolved in the early inflationary universe.
Earlier works based on decoherence [7] and the consideration that the initial vacuum state of the universe evolves into a highly squeezed state [8], led to a partially understanding of the issue, namely, that the predictions from the quantum theory are indistinguishable from those of a theory in which the random fluctuations are the result of a classical stochastic process [9]. Nevertheless, this argument by itself cannot address the fact that a single (classical) outcome emerges from the quantum the- * Electronic address: gleon@df.uba.ar † Electronic address: gabriel@iafe.uba.ar ory. In other words, decoherence (and the squeezing of quantum states) cannot solve the quantum measurement problem [10,11], a complication that, within the cosmological context, is amplified due to the impossibility of recurring to the "for all practical purposes" argument in the familiar laboratory situation; i.e. it is not clear how to define in the primordial universe entities such as observers, detectors, etc.
A more precise formulation of the problem at hand can be stated as follows: the inflationary universe is characterized by a background space-time that is completely homogeneous and isotropic, equivalently, the quantum state of the matter field, namely the vacuum state, is also perfectly homogeneous and isotropic. On the other hand, the present universe is the result of the evolution of primordial density inhomogeneities. Hence, the problem is: how do the inhomogeneities and anisotropies originated from the initial highly symmetric (the symmetry being the homogeneity and isotropy) state of the universe, described by both the background space-time and the vacuum state of the matter fields, given that the dynamical evolution, provided by Schrödinger and Einstein equations, does not break the translational and rotational symmetries? One possible answer to this question is to invoke the collapse of the wave function but without relying on any external objects, e.g. [12,13]. Other approaches to this problem have been based on Bohmian quantum mechanics, see for instance [14,15]. The idea of invoking a self-induced collapse in order to generate the primordial perturbations has been explored in great detail in previous works, e.g. [16][17][18][19][20][21][22][23]. Moreover, in [24] some generic collapse schemes have been tested using observational data coming from the 7-year release of the WMAP collaboration [25] and the Sloan Digital Sky Survey [26]. Therefore, the attempt to solve the aforementioned issue is not just a matter of philosophical concern but a relevant problem from the theoretical point of view, yielding predictions that can be confronted with observational data.
Furthermore, the proposal of a self-induced collapse of the wave function has been an active line of research since the early ideas of L. Diósi [27] and R. Penrose [28], advocating for gravity as the main agent triggering the collapse. Also, the Ghirardi-Rimmini-Weber (GRW) [29] model was among the first attempts to introduce an objective collapse model. In past decades, the Continuous Spontaneous Localization (CSL) model, which can be viewed as a continuous version of the GRW model, has been regarded as a promising model that can provide a solution to the quantum measurement problem [30,31]. In particular, the CSL mechanism is based on a non-linear stochastic modification of the standard Schrödinger equation, in this way, spontaneous and random collapses of the wave function occur all the time, to all particles, regardless they are isolated or interacting. The idea behind the CSL model [32], sometimes referred to as the "amplification mechanism," is that the collapses must be rare for microscopic systems, in order to not alter their quantum behavior as described by the Schrödinger equation. At the same time, their effect must increase when several particles are hold together forming a macroscopic system. Moreover, the testable predictions made by the CSL model are now considered to be feasible within the current available technology [33]. On the other hand, the CSL model is, at this stage, a nonrelativistic model and therefore a complete generalization to quantum fields is still under development [34]. Nevertheless, the particular features of the CSL model, e.g. the absence to rely on external agents to solve the measurement problem, make it a viable candidate, at least from the phenomenological point of view, to address the issue concerning the emergence of the primordial perturbations.
In [35], the CSL model was applied to the inflationary universe for the first time. There, the authors followed an approach in which the collapse affected in the same way all modes of the inflaton field and, as a consequence, the amplification mechanism of the CSL model was lost. Additionally, their predicted scalar power spectrum contained some features that conflicted with the nearly scaleinvariant power spectrum, which is consistent the CMB data. In order to retain the amplification mechanism, other authors [36], proposed a phenomenological manner in which the CSL model would affect each mode of the inflaton field. Additionally, they also assumed that the modification to Schrödinger equation, provided by the CSL mechanism, was by introducing a time-dependent parameter. This last assumption resulted in a prediction for the scalar (and also tensor) power spectrum with the correct shape (for a particular combination of their three free parameters) but that was time-dependent. In order to overcome this shortcoming, the power spectrum was chosen to be evaluated at the end on inflation, however, this choice makes their CSL collapse parameter extremely small (∼ e −120 ) [37].
A shared feature of the works in Refs. [35,36] is that the authors worked in a joint metric-matter quantization of the perturbations characterized by the Mukhanov-Sasaki variable [38]. On the other hand, in Ref. [39] it was followed a different approach to the problem, this is, the authors successfully applied the CSL collapse mechanism to the inflationary universe but within the semiclassical gravity framework. Consequently, as pointed out in [22,40], the amplitude of the primordial gravitational waves is exactly zero at first-order in the perturbations. Therefore, a confirmed detection of primordial gravity waves, would rule out this approach.
Another distinction between the approaches in Refs. [35,36] and [39] is the specific role played by the selfinduced collapse. This difference is subtle but important. The authors of [35,36] employed the CSL mechanism with the intention of localizing the evolved inflaton vacuum state in an eigen-state of the Mukhanov-Sasaki variable. Meanwhile, the authors of [39] attributed to the self-induced collapse the action of generating the primordial curvature perturbation, irrespectively of which eigen-state the initial state evolves into. Therefore, if no self-induced collapse occurs, the predicted power spectrum in Refs [35,36] is exactly the same as the standard one; on the contrary, in Ref. [39] the absence of a collapse results in no perturbations of the space-time at all, consequently, the predicted scalar power spectrum is exactly zero.
In the present Letter, we adopt the same role for the collapse, provided by the CSL model, as the one presented in [39], but we apply it to the joint metric-matter quantization of the inflaton field. In this way, we obtained predictions that are different from the ones in [35,36] and are also quite consistent with the observational data. In particular, our predicted scalar and tensor power spectra are not time-dependent, and the implementation of the CSL mechanism retains the amplification mechanism. Also, our model introduces just one free parameter and in contrast with Ref. [39], it predicts a non-vanishing tensor power spectrum at first-order in perturbations.
The article is organized as follows: in Section II we review some basics about the CSL collapse model applied to the inflationary universe, and we obtain the scalar power spectrum within such a framework; in Section III we show our results for the power spectrum of tensor modes and the tensor-to-scalar ratio; in Sections IV and V we make a discussion of our results and compare them with previous works, and finally in Section VI we summarize our conclusions.

II. THE CSL MODEL AND THE SCALAR POWER SPECTRUM
In this Section, we will summarize some concepts regarding the CSL collapse model, and we will use it to obtain a prediction for the primordial power spectrum of scalar perturbations.

A. Classical description of the perturbations
The inflationary universe is described by Einstein equations G ab = 8πGT ab (c = 1) and the dynamics of the matter fields dominated by the inflaton. We will assume the simplest inflationary model, which is a single scalar field ϕ in the slow-roll approximation. The background space-time is accelerating in a quasi-de Sitter type of expansion, 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 field V , and during the slow-roll inflation is satisfied the condition ǫ ≃ M 2 P /2(∂ φ V /V ) 2 ≪ 1, 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.
The scalar metric perturbations in a FRW background space-time are generically described by the line element: Since we will focus on a joint quantization of the metric and matter perturbations, it is convenient to work with the gauge-invariant quantity known as the Bardeen potential [41] defined as Φ ≡ φ + 1 a [a(B − E ′ )] ′ . The matter sector dominated by the inflaton, is separated in an homogeneous part plus small perturbations ϕ(η, x) = ϕ 0 (η) + δϕ(η, x). In a similar way, the perturbations of the inflaton can be modeled by the gauge invariant fluctuation of the scalar field δϕ (GI) (η, x) = δϕ + φ ′ 0 (B − E ′ ). The two objects Φ and δϕ (GI) can be used to form a new quantity called the Mukhanov-Sasaki variable [38] i.e., Written in a gauge invariant way, and in the absence of anisotropic stress, the components 00 and 0i of the Einstein perturbed equations at linear order δG ab = 8πGδT ab , can be combined to yield where z ≡ aφ ′ 0 /H. Equation (3) is expressed in terms of gauge-invariant quantities. Nevertheless, in the longitudinal gauge, Φ represents the curvature perturbation and is related to v exactly in the same way as described in (3). Under the slow-roll approximation, it leads us to the following useful expressions: z ′ /z ≃ (1+2ǫ−δ)/(−η) and a ′ /a ≃ (1 + ǫ)/(−η), similarly z ′′ /z ≃ (2 + 6ǫ − 3δ)/η 2 and a ′′ /a ≃ (2 + 3ǫ)/η 2 .

B. Quantization of the perturbations
As mentioned in the Introduction, the CSL model is based on a non-linear modification to the Schrödinger equation. Therefore, it is convenient to begin by presenting the theory of the inflaton in the Schrödinger picture, where the relevant theory objects are the Hamiltonian and the wave functional. One starts with the action of a scalar field, i.e. the inflaton, minimally coupled to gravity, and then by expanding up to second-order in the scalar perturbations one can express the action in terms of the Mukhanov-Sasaki variable v(η, x) [42]. The resulting action, expressed in Fourier modes of the field v(η, x), reads: Since where the indexes R, I denote the real and imaginary parts of v k and p k . We now promote v k and p k to quantum operators, by imposing canonical commutations re- In the Schrödinger picture, the wave functional Ψ[v(η, x)] characterizes the state of the system. Moreover, in Fourier space, the wave functional can be factorized into mode 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 , with the Hamiltonian provided by (6). The usual assumption is that at an early time τ (i.e. the onset of inflation), the modes are in their adiabatic ground state, which is a Gaussian centered at zero with certain spread. Since the initial quantum state is Gaussian, its form is preserved during the time evolution. For reasons that will be evident in the following, it is convenient to work in the momentum representation; thus, the Gaussian state evolves according to Schrödinger equation, with initial conditions given by A k (τ ) = 1/2k, B k (τ ) = C k (τ ) = 0 corresponding to the Bunch-Davies vacuum, which is perfectly homogeneous and isotropic in the sense of a vacuum state in quantum field theory.
In this work, the main reason for invoking the selfinduced collapse of the wave function is precisely to break the homogeneity and isotropy of the initial state, and that the emergence of the seeds of cosmic structures can be achieved. We will accomplish this goal by using the CSL collapse mechanism, and in the following subsection we will provide a very brief general description of it, for a complete review of the CSL model see for instance [32].

C. CSL collapse mechanism
The CSL collapse mechanism is based on a stochastic non-linear modification of the Schrödinger equation. Ideally, this modification induces a collapse of the wave function toward one of the possible eigenstates of an op-eratorΘ, 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. (9) implies that the most probable state will be the one with the largest norm. From (8) and (9), it can be derived the evolution equation of the density matrix operatorρ, i.e.
The density matrix operator can be used to obtain the ensemble average of the expectation value of an operator Ô = Tr[Ôρ]. Consequently, from (10) it follows that:

D. CSL and inflation
In order to apply the CSL to the inflationary regime, we need to establish which are the appropriate observables that emerge from the quantum theory of inflation. A reasonable observable is the curvature scalar perturbation, which is directly related to the temperature anisotropies of the CMB, and in the longitudinal gauge corresponds exactly to the Bardeen potential Φ. A quantization of the Mukhanov-Sasaki variable v yields automatically a quantization of Φ; therefore, the question that arises here is: what is exactly the relation between the quantum and classical objects? In particular, what is the relation betweenΦ and Φ?
We will adopt the point of view that the classical characterization of Φ is an adequate description if the quantum state is sharply peaked around some particular value. In consequence, the classical value corresponds to the expectation value ofΦ [19]. More precisely, the CSL collapse mechanism will lead to a final state such that the relation is valid. Also, equations (3) and (12) allow us to relate each mode of the curvature perturbation to the quantum field variables. This is, where we used the definition of the momentum provided by (5). Equation (13) relates the classical curvature perturbation with the momentum variable of the quantum field, consequently, this strongly suggest that the collapse operatorΘ to be considered in the CSL mechanism is the momentum operator. Thus, if the initial state is the Bunch-Davies vacuum state ], then 0|p R,I k |0 = 0, and as a consequence of (13), the curvature perturbation is Φ k = 0, i.e. the space-time is perfectly homogeneous and isotropic. It is only after the state has evolved, according to the CSL mechanism, that generically p R,I k = 0 and the curvature perturbation is born. This illustrates how the selfinduced collapse provided by the CSL model can generate the primordial perturbations.
Furthermore, from the previous discussion, it is evident that we will adopt the CSL mechanism for each mode of the field and the corresponding real and imaginary parts. Therefore, the evolution of the state vector characterizing each mode of the inflaton field, written in conformal time, will be with H R,I k given in (6). Moreover, with the previous considerations, the motion equation (11) for the ensemble average of the expectation value of an operator is: This completes our treatment of the CSL model during inflation.

E. The scalar power spectrum
Having established the relation between the objectsΦ and Φ, we now focus on the scalar power spectrum. The scalar power spectrum in Fourier space is defined as where P s (k) is the dimensionless power spectrum. The bar appearing in (16) denotes an ensemble average over possible realizations of the stochastic field Φ k . In our approach, the realization of a particular Φ k is given by the self-induced collapse that results from the CSL mechanism.
where the expectation values are being evaluated at the (evolved) state provided by (14). More explicitly, where in the second line we assumed that the CSL model does not induce modes correlations. Thus, from now on we will focus on calculate the quantities p R,I k 2 . In particular, we will calculate only the real part since the computation for the imaginary part proceeds in the same fashion. In order to simplify the notation, we will omit the index R unless it can create confusion, in which case we will write it explicitly.
By using the Gaussian wave function in the momentum representation (7), and the probability associated to W (η) (9), it can be shown that [39], This is, (4Re[A(η)]) −1 is the standard deviation of the squared momentum. It is also the width of every packet in momentum space. Thus, to calculate p k 2 , we only need to find the two terms on the right hand side of (18). The second term will be found from the CSL evolution equation (14), and the first one by using equation (15). Let us focus on the second term. Using the general Gaussian state in momentum space (7) and the CSL evolution equation (14) for the wave function, the motion equation for A k (η) results The previous equation can be solved by performing the change of variable A k (η) ≡ f ′ (η)/[2ik 2 f (η)], resulting in a Bessel differential equation for f . After solving such a equation, and returning to the original variable A k , we obtain being q 2 ≡ k 2 (1−2iλ k ) and where the initial condition for the Bunch-Davis vacuum A k (τ ) = 1/2k was used (recall that τ corresponds to the onset of inflation, thus, τ → −∞). J νs corresponds to a Bessel function of the first kind of order ν s = 1/2 + 2ǫ − δ. Note that equation (20) is exact. However, from the observational point of view, one is interested in modes that are well outside the Hubble radius during inflation, i.e. modes with k ≪ aH or equivalently, −kη → 0. Therefore, we can expand A k (η) in the limit −qη → 0 (provided that λ k ≪ 1): with Γ(ν) the Gamma function. From (21), it is straightforward to obtain Next, we focus on the first term of (18), i.e. p 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 . Thus, the evolution equations for Q, R and S are obtained using (15): (23) 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: where 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, S(τ ) = 0. Equation (24) is exact; expanding it again around −kη → 0 yields with β k ≡ −kτ − ν s π/2 − 3π/4. Also, recall that R(η) ≡ p 2 k . At this point, we have the two terms of (18), i.e. p 2 k from equation (25), and 1/(4Re[A k (η)]) from equation (22). Then, we can write where we have defined Equation (26) is valid for both, the real and the imaginary part ofp k . Therefore, substituting (26) into (17), and taking into account (13), we obtain: Finally, with the definition of the power spectrum (16), and the result in (28), the scalar power spectrum within the CSL model, results 1 We will discuss some physical implications from our prediction in Section IV.

III. THE CSL MODEL AND THE TENSOR POWER SPECTRUM
Once we have successfully applied the CSL model to the primordial scalar perturbations, we now proceed to focus on the tensor perturbations. As it is well known, these perturbations represent gravitational waves characterized by a traceless, transverse and symmetric tensor field. These properties imply that the gravitational waves are polarized in two ways. As it is traditional, we will consider one type of polarization and at the end we will just multiply by a factor of 2 the final results.
The action for the tensor perturbations is obtained from the Einstein-Hilbert action by expanding the tensor perturbations h ij (x, η) up to second-order [42]. The resulting action for the tensor field h ij (x, η) can be expressed in terms of its Fourier modes h ij (k, η) = h k (η)e ij (k), with e ij (k) representing a time-independent polarization tensor. Performing the change of variable 1 Note that in Eq. (28) the slow-roll parameter ǫ appears in the numerator, while in the expression for the scalar power spectrum (29) appears in the denominator. The reason for this difference is that, in the longitudinal gauge, the scalar curvature perturbation Φ becomes amplified by a factor of 1/ǫ during the transition from inflation to the radiation dominated stage [43,44], in which the CMB is created (the decoupling era). Therefore in order to obtain a consistent prediction to be compared with the observational data, we must multiply by a factor of 1/ǫ 2 the scalar power spectrum obtained during inflation associated to Φ k Φ ⋆ k ′ .
the action can be written as: By comparing the action in (4) with the one obtained for the tensor perturbations (31), we see that they are equivalent as expected. The only difference is that the action for the scalar perturbations contains a term z ′ /z while the action for the tensor modes contains the term a ′ /a. Therefore, the quantization procedure of Section II B remains exactly the same by simply replacing z ′ /z → a ′ /a. Now, as it was done in Section II D, we need to define which are the appropriate observables that emerge from the quantum theory of the tensor modes. By the same argument that allowed us to identify Φ k = Φ k , we will assume that is also valid. Hence, it is evident that a quantization of v k from the action (31) yields a quantization ofĥ k . In other words, equations (30) and (32) imply that: Thus, (33) relates the quantum field variablev k to the amplitude of the tensor mode h k . Furthermore, this relation is analogous to (13), in which the scalar curvature perturbation is related to the momentum field variablê p k .
On the other hand, as a matter of physical consistency in the operation of the CSL model, we will keep the assumption that the collapse operator corresponds to the momentum operator, even if the relation between h k and the quantum matter fields is in terms of v k . In other words, the quantum theory described by actions (31) and (4) are physically equivalent, i.e. they correspond to the action of a scalar field with time dependent mass. Thus, the fields 'do not know' how they will be related to the curvature perturbation. We will deepen this discussion in the next Section. Therefore, since we retained the momentum operator as the collapse operator, the treatment of the CSL model applied to inflation, remains the same as the one presented in (14) and (15). Moreover, given that the quantization procedure of Section II B stayed unchanged, we will continue working with the wave function shown in (7).
With the relation betweenĥ k and h k already established, we can now focus on the dimensionless tensor power spectrum defined as Furthermore, by using (33) it results in Proceeding in a similar manner with which we arrive at (18), we now find that v R,I The second term of this expression is obtained from Eq. (21), this is, in the limit −kη → 0 we can write, being ν t = 1/2+ǫ. This is different from ν s = 1/2+2ǫ−δ. This distinction arises because we have replaced z ′ /z → a ′ /a in the whole computation.
Next, in order to obtain the first term in (36) (v R,I k ) 2 , we note that its evolution equation is already shown in the system of differential equations given by (23), recalling that Q ≡ (v R,I k ) 2 (and also that z ′ /z → a ′ /a). The solution for Q can be found analytically. However, since we are interested in the modes such that −kη → 0, then Q(η) becomes with β k being the same as in (25) but replacing ν s → ν t . By using (37) and (38) into (36) Consequently, equation (39) along with (35) imply that Now, by comparing (40) with the definition of the tensor power spectrum (34), we finally obtain where we have multiplied by a factor of 2 due to the polarization of the gravitational waves, and used that a(η) ≃ −1/(Hη). With the expressions of the scalar and tensor power spectrum at hand, Eqs. (29) and (41) respectively, it is now straightforward to calculate the tensor-to-scalar ratio defined as r ≡ P t (k)/P s (k). Under the approximation ν s ≃ ν t ≃ 1/2, it leads to r ≃ 16ǫ, which is exactly the same prediction as in the traditional inflationary scenario. This result is also consistent with the one obtained in [45] in which the tensor-to-scalar ratio remained unchanged when considering a generic collapse scheme.

IV. DISCUSSION
With the power spectra within the CSL framework calculated, we will make a few remarks. Let us rewrite the scalar and tensor power spectra in the following suggesting forms: (43b) The quantities, A s and A t correspond to the amplitude of the scalar and tensor power spectrum respectively. From the definition of A s , it can be shown [46] that this amplitude is practically a time-independent quantity [i.e. d/dη{A s } = O(ǫ 2 , δ 2 )]. The same argument applies for A t . Actually, the amplitudes A s and A t coincide exactly with the ones from standard predictions. Moreover, since the amplitudes are constants, it is customary to evaluate the power spectra at some conformal time η * , which is usually taken to be the conformal time at the horizon crossing −kη * = 1. However, from our point of view, the self-induced collapse of the wave function generates the primordial perturbations. Therefore, we cannot evaluate the power spectra at some conformal time, e.g. at the horizon crossing, if the collapse has not taken place yet (or more precisely the CSL mechanism has not concluded), but this evaluation is purely conventional, our prediction for the power spectra is time independent.
The shapes of the power spectra (this is, its dependence on k) are of the form k ns−1 F (λ k , ν s ) and k nt F (λ k , ν t ) for the scalar and tensor power spectrum respectively. Let us focus on the function F (λ k , ν) shown in (27) (where ν = ν s , ν t ), and make the approximation ν s ≃ ν t ≃ 1/2, i.e.
where we used the definitions of β k , ζ k and θ k . Furthermore, if the following assumptions are valid i) −kτ ≫ 1 and ii) λ k ≪ 1, then Henceforth, if λ k = λ 0 /k, then F becomes essentially scale-invariant and the only dependence on k in the power spectra is of the form k ns−1 and k nt . In fact, from (43) and the definitions of ν s and ν t , we have n s −1 = −4ǫ+2δ and n t = −2ǫ. This is, the spectral indexes are exactly the same as in the standard approach. Thus, if λ k = λ 0 /k we recover the standard prediction, what was also noted in [39] but within the semiclassical gravity and just for the scalar case. In Section V A we will say more about the validity of i) and ii). Finally, it is interesting to observe that if λ k = 0, i.e. the vacuum state is evolving according to the unmodified Schrödinger equation and no self-induced collapse occurs, then F (0, ν) = 0 (with ν = ν s , ν t ). Consequently, P s = P t = 0. This implies that there are no primordial perturbations and the space-time is perfectly homogeneous and isotropic. On the other hand, the primordial perturbations are generated by "switching on" the self-induced collapse λ k = 0. In this case, Φ k is created randomly, and its power spectrum is obtained from Φ k Φ * k ′ . Note that one could argue that maybe the cancellation in the case λ k = 0 occurs only in the limit −kη → 0, i.e. only affects the super-horizon modes. We have checked that even if one considers the exact analytic expressions of p 2 k and 1/(4Re[A k (η)]), as well as (v R,I k ) 2 and |A k (η)| 2 /(Re[A k (η)]) in the tensor case, one arrives at the same result, this is, the space-time does not contain perturbations of any scale.

A. CSL parameter
We begin this Section by giving a rough estimate for the value of the CSL parameter λ k . As mentioned in the previous Section, if assumptions i) and ii) are valid, then by taking λ k = λ 0 /k the standard prediction for the power spectra is recovered. Assumption i) reads −kτ ≫ 1. An estimate for τ can be obtained by assuming that the energy scale at the onset of inflation is 10 16 GeV, and that inflation lasts ∼ 70 e-folds. Hence, τ ≃ −10 7 Mpc. Since the observational range for k is 10 −6 Mpc −1 ≤ k ≤ 10 −1 Mpc −1 , then assumption i) is verified. On the other hand, assumption ii) is λ k = λ 0 /k ≪ 1. Thus, by using again the range for k, the CSL parameter must satisfy λ 0 ≪ 10 −6 Mpc −1 . We can check that this is indeed the case form (42) along with Friedmann equation. This yields P s (k) ∝ V M 4 P ǫ (−λ 0 τ ) and as a consequence, if λ 0 ≃ −1/τ , our model prediction is consistent with the standard accepted result. Moreover, since τ ≃ −10 7 Mpc, then λ 0 ≃ 10 −7 Mpc −1 ≃ 10 −21 s −1 , in which case, assumption ii) is valid.
Thus, if the CSL parameter is of the form λ k = λ 0 /k, with λ 0 ≃ 10 −21 s −1 , the model is compatible with the CMB observational data. This result is also consistent with findings of the CSL collapse mechanism applied to situations different from the cosmological context. In particular, the fact that λ k depends on the wave number k captures the spirit of the 'amplification mechanism' of the CSL model, which states that micro and macroobjects do not behave the same way. For example, microobjects tend to exhibit quantum features (e.g. superposition in their wave functions) while macro-objects do not. In our model, the particular dependance on k of the CSL parameter (which generically sets the strength of the self-induced collapse) implies that the self-induced collapse affects the large-scale modes more than the lowscale modes; in some sense, the low-scale modes behave quantum mechanically while large-scale modes can be treated classically. It is also interesting to compare our value λ 0 ≃ 10 −21 s −1 with the one suggested by the GRW model (λ GRW = 10 −16 s −1 ), adopted later in the CSL model with the mass-density operator as the collapse operator [30,31]. Our estimated value for λ 0 makes the strength of the self-induced collapse five orders of magnitude weaker than the one in the GRW model. However, our estimate was based on very robust assumptions regarding the inflationary era, particularly, the energy scale of inflation and the number of e-foldings. Note that the value of the parameter ǫ can also change the value of λ 0 by a couple of orders of magnitude. Thus, a more careful analysis including the full angular power spectrum of the CMB data is required in order to constrain the value of the CSL parameter. This subject is left for a future work.

B. Squeezing of the modes
It is usually argued that during inflation, the Bunch-Davies vacuum evolves towards a squeezed state; this is, a state in which its uncertainty is increased in one variable and decreased in another one, such that the product of the uncertainties satisfy the minimum value allowed by the Heisenberg uncertainty principle. In the Schrödinger picture, the uncertainties of the real and imaginary parts ofp k andv k (associated to scalar perturbations) are given by Henceforth, by using (22) and (37), we see that ∆ 2v R,I k (η) ∼ (−kη) −2νs−1 and ∆ 2p R,I k (η) ∼ (−kη) −2νs+1 as −kη → 0. Thus, the uncertainty inv R,I k increases while the uncertainty inp R,I k decreases as inflation takes place. This is an expected result; in fact, it is normally used to show that inflation induces squeezing in the momentum variable [7,8,35,36]. Furthermore, recall that in our approach, the self-induced collapse gen-erates the primordial curvature perturbation. In the longitudinal gauge, our point of view was reflected in Eq. (13); therefore, the uncertainty in the momentum and in the curvature perturbation operator Φ are directly related, and since during inflation ∆ 2p R,I k (η) decreases, then ∆ 2Φ R,I k (η) also decreases. On the other hand, if we had chosen the comoving gauge, the curvature perturbation would be given by R = v/z, where z ≃ a. Consequently, in our approach, there would be a relation between the uncertainties ∆ 2R R,I k (η) = ∆ 2v R,I k (η)/a 2 , and since during inflation a(η) ≃ −1/η then ∆ 2R R,I k (η) ∼ (−kη) −2ν+1 . Thus, it decreases in the same fashion as ∆ 2p R,I k (η). The tensor modes follow this behavior, i.e. sinceĥ k ≃v k /a, then the uncertainty inĥ R,I k decreases even if the uncertainty inv R,I k increases. One could argue that this behavior occurs irrespectively if there is a quantum collapse or not. However, it should be borne in mind that, in our approach, if there is no quantum collapse then the classical curvature perturbation is exactly zero in spite of its quantum uncertainty ∆ 2Φ R,I k (η) (in the longitudinal gauge) or ∆ 2R R,I k (η) (in the comoving gauge) decreases as inflation takes place. We should also mention that even if the Bardeen potential Φ, the Mukhanov-Sasaki variable v and R are gauge-invariant quantities, the curvature perturbation, characterized by the spatial Ricci scalar, is not [46]. Hence, when a physical interpretation of these mathematical entities is needed, one must choose a specific gauge: the Bardeen potential represents the curvature perturbation in the longitudinal gauge, the Mukhanov-Sasaki variable represents the inflaton perturbation in a spatially flat gauge, and the quantity R represents the curvature perturbation in the comoving gauge.
In summary, our conceptual point of view regarding the collapse states the following: if there is no quantum collapse (λ k = 0), then Φ k = Φ R k + iΦ I k ∝ 0|p R k |0 + i 0|p I k |0 = 0, with the uncertainty inΦ R,I k decreasing as (−kη) −2νs+1 , and therefore the space-time is perfectly homogeneous and isotropic. It is only by introducing a self-induced collapse, provided in this case by the CSL mechanism, that λ k = 0, then Φ = Φ R k + iΦ I k ∝ p R k + i p I k = 0 and the uncertainty of the curvature perturbation operator decreases as ∆ 2Φ R,I k (η) ≃ (−kη) −2νs+1 /[ζ k cos(πθ k )]. Also, note that since we are treating the self-induced collapse as a sort of 'measurement process', then it must act on the hermitian operators, sayΦ R k andΦ I k . Thus, the uncertainties that are important from the 'measurement like process' are those corresponding to ∆ 2ΦR k and ∆ 2ΦI k .
are no curvature perturbations at all. As a consequence, P s = 0 = P t . This difference ultimately boils down to the quantum objects that are going to be identified with the power spectrum. In the framework of [35,36], the power spectrum corresponds exactly to the Fourier transform of the two-point quantum correlation function, this is, P s (k) ≃ Ψ|v kvk ′ |Ψ , where the state |Ψ corresponds to the post-collapse state, or to the Bunch-Davies vacuum if λ = 0. In our approach, we have assumed that, once the collapse took place, then it can be assumed that Φ k = Ψ|Φ k |Ψ . Therefore, the quantum collapse is the mechanism by which the curvature perturbation emerges. After the curvature perturbation was generated, as a result of a stochastic process, one can calculate the power spectrum associated to a classical stochastic field, i.e.
Another aspect that the authors in [35,36] share is that they treat the field variable v as their collapse operator, and then argue why it would be a desirable feature to obtain a squeezing in that variable. From our point of view, within the cosmological context, the choice of the collapse operator as the field variable, is equally valid as the choice of the momentum operator as we have done in the present paper. In fact, we motivated this choice by Eq. (13), but we could have chosen to work in the comoving gauge, in which R = v/z represents the curvature perturbation, and then have argued to set v as the collapse operator. We think that in the absence of a complete relativistic version of the CSL mechanism, it is not entirely clear which operator plays the role of the collapse operator. On the other hand, we also think that the squeezing of a variable does not necessarily denote that the system has become classical. For instance, squeezed states are regarded as highly non-classical states in Quantum Optics [47]. Therefore, in our approach neither the squeezing of the momentum nor the squeezing of the field variable are crucial for the classicalization of the perturbations. Moreover, as we have argued, if the squeezing occurs in the momentum variable, which is also our choice for the collapse operator, then by (13), the "squeezing" also occurs in the curvature perturbation operator (because the longitudinal gauge was selected). Additionally, (13) allows us to establish the relation P s (k) ∝ Ψ|p k |Ψ Ψ|p * k ′ |Ψ . Likewise, here the tensor modes are directly related to the field variable through (33), and in this case the collapse operator is still the momentum operator. Thus, even if the uncertainty inv R,I k increases, the uncertainty inĥ R,I k (which is associated to the tensor curvature perturbation), does not become larger. This is because of the relation ∆ 2ĥ R,I k ≃ ∆ 2v R,I k /a 2 . In both, the scalar and the tensor cases, we have focused on the uncertainties of the quantum operators that are directly related to the observables of the CMB; namely, the operators associated to the curvature perturbations.
Finally, it is also worthwhile to mention that the conceptual point of view adopted in this paper, regarding the role of the collapse during inflation, was the same as the one followed in [39]. However, in such a work, the authors applied the CSL model to inflation within the semiclassical gravity framework, and one important difference with our work is that, within their approach, the amplitude of the tensor modes is exactly zero at first-order in the perturbations. Additionally, the authors in [39] did not consider a full quasi-de Sitter background space-time and therefore their final expression for the scalar power spectrum does not contain a specific prediction for the scalar spectral index.

VI. CONCLUSIONS
In this work, we have adopted an important role for the self-induced collapse of the wave function characterized by the CSL model, i.e. the collapse acts as the main agent for the emergence of the primordial curvature perturbations. By focusing on the curvature perturbation operators, we have presented a prediction for the scalar and tensor power spectra (42) that are quite consistent with the standard prediction. In particular, the scalar and tensor spectral indexes are exactly the same, as well as the scalar A s and tensor A t amplitudes. The tensor-toscalar ratio is also of the same form as in the traditional approach. On the other hand, the addition of the CSL mechanism, is reflected in the appearance of the function F (λ k , ν) [defined in (27)] multiplying the standard prediction for the scalar and tensor power spectra. However, if ν ≃ 1/2 and the CSL parameter is of the form λ k = λ 0 /k, then F (λ k , ν) becomes practically independent of k, and we recover the standard prediction for the power spectra. This allows us to give a roughly estimate for the CSL parameter λ 0 ≃ 10 −21 s −1 , which can be refined by using the full CMB data set, what is left for a future work. On the other hand, small variations of the recipe λ k = λ 0 /k will result on different predictions from the standard case and that could be confronted with the observational data.
Our prediction for the power spectra differs from the ones obtained in previous works [35][36][37] in which the CSL mechanism was also applied to the inflationary universe. In particular, in contrast to Ref. [35], we have captured the amplification mechanism by assuming a k dependence of the CSL parameter. Additionally, our final result for the power spectra is not time dependent and we have only one free parameter λ 0 to be adjusted using the observational data. This is different form the findings of Ref. [36], in which their final predictions for the power spectra depend on the conformal time and their model introduces more than one free parameter. The reasons for these differences are traced back to the conceptual point of view regarding the collapse of the wave function during inflation. Specifically, if λ k = 0, within the framework of Refs. [35,36], then their model recovers the standard prediction. On the other hand, in our approach if λ k = 0, i.e. there is no quantum collapse. Then, the tensor and scalar power spectra are exactly zero, which means that the space-time does not contain any perturbations at all. Another important (conceptual) difference from previous works is that the squeezing of the field and momentum variables seems to play a relevant role in the approach of Refs. [35][36][37]. In our approach, we have shown that the uncertainty of the curvature perturbation operator (either in the comoving or in the longitudinal gauge) decreases as inflation goes on; however, the self-induced collapse is still needed for generating the primordial perturbation.
Our result, unlike the case of semiclassical gravity, predicts non-null tensor modes at first-order in perturbations, and it is also consistent with the findings of [45], where the tensor spectrum and the tensor-to-scalar ratio were calculated under a generic collapse scheme.
We conclude that, from the phenomenological point of view, the CSL mechanism can successfully be used to explain the origin of the primordial perturbations. However, there are still some aspects concerning the CSL model that need to be addressed further. For exam-ple, it is well known that the CSL model violates energy conservation, which can lead to divergences in the energy-momentum tensor. Also, the CSL model is nonrelativistic, hence, the choice of the collapse operator and the k dependence in the λ k parameter is at this point purely phenomenological. Therefore, we think that the early universe provides a natural laboratory to test new ideas that would allow us to reach a deep understanding of Nature.