Eternal inflation and the quantum birth of cosmic structure

We consider the eternal inflation scenario of the slow-roll/chaotic type with the additional element of an objective collapse of the wave function. The incorporation of this new agent to the traditional inflationary setting might represent a possible solution to the quantum measurement problem during inflation, a subject that has not reached a consensus among the community. Specifically, it could provide an explanation for the generation of the primordial anisotropies and inhomogeneities, starting from a perfectly symmetric background and invoking symmetric dynamics. We adopt the continuous spontaneous localization model, in the context of inflation, as the dynamical reduction mechanism that generates the primordial inhomogeneities. Furthermore, when enforcing the objective reduction mechanism, the condition for eternal inflation can be bypassed. In particular, the collapse mechanism incites the wave function, corresponding to the inflaton, to localize itself around the zero mode of the field. Then the zero mode will evolve essentially unperturbed, driving inflation to an end in any region of the Universe where inflation occurred. Also, our approach achieves a primordial spectrum with an amplitude and shape consistent with the one that best fits the observational data.


I. INTRODUCTION
The theory of the early Universe enjoying the most recognition amongst cosmologists is inflation [1][2][3][4]. The inflationary epoch is characterized by an accelerating expansion of the Universe, and its current success is based on the power to explain the primordial inhomogeneities generation that represent the seeds of cosmic structure [5][6][7][8][9]. Moreover, the recent Planck satellite data release reveals that the description of the primeval Universe is consistent with the inflationary paradigm [10][11][12]. These data suggest that the primordial perturbations spectrum is very close to scale invariance, favoring the simplest inflationary models [12,13].
In spite of the strong agreement between the inflationary predictions and the observational data, there are still some unsettled conceptual issues; among these is the usually called quantum-to-classical transition of the primordial perturbations. 1 As initially presented in [16], a precise formulation of that issue is as follows: Starting from an initial situation which is taken to be perfectly isotropic and homogeneous, both in the background spacetime and the quantum state characterizing the matter fields, and based on a dynamics that supposedly preserves those symmetries, one ends with a non-symmetric situation corresponding to the inhomogeneities and anisotropies observed in the Universe.
Most of the accepted proposals to address the quantum-to-classical transition of the perturbations during inflation are based on the decoherence framework [17][18][19][20][21][22][23][24], with a varying degree of reliance on non-orthodox interpretations of quantum mechanics, e.g. many-worlds [25,26], and the consistent histories formulation [27] (although we do not subscribe to such proposals because of the arguments exposed in [28][29][30]). In fact, in recent works [31,32], it has been shown how decoherence in de Sitter spacetime can lead to a violation of the time translation invariance (nonetheless, the spatial and vacuum symmetries remain intact). In any regard it is fair to say that the quantum measurement problem, and in particular in the cosmological context, is still very much an unsolved issue. Therefore, it is reasonable to explore all consistent proposed solutions. One possible solution to the aforementioned shortcoming is to consider an objective dynamical reduction mechanism 2 [16,28,29,36]. Such a process can break the homogeneity and isotropy of the inflaton vacuum state and, in turn, generates the primordial curvature perturbation.
In the present work, we employ the continuous spontaneous localization (CSL) model as the objective reduction mechanism. The CSL model belongs to a large class of models simply called collapse models, which attempt to provide a solution to the measurement problem of Quantum Mechanics [37][38][39][40][41]. The common idea in these models is to modify the Schrödinger equation by introducing nonlinear stochastic corrections that spoil its linearity, inducing a spontaneous random localization of the wave function in a sufficiently small region of space. The model parameters are chosen in a way that ensures that micro-systems evolve following closely the dynamics pro-vided by the Schrödinger equation; conversely, macrosystems are extremely sensible to the nonlinear effects resulting in a sharply localized wave function. Nevertheless, the main aspect of the CSL model is that there is no need to mention or to introduce the notion of an observer/measurement device as in the orthodox interpretation. This is a desired element in the context of the early Universe and cosmology in general. The implementation of the CSL model into the inflationary scenario has been analyzed in previous works [42][43][44]. Some of them resulting in theoretical predictions consistent with the observational data, in particular, with the observed amplitude and shape of the primordial spectrum (scalar and tensor) [42,[45][46][47][48].
Our view to implement the CSL model into inflation relies on the semiclassical gravity approximation, described by Einstein equations G ab = 8πG T ab . In this approximation, the matter degrees of freedom are treated quantum mechanically, while the gravitational sector, characterized by the spacetime metric is always classical. The validity of the semiclassical gravity approach is limited. In particular, it requires that the scalar curvature of the spacetime to be small compared to l −2 P (l P is the Planck length). It also requires that the fluctuations in the energy-momentum tensor are small compared to its absolute value [49]. The inflationary epoch is assumed to be associated with energy scales that are way below the Planck regime. Regarding the second constraint, we expect that once the collapse mechanism has ended, the quantum fluctuations of the energy-momentum tensor should be small, although we cannot be completely sure until a full covariant reduction mechanism is developed (see Ref. [50][51][52][53] for recent advances in that direction).
In this paper, we are concerned with the subject of eternal inflation [54][55][56][57][58][59][60]. In particular, we will focus only on the eternal inflation mechanism provided by potentials of the slow-roll/chaotic type [61]. At this point, it is worthwhile to differentiate the eternal chaotic inflation from the one produced by false vacuum tunneling [58,61]; the analysis of the latter type of inflationary models will be left for future work. Therefore, unless we state it explicitly, when referring to eternal inflation throughout this work, we have in mind the eternal chaotic inflation kind.
The traditional accepted idea, regarding eternal inflation, is to consider the vacuum quantum uncertainty of the inflaton as a true dynamical event, or commonly called quantum fluctuations. In the textbook presentation of inflation, the classical trajectory of the field is governed by the dynamics of the zero mode. Inflation unfolds as the zero mode slowly rolls down the potential, ending when the field lies at the bottom of the potential. However, the zero mode can fluctuate due to the quantum fluctuations of the inflaton. 3 The fluctuations can 3 Although, strictly speaking, we point out that there are no fluc-be sufficiently large that the field spontaneously goes upward the potential, starting a new cycle of inflation. This behavior can arise in any region of the Universe, thus, achieving a self-reproducing state of inflation, ordinarily called eternal inflation. It is known that eternal inflation is a typical feature of inflationary models with potentials that best fit the observational data [62].
Perhaps one of the most striking consequences of eternal inflation is its connection with a multiverse [58,63,64], in which "anything that can happen will happen, and it will happen an infinite number of times" [61]. For instance, the constants of nature can take a wide range of values depending on the region of the multiverse [65,66]. It is important to remark that the connection between eternal inflation and the multiverse is more rigorous in the case of eternal inflation via tunneling between false vacua than the eternal inflation mechanism considered in the present work. Some cosmologists view the multiverse as an issue questioning the validity of the inflationary theory [67][68][69], while others consider the multiverse hypothesis admissible and positive [64,65,70,71]. Nevertheless, as the multiverse is a widely accepted consequence of eternal inflation, we think it is worthwhile to meticulously analyze all the underlying assumptions of the possible eternal inflation scenarios (see Refs. [72,73] for other related works).
One of the basic premises of eternal inflation is that it is appropriate to treat the vacuum quantum uncertainty of the inflaton as an actual dynamical object, which literally disturbs the homogeneous part of the field, i.e. the zero mode. We find this picture misleading. For example, consider the one-dimensional quantum harmonic oscillator in its ground state. Ehrenfest's theorem guarantees that the expectation value of the momentum and position operators will follow the same dynamics as its classical counterparts. Nonetheless, there is nothing in ordinary quantum mechanics suggesting the need to modify the classical motion as a result of the quantum uncertainties. Also, the postulates of quantum mechanics indicate that, after a measurement of, say, the position, the amplitude of the wave function will change into a wave packet narrowly centered around the measured value. The peak of the wave packet, i.e. the expectation value, will continue to evolve according to the oscillator classical equations of motion. Moreover, before the measurement, the quantum state of the harmonic oscillator possessed a reflection symmetry. The dynamics of the vacuum quantum uncertainty does not lead to a breakdown of said symmetry. Only after the measurement process, the reflection symmetry of the wave function is lost.
In the case of the inflatonφ(x, t), the vacuum expectation value is precisely equal to the zero mode φ (x, t) = φ 0 (t) (the expectation value of the inhomogeneous part δ φ(x, t) is zero). Therefore, the evolution of tuations, only quantum uncertainties.
the field expectation value follows the zero mode equation of motion. Although, as in the case of the harmonic oscillator, there is nothing in the postulates of the quantum theory that signals a modification of the classical trajectory of the field; in particular, not an alteration that includes the quantum uncertainty. Furthermore, even if somehow one incorporates the quantum uncertainty into the inflaton dynamics, it is not clear that δ φ(x, t) 2 1/2 , taken in the vacuum state, should be associated to a classical stochastic field, which alters the evolution of the zero mode. At most it can be associated to the amplitude of a stochastic field, but the random nature of its phase will remain unknown. Another issue is to identify the element that plays the role of the measurement device, in the inflationary Universe, as in an ordinary laboratory situation. Thus, eternal inflation is deeply linked with the quantum-to-classical transition and the quantum measurement problem in general.
On the other hand, incorporating the CSL model into inflation could help to clarify the aforementioned issues (although open questions will remain since at the present there is no formal version of the CSL model for quantum fields). In particular, it can successfully change the symmetries of the vacuum state and, at the same time, be responsible for the birth of the primordial curvature perturbation [42,74]. In addition, since the evolution provided by the CSL model transforms spontaneously the vacuum state, it is possible to formally identify the expectation value δ φ(x, t) with a classical stochastic field δφ(x, t). The localization of the quantum state of the field, due to the CSL dynamics, implies that the quantum uncertainty of the field is decreasing during inflation, and becomes centered around a value that coincides with φ(x, t) = φ 0 (η)+δφ(x, t), which in turn avoids the condition for eternal inflation. Moreover, given that δφ ≪ φ 0 , in the usual interpretation, the zero mode does not have very large fluctuations.
The paper is structured as follows: In Sec. II, we present a brief review of slow-roll inflation focusing on the main aspects leading to the eternal inflation picture. Then, in Sec. III, we focus on the important features of the CSL model concerning its implementation into the inflationary context. Afterwards, in Sec. IV, we show how the CSL inflationary model modifies the usual conclusions drawn from slow-roll inflation resulting in the eternal inflation scenario. Finally, Sec. V contains discussion and conclusions.

II. THE STANDARD PICTURE
We will begin by reviewing the standard picture of eternal inflation within the framework of slow-roll inflation [64,66].
The inflationary Universe is described by Einstein equations G ab = 8πGT ab (c = 1). In slow-roll inflation, the matter sector is characterized by a single scalar field, the inflaton φ. Meanwhile, the background spacetime is described by a FRW Universe with line element (in conformal coordinates) ds 2 = a 2 (η)[−dη 2 + δ ij dx i dx j ]. The scale factor can be approximated by a(η) ≃ −1/[Hη(1 − ǫ 1 )], with H the Hubble factor, which during inflation is approximately constant and η the conformal time η ∈ (−∞, 0). The deviation from a perfect de Sitter expansion is parameterized by the first slow parameter defined as ǫ 1 ≡ 1 − H ′ /H 2 , where H ≡ a ′ /a = aH is the conformal expansion rate and the prime denotes derivative with respect to η. The energy density of the Universe is dominated by the inflaton potential V (φ), and during slow-roll inflation ǫ 1 ≃ (M 2 P /2)(∂ φ V /V ) 2 ≪ 1, with M 2 P ≡ (8πG) −1 the reduced Planck mass. The action of the inflaton with mass m minimally coupled to gravity leads to a quantum field theory for φ in a quasi-de Sitter spacetime. In particular, one can decompose the quantum field in Fourier modeŝ In order to avoid infrared divergences, we have introduced a regularization and consider the field in a box of side L. The sum is over the wave vectors k satisfying k i L = 2πn i for i = 1, 2, 3 with n i integers, and the field operator operator given bŷ The mode functions φ k satisfy the evolution equation and the normalization condition It is known that the normalization condition (4) does not determine the set of mode solutions unequivocally; a particular choice of the solutions of Eq. (3) implies an election of the vacuum. For k = 0, the general solution to Eq. (3) is a linear combination of η (3+2ǫ1−2ν)/2 and η (3+2ǫ1+2ν)/2 , with ν = 3/2 + ǫ 1 − m 2 /(3H 2 ). For k = 0, the general solution is a linear combination of the Hankel functions of first and second kind of order ν. We now choose the Bunch-Davies vacuum, i.e. we choose the modes such that at earlier times η → −∞ they behave as "positive frequency solutions", normalized according to Eq. (4). In fact, we highlight that the Bunch-Davies vacuum in an expanding space is more than simply a choice: for any given initial condition with certain quite natural assumptions for the UV behavior, the Bunch-Davies vacuum is always approached as it is a late-time attractor solution, see Refs; [75,76]. Therefore the Bunch-Davies vacuum corresponds to select the functions H (1) ν . The operatorsâ k andâ † k are the annihilation and creation operators respectively, which satisfy the commutation rules [â k ,â † k ′ ] = δ k,k ′ . As usual the vacuum is defined byâ k |0 = 0 for all k.
The inflaton can be split into an homogeneous part denoted byφ 0 plus small inhomogeneities, i.e.
where ε ≪ 1 is a parameter (do not confuse with the first slow-roll parameter ǫ 1 ) that helps to quantify the smallness of the inhomogeneities. The homogeneous part φ 0 (η) corresponds to the zero mode of the field. Therefore, we can rewrite Eq. (1) as followŝ whereδφ k denotes the modes of the fieldφ k with k = 0. The initial state of the fieldφ(x, η) is a featureless vacuum state, except for the zero mode that is excited. Explicitly the initial state is with k = 0. For k = 0, the state is the Bunch Davies vacuum. For the zero mode we assume a coherent state, i.e. a state such thatâ 0 |ξ 0 = ξ 0 |ξ 0 . The assumption that |ξ 0 corresponds to a coherent state is justified by noting that such states are sharply peaked around a value that coincides with the homogeneous part of the classical field φ 0 (η). In other words, the state |ξ 0 satisfies Additionally, expectation values satisfy their counterpart classical equations of motion. In particular, the homogeneous classical field φ 0 (η) satisfies which is the well known equation of motion corresponding to the (classical) homogeneous part of the inflaton. We now recall the definition of the second slow-roll parameter, ǫ 2 ≡ ǫ ′ 1 /Hǫ 1 . In terms of the slow-roll parameters, Eq. (9) is Assuming the slow-roll approximation ǫ 1 ≪ 1 and ǫ 2 ≪ 1, we recover the usual slow-roll equation for the homogeneous field 3Hφ ′ 0 ≃ −a 2 m 2 φ 0 . On the other hand, the quantum uncertainty 4 of the fieldφ(x, η), defined as operatorX is defined as δ QX ≡ X2 − X 2 ; however, we will abuse the language and refer to both δ 2 QX and δ QX as the quantum uncertainty.
is not zero. In particular, it is given by If we assume a renormalized energy-momentum tensor, then we can consider normal ordering. Therefore, for the coherent state |ξ 0 , one has ξ 0 |φ 0 (η) 2 |ξ 0 = (φ 0 (η)2Re[ξ 0 ]) 2 , i.e. the quantum uncertainty of the zero mode vanishes. Thus, the quantum uncertainty of the fieldφ is In summary, the quantum fieldφ(x, η) is in the initial state |in . The expectation value φ (x, η) is proportional to the homogeneous (classical) field φ 0 . The quantum uncertainty of the fieldφ(x, η) is given by 0|δφ(x, η) 2 |0 , that is, by the quantum uncertainty of the inhomogeneous part characterizing the inflaton.
In the traditional picture, eternal inflation is achieved when the variation of the quantum uncertainty of the field becomes large enough, so that the value of the field φ cannot be properly localized. In other words, the standard approach considers the quantum uncertainty as an actual dynamical fluctuation, referred to as quantum fluctuations. These quantum fluctuations will be superimposed onto the deterministic slow-roll trajectory obeyed by φ 0 . The fluctuations will lead the value of φ away from φ 0 in some regions, originating a new cycle of inflation in those regions. This process will occur endlessly, resulting in an infinite number of regions inflating, sometimes called "pocket universes" [61]. In some regions inflation may end, while others continue inflating.
Formally, eternal inflation arises when the variation of φ 0 over one Hubble time is smaller than the variation of 0|δφ(x, η) 2 |0 over the same time period. The mathematical expression can be deduced as follows.
The variation of φ 0 over one Hubble time is simply where we have switched from conformal to cosmic time t (the dot over the function represents derivative with respect to t).
Next, one needs an estimate for the amplitude associated with the variation of the quantum uncertainty 0|δφ(x, η) 2 |0 over one Hubble time. That estimate is usually obtained by introducing an UV cut-off [77] k H (η), which is the comoving wave number k H corresponding to Hubble radius crossing at a time η, so k H (η) = −1/η. Physically, the cut-off implies one is considering only super-Hubble modes since these modes contribute most to the quantum uncertainty amplitude. Therefore, the variation of the quantum fluctuations, in terms of the Fourier modes with k = 0, is Approximating the mode functions At this point, we can provide the condition for eternal inflation, using Eqs. (14) and (16), the explicit condition is [61,70] Equivalently, in terms of the inflaton potential and the first slow-roll parameter, the former condition is It is usually claimed that essentially any slow-roll potential satisfies the condition for eternal inflation [78]. For instance, let us consider the potential used in R 2 inflation, also known as Starobinsky's inflation [79], which is in very good agreement with the latest observational bounds from the Planck mission [12,13,80]. The first slow-roll parameter for this potential is where x ≡ 2/3φ/M P . Therefore, the condition for eternal inflation using this type of potential is If the inflaton is well enough in the slow-roll region x ≫ 1, which means the field takes super-Planckian values, then the condition for eternal inflation, Eq. (21), can be satisfied because e 2x ≫ 1. Furthermore, even if the energy scale is sub-Planckian m 4 < M 4 P , the field generically fulfills the condition for eternal inflation. In fact, eternal inflation commonly occurs for potentials of the hilltop type [78].
In the following, we will show how the traditional claims regarding eternal inflation can be modified if one takes into account the dynamical reduction of the wave function.

III. THE CSL MODEL DURING INFLATION
As we mentioned in the Introduction, we will use the CSL model as the particular model characterizing the dynamical reduction of the wave function. The implementation of the CSL model into the picture of slow-roll inflation has been analyzed in previous works [42,[46][47][48] (in Refs. [43][44][45] the CSL model is also considered during inflation but based on a different conceptual approach). Here, we present only the main features and results that will be of interest for the present work.
Before addressing the CSL inflationary model, we present our view regarding the relation between the spacetime description in terms of the metric and the degrees of freedom of the inflaton (see [16,28] for a detailed presentation). This particular view is based on the semiclassical gravity framework, which treats gravitation classically and all other fields quantum mechanically. We assume such a framework to be a valid approximation during the inflationary era, which is well after the full quantum gravity regime has ended. Note that this is a major difference between our approach and the standard one, since in the latter all scalar degrees of freedom (metric and matter fields) are quantized. On the other hand, it is not entirely settled that one should quantize the metric degrees of freedom. There are numerous arguments suggesting that the spacetime geometry might emerge from deeper, non-geometrical and fundamentally quantum mechanical degrees of freedom [81][82][83][84]. Therefore, the framework we will use is based on Einstein's semiclassical equation Let us discuss now how the dynamical reduction of the wave function fits into our general understanding of gravity and its, still incomplete, quantum mechanical characterization (see [85,86] for a more detailed explanation). According to general relativity, gravitation reflects the structure of spacetime itself, whereas quantum theory seems to fit most easily in contexts where this structure is a given one. That is, quantum states are associated with objects that "live" in spacetimes. For instance, the standard Schrödinger equation specifies the evolution of a system; the quantum states of fields characterize the system in connection to algebras of observables associated with predetermined spacetime regions. One of many complications that become manifest in the merging of gravity and quantum mechanics, following the canonical approach, is the appearance of a timeless theory, where the recovery of fully covariant spacetime notions is nontrivial [87]. Therefore, it seems natural to speculate that it is precisely in this setting where something that departs from quantum orthodoxy, as the dynamical collapse of the wave function, might find its origin. That is, the collapse process could reflect some remanent signature from a fundamental quantum gravity regime. It is worthwhile to mention, that Penrose [88] and Diósi, [89][90][91] have suggested similar ideas regarding the connection of the collapse of the wave function and its possible relation to gravity. According to this view, it seems reasonable to conjecture that a reduction of the wave function could correspond to lingering features of a quantum gravity theory. If that is the case, the emergence of spacetime itself would be tied to the incorporation of such effective quantum description of matter fields living on the spacetime, and evolving (approximately) according to standard quantum field theory on curved spaces, with some small deviations which might include our hypothetical collapse.
A specific example of this particular outlook is to analyze the departures from Eq. (22) that arise when considering a single and instantaneous collapse as described by the state |ψ(t) = Θ(t 0 − t)|0 + Θ(t − t 0 )|ξ (where Θ(x) is the Heaviside function). Using that state, Eq. (22) can be written as where the object ξ ab = Θ(t − t 0 )( ξ|T ab |ξ − 0|T ab |0 ) would represent an individual stochastic step. If we extrapolating this example, then the stochastic gravity framework [92,93] might correspond to the continuous version of several discrete stochastic steps, e.g. the one given by the CSL dynamical collapse model. Furthermore, we accept that at the fundamental level, the spacetime would have a quantum description in terms of some unspecified variables (e.g. those of loop quantum gravity, the causal set or dynamical triangulations approaches), but during the inflationary regime the spacetime can be viewed in effectively classical terms. A useful analogy is the hydrodynamical characterization of a fluid, which at certain scales is a useful description. However, at subatomic scales one should have a radically different depiction. Einstein equations would correspond to the Navier-Stokes equations, the spacetime metric to the density and velocity fields; whereas, the subatomic characterization of the fluid constituents would correspond to the fundamental degrees of freedom of quantum gravity. Moreover, this analogy also serves to picture a situation in which departures from the semiclassical gravity approximation occur. Specifically, when a phase transition (from say liquid to gas) is taking place, one might expect not only the violation of Navier-Stokes equations but also of the hydrodynamic characterization. In our case, small violations of Eq. (22) might occur when the dynamical collapse of the wave function is going on, much like in the fluid analogy during a phase transition. From now on, we will assume the validity of the semiclassical Einstein equations by assuming that the dynamical process, associated with the collapse of the wave function, has reached a final (different from the initial vacuum) state. However, we should bear in mind all the open issues and subtleties that we have indicated.
Our starting point will be the same as the standard inflationary account. We assume the initial state of the Universe characterized by the homogeneous and isotropic classical FRW spacetime and the equally symmetric Bunch-Davies vacuum. The dynamical reduction mechanism will drive the initial state of the matter field to a final state that does not need to share the symmetries of the Bunch-Davies vacuum; this mechanism acts as an effective spontaneous and stochastic quantum collapse of original the wave function. As a consequence, T ab will not have the symmetries of the initial state leading to a geometry, through Einstein semiclassical equation, that generically will no longer be homogeneous and isotropic. We would like to mention here that it might be valuable to consider the decoherence framework, plus coarse graining, since as found in Refs. [31,32], the time translation invariance is lost when considering said framework. Therefore, it would be interesting to investigate whether coarse graining could potentially lead to a breakdown of the symmetries of the quantum state.
Focusing on the scalar metric perturbations in the longitudinal gauge, and assuming no anisotropic stress, the semiclassical equations, in Fourier space, yield where Φ is the scalar metric perturbation (or the Newtonian potential in the longitudinal gauge) andδφ characterizes the inhomogeneities of the inflaton. From the former equation, we observe that in the vacuum state, the expectation value 0|δ φ ′ k |0 vanishes, which implies that Φ k is also zero for all k. Strictly speaking, there are no perturbations and the spacetime is perfectly symmetric. It is only after the CSL mechanism has evolved the vacuum state into the post-collapse state that generically δ φ ′ k = 0 and, by Eq. (24), the curvature perturbations are born. In this way, the self-induced collapse leads to a configuration in which neither the state of the field nor the spacetime is homogeneous and isotropic.
Next, we focus on the collapse mechanism. The CSL model is based on a modification of the Schrödinger equation. This alteration induces a collapse of the wave function towards one of the possible eigenstates of an operator called the collapse operator, with certain rate λ. The selfinduced collapse is due to the interaction of the system with a background noise W(t) that is characterized as a continuous-time stochastic process of the Wiener kind. The parameter λ sets the strength of the collapse process; the greater the value, the stronger the collapse (see [41] for a throughly review).
Since the CSL model modifies the Schrödinger equation, it is convenient to describe the theory of the inflaton in the Schrödinger picture, where the relevant objects of the theory are the wave function and the Hamiltonian.
The initial wave functional, corresponding to the quantum state given in Eq. (7), is where Ψ ξ0 denotes the wave function of the homogeneous part of the field, which is a coherent state. Meanwhile, Ψ BD denotes the wave function associated to the inhomogeneous sector of the inflaton, which corresponds to the Bunch-Davies vacuum. It will be suitable to separate the field into its real and imaginary parts; thus, in Fourier space, the initial wave function is with k = 0. On the other hand, the Hamiltonian characterizing the inhomogeneous sector of the field is H = 1 where we have defined the rescaled field variables (28) Additionally, we have ignored terms containing firstorder slow-roll parameters since we will be mostly interested in the amplitude of the fields rather than the shape of their spectrum. Note also that p R,I k = aδφ ′ R,I k ; thus, the right hand side of Eq. (24) involves the expectation value of the momentum operator.
Promoting y R,I k and p R,I k to quantum operators, the commutations relations are [ŷ R,I k ,p R,I k ] = iδ(k − k ′ ). The CSL mechanism will drive the initial state vector, represented by the wave function Ψ in [φ], Eq. (26), towards a final state that lacks the original symmetries of the system. However, the zero mode part of the wave function (recall that Ψ ξ0 corresponds to a coherent state) will remain unchanged, this is needed for the model to be self-consistent [85]. On the other hand, the collapse mechanism will change the part of the wave function corresponding to the modes k = 0. In our framework, these modes, together with the collapse of their wave function, are responsible for the birth of the primordial inhomogeneities.
We now apply the CSL model to each (non-zero) mode of the field independently. The initial state vector |Ψ R,I k , τ , with τ the conformal time at the beginning of inflation, will evolve according to T denotes the time-ordering operator. Note that we have chosen the momentum operatorp R,I k as the collapse operator. This is justified because the metric perturbation Φ is directly related to the expectation value of the momentum operator, Eq. (24). However, we could also have chosen the field variable as the collapse operator and perform a similar analysis with our conclusions unchanged; see Ref. [42]. The collapse parameter λ k depends on each mode. From dimensional analysis, one finds that λ k is dimensionless. Additionally, λ k sets the strength of the collapse, and since in principle all non-zero modes are subjected to it, we expect that λ k ≫ 1 for all k.
In addition, we want to highlight some open issues regarding the CSL inflationary model. First, the CSL model is a non-relativistic model; hence, in order to use it in the cosmological setting, it is required to introduce suitable modifications; see [94] for a detailed discussion. However, because of the perturbative nature of our analysis, such modifications will not be of importance in the present work. A second related issue is that, as is known [95], the CSL model in general violates energy conservation, although the energy increase is minimal (for example a particle of mass m = 10 −23 g, would take 10 10 years for an energy increase of 10 −8 eV [41]). Nevertheless, even if the energy violation can be neglected, at the fundamental level a more realistic model should remove this issue. In the case of the CSL inflationary model, the energy violation issue could lead to divergences in the energy-momentum tensor (see a recent work in which these violations in the energy momentum tensor are explicitly computed [53] and shown to be small). Third, let us recall that the amplification mechanism of the CSL model is the idea that the collapses must be rare for microscopic systems, but the effect of the collapses must increase when several particles are hold together forming a macroscopic system. This amplification mechanism is included in the CSL inflationary model in the k dependence of the λ k parameter. Nonetheless, the explicit k dependence, as well as the choice of the collapse operator, is purely phenomenological at this point. In our opinion, this third open issue is the most pressing one at the moment for the collapse models in the inflationary setting. We will say more about this at the end of the present section.
Returning to the CSL model during inflation, characterized by Eq. (29), we note that the initial wave function of each non-zero mode of the field is represented by a Gaussian centered at zero with certain spread. We expect that the CSL mechanism will not change the Gaussian nature of the evolved state (clearly, the spread and the peak of the Gaussian will not be the same as the original). Also, since we have selectedp R,I k as the collapse operator, it will be suitable to work in the momentum representation; thus, the wave function is characterized by The dynamical CSL evolution of the wave function follows Eq. (29), provided by the initial conditions A k (τ ) = (2k) −1 , B R,I k (τ ) = 0, C R,I k (τ ) = 0, corresponding to the Bunch-Davies vacuum. In fact, for the purpose of the present paper, we will be only interested in A k (η). The evolution equation for such a quantity is

This equation can be solved by performing a change of variable
, which results in a Bessel differential equation for f . After solving such an equation, and returning to the original variable A k , we obtain where q 2 ≡ k 2 (1 − 2iλ k ), and J n represents a Bessel function of order n.
One of the reasons we are interested in the quantity A k (η) is because it is related to the prediction for the scalar power spectrum within the CSL inflationary model. Recall that the dimensionless power spectrum in Fourier space is defined as where E{·} denotes an ensemble average over possible realizations of the stochastic field Φ. In our approach, each realization of Φ corresponds to a particular postcollapse state; the stochasticity of the post-collapse state is generated from the noise function W. Equation (24) implies that as a matter of fact, E{ p R k 2 } = E{ p I k 2 }. In can be shown that [42] E{ p R,I The predicted scalar power spectrum, within the CSL framework, is explicitly calculated in [42]. The result is which is scale invariant if λ k = λ/k (recall that λ is the universal CSL parameter). At this point we want to discuss some important aspects that led to the latter result. First, we can observe that the explicit dependence on k of the λ k parameter is purely phenomenological, and at the time of writing we cannot provide a more fundamental justification. Additionally, in order to obtain that result one had to fix a particular gauge (the longitudinal gauge) even before the quantization process. This gauge fixing is inevitable, since the use of the semiclassical gravity framework implies that matter and geometry degrees of freedom should be treated differently (the former are quantum mechanical while the latter are always classical). Also, the choice of gauge implies that the time coordinate is attached to some specific slicing of the perturbed spacetime, and thus our identification of the corresponding hypersurfaces (those of constant time) as the ones associated with the occurrence of collapses-something deemed as an actual physical change-turns what is normally a simple choice of gauge into a choice of the distinguished hypersurfaces, tied to the putative physical process behind the collapse. This aspect of the collapse process was also found in Ref. [96]. Naturally this is an issue that leads to tensions with the expected general covariance of a fundamental theory, a problem that afflicts all known collapse models, which in the non-gravitational settings becomes the issue of compatibility with Lorentz invariance of the proposals. We must acknowledge that this generic problem of collapse models is indeed an open issue for the present approach. One would expect that its resolution would be tied to the uncovering of the actual physics behind what we treat here as the collapse of the wave function (which we view as a merely an effective description). As has been argued in related works, and in following ideas originally exposed by Penrose [88], we hold that the physics that lies behind all this, ties the quantum treatment of gravitation with the foundational issues afflicting quantum theory in general, and in particular those with connection to the quantum measurement problem.
Finally, it will be convenient to write Eq. (36) reexpressing the potential as V (φ) = m 4 V 0 (φ/µ) (that is V 0 is dimensionless). The predicted spectrum is then where V * 0 and ǫ * 1 indicate that such functions are being evaluated at the time η * when the amplitude of the curvature perturbation is equal to the temperature of the anisotropies observed in the CMB, (δT /T ) 2 CMB ≃ 10 −9 .

IV. MODIFICATION OF THE STANDARD PICTURE
Having established the main features of the CSL inflationary model, we proceed to show how this model modifies the standard results regarding eternal inflation.
We consider first the expectation value of the inflaton in the state that resulted from the CSL mechanism, i.e. Ψ CSL |φ(x, η)|Ψ CSL = Ψ CSL |φ 0 (η)|Ψ CSL + ε Ψ CSL |δφ(x, η)|Ψ CSL . (38) As we mentioned in the previous section, the CSL mechanism only affects the non-zero modes of the field; therefore, the initial state of the zero mode, which is a coherent state, remains unchanged. On the other hand, the expectation value ofδφ can be characterized as Furthermore, since the CSL mechanism is based on a stochastic process in which each realization of the noise W α , with α denoting the particular realization, corresponds to a different state |Ψ CSL , the field δφ(x, η) is a classical stochastic field. The expectation value given in Eq. (38) is then Switching back to cosmic time, the variation ∆ Ψ CSL |φ(x, t)|Ψ CSL over one Hubble time ∆t = H −1 , using Eq. (40), is Moreover, εδ φ must be smaller thanφ 0 in order to maintain the consistency of the perturbative analysis. Consequently, the variation of the expectation value of the field over one Hubble time is The previous analysis contrasts with the corresponding one from the standard picture. As we saw in section II, the expectation value of the fieldφ(x, η) is taken with respect to the initial state, Eq. (7), during the whole inflationary regime; henceforth, in|φ(x, η)|in = φ 0 (η). In other words, in the standard approach, if one wishes to characterize classically the field φ(x, η) from the quantum fieldsφ 0 ,δφ and using only the |in state, then there is an issue because the expectation value ofφ can only provide homogeneous fields. In order to overcome such an issue, the traditional posture postulates the association δφ(x, η) ∼ δ φ(x, η) 2 1/2 , which, in principle, it is not clear where the stochasticity of the field δφ(x, η) comes from. In fact, after renormalization [57], the quantum expectation value δ φ(x, η) 2 leads to a result that is independent of x, and thus it is homogeneous.
Conversely, in our approach, the classical nature of the fields φ 0 and δφ, as well as their stochastic origin, are clearly identified in Eqs. (38) and (40).
We now focus on the quantum uncertainty of the field φ(x, η) using the CSL evolved state, i.e.
In discrete Fourier space, the latter expression is Note that only when k = −k ′ , one obtains a nonvanishing result. Using the wave function given in Eq. (30), and taking the continuum limit L → ∞, we find We have omitted the ε 2 factor as the only role played by ε is to ensure that the scalar perturbations are small compared to the background. Equation (45) implies that the uncertainty δ 2 Qφ (x, η) is essentially equal to the sum (integral) of the uncertainty of each mode of the field, that is, ]. An estimated value for the variation of δ 2 Qφ over one Hubble time can be found following the same procedure as in the standard approach. Namely, we introduce an UV cut-off k H (η) that selects only the super-Hubble modes. Using Eqs. (32) and (45), the variation ∆δ 2 where we have defined ζ k e iθ k ≡ √ 1 − 2iλ k . Furthermore, since k H (η) = −1/η, and dη/dt = a −1 = −Hη, we obtain ∆δ 2 Qφ (x, t) = ∆t with and ∆N the number of e-folds from the beginning of inflation to some time η. From Eq. (47), we observe that if there is no collapse, i.e λ = 0, which implies that the state is the Bunch-Davies vacuum, then we recover the estimate of the standard treatment; see Eq. (16).
Equation (47), also serves to illustrate the effect of the collapse. As long as λ H ≫ 1, the uncertainty of the field δ 2 Qφ will start to decrease with respect to the uncertainty in the vacuum state, which is expected from a reduction of the wave function. In fact, let us consider a super-Hubble mode k + , that is, k + |η| < 1. If λ H ≫ 1 then λ/k + ≫ 1; henceforth, the collapse affects the modes whose physical wavelength is greater than the Hubble radius during inflation, which are also the ones that originate the seeds of the observed structure in the Universe. Contrarily, for a sub-Hubble mode k − , i.e. k − |η| > 1, if λ H ≪ 1 then λ/k − ≪ 1; this implies that the uncertainty depicted in Eq. (47) is essentially the same as the one obtained in the vacuum state. Thus, in the CSL reduction model, the collapse affects the super-Hubble modes, during inflation, more than the sub-Hubble modes.
The main information contained in Eq. (47) is that the CSL evolution strongly localizes the wave function of the field if λ H ≫ 1. The localization is achieved around the expectation value of the field, i.e. φ 0 + δφ. Hence, the wave function of the inflaton is effectively centered around φ 0 , which obeys its classical equation of motion. Given that δφ < φ 0 , the evolution of φ 0 continues essentially undisturbed.
Next, we consider whether the condition for eternal inflation (17) is fulfilled within the CSL inflationary framework. The condition for eternal inflation occurs when the variation of the width of the wave function amplitude is bigger than the variation of its peak over a Hubble time.
The ratio of the variation of the quantum uncertainty with respect to the variation of the expectation value of the field, over a Hubble time, can be obtained from Eqs. (42) and (47). This ratio is Equation (49), is the main result of this section. Consequently, if λ H ≫ 1 then, in terms of the potential and the first slow-roll parameter, expression (49) is As we will show next, using a particular potential, the assumption λ H ≫ 1 implies that which is the opposite of the condition for eternal inflation to occur. On the other hand, the prediction for the scalar power spectrum is directly proportional to λ as illustrated in Eq. (37). If λ H ≫ 1 then we must adjust m such that the overall amplitude of the power spectrum fits the observational data, i.e. it must be consistent with P(k) ≃ 10 −9 . In particular, since λ|τ | is a large number, the factor m 4 /M 4 P must be quite small in order to satisfy Moreover, given the smallness of m 4 /M 4 P , Eq. (50) is ensured to be quite small (i.e. ≪ 1).
To further illustrate our picture, let us consider once again the R 2 inflationary model; the potential is given in Eq. (19). In this model, the number of e-folds from the beginning of inflation to some time η, is [80] recall that x ≡ 2/3φ/M P , so x τ denotes the value of the field at the conformal time τ corresponding to the beginning of inflation. From Eq. (50), we obtain But given Eq. (53) it is clear that ∆N ≫ x, hence Therefore, the assumption λ H ≫ 1, and considering sub-Planckian energy scales, implies that the condition for eternal inflation is bypassed, Eq. (51). In order to fit the observed amplitude of the spectrum, the mass of the inflaton must be where x * refers to the value of the field at the time η * .

V. DISCUSSION AND CONCLUSIONS
In this paper, we have considered an additional element to the standard picture of eternal inflation of the slow-roll/chaotic type. This element is the dynamical reduction of the wave function. The initial wave function of the fieldφ is characterized by a Gaussian, centered around φ 0 , with certain spread δ Qφ . The value of φ 0 is associated with the expectation value ofφ.
The traditional condition for eternal inflation is On the other hand, when we apply the self-induced collapse of the wave function, provided by the CSL mechanism, to the inflationary Universe we find that it is possible to overturn the condition for eternal inflation. In particular, we obtain as long as λ H ≡ λ|τ |e −∆N ≫ 1 holds. The parameter λ is the universal CSL collapse parameter, τ is the conformal time at the beginning of inflation, and ∆N the number of e-folds from τ to some time η during inflation. Equation (57) indicates that, during one Hubble time, the variation of the quantum uncertainty of the field becomes larger than the variation of the homogeneous field φ 0 . In other words, in the standard setting, one claims that the quantum uncertainties become actual dynamical fluctuations that dominate over the classical dynamics of φ 0 . Therefore, in the standard scenario, as the Universe inflates, the spread of the wave packet, corresponding to the wave function amplitude, increases more rapidly than the variation of its peak, failing to localize the center of the wave packet. And since the dynamics of the background is driven by φ 0 , one cannot assure that inflation will end because the peak of the wave packet, which follows the classical trajectory, has been diluted inside its width. Conversely, Eq. (58) shows that if the collapse is strong enough, i.e. λ H ≫ 1, then the width of the wave function will decrease, localizing itself effectively around φ 0 , which will drive inflation to an end.
In the traditional picture, the failure to localize φ 0 along the wave function is sometimes used to claim that the field is as likely to roll up as it is to roll down the potential, leading to an endless cycle of inflation for different regions of the Universe since the value of the field, in distinct regions, varies according to the spread of the Gaussian characterizing the initial wave function of the field. In our model, the CSL reduction mechanism strongly localizes the wave function around the expectation value of the field φ ≃ φ 0 in any sector of the potential. Once properly localized, the field φ 0 will evolve essentially unperturbed according to its classical equation of motion; this occurs in any region of the Universe indistinctly.
Another important difference between the standard approach and the CSL inflationary model involves the prediction for the power spectrum. In the traditional setting, the theoretical prediction for the power spectrum is of the same structure as the left-hand side of Eq. (57). On the other hand, in our approach, the predicted expression for the power spectrum is ∼ H 2 λ|τ |/φ, which is actually a different expression from the ratio [∆δ 2 Qφ (x, t)] 1/2 /∆φ 0 that results from the CSL model; see Eq. (58). The fact that two latter expressions are different, together with the additional parameter λ, provides our model with more freedom to achieve the observed amplitude of the spectrum and to evade the condition for eternal inflation. Concretely, the CSL inflationary model contains two main parameters: the collapse parameter λ and the overall energy scale of inflation m 4 . As seen from Eq. (58), the condition for eternal inflation can be avoided if λ|τ |e −∆N ≫ 1; additionally, the energy scale of inflation m 4 has to be adjusted to fit the observed amplitude of the spectrum. In the standard inflationary picture, the parameter m 4 is also fixed by the amplitude of the spectrum. However, once this parameter is fixed there is no other parameter that helps to bypass condition for eternal inflation resulting in a generic feature of most inflationary models.
An additional aspect worth considering is the implication when λ H ≪ 1. This situation virtually entails that no collapse has occurred for a particular set of modes, namely, the sub-Hubble modes. Moreover, the uncer-tainty of the fieldφ(x, t) is essentially a sum (integral) of the uncertainties of all mode operatorsφ k (t). In particular, δ 2 Qφ (x, t) is a sum of two kinds of quantum uncertainties. One kind is associated to the sub-Hubble modes, which we know are roughly equal to the uncertainty of the vacuum state. The other type encloses the uncertainties corresponding to the super-Hubble modes, which are decreasing due to the spontaneous collapse process. Henceforth, in the CSL inflationary model, the localization ofφ(x, t) is not completely perfect because the uncertainty of the sub-Hubble modes is bigger than the uncertainty of the super-Hubble modes. We think this is a trait inherited by considering the collapse acting on each mode of the field. We expect that a covariant reduction mechanism of the wave function, that acts on the field variablesφ(x, t) and/or its conjugated momentum field, will alter this picture leading to a sharp localization ofφ(x, t), or, equivalently, to a collapse that affects all the field modes equally. However, even if the localization ofφ provided by the CSL model is not perfect, the field expectation value consists of the zero mode plus δ φ , which is always smaller than φ 0 . Thus the zero mode does not suffer from large disturbances.
In summary, we have considered whether or not the standard arguments, regarding eternal inflation, are modified by considering the self-induced collapse of the wave function. Although there might be some issues related to the initial conditions for inflation [67]; once inflation has started in a region of the spacetime, the dynamical reduction of the wave function will generate the primordial inhomogeneities corresponding to all the observed structure in the Universe. After the primordial perturbations are born, the wave function of the field will be narrowly centered around φ 0 , disrupting the condition for eternal inflation. Note that the analysis shown in the present work, which considered the self-induced collapse proposal, only focused on the eternal inflation process of the slow-roll/chaotic type. Henceforth, our analysis would modify the usual conclusions regarding a never-ending inflationary phase and all its ramifications concerning the multiverse, only to those kind of inflationary models. On the other hand, we left for future work the analysis concerning the eternal inflation mechanism based on false vacuum decay [58,61], which many consider the essential concept leading to the multiverse landscape.