Inflationary quantum dynamics and backreaction using a classical-quantum correspondence

We study inflationary dynamics using a recently introduced classical-quantum correspondence for investigating the backreaction of a quantum mechanical degree of freedom to a classical background. Using specifically a coupled Einstein-Klein-Gordon system, an approximation that holds well during the very early inflationary era when modes are very deep inside the Hubble horizon, we show that the backreaction of a mode of the quantum field will renormalize the Hubble parameter only if the mode's wavelength is longer than some threshold Planckian length scale. Otherwise, the mode will destabilize the inflationary era. We also present an approximate analytical solution that supports the existence of such short-wavelength threshold and compare the results of the classical-quantum correspondence with the traditional perturbative-iterative method in semiclassical gravity.


Introduction
The existence of an inflationary era that precedes the standard radiation-(nonrelativistic) matter-dark energy expansion history is a powerful hypothesis that simultaneously solves the horizon, flatness, and monopole problems in cosmology [1][2][3][4]. Most importantly, inflation brings into the overall cosmological picture a quantum mechanical factory of inhomogeneities in the early Universe that could later on seed structure formation [5][6][7][8][9][10][11][12]. Observational signatures of primordial inflation can be found in the cosmic microwave background and a e-mail: rbernardo@nip.upd.edu.ph possibly in the future through gravitational wave observations [13][14][15][16][17][18][19][20][21]. The analysis of inflationary particle production is typically done within the confines of semiclassical gravity, wherein the modes, or rather particles, of a quantum field are drawn out of the vacuum by the expanding spacetime. However, the eventual response of the background spacetime to an accumulated cosmic fluid -the backreaction -is mostly restricted, if not completely overlooked, because of understandable technical reasons. Simply put, at the linear order formalism, the classical background merely acts as a fixed stage on which the perturbations behave as test fields. The backreaction effect causing the stage to respond to the fields can therefore be found at the higher order, nonlinear interactions. Several breakthroughs have been made even with the existing challenges and the analysis of quantum backreaction during inflation can be expected to remain as an active field in forthcoming years [22][23][24][25][26][27][28][29][30][31][32][33][34]. Inflationary particle production and backreaction must be cautiously dealt with all the same and any method for shedding light to this puzzle even in simpler analogphysical systems should be worth investigating.
A classical-quantum correspondence (CQC) was recently introduced to study the backreaction of a quantum harmonic oscillator (QHO) to a classical background [35]. Indeed, the dynamics of a quantum mechanical degree of freedom (d.o.f.) could be treated simply by using the classical Bogoliubov coefficients, leading to Hawking radiation and inflationary particle production. The CQC, on the other hand, shows that a complexified classical harmonic oscillator (CHO) could be used, instead of the Bogoliubov coefficients, to look at the backreaction of a QHO to a classical background. This complexified oscillator retains the quantum mechanical d.o.f.s just as much as the Bogoliubov coeffi-cients but it is arguably simpler, and more pedagogical, because it can be treated using integration methods familiar for dealing with the CHO. We clarify that the above methodology for dealing with the backreaction, i.e., utilizing a complexified CHO to deal with the coupled QHO-background dynamics, is what is meant by classical-quantum correspondence throughout this paper. In particular, we emphasize that this is different from perhaps synonymous terms pertaining to the important transition of quantum excitations to classical modes due to an expanding classical background [36][37][38]. Physical insights on backreaction have already been drawn by using the CQC in its early applications to a rolling ball-QHO system [35] and Hawking radiation [39] but it is the elegance by which these important results have been obtained that speaks loudly about the potential of the method. Moreover, the CQC has been extended to cover field theory applications [40] and applied for instance to a coupled scalar field system [41], one of which is the rolling inflaton familiar in the inflationary context, and the evaporation of breathers [42]. It is important to stress out that the CQC is an exact solution to the field equations and that perhaps its only limitation is whenever the coupled classical-quantum d.o.f.s picture no longer holds. Applying the CQC to more systems in order to obtain physical insights on quantum mechanical backreaction and to assess CQC's regime of validity is therefore of considerable importance. In this work, we apply the CQC for the first time in a cosmological background, specifically on inflationary dynamics.
We consider a coupled Einstein-(massless) Klein-Gordon system (Eqs. (48) and (49)), wherein the scalar field is quantized and driven by the cosmological dynamics. This model is admitedly not a faithful representation of the entire history of the early Universe but it holds well during the very early stages of inflation when the modes are deep inside the Hubble horizon [43]. At the same time, it appears that this model is one that is tailored for an exhibition of the CQC in the inflationary setting. Two specific reasons make it particularly appealing. First, it is the simplest and most natural testbed for analyzing the semiclassical effects of gravity. It is also quite flexible in the sense that the case of a massive scalar field can be accommodated by replacing the wavenumber k by √ k 2 + m 2 where m is the bare mass of the quantum field. Second, the d.o.f. of a massless scalar field can be used to describe the d.o.f.s of higher spin fields, e.g., the two independent polarizations of the massless tensor field. The results obtained in this paper could therefore be used to also gain insight on the backreaction of the gravitational waves produced during inflation. The same model, however, can-not be used to describe the production of fermions. This calls for an independent work using a coupled Einstein-Maxwell-Dirac system.
The main objective of this work is to demonstrate the relative ease of implementing the CQC and to gain insight on backreaction in an inflationary setting. The rest of this work proceeds as follows. We start by revisiting the CQC in detail and its application to the rolling ball-QHO system (Section 2). This prepares the application to inflationary particle production (Section 3). We recast the field equations of the coupled Einstein-Klein-Gordon system into a form that can be objectively studied using the CQC (Section 3.1) and present the results of various numerical integrations (Section 3.2). We furthermore compare the CQC to the perturbativeiterative method (Section 4) and discuss the relevant insights and inherent limitations of the present work (Section 5).
We work with the mostly-plus metric signature (−+ ++) and Planckian units c = 8πG = = 1 where c, G, and are the speed of light in vacuum, Newton's gravitational constant, and Planck's constant, respectively (Appendix A).

