Puzzling initial conditions in the R_h=ct model

In recent years, some studies have drawn attention to the lack of large-angle correlations in the observed cosmic microwave background (CMB) temperature anisotropies with respect to that predicted within the standard $\Lambda$CDM model. Lately, it has been argued that such a lack of correlations could be explained in the framework of the so-called $R_h=ct$ model without inflation. The aim of this work is to study whether there is a mechanism to generate, through a quantum field theory, the primordial power spectrum presented by these authors. Specifically, we consider two different scenarios: first, we assume a scalar field dominating the early Universe in the $R_h=ct$ cosmological model, and second, we deal with the possibility of adding an early inflationary phase to the mentioned model. During the analysis of the consistency between the predicted and observed amplitudes of the CMB temperature anisotropies in both scenarios, we run into deep issues which indicate that it is not clear how to characterize the primordial quantum perturbations within the $R_h=ct$ model.


Introduction
In addition to solve the horizon and flatness problems of the standard Big Bang model, inflation generates a nearly scale invariant power spectrum for density perturbations, which has been exquisitely tested with observations of the cosmic microwave background (CMB) angular spectrum [1,2,3,4,5,6].
Starting with the Cosmic Background Explorer observations [7], it was noted that the angular two-point a e-mail: gabriel@iafe.uba.ar b e-mail: gleon@df.uba.ar correlation function at angular scales larger than 60 • is unexpectedly close to zero, contrary to what the standard ΛCDM model predicts. Shortly after, it was rediscovered with the Wilkinson Microwave Anisotropy Probe (WMAP) data [8] and later by Planck mission [2,9]. This feature at large scales was studied in detail by several authors, e.g. [10,11,12,13,14]; it has been the source of some controversy (see for instance [15]) and today constitutes one of the persistent large-angle anomalies in the CMB data [16].
Recently, a series of theoretical and observational motivations exposed in [17,18,19] finished merging into what today is known as the R h = ct model [20]. This model has received considerable attention over the last few years, since it has been claimed to be favored over the standard ΛCDM by most observational data [21,22,23,24,25,26,27,28]. Even the authors argue that the mentioned horizon problem could be solved in the framework of this model without an inflationary epoch at the beginning of the Universe [29]. Basically, they hold that the Universe can be described by a FLRW cosmology, where the cosmic fluid filling the Universe satisfies, at all times, the overall equation of state ρ + 3P = 0, where ρ and P are the total energy density and pressure of the cosmic fluid, respectively. According to the authors, the condition w = −1/3 at all times is apparently required by the simultaneous application of the Cosmological Principle and Weyl's postulate [30,31]. We remind the reader that in the standard ΛCDM model, the equation of state ρ + 3P = 0 would lead to a Universe with negative curvature.
However, some observational objections were raised, and also the validity of the physical arguments underlying the R h = ct model have been criticized by a number of authors. Some of them can be found for instance in [32,33,34,35,36,37,38,39,40]. In particular, the claim made in Ref. [41], regarding that the analysis of the CMB anisotropies in the R h = ct model is preferred over the ΛCDM, appears to be incorrect; actually the formal computation of the angular power spectrum, i.e. the C l 's, is absent (in fact, in Sect. 4 we will show explicitly that it is very unlikely that the R h = ct model can be made consistent with the CMB observational data). Furthermore, the explanation, within this model, concerning how w is kept at −1/3 through the transitions from known matter to radiation sounds at least questionable, and the idea that ρ ∝ a −2 throughout nucleosynthesis, recombination, structure formation and today seems impossible to reconcile with all the observations put together. Some of these criticisms are claimed to have been answered in [42,43,44,45,46]. Nevertheless, after pointing out a number of objections to the R h = ct model based on recent observational data, in Ref. [47] the authors analyzed the central assumption underlying the original theoretical argument for the model, namely that the comoving Hubble radius should be constant, and showed that it is not required.
In the present manuscript, we will focus on the results presented in [41,48]. There, the authors analyzed the CMB angular correlation function for a fluctuation spectrum expected from growth in a Universe, whose dynamics is constrained by the equation of state w = −1/3. To accomplish this, they mention that since the exact form of the power spectrum emerging from the non-linear growth prior to recombination is unknown, a parameterization for this spectrum can be performed, for example, by assuming a scale-free initial power law spectrum and incorporating in its shape other relevant effects. Then they ensure that it is possible to obtain a better fit than the ΛCDM model to the data corresponding to the angular correlation function, and conclude stating that the absence of power on large scales exhibited by the angular correlation function might be evidence in support of the R h = ct model simply because it does not require inflation.
In this article, we perform a critical analysis whether there might be a mechanism for generating, through a quantum field theory, the primordial power spectrum presented by those authors in [41]. To do so, we are going to consider two different scenarios: first, we will assume a scalar field dominating the early Universe in the R h = ct cosmological model, and second, we will deal with the possibility of adding an early inflationary phase to the mentioned model. After that, we will analyze the consistency between the predicted and observed amplitudes of the CMB temperature anisotropies in both scenarios.
The article is organized as follows: in Sect. 2, we review some basics about the R h = ct model and how to describe classical perturbations in that framework; in Sect. 3, we search for a quantum mechanism to generate the primordial curvature perturbation, and we obtain the primordial power spectra within R h = ct model with and without an inflationary phase. Later, in Sect. 4, we analyze the amplitudes of the primordial power spectra and the consistency with the amplitude of the CMB temperature anisotropies. In Sect. 5 we make a discussion of our results, and finally in Sect. 6 we summarize our conclusions.
2 Classical perturbations in the R h = ct model In this section, we provide a summary of the main characteristics of the R h = ct Universe. Our main focus is the cosmological perturbations as presented in Refs. [20,41,48]. The first two subsections will be heavily based on the results presented in those references. However, in the last subsection, we will show how to relate the curvature power spectrum with the matter power spectrum proposed in Refs. [41,48]. We will use units in which c = = 1. We will make use of the reduced Planck mass M 2 P = 1/(8πG) and the "West Coast" signature (+ − −−) for the metric.

