Time in Quantum Cosmology of FRW f(R) Theories

Abstract

The time problem is a problem of canonical quantum gravity that has long been known about; it is related to the relativistic invariance and the consequent absence of an explicit time variable in the quantum equations. This fact complicates the interpretation of the wave function of the universe. Following proposals to assign the clock function to a scalar field, we look at the scalar degree of freedom contained in $f\left(R\right)$ theories. For this purpose we consider a quadratic $f\left(R\right)$ theory in an equivalent formulation with a scalar field, with a FRW metric, and consider its Wheeler-DeWitt equation. The wave function is obtained numerically and is consistent with the interpretation of the scalar field as time by means of a conditional probability, from which an effective time-dependent wave function follows. The evolution the scale factor is obtained by its mean value, and the quantum fluctuations are consistent with the Heisenberg relations and a classical universe today.

1. Introduction

Since its formulation, general relativity has been a successful theory, verified in many ways and at any scale. However, there are instances where it does not reproduce in a precise way the results of observations, in particular the origin of the universe and the early and present inflationary phases, as well as the present matter distribution. The other issues are quantum gravity, for which there are well known approaches, string theory, loop quantum gravity, and supergravity, which has been attracting attention recently. In all these approaches, the problem of time is present and there are diverse proposals for its solution, see [1,2,3]. In the proposal of Page and Wooters [4], it is argued that all the observables are stationary, but dynamics arises by the behavior of the observables relative to one which plays the role of clock.

One way to approach these issues is by $f\left(R\right)$ modified gravity theories, proposed by Starobisnky [5] as an effective action of gravity obtained by coupling it to quantum matter fields, which explains inflation and reheating, and recently it has being used to explain the effects of dark matter and dark energy [6,7], see also [8,9,10,11,12]. Even if these theories appear as effective theories, one appealing feature of them is that they are pure gravity theories, with additional structure as the action is higher order. In fact it is possible to give equivalent actions [13,14], with scalar degrees of freedom.

Quantum cosmology of $f\left(R\right)$ theory has been considered in [15], as a proposal for the origin of the universe from a tunneling from “nothing” to the de Sitter phase of the Starobinski model. In this work were also computed, in the WKB approximation, the subsequent curvature fluctuations and the duration of the inflationary phase. Quantum cosmology of $f\left(R\right)$ theories has been studied also in [16,17,18,19].

In this work we consider the FRW quantum cosmology of $f\left(R\right)=R+\alpha R^2$ theory in the form of a second order theory with a scalar field given by O’Hanlon [13,20]. In this approach, the scalar field is auxiliary and has an apparent dynamics which corresponds to its expression in terms of the Ricci scalar. As we do not implement conformal transformations, we do not expect that there are equivalence problems as in the cases discussed in [21].

We solve numerically the Wheeler-deWitt equation, and consider the scalar field as a clock. Considering the positivity of the scale factor a, for the hermiticity of the conjugate momentum of a, we take a vanishing wave function of the universe at $a=0$. For the numerical computation, we consider a compact domain in a and $\varphi$. The solution is consistent to the vanishing condition on all boundaries, hence it is expected to be normalizable. Moreover, it ensures that the conjugate momentum of a is hermitian [22]. For the time interpretation fo the scalar fields, we follow [22], where an effective wave function corresponding to a conditional probability is proposed. In the second section we make a short analysis of the classical solutions and in the third section we show the numerical solution of the WdW equation considering values of the parameter $\alpha$ for which classically the solutions are qualitatively different. In the last section we draw some conclusions.

2. Lagrangian Analysis

Let us consider the $f\left(R\right)$ model for gravity without matter, with action

$$A=\frac{1}{2{\kappa }^2}\int dt\sqrt{-g}f\left(R\right).$$

(1)

The variation of this action leads to the equations of motion [11] $F\left(R\right)R_{\mu \nu }-\frac{1}{2}f\left(R\right)g_{\mu \nu }-{\nabla }_{\mu }{\nabla }_{\nu }F\left(R\right)+g_{\mu \nu }\square F\left(R\right)=0$, where $F\left(R\right)=f^{\prime }\left(R\right)$, see e.g., [8,10]. For a FRW geometry $g_{\mu \nu }=\mathrm{diag}\left(-N^2,\frac{a^2}{1-k{r}^2},a^2{r}^2,a^2{r}^2{sin}^2\theta \right)$, the scalar curvature is $R=\frac{6}{a^2}\left(\frac{a\ddot{a}}{N^2}+\frac{{\dot{a}}^2}{N^2}-\frac{a\dot{a}\dot{N}}{N^3}+k\right)$. In reference [5], for the quadratic action

$$f\left(R\right)=R+\alpha R^2$$

(2)

