Scalar Product for a Version of Minisuperspace Model with Grassmann Variables

Abstract

Grassmann variables are used to formally transform a system with constraints into an unconstrained system. As a result, the Schrödinger equation arises instead of the Wheeler–DeWitt one. The Schrödinger equation describes a system’s evolution, but a definition of the scalar product is needed to calculate the mean values of the operators. We suggest an explicit formula for the scalar product related to the Klein–Gordon scalar product. The calculation of the mean values is compared with an etalon method in which a redundant degree of freedom is excluded. Nevertheless, we note that a complete correspondence with the etalon picture is not found. Apparently, the picture with Grassmann variables requires a further understanding of the underlying Hilbert space.

1. Introduction

There is a principal possibility to construct the theory of quantum gravity (QG) from the point of view that gravity is a usual physical system with constraints [1,2]. As a result, it has to be quantized using Dirac brackets [3].

The physical question arises: what type of gravity theory must be quantized? It is hardly general relativity (GR), because GR suffers from the loss of information (unitarity) in black holes [4,5] (however, see, e.g., [6,7,8]) and from the vacuum energy problem [9,10,11]. It seems possible [12] to repair GR by restricting it to a class of submanifolds without black holes [13,14]. Simultaneously, the possibility of arbitrarily choosing an energy density level appears [12,15], which removes the vacuum energy problem, at least for massless particles. Contributions of the masses of particles into vacuum energy density $\sim M_p^2{m}^2$ have to be mutually compensated for [16]. Contributions of the order of $\sim m^4$ also have to be compensated for but by taking condensates into account [15]. The resulting theory, which has features such as partial gauge fixing, preferred reference frame, and the existence of æther [17,18] (in the form of condensates and vacuum polarization) could be a suitable candidate for quantization.

Another (mathematical) question is how to realize the commutation relations corresponding to Dirac brackets, which again implies gauge fixing by defining auxiliary conditions to convert a system from the first class to the second class. So far, there is no constructive way to do so generally [19]. In particular, the implementation of Dirac brackets in 3 + 1 GR has not yet been achieved. Even for 2 + 1 gravity, Dirac brackets have a rather complicated structure [20]. Moreover, in the case of gravity quantization, the gauge-fixing conditions must be time-dependent to introduce time into the theory. Instead, one could use the quasi-Heisenberg picture [21,22,23,24], where the commutator relations corresponding to Dirac brackets are determined at some fixed moment of time, e.g., near a small-scale factor that simplifies a problem.

A more radical method is introducing Grassmann variables [25,26,27,28,29,30], which formally reduces a system with constraints to an unconstrained one. However, if one applies Grassmann variables to calculate not only the scattering amplitudes but also the mean values of operators, the questions about Hilbert spaces and scalar products arise [31,32,33].

For simplicity, the question about scalar products could be considered in a minisuperspace model example. Minisuperspace models represent examples of simple systems with constraints and are widely used [34,35,36,37,38,39,40,41,42,43] to understand the main features of gravity quantization. Without experimental data for the minisuperspace model, one would not be able to straightforwardly check different approaches to gravity quantization. Fortunately, an etalon quantization method for the minisuperspace model that “could not be wrong” exists. It consists of the explicit exclusion of the redundant degree of freedom, initially, to obtain a physical Hamiltonian [44,45] through explicit gauge fixing. Considering this “etalon” method implies that QG has to violate gauge invariance. Certainly, opposite points of view exist; e.g., loop quantum gravity promises “covariant with respect to diffeomorphisms” QG [46] (however, see some comments in [47]).

2. The Etalon Picture with the Exclusion of the Redundant Degrees of Freedom

Let us consider the action functional of gravity minimally coupled to a massless scalar field:

$$\begin{array}{c}\hfill S=\frac{1}{16\pi G}\int R\sqrt{-g}{d}^{4}x+\frac{1}{2}\int {𝜕}_{\mu }\varphi g^{\mu \nu }{𝜕}_{\nu }\varphi \sqrt{-g}{d}^{4}x,\end{array}$$

(1)

where R is the scalar curvature; Greek indexes run from 0 to 3; G is the Newton constant; $g_{\mu \nu }$ is the metric tensor, with g being its determinant. Considering a uniform, isotropic, and flat universe

$$d{s}^2=g_{\mu \nu }d{x}^{\mu }d{x}^{\nu }=a^2(N^{2}d{\eta }^2-d^2\mathit{r}),$$

