Stochastic dark energy from inflationary quantum fluctuations

We study the quantum backreaction from inflationary fluctuations of a very light, non-minimally coupled spectator scalar and show that it is a viable candidate for dark energy. The problem is solved by suitably adapting the formalism of stochastic inflation. This allows us to self-consistently account for the backreaction on the background expansion rate of the Universe where its effects are large. This framework is equivalent to that of semiclassical gravity in which matter vacuum fluctuations are included at the one loop level, but purely quantum gravitational fluctuations are neglected. Our results show that dark energy in our model can be characterized by a distinct effective equation of state parameter (as a function of redshift) which allows for testing of the model at the level of the background.


Introduction
The origin of dark energy (DE) is one of the most fascinating unsolved problems of modern science. In literature traditionally two main classes of solutions have been proposed [1][2][3][4][5][6][7][8][9]: • matter condensates or physical DE, of which the simplest representatives are scalar condensates (quintessence models); • modified gravity or geometrical DE, which mimics dark energy by changing the relation between geometry and matter or by supplying additional geometric fields to general relativity.
However, there is no impenetrable barrier between these two possibilities, and more generically DE can be both physical a e-mail: drazen.glavan@fuw.edu.pl b e-mail: t.prokopec@uu.nl c e-mail: alstar@landau.ac.ru and geometrical, i.e. a new matter field has to be introduced and gravity becomes modified, too. This just happens in the case of DE described by a non-minimally coupled scalar field considered in this paper. Recently the effective field theory (EFT) approach to dark energy [2,5,10] has been developed. Its beauty is in that it presents a unified framework for both approaches, but its drawback is in that it does not immediately select the fundamental theory that lies behind some EFT. Nevertheless, different theories can be mapped onto the same class of EFTs, such that one can think of EFTs as identifying universality classes associated with DE models. The question of naturalness of initial conditions is not addressed in traditional approaches. For example, in quintessence models typically a quintessence field starts running from a value which is not a (local or global) minimum of the potential. Criticisms are often brushed away by noting that similar malady plagues most of inflationary models. Arguably the main benefit of this work is in that we construct a theory that naturally explains the initial field value -which is accounted for by the calculable amplitude of infrared field fluctuations during inflation -thus addressing this fundamental criticism. The program advocated here can be thought of as a third way for understanding DE, in that in our class of models a link is established between primordial inflation and dark energy. This link, among other things, can be exploited when designing tests for these models.
Observers have devoted a lot of effort (and observational time) to nail down as accurately as possible the amount (and distribution) of dark energy. Since its discovery in 1998 [11,12] a lot of progress has been made in improving the accuracy of DE measurements [13][14][15]. At this moment the Planck satellite [16,17] and the Dark Energy Survey (DES) [19] collaborations provide the most stringent bounds on dark energy. Presently ΛCDM, which assumes a cosmological constant equation of state p = − ρ for dark energy, is consistent with all astronomical data. Assuming a general constant equation of state parameter p = wρ (the so-called wCDM model), combined Planck and Type Ia supernovae data provides constraints w =− 1.006 ± 0.045 (68% CL).
Considering dynamical dark energy models -where w that vary with redshift z -yields significantly relaxed constraints [17]. Still no statistically significant deviation of DE from an exact cosmological constant has been found. 1 The upcoming measurements of Large Synoptic Survey Telescope (LSST) [20], ESA's Euclid satellite mission [21] and the European Extreme Large Telescope (E-ELT) [22] will further constrain dynamical DE models. For example, the accuracy of Euclid [21] is projected to be a few percent for simple (1 or 2 parameter models) and weaker for more involved models. Furthermore, these missions will be able to test some prominent DE models, including clustering of DE, growth of Universe's structure, interactions of DE with itself and with other cosmological fluids, and -last but not least -the class of models presented here.
However, if the present DE traces its origin to a very early stage of the Universe, initial conditions for its subsequent evolution were quantum. As a result, it can remain quantum even up to the present time. Specifically in this work we investigate the influence of the quantum backreaction of a very light, nonminimally coupled spectator scalar field on the expansion dynamics of the Universe at late times. Our interest is to investigate what is the effect of inflationary quantum fluctuations that survive until late times and become comparable to -and eventually dominate overthe background nonrelativistic matter energy density driving the matter-dominated expansion. This idea goes back a while ago to [23,24] where it was noted that quantum fluctuations of a very light scalar could provide a CC-like contribution to the Friedmann equations. In recent years, similar ideas were examined [25][26][27][28][29][30]. Most recently in [30] it was shown that, for certain ranges of model parameters, the backreaction remains small throughout the expansion of the Universe, and becomes relevant only at late times in matter era where it behaves approximately like a cosmological constant. Three parameters were introduced in this model: the total number of e-foldings of inflation N I , the nonminimal coupling ξ , and the scalar field mass m. The conditions these parameters need to meet for the scenario to unfold are [30] (1) 1 Recent search for different channels of radioactive-type decay of DE with time-independent decay rates also resulted in upper limits on these rates less than the inverse present age of the Universe only [18].
The first condition guarantees that the scalar field stays light throughout the expansion, the second one comes from requiring the backreaction to remain small (perturbative) during inflation, and the third one ensures that the leading term in the backreaction at very late times is of the cosmological constant type. Here H I is the inflationary Hubble rate, H 0 the Hubble rate today, and H DE the Hubble rate at the onset of dark energy domination. The magnitude of the cosmological constant today implies another relation that determines the total number of e-folding of inflation N I in terms of the remaining two parameters, It turns out this limits the range of the nonminimal coupling to 0 > ξ ≥ −10 −2 , and the number of total e-foldings of inflation is N I 10 3 .
In the regime where the backreaction becomes comparable to the background it can no longer be treated perturbatively, but rather its effects have to be taken into account properly by solving the semiclassical Friedmann equations self-consistently, with the quantum backreaction as a source. In general, these are rather complicated integro-differential equations. In principle these equations can be solved numerically by discretizing both the time evolution and the integral over the modes for the backreaction. An example of such computation was done in the context of reheating in [31], but it still presents a challenging numerical problem. Instead of the numerical approach we make use of the observation from the analytic computations [26,28,30] that the quantum backreaction for very light fields is dominated by the infrared (IR) modes not only during inflation, but throughout the history, and that the spectrum of these modes is inherited from the inflation era. The evolution of IR modes (which are largely amplified) during inflation (for small enough nonminimal coupling) is very accurately described by the formalism of stochastic inflation fully formulated in [33] (though some of its simplest applications were already used in [32]). Here we adapt this method so that it becomes applicable to the evolution of IR modes in subsequent radiation and matter dominated eras of the Universe, and reproduce all the perturbative results of [30]. Such an approach to the problem was already advocated long ago in [3], just after the original inception of the idea [23,24], but -to our knowledge -it has not been carried out. We then use this method to approximate the semiclassical Friedmann equations in the regime where quantum backreaction is very large and write them as a closed set of differential equations which are much simpler to solve than the original integro-differential equations of semiclassical gravity. This approximation captures accurately the backreaction effects. These equations are then solved numerically. We find that the backreaction indeed accelerates the Universe, driving it towards a late de Sitter phase, very much like the observed behaviour of the Universe today, representing thus a novel dark energy candidate. This paper is organized as follows. Section 2 presents the definition of the scalar field model and the standard cosmological history. In Sect. 3 the stochastic formalism is derived for the model at hand, and in Sect. 4 the stochastic equations are solved in the regime of small backreaction and the results of [30] are reproduced. Section 5 presents the numerical solution of the full semiclassical Friedmann equations at late times, when backreaction is no longer small. In Sect. 6 we summarize and discuss the principal results.