Classical-quantum correspondence
In this section we provide a detailed overview of the classical-quantum correspondence (CQC) and revisit its application to a toy problem. This prepares the application to the cosmological background.

An overview of the CQC
The Hamiltonian for a driven QHO with time-dependent frequency is given bŷ where ω (t) is the classical driving frequency or the background. As is well-known, this Hamiltonian can be diagonalized in the Heisenberg picture by writing it aŝ whereâ andâ † are the annihilation and creation operators, respectively, i.e., a, a † = 1, given bŷ Solving operator equations of motion, however, is often extremely challenging 1 . For this reason, the quantum dynamical problem becomes more manageable after using the Bogoliubov transformation where α and β are classical variables andâ 0 =â (0) and a † 0 =â † (0) are the initial values of the annihilation and creation operators. The Heisenberg equations of motion for the operatorsâ andâ † becomė and it can be shown that the Bogoliubov coefficients α and β satisfy Furthermore, in terms of the Bogoliubov coefficients, the vacuum expectation value of the Hamiltonian can be written as where |0 is the vacuum state at t = 0, i.e.,â 0 |0 = 0. This quantifies the excitation driven by the background ω (t). Quite remarkably, it was shown in Ref. [35] that an improvement, and arguably a simplification, to the traditional Bogoliubov transformation method can be further obtained by using the following transformation which trades the Bogoliubov coefficients α and β for a single complex variable z. This variable z describes a complex classical harmonic oscillator (CHO), and in terms of which the original QHO problem can be shown to be equivalent to that of two CHOs with constrained initial conditions. Indeed, it can be shown 1 The dynamics of an operatorÂ is given by i dÂ dt = Â ,Ĥ + i ∂Â ∂t whereĤ is the system's Hamiltonian.
that the real and imaginary parts of z independently satisfy the CHO equation of motion, ξ + ω 2 (t) ξ = 0 (16) χ + ω 2 (t) χ = 0, (17) but with the constrained initial conditions (coming from α (0) = 1 and β (0) = 0) ξ 0 = 0 (18) ξχ − χξ = 1/2 (22) where ω 0 = ω (0). Obviously, ξ and χ are also CHOs where ξ is initially pushed to the right away from the origin and χ is initially maximally stretched to the left and released from rest. In other words, the CHOs ξ and χ are out of phase by π/2. Furthermore, in terms of the CHOs, the vacuum expectation value can be written simply as the sum of the energies of two driven CHOs, i.e., where E ξ and E χ are given by Shortly, we shall revisit the application of the above transformation to the problem of a QHO coupled to a classical rolling ball (Section 2.2). It is possible work out the direct transformation from the original QHO operatorsx andp to the CHOs ξ and χ 2 . The resulting transformation is given bŷ wherex 0 andp 0 are the initial values of the position and momentum operators, respectively, at t = 0.
2.2 Application to a toy problem: classical rolling particle and QHO Consider the toy problem described by the Hamiltonian 2 To see this, one should simply track the transformations used throughout. For instance, this procedure leads tox = z * â 0 + zâ † 0 andp = m ż * â 0 +żâ † 0 which leads straight to Eqs. (26) and (27).
In this problem, the QHO d.o.f. is given byx and the classical rolling particle d.o.f. is given by y. Also, ω 0 and a are constants that would be the QHOs natural frequency and the rolling ball's acceleration in the uncouplied limit. It should be stressed out that the classical terms in Eq. (28), e.g., p 2 y /2, should always be viewed besides the identity operatorÎ for consistency with the explicitly operator terms.
To proceed using the CQC, we write down the Hamiltonian in terms of CHOs described by Eqs. (26) and (27). This straightforward calculation leads tô The classical equation of motion for y can be derived from this simply first by taking the vacuum expectation value and using the Hamiltonian equations, This leads tö Clearly, this shows that the rolling ball would be uninterupted and fall with constant acceleration a provided that there's no coupling with the QHO (λ = 0). The backreaction of the QHO to the background y is therefore quantified by the terms attached to λ in Eq. (32). Before this can be used, however, it is important to specify the initial shape of the wave packet as as this is a quantum dynamics problem 3 . Using Eq. (33), it can be shown that 0|p 0x0 +x 0p0 |0 = 0 (36) and Eq. (32) reduces tö The complete dynamics incorporating the backreaction of the QHO to the background can finally be obtained by solving the coupled classical system given by Eq. (37) and the CHO equations of motion together with the initial conditions The vacuum expectation value could also be simply calculated by adding the energies of the rolling ball and the CHOs, i.e., In Ref. [35] it was shown using this system of equations that backreaction renormalizes the bare acceleration a of the background.