(2)

where functions a and N depend only on $\eta$, reduces action (1) to

$$S=\frac{1}{2}\int \frac{1}{N}\left(-M_p^2{a}^{\prime 2}+a^2{\varphi }^{\prime 2}\right)d\eta =\int \left(-p_a{a}^{\prime }+{\pi }_{\varphi }{\varphi }^{\prime }-N\left(-\frac{1}{2}{p}_a^2+\frac{{\pi }_{\varphi }^2}{2a^2}\right)\right)d\eta ,$$

(3)

where the reduced Planck mass, $M_p=\sqrt{\frac{3}{4\pi G}}$, is used, which is set to unity everywhere further for simplicity. The Hamiltonian

$$H=N\left(-\frac{1}{2}{p}_a^2+\frac{{\pi }_{\varphi }^2}{2a^2}\right),$$

(4)

also determines the Hamiltonian constraint

$${\mathsf{\Phi }}_1=-\frac{1}{2}{p}_a^2+\frac{{\pi }_{\varphi }^2}{2a^2}=0,$$

(5)

by virtue of $\frac{\delta S}{\delta N}=0$. The time evolution of an arbitrary observable A is expressed with the Poisson brackets

$$\frac{dA}{d\eta }=\frac{𝜕A}{𝜕\eta }+\{H,A\},$$

(6)

which are defined as

$$\{A,B\}=\frac{𝜕A}{𝜕{\pi }_{\varphi }}\frac{𝜕B}{𝜕\varphi }-\frac{𝜕A}{𝜕\varphi }\frac{𝜕B}{𝜕{\pi }_{\varphi }}-\frac{𝜕A}{𝜕p_a}\frac{𝜕B}{𝜕a}+\frac{𝜕A}{𝜕a}\frac{𝜕B}{𝜕p_a}.$$

(7)

The full system of equations of motion has the form

$$\begin{array}{c}\hfill {\varphi }^{\prime }=\frac{𝜕H}{𝜕{\pi }_{\varphi }}=N\frac{{\pi }_{\varphi }}{a^2},{\pi }_{\varphi }^{\prime }=-\frac{𝜕H}{𝜕\varphi }=0,⟹{\pi }_{\varphi }=k=const,\end{array}$$

(8)

$$\begin{array}{c}\hfill \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}{a}^{\prime }=-\frac{𝜕H}{𝜕p_a}=N{p}_a,p_a^{\prime }=\frac{𝜕H}{𝜕a}=-\frac{N{k}^2}{a^3}.\end{array}$$

(9)

The additional time-dependent gauge-fixing condition

$${\mathsf{\Phi }}_2=a-\sqrt{2|{\pi }_{\varphi }|\eta }$$

(10)

can be introduced as the constraint ${\mathsf{\Phi }}_2$, which fixes N to be equal unity. The condition (10) also fixes the direction of time because the scale factor increases with time. The solutions of Equation (9) are

$$a=\sqrt{2|{\pi }_{\varphi }|\eta },\varphi =\frac{{\pi }_{\varphi }}{2|{\pi }_{\varphi }|}ln\eta +const.$$

(11)

The constraints ${\mathsf{\Phi }}_1$, ${\mathsf{\Phi }}_2$ allow reducing this simple system to a sole degree of freedom. Let us make some general notes about the exclusion of the variables and coming to the physical Hamiltonian. Let system variables be $Q_1,\dots ,Q_m,q_{m+1},\dots ,q_n$, $P_1,\dots ,P_m,p_{m+1},\dots ,p_n$ and one would like to exclude $n-m$ coordinates $q_k$ and momentums $p_k$ using the constraints and additional conditions ${\mathsf{\Phi }}_k$, to have only the variables $Q_k$ and momentums $P_k$ finally. Equating the Lagrangians before and after exclusion leads to

$$\sum _{k=1}^m{P}_{k}d{Q}_k+\sum _{k=m+1}^n{p}_{k}d{q}_k-\sum _{k=m+1}^m{\lambda }_k{\mathsf{\Phi }}_k=\sum _{k=1}^m{P}_{k}d{Q}_k-H_{phys}(\mathit{Q},\mathit{P})d\tau +dF(\mathit{Q},\mathit{P},\tau ),$$

(12)

