Quantum randomness in the Sky

In this article, we study quantum randomness of stochastic cosmological particle production scenario using quantum corrected higher order Fokker Planck equation. Using the one to one correspondence between particle production in presence of scatterers and electron transport in conduction wire with impurities we compute the quantum corrections of Fokker Planck Equation at different orders. Finally, we estimate Gaussian and non-Gaussian statistical moments to verify our result derived to explain stochastic particle production probability distribution profile.

It is a well known fact that the particle production scenario in the early universe cosmology (during reheating) follows the dynamical master equation, aka Klein-Gordon equation. On the other hand, transport phenomena of electron through a conduction wire with impurities follow time independent Schrödinger equation. Both of this dynamical time dependent phenomena have structural one to one correspondence [1,2]. Anderson Localization and saturation of the chaos are some well studied phenomena in the context of scattering problem can be extended to describe the quantum randomness phenomena during cosmological particle production. From their inherent stochastic nature quantum chaos can be related to them and chaos bound can be defined either by Lyapunov exponent [3] or by Spectral Form Factor [4,5]. The possible quantum effects arising from higher order corrections in dynamical master equation aka Fokker Planck equation for particle production scenario in the early universe cosmology (during reheating) can be achieved from the present discussion. For comparing scattering event with stochastic particle production Dirac Delta profile of time dependent coupling (mass function) is chosen, m 2 (τ ) = N j=1 m j δ D (τ − τ j ), localized at time scale τ = τ j (where j represents the number of non-adiabatic events). Further using the concept of transfer matrices occupation number can be computed from this set up. To model a phenomenological situation where width (w j ) of the profile of the time dependent coupling is finite and the scattering event is relevant, we consider sech scatterers. It is important to note that, in the limit w j → ∞ the Dirac Delta profile can be recovered from this phenomenological profile.
In the context of disippative system, Fokker Planck equation explains the probability density for particle position of Brownian motion in a random system. For a Markovian process this situation can be expressed by Chapman-Kolmogorov equation [1]. Now considering Maximum Entropy Anstaz we can derive the Fokker Planck equation from Smoluchowski equation when we integrate the probability density over the angular coor-dinate θ: P (n, θ, φ; τ + δτ ) ≡ P (n, θ; τ + δτ ) → P (n + δn; τ ) δτ (1) where we consider an infinitesimal change (δθ) is not functionally dependent on θ. Further Taylor expansion of P (n + δn; τ ) δτ with respect to the infinitesimal occupation number (δn) with the constraint in this context P (n; τ ) δτ = P (n, τ ) gives the following result: (2) This gives the following general structure of Fokker-Planck equation which we will use for our all calculations: Using Smoluchowski equation the occupation number can be expressed as: δn ≡ n 2 (1 + 2n) − 2 (1 + n 2 )(1 + n)n 2 n cos 2(φ 2 − θ) (4) which help us to further define various statistical moments from the probability density function. Assuming that the particle production rate is small locally (µδτ < 1) we have the truncation in Taylor expansion. With primary truncation in first order (δn) δτ Fokker-Planck equation is derived as: Here the mean particle production rate have Fourier mode dependence (µ k ). By Fourier transformation with respect to the occupation number n of the distribution function: P (n; τ ) = (2π) −1 dk e iknP (k; τ ).
Which simplifies the Fokker Planck equation in Fourier space : Imposing initial condition for probability distribution function at time τ is given by the Dirac Delta profile or its derivatives in different orders we get: where J denoting the order of quantum corrected Fokker Planck Equation.
For J = 1 we get the following solution of the probability density function: Comparing the coefficient of δτ from the both sides of the Taylor expansion we get quantum corrected Fokker Planck equation at different order. Without truncation on both sides of this expression additional contributions in δτ and in its higher order can be obtained and generate quantum corrected version of the Fokker Planck equation valid upto higher orders. All such higher order corrections justify non-Gaussian effects appearing during cosmological stochastic particle production in reheating phase. In another words origin of higher order contributions describe the quantum effects from its non vanishing statistical moments originating from quantum correlations.
FIG. 1. Variation of fourth order quantum corrected probability distribution with respect to the particle number density.
From Fig. (1) the P i (i=1,2,3,4) denote the i-th order probability distribution. The order by order small corrections (fluctuations) from Gaussian profile support the quantum effects in stochastic particle production. From the quantum corrected probability distribution we can further calculate different statistical moments using Eq (3). Calculating expression for n , n 2 , n 3 and n 4 and standard deviation, skewness and kurtosis for a given time solidify the quantum nature as predicted earlier.
To compute the first moment of the occupation number we use the first order master evolution equation: To compute the second moment we use first and second order master equations in two different orders: 1st order : µ −1 k ∂ τ n 2 = 4 n + 6 n 2 , (17) 2nd order : µ −2 k ∂ 2 τ n 2 = 12 n + 12 n 2 + 2. (18) Continuing in the same way one can similarly calculate third and fourth moments corrected upto different orders. From Fig. (2(a)) we obtain the large variance with increasing τ . But the quantum corrected and uncorrected distribution have same variance at all time signifying that width of the peak is unchanged by the quantum effects.
Additionally, it is important to note that the computed probability distribution function has a long right tail a specific effect of positive skewness. Considering different order correction in n 3 and standard deviation we calculate skewness with and without correction. Now from Fig. (2(b)), we can say that the corrected Skewness deviate significantly from the uncorrected one at low τ limit. But we can see that at higher time scale they overlap. So for particle production at initial time the skewness is dominant over uncorrected skewness. So the effects of quantum corrections are more clearly visible for initial time scale. Using the corrected n 4 and standard deviation we calculate the kurtosis for particle production event which we have shown in Fig. (2(c)). Here we have shown the quantum corrections are dominant at large time scale, but at low time scale both the corrected and uncorrected kurtosis overlap with each other. Itô : ∂ τ P (n; τ ) = ∂ 2 n (D(n)P (n; τ )), (19) Stratonovitch : ∂ τ P (n; τ ) = ∂ n ((D(n)) 1−Q ∂ n ((D(n)) Q P (n; τ ))) (20) where D(n) = n(n + 1). Using this we get the following solution of probability distribution: Stratonovitch : P (n, τ ) = 1 2 µ k πτ (n(n + 1)) Q × exp − n 2 (n + 1) + µ k τ (2n + 1)Q(Q + 1)(n(n + 1)) Q 2 4µ k τ (n(n + 1)) Q+2 .
The probability distribution function obtained from this have the same the log normal form. From General perspective the Fokker Planck equation with effect of potential(U (n)) can be expressed as: ∂ n (D(n)∂ n (exp(βV (n)/2)P (n; τ )) −U (n)(exp(βV (n)/2)P (n; τ )) = ∂ τ (exp(βV (n)/2)P (n; τ )), where Choosing D(n) = n(n + 1) and V (n) = n 2 we get the following probability distribution function: Fig. (4) of the probability distribution function we observe that for large value of occupation number the distribution function decays to a finite saturation value. On the other hand for small occupation number we get peak in the distribution function for different values of µ k τ .
With different order solutions of Fokker Planck equation we construct probability density function which explain the quantum nature in stochastic particle production scenario in early universe cosmology. Also the existance of higher order statistical moments of the probability density function. The present approach can extend to explain the semi-classical behaviour of particle production event and relating chaos to this approach eventually set a bound to the quantum randomness [4].
From this non-gaussianity in stochastic particle production during inflation period we can connect it with our idea of non-gaussianity in finite universe. Considering interacting background field it will be possible to introduce the other non-linear and dissipative effects into the system introduced by the background itself and can be studied as open quantum system interacting with the defined background set-up [6].Using more general statistical field theory { Itôo,Stratnovitch,General background potential }this results can be genralized for any system with randomness within it.