Backreaction of quantum scalar field on a cosmological background
In this section, we use the CQC to study the evolution of a mode in the cosmological background, where a is the scale factor. We setup the field equations in Section 3.1 and present the results in Section 3.2.

Setup: Einstein-Klein-Gordon system
We rely on the coupled Einstein-Klein-Gordon theory. In this way, the quantum fieldφ satisfies the (massless) Klein-Gordon equation and the metric (the background) satisfies the Einstein equation where variables with hats, e.g.,φ, are (quantum mechanical) operators. Also, in Eq. (49), the quantity 0|T αβ |0 is the vacuum expectation value of the scalar field's stress-energy tensor (SET),T (φ) ∼φ 2 , ∂φ 2 , and T (M ) αβ is the matter sector's SET. To support the inflationary era, we assume that the matter sector is completely dominated by vacuum energy, a barotropic perfect fluid with an equation of state w = P/ρ = −1, where ρ and P are the energy density and pressure, respectively 4 . On the other hand, the SET of the quantum scalar field is given bŷ which could be identified as a fluid with an energy densityρ (φ) and pressureP (φ) given bŷ This identification also shows that particles driven out of the vacuum by the expansion could be described by an effective equation of state, w φ = P (φ) / ρ (φ) = 1, where the averaging, · · · , is taken over a quantum state. Thus, because cosmic fluids with an equation of state w > −1/3 contribute to deceleration, it is reasonable to expect the production ofφ-particles and their backreaction to the spacetime would slow down the ongoing expansion. An important comment is in order. First of all, the object 0|T αβ |0 sourcing the response of the background (encoded in the Einstein tensor G αβ ) in Eq. (49) is a generally divergent quantity due to the sum of albeit finite but infinitely many zero point energies coming from each of the modes ofφ(x). This divergence can be dealt with using renormalization and has been the topic of many important papers in the field [44][45][46]. On the other hand, we shall deal with only a mode of the quantum field in what follows to obtain finite predictions of the backreaction per mode. This goes well by the scope of the paper which is to demonstrate the application of the CQC for the first time in an expanding spacetime background. However, the renormalization of the zero point fluctuations within the context of backreaction and the CQC will have to be confronted in a future work.
We proceed to setting up the dynamical system. To achieve this our first goal is to rewrite Eqs. (48) and (49) into a QHO driven by an expanding spacetime. In the cosmological background (Eq. (47)), the Klein-Gordon equation reduces tö and the Einstein equation (Eq. (49)) turns into the Friedmann equations where H is the Hubble parameter, H =ȧ/a, and an overdot indicates differentiation with respect to the comoving time t. Recognizing that the Klein-Gordon equation (Eq. (52)) must be restructured as a driven QHO problem to use the CQC, we deliberately express the field equations in terms of the conformal time η and conformal scalar fieldψ defined by 5 and respectively. This leads tô where a prime denotes differentiation with respect to the conformal time η. The above equation can easily be recognized as the QHO's equation of motion for the 5 In terms of the conformal fieldψ = aφ and the comoving time t, the Klein-Gordon equation is given by On the other hand, the transformation from comoving time t to conformal time η is given bÿ where H =ȧ/a = a /a 2 .
fieldψ at a given location x = X. We must also express the Friedmann equations in terms of the conformal variables η andψ. Straightforward calculation leads to where the energy density and pressure of the scalar field becomê This achieves our first goal. Eqs. (60), (61), and (62) indeed make up a system that could be solved using the CQC. The above set of equations, however, while it could be considered already as a QHO problem coupled to a dynamical background, contains infinitely many QHOs defined at each point x = X. To simplify the analysis, we instead focus only on a single k-mode,ψ k (η), of the fieldψ (η, x) defined by the Fourier transform In this way, we work instead with modesψ k (η) which satisfy the driven QHO equation of motion with the time dependent frequency The energy density and pressure of a mode are given bŷ This sets the stage for the CQC. Noting that each modê ψ k is described by the Lagrangian the Hamiltonian is therefore given bŷ where is the field conjugate to ψ k . Analogous to Eqs. (26) and (27), we implement the CQC by performing the transformation where ψ k0 ,p k 0 = iδ kk and ψ k0 ,p k0 are the operators ψ k ,p k at η = η 0 . This converts the driven QHO dynamical problem (Eq. (65)) into that of two driven CHOs ξ and χ described by where ω 0 = ω (η 0 ). In terms of the CHOs, the vacuum expectation value of the mode HamiltonianĤ k can be written as where E ξ and E χ are the total energies of the CHOs. The Friedmann equations (Eqs. (61) and (62)) and Eqs.
(73) and (74) form the CQC dynamical system for the vector (a, ξ, χ) which describes cosmological dynamics incorporating single k-mode particle production and backreaction. We complete the setup by expressing the vacuum expectation values 0|ρ (ψ) |0 and 0|P (ψ) |0 in terms of the CHOs (ξ, χ) and preparing the system in some initial state. To do so, we note that Eqs. (71) and (72) straighforwardly lead to and  61) and (62)) in terms of the CHOs ξ and χ. Lastly, we prepare the initialψ k -configuration in a Gaussian with a width σ, i.e., In this way, Eqs. (80), (81), and (82) further reduce to and the vacuum expectation values ofρ Some final notes are in order to guarantee a wellposed initial value problem. First, by using Eq. (66) and the Friedmann equations (Eqs. (61) and (62)), we can express ω 0 as where a 0 and 0|ρ (ψ) |0 0 are the initial scale factor and the scalar field's energy density, respectively. The above equation also shows that the scalar field sources its own frequency -a backreaction effect -and that this effect is opposite to one where only the classical expansion drives the quantum field. It can also be shown that the initial value of the vacuum expectation values of the energy density and pressure of the modes are given by 0|ρ (ψ)