where the left-hand part of (12) contains all the constraints, including the Hamiltonian one, whereas the right-hand side contains the full differential of some function $F(\mathit{Q},\mathit{P},\tau )$. At the constraint surface ${\mathsf{\Phi }}_k$, Equation (12) reduces to

$$H_{phys}=-\sum _{k=m+1}^n{p}_k{q}_k^{\prime }+\frac{dF(\mathit{Q},\mathit{P},\tau )}{d\tau }.$$

(13)

The function F has to be chosen in such a way that $H_{phys}$ reproduces correct equations of motion by

$${\mathit{Q}}^{\prime }=\frac{𝜕H_{phys}}{𝜕\mathit{P}},{\mathit{P}}^{\prime }=-\frac{𝜕H_{phys}}{𝜕\mathit{Q}}.$$

(14)

Let us take ${\pi }_{\varphi }$ and $\varphi$ as physical variables, then a and $p_a$ have to be excluded by the constraints. Substituting $p_a$, $a^{\prime }$ and a into (3) results in

$$L=\int \left({\pi }_{\varphi }{\varphi }^{\prime }-H_{phys}(\varphi ,{\pi }_{\varphi },\eta )\right)d\eta ,$$

(15)

where

$$H_{phys}(\varphi ,{\pi }_{\varphi },\eta )=p_a{a}^{\prime }=\frac{|{\pi }_{\varphi }|}{2\eta }.$$

(16)

Hamiltonian (16) reproduces the equations of motion (8) correctly; thus, function F in (13) equals zero for this particular case1.

The most simple and straightforward way to describe quantum evolution is to formulate the Schrödinger equation

$$i{𝜕}_{\eta }\mathsf{\Psi }={\widehat{H}}_{phys}\mathsf{\Psi }$$

(17)

with the physical Hamiltonian (16). In the momentum representation, the operators become

$${\widehat{\pi }}_{\varphi }=k,\widehat{\varphi }=i\frac{𝜕}{𝜕k}.$$

(18)

The solution of Equation (17) is written as

$$\mathsf{\Psi }(k,\eta )=C\left(k\right){\left(\frac{2}{e}|k|\eta \right)}^{-i|k|/2}.$$

(19)

It is possible to calculate the mean values of an arbitrary operator $\widehat{A}(k,i\frac{𝜕}{𝜕k})$ built from $\widehat{\varphi }=i\frac{𝜕}{𝜕k}$ and $a=\sqrt{2|k|\eta }$ with respect to a wave packet $C\left(k\right)$ in the following way

$$<C|\widehat{A}|C>=\int {\mathsf{\Psi }}^{\ast }(k,\eta )\widehat{A}\mathsf{\Psi }(k,\eta )dk.$$

(20)

Since the basic wave function ${\left(\frac{2}{e}|k|\eta \right)}^{-i|k|/2}$ contains a module of k, a singularity may arise at $k=0$ if $\widehat{A}$ contains degrees of the differential operator $\frac{𝜕}{𝜕k}$. That may violate hermicity. To avoid this, the wave packet has to be turned to zero at $k=0$. For instance, it could be taken in the Gaussian form

$$C\left(k\right)=\frac{2{\sigma }^{-5/2}}{\sqrt{3}{\pi }^{1/4}}{k}^{2}exp(-k^2/\left(2{\sigma }^2\right))$$

(21)

with the multiplier $k^2$ in the front of the exponent.

Let us come to the calculation of some mean values taking the parameter $\sigma =1$. The mean value of $a^2$ is2

$$<C|a^2|C>=\frac{16}{3\sqrt{\pi }}\eta {\int }_0^{\infty }{e}^{-k^2}{k}^{5}dk=\frac{16}{3\sqrt{\pi }}\eta .$$

(22)

The next quantity is expressed as

$$<C|a^4|C>=\frac{16}{3\sqrt{\pi }}{\eta }^2{\int }_{-\infty }^{\infty }{e}^{-k^2}{k}^{6}dk=10{\eta }^2.$$

(23)

Other mean values for this wave packet were calculated in [45].

3. Evolution in Extended Space

Another derivation of the physical Hamiltonian (16) is given by the continual integrals considering the transition amplitude from in to out states. Using canonical gauge fixing condition ${\mathsf{\Phi }}_2=a-\sqrt{2\eta |{\pi }_{\varphi }|}$ leads to