Scalar field in FLRW
In this section we define the cosmological background with the standard expansion history, on which our spectator scalar field model is defined and quantized.

FLRW background
The metric of the D-dimensional spatially flat Friedman-Lemaître-Robertson-Walker (FLRW) spacetime is given by where a is the scale factor, and the derivative of its logarithmic with respect to time defines the Hubble expansion rate, H =ȧ/a. The conventions that we use for the geometric quantities are Γ α μν = 1 2 g αβ (∂ μ g νβ + ∂ ν g μβ − ∂ β g μν ) for the Christoffel symbol, R α μβν = ∂ β Γ α μν −∂ ν Γ α βμ +Γ ρ μν Γ α ρβ − Γ ρ μβ Γ α ρν for the Riemann tensor, and R μν = R α μαν for the Ricci tensor. A convenient quantity we often use instead of time is the number of e-foldings N , which for the period between some t 0 and t is equal to N = ln(a/a 0 ).
The classical evolution of a spatially homogeneous, spatially flat universe is governed by the Friedmann equations, where M P = (8π G N ) −1/2 ≈ 2.45×10 18 GeV is the reduced Planck mass, 2 G N is the Newton's constant, and ρ c = ρ c (t) and p c = p c (t) are the classical pressure and energy density of the (dominant) cosmological fluid, respectively. Here we are primarily interested in solving for the Universe's dynamics which also includes the quantum backreaction from inflationary quantum fluctuations of matter fields at the one- where ≈ 0 during inflation (where it is known as the slowroll parameter), = 2 during radiation-dominated period, and = 3/2 during matter dominated period (see Fig. 1 for a schematic depiction of the evolution of parameter). During = const. periods the scale factor and the Hubble rate evolve as where a 0 and H 0 refer to the values at some particular time t 0 , conveniently taken to be the beginning of a given constant period. The evolution of the Hubble rate is depicted schematically in Fig. 1.