and for $k=0$, $F\left(R\right)$ is regarded as a scalar degree of freedom, the “scalaron” $\phi$. In this case the equations of motion are second order, and the dynamics of the scale factor is simply given by $a=exp\int Hdt$.

Instead of it, here we adopt the O’Hanlon action [13]

$$A=\frac{1}{2{\kappa }^2}\int \sqrt{-g}\left[R+\varphi (\beta \varphi +R)\right],$$

(3)

where $\beta$ is a free parameter. This action resembles the action used in [15], where the definition of the scalar curvature is regarded as a constraint. A variation with respect to $\varphi$ gives $\varphi =-\frac{1}{2\beta }R$, which leads to $A=\frac{1}{2{\kappa }^2}\int \sqrt{-g}\left(R-1/4{\beta }^{-1}{R}^2\right)$, i.e., $\alpha =-1/4{\beta }^{-1}$. Thus, with the FRW metric, action (3) can be written as

$$\begin{array}{c}\hfill A=\frac{1}{\kappa }\int dt\left[3(1+\varphi )(Nak-N^{-1}a{\dot{a}}^2)-3N^{-1}{a}^2\dot{a}\dot{\varphi }+\frac{1}{2}\beta N{a}^3{\varphi }^2\right].\end{array}$$

From the equation of motion of $\varphi$ follows

$$\varphi =-\frac{3}{\beta a^2}\left(a\ddot{a}+{\dot{a}}^2+k\right),$$

(4)

as well as the higher order equation for a [11]

$$3k^2+{\dot{a}}^2\left(6a\ddot{a}-9{\dot{a}}^2-2\beta a^2-6k\right)+a^2\left(6\dot{a}\stackrel{⃛}{a}-3{\ddot{a}}^2-2k\beta \right)=0.$$

(5)

A numerical solution of the last equation with vanishing initial scale factor, at $t=0$, leads as shown in [20], to an inflationary stage with exit, see Figure 1. See [11] for an analytical solution.

Substituting this solution in the expression for $\varphi \left(t\right)$ (4), we get Figure 2.

Thus, as can be also deduced from the solution of $\varphi$ in terms of the scalar of curvature, the scalar $\varphi$ take positive values. From these graphics we see that this scalar could be taken as a clock inverting the time, omitting far times, as shown in Figure 3.

3. Quantum Cosmology

The canonical momenta of the action (3) are ${\pi }_a=-\left[\frac{12a\varphi \dot{a}}{N}+\frac{12a\dot{a}}{N}+\frac{6a^2\dot{\varphi }}{N}\right]$, ${\pi }_{\varphi }=-\frac{6a^2\dot{a}}{N}$, and ${\pi }_N=0$, and the Hamiltonian is $H=N{H}_0$, where $H_0=-6ka\varphi -6ka-\beta a^3{\varphi }^2+\frac{\varphi {\pi }_{\varphi }{}^2}{6a^3}+\frac{{\pi }_{\varphi }{}^2}{6a^3}-\frac{{\pi }_a{\pi }_{\varphi }}{6a^2}$. After choosing Weyl ordering for the ambiguous operator products, the Wheeler-DeWitt equation follows

$$\left[-{\hslash }^2\left(\varphi +1\right)\frac{{\partial }^2}{\partial {\varphi }^2}+a{\hslash }^2\frac{{\partial }^2}{\partial \varphi \partial a}-2{\hslash }^2\frac{\partial }{\partial \varphi }-6\beta a^6{\varphi }^2-36k{a}^4\left(\varphi +1\right)\right]\psi (a,\varphi )=0.$$

(6)

This equation does not admit analytic solutions, unless approximation methods are considered [15]. Further, the scale factor satisfies $a\ge 0$, and as shown in the previous section, the scalar field satisfies the same condition $\varphi \ge 0$. Thus, to ensure hermiticity of the momenta ${\pi }_a$ and ${\pi }_{\varphi }$, it is required that the wave function vanishes at $a=0$ [22] and $\varphi =0$. Furthermore, for $a\to \infty$ the wave function must vanish, but if we allow that the universe does not collapse at future, then we do not impose conditions at an upper value of $\varphi$.

An analysis in this case can be done with an ansatz with a power series solution of the form

$$\psi (a,\varphi )={\displaystyle \sum _{n=0}^{\infty }}{f}_n\left(\varphi \right)a^n,$$

(7)

which entered in (6) leads to a system of differential equations for the coefficients $f_n\left(\varphi \right)$

$$\begin{array}{cc}\hfill \left(1+\varphi \right)f_0^{″}+2f_0^{\prime }& =0,\hfill \\ \hfill \left(1+\varphi \right)f_1^{″}+f_1^{\prime }& =0,\hfill \\ \hfill f_2^{″}& =0,\hfill \\ \hfill \left(1+\varphi \right)f_3^{″}+f_3^{\prime }& =0,\hfill \\ \hfill \left(1+\varphi \right)f_4^{″}-2f_4^{\prime }+36k\left(1+\varphi \right)f_0& =0,\hfill \\ ⋮\hfill & \end{array}$$

