Stochastic inflaton wave equation from an expanding environment

We discuss the inflaton ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document} in an environment of scalar fields χn\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi _{n}$$\end{document} on flat and curved manifolds. We average over the environmental fields χn\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi _{n}$$\end{document}. We study a contribution of superhorizon k≪aH\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k\ll aH$$\end{document} as well as subhorizon k≫aH\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k \gg aH$$\end{document} modes χn(k)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi _{n}(\mathbf{k})$$\end{document}. As a result we obtain a stochastic wave equation with a friction and noise. We show that in the subhorizon regime in field theory a finite number of fields is sufficient to produce a friction and diffusion owing to the infinite number of degrees of freedom corresponding to different k\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{k}$$\end{document} in χn(k)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi _{n}(\mathbf{k})$$\end{document}. We investigate the slow roll and the Markovian approximations to the stochastic wave equation. A determination of the metric from the stochastic Einstein–Klein–Gordon equations is briefly discussed.


Introduction
Inflation is generated by an interaction of gravity with a scalar field (inflaton) [1][2][3][4][5]. In a purely Hamiltonian system it is impossible to reach a thermodynamic equilibrium exhibited in the CMB radiation. The equilibration can be achieved by an interaction with an environment of other fields (heat bath). Then, the temperature of the inflaton subsystem is decreasing during inflation. In the standard (cold inflation) approach a mechanism of reheating [6,7] is applied to raise the temperature. The presence of the environment limits the decrease of temperature [8][9][10]. Berera [11] described the effect of the environment by a stochastic modification of the inflaton equation. He derived a stochastic equation from a linear interaction with an infinite set of scalar fields in a heat bath. The model follows the well-known derivation of the Brownian motion from dynamical systems [12,13]. We have obtained the same stochastic equation in the low momentum regime (k 0) by different mathematical methods in [14]. In [11,14] the expansion of the environment has been neglected in the derivation (as we set a −2 k 2 0). What was essential for the a e-mail: zbigniew.haba@uwr.edu.pl (corresponding author) result was an infinite set of fields with masses proportional to couplings. Such a relation allows to apply the Markov approximation. In Appendix B of Ref. [11] Berera derives the same diffusion equation if there is only one field in the environment but k a H. In this paper we investigate models of an interaction with the environment in more detail starting from models in the Minkowski space. We consider some generalizations of the models of [11,14]. The main new result concerns the calculation of the effect of the subhorizon modes k a H on the wave equation of the inflaton. We show that if the momentum and scale-dependent term in the scalar-field equations is dominating (in the subhorizon regime k a H) then we obtain a different diffusive inflaton equation than the one resulting from averaging over superhorizon modes (in disagreement with [11]). The appearance of the diffusive behaviour in this model is similar to the one discussed by Starobinsky and Vilenkin [15][16][17][18] when modes with large k of the quantum scalar field lead to a diffusive behaviour of the remaining superhorizon modes k a H in an expanding universe. In such a case those authors were able to derive a stochastic wave equation describing quantum fluctuations in an arbitrary inflaton potential. Quantum as well as thermal fluctuations determine the CMB spectrum which can be compared with observations [19,20].
The plan of this paper is the following. In Sect. 2 we review a minor generalization of the model of Refs. [11,14]. In Sect. 3 we discuss in detail the model in Minkowski spacetime. In Sect. 4 we study the environmental fields for subhorizon momenta in de Sitter space. In Sect. 5 we explore an expansion in a fixed homogeneous metric which is close to exponential. We eliminate the environmental fields, average over the initial values and obtain a stochastic wave equation for the inflaton with a friction and noise. We briefly discuss the slow roll conditions resulting from our modification of the stochastic equation for warm inflation. We show that in some approximations the resulting wave equation for the inflaton can be treated as a Markovian stochastic wave equation. In Sect. 6 we briefly discuss the Friedman equation which determines the homogeneous metric ( the Hubble variable as a function of the inflaton field). We believe that the averaging over the initial values simulates some quantum effects in cosmological models as the quantum field theory at finite temperature tends to the classical field theory with random initial conditions distributed according to the classical Gibbs law. The thermal noise in the inflaton equation can play a similar role as the scalar quantum fluctuations (discussed in [4]) which together with the quantum fluctuations of the gravitational field during de Sitter expansion (calculated earlier in [21]) lead to the correct evaluation of the power spectrum (this can be done by means of the methods developed in [22,23]).