Nonminimally coupled massive scalar
The action for the nonminimally coupled scalar field in a curved background is where m is the scalar field mass, R is the Ricci scalar, and ξ is the nonminimal coupling with the sign convention − −− → 1/6 corresponds to the conformal coupling. We consider this model on an FLRW background.
In order to quantize the scalar field we first need to define the canonical formalism. The canonical momentum is and accordingly the Hamiltonian, Canonical quantization is now accomplished by promoting the field and its conjugate momentum to operators, and their Poisson brackets to equal-time commutation relations, The Heisenberg equations of motion for the field operators are now, where the effective mass-squared is It is convenient to expand the field operators in Fourier modes,Φ where the expansion of the conjugate momentum operator is such that it automatically satisfies Eq. (10). Observations tell us that the Universe is nearly spatially homogeneous on large scales. This motivates a spatially homogeneous Ansatz for the mode function ϕ(t, k), where k denotes the modulus of the comoving momentum k = k , and the annihilation and creation operators,b andb † , satisfy canonical commutation relations This puts a constraint on the mode function in the form of a Wronskian normalization, The equation of motion for the mode function inferred from (10-11) is The Hilbert space of states is constructed in the standard manner by defining a vacuum state |Ω , which is annihilated by all the annihilation operators,b(k)|Ω = 0, ∀k. Then the rest of the state space is constructed by creation operators acting on |Ω . The energy-momentum tensor operator of the scalar field is given bŷ Its expectation value · with respect to a homogeneous and isotropic state |Ω defined on a FLRW background (4) is diagonal and takes the form of ideal fluid, where ρ Q = ρ Q (t) and p Q = p Q (t) are spatially homogeneous.

Stochastic formalism
Here we briefly introduce the stochastic formalism of [33] for scalar fields in expanding cosmological space. The dominant contribution to the light scalar field correlators in inflation comes from the superhorizon modes (of wavelengths k > 1/a H). This is also true for the very light or massless scalar fields in subsequent evolution of the Universe, namely during radiation and matter periods [28,30]. The slow-roll approximation usually employed to derive the stochastic equations need not be correct in situations when the scalar becomes very massive, meaning m H . Since during decelerating periods of expansion the Hubble rate decays, this can eventually become true. This is why we do not use the slow-roll approximation in the stochastic formalism, but rather derive the equations for all three IR correlators. These reduce to the standard equations when the slow-roll hierarchy between the correlators is present.

Equations of motion for field operators
We split the field operators into the contributions from the long wavelength modes and from the short wavelength modes as follows, separated by a comoving scale μa H, where 0 < μ 1 is the control parameter of the splitting. This splitting is most conveniently imposed in Fourier space where the long wavelength parts are and the short wavelength partŝ where, for simplicity, we took the Heaviside θ -function for the window functions. We want to derive the analog of equations of motion (10) and (11), but just for the long wavelength parts of the field operators (22) and (23). Making use of the equation of motion (17) for the mode function the resulting equations are where the sources -which originate from the coupling between the short and long wavelength fields -arê The two sources (28) and (29) are to be considered as stochastic forces acting on the long wavelength fields, and here they originate from the effect of modes leaving (entering) the Hubble sphere due to accelerating (decelerating) expansion.