which can be solved iteratively starting from $f_0$

$$\begin{array}{cc}\hfill f_0& =c_2-\frac{c_1}{\varphi +1},\hfill \\ \hfill f_1& =c_{3}log(\varphi +1)+c_4,\hfill \\ \hfill f_2& =c_5\varphi +c_4,\hfill \\ \hfill f_3& =c_5\left(\frac{{\varphi }^2}{2}+\varphi \right)+c_6,\hfill \\ \hfill f_4& =\frac{18c_{2}k{\varphi }^2}{{\hslash }^2}+\frac{36c_{2}k\varphi }{{\hslash }^2}-\frac{18c_{1}k\varphi }{{\hslash }^2}+\frac{c_6{\varphi }^3}{3}+c_6{\varphi }^2+c_6\varphi +c_7,\hfill \\ ⋮\hfill & \end{array}$$

It can be seen that imposing the previous conditions $\psi (0,\varphi )=0$ and $\psi (a,0)=0$ determines the coefficients in function of $c_1$, $c_3$, and $c_5$. In Figure 4 we illustrate this wave function obtained by numerical analysis of (6) for small values of a and $\varphi$.

For a solution around $a=\infty$ we write (6) for $b=1/a$

$$\left[{\hslash }^2{b}^6\left(1+\varphi \right)\frac{{\partial }^2}{\partial {\varphi }^2}+{\hslash }^2{b}^7\frac{{\partial }^2}{\partial \varphi \partial b}+2b^6{\hslash }^2\frac{\partial }{\partial \varphi }+36b^{2}k\left(1+\varphi \right)+6\beta {\varphi }^2\right]\psi =0,$$

it is consistent with $\psi =0$ at $b=0$. Regarding the $\varphi$-direction, we can see that as $\varphi \to \infty$ the ${\varphi }^2$ term dominates in Equation (6), with a vanishing solution in this limit.

4. Time

The problem of time in quantum cosmology amounts to the impossibility to implement dynamics by the equation of Wheeler-deWitt, as the Hamiltonian operator vanishes and a time dependent Schrödinger equation implies that the states are time independent [1,2]. As there are no external observators, it has been argued that the universe must contain its own clock. In [2] a scalar field has been proposed as a clock, provided its classical dependence of time is monotonically increasing.

In [22], we proposed an “effective” wave function

$$\Psi (a,\tau )={\left.\frac{1}{\sqrt{{\int }_0^{\infty }da{\left|\psi (a,\varphi )\right|}^2}}\psi (a,\varphi )\right|}_{\varphi =\tau }$$

(8)

$${\left|\Psi (a,\tau )\right|}^2={\left.\frac{{\left|\psi (a,\varphi )\right|}^2}{{\int }_0^{\infty }da{\left|\psi (a,\varphi )\right|}^2}\right|}_{\varphi =\tau }$$

(9)

is the conditional probability for the universe [1] to be at a, if the universe is at $\varphi =\tau$ regardless of a. The mean value of the scale factor is

$$a\left(\tau \right)={\int }_0^{\infty }a{|\Psi (a,\tau )|}^{2}da$$

(10)

In Figure 5 we show the result for this mean value for the wave function given in Figure 4, where the initial values are due to the vanishing boundary condition of the wave function at $\varphi =0$, which requires a closer analysis in the context of the the ansatz (8).

Consistency requires verification of the uncertainty relation $\Delta a\left(\tau \right)\Delta {\pi }_a\left(\tau \right)\ge \frac{1}{2}$. Further, the quantum fluctuations of a scale factor measurement must be evaluated considering that the scale factor is determined from its velocity by red shift measurements. Thus, these fluctuations must be computed from $\dot{a}=-\frac{1}{6}{a}^{-2}{\pi }_{\varphi }$. As shown in [22], the fluctuations are consistent with a classical universe today.

5. Conclusions

We have studied the classical and quantum formulation of a quadratic $f\left(R\right)$ modified theory of General Relativity based on the Starobinsky model, with a scalar field in an equivalent action, in a cosmological setting with a FRW metric. The lagrangian and hamiltonian formulations are straightforward. We consider the numerical solutions for the exact equations in both scenarios, classical and quantum, taking a compact domain for the numerical computation in the second case. With suitable boundary conditions, these solutions tend to zero at the boundaries, pointing to normalizability of the wave function, consistently with the probabilistic interpretation. The wave function is interpreted by a conditional probability as in [22], where the scalar field plays the role of time.