Stochastic quantum inflation for a canonical scalar field with linear self-interaction potential

We apply Starobinsky's formalism of stochastic inflation to the case of a minimally coupled scalar field with linear self-interaction potential. We solve the corresponding Fokker-Planck equation exactly, and obtain analytical expressions for the stochastic expectation values.


Introduction
Inflation [1,2] is widely accepted as the standard paradigm of the early Universe. The first reason is due to the fact that several long-standing puzzles of the Hot Big-Bang model, such as the horizon, flatness, and monopole problems, find a natural explanation in the framework of inflationary Universe. In addition, and perhaps the most intriguing feature of inflation, is that it gives us a causal interpretation of the origin of the Cosmic Microwave Background (CMB) temperature anisotropies [3], while at the same time it provides us with a mechanism to explain the Large-Scale Structure (LSS) of the Universe, since quantum fluctuations during the inflationary era may give rise to the primordial density perturbations [4].
Despite its success a theory of inflation is still missing, since as of today we do not know neither what the inflaton is nor why the inflaton potential is so flat. Furthermore, the inflaton potential cannot be derived from a fundamental theory in a unique way. All we have is a large collection of inflationary models (see e.g. [5], and for a classification of inflationary models see e.g. [6]). In a e-mail: grigorios.panotopoulos@tecnico.ulisboa.pt single-field inflationary models with a canonical scalar field in Einstein's General Relativity (GR) there are basically four classes of models, namely large field models, small field models, hybrid models and linear inflation that lies between large and small field models. After the latest Planck results [7] many inflationary models based on monomial potentials have been ruled out or are disfavored by data. There are of course the Starobinsky's model [1] as well as the Higgs inflation model [8] which are in agreement with the data, but they are extensions of GR. The linear inflaton potential is still in agreement with observations, but difficult to realize in particle physics models. However, recently it was shown that the Coleman-Weinberg potential [9] together with a nonminimal coupling of the inflaton to gravity leads to attractor solutions that interpolate between quadratic inflation, which is ruled out by recent measurements, and linear inflation, which lies within the allowed region [10,11]. In particular, first in [10] it was shown that the predictions of linear inflation for the observables can be achieved in the context of well-defined quantum field theory, without introducing complicated interactions by hand. Then in [11] the authors extended the previous discussion, presenting a more detailed study of the parameter space, and they also added a discussion on reheating.
During inflation infrared logarithms arise in the expectation values of operators of quantum field theories that contain massless minimally coupled scalar fields or gravitons. For an incomplete list see e.g.
[12] and references therein. These terms are powers of the logarithm of the inflationary scale factor a = exp(Ht), with H being the Hubble constant, and are very exciting because they may signal important quantum effects [13]. However, their continued growth implies that when inflation has proceeded for a long time, the large 2 logarithms eventually overcome the small coupling constants. Therefore the effects become non-perturbative and perturbation theory breaks down. A natural approach to obtain non-perturbative information is the leading logarithm approximation [14], in which one attempts to sum the series comprised of just the leading infrared logarithms at each order. For scalar fields with non-derivative interactions Starobinsky's technique of stochastic inflation [15,16] recovers the leading infrared logarithms at each order, and the series of these leading effects at all orders can be resummed [17] to give non-perturbative predictions.
Given the observational bound on the tensor-toscalar ratio r < 0.07 [18], the slow-roll parameter ǫ = r/16 is extremely close to zero, ǫ < 0.0044, and therefore the de Sitter spacetime is an excellent approximation to inflation. The goal of this work is to study the effect of a minimally coupled scalar field with a linear self-interaction potential using Starobinsky's formalism. Our work is organized as follows: After this introduction, we briefly summarize stochastic processes in the Langevin and Fokker-Planck approach in the second section. Then we apply this formalism to the case of a canonical scalar field with a linear potential following the Starobinsky's technique in section three. Finally we conclude in the last section.