Scalar fields interacting linearly with an environment
We consider the Lagrangian which is an extension of the well-known oscillator model discussed in [12,13] [11,14]) but the results of averaging over χ n do not depend essentially on U . The number of particles in the classical mechanics of [12,13] is infinite. Only in the limit of an infinite number of degrees of freedom the irreversible diffusive behaviour can appear. In field theory even with a finite number of fields we have an infinite number of degrees of freedom (the spatial Fourier modes). Averaging over those modes can lead to a diffusive dynamics. Equations of motion read where g μν is the metric tensor and g = | det[g μν ]|. We can consider classical as well as quantum systems (2)-(3). Eliminating the quantized fields χ we obtain a quantum version of the environmental noise. The quantum noise can approximate the quantum fluctuations of scalar and gravitational fields as it does in the e-fold time [22,23].
In the flat expanding metric Equation (3) reads where H = a −1 ∂ t a. We can solve Eq. (3) for χ n where we denote U (x) = U (φ(x)), G n is the Green function of the operator M n , A n is any operator commuting with M n and χ cl n are solutions of the homogeneous equation M n χ cl n = 0. When we insert χ n of Eq. (6) in Eq. (2) then it takes the form where In the homogeneous metric (4) we can take the spatial Fourier transform of Eq. (3). Then, the operators A n become just functions A n (k). In classical field theory (with the energymomentum tensor T μν ) in Minkowski space T 00 (χ )dx is a constant of motion as long as there is no coupling to φ. Then, the Gibbs distribution at temperature β −1 is dχ exp(−β T 00 (χ )dx). A natural covariant generalization to a manifold is the weight factor for classical configurations of the form where the integral d 4 x is over a volume β in space-time.
There is an analog in quantum statistical mechanics to the formula (9) which is applied as a statistical operator in quantum non-equilibrium statistical mechanics. For a homogeneous metric (4) the volume β in Eq. (9) can be expressed as a time interval [0, β] times the space volume. After averaging over the solutions A n χ cl n (t, k) the correlation function of the noise is In this way we obtain a stochastic wave equation (7)with a friction δφ and the noise η. There is some arbitrariness in the choice of the potentials U and V , the number of fields χ , the masses and the couplings. In [11,14] U (φ) = φ, an infinite set of fields is chosen and the masses are proportional to the couplings. As shown in [14] if we ignore the k dependence of M n (superhorizon domain) then we obtain a wave equation of warm inflation [8,9] with the friction proportional to ∂ t φ and the noise η as the white noise. In the Appendix B of Ref. [11] one field χ with large momenta is discussed. It is claimed that its effect is the same as an infinite number of χ fields leading to an alternative derivation of the same stochastic equation for inflaton. We disagree with this claim. For low momenta of the environmental fields another stochastic equation (studied in [30]) appears than for large momenta. We study these equations in detail in this paper beginning with the simple case of the Minkowski space in the next section.
The Hamiltonian for the wave fields χ n is ( = ∂ τ χ ) Then, the expectation value with respect to the Gibbs density (9) is In quantum theory A * n (k) and A m (k ) become creation and annihilation operators. Then, the expectation value is When we calculate the expectation value of η with respect to the measure d A exp(−βH) then we obtain In quantum theory we should consider an anticommutator on the lhs of Eq. (16). Then, β −1 ω −2 n → ω −1 n 1 2h coth(h 1 2 βω n ) on the rhs. So that the quantum thermal expectation values tend to classical expectation values (in the Gibbs state) when hβ → 0.
Let us distinguish two regimes: small k and large k. In the first case we neglect k in ω n . Then Let (2π) 3 2 λ n =γ m n then If m n are uniformly distributed then (because of oscillations of the cosine) the sum (18) is concentrated at small t − t .
We could e.g. apply the Poisson summation formula (for m n n) or approximate the sum in Eq. (18) by an integral over m n leading to the δ function. In both cases (for a large L or in the approximation nγ 2 · · · = γ 2 . . .) we obtain with a certain γ γ . Next let us consider the friction δφ in Eqs. (7) and (13). Writing sin(ω n (t − s))ω −1 n = ∂ s (cos(ω n (t − s))ω −2 n ) and integrating by parts we get the term which under the assumptions leading to Eq. (19)can be considered as modification of the potential V → V + qU 2 with q δ(0) (this is the mass renormalization if U φ as in the Leggett model [24]). The second term from the integration by parts under the assumptions used at Eq. (19) is proportional to δ(t) and to the initial values of the field. We will neglect this term for t > 0. After the integration by parts the friction term takes the form We have got the same kernel as the one in Eq. (16) which we approximated by δ(t −s) in Eq. (19) for the correlation of the noise. In fact these kernels are related by a version of the fluctuation dissipation theorem [25]. With these approximations Eq. (2) in Minkowski space reads with the noise (19). This is the equation which would have been derived in [11] and [14] when a = 1. Note that in the limit of a strong friction we obtain Hence, U (φ) behaves like a Brownian motion in classical field theory and in the high temperature limit of the quantum field theory. If V is of the same order as γ 2 U ∂ t U (this may be the case because of the U 2 renormalization mentioned below Eq. (19)) then Eq. (20) after the neglect of the second order derivatives reads U has a large time asymptotic distribution where dU V U (U ) means that we first express V U −1 as a function of U and subsequently calculate its indefinite integral. So, U = 0 may change the form of the asymptotic behaviour of the inflaton modifying the discussion in [26].
We can derive the slow roll conditions leading to the approximation (21). They require The slow roll approximation can be applied to quantum fields as well. We can see that with the friction term the quantum field behaves like a diffusion process. In a linear model when V = m 2 2 φ 2 and U = 1 then the slow roll requirement is m 2 γ 4 which does not depend on φ. The correlation functions in this model are the same as the ones for the Ornstein-Uhlenbeck process [27]. The model (21) can be considered as a limit H γ 2 of the model in de Sitter space discussed in the next section (in this sense ∂ t U brings the friction in the Minkowski space as the Hubble constant does in the de Sitter space, in both cases the strong friction approximation is applicable).
Next, assume that k is large in comparison to m n . Then, neglecting the masses in Eq. (16) we have where In general, from the integral ( [28], formula 3.721) we have Note that the covariance (25) is vanishing for the time-like separations Let us still discuss various forms of the friction term δφ in Eq. (7). If we do not integrate by parts then the friction in Eq. (7) has an explicitly Lorentz covariant form is the commutator function for the massless quantum scalar field. The friction term in Eq. (13) can also be expressed in a Lorentz invariant way as Now, Eq. (7) takes the form The kernel in Eq. (13) with m n = 0 after an integration by parts is This kernel is the same as the covariance (25) of the noise and is vanishing for time-like separations. The stochastic equation after an integration by parts reads with the friction kernel K and the Gaussian noise (25). That the variance of the noise and the kernel of the friction are related follows from a version of the fluctuation dissipation theorem [25]. Let us note that the result (27)-(28) is exact if all m n = 0 and if the initial values of the χ n fields are distributed according to the Gibbs weight factor exp(−βH n )(the number of fields is irrelevant). If then we can approximate the kernel K by a δ-function leading to the stochastic equation where η in this approximation has the correlation (23). The condition (29) can be interpreted as a negligence of the spatial dependence of φ. If the spatial dependence of φ is ignored then we could identify Eqs. (20) and (30) although we have derived them on a basis of different assumptions concerning the χ fields (this could explain the "alternative" derivation of the stochastic equation in Appendix B of [11]). In the strong friction limit we obtain In Eq. (30) we can take V into account and integrate both sides over x. Then, we could conclude that dyU has an equilibrium distribution. In contradistinction to the strong friction limit of Eq. (20) only the space average tends to the Brownian motion in Eq. (30) (and eventually to an equilibrium). However, there is a distinction between the results (20) and (30). Equation (20) is derived under an assumption of an infinite number of χ fields with properly chosen masses, whereas Eq. (30) follows from Eq. (27) which is exact if m n = 0. In conclusion, in the strong friction limit the classical thermal field as well as the thermal quantum fields behave as a stochastic process. As an example, if in Eq. (30) we neglect second order time derivatives and assume V = U 2 then Hence, dyU (t, y) is the Ornstein-Uhlenbeck stochastic process.