Equations of motion for IR correlators
Here we derive the equations for motion for the coincident long wavelength (infrared, IR) two-point functions, which we define to be conveniently rescaled, The convenience of time dependent factors in the definitions above is in that all of these correlators become of the same dimension and the way they appear in the equations allows us to compare them directly which makes comparing their magnitudes straightforward. From now on we set D = 4 since the UV divergences are captured by the short wavelength part of the fields. The equations of motion for these correlators follow from the equations of motion (26) and (27) for the long wavelength parts of the field operators (in D = 4), In deriving these equations we have thrown away gradients, in particular terms of the form (∇φ) 2 , since they are suppressed by a factor of μ 2 1 compared to the rest. This allows one to close the set of equations. Also we have switched to the number of e-foldings N = ln(a) as the time variable (with the choice of time (gauge) that corresponds to a = 1 at the beginning of inflation), and = −Ḣ /H 2 . The stochastic sources on the right hand side of (33-35) are the coincident field-noise correlators, These are straightforwardly computed from (22,23) and (28,29) to be f φ (t, x),π(t, x) f π (t, x),π(t, x) where the anticommutator is defined as {Â,B} ≡ÂB +BÂ, and where we have used that θ(0) = 1/2. 3 The sources (36)(37)(38) are referred to as the stochastic sources in analogy with the stochastic inflation formalism, though their physical interpretation during decelerating cosmological eras is not necessarily the one of standard stochastic sources. During decelerating eras ( > 1) these sources have a negative overall sign, as evident from (39)(40)(41)(42), and in fact contribute towards decreasing the correlators. This is not at all surprising if one notes that these sources arise as a time-dependent split of the modes into super-Hubble and sub-Hubble ones. During decelerated eras modes leak from the super-Hubble to sub-Hubble phase space (as opposed to accelerating periods where the leak is reversed) which accounts for the decrease in the IR correlators. The energy density and pressure expectation values can be expressed in terms of the coincident IR correlators, The gradients were thrown away here on the same grounds as in the equations of motion (33)(34)(35) of the previous subsection. The short wavelength contributions were neglected sinceafter renormalization -except possibly during inflation, they contribute negligibly. 4 These expressions for energy density and pressure can now be used to quantum-correct the Friedmann equations (in the sense that one should exact the replacements, ρ c → ρ c + ρ Q and p c → p c + p Q in (4)) and together with equations of motion (33-35) form a closed set of stochastic differential equations. It turns out that at late times (during radiation and matter era) stochastic sources contribute negligibly which will greatly simplify solving these equations.

Comparison with field theoretic calculations
Here we employ the stochastic formalism of the preceding section to calculate the backreaction of the scalar field quantum fluctuations during the three relevant cosmological periods and solve them perturbatively. These solutions are applicable in the regime where the quantum backreaction (43,44) is still negligible in comparison with the classical sources ρ c and p c in (4). The sources for the equations of motion (33)(34)(35) are determined using the mode functions from [30]. This section serves to reaffirm the results of [30] as well as to test the (perturbative) correctness of the stochastic formalism.

De Sitter inflationary period
The results presented in this subsection were already derived in [34,35], both for test fields and for cosmological perturbations, and on the more general slow roll inflationary background.
The Chernikov-Tagirov-Bunch-Davies (CTBD) mode function of the scalar during the de Sitter inflation is where H (1) ν stands for the Hankel function of the first kind, and H I = const. is the inflationary Hubble rate. This mode function coincides with the one of the adiabatic vacuum state for the subhorizon modes. The leading order IR expansion of the mode function (45) (for k <μa H I a H I ) is With this we can determine stochastic sources (36)(37)(38) for the equations of motion, which are just constant. To leading order in nonminimal coupling and mass, 0 < |ξ |, m 2 /H 2 I 1, they are independent of the separation scale μ, where the shorthand notation is X = 12ξ + (m/H I ) 2 , and it is clear that |X | 1, which can be used as an expansion parameter, encompassing both minimally coupled and massive case. Parameter X plays the role of an effective mass squared (in units of H 2 I ) of the field during inflation. It can take negative values X < 0, signaling tachyonic particle production. This is not troubling, since the tachyonic nature of X depends on the background spacetime, and is not present e.g. in flat space or in radiation era. 5 The equations of motion for the coincident correlators (33)(34)(35) are now, The natural initial conditions for these equations are These can be seen as arising from the set of modes which were initially super-Hubble (k < μH I ), meaning that initially the IR phase space is zero (this point we identify with the beginning of inflation at which a = 1). This is not strictly speaking so, we rather assume that the super-Hubble modes at the beginning of inflation are suppressed due to some physical mechanism (for examples of how this can be done in practice see [36]). In practice it is sufficient to assume that Δ φφ < H 2 I at the beginning of inflation. This results in them contributing subdominantly at late times. Now the leading order solution in |X | 1 for these IR correlators is The subdominant contributions are suppressed by additional factors of X 3 compared to the leading order one, 6 and also decay much faster with the number of e-foldings. Given the initial conditions (52), the stochastic sources determine the amplitude of the IR correlators. These sources also induce the hierarchy between the correlators which survives until very late times. There are two interesting limits we discuss briefly in the two following subsections. 5 The tachyonic nature of the IR modes during inflation is not problematic, since one can use e.g. pre-inflationary radiation era to regulate the infrared sector of the theory by a smooth matching of inflationary modes onto the radiation era modes [36]. 6 The exact form of the exponent in solutions (53) is −3 1 − √ 1−4X/9 N , but we have expanded it under the assumption |X | 2 N I 1, where N I is the total number of e-folding. This will be true in our case.