The background
The R h = ct Universe is characterized by a spatially flat FLRW spacetime, which in comoving coordinates is represented by the line element Additionally, the authors of the R h = ct Universe claim that the total matter components in the Universe combined (dark matter, ordinary matter, radiation and dark energy) behave as a perfect fluid with the overall equation of state 1 where ρ and P represent the total energy density and pressure of the Universe, respectively. Therefore, the Friedmann equation H 2 ≡ (ȧ/a) 2 = ρ/(3M 2 P ) and the continuity equationρ + 3H(ρ + P ) = 0 (with the dot over functions representing derivative with respect to cosmic time t) lead to a scale factor of the form where we have normalized the scale factor to a(t 0 ) = 1 at the present cosmic time t 0 . Consequently, the Hubble radius evolves as R h ≡ H −1 = t. This is one of the main features of the R h model, i.e., the Hubble radius satisfies the relation H −1 = t during the whole cosmic evolution and not "just today" as in the standard ΛCDM model. Therefore, the total energy density of the Universe evolves as ρ ∝ 1/a 2 .

Cosmological perturbations
The dynamical evolution of the cosmological perturbations in the R h = ct Universe follows from Einstein equations δG ab = δT ab /M 2 P . In particular, by using the Newtonian (longitudinal) gauge, the Fourier modes associated to the density contrasts defined as δ k (t) ≡ δρ k (t)/ρ(t) (whereρ(t) is the background energy density) satisfÿ Therefore, by considering the equation of state associated to the R h = ct Universe, it is assumed that w = c 2 s = −1/3. Thus, the motion equation for δ k is given bÿ where ∆ k ≡ k/(aH). Note that aH = H 0 is a constant (H 0 denotes the Hubble parameter today). It is not hard to find the solutions of (4); nevertheless, the solutions depend on whether k is greater or less than aH, i.e. ∆ k > 1 or ∆ k < 1. If ∆ k > 1, then the solutions are a growing mode δ k ∼ t α (with α > 0) and a decaying mode. On the contrary, if ∆ k < 1, then the solutions are a constant and a decaying mode. In the R h = ct Universe one is primarily interested in the modes such that k > aH since these are the modes that can grow into the large scale structure. As a matter of fact, motivated by the angular correlation of the CMB and the solution corresponding to the growing modes of Eq. (4), the authors of Refs. [41,48] proposed that the initial matter power spectrum is of the form where b is an unknown constant to be adjusted, and R e is the proper distance to the last scattering surface at time t e , which corresponds to the cosmic time at the decoupling epoch. The power spectrum (5) can be recast as where A is the amplitude of the power spectrum, H ≡ aH and θ max is the maximum angular size of any fluctuation associated with the CMB emitted at t e , that is, θ max = [2πa(t e )]/[k max R e (t e )]; also k max /H = 1.
2.3 Curvature and matter power spectra in the R h = ct Universe Our next step is to relate, also through a classical analysis, the matter power spectrum with the curvature power spectrum in the R h = ct Universe. Later, we will investigate whether it is possible to find a quantum mechanism for generating the curvature perturbation. If possible, we will relate that spectrum with the matter power spectrum, and then we will compare it with the one proposed in (6). We start the discussion by switching to conformal time η, i.e. dt 2 = a 2 dη 2 . In these coordinates H ≡ aH = a ′ (η)/a(η), where a prime denotes derivative with respect to conformal time. As a matter of fact, using the equation of state P = −ρ/3, the continuity equation ρ ′ + 3H(ρ + P ) = 0 and Friedmann equation H 2 = a 2 ρ/3M 2 P , we arrive at the important result That is, H is a constant of motion in the R h = ct Universe and has the value of the Hubble parameter today. For the sake of completeness, we present the explicit form of the scale factor in conformal time coordinates: where η 0 corresponds to the conformal time today. The most generic metric associated to a flat FLRW Universe with linear scalar perturbations is where ϕ, ψ, E, B are scalar functions of the spacetime. In the Newtonian gauge, ϕ = Φ, ψ = Ψ and E = B = 0. In the absence of anisotropic stress, Einstein equations (EE) δG ab = δT ab /M 2 P lead to Φ = Ψ . Moreover, considering once again that in the R h = ct Universe c 2 s = w = −1/3, the equation of motion for the Fourier mode Φ k (η) that results from combining EE is independent solutions of (10) can be approximated by a constant and a decaying mode exp(−2Hη) ∝ a(η) −2 .
On the other hand, if k ≫ H then q ∼ k/ √ 3, and the linearly independent solutions of (10) are approximately given by a growing mode Φ k ∼ exp[(k − H)η] and a decaying mode exp[−(k + H)η] ∝ exp(−kη)/a(η) (note that the conformal time η is an increasing variable).
The EE with component δG 00 = δT 00 /M 2 P is useful to relate the density contrasts with the metric perturbation Φ. That is, As we mentioned in the previous subsection, in the R h = ct Universe one is interested in the modes such that k > H; that is, the modes whose associated proper wavelength is less than the Hubble radius. These are the modes that evolve as Φ k ∼ exp (k − H)η; consequently Φ ′ k = (k − H)Φ. By using that result, Eq. (11) becomes We emphasize that Eq. (12) is valid only for k > H and w = c 2 s = −1/3. At this point we have to do a technical digression. The quantum analysis of the field perturbations usually involves the so called Mukhanov-Sasaki variable, and then one relates that variable with the comoving curvature perturbation R. We will follow such an analysis in the next section; however, Eq. (12), which will help us to relate the matter power spectrum with the curvature one, was obtained in the Newtonian gauge. Therefore, it will be useful to change from the Newtonian gauge to the comoving gauge. That relation is generically given (for constant w) by [49] Thus, for w = −1/3, Eq. (13) leads to R = −2Φ − H −1 Φ ′ . Moreover, if we focus on the modes such that k > H and recall that for such modes Φ ′ k = (k − H)Φ k , then we arrive at With Eqs. (12) and (14) at hand, it is straightforward (in the comoving gauge) to relate the corresponding matter and curvature spectra, namely where we have retained only the first three dominant terms in powers of k/H.
Equation (15) is the main result of this subsection. One can immediately observe that if P R ∝ k −1 , then the resulting matter power spectrum will be of a similar structure as the one shown in Eq. (6), except for a constant term.
In the following section, we will attempt to construct a mechanism for generating the curvature power spectrum P R (k).

