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 one 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.


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, a e-mail: gleon@df.uba.ar b e-mail: gabriel@iafe.uba.ar 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 the technical level. Therefore, it is expected that when these theories are used in the same footing, some difficulties might arise. Specifically, the quantum-to-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.
A more precise formulation of the problem at hand can be stated as follows: the inflationary universe is characterized by a background spacetime that is completely homogeneous and isotropic. Equivalently, the quantum state of the matter field, namely the vacuum state, is also perfectly homogeneous and isotropic (i.e. it is an eigenstate of the operators associated to the generators of spatial translations and rotations). On the other hand, the present universe is the result of the evolution of primordial density inhomogeneities.
Henceforth, the problem is: how do the inhomogeneities and anisotropies originate from the initial highly symmetric (the symmetry being the homogeneity and isotropy) state of the universe, described by both the background spacetime and the vacuum state of the matter fields, given that the dynamical evolution, provided by the Schrödinger and Einstein equations, does not break the translational and rotational symmetries? In other words, if one considers quantum mechanics as a fundamental theory applicable, in particular, to the universe as a whole, then one must regard any classical description of the state of any system as a sort of imprecise characterization of a complicated quantum-mechanical state. The universe we observe today is clearly well described by an inhomogeneous and anisotropic classical state; therefore, such description must be considered as an imperfect descrip-tion of an equally inhomogeneous and anisotropic quantum state. Consequently, if we want to consider the inflationary account as providing the physical mechanism for the generation of the seeds of structure, such account must contain an explanation for why the quantum state that describes our actual universe does not possess the same symmetries as the early quantum state of the universe, which happened to be perfectly symmetric. Since there is nothing in the dynamical evolution (as given by the standard inflationary approach) of the quantum state that can break those symmetries, the traditional inflationary paradigm is incomplete in that sense (sometimes, an usual argument is that since the vacuum fluctuations are not zero, somehow the system contains inhomogeneities. See Appendix A: for a discussion). Here, we also refer the reader to Refs. [7,8] where the relation between symmetry-breaking vacuum and the reduction of the wave function is discussed.
Earlier works based on decoherence [9][10][11][12] and the consideration that the initial vacuum state of the universe evolves into a highly squeezed state [13] led to a partial 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 [14]. Nevertheless, this argument by itself cannot address the fact that a single (classical) outcome emerges from the quantum theory. In other words, decoherence (and the squeezing of quantum states) cannot solve the quantum measurement problem [15,16], 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. Other cosmologists seem to adopt the Everett "many-worlds" interpretation of quantum mechanics plus the decoherence process when confronted with the quantum-toclassical transition in the inflationary universe [17]. In the Everettian formulation, reality is made of a connected weave of ever-splitting worlds, each one realizing one of the alternatives that is opened by what we would call a quantummechanical measurement. Regarding this point, we would like to refer the reader to Refs. [18][19][20] where arguments against the Everett interpretation are presented.
One possible way to address the mentioned issues or, in other words, to avoid the standard "measurement problem" of quantum theory, is to invoke the collapse of the wave function but without relying on any external objects, e.g. [18,21]. Other approaches to this problem have been based on Bohmian quantum mechanics; see for instance [22,23]. 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. [24][25][26][27][28][29][30][31][32][33]. Moreover, in [34] some generic collapse schemes have been tested using observational data coming from the 7-year release of the WMAP collaboration [35] and the Sloan Digital Sky Survey [36]. 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 Diósi [37,38] and Penrose [39], advocating gravity as the main agent triggering the collapse. Also, the Ghirardi-Rimmini-Weber (GRW) [40] 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 [41,42]. In particular, the CSL mechanism is based on a nonlinear 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 whether they are isolated or interacting. The idea behind the CSL model [43], sometimes referred to as the "amplification mechanism", is that the collapses must be rare for microscopic systems, in order not to 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 [44]. On the other hand, the CSL model is, at this stage, a non-relativistic model and therefore a complete generalization to quantum fields is still under development [45,46]. Nevertheless, the particular features of the CSL model, e.g. the absence to rely on external agents to solve the measurement problem, would make it a viable candidate to address the issue concerning the emergence of the classical primordial perturbations if it predicts a spectrum consistent with the observational data.
In [47], the CSL model was applied to the inflationary universe for the first time. In this reference, 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 scale-invariant power spectrum, which is consistent with the CMB data. In order to retain the amplification mechanism, other authors [48], 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 ) [49].
A shared feature of the works in Refs. [47,48] is that the authors worked in a joint metric-matter quantization of the perturbations characterized by the Mukhanov-Sasaki variable [50]. On the other hand, in Ref. [51] one followed a different approach to the problem; the authors successfully applied the CSL collapse mechanism to the inflationary universe but within the semiclassical gravity framework. Consequently, as pointed out in [30,52], 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 make this approach face serious issues.
Another distinction between the approaches in Refs. [47,48] and [51] is the specific role played by the self-induced collapse. This difference is subtle but important. The authors of [47,48] employed the CSL mechanism with the intention of localizing the evolved inflaton vacuum state in an eigenstate of the Mukhanov-Sasaki variable. Meanwhile, the authors of [51] attributed to the self-induced collapse the action of generating the primordial curvature perturbation, irrespectively of which eigenstate the initial state evolves into. Therefore, if no self-induced collapse occurs, the predicted power spectrum in Refs. [47,48] is exactly the same as the standard one. On the contrary, in Ref. [51] the absence of a collapse results in no perturbations of the spacetime at all, consequently, the predicted scalar power spectrum is exactly zero (see Sect. 2.4 and Appendix B: for a more detailed discussion).
In the present work, we adopt the same role for the collapse, provided by the CSL model, as the one presented in [51], but we apply it to the joint metric-matter quantization of the inflaton field, to improve the standard inflationary treatment providing a mechanism for the emergence of the primordial curvature perturbations, and also for driving the quantum-to-classical transition. In this way, we obtained predictions that are different from the ones in [47,48] 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. [51], it predicts a non-vanishing tensor power spectrum at first order in the perturbations.
The article is organized as follows: in Sect. 2 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 Sect. 3 we show our results for the power spectrum of tensor modes and the tensor-to-scalar ratio; in Sects. 4 and 5 we make a discussion of our results and compare them with previous works, and finally in Sect. 6 we summarize our conclusions.

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.

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 spacetime 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 2 1, with M 2 P ≡ (8π G) −1 the reduced Planck mass. Since we will work in a full quaside 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 spacetime 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 [53] . 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 [50] i.e., v(η, x) ≡ a δϕ (GI) 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 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 .

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) [54]. The resulting action, expressed in Fourier modes of the field v(η, x), reads δ (2) The canonical conjugated momentum associated to v k is p k = ∂L/∂v k , i.e. Since where the indices 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 relations 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 com- . 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 (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 become 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 (see Appendix A:).
In this work, the main reason for invoking the self-induced collapse of the wave function is precisely to break the homogeneity and isotropy of the initial state, and as a consequence, 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 [43].

CSL collapse mechanism
The CSL collapse mechanism is based on a stochastic nonlinear modification of the Schrödinger equation. Ideally, this modification induces a collapse of the wave function toward 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 continuoustime 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

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 ?
To illustrate the point of view we will adopt, let us focus on the temperature anisotropies of the CMB observed today on the celestial two-sphere and its relation to the scalar metric perturbation . 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 expansion δT /T 0 (θ, ϕ) = lm a lm Y lm (θ, ϕ), hence here θ and ϕ are the coordinates on the celestial two-sphere, with Y lm (θ, ϕ) as the spherical harmonics. Given Eq. (12), the coefficients a lm can be further reexpressed in terms of the Fourier modes associated to , i.e.
with j l (k R D ) the spherical Bessel function of order l of the first kind and R D is the comoving radius of the last scattering surface. We have explicitly included the modifications associated with late-time physics encoded in the transfer functions (k). The metric perturbation k is the primordial curvature perturbation. Now, how to relate k , which appears in Eq. (14), with the quantum operatorˆ k coming from the quantum theory? Evidently, if we compute the expectation value ˆ k in the vacuum state |0 and identify it exactly with k , then we obtain exactly 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 (clearly the concepts of measurements, observers, etc. are not well defined in the early universe).
Nevertheless, in the standard approach (by invoking decoherence, squeezing of the vacuum, many-world interpretation of quantum mechanics, etc. although we do not subscribe to such postures for the reasons exposed in Ref. [18]) somehow one can make the association k = Ae iα k with α k a random phase and A being identified with the quantum uncertainty of But the random nature of k , codified in the random number e iα k , remains unclear. In fact, the relation k = Ae iα k is not valid at all times in the standard approach. It is only valid after the proper wavelength of the mode is larger than the Hubble radius (or when "the mode has crossed the horizon"). That is, the traditional statement is that during inflation the modes become super-horizon and then occurs the transitionˆ → . Furthermore, the relation k = 0|ˆ 2 k |0 e iα k could generate some misunderstandings, in the sense that one might think that, since the vacuum fluctuations 0|ˆ 2 k |0 are not zero, k should be different from zero when the mode is super-horizon, and the spacetime perturbations are "born" (see Appendix A: for clarification).
On the other hand, in our approach, the random nature of k will come directly from the stochastic aspects of the quantum dynamical reduction, i.e. from the CSL mechanism. Moreover, 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ˆ [27]. More precisely, the CSL collapse mechanism will lead to a final state such that the relation is valid.
Therefore, in our approach, the coefficients a lm in Eq. (14), will be given by (16) where | corresponds to the evolved state according to the non-unitary modification of the Schrödinger equation provided by the CSL mechanism. At this point, we encourage the reader to consult Appendix B: for a discussion of other possible ways to relateˆ and within the CSL framework. Additionally, Eqs. (3) and (15) allow us to relate each mode of the curvature perturbation to the quantum field variables. That is, where we used the definition of the momentum provided by (5). Equation (17) 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 ], then 0|p R,I k |0 = 0, and as a consequence of (17), the curvature perturbation is k = 0, i.e. the spacetime 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 d dη This completes our treatment of the CSL model during inflation.

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 (20) 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. Henceforth, Eq. (17) implies that k k ∝ p k p k , where the expectation values are being evaluated at the (evolved) state provided by (18). 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 [51], That 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 (22). The second term will be found from the CSL evolution Eq. (18), and the first one by using Eq. (19). Let us focus on the second term. Using the general Gaussian state in momentum space (7) and the CSL evolution Eq. (18) for the wave function, the equation of motion for A k (η) results, where 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 Eq. (24) 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 a H or equivalently, −kη → 0. Therefore, we can expand A k (η) in the limit −qη → 0 (provided that λ k 1): with (ν) the Gamma function. From (25), it is straightforward to obtain where we have defined ζ k e iθ k ≡ √ 1 − 2iλ k . Next, we focus on the first term of (22), i.e. p 2 k . It will be useful to define the 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 (19): 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 (28) 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 (22), i.e. p 2 k from Eq. (29), and 1/(4Re[A k (η)]) from Eq. (26). Then we can write where we have defined Equation (30) is valid for both the real and the imaginary part ofp k . Therefore, substituting (30) into (21), and taking into account (17), we obtain Finally, with the definition of the power spectrum (20), and the result in (32), the scalar power spectrum within the CSL model, results 1 We will discuss some physical implications from our prediction in Sect. 4.

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 1 Note that in Eq. (32) the slow-roll parameter appears in the numerator, while in the expression for the scalar power spectrum (33) 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 [55,56], 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 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 is traditional, we will consider one type of polarization and at the end we will just multiply the final results by a factor of 2.
The action for the tensor perturbations is obtained from the Einstein-Hilbert action by expanding the tensor perturbations h i j (x, η) up to second order [54]. The resulting action for the tensor field h i j (x, η) can be expressed in terms of its Fourier modes h i j (k, η) = h k (η)e i j (k), with e i j (k) representing a time-independent polarization tensor. Performing the change of variable the action can be written as δ (2) By comparing the action in (4) with the one obtained for the tensor perturbations (35), 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 Sect. 2.2 remains exactly the same by simply replacing z /z → a /a. Now, as was done in Sect. 2.4, 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 (35) yields a quantization ofĥ k . In other words, Eqs. (34) and (36) imply that Thus, (37) relates the quantum field variablev k to the amplitude of the tensor mode h k . Furthermore, this relation is analogous to (17), in which the scalar curvature perturbation is related to the momentum field variablep 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 (35) 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 (18) and (19). Moreover, given that the quantization procedure of Sect. 2.2 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 (37) it results in Proceeding in a similar manner with which we arrive at (22), we now find that v R,I The second term of this expression is obtained from Eq. (25), that is, in the limit −kη → 0 we can write where ν 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 (40), i.e. (v R,I k ) 2 , we note that its evolution equation is already shown in the system of differential equations given by (27), 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, Q(η) becomes with β k being the same as in (29) but replacing ν s → ν t . By using (41) and (42) Consequently, Eq. (43) along with (39) imply that Now, by comparing (44) with the definition of the tensor power spectrum (38), we finally obtain where we have multiplied by a factor of 2 due to the polarization of the gravitational waves, and we used a(η) −1/(H η). With the expressions of the scalar and tensor power spectrum at hand, Eqs. (33) and (45), 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 [57] in which the tensor-to-scalar ratio remained unchanged when considering a generic collapse scheme.

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 suggestive forms: where we define 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 [58] 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 (that is, its dependence on k) are of the form k n s −1 F(λ k , ν s ) and k n t F(λ k , ν t ) for the scalar and tensor power spectrum, respectively. Let us focus on the function F(λ k , ν) shown in (31) (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 scaleinvariant and the only dependence on k in the power spectra is of the form k n s −1 and k n t . In fact, from (47) and the definitions of ν s and ν t , we have n s −1 = −4 +2δ and n t = −2 .
That is, the spectral indices are exactly the same as in the standard approach. Thus, if λ k = λ 0 /k we recover the standard prediction, what was also noted in [51] but within the semiclassical gravity and just for the scalar case. In Sect. 5.1 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 spacetime 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. it only affects the superhorizon 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, that is, the spacetime does not contain perturbations of any scale.

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 , 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 (46)  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 macro-objects do not behave the same way. For example, micro-objects tend to exhibit quantum features (e.g. superposition in their wave functions), while macro-objects do not. In our model, the particular dependence 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 low-scale modes; in some sense, the low-scale modes behave quantum mechanically while largescale 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 [41,42]. 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.

Squeezing of the modes
It is usually argued that during inflation, the Bunch-Davies vacuum evolves toward a squeezed state; that 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 (26) and (41), we see that 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 [9,13,47,48]. Furthermore, recall that in our approach, the self-induced collapse generates the primordial curvature perturbation. In the longitudinal gauge, our point of view was reflected in Eq. (17); therefore, the uncertainty in the momentum and in the curvature perturbation operator are directly related, and since during inflation 2p R,I k (η) decreases, 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/η, we have 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, the uncertainty inĥ R,I k decreases even if the uncertainty in v R,I k increases. One could argue that this behavior occurs irrespectively of whether 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 [58]. 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 spacetime 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 . Also, note that since we are treating the self-induced collapse as a sort of 'measurement process', 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 .

Comparison to previous works
As we have mentioned in the Introduction, the CSL model has been applied to the inflationary universe in previous works by several authors, and they have reached different conclusions with each other. Here we will present a brief summary of their results and a comparison to our findings. The first work we will address is the one presented in [47]. There, the authors applied the CSL model to inflation using the Mukhanov-Sasaki variable as the collapse operator, and assumed that the CSL parameter λ to be the same for all modes, i.e. independent of k. Those assumptions led to a scalar power spectrum with two branches, one that is scaleinvariant, and another one with a spectral index n s = 4 in an evident conflict with the CMB data. In order to suppress the conflicting branch, the CSL parameter had to be severely constrained. Furthermore, the fact that λ is independent of k fails to incorporate the amplification mechanism that is crucial in the CSL model. In some sense, if λ is the same for all modes, then small and large-scale modes exhibit a classical behavior. Another issue mentioned by the authors is that the uncertainty in the Mukhanov-Sasaki variable increases even if such a field variable corresponds to the CSL collapse operator.
In order to address some of the issues of [47], other authors proposed that the CSL parameter is of the form λ(k, η) = λ 0 (k/k 0 ) β /(−kη) α (with k 0 a pivot scale), i.e. a function of both, the conformal time and k, while retain-ing the Mukhanov-Sasaki variable as the collapse operator [48,49]. Hence, since the CSL parameter depends on k, the amplification mechanism of the CSL model is recovered. On the other hand, their model introduced three free parameters: λ 0 , α and β. However, by taking into account the tensor modes [49], the modification to the scalar and tensor spectral indices could be captured in just one effective free parameter δ = 3 + α − β and λ 0 . In fact, if δ = 0, then their model predicts a scale-invariant scalar power spectrum. Furthermore, if 1 < α < 2, then, as inflation goes on, the squeezing occurs in the Mukhanov-Sasaki variable rather than its conjugated momentum. This is because, in their approach, the CSL parameter depends on the conformal time. Nevertheless, this dependence on η is inherited in their final prediction for the scalar and tensor power spectra; thus, they choose to evaluate their power spectra at the end of inflation. This choice, in turn, translates into an extremely small value for the CSL parameter λ 0 e −120 , however, a higher value can be achieved by bringing down the energy scale of inflation.
References [48,49] can be considered as an improved version of [47]; nevertheless, they share some features. First, in both approaches, if λ = 0, which means there is no selfinduced collapse, their prediction for the (scalar) power spectrum is exactly the same as in the standard approach. 2 This contrasts with the result presented in this letter, since in our approach, if λ k = 0 there 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 [47,48], the power spectrum corresponds exactly to the Fourier transform of the two-point quantum correlation function, that 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, 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. P s (k) = k * k = |ˆ k | |ˆ * k | . Another aspect that the authors in [47,48] 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. (17), 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 [59,60]. 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 (17), the "squeezing" also occurs in the curvature perturbation operator (because the longitudinal gauge was selected). Additionally, (17) 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 (37), 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 [51]. 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 [51] did not consider a full quasi-de Sitter background spacetime and therefore their final expression for the scalar power spectrum does not contain a specific prediction for the scalar spectral index.

Conclusions
In this work, we have adopted an important role for the selfinduced 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. We are aware that other people in the scientific community might prefer to adopt the Everett "many-world" interpretation of quantum theory plus decoherence in the case of the inflationary universe. However, we have presented the reason that a dynamical reduction of the wave function, provided by the CSL model, can be a viable alternative explanation in the sense that it is consistent with observations and also addresses the standard "measurement problem" of quantum theory.
Specifically, by focusing on the curvature perturbation operators, we have presented predictions for the scalar and tensor power spectra (46) that are quite consistent with the standard prediction. In particular, the scalar and tensor spectral indices are exactly the same, as well as the scalar A s and tensor A t amplitudes. The tensor-to-scalar 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 (31)] 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 (similar to the one suggested by the GRW model, λ GRW = 10 −16 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 in 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 [47][48][49] in which the CSL mechanism was also applied to the inflationary universe. In particular, in contrast to Ref. [47], 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 from the findings of Ref. [48], 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. [47,48], then their model recovers the standard prediction. On the other hand, in our approach, if λ k = 0, there is no quantum collapse. Then the tensor and scalar power spectra are exactly zero, which means that the spacetime does not contain any perturbations at all. Another important (conceptual) difference from previous research is that the squeezing of the field and momentum variables seems to play a relevant role in the approach of Refs. [47][48][49]. 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 the perturbations, and it is also consistent with the findings of [57], where the tensor spectrum and the tensor-to-scalar ratio were calculated within a generic collapse scheme.
We conclude that, from the phenomenological point of view, the CSL mechanism can successfully be used to improve the standard inflationary paradigm by explaining the origin of the primordial perturbations and driving the quantum-to-classical transition. However, there are still some aspects concerning the CSL model that need to be addressed further. For example, 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.
is plausible to think that if λ k << 1, then the scalar spectrum of Refs. [47][48][49] will be consistent with the observations too.