Minimally coupled limit
The minimally coupled limit corresponds to taking X = m 2 /H 2 I . There are two different cases to consider, the first one of which is the "long" inflation limit N I (H I /m) 2 , where the correlators eventually saturate to constant (de Sitter invariant) values, The other case is the limit of "short" inflation In both cases the hierarchy between the correlators is dictated by the small mass.
In the case of "long" inflation the energy density and pressure by the end of inflationary period are and in the case of "short" inflation they are Note that the contribution from the conformal anomaly ρ C A = −p C A = −119H 4 I /(960π 2 ) is included in both expressions since it is not negligible in these cases. In subsequent periods though, its contribution will be negligible. It should be noted that for our purposes the "long" inflation limit is extremely long (about 10 13 efoldings), and the result (56) is computed under the assumption of exact de Sitter inflation (exactly flat plateau inflation). In cases of small deviations from such assumptions on inflation the scalar spectator does not necessarily reach stationary values for its correlators as pointed out in [44], and one needs to be more careful when relating the correlator values of the spectator to the inflationary period history. What is of primary interest for us in this work is to establish that the scalar can attain super-Planckian correlator values by the end of inflation. The details of precisely how it attains these values will introduce model dependent corrections to the relation (2) relating the number of e-foldings of inflation to the properties of the late time Universe.

Massless limit
In the limit of negative nonminimal coupling which dominates over the mass term, 1 |ξ | (m/H I ) 2 , to leading order the field is effectively massless, and the correlators grow exponentially with the total number of e-foldings of inflation, This growth of correlators during inflation is important for us for the model building of DE. Note that in this case the leading order behavior of the correlators at late times in inflation is independent of the stochastic source, namely their growth is dominated by the instability of IR modes during inflation. Nevertheless, the stochastic sources, which in this case contribute significantly only during the beginning of inflationary period, fix the amplitude and the hierarchy of the correlators, accounting for the fact that all the contributing modes are of UV origin.
The energy density and pressure of the backreaction by the end of inflationary period in this case are

Radiation period
In radiation period we have R = 2, so a 2 H = a 2 1 H I , where a 1 refers to the value of the scale factor at the end of inflation/beginning of radiation. From [30] we have for the IR limit (k < μa H a H) of the modulus squared of the scalar mode function, where the expression is valid up to the first subleading order in small mass. The other combinations of the mode function appearing in the stochastic sources can be determined from taking derivatives of the expression above, and by using the equation of motion (17). In particular, where the gradient terms were dropped. The stochastic sources during radiation period, to leading order in 0 ≤ |ξ |, (m/H ) 2  The order (m/H ) 2 to which the mode function was computed in [30] does not allow us to compute the leading order contribution to the last stochastic source above, but this is not important. Stochastic sources are relevant in inflation since they determine the amplitude of the correlators, but can actually safely be neglected in the subsequent periods of expansion. This argument can easily be made by examining the equations of motion, and by comparing the relative size of various terms. The final values of correlators at the end of inflation (55) or (58) serve as the initial conditions for the radiation period, and we note that both in minimally coupled, or massless limit they contain enhancement factors of N I or e 8|ξ |N I /|ξ |, respectively. These derive partly from expanding the IR phase space of modes included in long wavelength correlators. On the other hand, none of the sources on the right hand side of (63)-(65) has this enhancement factor, making them negligible in comparison. Hence, to a very good approximation, the equations for IR correlators during (late) radiation era are the same as the classical ones, These equations we may solve perturbatively in small mass (and small nonminimal coupling contained in the initial conditions). To do this properly we need to take into account that the initial conditions of the radiation period satisfy the hierarchy inherited from inflation, which is dictated by the same perturbative parameters. Therefore the equations that correctly capture the leading order correlators (denoted by superscript (0)) are which we can readily solve, This result encompasses both minimally coupled and massless limits. At the end of radiation period (when N R ∼ 50) we can further simplify (73-75) by dropping the terms that have exponentially decayed, where N R ≈ 50 is the total number of e-foldings during radiation, 7 and H eq is the Hubble rate at the end of radiation and beginning of matter period (radiation-matter equality). The energy density and pressure to leading order in radiation period are which is exactly what was obtained in [30]. There are two competing contributions in (79, 80): a radiation-like fluid with negative energy and pressure, and a cosmological constant type contribution.