Generation of the primordial curvature perturbation
In this section, we will consider two possibilities for generating the primordial curvature perturbations: a scalar field dominating the early R h = ct Universe, and a preceding inflationary era in the R h = ct Universe.
Since in the R h = ct model the combination of different types of matter is such that it mimics a perfect fluid with an overall equation of state P = −ρ/3 (which involves a negative pressure), we will make the standard assumption that the early Universe was dominated by a scalar field φ(x, t), with some potential V (φ), such that P and ρ associated to φ satisfy P (φ) = −ρ(φ)/3 at all times. Afterwards, the scalar field should decay into particles of the standard model and possibly into dark matter particles, and the evolution of the Universe then follow the R h = ct model. Since we are considering a canonical scalar field, the action is given by In contrast with the standard ΛCDM model (plus inflation) in which the end of a different cosmological era is linked to a change in the equation of state, in the R h = ct Universe the equation of state P = −ρ/3 should be satisfied at all times during the evolution of the Universe. As a consequence, we need to provide a condition that marks the end of the early cosmological era dominated by the field φ. We propose that the value of the adiabatic speed of sound c 2 s will help to provide such condition.
For ordinary matter and constant equation of state we know that w = c 2 s . However, for a scalar field generically c 2 s = w. In particular, for a canonical scalar field (a field with canonical kinetic term), c 2 s = 1 [49]. In fact, in standard slow-roll inflation c 2 s = 1 and w ≃ −1. Therefore, in the R h = ct Universe, we will consider a canonical scalar field that dominates the matter content of the early Universe, and such scalar field will be characterized by c 2 s = 1. Then, at some point during the evolution, the scalar field will decay in such a way that c 2 s will decrease from c 2 s = 1 to c 2 s = −1/3. Note from Eq. (3) that it is crucial to have c 2 s = w = −1/3 in order to obtain Eq. (4), which results in a solution for the growing modes. It is important to mention that other combinations of w and c 2 s would lead to a solution of Eq. (3) with a growing mode; in particular, the condition for the R h = ct model is w = −1/3. Hence, other values of c 2 s , but maintaining w = −1/3 could lead to a growing mode in the R h = ct model. On the other hand, Eq. (4) is the main equation used by the authors of the R h = ct model to analyze the growth of structure in [41,48]; and to obtain Eq. (4) from Eq. (3), one must satisfy c 2 s = w = −1/3. To continue, we split the scalar field into an homogeneous part plus small inhomogeneities, i.e. φ(x, t) = φ 0 (t)+δφ(x, t). The homogeneous part of the field drives the background evolution, that is, the one characterized by the R h = ct Universe, and the quantum theory of δφ(x, t) will result in the primordial power spectrum of the perturbations. In the following, we will attempt to construct a quantum theory for δφ, but first we will derive some useful quantities to describe the background.
Since the background field, φ 0 , drives the evolution of the R h = ct Universe, we can associate the standard energy-momentum tensor T α β to the field φ 0 . In particular, from the time component . Using the fact that H is constant [see (7)], and from the continuity equation ρ ′ + 3H(ρ + P ) = 0 applied to the scalar field φ 0 , one obtains Consequently, from Eqs. (7) and (17), it is clear that in the R h = ct Universe, H and φ ′ 0 are exactly constants of motion. As a matter of fact, using the Friedmann equations, it can be shown that Furthermore, using that φ ′′ 0 = 0 and that H is a constant, we can find a potential that is consistent with P (φ) = −ρ(φ)/3. This potential turns out to be Now, let us focus on the linear scalar perturbations. The field perturbations δφ induce metric perturbations δg µν via EE. As we mentioned in Sect. 2.3, Eq. (9) represents the most generic metric associated to a FLRW Universe with scalar perturbations. As is well known, the relativistic perturbation theory has the issue of gauge redundance [50,51]. However, the gravitational part can be characterized by a single, gauge-invariant object known as the Bardeen potential defined as [52] In the same manner, the matter sector can be described by the gauge-invariant field perturbation The Einstein equations relate Φ B and δφ (gi) through a constraint equation. That implies that the scalar sector can be characterized by a single object; this object is the so-called Mukhanov-Sasaki variable, defined by All other relevant quantities can be expressed in terms of v(η, x), i.e. it fully characterizes the scalar sector. Moreover, we can expand the action of our theory, that is, the action of a scalar field minimally coupled to gravity, up to second order in the scalar perturbations, obtaining From the action (23), the equation of motion is Notice that Eq. (24) implies that z ′′ /z = a ′′ /a. Additionally, the fact that H ′ = 0 implies that a ′′ /a = a ′2 /a 2 = H 2 . Thus, the equation of motion can be rewritten as where we have defined the operator (26) is a Klein-Gordon type of equation with the "wrong" mass sign, that is, the motion equation of a free tachyon field. This can also be read directly from action (23), which is δS (2) Thus, quantizing the scalar field v(η, x) in the R h = ct Universe, is equivalent of quantizing a free tachyon with constant mass given by m 2 = −H 2 < 0. There are various methods proposed for constructing a quantum field theory of a free tachyon in the past 6 [53,54,55,56,57,58]. Nevertheless, there are some issues that seem to be always present in such theories [59]. Among them, we can mention the non-locality of the tachyonic field, represented in the present paper by the field v(x) [the short-hand notation x refers to a point in spacetime (x, η)], in the sense that the commutator (as well as the anti-commutator in some methods) does not vanish for spacelike arguments. Another puzzle is that the energy operator, normally associated to the Hamiltonian, does not have a lower bound on its spectrum, i.e. there are infinitely negative energy states, which requires some reinterpretation principle [60]. But perhaps the most serious difficulty in formulating a theory of tachyons is that the resulting S-matrix is non-unitary. Thus, it is unknown how to describe interactions within the theory of a tachyonic field [59].
In spite of the aforementioned issues, we could proceed in a pragmatic way, and construct a quantum theory of the field v(x) but only considering the modes such that k > H, i.e. modes with a proper wavelength less than the Hubble radius λ p < H −1 . Also, according to Ref. [41] those modes are the ones that can grow and evolve into large scale structure. 2 Afterwards, we could compute the quantum two-point correlation function and extract its corresponding power spectrum.
There are various known methods for constructing a quantum theory for a field with the Lagrangian (27) that ignores the "problematic modes". Among those, we can mention the one proposed by Feinberg [56] and another one developed by Arons and Sudarshan (AS) [57]. We will focus on those methods as they illustrate the kind of puzzles one encounters when trying to com- Both methods assume that the field v(x) possesses non-vanishing Fourier components only for k ≥ H and is expanded as where Then one promotes the field v into an operatorv. The difference between the AS and Feinberg's method is the quantum interpretation of the coefficients c + (k) and c − (k).
Feinberg's method follows the traditional approach of promoting c + (k) →ĉ(k) and c − (k) →ĉ(k) † into annihilation and creation operators respectively. Furthermore,ĉ(k) andĉ(k) † satisfy anti-commutator relations The anti-commutator replaces the commutator since the former is compatible with Lorentz invariance, within the quantization of a free tachyon. However, under a suitable Lorentz transformation,ĉ(k) can be converted intoĉ(k) † . Thus, the vacuum state defined asĉ(k)|0 = 0 is not an invariant vacuum state since in another frame of reference it takes the formĉ(k) † |0 = 0. For this reason, we find Feinberg's method not to be suitable for the problem at hand.
On the other hand, in the AS method both coefficients are promoted to annihilation operators. The fact that both operatorsĉ + (k) andĉ − (k) are annihilation operators is needed in this approach in order to preserve the Lorentz invariance symmetry of the vacuum state [53,59]. Moreover, one also has anti-commutation relations {ĉ ± (k),ĉ ± (k ′ ) † } = δ 3 (k − k ′ ) and the vacuum state defined as c ± (k)|0 . Consequently, we can com- In the comoving gauge, the curvature perturbation is given by R = v/z. That is, from (29) we can extract the primordial power spectrum P R (k, η) = P v (k)/z(η) 2 which, using Eq. (24), results in The previous approximated expression is valid only for k > H. As a matter of fact, P R (k, η) = 0 for k < H; i.e. there are no "super-Hubble" modes (see footnote 2). Substituting Eq. (30) into Eq. (15) yields the matter power spectrum, which is valid for k > H, while P δ (k) = 0 if k < H. The quantum theory proposed above resulted in a matter power spectrum (31) of the same structure in k, plus a constant term, as the one in Eq. (6), whose form was proposed by the authors of [41,48] motivated by observational data. It may be the case that the spectrum (31), including the constant term, could reproduce the results obtained from the one proposed heuristically in Refs. [41,48], Eq. (6), for some values of the parameters considered in those references. However, the quantum theory of the primordial perturbation in the present section contains at least two fundamental issues: (i) The theory describes a free tachyon field and (ii) the final primordial spectrum, Eq. (30), depends on the scale factor. We will study the implications of the second issue in the next section. Here, let us focus on the first issue. The fact that the spectrum obtained involved the quantum theory of a free tachyon field could discourage some readers to consider the quantum theory of the field v(x) as a serious mechanism for generating the primordial spectrum in the R h = ct Universe. The reasons are vast and we entirely subscribe to most of them. However, a possible way to deal with that issue is to abandon the R h = ct model framework for the early Universe and instead use the standard inflationary paradigm. In other words, we can assume that inflation did occur in the early Universe, but then, after the reheating era, the Universe followed the evolution described by the R h = ct Universe.
In slow-roll inflation, one has the standard theory of the inflaton field, and the end of the inflationary era is achieved when the slow-roll parameters are close to unity. As is well known, the quantum theory of inflation leads to the following expression for the Mukhanov-Sasaki variable: and z = √ 2ǫM P a, where ǫ is the standard Hubble slow-roll parameter defined as ǫ ≡ 1 − H ′ /H and during inflation ǫ ≪ 1. As a consequence, the primordial spectrum for slow-roll inflation is P R ≃ |v k (η)| 2 /z 2 , that is, where we used that H ≃ −1/η during inflation. For the "super-Hubble" modes, i.e. modes that k ≪ H during inflation, one has the familiar result (ignoring the numerical factors) that is, the scale invariant primordial spectrum, which remains constant after the "horizon crossing." On the other hand, the "sub-Hubble" modes, which satisfy k ≫ H during inflation, lead to a spectrum of the form In the standard ΛCDM model, the sub-Hubble modes are ignored since they decay as ∼ a −2 . However, in the R h = ct Universe these modes need not to decay so much because in the R h = ct model there is no "horizon problem" (see Ref. [29]). Consequently, there is no minimum value of e-folds needed for inflation to solve the horizon problem. Thus, the inflationary era could end after a few e-folds and the sub-Hubble modes still contribute to the observable modes in the CMB.
Additionally, in the inflation plus R h = ct Universe scenario, the dynamical evolution of H, ignoring the reheating era, is where η f denotes the conformal time at which inflation ends, thus, η f = −1/H 0 . Equation (36) implies that super-Hubble modes k < H during inflation, stay super-Hubble at all times, namely they do not re-enter the horizon as in the standard ΛCDM model, instead, if a mode satisfies k < H during inflation, then it also satisfies k < H during the whole evolution of the R h = ct Universe. Therefore, taking into account that Φ k is a constant for super-Hubble modes and that w = −1/3, one obtains from Eqs. (11) and (13) the matter power spectrum, where H * and ǫ * are valuated at the time −kη * = 1 during inflation. On the other hand, Eq. (36) implies that if a mode is sub-Hubble during the R h = ct Universe expansion, k > H = H 0 , then it is also sub-Hubble during inflation −kη > 1 (or equivalently k > H inf ). As a consequence, the primordial spectrum associated to these modes, Eq. (35), should be evaluated at some conformal timeη after inflation ends, namely when ǫ = 1 and a(η) > a(η f ). Therefore, after substituting Eq. (35) into Eq. (15), the matter power spectrum associated to the super-Hubble modes, at leading order in k/H, is Thus, adding a standard inflationary era to the R h = ct model, results in a matter power spectrum of the form The spectrum (39) has a resemblance to the standard prediction of the ΛCDM model but with an important difference. In the traditional model, one divides the parts proportional to k and k −3 using the value k eq. , instead of the value H 0 as in (39), where we have defined k eq. ≡ 1/η eq. and η eq. denotes the conformal time at the epoch of matter-radiation equality. If k ≪ k eq. , then the mode enters the horizon (or becomes sub-Hubble) during the matter dominated epoch and P δ (k) ∝ k. On the contrary, if k ≫ k eq. , then the modes becomes sub-Hubble during the radiation era and P δ (k) ∝ k −3 .
Also, note that the matter power spectrum shown in (39) is not of the form proposed by the authors of the R h = ct Universe [see (6)], which we obtained by adding a previous inflationary phase to the R h = ct model. The only similarity between the two expressions, Eqs. (39) and (6), is in the term that goes as k, the rest of the terms are not equivalent. The spectrum of Eq. (6) was proposed heuristically (not deduced from a physical mechanism) by the original authors of the R h = ct model in order to provide a solution for the low correlation observed at large angles. Given that Eq. (39) is not the same as (6), we cannot say if the analysis made by the authors of R h = ct model still is valid for Eq. (39), i.e. we cannot claim that the spectrum (39) solves the low large-angle correlation. In the next section, we will deepen the discussion regarding the viability of the primordial spectra obtained when considering the amplitude of the temperature anisotropies in the CMB.