$$\begin{array}{c}\hfill <out|in>=Z=\int e^{i\int \left({\pi }_{\varphi }{\varphi }^{\prime }-p_a{a}^{\prime }-N\left(-\frac{1}{2}{p}_a^2+\frac{{\pi }_{\varphi }^2}{2a^2}\right)\right)d\eta }{\mathsf{\Pi }}_{\eta }{\Delta }_{FP}{\mathsf{\Pi }}_{\eta }\delta \left({\mathsf{\Phi }}_2\right)\mathcal{D}N\mathcal{D}{p}_a\mathcal{D}a\mathcal{D}{\pi }_{\varphi }\mathcal{D}\varphi =\\ \hfill \int e^{i\int \left({\pi }_{\varphi }{\varphi }^{\prime }-p_a{a}^{\prime }\right)d\eta }{\mathsf{\Pi }}_{\eta }{p}_a{\mathsf{\Pi }}_{\eta }\delta \left(-\frac{1}{2}{p}_a^2+\frac{{\pi }_{\varphi }^2}{2a^2}\right){\mathsf{\Pi }}_{\eta }\delta \left(a-\sqrt{2\eta |{\pi }_{\varphi }|}\right)\mathcal{D}{p}_a\mathcal{D}a\mathcal{D}{\pi }_{\varphi }\mathcal{D}\varphi =\\ \hfill \int e^{i\int \left({\pi }_{\varphi }{\varphi }^{\prime }-|{\pi }_{\varphi }|/\left(2\eta \right)\right)d\eta }\mathcal{D}{\pi }_{\varphi }\mathcal{D}\varphi .\end{array}$$

(24)

Under the derivation of (24), it was used that ${\Delta }_{FP}=\{{\mathsf{\Phi }}_1,{\mathsf{\Phi }}_2\}=p_a$, and ${\mathsf{\Pi }}_{\eta }\delta \left(-\frac{1}{2}{p}_a^2+\frac{{\pi }_{\varphi }^2}{2a^2}\right)$$=\frac{{\mathsf{\Pi }}_{\eta }\delta \left(p_a-|{\pi }_{\varphi }|/a\right)}{{\mathsf{\Pi }}_{\eta }{p}_a}$. From (24) follows the formula (16) for a physical Hamiltonian. A principle of derivation is clear: obtaining an expression of a kind $\int e^{i\int \left({\sum }_i{p}_i{q}_i^{\prime }-H(\mathit{p},\mathit{q})\right)d\eta }\mathcal{D}\mathit{p}\mathcal{D}\mathit{q}$ without any pre-exponential factors and extracting $H(\mathit{p},\mathit{q})$ to use in the Schrödinger equation.

However, the etalon picture with $H_{phys}$ cannot be applied in the general case to QG because one cannot resolve the constraints. It is believed that the Grassmann variables allow writing the Lagrangian in a form where there are no constraints [27,28,50,51,52]. Using non-canonical gauge fixing [31,53] leads to

$$Z=\int e^{i\int \left({\pi }_{\varphi }{\varphi }^{\prime }-p_a{a}^{\prime }-N\left(-\frac{1}{2}{p}_a^2+\frac{{\pi }_{\varphi }^2}{2a^2}\right)\right)d\eta }{\mathsf{\Pi }}_{\eta }\frac{\delta F}{\delta \epsilon }{\mathsf{\Pi }}_{\eta }\delta \left(F\right)\mathcal{D}{p}_a\mathcal{D}a\mathcal{D}{\pi }_{\varphi }\mathcal{D}\varphi \mathcal{D}N,$$

(25)

where $F\left(N\right)$ is a gauge-fixing function.

The action (3) is invariant under the infinitesimal gauge transformation:

$$\begin{array}{c}\hfill \tilde{a}=a+\delta a=a+\epsilon a^{\prime },\end{array}$$

(26)

$$\begin{array}{c}\hfill \tilde{\varphi }=\varphi +\delta \varphi =\varphi +\epsilon {\varphi }^{\prime },\end{array}$$

(27)

$$\begin{array}{c}\hfill \tilde{N}=N+\delta N=N+{\left(N\epsilon \right)}^{\prime },\end{array}$$

(28)

where $\epsilon$ is an infinitesimal function of time. If one takes the differential gauge condition $F=N^{\prime }=0$, then (28) follows in

$$\delta F=\delta N^{\prime }={\left(N\epsilon \right)}^{\prime \prime },$$