Matter period
For the matter period the correlator values at the end of radiation (76)-(78) serve as initial conditions. Applying the same reasoning as in the radiation period case, the leading order equations of motion for the IR correlators in matter period are d There are two competing contributions in (90), one behaving like nonrelativistic matter (with negative energy), and the other one like a cosmological constant (CC). If the parameters of the model are such that the CC-like term starts dominating over the matter term before late times in matter era it should also eventually dominate over the background fluid and accelerate the expansion. This regime is beyond the (perturbative) regime discussed in the present section, where the backreaction on the expansion rate was neglected. It requires solving the set of semiclassical Friedman equations (4) (with ρ c → ρ c + ρ Q and p c → p c + p Q ) together with the matter era correlator equations self-consistently, as described at the end of Sect. 3. The expansion dynamics during that regime are determined in the next section.

Self-consistent solution and results
Here we present solutions of the semiclassical Friedman equations with the initial conditions set during matter era at time t = t * when still H (t * ) = H * H DE . The quantum backreaction energy density and pressure (43,44) represent sources for the Friedman equations (4) on top of the classical matter fluid. The equations of motion for the correlators (33)(34)(35) with the sources neglected close the full set of equations.
More specifically, initial conditions are set at redshift z = z in = 10 (corresponding to about 2 e-foldings in the past from today, N (z) = − ln(1+z)) in such a way that the energy density and pressure of the would-be cosmological constant is substituted for the quantum backreaction. The equations are then solved numerically. The solutions for different choices of parameters are given in Figs. 2 and 3 which show how geometric quantities such as the Hubble rate H and = (3/2)(1+w) (vertical axis) depend on the number of e-foldings N = − ln(1+z) (horizontal axis).
In Fig. 2 we see that, when ξ = 0 and m = 0.20H 0 [m = 0.01H (z = 10)], the relative difference between the Hubble expansion rate and in our model and ΛCDM is at a sub-percent level, which is not observable by current observations and at the verge of being observable by the near future observatories. However the differences become more pronounced when the scalar mass approaches the Hubble rate today. Indeed, when m 0.61H 0 [m 2 = 10 −3 H 2 (z = 10)] the difference in the expansion rate and (or equivalently w) reaches several percent (see right panels in Fig. 2), which is testable by the next generation of observatories such as Euclid [21] and ELT [22].
In order to study the effect of the nonminimal coupling ξ on the geometrical quantities H (z) and (z), in Fig. 3 we show the evolution of the Hubble parameter and as a function of the number of e-folding N = − ln(1+z) for ξ = −10 −3 . As in Fig. 2 we see that for small masses, when m H 0 (m 0.20H 0 ) the difference between ΛCDM and Note that the uncertainty is larger at larger redshift. However, that does not reduces in any significant way testability of the model, because at large redshifts the relative contribution of dark energy is small our model is at a sub-percent level and hence barely reachable by planned observatories. However, for a mass comparable to the Hubble rate today (m 0.61H 0 ), the difference becomes significant and well within reach of Euclid and ELT. Comparing Figs. 2 and 3 reveals that increasing |ξ | seems to increase the difference between our model and ΛCDM.
However, for small values of ξ 's (ξ = −10 −3 in Fig. 3) this increase is rather moderate. Note that a non-vanishing mass term (or a potential term) is crucial for the model to exhibit DE behaviour. (The case m = 0 and ξ <0 is examined in [28], and that model does not exhibit a late time dark energy.) The nonminimal coupling is a relevant parameter here more from the point of view of relating the DE period to the primordial inflationary period, more precisely to the total number of e-foldings of inflation N I . In case of zero nonminimal coupling this number has to be huge [29], N I ∼ 10 12 for both cases of Fig. 2, while for even rather small nonminimal couplings of Fig. 3 the total number of e-foldings of inflation is much smaller [30], N I ∼ 10 5 , which can be determined from (2). Larger values (in magnitude) for nonminimal cou-pling are problematic in the sense that one has to worry about large backreaction much earlier, during primordial inflation. Although this is an interesting possibility to study in the context of graceful exit from inflation, it is not considered in our DE model (for a recent study of the backreaction of a heavy nonminimally coupled scalar in de Sitter inflation see [37]).
While Figs. 2 and 3 nicely show the recent Universe dynamics in our model and compare them to those of ΛCDM, it is instructive to show how the dark energy component in our model evolves, and that is illustrated in Fig. 4. The upper two plots show that for the vanishing nonminimal coupling (ξ = 0), the backreaction behaves like a cosmological constant to a very good approximation in the late Universe, causing the transition from the matter-dominated to the quasi-de Sitter phase. This behaviour tends to persist until the Hubble rate drops below the mass scale of the scalar field (m ∼ H ). After this point the Universe transitions to a decelerated phase where oscillates between 0 and 3.
The lower two plots in Fig. 4 demonstrate the effect of non-minimal coupling on the dark energy equation of state.
The shaded regions correspond to the uncertainty we have in estimating its equation of state. These are introduced for the following reason -due to nonminimal coupling the backreaction contains a contribution (whose energy density and pressure we denote as ρ m Q and p m Q ) that to a good approximation scales like nonrelativistic matter at higher redshifts (while the backreaction is still small), as is evident from Eq. (90). It is more appropriate not to consider this contribution a part of dark energy, but rather to group it together with the background matter fluid. Hence, we have subtracted the contribution from the backreaction that scales away exactly like mater, and determined the equation of state for the remaining contribution in the backreaction fluid. The uncertainty in w Q is due to the fact that the amplitude of this matter like contribution is computed approximately to order O(ξ ), and that it does not scale exactly like matter, but its behaviour has corrections of the order of O(|ξ |N M ), where N M is the number of e-foldings of the matter-dominated period. This means we do not know exactly how much of this matter-like fluid to subtract, which introduces sizeable uncertainties in the equation of state only for higher redshifts, but for smaller redshifts (z 3) we get -up to sub-percent uncertainties -unique predictions. This implies that our model of dark energy is in principle testable as it gives a definite prediction for the equation of state parameter w Q as a function of the redshift.
Even though we know the subtracted matter-like contribution redshifts like ordinary matter on cosmological scales, we do not know how it clusters. If its clustering turns out to be insignificant, its principal effect would be a finite renormalization (reduction) of the Newton constant. In the opposite case the effect would fall under the class of some interacting dark matter/dark energy models, the precise nature of which requires further investigation, which is left for future work.
As of yet dark matter has not been directly seen. However, in the near future dark matter may be seen in the laboratories and if that happens, we may have enough to infer its energy density. Since our model predicts a negative component that redshifts as dark matter (see Eq. (90)), if it also clusters like dark matter, there should be more dark matter than what indirect observations (such as lensing) indicate. More concretely, for [m/H (z in )] 2 = 10 −4 and ξ = −10 −3 the ratio is: ρ m Q /ρ c =− 0.0687 and for [m/H (z in )] 2 = 10 −3 , ξ = −10 −3 the ratio is ρ m Q /ρ c = − 0.0085, which is large enough to be potentially observable (both evaluated at redshift z 10). As far as we know, our model is the only one in the literature that exhibits this feature, which may become a smoking gun for testing our model.
Let us now compare with current observations. The Planck collaboration papers [16,17] [17]), which shows that the constraints w 0 and w a depend significantly on the data sets used and hence the quality of the data is still not good enough to place strong bounds on the model. Since these bounds are only slightly modified when the more recent 1 year DES data are accounted for [19], they suffice for our purpose. It is clear that these bounds are not precise enough to place meaningful constraints on the dark energy models shown in Fig. 4. This can be seen even clearer in Fig. 5, in which we observe that the contours of our dark energy model lie comfortably within one standard deviation contours of the w 0 w a CDM model (which meaningfully constrains dark energy only at relatively small redshifts, z 2, or N −1).
However, with the incipience of upcoming high precision observatories such as Euclid and LSST, the constraints may become strong enough to curb our model, especially when combined with inhomogeneous probes.