Amplitude of the primordial spectra and the CMB temperature anisotropies
In the previous section, we proposed a mechanism for deriving the primordial spectrum from the quantum fluctuations of a field φ that dominated the early Universe, but with the condition that the field must satisfy P (φ) = −ρ(φ)/3. That procedure resulted, albeit the need of a quantum theory for a free tachyon, in a prediction with a similar shape to the one proposed by the authors of the R h = ct Universe [41,48], but with the difference that the final primordial spectrum showed in Eq. (30) depends on the scale factor. In the present section, we return to this subject by analyzing the amplitude of the spectrum, which is tightly related to the amplitude of the CMB temperature anisotropies.
We will consider the Sachs-Wolfe effect on the temperature anisotropies. That effect is the dominant source for the anisotropies at large angular scales (l ≤ 20). It also relates the anisotropies in the temperature observed today on the celestial sphere to the inhomogeneities in the Newtonian potential on the last scattering surface, Here, η D is the conformal time of the decoupling era and x D = R D (sin θ sin ϕ, sin θ cos ϕ, cos θ), with R D the comoving radius of the last scattering surface. It is useful to perform a multipolar series expansion δT T = l,m a lm Y lm (θ, ϕ).
Using the Fourier expansion of Φ(η D , x D ) and the expression for the Sachs-Wolfe effect, the coefficients a lm can be expressed as with j l the spherical Bessel functions of order l. The observed data is presented in terms of the angular power spectrum defined as C l ≡ 1/(2l + 1) m |a lm | 2 , that is, Eq. (41) yields It is straightforward to check that if |Φ k (η D )| 2 = A/k 3 , with A some constant, then C l = A/[l(l + 1)]. In other words, for l ≤ 20, the quantity l(l + 1)C l is a constant A and is equal to the amplitude of the squared temperature anisotropies ∼ 10 −9 [5]. Thus, it is a necessary condition that the Newtonian potential at the time of decoupling should scale as ∼ k −3 if C l is to be consistent with the temperature anisotropies of the CMB.
In the standard ΛCDM model, the value of |Φ k (η D )| 2 is determined by the modes that became super-Hubble during inflation. Those modes behave as |Φ inf k | 2 = A/k 3 and remained constant during the whole cosmological evolution up until they became sub-Hubble at some point. If the modes became sub-Hubble during the matter dominated epoch, then they remain constant even for k > H. On the other hand, the modes that became sub-Hubble during the radiation dominated epoch decayed as ∼ a 2 . Thus, in the traditional scenario once |Φ inf k | 2 is generated during inflation, it remains fixed at that value and then one simply relates |Φ k (η D )| 2 ∝ |Φ inf k | 2 . In other words, the amplitude A = k 3 |Φ inf k | 2 is fixed during inflation and is the same for all modes up to the decoupling epoch. Now, let us focus on the value of |Φ k (η D )| 2 in the R h = ct Universe. As mentioned previously, the potential Φ k corresponds to the general solution of Eq. We have defined α ≡ q/H; similarly, the second solution is given by a −α−1 . Since α > 0 the second solution corresponds to a decaying mode; on the other hand, the first solution is explicitly If k < H then α ∼ 1; on the contrary, if k > H then α > 1. Therefore, depending on whether k < H or k > H, the first linearly independent solution of (10) can be approximated by a constant or a growing mode [which is consistent with the discussion after Eq. (10)].
In the R h = ct model one is interested in the growing mode, hence k > H and α > 1. The primordial spectrum P R (k) obtained in Eq. (30) can be related to the amplitude of the Newtonian potential |Φ k | 2 through Eq. (14), which results in Note that we have evaluated the scale factor, and consequently the power spectrum, at some conformal time η p , i.e. a(η p ) = a p . At this point, we will make the assumption that the value of |Φ k | 2 obtained during the period dominated by the scalar field φ [Eq. (44)], when c 2 s = 1, is the same as the one given by (43), when c 2 s = −1/3 at the time η p . Note, however, that in both cases w = −1/3, hence the R h = ct Universe expansion remains unchanged. In particular, we are assuming that the following condition is satisfied: but w = −1/3 in both situations. In other words, we are assuming that the "reheating" period in the R h = ct Universe is practically instantaneous. Furthermore, with the condition (45) and expression (44), we can find the explicit value of the integration constant C k in (43): Consequently, the expression for |Φ k (η)| 2 in the R h = ct model is given by The Newtonian potential obtained, Eq. (47), has a scale dependence k −3 , which could guarantee the same amplitude for all the modes. But, unfortunately, it also carries an additional k-dependence through α; specifically, for modes k > H, we can approximate α ≃ k/H. Therefore, different modes, grow at a different rate, and for that reason, the amplitude of each mode at the time of decoupling would be different for each mode. Nevertheless, we could make use of the fact that up to this point a p has remained unspecified. Then, to avoid the mentioned issue, we must have a 2α p = N 2 a 2α−2 D , with N 2 some normalization constant and a D the scale factor at the time of decoupling. In other words, we are adjusting the value of a p for each mode in order to achieve that all modes arrive with the same amplitude at the time of decoupling, and thus we obtain a nearly scale invariant spectrum as observed in the CMB.
By using the expression a 2α p = N 2 a 2α−2 D and Eq. (47), we obtain the value of |Φ k | 2 evaluated at the time of decoupling, Afterwards, we could simply adjust N 2 so that can be rewritten as which implies that if α ≫ 1 (or equivalently k ≫ H) then a p ≃ a D . That is, for these modes, the epoch dominated by the scalar field φ should last up until the decoupling epoch in order to the primordial spectrum obtained can have a consistent amplitude with the corresponding observed in the temperature anisotropies of the CMB. Another way to show that the spectrum (47) presents some issues with the CMB is to calculate the angular power spectrum l(l + 1)C l from Eq. (42) using precisely Eq. (47). It is known that, for low l, say l < 20, the shape of the angular power spectrum must be essentially a constant, independent of l, which results from a nearly scale invariant primordial power spectrum. That is, the region of the angular spectrum where the Sachs-Wolfe effect is dominant must not depend on l in order to be consistent with the CMB. Thus, our next task will be to compute the angular power spectrum using the spectrum (47) for l < 20.
Assuming that a p = a D /γ, with γ > 1 a constant and evaluating Eq. (47) at the time of decoupling, we have Substituting Eq. (50) into Eq. (42), we find Since the spectrum (47) was obtained for the modes k > H = H 0 note that the lower limit of integration is H 0 . As a matter of fact, the modes with k < H vanished when we considered the quantum theory of a free tachyon. The value of R D , which corresponds to the comoving radius of the last scattering surface, was calculated in Refs. [41,48] in the context of the R h = ct model resulting R D ≃ 10/H 0 . Additionally, for k > H we can approximate α ≃ k/H = k/H 0 . Performing the change of variable in the integral (51) x ≡ k/H 0 yields Noting that the asymptotic form of the functions j l (10x) 2 contribute with a factor of 1/x 2 as x → ∞ and that γ > 1, we can conclude that the integral in Eq. (52) diverges.
We remind the reader that the standard ΛCDM prediction, which is consistent with the CMB data, would have resulted in l(l + 1)C l = const. for the lowest multipoles (approximately for l < 20). On the other hand, the result of Eq. (52) diverges even for the lowest values of l. Thus, in addition to the aforementioned problems, e.g. the need of a quantum theory of a free tachyon and that the epoch dominated by the scalar field φ should last up until the decoupling epoch, the resulting angular power spectrum is divergent.
Furthermore, the condition (49) also applies to the situation in which one drops the R h = ct model in the early Universe in favor of an inflationary era. As we described in Sect. 3, during inflation, the dynamical evolution of the Mukhanov-Sasaki variable v leads to the following expression for the comoving curvature perturbation: which for sub-Hubble modes k > H is approximated by |R k (η)| 2 ≃ (4M 2 P ǫa 2 k) −1 . Next, we evaluate |R k (η)| 2 at some time η p near the end of inflation, i.e. when ǫ = 1, and we make the assumption that once again neglecting the reheating era. Using Eq. (14), which allows us to relate |R k (η)| 2 with |Φ(η)| 2 when w = −1/3, and since we are connecting the inflationary regime with the R h = ct Universe expansion, we have H = H 0 . Consequently, the amplitude of the primordial Newtonian potential is which is essentially the same as the one in (44). Therefore, all the mathematical steps that lead from (44) up to (49) are equivalent, including Eq. (47). And, from condition (49), one is led to conclude that inflation should last until the decoupling epoch in order to the primordial spectrum obtained can have a consistent amplitude with the temperature anisotropies of the CMB.
Moreover, the discussion regarding the shape of the angular power spectrum also remains the same since the spectrum obtained in Eq. (47) will be exactly the same in the present scenario of adding an inflationary era to the R h = ct model. Given that both approaches, for generating the primordial perturbation, require very unlikely conditions to be compatible with the observed shape and amplitude of the temperature anisotropies, we could do a search for the initial value of |Φ k (η p )| 2 so that it is consistent with the CMB temperature anisotropies based solely in the dynamics of the R h = ct model.
The equation of motion for Φ k , Eq. (10), implies that That is, Φ k (η)a(η) 1−α is a constant of motion in the R h = ct Universe. Consequently, we have the relation We now select the value of the scale factor at the initial time η p . Hence, we assume that a P = Ca D , with C some normalization constant and C < 1, which from Eq. (57) yields Therefore, if the initial amplitude of the Newtonian potential is of the form then the amplitude of the Newtonian potential at the time of decoupling, Eq. (58), is |Φ k (η D )| 2 = C 2 /k 3 , which will be consistent with the amplitude of the temperature anisotropies if C 2 ≃ 10 −9 . It is evident that Eq. (59) is not equivalent to (44) and/or (55), which corresponds to the primordial amplitude obtained in the two approaches described in the previous section. Thus, any mechanism proposed for generating the primordial curvature perturbation in the R h = ct Universe must be of the form of Eq. (59) in order to be consistent with the CMB temperature anisotropies. However, neither a scalar field dominating the early Universe satisfying an equation of state P (φ) = −ρ(φ)/3 nor the inflaton yield a primordial spectrum compatible with Eq. (59).