(29)

and the Faddev–Popov determinant [53] takes the form of ${\Delta }_{FP}=\frac{\delta F}{\delta \epsilon }=\frac{\delta N^{\prime }}{\delta \epsilon }=N^{\prime \prime }+2N^{\prime }\frac{𝜕}{𝜕\eta }+N\frac{{𝜕}^2}{𝜕{\eta }^2}$. The functional (25) could be rewritten as

$$\begin{array}{c}\hfill Z=\int e^{i\int \left({\pi }_{\varphi }{\varphi }^{\prime }-p_a{a}^{\prime }-N\left(-\frac{1}{2}{p}_a^2+\frac{{\pi }_{\varphi }^2}{2a^2}\right)-\overline{\theta }{\left(N\theta \right)}^{″}\right)d\eta }{\mathsf{\Pi }}_{\eta }\delta \left(N^{\prime }\left(\eta \right)\right)\mathcal{D}N\mathcal{D}{p}_a\mathcal{D}a\mathcal{D}{\pi }_{\varphi }\mathcal{D}\varphi \mathcal{D}\theta \mathcal{D}\overline{\theta }=\\ \hfill \int e^{i\int \left({\pi }_{\varphi }{\varphi }^{\prime }-p_a{a}^{\prime }-N_0\left(-\frac{1}{2}{p}_a^2+\frac{{\pi }_{\varphi }^2}{2a^2}-{\overline{\theta }}^{\prime }{\theta }^{\prime }\right)\right)d\eta }\mathcal{D}{p}_a\mathcal{D}a\mathcal{D}{\pi }_{\varphi }\mathcal{D}\varphi \mathcal{D}\theta \mathcal{D}\overline{\theta }=\\ \hfill \int e^{i\int \left({\pi }_{\varphi }{\varphi }^{\prime }-p_a{a}^{\prime }+{\overline{\theta }}^{\prime }{\pi }_{\overline{\theta }}+{\pi }_{\theta }{\theta }^{\prime }-N_0\left(-\frac{1}{2}{p}_a^2+\frac{{\pi }_{\varphi }^2}{2a^2}+{\pi }_{\theta }{\pi }_{\overline{\theta }}\right)\right)d\eta }\mathcal{D}{p}_a\mathcal{D}a\mathcal{D}{\pi }_{\varphi }\mathcal{D}\varphi \mathcal{D}{\pi }_{\theta }\mathcal{D}{\pi }_{\overline{\theta }}\mathcal{D}\theta \mathcal{D}\overline{\theta },\end{array}$$

(30)

where using the Grassmann variables [53] in the first equality of (30) raises the Faddeev–Popov determinant into an exponent. Here, a Grassmann number $\overline{\theta }$ is considered as a complex conjugate to $\theta$. Integration over N has been performed explicitly in (30). If $N\left(\eta \right)$ is discretized over the interval $\Delta \eta$, the term containing a product of the delta functions ${\mathsf{\Pi }}_{\eta }\delta \left(N^{\prime }\left(\eta \right)\right)$ takes the form

$$\int \cdots \delta \left(\frac{N_0-N_1}{\Delta \eta }\right)\delta \left(\frac{N_1-N_2}{\Delta \eta }\right)\dots \delta \left(\frac{N_{k-1}-N_k}{\Delta \eta }\right)d{N}_1\cdots d{N}_{k-1}\sim \Delta {\eta }^{k-1}\delta (N_0-N_k),$$

(31)

i.e., an initial value of $N_0$ has to equal a final value $N_k$, for instance, one may take $N_0=1$. For deducing third equality of (30), see Appendix A. Extracting the Lagrangian from (30) gives

$$L={\pi }_{\varphi }{\varphi }^{\prime }-p_a{a}^{\prime }+{\overline{\theta }}^{\prime }{\pi }_{\overline{\theta }}+{\pi }_{\theta }{\theta }^{\prime }-\left(-\frac{1}{2}{p}_a^2+\frac{{\pi }_{\varphi }^2}{2a^2}+{\pi }_{\theta }{\pi }_{\overline{\theta }}\right).$$

(32)

The action (32) is a fixed gauge action with no Hamiltonian constraint, but instead, the ghost (Grassmann) variables arise in (32).

Following Vereshkov and Shestakova et al. [27,28], one may consider the Hamiltonian

