Enlightening the CSL model landscape in inflation

We propose a novel realization for the natural extrapolation of the continuous spontaneous localization (CSL) model, in order to account for the origin of primordial inhomogeneities during inflation. This particular model is based on three main elements: (i) the semiclassical gravity framework, (ii) a collapse-generating operator associated to a relativistic invariant scalar of the energy-momentum tensor, and (iii) an extension of the CSL parameter(s) as a function of the spacetime curvature. Furthermore, employing standard cosmological perturbation theory at linear order, and for a reasonable range within the parameter space of the model, we obtain a nearly scale invariant power spectrum consistent with recent observational CMB data. This opens a vast landscape of different options for the application of the CSL model to the cosmological context, and possibly sheds light on searches for a full covariant version of the CSL theory.


I. INTRODUCTION
One of the most interesting challenges when trying to combine General Relativity with Quantum theory to explain the very early epoch of the universe, has to do with giving a satisfactory explanation of the quantumto-classical transition of perturbations originated during inflation. The proposal that, at the very beginning, the universe went through an accelerated phase, called inflation, is extremely successful since it allowed explaining not only some well-known problems of the hot Big Bang model (e.g. the horizon and flatness problems), but it has also offered a way to address the quantum origin of largescale structures that we see today in the sky, from primordial cosmic seeds [1][2][3][4][5][6][7][8][9][10]. The inflationary paradigm predictions for primordial power spectra are strikingly consistent with the most recent accurate observations of cosmic microwave background (CMB) radiation [11][12][13].
Notice that, the aforementioned quantum to classical transition is closely related to the so-called measurement problem in Quantum physics [14][15][16][17][18][19][20][21] 1 . Let us mention why this problem is notoriously enhanced in the cosmological case [22][23][24][25] 2 . The central point here is that, according to standard Quantum theory, the evolution of any quantum state is always unitary, dictated by the Schrödinger equation, which does not break any initial symmetry of the system or destroy quantum superpositions. If at the beginning of inflation the spacetime is assumed to be spatially isotropic and homogeneous, and the perturbations of matter fields (i.e. the inflaton field) were in a quantum vacuum state also invariant under spatial rotations and translations (i.e. the so-called Bunch-Davies vacuum); then one question arises: how do we arrive at an anisotropic and inhomogeneous state (with cosmic seeds of structures) from a quantum vacuum state that is perfectly isotropic and homogeneous in space? The traditional early universe paradigm seems to be incomplete in that sense 3 . And any model that claims to provide a mechanism for the emergence of the seeds of structure must be able to give a convincing answer to the mentioned scenario, without appealing, of course, to the existence of observers or measuring devices which cannot constitute fundamental elements in the very early universe [28]. There are various proposals that try to explain the quantum-to-classical transition from different approaches. For instance, by invoking the decoherence framework [29][30][31][32], it could provide a partial understanding of the issue. However, it was shown that decoherence does not fully addresses the problem [24,[33][34][35][36][37]. Other authors, within the framework of the de Broglie-Bohm theory [38], have also explored applications to this cosmological puzzle, e.g. [39][40][41][42][43][44].
Known as objective collapse theories, they constitute another approach that seeks to address the aforementioned issue, since they were specifically designed to solve the measurement problem in Quantum Mechanics [45][46][47][48][49][50]. Reviews on these sort of theories can be found, for instance, in [51,52]. The idea of invoking a self-induced collapse in cosmology in order to generate the primordial perturbations and/or accomplish the quantum-to-classical transition has been explored in great detail since 2006 [23], and it has led to numerous investigations in recent years with varied proposals and results, e.g. [24,25,. In several of those works, a specific version based on the continuous spontaneous localization (CSL) model [45,47] adapted to the situation of interest, was chosen to analyze this subject.
The CSL constitutes a particular proposal, based on a non-linear stochastic modification of the standard Schrödinger equation, where spontaneous and random collapses occur, resulting in an objective dynamical reduction of the wave function. For the first time, in Ref. [56] the CSL model was applied to the inflationary universe and in Ref. [57] such exploration was done within the semiclassical gravity framework. The debate on the particular details involved in implementing the CSL model in the cosmological context is still open, and contains an extensive landscape of possibilities constituting an active research area at the moment [80][81][82][83][84][85].
As was already pointed out in [86], objective collapse models have the appealing feature of connecting plausible resolutions of other open problems in a single unified picture; not only the measurement problem in Quantum Mechanics and the emergence of the primordial inhomogeneities during inflation, but also the problem of time in canonical quantum theories of gravitation, the black hole information paradox [87][88][89] and other applications to cosmology, for instance to account for the dark energy and the late-time accelerated expansion of the universe [77,90,91].
Recently, in Ref. [82], a wide variety of open alternatives were discussed when applying the CSL dynamical collapse model to the inflationary era. That exploration included the important aspects one must face in order to address such a problem. These are: (i) the two different approaches to deal with quantum field theory and gravitation, (ii) the identification of the collapse-generating operator, and (iii) the general nature and values of the CSL model parameters. In that same work, it was also emphasized that all the choices connected with those issues have the potential to dramatically alter the conclusions one can draw.
In this article, we will explore a particular CSL model in the inflationary context, within the vast theoretical landscape available for the extrapolation of the standard CSL (as constructed to deal with non-relativistic many particle Quantum Mechanics) into the realms of relativistic quantum field theory in curved spacetimes. To accomplish this task, we will follow some options raised in Ref. [82]. In this way, our model possess three main elements which we mention below.
First, one must choose the setting within which one is going to combine quantum field theory (QFT) with gravitation. While in some previous works, we have explored the use of the CSL using the standard framework of quantization (in which metric and matter fields are quantized simultaneously). Here, as in other previous works, we will adopt the semiclassical gravity framework where grav-ity is always classical while the matter fields are treated quantum mechanically. This approach has received some known criticisms [92,93], but those arguments have been discussed and refuted [94][95][96][97][98][99][100][101]; currently, the theory continues to be of great interest. Moreover, it is not often emphasized that most of the conclusions against are only valid in those contexts where quantum mechanical evolution does not include any sort of measurement-related or spontaneous collapse of the quantum state. We will assume such a framework to be a valid approximation during the inflationary era, which is well after the full quantum gravity regime has ended. This choice is due to the fact that it appears favored from a theoretical and conceptual point of view, in particular, when one wants to incorporate collapse models [54,69,70,82,102]. In fact, it was shown that, within semiclassical gravity, the CSL model predicts generically a very small amplitude of the tensor power spectrum [72,74,75], which is very consistent with the current non-detection of B-modes in the CMB data [12].
Second, as mentioned in Ref. [82], there are a large number of operators that are relativistically invariant and do reduce, in the simple non-relativistic regimes encountered in laboratory situations, to the required "mass density" or "energy density". Some simple examples we might consider are scalars constructed, for example, using the energy-momentum tensor T ab . In the present work, for the first time in this line of research, we will explore the consequences of assuming a collapse generating operator composed of a contraction of the energymomentum tensor. Specifically, we will choose the scalar (T ab T ab ) 1/2 as the collapse operator in our proposal to extrapolate the CSL model into the cosmological context.
Third, we will make a particular choice regarding the parameters of the CSL model. The main idea is the following.
One of the central features of the CSL model, 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. But, on the other hand, their effects must increase when several particles are hold together forming a macroscopic system. The amplification mechanism can be characterized through the collapse rate parameter λ of the CSL model.
As was discussed at length in [82], such a parameter used in the non-relativistic versions of the CSL model, let us call it λ 0 , might not be a fundamental constant after all. The original motivation [45-47, 51, 103] for choosing the particular value of λ 0 ≈ 10 −16 s −1 was based on some phenomena occurring at laboratory scales, e.g. the number of atoms in a pointer apparatus, the matter density of nuclear matter, etc. Therefore, it is not obvious at all that there should be any comparison between the value of λ 0 , informed by the existence of hadrons, nuclei and the atomic densities of solids, with parameters at an era when there were no nucleons, no atoms and certainly no solids.
On the other hand, the second parameter of the nonrelativistic CSL model, r c , is a length scale characterizing the level of de-localization associated with the onset of the spontaneous collapse. For instance, for a single particle, the effect of the collapse aspect of the modified dynamics is negligible if the width of its wave function is smaller than r c . But, if it is larger than r c , the collapse mechanism becomes very relevant. The actual numerical value of r c comes from experimental constraints using a non-relativistic version of the CSL model. Additionally, for that particular CSL model, the parameter r c serves to characterize physical processes at laboratory scales. For example, r c takes into account notions such as atomic dimensions, mean spreads around the equilibrium positions of the lattice points of a crystal, etc. (see Sec. 6.5 of Ref. [51]). Thus, there are ample grounds to doubt any simple one-on-one connection between the value of λ 0 and r c , relevant for one regime (involving laboratory scales), with parameters characterizing the model in a completely different one, such as the inflationary universe.
Previous works [57,62,71,76], where a particular version of the CSL model was applied to inflation, were able to recover the correct shape and amplitude of the CMB spectrum (in addition to being compatible with laboratory constraints). Those results required making suitable assumptions regarding the dependence of the model parameters on the Fourier mode wavelengths, which in turn, might be viewed as tied to a more fundamental dependence on curvature. However, in the present work, we will make explicit this dependence on the curvature in the CSL parameters, and make no assumption involving a particular decomposition used in the field's modes (e.g. Fourier modes).
As a matter of fact, as suggested by L. Diosi and R. Penrose in e.g. [48][49][50]104], it seems quite natural to think that, at a fundamental level, the spontaneous collapse dynamics might be intimately tied with gravitation. Therefore, accepting such a hypothesis, one reasonable option for exploring the relation between the collapse of the wave function and gravitation, is to replace the collapse parameter of the ordinary CSL model by a function of some geometrical scalar associated with curvature (see e.g. [86]). In particular, in a series of works devoted to the examination of the famous black hole information puzzle [87,88], it was found that, during the black hole evaporation, information is lost as a result of the modified evolution described by the CSL model with the additional hypothesis that the collapse rate λ is enhanced by the curvature of the spacetime, through the Ricci scalar R.
In this paper, we will adopt the same parameterization of λ as the authors of Refs. [87,88], i.e. we will assume that λ is a function of R (a more detailed discussion about this point will be given in Sect. III C). This idea was first suggested in [103,105], and was also one of the options mentioned in [82]. We acknowledge that there are other possibilities for including curvature scalars in the parameterization of λ. For instance, one could consider the Kretschmann scalar R αβµν R αβµν , or the scalar W αβµν W αβµν , where W αβµν denotes the Weyl tensor. In fact, as argued in [86], the latter choice seems to have the characteristics that could lead to an association of low entropy with the early state of the universe and a large entropy with its late-time state. That is, the adoption of an objective collapse model with a curvature-dependent collapse parameter involving the Weyl tensor, might result in dynamical-based explanation of the Weyl curvature hypothesis (see Ref. [106] and references within for more details). The point is that including an explicit dependence of some curvature scalar in λ, appears to be a valid option if one accepts the premise that the objective collapse is enhanced by the gravitational degrees of freedom.
The dependence on curvature could result, of course, in an effective temporal dependency of the model parameters. In fact, this kind of effective time dependence (i.e. a time dependence that appears in a certain regime to be a constant, but is in fact a coupling with some other dynamical variable) is one we have already encountered in cosmology; e.g., the masses in the standard model of particles depend on a vacuum expectation value in the Higgs sector, and that field experienced dramatic changes between the inflationary epoch and the present day. So an effective time dependence of the CSL parameters, in the cosmological context, might very well be considered likely.
With respect to the second parameter of the nonrelativistic CSL model r c , which is related to the localization of the wave function corresponding to a system of (non-relativistic massive) particles, in this article we will propose a novel approach to generalize such element when analyzing the inflationary universe. In particular, we will assume that, when taking into account the collapse of the wave function corresponding to the inflaton, the quantum two-point correlation function of the field variables is a straightforward way to extrapolate the notion of "localization" as given by the ordinary CSL model (we will be more precise in Sec. III D). Recall that, in the latter, the level of localization is encoded in the quantum uncertainty of the position operator.
Finally, let us mention that we will base our treatment on a model of spontaneous collapses that is strictly speaking not covariant. We acknowledge that a truly satisfactory proposal, to deal with the subject considered here, should be based on a fully covariant theory. However in absence of such a theory, we can move forward by making reasonable assumptions that could help in finding a complete and workable theory. Fortunately, recent proposals for special relativistic versions of these type of theories indicate that we should be able to address this shortcoming in near future [89,[107][108][109].
Our goal in this paper can be summarized as follows: based on the three main elements mentioned above, and employing standard cosmological perturbation theory at linear order, we will obtain the primordial scalar power spectrum. Also, we will analyze the characteristics needed for the spectrum to be consistent with both the laboratory constraints and those from the CMB observations.
This manuscript is divided in 6 sections plus two Appendices, where in the latter we have included most of the calculations. We start in section II with a brief review of the standard CSL model by discussing the mass proportional CSL model. Next, in section III we describe our proposed framework, in section IV we find an equivalent scalar power spectrum, and in section V we analyze the features of the obtained spectrum. Finally, in section VI, we present our conclusions. Regarding conventions and notation, we use a (−, +, +, +) signature for the spacetime metric and units where c = 1 = .

II. BRIEF REVIEW OF THE CSL MODEL
In this section, we provide a very brief review of the non-relativistic CSL model, so there is no original work here. Specifically, we will focus on the mass proportional CSL model. Those readers familiar with the subject can safely skip to the next section. For a pedagogical review, including the main technical aspects, we suggest Ref. [110], and for a general overview of collapse models see [51,52].
The CSL model is a generalization of the spontaneous collapse model originally proposed by Ghirardi, Rimini and Weber [46] to a system of identical particles. In the CSL model the collapse of the wave function occurs continuously in time, while in the GRW one the collapse happens discretely. The mass proportional CSL model has been probed in experiments, where relativistic effects can be ignored (see e.g. Ref. [111]). It consists of two main equations. The first one is a modified version of the Schrödinger equation, whose general solution is Here,T is the time-ordering operator and w(x, t) characterizes a stochastic process of white noise type. The probability for the noise w(x, t) is given by the second main equation, which represents the Probability Rule: The state vector norm evolves dynamically, i.e. it is not equal to 1. In particular, Eq. (2) implies that state vectors with largest norm are the most probable. The total probability is then P (w)dw = 1, provided that ψ, 0|ψ, 0 = 1. The parameter λ 0 is the collapse rate for a reference particle in a spatially superposed state, and usually the reference particle is chosen as a nucleon. On the other hand, the operatorÂ(x) is called the collapse generating operator, and for the mass proportional CSL model it is defined aŝ whereN i (y) =χ † i (y)χ i (y) is the particle number density operator, constructed from (non-relativistic) creation an annihilation operators for a particle of type i at y. In addition, m i is the mass of the corresponding particle species and m 0 is the mass of the reference particle, e.g. a nucleon. The second parameter of the model is the smearing length r c , a length scale characterizing the level of de-localization associated with the onset of the spontaneous collapse. In this way,Â(x) is considered as a smeared mass density operator at x.
Equation (1) describes the complete evolution of an individual system prepared in a given initial state. However, the presence of the noise w(x, t) clearly makes the evolution highly random. The only certainty is that, by construction, the system will be driven towards one of the eigenstates ofÂ(x). Therefore, it is useful to investigate what happens with an ensemble of identically prepared systems. Because of the CSL dynamics, the ensemble will be constituted of a set of different evolved state vectors, each characterized by a different realization of w(x, t), and the density matrix is the object that best serves to describe the ensemble's evolution.
The density matrix is defined from Eqs. (1) and (2) aŝ From Eq. (4), it can be shown that the corresponding density matrix evolution satisfies As a consequence, the ensemble expectation value of any observable corresponding to an operatorÔ can be obtained as Ô = Tr {Ôρ(t)}. Note that, from Eq. (1) [or equivalently (5)], one can also identify an effective collapse rate for each particle species λ i ≡ λ 0 (m i /m 0 ) 2 , i.e. the strength of the collapse scales with the mass of the particle. This leads to one of the most important features of collapse models: the amplification mechanism. For microscopic systems (for instance, a single electron) the collapse is very weak, so superpositions can occur. For macroscopic objects (for example, 10 24 atoms) the collapse is very strong, implying a rapid and efficient collapse of the wave function.
If the model is correct, accurate values of λ 0 and r c have to be determined by empirical evidence. It should be noted that as the way the model was constructed, relativistic effects have not been taken into account; also, in ordinary matter, the collapse is mostly caused by nucleons. Until experiments can provide precise values of the parameters, one can adopt some estimated values based on reasonable arguments. For example, the choice r c ≃ 10 −5 cm, originally suggested by GRW, was based on the assumption that r c should be large with respect to the atomic dimensions and to the mean spreads around the equilibrium positions of the lattice points of a crystal. Similar arguments, which involve human perception time and the number of nucleons marking the quantumclassical threshold, lead to the choice λ 0 ≃ 10 −16 s −1 [45-47, 51, 103].

III. THE PROPOSED FRAMEWORK
In this section, we present one of the possible extrapolations of the non-relativistic CSL model, described in the previous section, to more general situations in which relativistic effects cannot be ignored. In our view, the proposed model is a natural extension of the original one, since it captures the essence and main features of the mass proportional CSL model. Also, it can be applied directly to the cosmological case; and in particular, to the inflationary universe. To accomplish this, we will follow some options discussed in [82].
The background metric will be given by a spatially flat FLRW spacetime. In conformal comoving coordinates the metric components are g µν = a 2 (η)η µν , with η µν the Minkowski metric with signature (−, +, +, +). We also introduce the Hubble parameter H, and then the comoving Hubble parameter is given as H = a ′ /a = aH, with primes over functions denoting derivative with respect to conformal time η. Inflation is characterized by a shrinking comoving Hubble radius d/dη[H −1 ] < 0, with H ≃ constant. Also, if needed, we can fix the normalization of the scale factor today, so a 0 = 1.
We proceed now to characterize the metric perturbations. In particular, we choose to work in the longitudinal gauge where the line element associated to the metric is with E and Ψ being scalar fields corresponding to scalar perturbations at first order. In fact, by assuming that there is no anisotropic stress, Einstein's equations lead to Ψ = E. From now on, we will use this result and refer to Ψ as the Newtonian potential. Moreover, in this gauge, Ψ represents the curvature perturbation (i.e. the intrinsic spatial curvature on hypersurfaces on constant conformal time for a flat universe). During inflation, the matter degrees of freedom are represented by a single scalar field φ(x, η), called the inflaton. Furthermore, we separate φ in its homogeneous part φ 0 (η) plus a small inhomogeneous perturbation |δφ(x, η)| ≪ |φ 0 |. The field φ 0 will be treated in a classical way, while the perturbation δφ can be described as a QFT in a flat FLRW background metric. The condition for an inflationary phase requires that the potential energy V of the scalar field dominates over the kinetic term φ ′2 /2a 2 . Furthermore, we will focus on slow roll inflation, where the slow roll parameter is defined as ǫ ≡ 1 − H ′ /H 2 . During slow roll inflation ǫ ≪ 1, and the inflationary phase ends when ǫ ≃ 1.

A. Semiclassical Gravity
Given the fact that we have not yet at our disposal a complete and satisfactory theory of quantum gravity, we will employ the semiclassical gravity (SCG) formulation. This framework describes the way in which the gravitational and matter degrees of freedom are related to each other. We consider the SCG approach as an effective setting rather than something akin to a fundamental theory. The SCG framework is characterized by Einstein semiclassical equations where R ab and R are the Ricci tensor and scalar respectively; also we are neglecting the contribution of a cosmological constant. The expectation value of the energymomentum tensor T ab , characterizing the matter fields, acts as a source of the spacetime geometry. In this way the description of gravity in terms of the metric is always taken as classical.
The inflationary epoch is presumed to develop at energy scales smaller than the Planck mass, therefore, one expects that the SCG approach can be considered as a suitable framework. For our case of interest, the use of SCG possess two main conceptual advantages. First, the spacetime, and thus the metric, is always classical. There is no issue with the "quantum-to-classical transition" in the characterization of the spacetime. We will not need to justify the transition from "metric operators" (e.g. Ψ) to classical metric variables (such as Ψ). The fact that the spacetime remains classical is of great importance when including the CSL model. This is because the "localization" (or collapse) of the wave function is regarded as a physical process occurring in time. Consequently, it is preferred to have a setting that admits consideration of full spacetime notions [54,69]. The second advantage is that it allows to present a transparent picture of how the primordial perturbations are born from the wave function collapse: the initial state of the universe (i.e. the one resulting after a few e-folds since the beginning of inflation) is described by the completely symmetric Bunch-Davies (BD) vacuum and the equally symmetric FLRW spacetime; the symmetry being spatial homogeneity and isotropy. Then, after the CSL mechanism has ended, i.e. when the wave function has been effectively localized/collpased, the state emanating from the collapse needs not to share the symmetries of the initial state. After the collapse, the gravitational degrees of freedom are assumed to be accurately described by Eq. (7). Thus, in this post-collapsed state, in which the symmetries of the original state have been lost, the expectation value T ab yields a geometry that generically will no longer be homogeneous and isotropic.

B. Main equations of the extended CSL model
We now introduce the two main equations that characterize our version of the CSL model. The first is the evolution equation adapted to the cosmological context: whereT is the time-ordering operator, dηd 3 x |g| is the 4-volume associated to the background metric g µν , and τ → −∞ is the conformal time at the beginning of inflation. The operatorĤ corresponds to the Hamiltonian density of the system. In addition, the CSL parameter λ and the collapse generating operatorĈ appear here. We will be more specific about these objects later on. The classical scalar field W (x, η) is of white noise type, whose probability is given by the second equation which constitutes the Probability Rule: Note that we have chosen the time coordinate to be the conformal time η, which implies |g| 1/2 = a 4 . The justification of such an election is based on the widely known result indicating that, in these coordinates, the "free" Hamiltonian of a scalar field in a flat FLRW spacetime is mathematically equivalent to the one of a scalar field in Minkowski spacetime with time dependent mass. Furthermore, from Eqs. (8) and (9), and assuming that the initial states are normalized Φ, τ |Φ, τ = 1, one obtains P (W )dW = 1 as expected.
The density matrix operator, which is an useful expression describing the ensemble's evolution, is constructed from the state vectors and their associated probabilities, here given by Eqs. (8) and (9): According to Eq. (10), the corresponding density matrix evolution is (11) whereĤ is the total "free" Hamiltonian 4 , i.e.Ĥ ≡ d 3 x |g|Ĥ. As a first step, let us ignore the CSL term in the evolution equations and concentrate on finding the "free" Hamiltonian of the theory during inflation. Given that we are working within the semiclassical gravity framework, only the matter degrees of freedom will be quantized. Consequently, the field variable operator can be taken as y(x, η) ≡ a(η)δφ(x, η). Expanding the action of the standard inflationary model (i.e. a single scalar field-the inflaton-with canonical kinetic term and minimally coupled to gravity) up to second order in the perturbations, one can find the action associated to the matter perturbations. Thus, the second order action is and V ,φ indicates partial derivative with respect to φ.
The previous Lagrangian can be used to obtain the Hamiltonian density involved in Eq. (8). We define the canonical momentum p(x, η) ≡ ∂L (2) y /∂y = y ′ , in this way, the Hamiltonian H y ≡ py ′ − L (2) y is where we have used φ ′ 0 = aHM P √ 2ǫ and the slow roll approximation −V ,φ ≃ 3Hφ ′ 0 /a. At this point we promote y and p to quantum variables, with equal time commutator Also, recall that we are working within the SCG framework, where the metric variables are always classical and only the matter fields are quantized. In particular, in our approach, because of the CSL mechanism together with SCG, the metric would be changing from a homogeneous and isotropic spacetime to another one with actual inhomogeneities/anisotropies. The latter in fact would become manifested in the metric perturbations such as Ψ. As a consequence, the terms involving Ψ in Hamiltonian (13), could be treated as reflecting the effect on φ of the metric response to the "fluctuations" of the field itself δφ. In other words, the Ψ-terms can be considered as a back-reaction on the evolution of the field. However, in [78] was shown that these terms are of second order in the slow-roll parameters. Therefore, the Ψ-terms in the Hamiltonian (13) can be ignored at the leading order considered in this work. Thus, the total "free" Hamiltonian during inflation, turns out to be [see the first term in the r.h.s. of Eq. (11)] Let us turn our attention to the CSL terms in the evolution equation, Eq. (8) [or equivalently Eq. (11)]. As mentioned above, the parameter λ is the collapse rate. Notice that, in the original GRW collapse model (as well as in the mass proportional CSL model), there is an effective collapse rate that depends on the mass of the particle involved. This feature leads to the amplification mechanism, which basically states that as more massive particles are involved, their collapse rates become strengthened. On the other hand, R. Penrose and L. Diosi have, for a longtime, advocated that the collapse of the wave function might be a dynamical process, in which the underlying mechanism might be related to gravitational interaction [48][49][50]104]. Driven by these ideas, and as extensively discussed in Ref. [82], we think it is natural to consider that the collapse rate should incorporate some aspects tied to the spacetime curvature for systems where this factor is important. A possible extrapolation of λ to more general regimes is the one that assumes an explicit dependence of λ upon the spacetime curvature. That is, in laboratory experiments where one can safely presume a flat spacetime, the value of the CSL parameter is approximately λ 0 . While in a regime where the spacetime curvature is strong, such a value would completely change.
As we have mentioned in the Introduction (as well as in Ref. [82]), there are several ways to include a spacetime curvature dependence in the collapse rate parameter. For example, if one were restricted only to include the geometrical features of the spacetime, a reasonable parameterization would be λ(W 2 ), where W 2 ≡ W αβµν W αβµν with W αβµν the Weyl tensor. However, in the case of a FLRW spacetime W = 0, hence one would need to include the first order perturbations of the metric so W = 0. The corresponding analysis, although valid, would require second-order perturbation theory [see e.g. Eq. (11)] which is beyond the scope of the present work. On the other hand, if it is assumed that the relevant curvature scalar, present in the parameterization of λ, is that due to the presence of a local matter distribution, then a dependence of scalars constructed using the Riemann tensor seems to be the reasonable option (because of Einstein's field equations). Once again, there are many choices for expressing λ as a function of a curvature scalar constructed using R α βµν , e.g. R, R µν R µν , R αβµν R αβµν , etc. 5 In Refs. [87,88], motivated by finding a resolution to the black hole information loss paradox based on the CSL 5 In Appendix A of Ref. [82] we introduced a toy model in which the CSL parameter λ depends on the Kretschmann curvature scalar R αβµν R αβµν (the case for λ(W 2 ) was also considered). This model, however, did not involve a FLRW spacetime; instead, we analyzed the spacetime corresponding to a static spherical object of certain radius and a homogeneous mass-density distribution.
mechanism, the authors assumed that λ is determined by the Ricci scalar R as with λ 0 the value of the CSL parameter corresponding to the non-relativistic CSL models, α > 0 is a free parameter and µ corresponds to an appropriate physical scale with dimensions of length −2 . We point out that, in Refs. [87,88], the parameter α was constrained to α ≥ 1. In flat spacetime (or in laboratory experiments perfomed in regions of the spacetime where the curvature can be safely neglected) Eq. (16) reduces to λ 0 , i.e. the standard (non-relativistic) CSL parameter. The advantage of Eq. (16), is that it shows explicitly how the collapse rate is now sensitive to the local curvature. In other words, the local curvature scalar R serves to adjust the dominance of the non-linear, stochastic term over the linear term in Eq. (8) (or Eq. (11) equivalently). Note the former term is the one responsible for breaking the linear superposition of various basis vectors. Of course, at the present moment, one cannot verify or falsify the hypothesis that led to Eq. (16). However, in the future, it is quite possible that direct observations of gravitational effects in the quantum regime could be achieved at laboratory scales, which, in turn, will provide strong empirical evidence in favor or against of the proposal that resulted in Eq. (16). Motivated by the aforementioned discussion, we will assume that the collapse rate parameter λ, used in the CSL model considered in the present paper, is also given by Eq. (16) identically. In order to provide a reasonable value of µ, we note first that the CSL model, shown in Eq. (8), evolves the state of the quantum field with respect to the conformal time η. Therefore, we are explicitly considering a particular foliation of the spacetime, in which the states evolve from one space-like surface of constant η to another one. Additionally, in our proposal, the rate of collapse is enhanced by the curvature of the spacetime, according to Eq. (16). Hence, in principle, the collapse rate depends on the time parameter defined by the foliation λ = λ(η).
Assuming that the matter content in the universe can be described by the energy-momentum tensor of a perfect fluid T µν = (ρ+P )u µ u ν +P g µν , where ρ, P correspond to the energy density and pressure in the rest frame of the fluid respectively, and u µ represents its 4-velocity (relative to the observer). Then, from Einstein equations (in trace form) R = −T /M 2 P , together with equation of state P/ρ = ω ≃ constant, one obtains At this point, we can discuss the specific physical scale µ employed in our case of interest. During slow roll inflation, ω inf ≃ −1 and ρ inf = 3M 2 P H 2 inf ≃ constant, so R ≃ 12H 2 inf . Let us introduce a spatial physical length of size ℓ phys = a|τ |. This scale 6 induces the particular choice µ = ℓ −2 phys . In this case, we can estimate |τ | ≃ 1/(a τ H inf ), and by taking into account the number of e-folds N since the beginning of inflation a = e N a τ , we conclude that R/µ ≃ e 2N ≫ 1. Therefore, the choice µ = ℓ −2 phys ensures that, during the full inflationary regime, the collapse rate λ in Eq. (16), is enhanced by the term R/µ. In particular, one has and from now on we will use this expression for the collapse rate when analyzing the CSL model during inflation.
It is also relevant to analyze the behavior of λ in the subsequent cosmic stages after inflation. For pure radiation ω rad = 1/3, one obtains that R rad vanishes exactly. For pressure-less matter (dust) ω ≃ 0, the Ricci scalar is R matt ≃ ρ matt /M 2 P . After inflation ends, we can approximate the matter content in the universe by a twocomponent mixture of radiation and pressure-less matter. The Ricci scalar is then R = i ρ i (1−3ω i )/M 2 P = R matt , where the index i labels each component of the matter content, i = {rad, matt}. Therefore, after the inflationary regime, the evolution 7 of the quantity R/µ is given as where we used the choice µ = ℓ −2 phys . The quantity a eq denotes the scale factor evaluated at the matter-radiation equality epoch. In particular, one has a eq = e Neq a τ , where N eq is the total number of efolds from the beginning of inflation up to a eq . Moreover, assuming an inflationary energy scale of ρ 1/4 inf ≃ 10 15 GeV and that the energy scale at the radiation-matter equality 6 We note that |τ | is, to a very good degree of approximation, the size of the comoving particle horizon at any time η after the end of inflation. The particle horizon at time η can be envisaged as the intersection of the past light cone of an observer at point p with the spacelike surface at η = τ . Causal influences must originate inside this region. Only comoving particles whose worldlines intersect the past light cone of p can send a signal to an observer at p. Therefore, the physical particle horizon a|τ | is a good candidate for a physical length to be compared with R during inflation. Conversely, other natural choices do not enhance the term R/µ during inflation, e.g. assuming µ = ℓ −2 P , where ℓ P is Planck's length, implies R/µ ≪ 1 because inflation involves energy scales less than Planck scale. So, the term R/µ does not affect the collapse rate, i.e. λ ≃ λ 0 . Evidently, a similar situation occurs if one chooses µ = H 2 inf , which is essentially R during inflation, hence R/µ ≃ 1, so λ ≃ λ 0 . That is, once again the curvature term does not dominate in Eq. (16). 7 The total energy density ρ ≡ ρmatt + ρ rad evolves as epoch is ≃ 1 eV, we have ρ eq ≃ 10 −96 ρ inf . Also, if we consider that inflation lasts approximately 60 e-folds and that a eq ≃ 10 −4 a 0 , then N eq ≃ 115. Taking into account all these estimates, Eq. (19) yields The above equation implies that at the onset of the matter-radiation equality era R eq /µ ≃ 10 4 , and then it decays 8 as a −1 . Specifically, assuming a matter dominated regime, the quantity R/µ evaluated today is approximately of order 10 0 . Consequently, in the matter dominated epoch, a > a eq , from expression (16), we find that the CSL parameter λ will be constrained between The experimental constraints on λ 0 are between [111] 10 −10 s −1 ≥ λ 0 ≥ 10 −19 s −1 .
Thus, for α ≤ 2 the variation of λ within the interval (21) due to the matter dominated cosmological epoch, can fit the allowed experimental range. In particular, if we consider values close to the lowest empirical bounds, e.g. λ 0 = 10 −19 s −1 .
The previous analysis serves to illustrate that, as a result of the choice µ = ℓ −2 phys , the term R/µ becomes very dominant in Eq. (16) during inflation. However, after the end of inflation R/µ decays as a −1 . In particular, in the matter dominated epoch, the CSL rate parameter is essentially the same as the non-relativistic version of the CSL model. This is a desirable feature to be included in our model because all results and observational bounds, regarding the CSL model, obtained from contemplating this cosmic stage, remain intact [90,91,112].

D. The collapse generating operator
The next important element for the CSL cosmological model is the collapse generating operator. As we have argued in the Introduction, we are interested in a collapse operator that is obtained from covariant objects that reduce, in the simple non-relativistic regimes encountered in laboratory situations, to the "mass density" or the "energy density". Specifically, we would like to recover something that can be interpreted as the mass density operator used in the mass proportional CSL model [see Eq. (3)]. One possible covariant object, as mentioned in [82], is the scalar Γ ≡ (T µν T µν ) 1/2 . Moreover, given that our main focus is to implement the CSL model during inflation, and obtain the primordial spectrum, we will be working within the framework of cosmological perturbation theory at linear order. As a consequence, our choice for the collapse operator C will be motivated by the first order perturbations associated to the scalar Γ. In particular, the collapse operator will be derived from: where, d 3 y |h| is the 3-volume associated to the spatial metric h ij , the integration is over one of the hypersurfaces of constant time η, and δΓ = 1 2(T µνT µν ) 1/2 (T µν δT µν +T µν δT µν ) (23) withT µν denoting the background part. The function s(x, y; η) is a smearing function defined over the spatial hypersurfaces of constant η with dimensions of length −3/2 (as it occurs in the non-relativistic CSL model). The factor M −1 P in Eq. (22), fixes the correct dimensions of C and can also be considered to play a similar role as the reference mass usually introduced in non-relativistic CSL models (in which the mass of a nucleon is commonly used).
With the same spirit as the original CSL model, where the smearing function helps to localize a macroscopic body in position space, here we adopt the view that s(x, y; η), induces a "localization" of the fields associated to δΓ. In other words, at first sight s(x, y; η) should depend on the "quantum uncertainties" of the field variables rather than the uncertainties of the position operator as in the non-relativistic CSL model.
In the inflationary case, one can find the explicit dependence on the field variables for the operator δΓ. Using Eq. (23), and taking into account the canonical quantum variablesŷ,p, while neglecting the terms proportional to Ψ [see discussion after Eq. (14)], one has In the non-relativistic CSL model, the smearing function introduces the second parameter of the model, namely the smearing length r c . The parameter r c is a length scale characterizing the localization of the wave function through the collapse mechanism. For a single particle, if the width of its wave function is much less than r c , the effect of the collapse becomes negligible. On the contrary, if the width is larger than r c the collapse effects become significant. Therefore, as we have argued, the quantum uncertainties of the position operator, associated to the width of the wave function in position space, play an important role in defining the smearing function and in characterizing the parameter r c .
In our view, it is not at all obvious that there is a simple one-on-one connection between the values of r c , characterizing laboratory experiments effectively described by non-relativistic Quantum Mechanics, and the parameters of a CSL-like model describing a complete different regime such as the inflationary universe [81,82]. As a matter of fact, in this regime certain notions are lost, like the position operator of a particle; but, more fundamental ones emerge, e.g. the quantum fields.
In order to move forward, we will adopt as fundamental the relation between "quantum uncertainties" and the smearing function. Intuitively, we will think the smearing function s as strongly dependent on the "quantum uncertainties" of the system under consideration. For example, if the system can be effectively described by non-relativistic Quantum Mechanics, then "quantum uncertainties" in the position operator should serve to define the size of the smearing region, and thus set the value of r c . On the other hand, if the system we want to analyze is, e.g. the radiation emitted by Black Holes, then s should depend on the "quantum uncertainties" of the matter fields and, quite possibly, on the background spacetime curvature. However, it is a well-known fact that, strictly speaking, the uncertainty of a quantum field is not well-defined generically because vacuum expectation values like φ 2 (x) vac. diverge. In order to take the next step, we can focus on the two-point correlation function(s). In fact, until a full covariant CSL theory is formulated, we can only make progress by considering particular smearing functions with the desired characteristics.
Returning to our case of interest, the operator δΓ in Eq. (24) depends linearly on the field variablesp/a 2 and y/a. Without loss of generality, we will analyze the 2point function ofŷ/a, and the corresponding analysis of p/a 2 will proceed in an equivalent manner. Let us recall thatŷ/a =δφ, and concentrate on very small proper wavelengths (high-frequencies) λ phys ≪ H −1 . In the vacuum state, characterized by the standard Bunch-Davies vacuum [113], the 2-point function ofŷ/a is where r ≡ |x 1 − x 2 |. On the opposite regime, λ phys ≫ H −1 , the 2-point function in the BD vacuum becomes Therefore, as inflation takes place, the 2-point function decreases because the physical distance (ar) increases. But, as ar becomes larger than H −1 , the 2-point function remains basically a constant of order H 2 . Thus, we can make use of H as characterizing the size of the smearing region for the field variableŷ/a whose "uncertainty" the CSL model will attempt to decrease. In view of the previous discussion and the analysis of the 2-point function, we propose the following smearing function: Here, H acts akin to the r c parameter of the nonrelativistic CSL model. To further illustrate this point, we can write the explicit form of the collapse operator obtained from Eq. (22) by gathering all its elements. Using the smearing function in Eq. (27), the collapse operator is then Switching to Fourier space, and using Eq. (24), the collapse operator can be expressed as: We observe that, for modes such that k ≫ aH, the op-eratorĈ k (η) becomes exponentially suppressed. That is, the effect of the CSL term on the evolution of the quantum state is negligible. On the contrary, for modes such that k ≪ aH, the effect of the collapse (induced by the CSL evolution), grows stronger. This is the reason why in our proposed model H acts in a similar manner as the r c parameter. Furthermore, this is consistent with the fact that super-Hubble modes (k ≪ aH) contribute the most to the observed spectrum. Namely, we expect the collapse affects the most to the state characterizing the modes which have the strongest presence in the primordial inhomogeneities, and which later become imprinted in the CMB. For a more precise analysis regarding this subject, we invite the reader to check B. After presenting our version of the CSL model, and how it is adapted to inflation, in the next section, we will focus on obtaining one of the main predictions of the inflationary paradigm: the primordial power spectrum. However, before ending this section, let us discuss how the modified Schrödinger equation, corresponding to the mass proportional CSL model employed in ordinary laboratory circumstances, can be recovered from our proposal.
The generalized collapse operator, C(x, η), as given by Eq. (22), involves three distinctive elements: (i) the perturbation associated to the quantity Γ, (ii) the smearing function s and (iii) the reference mass M P . Regarding (i), given the definition Γ ≡ (T µν T µν ) 1/2 , it is clear that in the non-relativistic limit, Γ → ρ where ρ denotes the mass density function. 9 In particular, for several nonrelativistic massive particles ρ(x) = i m i N i (x), where 9 This case is analogous to that encountered in General Relativity when one is interested in recovering the Poisson equation for the Newtonian potential from the Einstein field equations, ∇ 2 Φ = 4πGρ . The matter density source ρ comes from the energymomentum tensor N i (x) is the particle number density for a particle of type i at x and m i is the mass of the corresponding particle species. With respect to element (ii), from the arguments presented in this section [recall discussion after Eq. (24)], in the regime where non-relativistic Quantum Mechanics can be used, the smearing function takes the form s → (πr 2 c ) −3/4 e −|x−y| 2 /2r 2 c . And as regards to (iii), the reference mass suitable at laboratory scales corresponds to the mass of a nucleon m 0 instead of the reduced Planck mass M P . Furthermore, since ordinary experiments are performed in regions of the spacetime where the local curvature can be ignored, from Eq. (16), we have that λ(R) → λ 0 . In view of the previous arguments, in the non-relativistic quantum regime, the CSL term proposed in Eq. (8), which modifies the Schrödinger equation, reduces to whereÂ(x, t) is given in Eq. (3). Note that we have deliberately made a distinction between w and W due to the fact that there is very little knowledge about the physical origin of this term, so, in principle we cannot assure that the noise term is the same for all physical systems and for all scales (however its statistical features are the same).

IV. THE NEWTONIAN POTENTIAL AND AN EQUIVALENT POWER SPECTRUM
In this section we present the derivation of an important equation, relating the primordial curvature perturbation and the proposed CSL model of the previous section. Also, we provide an expression for the primordial power spectrum, which is equivalent but not exactly equal to the standard one.
The perturbed Einstein equations (EE) at linear order in the longitudinal gauge can be combined into [114,115]: with the following definitions: By using the definition of the slow roll parameter ǫ and z = aM P √ 2ǫ, we can combine equations (31) to obtain where we have made use of the SCG approach. Also, Eq. (33) is exact, i.e no approximations were made. This equation relates the quantum expectation value of the matter degrees of freedom, obtained from the CSL collapse mechanism, and the primordial curvature perturbation represented by the Newtonian potential, which is always classical. In particular, in the vacuum state ŷ vac = 0, meaning that no curvature perturbations are present. It is only after the collapse has taken place that ŷ = 0, and the primordial perturbation is "born." We introduce a well-known quantity, defined generically as where ρ and P are associated to the type of matter driving the expansion of the universe. A main feature of the quantity R is that, for adiabatic perturbations, it is conserved for super-Hubble scales, irrespective of the cosmological epoch one is considering [114]. From the components of the energy-momentum tensor of a single scalar field, we have ρ + P = φ ′2 0 /a 2 = M 2 P H 2 2ǫ/a 2 , and because of Friedmann's equation H 2 = a 2 ρ/3M 2 P , one obtains The above equation is exact, but during slow roll inflation ǫ ≪ 1; consequently, at the lowest order in ǫ, we can approximate the latter expression as where in the last equality we have used our main equation (33). Another important aspect of the quantity R is that, in the comoving gauge, it represents the curvature perturbation. In fact, the primordial power spectrum usually shown in the literature is associated to R. The scalar power spectrum (associated to the curvature perturbation in the comoving gauge and in Fourier space) is defined as where P s (k) is the dimensionless power spectrum. The bar appearing in (37) denotes an ensemble average over possible realizations of the stochastic field R k . In the CSL inflationary model, each realization will be associated to a particular realization of the stochastic process characterizing the collapse mechanism, which in turn is related to a specific realization of the noise W . On the other hand, our main equation (33), was obtained in the longitudinal gauge. Fortunately, Eq. (34) relates Ψ, Ψ ′ and R exactly. That is, we can compute the curvature perturbation in the longitudinal gauge (where the SCG framework together with the CSL collapse mechanism generate the primordial inhomogeneities as given in Eq. (33)), and then, we can switch to the comoving gauge in order to compare the primordial spectrum obtained in our model with the standard one. In particular, we can use approximation (36) to compute the scalar power spectrum, associated to R k , that results from our main equation (33). This is, Finally, from definition (37) and Eq. (38), we can identify an equivalent scalar power spectrum as:

V. ANALYSIS OF THE SCALAR POWER SPECTRUM
The next step is to use the CSL model described in section III, to compute the expectation values in Eq. (39). Furthermore, one needs to evaluate P s (k) in the super-Hubble regime −kη → 0 to be able to compare the corresponding theoretical prediction with the observational data. The calculations are long but straightforward; all the technical details employed are presented in A.
At the lowest order in the slow roll parameter ǫ, the primordial power spectrum obtained is where N is the total number of e-folds since the beginning of inflation, and γ ε is the Euler-Mascheroni constant. We also introduce the coefficient and the function F (σ k , −kη)| −kη→0 is defined in (A53) [see also Eq. (A52)]. One of the main assumptions used to obtain the primordial spectrum (40) was to fix the free parameter α = 3/2, see Eq. (18). This is consistent with the previous analysis, in which for α ≤ 2 and the matter dominated cosmological epoch, the parameter λ can fit the experimental bounds coming from the non-relativistic CSL model. We have analyzed other values of the free parameter, for example α = 1, 2 , and we found that these values do not yield the correct shape of the spectrum, i.e. the one favored by observational data.
We proceed to examine the shape and amplitude of the predicted power spectrum (40). As shown in A, in the super-Hubble regime, the function F appearing in the expression for the power spectrum behaves as which means that F (σ k , −kη)| −kη→0 ≪ 1. The latter result implies that, in Eq. (40), the last term of P s (k) can be neglected. Additionally, we note that it is possible to express F (σ k , −kη) as Hence, the fact that F (σ k , −kη) ≪ 1, when −kη → 0, means that the width of the wave functional associated to y k , which in the BD vacuum state is a Gaussian centered at zero, is decreasing during inflation due to the CSL mechanism. In other words, the CSL process is effectively localizing (or "collapsing") the initial wave function.
On the contrary, if we turn off the collapse mechanism (i.e. λ = 0), then σ k = 0 and F (0, −kη) = 1 [see also Eq. (A52)]. Thus, in the absence of the CSL collapse term, P s (k) = 0. This is an expected result because if there is no collapse, then the quantum state of the field remains homogeneous/isotropic, and, as a consequence, there are no inhomogeneities/anisotropies in the spacetime either.
The expression for P s (k) was obtained working at the lowest order in the slow roll parameter ǫ. Henceforth, at this order, the predicted spectrum must be nearly scale invariant for it to be consistent with the observational data. At the next leading order in ǫ, one could find the standard scale dependence reflected in the so called scalar spectral index. It is thus clear that any significant departure from a perfect scale invariant spectrum in Eq. (40) would mean a fatal failure for the CSL model proposed in this work. We have argued that the last term in Eq. (40) can be neglected if λ = 0. Therefore, the other term that might induce a strong departure from scale invariance is the one with the coefficient σ k . Let us focus on this term.
The observed amplitude of the temperature anisotropies in the CMB constrains the value of H during inflation, and this value is approximately H ≃ 10 −5 M P ǫ 1/2 . Moreover, we can choose a particular value of λ 0 consistent with laboratory experiments, for instance λ 0 ≃ 10 −17 s −1 , or in Planck units λ 0 ≃ 10 −61 M P . From the definition of σ k , Eq. (41), it then follows the estimate Considering a typical energy scale for inflation V 1/4 ≃ 10 −3 M P , implies ǫ ≃ 10 −3 . Assuming also a total duration for inflation of N = 65, yields |τ | ≃ 10 5 Mpc. In addition, by taking into account k in the range of observational interest [12], one has k ∈ [10 −4 , 10 −1 ] Mpc −1 . Putting together all these values, we find that the coefficient is within the following estimated range: Consequently, σ k ≪ 1 for the range of k that contributes the most to the observed spectrum. We point out that other values of λ 0 , within the window allowed by laboratory experiments, also lead to a strong suppression of the coefficient σ k . The former analysis indicates that the primordial power spectrum (40) is scale invariant to a very good degree of approximation. Thus, we have found that, at the lowest order in the slow roll parameter, the predicted spectrum is essentially the same as the traditional one, including its amplitude, i.e.

VI. CONCLUSIONS
The adequate implementation of the CSL model during inflation has been critically analyzed in recent works [80][81][82][83][84][85]. In particular, in Ref. [82] the most important elements to be considered were discussed, according to our point of view. These are: (i) a suitable framework for treating the quantum degrees of freedom and gravity, which we think is provided by the semiclassical gravity approximation, (ii) a collapse-generating operator that can be constructed from covariant objects and, at the same time, can be interpreted as the "mass density" operator usually encountered in the laboratory applications of the CSL model, and (iii) a generalization of the CSL parameters that includes a dependence on the spacetime curvature. In this work, we have provided a concrete realization based on these three hypotheses. We have shown that the proposed model yields a primordial scalar power spectrum with the characteristics required to be consistent with observational CMB data, i.e. we have found a nearly scale invariant power spectrum. Moreover, there is no conflict between our proposal and empirical constraints, found in laboratory experiments, on the parameters of the non-relativistic version of the CSL model.
On the other hand, given that our model is not fully covariant, some potential problems might arise. However, we have taken a first step in exploring the theoretical landscape corresponding to the extrapolation of the CSL model suitable to more complex systems, in which gravity and quantum fields play an important role. Evidently, open issues remain to be analyzed, and we hope to contribute in addressing some of them in the future ahead.
Plata I+D G175 and PIP 112-2020-0100729CO of CON-ICET (Argentina). G.R.B. is supported by CONICET (Argentina) and he acknowledges support from grant: PIP 112-2017-0100220CO of CONICET (Argentina). We thank Philip Pearle for encouraging this study. We are also especially grateful to Daniel Sudarsky for his constructive comments and suggestions.
Appendix A: Computation of the primordial spectrum The standard primordial power spectrum is expressed in Fourier space. So, we start by defining the Fourier transform of any scalar field as: In Fourier space, Eq. (15) takes the form The quantization procedure will be carried out in the Schrödinger picture, so it is more convenient to work with real variables, which later can be associated to Hermitian operators. In this way, we separate the canonical variables in their real and imaginary parts. In Fourier space, this is given aŝ The quantum commutator in Eq. (14) implies where s =R,I and δ ss ′ is Kronecker's delta. The quantum "free" Hamiltonian term corresponding to (A2) is then with the following definitionŝ In the Schrödinger approach, the quantum state of the system is described by a wave functional, Φ[y(x, η)]. In Fourier space (and since the theory is still free in the sense that it does not contain terms with power higher than two in the Lagrangian), the wave functional can also be factorized into mode components as In the field representation, the operators would take the form: The decomposition shown in Eq. (A7) and the linear evolution equation (8) imply that the quantum state of each mode evolves independently (see also [85]). Thus, in Fourier space, the CSL evolution equation of the quantum state, corresponding to the wave functional(s) of each mode Φ R,I k (y R,I k ), is where the collapse generating operator (in Fourier space) is given aŝ and the Probability Rule is 2πλ/dη (A11) We will adopt the standard assumption regarding the vacuum state. That is, at an early conformal time τ , the modes are in their adiabatic ground state, which is a Gaussian centered at zero with certain spread. This ground state is commonly referred to as the Bunch-Davies (BD) vacuum. Taking into account that the initial quantum state is Gaussian, and the CSL evolution equation (A9) is quadratic in the canonical variablesŷ R,I k ,p R,I k , the form of the state vector in the field basis at any time is (A12) Therefore, the wave functional evolves according to the CSL equation (A9), with initial conditions given by Those initial conditions correspond to the BD vacuum, which is perfectly homogeneous and isotropic in the sense of a vacuum state in quantum theory.
Considering that our main goal is to obtain the primordial power spectrum, see Eq. (39), we focus our efforts in computing the quantity ŷ k ŷ q * . The expectation values, of course, will be evaluated at the evolved state provided by (A9). In terms of the real and imaginary parts we have Also, from (A14), it is clear that we are interested in computing the quantities ŷ R,I k 2 . In fact, the calculation of the real and imaginary part are exactly the same, so we will only focus on one of them. Additionally, we will simplify the notation by omitting the indexes R,I from now on.
Using the Gaussian wave function (A12), and the two main CSL equations for the mode's quantum state [Eqs.

First term of (A15)
We now turn our attention to the first term on the right hand side of (A15), i.e ŷ 2 k . The two main CSL equations for the quantum state of each mode, Eqs. (A9) and (A11), can be used to obtain the equation of evolution for the density matrix. This yields one equation per Fourier mode, which can be written as (A16) The latter equation can also be obtained form the equation of evolution of the density matrix (11) in Fourier space. Thus, it serves as a consistency check.
The importance of Eq. (A16) is that one can use it to derive the equation of evolution for the stochastic mean of the quantum expectation value of any operatorÔ k . This is, We take into account that λ during inflation is given by Eq. (18). This is, Moreover, by choosing α = 3/2 and estimating a τ ≃ 1/(H|τ |), we obtain We define the quantities Q ≡ ŷ 2 k , R ≡ p 2 k and T ≡ p kŷk +ŷ kpk , notation warning: do not confuse this R with the Ricci scalar. The evolution equations of Q, R and T obtained from (A17) are: and We combine Eqs. (A20) to obtain a single differential equation for Q, where the source function S is In Ref. [80] it is shown that Eq. (A23) can be solved by introducing the Green function of the free theory, where y k (η) is a solution of y ′′ k (η) + k 2 − w k (η) y k (η) = 0 (A26) and W = y ′ k y * k − y * ′ k y k is the Wronskian, which for Eq. (A26) is a constant. The solution to Eq. (A23) is then where τ → −∞ is the conformal time at the beginning of inflation.
Considering the approximated form of the scale factor during inflation a(η) ≃ −1/Hη, with H ≃ constant; one has w k (η) ≃ k 2 − 2/η 2 . The solution of (A26) is then Consequently, the Green function (A25) is Moreover, in the super-Hubble limit, we have On the other hand, the amplitude of the source function S(η) is exponentially suppressed in the regime kτ → −∞, while it reaches its maximum in the opposite regime −kη → 0. Therefore, when computing the integral in (A27), we can consider the approximated expression (A30); i.e. we use the form of G(η,η) in the super-Hubble limit.
Furthermore, we are interested in computing the power spectrum that is used to compare the theoretical predictions with the observational data. Specifically, we have to evaluate Q in the regime of observational interest, i.e. in the super-Hubble limit. After taking this limit, at the leading order we finally obtain where N is the total number of e-foldings from the beginning of inflation, γ ε is the Euler-Mascheroni constant and the coefficient σ k is defined in (41).

Second term of (A15)
Continuing with the second half of the calculations, we focus now on the second term on the right hand side of (A15); i.e. the term [4Re(A k )] −1 . We take the time derivative of (A9), obtaining We apply the CSL evolution operator, as characterized by Eq. (A32), to the wave function (A12) and regroup terms of order y 2 , y 1 and y 0 ; the evolution equations corresponding to these terms are thus decoupled. Fortunately, the evolution equation corresponding to y 2 only contains A k (η), which is the function we are interested in. The evolution equation is then where w k and λ 1 are defined in Eqs. (A21) and (A22) respectively. Performing the following change of variable with N ≡ 2(i + 2λ 1 a 3 e − k 2 a 2 H 2 ) (A36) and g ≡ − 1 2 a 2 H 2 (A37) equation (A34) can be recast in a more useful way: We note that if λ 1 = 0 (i.e. no collapse) f k satisfies the usual equation f ′′ k = −w k f k . The exact solution of Eq. (A38) is difficult to find, but it can be done perturbatively in λ 1 . Therefore, at leading order in λ 1 , the function P (η) is explicitly given as The perturbed solution to Eq. (A38) can be written as: The function f (0) k is the solution of the mode equation for λ 1 = 0. Specifically, we can choose the normalized mode function in the Bunch-Davies vacuum, which at leading order in slow roll is given by Inserting expansions (A40), (A42) in (A38), one finds that the function f In order to solve Eq. (A44), we will work at the leading order in slow roll, this is, a(η) ≃ −1/Hη, which implies w k = k 2 −2/η 2 . Therefore, the function Σ(η) is explicitly Σ(η) = 2ie −k 2 η 2 η 5 H 3 9 − 4k 2 η 2 + 2η 4 k 4 (A45) The solution to Eq. (A44) is where G(η,η) is the Green function introduced previously, Eq. (A29). As before, the amplitude of the source term is exponentially suppressed in the regime kτ → −∞, so the mayor contribution to the integral in Eq. (A46) comes from the super-Hubble limit. This means that we can use the approximated expressions for f (0) k (η) and G(η,η), Eq.(A30) in such a limit. The result is where we have used the definition x ≡ −kη. The function Ei(z) is the exponential integral function defined as and Erf(z) is the Gaussian error function Returning to the original variable A k , Eq. (A35), with the perturbed solution in terms of Eqs. (A43) and (A47), we obtain From the observational point of view, we are interested in taking the limit x → 0. Retaining the dominant terms in the perturbed solutions, we find that the real part of A k can be approximated by In fact, we rewrite the previous expression using the coefficient σ k defined in Eq. (41), which yields Re(A k )| x→0 ≃ kx 2 2 1 + 16 σ k x 6 (A52) We also introduce the function As we can see from Eq. (A52), in the limit x → 0, the function F (σ k , x)| x→0 goes to zero as x 6 , i.e.
|F (σ k , x)| x→0 ≪ 1 On the other hand, if there is no collapse, this is, if λ = 0, then σ k = 0 and F (0, x)| x→0 = 1 It is now straightforward to write the explicit expression for the equivalent power spectrum. Substituting Eqs. (A31) and (A53) in (A15), and making use of Eq. (39), we obtain the primordial power spectrum as shown in Eq. (40). In this Appendix, we provide an analysis regarding the evolution of the subhorizon modes as dictated by the CSL evolution equation. For ease of the calculations, we will rely on the evolution equation for the density matrix Eq. (A16), the analysis of course would be equivalent if the CSL master equation (A9) was used instead. Therefore, from Eq. (A16), we observe that the term affecting the dynamical evolution of any mode, due to the CSL process, isD We proceed to analyze the last term in the previous equation, i.e. the one with the highest power of a. The analysis corresponding to the other terms is equivalent. Thus, the dynamical evolution for the last term is dictated by the function ∆ k (a) ≡ a 5 e −k 2 /a 2 H 2 . (B4) We note that ∆ k (a) < 1 if and only if κ 2 > 5 ln a, where κ ≡ k/aH. Assuming that inflation lasts approximately 65 e-folds, then a ∼ e 65 . As a consequence, ∆ k (a) < 1 if and only if κ > O(10). In particular, for a subhorizon mode κ ≫ 1, one has that ∆ k (a) ≪ 1. This demonstrates that, for subhorizon modes, the additional dynamical evolution introduced by the CSL mechanism is exponentially suppressed.
On the other hand, an analogous analysis for superhorizon modes κ ≪ 1, indicates that ∆ k (a) ≫ 1, i.e. the superhorizon modes are the most affected by the CSL evolution term.