Discussion
In this work, we began by proposing that the small inhomogeneities of a scalar field δφ could generate the primordial spectrum in the same fashion as the one in the standard inflationary scenario, but with the important difference that during the period dominated by the scalar field, the equation of state P (φ) = −ρ(φ)/3 should be satisfied at all times. That is an important condition in the R h = ct cosmological model.
Under that proposal, the quantum theory of perturbations led to a theory of a free tachyon field with mass m 2 = −H 2 = −H 2 0 < 0. We pushed forward and followed a suitable method for dealing with that kind of theory, which resulted in a matter power spectrum (31) that is similar in structure in k, plus a constant term, to the one proposed in an empirical manner in Refs. [41,48], Eq. (6). It might be the case that, for some values of the parameters corresponding to the spectrum proposed by the authors of the Refs. [41,48], the spectra (31) and (6) coincide and the analysis of Refs. [41,48] continues to be valid for the matter power spectrum given in Eq. (31). Nevertheless, as we will see in the rest of this discussion there are other important problems associated to the spectrum, Eq. (31).
The fact that the quantum theory of the perturbations in the R h = ct Universe resulted in that of a free tachyon carries deep issues. Among them, perhaps the most important issue in the cosmological context is that the corresponding S-matrix is non-unitary. Therefore, there is no clear way how to describe interactions. The interactions with other fields are important since at some point the scalar field φ, dominating the early Universe in the R h = ct model, should decay into the particles of the Standard Model, and in the absence of a well defined S-matrix, it is a puzzle how to describe such interactions. Another related issue with the non-unitarity of the S-matrix is that the self-interactions of the scalar degree of freedom, characterized here by the Mukhanov-Sasaki variable, result in primordial non-Gaussianities; however, since there is no way to characterize the interactions, one cannot quantify the amount of primordial non-Gaussianities generated by the perturbations in the R h = ct Universe.
Since in the R h = ct framework z ′′ /z = a ′′ /a, the action given in Eq. (23) will be identical for a hypothetical analysis of tensor modes. Thus, if one wanted to study the tensor case, the quantum field theory used here for scalar perturbations would be equivalent. Therefore, in the light of our results, it is not obvious how tensor modes could be generated within the framework of a quantum theory in the R h = ct model. In case of adding an inflationary epoch prior to the R h = ct evo-lution, tensor modes could be generated but since the modes satisfying k > H are relevant, their amplitudes will be exponentially suppressed and the suppressing will continue also during the R h = ct evolution.
Those issues added to the technical and conceptual problems raised by a quantum theory of a free tachyon might suggest that we should abandon the idea of describing the quantum perturbations in the R h = ct Universe.
The aforementioned problems led us to consider the quantum theory of perturbations during an inflationary era preceding the R h = ct cosmological expansion. However, given the nature of the R h = ct Universe, we needed to focus on the primordial spectrum of the sub-Hubble modes. The primordial spectrum that resulted from inflation for the sub-Hubble modes is not consistent with the matter power spectrum proposed by the authors of the R h = ct Universe.
In fact, the matter power spectrum that we obtained by adding an early inflationary regime in the R h = ct model, resembles to the traditional one from the ΛCDM model, but it is not exactly the same. More precisely, the matter power spectrum obtained by adding an early inflationary regime to the R h = ct model, can be separated into two cases. In one case the spectrum goes as k, and in the second case the spectrum goes as k −3 . In the ΛCDM model the matter power spectrum can also be separated into two cases, one that goes as k and a second case where the spectrum goes as k −3 . However, as mentioned after Eq. (39), the conditions for the separation into two cases in the R h = ct model are not the same as the conditions in the standard model. Consequently, the functional form of the ΛCDM matter power spectrum is completely different from the R h = ct model. Moreover, the matter power spectrum obtained by adding an early era of inflation to the R h = ct model leads to a matter power spectrum that is different from the one proposed empirically by the authors of the R h = ct model. The main motivation, as stated by those authors when proposing such a spectrum, was to solve the observed low correlation at large angles in the angular correlation, and since the spectrum that we have obtained, Eq. (39), is not equal to the one heuristically proposed, Eq. (6), we cannot claim that the spectrum in Eq. (39) explains the observed low angular correlation at large angles.
Finally, we investigated the predicted amplitude of the CMB temperature anisotropies following the two approaches described. In the first approach, we considered an early phase dominated by a scalar field φ satisfying the equation of state of the R h = ct model; in the second framework, we assumed an inflationary stage preceding the R h = ct cosmological evolution. In both approaches, the amplitude of the Newtonian potential at the time of decoupling, which is the main source of the temperature anisotropies at low angular multipoles, depends on the wavenumber k in a nontrivial way. The reason is that in the R h = ct Universe the evolution of each mode associated to the Newtonian potential evolves as ∼ a α−1 , with α ≃ k/H. As a consequence, we were forced to choose a particular initial condition for the evolution of the modes that translates into adjusting the value of the scale factor for each mode at the initial time η p . Specifically, we had to choose a 2α p = N 2 a 2α−2 D , with N some normalization constant and a D the value of the scale factor at the time of decoupling. Such an election implies that, for modes k ≫ H, the initial value of the scale factor is a p ≃ a D . That is, the era dominated by the scalar field φ, in the first approach, or the inflationary era preceding the R h = ct evolution, should last up to the decoupling epoch in order to the primordial spectrum obtained can have a consistent amplitude with the corresponding observed in the temperature anisotropies of the CMB. Another related problem that emerges from considering the Newtonian potential Eq. (47) is that the predicted angular power spectrum Eq. (52) diverges in the region where the Sachs-Wolfe effect is dominant. The expected behavior for the angular power spectrum, which is consistent with the CMB data in the Sachs-Wolfe region, is l(l + 1)C l ≃ constant.
We ended our analysis by obtaining the desired form of the initial amplitude of the Newtonian potential based solely on the dynamics of the R h = ct Universe, and that it could be consistent with the CMB temperature anisotropies. That primordial amplitude should be |Φ k (η p )| 2 = C 2α /k 3 , which is not the one obtained from the two approaches considered so far. Thus, we think that any proposed mechanism for generating the primordial spectrum should predict that particular initial amplitude to be consistent with the observed CMB anisotropies. In fact, if some physically motivated mechanism, within the R h = ct model, can reproduce the primordial spectrum |Φ k (η p )| 2 = C 2α /k 3 then all the concerns raised in our paper would possibly disappear. Nevertheless, neither the inflaton nor a scalar field dominating the R h = ct Universe can produce that kind of spectrum.