Summary
The Einstein-(massless) Klein-Gordon system for a kmode of the quantum scalar field that is initially Gaussian distributed with a width σ is equivalent to the coupled system and where the driving frequency ω (η) is given by and the initial conditions are

Numerical investigation
The numerical execution is described as follows. We combine the Friedmann constraint (Eq. (92)) with the Hubble equation (Eq. (93)) to obtain Also, we fix σ = ω 0 throughout, thus, locking the initial state of the QHO to its instantaneous ground state. Previous implementations of the CQC have always invoked this. Remarkably, the main conclusion of this paperthe existence of a short wavelength threshold for inflationary stability -is valid even for σ = ω 0 . Eqs. (94),   (102)) can be computed using these parameters, i.e., calculate ω 0 using Eqs. (90) and (91).
tion of this shows that the scalar field's energy density, ρ φ = 0|ρ (ψ) |0 , has a minimum value of ρ φ ≈ 1.00 at k ≈ 1.41 and that ρ φ < 2.02 for 0 ≤ k < 3.90. It is also worth pointing out that the vacuum energy, ρ, assumes negative values for k 5.92. The destabilization of the inflationary era might therefore be naively associated with this. However, it will be shown for the nearly critical wavelengths (k ∼ 5.78) that this does not appear to be a straightforward conclusion. The succeeding discussion is hereby divided into long wavelength-modes (Section 3.2.1), nearly critical wavelength-modes (Section 3.2.2), and short wavelength-modes (Section 3.2.3).