Stochastic processes
The observation that small pollen grains, when suspended in water, are found to be in a very animated and irregular state of motion, was first systematically investigated by the Scottish botanist Robert Brown in 1827, and the observed phenomenon naturally was named after him. First A. Einstein [19] in 1905 and independently M. Smoluchowski in 1906 [20], and some time later Paul Langevin [21] in 1908 (see [22] for a translation of the original Langevin paper in English) explained Brownian motion using different but equally successful approaches. From the one hand, Einstein's analysis was based on the diffusion equation The solution of the diffusion equation is given by [23] f (t, and the mean of the square of displacement is given by x 2 = 2Dt. On the other hand, Langevin started from Newton's equation of motion assuming a Stokes's drag force and a random thermal force due to the continuous bombardment from the molecules of the liquid. Although he did not exploit all richness of his model, Langevin obtained in the late-time regime Einstein's result, namely where k B is Boltzmann's constant, T is the temperature, µ is the fluid viscosity, and a is the radius of the particle. Therefore, thanks to Langevin's approach the Brownian motion, diffusion as well as the random walk were all linked together, a view that soon quantified experimentally by Perrin [24]. Therefore the diffusion coefficient can be computed in terms of properties of the fluid and the Brownian particles, which is the Einstein-Stokes for- In modern times Langevin's approach is still in use. The Langevin equation for the process x(t) (let us call it the position of a moving particle) readṡ where A(x) is an external applied force, the dot denotes derivative with respect to time, and ξ(t) is assumed to be a Gaussian white noise Since the random force ξ(t) is not known, we can only compute mean values of powers of the position, or the moments, x n , once the density probability function is known. The density probability function u(t, x) satisfies the corresponding Fokker-Planck (FP) equation [25] ∂u(t, x) where the first term is due to the external force while the second term is the diffusion term with a constant diffusion coefficient D. If we ignore the external applied force, the FP equation reduces to the standard diffusion equation. That explains why the Einstein's and Langevin's approaches were equally successful. Solving the FP equation we then can compute the moments performing the integrals and they are functions of time. Although seemingly quite different equations, when the diffusion equation and the FP equation at hand have the same number of symmetries, or equivalently when the following condition is satisfied [26] 2DA the FP equation can be recast into the diffusion equation, and therefore admits an exact analytical solution. 3 3 Application of Starobinsky's formalism to a massless canonical scalar field with a linear potential In this section we apply Starobinsky's formalism [15,16] to a scalar field Ψ described by the Lagrangian where V (Ψ ) = −M 3 Ψ is the self-interaction potential of the scalar field taken to be linear in Ψ , with M being a mass scale. The resulting Klein-Gordon equation reads where V ,Ψ is the derivative of the scalar potential with respect to the scalar field. We take the onset of inflation to be at t = 0, and we work perturbatively around the non-dynamical de Sitter background with the inflationary scale factor being a(t) = exp(Ht), and H being the Hubble constant.
In the leading logarithm approximation the scalar field behaves like a stochastic variable φ satisfying Langevin's equatioṅ where the stochastic source f (t, x) has the properties of the Gaussian white noise (see however eq. (11) of [16] for the general case where the two space-like points x 1 , x 2 are taken to be different). Langevin's equation above can be rigorously derived from the underlying quantum field theory [15], or it can be easily obtained from the Klein-Gordon equation by applying the stochastic rules of [27], namely -During inflation the scale factor varies faster than the scalar field, and therefore the Hubble friction term dominates over the spatial derivative and the second time derivative terms.
-We replace the full field Ψ with its stochastic counterpart φ.
-We make the substitutionφ →φ − f , where f is the stochastic source.
Since the scalar potential is linear, the corresponding Fokker-Planck equation takes the form where we have put ∆ = H 3 /(4π 2 ). It is trivial to check that the condition (8) is satisfied, and therefore for this model we can find exact analytical solution of the FP equation. According to the previous discussion, the FK equation can be recast in the diffusion equation ω τ = ω yy and the solution is given by where h(t, φ), y(t, φ), τ (t, φ) are given by [26] h(t, φ) = exp Therefore, the final expression for the solution is given by Therefore it is now straightforward to compute the first stochastic expectation values φ , φ 2 , φ 3 and φ 4 , which are found to be where we have made use of the Gaussian integrals In the formulas obtained above, it is easy to check that when M = 0 we recover the known result for the pure diffusion case. Note that there are two contributions with different time dependence. In particular, the M dependent terms are due to the classical linear potential, while the M independent terms are due to quantum effects, and they are just the ones of massless scalar fields. Equivalently, we switch from the cosmological time t to the inflationary scale factor a by using t = ln(a)/H and we find the following expressions The factors ln(a) associated to the M independent quantum effects are the IR logarithms mentioned in the Introduction. The model is characterized by two mass scales, H, M , and we can view the ratio M/H as a dimensionless coupling constant. As already mentioned, even if M/H is small, when inflation has proceeded for a long time the large logarithms eventually overcome the small coupling constant.
Powers of H t due to quantum effects of massless gravitons slow inflation [28]. In this model, however, it is the classical effect of the scalar field rolling down the linear potential that slows inflation. As it has been shown in [29,30], massless canonical scalar fields (with vanishing potential) contribute negligibly to the energy density in de Sitter, and lead to negligible backreaction. On the contrary, the classical homogeneous field rolling down the linear potential discussed here will eventually develop an energy density comparable to the cosmological constant.
Since the probability density function is Gaussian, it is characterized by only two parameters, and therefore all moments can be given in terms of the first two. However, in the following we shall compute the generic stochastic expectation value φ 2n for even powers and φ 2n+1 for odd powers in closed form. To this end we change variable z = φ − (M 3 t)/(3H), use the binomial expansion and make use of the integral We obtain the final result φ 2n for even powers, and similarly for odd powers we obtain the formula φ 2n+1 where we have defined α −1 = 2∆t andφ = (M 3 t)/(3H), and it is straightforward to check that for n = 1, 2 one obtains the previous expressions for φ 2 , φ 3 and φ 4 . The last two expressions are the main result of this article. As most of the equations in realistic models are solved either numerically or approximately, it is always desirable to have exact analytical solutions. We find it remarkable that the linear scalar potential i) can be derived in the framework of well-established quantum field theory, ii) is still in agreement with observations, and iii) the corresponding Fokker-Planck equation can be recast into the diffusion equation and thus be solved exactly.
A final remark is in order here. Throughout this work we imagine that the scalar field is a spectator to inflation, and we have worked perturbatively around a non-dynamical de Sitter background, which is an excellent approximation to inflation. However, the scalar field can no longer be treated as a spectator when its energy density becomes comparable to the cosmological constant. This happens for a number of e-folds N * = H t given by Before finishing we show that it is possible to obtain the same results following another approach by solving the Langevin's equation directly. With the initial condition φ(0) = 0, the solution for φ(t) is given by From this one can immediately see that φ = (M 3 t)/(3H) since f (t) = 0. Then, by squaring the solution and using the property of the stochastic source, namely f (t 1 )f (t 2 ) = ∆δ(t 1 − t 2 ), one obtains the previous expression for the second moment. In a similar way the third and the fourth moments can also be computed.

Conclusions
In the present article we have applied Starobinsky's technique of stochastic inflation to the case of a minimally coupled scalar field with a linear self-interaction potential. This type of inflaton potential, although still in agreement with the latest Planck results, could not be obtained from a fundamental theory of particle physics. However, it recently became relevant since it has been shown that the linear potential can be obtained in the context of well-defined quantum field theory from a Coleman-Weinberg potential, provided that a nonminimal coupling to gravity is also present. We have obtained analytical expressions for the stochastic expectation values φ 2n (even powers) and φ 2n+1 (odd powers) in two ways. First by solving exactly the corresponding Fokker-Planck equation for the density probability function, and then directly by using the Klein-Gordon-Langevin equation and the properties of the stochastic source. The two approaches give us the same results as expected.

Acknowledgments
The author thanks the Fundação para a Ciência e Tecnologia (FCT), Portugal, for the financial support to the Multidisciplinary Center for Astrophysics (CEN-TRA), Instituto Superior Técnico, Universidade de Lisboa, through the Grant No. UID/FIS/00099/2013. In addition, he thanks the anonymous reviewer for her/his comments and suggestions that improved the quality of the manuscript, and R. Woodard for useful correspondence.