Conclusions
Some studies have drawn attention to the lack of largeangle correlations in the observed CMB temperature anisotropies with respect to that predicted within the standard ΛCDM model. Recently, some authors have suggested that this lack of correlations could be explained in the framework of the so-called R h = ct model without inflation, by selecting an explicit form for the matter power spectrum and showing that it could achieve a better fit than the ΛCDM model to the data corresponding to the CMB angular correlation function. The aim of this work was to critically investigate whether there may be a mechanism to generate, through a quantum field theory, the primordial power spectrum presented by these authors.
During this search we run into deep issues and, we also studied the possibility of adding an inflationary phase prior to the evolution given by the mentioned model. The resulting power spectrum for the relevant sub-horizon modes within this approach is not consistent with the matter power spectrum displayed by the mentioned authors; thus, it cannot explain the unexpectedly close to zero angular two-point correlation function observed at angular scales larger than 60 • .
Also, we analyze the consistency between the predicted and observed amplitudes of the CMB temperature anisotropies with and without the inflationary epoch added prior to the R h = ct evolution. We found that for modes satisfying k ≫ H, the epoch dominated by the scalar field (representing the matter field in the R h = ct Universe or the inflaton) should last up to the decoupling epoch in order to the primordial spectrum obtained can have a consistent amplitude with the corresponding observed in the temperature anisotropies of the CMB. That is an implausible condition. Additionally, we have performed a brief analysis by focusing on the lowest angular multipoles l < 20, where we expect the Sachs-Wolfe effect to be the dominant effect, and we obtained, for l < 20, the angular power spectrum l(l + 1)C l is divergent; clearly not consistent with the observations from the CMB.
Finally, we showed the generic form that a primordial curvature power spectrum should exhibit in the R h = ct framework to be consistent with the CMB temperature anisotropies observed. Neither a scalar field dominating the early Universe satisfying an equation of state P (φ) = −ρ(φ)/3 nor the inflaton yield a primordial spectrum compatible with this requirement. Based on the results obtained in this paper, we conclude that (in addition to the criticisms already raised by other authors) it is not clear how to characterize the quantum perturbations within the R h = ct Universe, rendering this model a very unlikely alternative to the standard ΛCDM model.