$$H=-\frac{1}{2}{p}_a^2+\frac{{\pi }_{\varphi }^2}{2a^2}+{\pi }_{\theta }{\pi }_{\overline{\theta }}$$

(33)

as describing the quantum evolution of a system.

To quantize the system, the anticommutation relation has to be introduced for the Grassmann variables

$$\{{\pi }_{\theta },\theta \}=-i,\{{\pi }_{\overline{\theta }},\overline{\theta }\}=-i.$$

(34)

In the particular representation $\alpha =lna$, ${\widehat{p}}_{\alpha }=i\frac{𝜕}{𝜕\alpha }$, $\widehat{\varphi }=i\frac{𝜕}{𝜕k}$, ${\widehat{\pi }}_{\varphi }=k$, ${\widehat{\pi }}_{\theta }=-i\frac{𝜕}{𝜕\theta }$, ${\widehat{\pi }}_{\overline{\theta }}=-i\frac{𝜕}{𝜕\overline{\theta }}$, the Schrödinger equation reads as

$$i\frac{𝜕}{𝜕\eta }\psi =\left(\frac{1}{2}{e}^{-2\alpha }\left(\frac{{𝜕}^2}{𝜕{\alpha }^2}+k^2\right)-\frac{𝜕}{𝜕\theta }\frac{𝜕}{𝜕\overline{\theta }}\right)\psi ,$$

(35)

where the operator ordering in the form of the two-dimensional Laplacian [54] has been used3. It should be supplemented by the scalar product

$$<{\psi }_1|{\psi }_2>=\int {\psi }_1^{\ast }(\eta ,k,\alpha ,\overline{\theta },\theta ){\psi }_2(\eta ,k,\alpha ,\overline{\theta },\theta )e^{2\alpha }d\alpha dkd\theta d\overline{\theta },$$

(36)

where the measure $e^{2\alpha }$ arises due to the hermicity requirement [26,54]. This measure is a consequence of a minisuperspace metric if the classical Hamiltonian is written in the form of $H=\frac{1}{2}{\mathcal{G}}^{ij}{p}_i{p}_j+{\pi }_{\theta }{\pi }_{\overline{\theta }}$ with $p_i\equiv \{p_{\alpha },{\pi }_{\varphi }\}$, ${\mathcal{G}}^{ij}=diag\{-e^{-2\alpha },e^{-2\alpha }\}$. Thus, the measure takes the form $\sqrt{\mathcal{G}}=e^{2\alpha }$, $\mathcal{G}=|det{\mathcal{G}}_{ij}|$ and the Laplacian $\frac{1}{\sqrt{\mathcal{G}}}\frac{𝜕}{𝜕q^i}\sqrt{\mathcal{G}}{\mathcal{G}}^{ij}\frac{𝜕}{𝜕q^i}=e^{-2\alpha }\left(-\frac{{𝜕}^2}{𝜕{\alpha }^2}+{\pi }_{\varphi }^2\right)$ is self-adjoined [54] with this measure. Formal solutions of the Equation (35) can be written as

$$\psi (\eta ,k,\alpha ,\overline{\theta },\theta )=(\overline{\theta }+\theta )u(\eta ,k,\alpha )+i(\overline{\theta }-\theta )v(\eta ,k,\alpha ),$$

(37)

where the functions u and v satisfy the equation

$$i\frac{𝜕}{𝜕\eta }u={\widehat{H}}_{0}u$$

(38)

with

$${\widehat{H}}_0=\frac{1}{2}{e}^{-2\alpha }\left(\frac{{𝜕}^2}{𝜕{\alpha }^2}+k^2\right).$$

(39)

Then, the scalar product (36) reduces to

$$<{\psi }_1|{\psi }_2>=-2i\int \left(u_1^{\ast }{v}_2-v_1^{\ast }{u}_2\right)e^{2\alpha }d\alpha dk.$$

(40)

Although the constraint $H_0=0$ formally disappears from the theory, one may think that the space of solutions of the Wheeler–DeWitt equation (WDW) still plays a role [33]. Otherwise, the question of correspondence with the classical theory, where the Hamiltonian constraint holds, arises. We would like to relate the space of the functions, satisfying the Schrödinger Equation (35) with the functions $\chi$ satisfying the equation $H_0\chi =0$, i.e., the WDW equation. The operator ${\widehat{H}}_0$ (39) has the Klein–Gordon form. Thus, the Klien–Gordon-type scalar product has to be used. According to this hypothesis, let us represent the functions u, v as