Long wavelength-modes
We present results of the integration for long wavelengthmodes. Figures 2(a) and (b) show that the background spacetime eventually ends up with an exponential expansion. This is more convincingly observed in figure  2(b) where the logarithm of the Hubble parameter (H = a /a 2 ) undoubtedly becomes constant, marking the onset of (a renormalized ) inflationary era, during the late stages of the evolution. The backreaction can indeed be confirmed to renormalize the Hubble parameter to a smaller value H eff < h, where H eff = H (η → η f ) with η f corresponding to t → ∞, for all k < 5.78. This defines the long wavelength regime. In addition, figures 2(a) and (b) explicitly show that k = 1 backreacts to the spacetime more strongly compared to k = 2, i.e., H eff,k=1 < H eff,k=2 < h. This can be explained by the minimum value of ρ φ for k ≈ 1.41 (figure 1) and in particular because ρ φ (k = 1) > ρ φ (k = 2), noting that ρ φ (k) can be taken to quantify the strength of backreaction of a mode k. The results of the integration, (a, ξ, χ), could of course be used to compute any desired physical quantity. Shown in figures 2(c), (d), (e), and (f) are therefore the corresponding QHO Hamiltonian 0|Ĥ|0 , frequency squared ω 2 , particle number |β| 2 , and scalar field's energy density 0|ρ (ψ) |0 , respectively. It can be seen that the QHO's Hamiltonian (figure 2(c)) as well as the square of the frequency (figure 2(d)) follow a decreasing trend in their evolution, even crossing to negative values for the longer wavelengths. The trail to negative Hamiltonian shows that the long wavelength modes are themselves unstable and so they must interact with other fields later through some unaccounted mechanism. Nonetheless, particle creation continues to occur as shown in figure 2(e). Mode production, marked by the eventual rise in the particle number, |β| 2 , occurs earlier for the longer wavelength cases, as they should be easier to pull out of the vacuum. This is clearly the case in figure 2(e) where |β| 2 begins to exponentially rise earlier for k a compared to k b > k a , e.g., for k = 1 the onset of mode production occurs near η = 0.04 while for k = 2 this occurs near η = 0.36. The scalar field's energy density, 0|ρ (ψ) |0 , declines despite this because of  the ongoing expansion and anticipates the upcoming inflationary era (figure 2(f)).