Discussion and outlook
This paper is one in a series [26,28,30], whose goal is to investigate viability of DE models based on quantum fluctuations during inflation. The attractiveness of these models stems from the fact that -when compared to more conventional DE models -they are amenable to additional tests. For example, since the origin of these models can be traced back to inflation, there is an intricate relationship between DE and inflationary observables. Furthermore, since the origin of DE in these models is in quantum fluctuations of matter fields, in general we expect these models to generate inhomogeneous dark energy with calculable amount of inhomogeneities, which can also be used to test this class of models.
In this paper we reanalyze the dark energy model proposed in Ref. [30]. The model utilizes an ultra-light, non-minimally coupled scalar field and its dynamics is governed by the action (6). Even though this field is a spectator field during inflation, quantum fluctuations naturally grow large such that, by the end of inflation, they typically reach super-Planckian values. In fact, due to the assumed negative non-minimal coupling, it can take as little as a thousand of e-foldings to achieve this feat. This mechanism addresses the question of naturalness of initial conditions -often raised in quintessence DE models -by dynamically generating super-Planckian fluctuations from initially subhorizon quantum fluctuations. In order to be able to calculate beyond the perturbative regime of Ref. [30], in Sect. 3 we develop a stochastic framework which permits us to solve the problem self-consistently, also in the regime where the quantum backreaction dominates the background dynamics. To establish the accuracy of our stochastic model, in Sect. 4 we use the stochastic framework to calculate the one loop contribution to the energy density and pressure of quantum fluctuations and demonstrate agreement with the perturbative calculations of Ref. [30] in all of the relevant epochs: inflation, radiation and matter era. This is sufficient to establish perturbative accuracy of the model. The main advantage of the stochastic formulation is in that a self-consistent solution of semiclassical gravity comes within reach of even modest programmers, as the original set of integro-differential equations of semiclassical gravity is reduced to a much simpler set of ordinary differential equations with rather facile stochastic sources.
A notable result of our work is that stochastic sources -which are of essential importance for capturing the correct dynamics during inflation -can be neglected during the post-inflationary radiation and matter epochs, such that the field dynamics becomes essentially classical. The memory of stochastic sources is still kept in the initial conditions for radiation era and for matter era.
Armed with this formalism we then study how quantum inflationary fluctuations in our model (6) evolve in time and how they affect the background evolution. In particular in Figs. 2 and 3 we compare the evolution of the Hubble rate H (t) and its dimensionless rate of change, = −Ḣ /H 2 (which is equal to the principal slow roll parameter during inflation) in our model with the same in ΛCDM. Our results show that, as the ultra-light mass m is taken to be closer to, but still remains smaller than, the Hubble parameter today, the differences become significant and observable by the near future observatories such as Euclid and ELT. From Figs. 4 and 5 we see that the current data are not constraining enough to meaningfully test our the backreaction DE models. However, the near future data (from Euclid, SST, ELT and other planned observatories) will penetrate deeper into the redshift space and will be collecting much more data and thus will be able to constrain the quantum backreaction DE models.
Next, we point out that the dependence of w Q = w Q (z) in our quantum backreaction model is quite different from that in typical quintessence models [6]. This owes to the rather unusual feature of the perturbative initial condition in matter era (79, 80), which is in turn dictated by the standard (IRregulated) Chernikov-Tagirov-Bunch-Davies initial condition in inflation (45). Namely, there is a significant negative energy component that initially scales as the dominant background (dark matter) component (cf. Eqs. (79, 80)). This feature is akin to negative energy scenarios studied in [38] and -once the properties of (particle-like) dark matter are mapped out -can be used to test this class of DE models.
While in this paper we only study the background evolution, keeping in mind the quantum inflationary origin of scalar field dark energy in our model, it would be of utmost importance to investigate the growth of structure in the visible and dark energy sector in our DE model, both at Gaussian and non-Gaussian levels. Namely, due to the non-minimal coupling the model is non-Gaussian and hence we expect it to leave distinct Gaussian and non-Gaussian footprints at late times on Universe's large scale structure. Determining these will require a suitable adaptation of the stochastic formalism for off-coincident correlation functions in inflation [39] or of further developments in the stochastic δ N formalism used in [40,41].
From the theoretical point of view, the model we presented here addresses the issue of naturalness of initial conditions, but there still remains the issue of fine-tuning in the form of requiring the scalar to have a tiny mass m ≤ H 0 . This issue is not specific to our model, but is a generic feature of all classical quintessence models. It would be of great interest to investigate whether a similar effect can be reproduced within different (interacting) models, via a mechanism that utilises spontaneous symmetry breaking [42] or 't Hooft's naturalness hypothesis [43].