$$\begin{array}{c}\hfill \hspace{1em}\hspace{1em}\hspace{1em}v(\alpha ,k)=e^{-i{H}_0\eta }{\widehat{D}}^{1/4}\chi (\alpha ,k),\end{array}$$

(41)

$$\begin{array}{c}\hfill u(\alpha ,k,\eta )=e^{-i{H}_0\eta }{\widehat{D}}^{-1/4}\delta (\alpha -{\alpha }_0)\frac{𝜕}{𝜕\alpha }\chi (\alpha ,k),\end{array}$$

(42)

where the operator $\widehat{D}=-\frac{{𝜕}^2}{𝜕{\varphi }^2}$, or $D=k^2$ in the representation (18) and

$$\chi (\alpha ,k)=\frac{e^{-i\alpha |k|-{\alpha }_0}}{\sqrt{2|k|}}C\left(k\right),$$

(43)

and, as in (19), only a half-space corresponding to the negative frequencies’ solutions of the WDW equation is taken because only in this case does the Klein–Gordon product imply a positive definite norm of a state. The operator $\widehat{D}$ (see Appendix in [57]) is a necessary attribute of the scalar product for the Klein–Gordon equation to obtain hermicity. It should be noted that in fact, the function v does not depend on the time $\eta$ because ${\widehat{H}}_0\chi =0$ and $\widehat{D}$ commutes with $H_0$. Thus, the time evolution arises only due to function u, or more accurately, due to the presence of the Dirac delta function in (42).

Thus, the scalar product (40) reduces to

$$<{\psi }_1|{\psi }_2>=-2i\int \left(\frac{𝜕{\chi }_1^{\ast }}{𝜕\alpha }{\chi }_2-{\chi }_1^{\ast }\frac{𝜕{\chi }_2}{𝜕\alpha }\right)e^{2\alpha }{|}_{\alpha ={\alpha }_0}dk.$$

(44)

The expression for the mean value of an operator $\widehat{A}$ has the form:

$$<\psi |\widehat{A}|\psi >=-2i\int e^{2\alpha }\left(u^{\ast }\widehat{A}v-v^{\ast }\widehat{A}u\right){|}_{\alpha ={\alpha }_0\to -\infty }dk,$$

(45)

where u, v are given by (41) and (42), and it is assumed that an operator $\widehat{A}$ does not contain the ghost variables $\theta$, $\overline{\theta }$, that is expected for physical operators. The limit $\alpha \to -\infty$ in (45) implies that an evolution begins at $\eta =0$ when $a=0$ and $\alpha =lna$ tends to $-\infty$.

The evaluation of the mean values is illustrated schematically in Figure 1. Initially, we start with the negative frequency functions $\chi$, $\frac{𝜕\chi }{𝜕\alpha }$ both satisfying $H_0\chi =0$. A direct product of these half-spaces is taken. Then, the function $\frac{𝜕\chi }{𝜕\alpha }$ is multiplied by $\delta (\alpha -{\alpha }_0)$ and runs into an extended space, where “evolution” occurs, and, thus, the mean values of the operators can be evaluated.

Both Schrödinger and Heisenberg pictures are possible with this scalar product. In the Heisenberg picture, the time-dependent operators have the form

$$\widehat{A}\left(\eta \right)=e^{i{\widehat{H}}_0\eta }\widehat{A}{e}^{-i{\widehat{H}}_0\eta },$$

(46)

while the functions u and v have to be used without multiplier $e^{-i{\widehat{H}}_0\eta }$.

4. Expectation Values of Scale Factor Degrees

The simplest way to test a theory is to compare it with the etalon picture by calculating the mean value of the squared scale factor, which has to be equal to $\frac{16}{3\sqrt{\pi }}\eta$ according to (22). To do that, it is sufficient to expand $e^{-i{\widehat{H}}_0\eta }\approx 1-i{\widehat{H}}_0\eta -\frac{1}{2}{\widehat{H}}_0^2{\eta }^2$ in (41) and (42) and perform the calculation according (45). It turns out that the mean value of $a^2=e^{2\alpha }$ actually coincides with that given by (22). The next test is the calculation of $<C|a^4|C>$. The result of the calculation is