Near threshold wavelength-modes
We present the results of the integration of the field equations for nearly critical wavelength-modes, i.e., modes characterized by wavelengths very close to the threshold. Figures 3(a) and (b) show the drastic shift in the behavior of the background, from inflationary to collapsing, upon crossing a certaint cutoff wavenumber, k max ≈ 5.77. This defines the threshold wavelength λ min = 2π/k max which conveniently differentiates the long wavelength-modes (λ > λ min ) which will only tend to renormalize the Hubble parameter from short wavelengthmodes (λ < λ min ) that will cause the spacetime to collapse. It should be stressed out that ρ > 0 for all of  the cases considered in figure 3 and so the eventual collapse for the short wavelength-cases (λ < λ min ) cannot be straightforwardly attributed to a negative value of vacuum energy. In particular, it can be shown by numerical inspection of figure 1 that the nondynam-ical vacuum energy, ρ, will become negative only for k 5.92. The near threshold wavelength-modes can also be seen to support a positive, slightly decreasing, Hamiltonian (figure 3(c)) and squared frequencies (figure 3(d)) throughout their evolution. They can there-fore be considered to be more stable compared to the long wavelength-modes which cascade to negative energies. Particle creation, however, cannot anymore be observed to continue because of the stronger backreaction (figure 3(e)); instead, the particle number, |β| 2 , just oscillates up and down for the near threshold-cases. The scalar field's energy density, 0|ρ (ψ) |0 , on the other hand, increases during the late stages of the evolution for the nearly critical wavelength-modes with λ < λ min as shown by careful inspection of figure 3(f) for k = 5.78 and k = 5.79. This is of course understood to be an effect of the collapsing era on the scalar field's energy density.
In a short while, it will be shown that the existence of the threshold wavelength can be crudely predicted analytically by assuming that the frequency ω is constant throughout (Section 3.2.4). This rough approximation (Eq. (105)) holds well for very short wavelengthmodes and supports the existence of a threshold.

Short wavelength-modes
We present the results of the integration for the short wavelength-modes. Figures 4(a) and (b) for the scale factor and the Hubble parameter, respectively, confirm that the spacetime collapses due to the backreaction of short wavelength-modes. For all cases in figure 4, however, the vacuum energy is negative and so this might be the reason why the integration leads to an eventually collapsing spacetime background. The discussion of the nearly critical wavelength-modes, on the other hand, have shown that even for ρ > 0 it can be inevitable to end up with a collapsing spacetime background because of the backreaction of the quantum field 6 . Figures 4(c) and (d), moreover, show that short wavelength-modes are supported by a roughly constant Hamiltonian and frequency and are therefore more stable compared to the long wavelength and nearly critical wavelength counterparts. At the same time, the scalar field's energy density (figure 4(f)) increases notably because of the collapse during the late stages of the evolution. This is accompanied by a production of modes ( figure 4(e)). It must also be mentioned that the larger k is, or rather the shorter λ < λ min is, the earlier the spacetime collapses due to backreaction.

An analytical solution for short-wavelength modes
The assumption of constant frequency, which is supported by the results for the short wavelength-modes 6 The precise solution to this apparent puzzle will be discussed in a different paper investigating the initial conditions of inflation in light of the CQC.
( figure 4(d)), can be used to obtain an approximate analytical solution and moreover deduce the existence of the threshold. This calculation can be done starting with Eq. (96) and taking the left hand side to be a constant that is completely determined by Eqs. (90) and (91). The vacuum energy can be determined by using the Friedmann constraint (Eq. (61)). This leaves us with a second-order linear ordinary differential equation for the scale factor, a + ω 2 0 − k 2 a = 0, which incidentally is the harmonic oscillator's equation of motion itself. By using the same initial conditions as with the semiclassical solution (Eq. (104)), we obtain the approximate short wavelength-mode solution This analytical expression (recall that sinh (x) , cosh (x) ∼ e x , x → ∞) shows that a de Sitter state, a.k.a., an inflationary era, would survive at late times even with backreaction provided that the quantity f (k) is real or rather that Fixing σ = ω 0 (k), this therefore predicts a cutoff wavenumber, k = k max , defined by or rather a threshold wavelength, λ min = 2π/k max , such that modes with k > k max would destabilize the background to collapse. This approximation understandably predicts k max ≈ 3.87 instead of k max ≈ 5.77 because the accuracy of this approximation is expected to improve only for large k, where the assumption of a constant frequency holds to a certain degree. Nonetheless, it supports the conclusion that there should be a threshold Planckian wavelength. Plots of the scale factor given by Eq. (105) and the CQC are shown in figure 5 for k = 6 and k = 8. This confirms that the approximation (Eq. (105)) only starts to become better for larger k, which is to be expected based on the validity of the constant frequency assumption.