Exponentially expanding environment
In an expanding universe we write Then, in the momentum space where First, let us consider low momenta a −1 k H so that for a large time (as in [14]) we can neglect a −2 k 2 term. Then, we assume that ω 2 n > 0, H is slowly varying and ω n are approximately constant. The approximate solutions of Eq. (32) at λ = 0 are (the same as in Eq. (11)) They are normalized as In de Sitter space H = const. If we neglect a −2 k 2 and assume m 2 n > 9 where the last two terms on the lhs of Eq. (34) result from the integration by parts as in Eq. (20) with an extra term coming from the transformation (31).The noise has the covariance At slow roll approximation with the neglect of second order derivatives we obtain with a = exp(Ht). The slow roll requirements are |∂ 2 t φ| |(3H +γ 2 (U ) 2 )∂ t φ| and |∂ 2 t φ| |V + 3 2 γ 2 HU U |. These conditions are satisfied if As an example, if H = const, U = 1 and V = m 2 φ then the slow roll condition is m 2 (3H + γ 2 ) 2 . It does not depend on φ and allows to solve the slow roll stochastic equation explicitly.If H is determined by Einstein equations then H depends on φ (as will be discussed in Sect. 6). Then, the slow roll conditions still must involve some requirements allowing to determine H as a function of φ.
The solution of Eq. (36)in the limit of large γ is Hence, we get an extra damping factor in comparison to the case of the Minkowski space (discussed at the end of Sect. 3). It is instructive to consider the soluble case of Eq. (36) with U = 1, V = m 2 φ. Then, the solution (with zero initial condition) is the Ornstein-Uhlenbeck process We have got the same result as the one for the de Sitter space with the Hubble constant H → H + γ 2 3 or in Minkowski space with a friction (and H = 0).
Another (complementary) regime is the one when the term a −2 k 2 is dominating (a −1 k H and a −1 k m n ). In such a case ω n a −1 k. Hence, in the solution (33) of the wave equation (32) where τ is the conformal time. In the conformal time (for the exponential expansion a(t) = exp(Ht)) Equation (5) reads and Eq. (2) The free Lagrangian of the χ fields is The statistical weight (9) with the Hamiltonian H = T 00 (with the momentum = τ 2 (∂ τ χ )) is For a large k we may ignore the masses m n . Then, the solution of the free wave equation is where H (1) 3 2 is the Hankel function [28].
We can write the general solution in the form If we insert the solution (43) in the Hamiltonian (41) and neglect the terms (decaying fast for large τ ;large τ means small a, hence close to the Big Bang) and τ −2 (∂ τ χ) 2 ,then The τ -dependence cancels for a large τ in Eq. (45) as expressed in Eq. (41) because τ −2 in √ gH cancels with τ 2 in (∇χ) 2 resulting from the terms τ exp(ikτ ) in Eq. (43). We can calculate the expectation value of the noise with respect to the weight factor (41). We obtain We can see that the last term of the sum in Eq. (47) is divergent at small k. This is a typical infrared divergence of a massless χ field theory in de Sitter space. It means that instead of the random variable η we should consider η(x) − η(z) with a fixed z = (s, z). Then, the variance Using the formula (24) we calculate the Fourier transform in Eq. (47) (neglecting the terms decaying fast for a large time) with the result where χ cl is the solution with λ n = 0 and Using the formulas 3.741 of [28] ∞ 0 du sin(zu) and the integral (24) we can calculate the Fourier integrals in order to express the kernel G in configuration space Note that where (x −x ) 2 is the Minkowski distance. Hence, the leading term (for a large time) is Lorentz invariant. After an insertion of the solution of Eq. (38) in Eq. (39) the stochastic equation (7) for φ in the configuration space reads For |x − y| |τ − τ | and large τ and τ we get the stochastic equation by an expansion of the kernel (51) in |x − y| and integration by parts over τ in Eq. (52) (as in an analogous derivation of Eq. (30)). The limit of small |x − y| of the noise follows from Eq. (48) In the cosmic time Eq. (53) reads With the neglect of ∂ 2 t φ (slow-roll) and φ ( which disappears after a space average in Eq. (53)) we obtain In the limit of the strong friction (only κ 2 terms on the lhs) The solution is an Ornstein-Uhlenbeck process We can see that the classical field in a thermal environment tends to a diffusion process. We could prove the same behaviour for the quantum field in a thermal state. It seems that in both cases the appearance of fluctuations and friction is crucial for the diffusive behaviour. This phenomenon is similar to the diffusive behaviour of zero temperature quantum fields in de Sitter space which also results from "friction" H ∂ t φ and from quantum fluctuations [15,18,29].