$$<C|a^4|C>=2{\eta }^2,$$

(47)

while the etalon model gives another value (23). The origin of this discrepancy could be better seen in the Heisenberg picture. Evolution equations for the Heisenberg operators follow from the operator commutators with the Hamiltonian (39)

$$\frac{d{\widehat{a}}^2}{d\eta }=i[{\widehat{H}}_0,{\widehat{a}}^2].$$

(48)

It is possible to guess a solution for this particular case:

$${\widehat{a}}^2\left(\eta \right)=e^{2\alpha }+2\eta e^{-\alpha }{\widehat{p}}_{\alpha }{e}^{\alpha }-2{\eta }^2{\widehat{H}}_0,$$

(49)

where we define ${\widehat{p}}_{\alpha }=i\frac{𝜕}{𝜕\alpha }$4.

Actually, the calculation of the commutator (48) using (39), (49) gives

$$i[{\widehat{H}}_0,{\widehat{a}}^2\left(\eta \right)]=2e^{-\alpha }{\widehat{p}}_{\alpha }{e}^{\alpha }-4\eta {\widehat{H}}_0,$$

(50)

which is exactly equal to the derivative of (49) over $\eta$. Under calculation of the mean value of $<C|{\widehat{a}}^2|C>$, the third term in (49) does not contribute, and the result coincides with that of the etalon method. However, under the calculation of $<C|{\widehat{a}}^4|C>$, the first and third terms in (49) play a role, and the discrepancy with the etalon method arises. One can calculate the mean values of the other degrees of a, which are presented in

Table 1. The expectation values $<C|a^{2n}|C>=k_{2n}{\eta }^n$ for the wave packet (21).

$2\mathit{n}$	2	4	6	8	10	12	14
$\sqrt[n]{k_{2n}}$ for the etalon model	$\frac{16}{3\sqrt{\pi }}$	$\sqrt{10}$	$\sqrt[3]{\frac{64}{\sqrt{\pi }}}$	$\sqrt[4]{140}$	$\sqrt[5]{\frac{1024}{\sqrt{\pi }}}$	$\sqrt[6]{2520}$	$\sqrt[7]{\frac{20480}{\sqrt{\pi }}}$
$\sqrt[n]{k_{2n}}$ for the model with the Grassmann variables	$\frac{16}{3\sqrt{\pi }}$	$\sqrt{2}$	$\sqrt[3]{\frac{512}{3\sqrt{\pi }}}$	$\sqrt[4]{876}$	$\sqrt[5]{\frac{7936}{3\sqrt{\pi }}}$	$\sqrt[6]{118280}$	$\sqrt[7]{\frac{1172480}{\sqrt{\pi }}}$

. It is interesting to plot the values of $\sqrt[n]{k_{2n}}=\frac{1}{\eta }\sqrt[n]{<C|{\widehat{a}}^{2n}\left(\eta \right)|C>}$, which is shown in Figure 2.

Recently, analog models of QG and minisuperspace have been discussed [58,59]. In Appendix B, quantization of a particle-clock is considered, which has some similar features to the minisuperspace model but does not have an operator ordering issue.

5. Discussion and Conclusions

A reasonable expression for the scalar product using the Grassmann variables is suggested. It establishes a relation of a picture with the Grassmann variables to the Klein–Gordon scalar product and allows calculating the mean values of operators in both Schrödinger and Heisenberg pictures, which give the same results. However, it is shown that the mean values of $a^{2n}$ are different for $n>1$ than those calculated in the etalon method, implying an explicit exclusion of the superfluous degrees of freedom. One may guess that the above methods could have different Hilbert spaces. That means that the different wave packets have to be taken for these methods to obtain the same set of operator mean values. Here, we cannot find a wave packet $\tilde{C}\left(k\right)$, which would give the same mean values as a wave packet $C\left(k\right)$ for the etalon method.

The possible influence of the Zitterbewegung phenomenon in extended space was investigated in Appendix C but without a breakthrough in the results achieved. It should be noted that the quasi-Heisenberg picture corresponds entirely with the etalon method [45].

One of the possible ways to correct the picture with the Grassmann variables is to assume that the operators of physical observables act not only in k and $\alpha$ space but also in the extended space of the Grassmann variables $\theta$, $\overline{\theta }$. This hypothesis needs further investigation5 as well as the general issue of the scalar product for the approach with the Grassmann variables.