CQC versus perturbative-iterative method
Prior to the CQC, a common way of obtaining semiclassical gravity effects is by solving for the Bogoliubov coefficients provided a fixed classical background. The solution to this, i.e., the Bogoliubov coefficients, could be used to obtain the leading order (LO) backreaction effect to the background which henceforth could be used to obtain the modes at next-to-leading order (NLO), and so on. This is the perturbative-iterative method.
It was shown that the CQC is the large iteration limit of the perturative-iterative method in the context of the rolling ball-QHO toy problem [35]. The same assertion was also argued to be the case for CQC's field theory counterpart [40]. In this section, we show that this is also the case for the Einstein-Klein-Gordon system. To obtain the perturbative-iterative results, we start with the de Sitter background given by Eq. (104) and from this solve for (ξ, χ) using Eqs. (94) and (95). We stress that this is the familiar elementary step of solving for the Bogoliubov coefficients, provided a fixed background, but expressed more conveniently in terms  (93)) for the scale factor a. This scale factor can then again be used to obtain (ξ, χ) which will determine a better approximation to a at next-to-leading order (NLO). This procedure can be continued for any practical order in iteration to obtain a better approximation to the exact solution. Figures 6 and 7 show the LO and NLO results for the scale factor, a, and the scalar field's energy density, 0|ρ (ψ) |0 , together with the CQC. All of the cases considered indeed strongly confirm that perturbativeiterative method approaches the CQC. Remarkably, even at NLO (N = 2 iterations) only, the agreement between the scale factor and the CQC counterpart can be considered to be satisfactory, most especially for long wavelength modes where it becomes impractical to visually distinguish the NLO result and CQC. In particular, the relative deviation of the NLO from the CQC for the scale factor peaks only at 0.25%, 0.10%, 0.80%, and 2.5% for k = 2, 4, 6, 8, respectively. The subpercent error level explains why the plots of the NLO and CQC cannot be visually distinguished, hence, the need to show the inset plots of relative error in figure  6. A similar result can be observed between the energy density, 0|ρ (ψ) |0 , obtained at NLO and the CQC, although with a notably larger deviation as compared with the scale factor results. This, nonetheless, continues to support the assertion that the CQC is the large iteration limit of perturbative-iterative method as the NLO can be observed to approach the CQC in all cases. The visual indistinguishability of the NLO and CQC results, especially for longer wavelengths, demands again the need to show the inset of relative error in figure 7. At NLO, thus, the peak relative error from the CQC is at 15%, 10%, 30%, and 60% for k = 2, 4, 6, 8, respectively. The relative deviation clearly increases for shorter wavelengths (larger k) where backreaction becomes more significant. This calls for larger iteration orders in the perturbative-iterative method for short wavelength-modes.
Building on the NLO, we demonstrate that the CQC is the undeniable limit of the perturbative-iterative method by comparing the tenth order (N = 10) iteration and the CQC for the scale factor, a, and the scalar field's energy density, 0|ρ (ψ) |0 . This result is shown in figure 8 for k = 8 which was previously where the NLO can be regarded as an awkward approximation. The subpercent agreement (δa, δρ ψ ∼ 10 −5 %) shown in figure 8 strongly supports the CQC as the limit of the perturbative-iterative method.