Environment in an almost exponential expansion
We consider a flat expanding metric which is close to de Sitter. For k a H we may repeat the calculations of Sect. 4 and Ref. [14] leading to an analog of the stochastic equation (34) (in the cosmic time) with the noise (19). It can be shown that without the friction terms Eq. (56) is invariant with respect to the general change of coordinates. With the friction terms even the Lorentz invariance is violated because in the derivation of Eq. (56) from Eqs. (2)-(3) we have neglected the spatial derivatives of χ fields. Next, we explore χ fields with large k a H in a conformal time (we follow to some extent our earlier paper [31]) where the conformal time τ is

With a slowly varying H we have approximately
(t is the cosmic time). If the expansion is close to exponential then with the relation (57) Eq. (3) is where 3η n = m 2 n H −2 Let Then We obtain another form of Eq. (3) if we write Then where With λ n = 0 the solution of Eq. (65) is the Hankel function [28]˜ n = H (1) ν n (kτ ). For a general expanding metric the solution of a homogeneous Eq. (59) (λ n = 0) is In general, we have a superposition of classical solutions with different k The noise is The solution of Eq. (62) with λ n = 0 on the rhs is ψ ν = τ −α+μ H ν . Then, the solution of the wave equation (62) for n with U (φ) on the rhs is where G is the Green function of the operator on the lhs of Eq. (62). The Green function can be constructed from the two independent solutions of the homogeneous equation (62) where from Eqs. (61) and (64) The Bessel functions J and Y [28] are defined by the Hankel function H where w(k) is the wronskian. If a is close to exponential then w(k) = 2k π . The Green function could also be expressed as We insert the Green function (74) in Eq. (69) in order to calculate δφ of Eq. (7) for a large time (we use the asymptotic expansion of H ν [28]). Then From Eq. (75) there follows the stochastic equation (7) ( In Eq. (77) we have a non-Markovian and non-local friction term. For a large time G can be approximated by the first term of Eq. (51). If we assume that only small |x − y| contribute and expand in this variable then As a result for a small |x − y| we obtain a stochastic equation generalizing Eq. (53) Inserting in Eqs.
(2)-(3) the Green function (92) and the wave functions exp(ik t a −1 ) we obtain in the limit |τ − τ | |x − x | the stochastic equation where the noise η has the covariance With Eq. (93) we can repeat the argument, which has been based on Eqs. (56), (83)-(89), in order to express H and a in terms of φ. The study of the Einstein-Klein-Gordon system without the assumption |τ − τ | |x − x | is more involved and requires a separate investigation because in this case H ,a and the noise depend on φ in a non-local way.

Summary and outlook
We considered a system of non-linear wave equations describing an interaction of the inflaton φ with other scalar fields χ n . In quantum mechanics if the fields χ n are not observable then we average over the states of χ n . As a result the inflaton is described by a density matrix resulting from the averaging. The natural state for averaging is the Gibbs (thermal) state of maximal entropy for the χ n fields. Then, further evolution does not depend on the initial time when the interaction between the inflaton and the χ n fields is switched in. We discussed the model in the classical limit when the probability distribution is described by the canonical Hamiltonian of the χ n fields. In an expanding metric the canonical Hamiltonian depends on time. However, for large conformal time in models with an expansion close to the exponential the time-dependence disappears. The average over the environmental fields is reduced to an average over the initial values of the χ n fields in the Gibbs state. There is some arbitrariness in the choice of the environment; the number of fields χ n , the couplings λ n , the masses m n and the interaction λ n χ n U (φ) (we choose an interaction linear in the environmental fields). We discuss evolution of the inflaton field in an external expanding metric in two regimes. First, we choose the increasing masses proportional to couplings and assume negligible (superhorizon) momenta k a H of an infinite set of the environmental fields. These assumptions after an averaging over χ n lead directly to a Markovian wave equation with a friction and white noise. In such a case the quantum evolution is described by a master equation of the Lindblad type. Another (subhorizon) limit discussed in this paper involves large momenta k a H and k m n . The number of fields χ n can be finite. After an averaging we obtain a wave equation for the inflaton which is non-local and non-Markovian in the friction and noise. If we restrict ourselves to correlations φ(τ 1 , x 1 ) · · · φ(τ n , x n ) such that |τ j − τ k | |x j − x k | then the stochastic wave equation has a Markovian limit as a diffusion wave equation with friction. This equation is different from the wave equation derived for low momenta of the environmental fields. Both equations may coincide in a linearized version after a spatial averaging. In the last section we discussed a determination of the expansion scale factor a from the Einstein-Klein-Gordon system of equations. In the superhorizon case it is possible to obtain a stochastic wave equation with a(φ) and H (φ) dependent on φ and the noise independent of φ. In such a case we can calculate all expectation values on the basis of the Fokker-Planck equation. An interesting problem for further studies is the dependence of slow roll conditions and spectral indices on the interaction U (φ). In the subhorizon case the noise and a(φ) depend in a non-local way on φ. Only in the limit of small spatial distances (|x j − x k | |τ j − τ k |) we are able to derive a Markovian stochastic wave equation with a(φ) and H (φ) determined by Einstein equations.
In this paper we concentrate our attention on the stochastic equations satisfied (in various regimes) by the inflaton interacting with an environment. The environment consists of superhorizon or subhorizon modes of the classical thermal environmental fields. In [15,16] the subhorizon modes of the quantum field have been treated as the environment. The environmental modes are approximated by a noise. If the noise is approximately Markovian then the quantum states satisfy the irreversible Lindblad equation. It is well-known that Lindblad equation leads to decoherence (meaning a destruction of interference of quantum states). The interference is a substantial obstacle to the classical limit. In the standard model the large scale structure is formed from quantum fluctuations [4,21] in the early universe. Hence, in the meantime quantum fluctuations should become classical. In [32,33] it is shown that if initially the universe starts from the Gaussian state then with the time evolution determined by a linear perturbation theory the resulting squeezed state becomes classical. The measured CMB fluctuations are indistinguishable from classical fluctuations. At present, all observations [19,20] show no departure from Gaussian correlations. However, it may be that prospective CMB measurements can reveal some non-Gaussian behavior. Such CMB results could show a destruc-tion of interference of states in the early universe and the role of the environment in their decoherence as discussed in [34]. The noise in Einstein-Klein-Gordon system plays at least two roles:it leads to decoherence and it determines the power spectrum (as calculated in [23] with the quantum noise taking into account quantum gravitational fluctuations). The detailed analysis of the contribution of the noise from the environment to the power spectrum and to higher order correlation functions of the inflanton could give important information about the quantum state in an early universe and its classical limit.

Data Availability Statement
This manuscript has associated data in a data repository. [Authors' comment: All data included in this manuscript are available upon request by contacting with the corresponding author (Zbigniew Haba).] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .