Emergent universe revisited through the CSL theory

In this work we analyze how the spectrum of primordial scalar perturbations is modified, within the emergent universe scenario, when a particular version of the Continuous Spontaneous Localization (CSL) model is incorporated as the generating mechanism of initial perturbations, providing also an explanation to the quantum-to-classical transition of such perturbations. On the other hand, a phase of super-inflation, prior to slow-roll inflation, is a characteristic feature of the emergent universe hypothesis. In recent works, it was shown that the super-inflation phase could generically induce a suppression of the temperature anisotropies of the CMB at large angular scales. We study here under what conditions the CSL maintains or modifies these characteristics of the emergent universe and their compatibility with the CMB observations.

There are renowned merits attributed to the inflationary paradigm. In addition to solving the horizon problem, predictions from the simplest slow-roll inflationary model, such as a spatially flat geometry characterizing the universe and a quantum origin of the spectrum of primordial perturbations (i.e. a nearly scale-invariant power law), are some of them. Moreover, the predictions are extremely consistent with recent observations from the cosmic microwave background (CMB) radiation [16]. However, the exploration of some alternatives seems to be interesting, in the light of small features in the CMB power spectra that remain unexplained and also open discussions in some recent works, which we mention below.
Some studies have drawn attention to the lack of largeangle correlations and a weak power in the low-ℓ multipole moments of the angular power spectrum, in the observed CMB temperature anisotropies with respect to that predicted within the standard ΛCDM model. This was first mentioned in the COBE results [17] and later confirmed in subsequent generations of satellites [18][19][20][21][22][23][24][25]. This feature on the largest angular scales was analyzed, with some controversies included, by several authors (see for instance [26][27][28][29][30][31]). As is well known, the largest observable angular scales contain direct information on primordial physical processes that occurred during the inflationary era (or prior to it), and could only have un- * gabriel@iafe.uba.ar † mpp@fcaglp.unlp.edu.ar ‡ gleon@fcaglp.unlp.edu.ar dergone modifications by the physics involved in the relatively recent past, e.g. through the late-time integrated Sachs-Wolfe effect. In fact, several authors have shown that some of the observed CMB anomalies could be explained in this way, e.g. [32,33]. However, there is still no verdict on whether it is just a statistical fluke or if there really is something new and interesting behind it. Today it constitutes one of the persistent large-angle anomalies in the CMB data that makes up the list of current challenges of the standard ΛCDM model [34,35]. In search of a convincing explanation, new theoretical ideas have been considered from different approaches. The results for the suppression of the low multipoles in the CMB spectrum can give us clues towards new physics that determined the initial conditions for slow-roll inflation or even tell us something about an earlier phase before slow-roll inflation .
As mentioned above, the prediction that our current universe is spatially flat is one of those that is often mentioned in relation to inflationary models, and the data seem to indicate that indeed the observed universe is very close to flat. However, it does not imply that the spatial sections are exactly flat. Cosmological models in the context of non-flat cases have also been explored with interesting results, see e.g. [59][60][61][62][63][64][65][66][67][68][69][70][71][72][73][74][75]. Recently, in light of the analysis of some observational data, the debate about what is the spatial curvature of the universe has resurfaced, typically quantified via its fractional contribution to the cosmic energy budget today parameterized by Ω K [76][77][78][79][80][81][82][83][84][85][86][87][88][89][90][91][92][93]. This parameter has an important role in determining the evolution of the universe and is closely related with the early universe physics, because if the spatial curvature is positive, then the curvature term will always dominate at early enough times in a slow-rolling inflationary epoch. And indeed, some of these recent analyzes suggest that a small positive spatial curvature could be present favoring, perhaps, the case of a closed universe. It is also interesting to mention that the question of spa-tial curvature participates in the so-called "H 0 tension" between sets of early time probes and a number of late time observational data (see e.g. [94][95][96][97]). Additionally, spatial curvature also affects the search for the nature of dark energy, an energy component to explain the current stage of accelerated expansion and that would constitute around 70% of the energy budget of the universe according to the ΛCDM model [87,[98][99][100][101]. In particular, it is known that by assuming flatness, when in fact Ω K = 0, would induce critically large errors in reconstructing the dark energy equation of state, showing that including curvature as a free parameter is imperative in any future analyses that attempts to determine whether dark energy is a cosmological constant or is something more exotic. And, on the other hand, closed universe models can generally relax the Hubble tension between supernovae observations and the CMB. However, at this stage, there is no conclusive evidence for a positive curvature, but this is at least an attractive possibility supported by the data that warrants further exploration.
An important feature that was highlighted in [53], and of particular interest for the present work, is that a phase of super-inflation (i.e. a period where the Hubble parameter increases with time) prior to slow-roll inflation could be related to the suppression of power in the low CMB multipoles. Furthermore, any mechanism that attempts to solve the cosmological singularity problem, within a semiclassical spacetime description, will naturally contain such a phase.
The emergent universe (EU) of [127,130] is one in which a spatially closed universe (based on General Relativity and dominated by a scalar field minimally coupled to gravity) emerges from an initially past-eternal Einstein static state (with a finite initial size), enters a phase of super-inflation and then evolves towards slowroll inflation to finally give rise to the standard hot-Big Bang. Because of how it is built, there is neither a horizon problem nor singularity. This model has been studied in recent years through different approaches and variants, establishing its stability conditions, analyzing fine-tuning issues, its viability both theoretically and observationally and making it clear which questions are still open; see, for instance . Among the alternative scenarios of the early universe [170,171], a recent realization of an emergent universe is the so-called string gas cosmology [172][173][174][175]. In this model the universe begins in a long hot and almost static phase, dominated by a thermal gas of closed fundamental strings. Recently, some conjectures such as Swampland [176] and Trans-Planckian Censorship (TCC) [177] have put very strong constraints on possible inflationary models [178][179][180][181]. However, alternative cosmologies such as bouncing and emergent models are consistent with these conjectures and trivially satisfied. A good understanding of the emerging phase is still missing, but there are some promising approaches [182].
The next point to consider is the fact that any model that claims to provide a mechanism for the generation of the seeds of cosmic structure, must be able to give a convincing answer to the following issue. In the early stages at the beginning of the universe (and well after the Planck era has ended), the spacetime is assumed to be spatially isotropic and homogeneous. In addition, a standard assumption is that the perturbations of matter fields (e.g. the inflaton field) were in a quantum vacuum state also perfectly symmetric (the symmetry being spatial isotropy and homogeneity), usually described by the so-called Bunch-Davies vacuum. Then, an important puzzle arises, namely to explain the transition from a perfect symmetric state portraying the early universe to the non-symmetric state that characterizes the current universe, which cannot be attributed to quantum unitary evolution. The standard approach to address that issue is based on the study of the role of quantum fluctuations during, for instance, the inflationary epoch. However, since the evolution of any quantum state, according to standard quantum theory, is always dictated by the Schrödinger equation (which does not break any initial symmetry of the system or destroy quantum superpositions), the traditional early universe paradigm is incomplete in that sense. In other words, the solely existence of vacuum fluctuations is in no way sufficient to claim that there are actual inhomogeneities of any kind present in the universe. The aforementioned puzzle is usually referred in the literature as the quantum-to-classical transition of the primordial perturbations. In fact, this subject is closely related to what is known as the measurement problem in quantum physics [183][184][185][186][187][188][189][190], and is notoriously exposed in the case of the quantum description of the primordial universe [191][192][193][194]. This is because measuring devices and observers who decide when and how to perform some kind of measurements cannot be fundamental notions in a theory which seeks to describe the early universe where neither existed [195]. Cosmologists have tried to account for the quantum-to-classical transition by different types of arguments. A recent critique of different efforts at explaining this subject can be found in [196] 1 .
In order to classify the possible alternatives to address the issue described above, in [185] the measurement problem was written in a elegant manner showing that there are three statements that are mutually inconsistent. In short: (A) the physical description provided by the quantum state is complete, (B) quantum states always evolve according to the Schrödinger equation, and (C) measurements always have definite results.
If statement (A) is denied, then the quantum state does not contain all the information necessary for the description of a quantum system. In this way, the addition of hidden variables and the equations that determine their evolution is required. The best known proposal for this case is the de Broglie-Bohm model [198]. Some applications to the cosmological case can be found in [199][200][201][202][203].
Works based on decoherence [204][205][206][207] led to a partial understanding of the issue. Nevertheless, this argument by itself (i.e. without any extra assumptions) cannot address the fact that a single (classical) outcome emerges from the quantum theory. In other words, decoherence alone cannot solve the quantum measurement problem [208][209][210]. Other cosmologists seem to adopt the Everett "many-worlds" interpretation of quantum mechanics [211] plus the decoherence process when confronted with the quantum-to-classical transition in the inflationary universe, e.g. [212]. Regarding this point, we would like to refer the reader to [193,213,214] where arguments against decoherence and the Everett interpretation are also presented. This proposal and the Everettian interpretations are some of the approaches that somehow discard statement (C).
The last remaining choice is to negate statement (B). This path leads to non-standard quantum theories, i.e. theories where the collapse of the wave function is selfinduced by some novel mechanism. Known as objective collapse theories, from the mid-1970s several authors began to develop modifications to the Schrödinger equation, with the aim to alter the evolution of the wave function. In this way, the collapse of the wave function would occur without any reference to external observers or devices present that should perform some sort of measurements [215][216][217][218][219][220]. One virtue of collapse models is that they have shown in recent years to have the attractive feature of connecting plausible resolutions of other open problems in a single unified picture [221]. Reviews on these sort of theories can be found, for instance, in [222,223].
In the present work, we will approach the emergent universe from the perspective of the Continuous Spontaneous Localization (CSL) model [215,217], which will be incorporated into the situation at hand as a mechanism to break the original symmetries of the quantum vacuum state of the field driving the expansion of the early universe, and generating the primordial cosmological perturbations. In this manner, the CSL model naturally provides an explanation of the quantum-to-classical transition of such perturbations. The incorporation of objective collapse schemes and theories in the cosmological context has been studied since 2006 [192]. This has led to numerous investigations in recent years with varied proposals, particularly in the framework of semiclassical gravity (but also with exploratory works in the frame-work of standard quantization), with very encouraging results [193,194,. Other authors have investigated similar ideas and some of these works can be seen, for instance, in [245][246][247][248][249]. The debate about the particular details involving the implementation of the CSL theory into the cosmological context is still open. In fact, there is an extensive landscape of possibilities, which constitutes an active line of research at the moment [250][251][252][253][254][255][256].
In particular, a relevant aspect to examine is the following: when one decides to combine quantum field theory (QFT) with gravitation, one must choose the setting within which such a link is to be made. At the time of writing there is no complete and finished program merging successfully both theories, so a couple of options arise. In Ref. [252] some of these approaches were analyzed, evaluating their pros and cons. There, it was argued that the semiclassical gravity framework appears favored from a theoretical and conceptual point of view when one wants to incorporate collapse models. Therefore, our present analysis will be based on the semiclassical gravity (SCG) framework, in which gravity is treated classically and the matter fields are treated quantum mechanically [192,225,227,235,236,257]. This approach accepts that gravity is quantum mechanical at the fundamental level, but considers that the characterization of gravity in terms of the metric is only meaningful when the spacetime can be considered classical. Namely, we will be dealing with the description of an epoch well after the full quantum gravity regime has ended (i.e. from when the emergent universe begins to evolve), where the energy scales involved allow one to suppose valid the consideration of the metric as classical and well described by semiclassical gravity equations. Therefore, semiclassical gravity can be treated as an effective description of quantum matter fields inhabiting a classical spacetime. While this approach has received some criticisms [258,259], those arguments have been refuted [260][261][262][263][264][265][266][267]. In the particular case of CSL theory, its first implementation into the primordial universe, based on the SCG framework, was done in [227], and some of us have continued to explore its consequences. For example, in [238,240,241] it was shown that a strong suppression of primordial B-modes in the CMB is predicted generically. In addition, within that same framework, observational constraints were analyzed in [237], and in [231] it was found that the condition for eternal inflation can be bypassed. In the next section, some additional motivations for this choice in the present case of analysis will become apparent.
To finish this Introduction, let us mention that in some previous works [229,233,234], where a spontaneous collapse of the wave function was implemented during a phase of slow-roll inflation, certain features on the low multipoles of the CMB were analyzed as consequence of the collapse. On the other hand, recently in [163] the emergent universe model, originally put forward in [127,130], was analyzed and it was shown that the superinflation phase (a characteristic shared by all emergent universe models) could be responsible for part of the anomaly in the low multipoles of the CMB; in particular, for the observed lack of power at large angular scales. Motivated by these results, and under the same assumptions of [163], here we calculate the primordial power spectrum of scalar perturbations, but incorporating a particular version of the CSL model for the situation at hand. We will analyze whether the super-inflation phase in the framework of the emergent universe plus CSL continues (or not) to produce the power suppression in the low multipoles. Lastly, we will study under what conditions the CSL maintains or modifies such characteristics and their compatibility with the CMB observations.
Our manuscript is divided as follows. We start in section II presenting the theoretical framework of the emergent universe plus the CSL proposal; we also obtain the predicted scalar power spectrum. Next, in section III we present and discuss our results; there, we also perform a further exploration of the corresponding parameter space. Finally, in section IV, we present our conclusions. Regarding conventions and notation, we use a (−, +, +, +) signature for the spacetime metric and units where c = 1 = .

II. EMERGENT UNIVERSE IN THE CSL FRAMEWORK
In this section, we present the implementation of the CSL model into the emergent universe (EU) model.

A. Theoretical background
As we mentioned in the Introduction, many models of the emergent universe have been studied in recent years. Since we will be closely following the results of Ref. [163], here we will do our analysis under the same assumptions considered there, which in turn are based on [127]. In particular, we assume the action of General Relativity with a scalar field φ, which represents the dominant matter driving the universe early expansion, minimally coupled to gravity, with canonical kinetic term. The scalar potential employed (inspired by R 2 -inflation) takes the form V (φ) = (4πG) −1 (e Cφ − 1) 2 as the one reconstructed in [130] (using techniques developed in [63]). In obtaining that potential, the evolution of the scalar factor a(t) ≃ a 0 + A e H0t was taken into account; where, a 0 > 0 is the (initial) radius of the Einstein static universe, C, A are positive constants, and H 0 is the Hubble parameter at the onset of slow-roll inflation 2 .
In the EU model, after leaving its initial static state, the universe enters a slow-rolling regime (at a few efoldings after leaving its initial static state) where the 2 The type of potential V is as the one shown in Figs. 1 of Refs. [130,163].
scale factor grows sufficiently quickly to mitigate neglecting the curvature effects. This period of a de Sitter type of inflation comes naturally to an end (as the scalar field starts oscillating around the minimum of the potential), it is then followed by a re-heating phase, and finally continues to the standard hot Big Bang expansion. In Ref. [166], it was shown that the temporal evolution, given by Friedmann equation along with the scalar field Klein-Gordon equation, leads the system towards an attractor where H tends to a constant andφ 2 → 0; this is, the system evolves from an Einstein static state to a de Sitter type of expansion.
The dynamics of this model is such that, prior to the traditional slow-roll inflation, there is a phase of super-inflation where the Hubble parameter increases with time, i.e.Ḣ > 0. The mechanism which generates this superinflationary period depends on the particular model under consideration, but it is a generic characteristic of the EU scenario. For example, in the models of Refs. [127,130], it is considered a FLRW closed universe where the spatial curvature is responsible for the superinflationary period. However, let us note that we could have chosen another model from those mentioned in the Introduction, provided that the evolution of the background (given by the mentioned scale factor a(t)) produces the phases of super-inflation and slow-roll inflation, in which the generation of curvature perturbations is analyzed here.
As discussed in depth in [130], even though the traditional emergent universe is with positive spatial curvature, it is quickly negligible in a few e-foldings and furthermore slow-roll inflation can always be made to end for some negative value of φ. On the other hand, it is also possible to find the (finite) N number of e-foldings for the total slow-roll phase within the emergent universe. Analogously to [163], we will now make a first approach to the problem at hand and therefore we will neglect the contributions of the space curvature to the primordial perturbation 3 .
We follow the standard procedure and separate the scalar field and the metric into a homogeneous background plus small perturbations, i.e. g µν = g (0) µν + δg µν and φ = φ 0 +δφ. We will now fix the gauge of the perturbations, and work in the so called longitudinal (or Newtonian) gauge. In this gauge, at first order in the scalar metric perturbations, and assuming no anisotropic stress components, the corresponding line element is In these coordinates, the scale factor can be modeled by We define H ≡ a ′ /a; the prime over variables denotes derivative with respect to conformal time η.
The metric degrees of freedom will remain classical because of the semiclassical gravity approach, these include the background and the perturbation Ψ. In the matter fields sector, the background scalar field φ 0 will be treated also in a classical fashion; however, the perturbed part δφ will be subjected to quantization. Taking into account that the CSL theory modifies the Schrödinger equation, it will be convenient to carry out the quantization in the Schrödinger picture. Therefore, we focus on finding the total Hamiltonian of the system.
We introduce the (re-scaled) field variable y = aδφ. Expanding the action of the system (i.e. a single scalar field minimally coupled to gravity) up to second order in the perturbations, one can find the action associated to y. In this way, the second order action is S (2) = d 4 xL (2) , where We define the canonical momentum p(x, η) ≡ ∂L (2) y /∂y = y ′ , in this way, the Hamiltonian density is given by We now promote the fields y and p to quantum operators satisfying the following equal time commutator relation [ŷ(x, η),p(y, η)] = iδ(x − y).
Our next step is to decompose the field and the conjugated momentum in Fourier modes. This is justified by the fact that we work with a linear theory and, hence, all the modes evolve independently. In Fourier space, the total Hamiltonian corresponding to Eq. (4) takes the formĤ Furthermore, it will be convenient to work with real variables. In this way, we separate the canonical variables into their real and imaginary parts, i.e.
The quantum commutator in Eq. (5), implies where s =R,I and δ ss ′ is Kronecker's delta. Using this separation the Hamiltonian becomesĤ = R 3+ d 3 k(Ĥ R k + H I k ), with the following definitionŝ In order to apply the CSL model into the EU scenario, we will follow the approach first introduced in [227] for the inflationary regime. There, it was found that with an appropriate selection of the field collapse operators and using the corresponding CSL evolution law, it is possible to attain a "collapse" in the relevant operators corresponding to the Fourier components of the field. Furthermore, we will assume linearity in the collapse generating operator, therefore, the reduction mechanism will act on each mode of the field independently, i.e. there will be no mode mixing because of the CSL process.
In view of the above, the evolution of the state vector characterizing each mode of the quantum field as given by the CSL theory is: whereT is the time-ordering operator, and τ denotes the conformal time at the beginning of the EU regime. Note that the stochastic field W k = W R k + iW I k depends on k and the conformal time. In other words, it is reasonable to introduce a stochastic function for each independent degree of freedom given that we are applying the CSL collapse dynamics to each mode of the field. Consequently, the stochastic field W k (η) might be regarded as a Fourier transform on a stochastic spacetime field W(x, η). The probability for the stochastic field is given by the second main CSL equation: the Probability Rule, i.e.
As can be seen from the CSL evolution equation (10), we have chosen the field variableŷ R,I k as the collapse generating operator. Technically, this means that the CSL process will drive the initial state vector towards an eigenstate ofŷ R,I k . The motivation for this choice is based on the fact that when one implements the SCG formulation into cosmological perturbation theory, one obtains, at linear order in Fourier space, the following relation: where µ ≡ φ ′ 2 0 /(2H 2 M 2 P ). We observe in the above expression that the quantum expectation value ŷ k acts as a source for the curvature perturbation. This might be interpreted as indicating that the collapse is tied with some aspect of the quantum matter that "gravitates" (i.e. that would characterize the interaction between gravitation and matter degrees of freedom). Moreover, this view is, in principle, consistent with the proposals by R. Penrose and L. Diosi suggesting that gravity might play a fundamental role in the so called collapse of the wave function [218,220].
We denote by Φ[y, η] the wave functional characterizing the quantum state of the field. In Fourier space, the wave functional can be factorized into mode components . It is known that the ground state of the Hamiltonian (9), characterized by a wave functional Φ R,I 0 [y R,I , η], is a Gaussian. Also, the Hamiltonian (9) and the CSL evolution equation (10) are quadratic in bothŷ R,I k andp R,I k ; consequently, the wave functional at any time can be written in the form: The initial state of the field |Φ R,I k , τ will be the chosen as the standard Bunch-Davies (BD) vacuum. The corresponding wave functional for the BD vacuum is characterized by the initial conditions

B. Power spectrum
After having introduced the theoretical basis of our model, here we focus on deriving a prediction for the primordial spectrum, which is the observational quantity of interest. The standard expression of the primordial spectrum, associated to the curvature perturbation, is normally expressed in the so-called comoving gauge, while our main equations were obtained in the longitudinal gauge. For a single scalar field, the relation between the curvature perturbation in the comoving gauge R and in the longitudinal gauge Ψ is given by [9,212] Let us note that, in cosmic time t coordinates, µ can be expressed as µ =φ 2 0 /(2M 2 P H 2 ). Given that the system has an attractor point such that H → constant andφ 2 0 tends to an infinitesimal small number, then it follows that µ → 0. Therefore, in Fourier space, Eq. (15) implies where in the last equality we have used our main equation (12). The scalar power spectrum associated to R k is defined as where P s (k) is the dimensionless power spectrum. The bar appearing in (17) denotes an ensemble average over possible realizations of the stochastic field R k . In the CSL model, each realization will be associated to a particular realization of the stochastic process characterizing the collapse. We can use approximation (16) to compute the spectrum associated to R k , i.e.
From definition (17) and Eq. (18), we can identify an equivalent scalar power spectrum as: In terms of the Real and Imaginary parts ofŷ k , the ensemble average in (19) is Furthermore, ŷ R k 2 = ŷ I k 2 , thus we will omit the indexes R,I from now on.
Using the main equations of the CSL model, Eqs. (10) and (11) one obtains: Substituting Eqs. (20) and (21) into Eq. (19), we find that the power spectrum can be expressed as: We observe that the prediction for the power spectrum depends on the terms ŷ 2 k and (Re[A k (η)]) −1 , which can be obtained from the CSL equations.
The quantity (Re[A k (η)]) −1 represents the variance of the field variable, which in turn is related to the width of the wave functional (13). The evolution equation for this quantity can be found by taking the time derivative of (10), and applying the resulting operator to the wave functional (13). Then, regrouping terms of order y 2 , y 1 and y 0 , the evolution equations corresponding to these terms become decoupled. In particular, the evolution equation associated to y 2 only contains A k (η), so it decouples from the other variables B k (η) and C k (η).

The evolution equation is then
Performing the change of variable A k ≡ f ′ /(2if ), Eq. (23) can be expressed as where: The solution to Eq. (24), which satisfies the BD initial condition corresponding to A k (τ ) = k/2, is where 2 F 1 is the hypergeometric function, and With solution f , one can return to the original variable A k and obtain the sought quantity (Re[A k (η)]) −1 . As a matter of fact, one has Note that if λ k = 0, then W = i for all η, and q = k. The other important term in the power spectrum is ŷ 2 k . In order to find this quantity, it will be useful to define the following objects: Q ≡ ŷ 2 k , R ≡ p 2 k and S ≡ p kŷk +ŷ kpk . The evolution equations for Q, R and S obtained from the CSL equations are: (30) with w k (η) ≡ k 2 − a ′′ /a. Therefore, we have a linear system of coupled differential equations, whose general solution is a particular solution to the system plus a solution to the homogeneous equation (with λ k = 0). In this way, the solution can be written as: where the constants C 1 , C 2 and C 3 are found by imposing the initial conditions corresponding to the Bunch-Davies vacuum state: Q(τ ) = 1/2k, R(τ ) = k/2, S(τ ) = 0. The functions y 1 and y 2 are two linearly independent solutions of y ′′ = −w k y, and the function Q p is a particular solution of The exact solutions y 1 and y 2 are (33) and y 2 = y * 1 , also k ± and b are defined in the same manner as in (27) (28) but replacing q → k.
On the other hand, the exact solution of Eq. (32) is difficult to find, but we can find approximate solutions in the regimes of interest. In particular, we are interested in the static regime, which corresponds also to the regime where the BD initial conditions are imposed. The other regime involved is the de Sitter phase, where the power spectrum is evaluated for the purpose of comparing it with the standard prediction. Therefore, in the static regime w k ≃ k 2 while in the de Sitter phase w k ≃ k 2 − 2/η 2 . It is remarkable that, in these two regimes, Q p can be approximated by the same solution, i.e.
Thus, the constants obtained from imposing the initial conditions are, Now that we have all the elements needed for obtaining the power spectrum (22), it is straightforward to check first that if λ k = 0 then P s = 0, because Q(η) = (4Re[A k (η)]) −1 exactly in that case. This result is also consistent with our view in which, if there is no collapse, then the metric perturbations are zero, i.e. there are no inhomogeneities in the spacetime.
On the other hand, considering the modes in the super-Hubble limit (−kη → 0), the power spectrum (22) can be approximated by where Note that in the definition of the amplitude A s , we have used that the quantity H 2 /φ ′2 0 = H 2 /φ 2 0 tends to a constant, given the existence of the aforementioned attractor point for the dynamical background variables in the limit −kη → 0.
Thus, we have found the main prediction of this section, namely the primordial scalar power spectrum originated by the CSL mechanism within the EU model.

III. RESULTS
In this section, we shall proceed to examine the observational effects of implementing the CSL model in the EU scenario. It will be useful to set as a reference model the one described in [163], we will refer to it as the original model.
Furthermore, as also argued in [163], we can generalize the power spectrum obtained in order to include the small scale dependence normally associated with the scalar spectral index n s . Thus, for the present section we will use the following expression for the primordial power spectrum (PPS): where k P is a pivot scale, which is traditionally set as k P = 0.05 Mpc −1 The next step in our analysis is to introduce a parameterization of λ k , in order to explore the possible observational features of our model. We propose the linear parameterization in k given by with λ 0 acting as a proportionality constant plus a parameter B ≥ 0. In fact, if B = 0, we recover a very similar expression for the PPS as the one obtained in the original model. In this way, B quantifies small deviations from the original model reflecting the inclusion of the CSL model. A very similar parameterization was also implemented in a recent work involving the CSL proposal during inflation [244]. Also, we will fix the value of the proportionality constant as λ 0 = 10 −14 s −1 . This choice is motivated by the fact that such a value is within the range allowed by laboratory experiments testing nonrelativistic versions of the CSL model [268], where λ 0 corresponds to the CSL parameter for these kind of models.
In the units used in the present paper, the former choice is equivalent to λ 0 = 1.029 Mpc −1 . Given that k has also units of Mpc −1 , the constant B has units of Mpc −1 too. For ease of notation, from now on we will neglect the units of these quantities with the understanding that the corresponding units have been well established. At this point, we would like to discuss the following issue. Numerical calculations set a restriction for implementing the exact formula (42). That is, as long as the value of k increases, the Gamma functions become exponentially small beyond the capability of machine representation. This demands a cutoff value k max , which in principle we set it to 0.015, to be consistent with the value chosen in the original model [163]. On the other hand, we can also approximate expression (42), obtaining The advantage of this approximation is that no numerical restrictions are imposed, hence it can be used in the numerical calculations up to the end of the observable window k = 1. Consequently, we have two options for performing the analysis. The first one is to use the exact expression for the PPS (42) up to k max , and then perform an analytic continuation such that it approaches smoothly to the standard expression corresponding to the canonical model, i.e. to the PPS of the standard cosmological model. The second option we can consider is to use the approximate formula (44) for the whole k range. Both alternatives can be seen in Fig. 1, which shows the PPS for different values of the B parameter comparing both the exact and approximate expressions. After an exhaustive exploration, we find that the values 10 −3 and 10 −4 are representative to show the behavior of the power spectrum curve. For each case, both calculations, exact (solid line) and approximate (dashed line), are shown.
Moreover, we will include in each Figure a plot called canonical model, which is used as a second reference (in addition to the original model). As its name suggests, the canonical model corresponds to the standard ΛCDM cosmological model, with the cosmological parameters determined by latest data from Planck collaboration [3]. In particular, we focus on the angular spectrum corresponding to the temperature and E-mode polarization autocorrelation and cross correlation functions. These data leads to the following set of parameters at the 68% confidence level: Ω b h 2 = 0.02236, Ω c h 2 = 0.1202, H today = 67.27 km s −1 Mpc −1 , A s = 2.101 × 10 −9 , n s = 0.9649, τ d = 0.0544 (called the optical depth parameter), and considering no running of the scalar spectral index. The canonical model plot will be useful not only as a guide to quickly spot the novel features introduced by our model, but also to confirm the fact that, for given a value of a 0 H 0 , a certain scale k max can be found where the PPS meets naturally the standard ΛCDM model. This was the spirit in which the value of k max = 0.015 was adopted, as mentioned above.
The next step is to analyze the effects of varying the model parameters in the angular power spectrum. To accomplish this, we perform the corresponding modifications in the Code for Anisotropies in the Microwave Background (CAMB) software [269]. Figure 2 depicts the resulting angular spectrum for different values of B.
In order to have a better visualization on the impact our model may have on observational signatures, we show the best-fit ΛCDM prediction for the angular spectrum and the corresponding Planck data 4 together with their error bars (shown in blue). The effect of varying B is mainly seen in the low multipoles (l < 50). If B tends to zero, the CSL spectrum approaches to the one of the original model, and at the same time, it separates from the canonical model at the bottom of the graph. On the contrary, if B is increased, the CSL spectrum splits from the canonical model at the top part of the plot corresponding to the lowest multipoles.
The previous analysis indicates that our model has the potential to exhibit different features in the low-l range, approaching the ΛCDM angular spectrum from upwards or below, and having the original model as the lower limit. In other words, including the CSL mechanism in the EU model could result not only in a suppression of the temperature anisotropies of the CMB at large angular scales, but also in an excess.

A. Further parameter space exploration
From the previous discussion, it is clear that the implementation of the CSL model into the EU scenario could be in good agreement with observational data. However, it is also important to check whether one can, in principle, truly distinguish the predictions between the canonical model and the one proposed in this work. Recall that we have introduced a new parameter B, while the combination a 0 H 0 comes from the original model. Therefore, we are interested in testing the robustness of B and a 0 H 0 in the predicted angular spectrum. In other words, we will vary B and a 0 H 0 enough to see how much our predicted spectrum deviates from the standard one.
At first glance, one could argue that any well fitted value of B would make the model consistent with observations, because any new features introduced in the spectrum could remain masked under the so called cosmic variance. However, exploring a wide range of the B parameter space, shows that, even though highest multipoles are unaffected, the first acoustic peak is half missed for increasing values of B, see Fig. 3. This suggests that B has an upper limit and its value can be constrained with observational data; therefore, the model has predictability. Up to this point, we have worked with the same fixed value of a 0 H 0 as the one considered in the original model [163]. Henceforth, it is interesting to explore the possibility to regard the combination a 0 H 0 as a free parameter. In Fig. 4, we vary a 0 H 0 along four orders of magnitude. Increasing a 0 H 0 , implies a strong suppression in the angular power spectrum. This fact indicates that not any value of a 0 H 0 would make the model consistent with the data. Furthermore, it also indicates that the suppression observed in the original model could be in part explained by a particular combination of the cosmological parameters that has been chosen there.
The final issue we want to address is: given a particular value of a 0 H 0 that is not compatible with observational data, could it be compensated in some way by varying the B parameter introduced by the CSL model?

IV. CONCLUSIONS
In recent years, the emergent universe (EU) has been studied from different perspectives as a viable cosmological model, which could not only avoid the horizon problem and the initial singularity of the standard approach, but could also account for some anomalies in the observations of the CMB at large angular scales.
Motivated by the recent results of [163], where it was shown that a phase of super-inflation prior to that of standard slow-roll inflation (originated in the context of the emergent universe presented in [127,130]) could explain the power suppression observed in the low multipoles of the CMB, we decided to explore in this work the same emergent universe but within the framework of a particular objective collapse theory, known as Continuous Spontaneous Localization (CSL). The CSL theory has been studied with encouraging results in cosmological applications for several years, as a mechanism to break the original symmetries of the quantum vacuum state of the fields, generating the primordial cosmological perturbations and giving an explanation to the quantum-toclassical transition of such perturbations.
We obtained a prediction for the primordial power spectrum. The computed spectrum can be consistent with CMB observational data, if a particular parameterization for the collapse rate λ k is assumed. The parametrization we considered was also proposed and analyzed in previous works but in the framework of standard inflation. Such a parameterization introduces an extra free parameter B, in addition to the parameter a 0 H 0 of the emergent model studied in [163]. We also found that the predictions and results of such a work can be recovered, when B = 0 is chosen.
From the analysis in Sec. III, we have found that implementing the CSL collapse proposal to the emergent universe scenario introduces extra modifications at the lowest multipoles. Specifically, through the CSLparameter B, the angular spectrum in the low multipoles sector, exhibits a suppression or an increment. This is a different feature from what is generically produced in models with a super-inflation phase which only decrease the curve spectrum at large angular scales. On the other hand, in the model proposed in this work, whether there is an excess or suppression, will be determined by a combination of the parameters derived from the original model and the novel one introduced by the collapse rate parameterization. This fact enables our model to introduce some new features in the angular spectrum, particularly in the sector of interest where the so-called low-l anomaly is located.
On the other hand, increasing the value of the combination a 0 H 0 produces a suppression in the lowest multipoles. We have also seen that the variation of the B parameter produces opposite (and similiar) effects to a 0 H 0 . Nevertheless, a suppression of the angular spectrum given by a 0 H 0 cannot be compensated with an increase of B .
One important aspect of our model is that the primordial power spectrum obtained smoothly approaches to the one from the standard ΛCDM model within the observational range of interest. This characteristic constitutes a natural modification of the (temperature) angular spectrum, affecting mainly the lowest multipoles without altering the highest ones (which are already constrained by observations to a high degree of accuracy).
Exploration of the free parameters shows that not any value will make the proposed model consistent with observational data. This suggests that there are good opportunities of effectively constraining the parameter space with the full machinery of statistical analysis based on Monte Carlo-Markov chains for cosmology. Such a scenario gives good predictability to the CSL collapse proposal in the emergent universe model. In fact, we expect that the extra free parameter B should take a small value (between 10 −3 and 10 −4 ) but not centered at B = 0, which clearly distinguishes our proposal from the one explored in [163]. ACKNOWLEDGMENTS G.R.B. is supported by CONICET (Argentina) and he acknowledges support from grant PIP 112-2017-0100220CO of CONICET (Argentina). G.L. and M.P.P. are supported by CONICET (Argentina), they also acknowledge support from the following project grants: Universidad Nacional de La Plata I+D G175 and PIP 112-2020-0100729CO of CONICET (Argentina). M.P.P. thanks Susana Landau for useful discussions regarding the impact of non-standard cosmological models on CMB angular power spectrum. The authors acknowledge the use of the supercluster MIZTLI of UNAM through project LANCAD-UNAM-DGTIC-132 and thank the people of DGTIC-UNAM for technical and computational support.