Discussion
We have reviewed the CQC in detail (Section 2) and, for the first time, successfully applied it to a quantum field that is coupled to a cosmological background (Section 3). In addition, we have shown that the CQC is the undeniable limit of the perturbative-iterative method in semiclassical gravity (Section 4), supporting previous assertions of this statement in the rolling ball-QHO system [35] and CQC's field theory extension [40]. This is also a good place to remind that we have not dealt with renormalization thus far and have considered only the individual dynamics of each of the modes of the scalar fieldφ(x) in an expanding spacetime background. Each of these modes were then taken to be initially in their lowest energy state which is a superposition of the growing and decaying modes familiar in inflationary studies focusing on the transition of quantum mechanical modes to classical degrees of freedom. The backreaction effects presented in this paper must therefore be carefully interpreted in a per mode basis.
The application of the results of the paper are physically restricted to the domain of the coupled Einstein- (massless) Klein-Gordon system (Eqs. (48) and (49)). The early Universe for all we know might be filled with more than just the inflaton and it is of interest to extend this work later to study the backreaction of spin-1/2 fields and massive vector fields to the spacetime. Nonetheless, the Einstein-(massless) Klein-Gordon system may also be able to draw information on the backreaction of massless higher-spin fields, e.g., each of the two polarization modes of gravitational waves generated during the inflationary era satisfies u + (k 2 − a /a)u = 0, which is the same equation as the massless scalar field (Eq. (65)). The case of a massive scalar field, with a bare mass m, could also be very simply accommodated in the model by transforming the wavenumber k as k → √ k 2 + m 2 . Through this simple transformation it can be realized that the existence of a short wavelength threshold translate to the existence of a heavy mass cutoff that would destabilize the inflationary era.
There is one more feature that dissociates this work from the realistic scenario but one which we heavily clinged on to for great simplification. We have assumed from Eq. (63) onwards that the quantum field lives only in one of its Fourier modes defined by Eq. (64). However, the more realistic picture will always be one which involves a nontrivial superposition of Fourier modes, say for instance, a Gaussian distribusion centered at k = K with a width σ K . To proceed in this way demands an extension of the CQC to a system with interacting QHOs each of which is also coupled to the same spacetime background. Alternatively, the construction of a spatial lattice could be established as in the field theory-extension of the CQC [40][41][42]. Either way is more computationally exhausting but is irrefutably a relevant future work.
A potential application of the CQC is on the preheating era where the quantum nature of the modes play a significant role [34]. The CQC can be used to improve the existing numerical recipes (see for instance Refs. [47,48]) by taking into consideration the quantum nature of the matter fields and could be used to study nonperturbatively the backreaction of quantum fields to both the inflaton and the metric. The backreaction can also be used to gain insight on the primordial power spectrum. The modes of relevance to this are those which start very deep inside the horizon but then eventually cross the horizon and become a constant. During the very early inflationary era, the perturbations of the metric can be shown to be subdominant when compared to that of the inflaton [43], and so the model used in this paper can be considered to be an approximate working model. Accounting for backreaction, the condition for horizon crossing of a comoving k-mode will therefore change from the naive k/a(η * ) = h to k/a(η * ) = H eff (η * ) where H eff (< h) is the renormalized Hubble parameter and η * is the conformal time at horizon crossing. The two-point correlation function of a quantum field,φ, which gives rise to the primordial power spectrum, will then become proportional to H eff (η * ) 2 rather than to h 2 . The detailed analysis of this is left for future work.

Conclusions
We have applied the CQC, for the first time, to a system described by a quantum field that is driven by a cosmological background (Eqs. (48) and (49)) in order to draw insigths on quantum mechanical backreaction during inflation. The results have shown that there exists a minimum Planckian wavelength that distinguishes between long wavelength-modes that would only renormalize the Hubble parameter and short wavelengthmodes that would destabilize the inflationary era. We have also obtained an approximate analytical solution for the short-wavelength modes which support the existence of such a threshold wavelength. The results of the CQC and the standard perturbative-iterative method was also compared. This lead to a strong support to the assertion that the CQC is the inevitable limit of the perturbative-iterative method.
The model used in this paper can be understood to be a working theory of the very early Universe when the modes are deep inside the Hubble horizon and the (b) δρ ψ,N =10 /ρ ψ = | (ρ ψ,N =10 − ρ ψ,CQC ) /ρ ψ,CQC | Fig. 8 Relative error between the tenth-order (N = 10) iterative-perturbative method and CQC for the scale factor, a, and the scalar field's energy density, 0|ρ (ψ) |0 for k = 8 with a 0 = 1 and h = 1. metric perturbations can be discarded compared to the inflaton perturbation. Nonetheless, it can be improved in many ways to be a more faithful representation of the early Universe. As a start, it would be desirable to extend this work to involve more than just a single mode backreacting to the metric. Also, we have relied so far on a fluid that is vacuum energy-like, to initially support the inflationary era. However, it is known that primordial inflation is more convincingly supported by a quantized slowly-rolling scalar field (a.k.a., the inflaton) in a potential that coexists with the spacetime and perhaps with some other fields in the early Universe. It is interesting to check out whether the addition of other fields could usher the inflationary era gracefully to the radiation era. The work must therefore be later on extended to viable inflationary theories (see for instance Refs. [49,50]) and compute the resulting primordial power spectrum and adiabatic invariants when the modes exit the horizon. Lastly, backreaction is a nontrivial effect even outside the context of inflation. The application of the CQC to simpler physical systems will therefore continue to be interesting and encouraged in future works.