Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Cosmic-Time Quantum Mechanics and the Passage-of-Time Problem
Next Article in Special Issue
Particles of Negative and Zero Energy in Black Holes and Cosmological Models
Previous Article in Journal
Solar Sail Trajectories to Earth’s Trojan Asteroids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

DeWitt Boundary Condition in One-Loop Quantum Cosmology

by
Giampiero Esposito
1,2
1
Dipartimento di Fisica “Ettore Pancini”, Complesso Universitario di Monte S. Angelo, Via Cintia Edificio 6, 80126 Napoli, Italy
2
INFN Sezione di Napoli, Complesso Universitario di Monte S. Angelo, Via Cintia Edificio 6, 80126 Napoli, Italy
Universe 2023, 9(4), 187; https://doi.org/10.3390/universe9040187
Submission received: 9 March 2023 / Revised: 11 April 2023 / Accepted: 13 April 2023 / Published: 14 April 2023
(This article belongs to the Special Issue The Friedmann Cosmology: A Century Later)

Abstract

:
DeWitt’s suggestion that the wave function of the universe should vanish at the classical Big Bang singularity is considered here within the framework of one-loop quantum cosmology. For pure gravity at one loop about a flat four-dimensional background bounded by a 3-sphere, three choices of boundary conditions are considered: vanishing of the linearized magnetic curvature when only transverse-traceless gravitational modes are quantized; a one-parameter family of mixed boundary conditions for gravitational and ghost modes; and diffeomorphism-invariant boundary conditions for metric perturbations and ghost modes. A positive ζ ( 0 ) value in these cases ensures that, when the three-sphere boundary approaches zero, the resulting one-loop wave function approaches zero. This property may be interpreted by saying that, in the limit of small three-geometry, the resulting one-loop wave function describes a singularity-free universe. This property holds for one-loop functional integrals, which are not necessarily equivalent to solutions of the quantum constraint equations.

1. Introduction

After the birth of relativistic cosmology thanks to Friedmann’s work [1], and the subsequent proof of the singularity theorems of Penrose, Hawking, and Geroch [2,3,4,5,6,7], it became well-accepted by the scientific community that classical general relativity leads to the occurrence of cosmological singularities (a spacetime being singular if it is timelike or null geodesically incomplete) in a generic way. Since then, several developments occurred, and in particular, we here mention what follows.
(i)
At the classical level, the work of Christodoulou and Klainerman [8] led to the discovery of asymptotically flat spacetimes, which are timelike geodesically complete.
(ii)
At the quantum level, DeWitt [9] proposed to look at the behavior of the wave function of the universe in correspondence with the classical Big Bang singularity. Such a proposal was assessed in the outstanding work in Ref. [10].
(iii)
Over many years, various concepts of singularity have been conceived, as can be seen in an important review of Kamenshchik [11].
Moreover, in the literature on quantum gravity and quantum cosmology, several approaches were developed to study the possible quantum origin of spacetime geometry. One-loop effects in the early universe were investigated in detail, especially with the help of ζ -function methods. It is the aim of our review to describe them and then discuss their relevance for the singularity issue in cosmology. The structure of the paper is as follows.
Section 2 presents in detail a ζ ( 0 ) calculation when only transverse-traceless perturbations are considered, with boundary conditions requiring the vanishing of linearized magnetic curvature on the three-sphere boundary. Section 3 discusses a one-parameter family of ζ ( 0 ) values obtained with mixed boundary conditions for metric perturbations and ghost fields. Section 4, Section 5, Section 6 and Section 7 outline the basic steps of the ζ ( 0 ) calculation with diffeomorphism-invariant boundary conditions. Open problems are discussed in Section 8.

2. Linearized Magnetic Curvature Vanishing on S 3

We study pure gravity at one loop about flat Euclidean four-space with a three-sphere boundary of radius a, because when a 0 , this is the limiting case of a four-sphere geometry bounded by a three-sphere [12]. The prefactor of the semiclassical wave function is given by the following (with I 2 denoting the part of the action quadratic in metric perturbations)
P ( a ) = e I 2 [ γ ] D γ ,
which is a functional integral over all metric perturbations γ a b that are regular at the origin τ = 0 and satisfy a given boundary condition at τ = a . Integration is here restricted to the physical degrees of freedom, which are found by using the Hamiltonian formulation with the following transverse-traceless choice of gauge condition:
i ( D i γ ) i j = 0 , k γ k k = 0 .
These relations pick out the transverse-traceless tensor hyperspherical harmonics G i j ( n ) ( ϕ k ) multiplied by functions of the radial coordinate τ . Hence, we write
γ i j = γ i j T T = n = 3 q n ( τ ) G i j ( n ) ( ϕ k ) .
Our work in Ref. [13] studied the Breitenlohner–Freedman–Hawking [14,15] local boundary conditions for fields of spin 0 , 1 2 , 1 , 3 2 , 2 . For gravity, these imply that the linearized magnetic curvature should vanish at the boundary. Our detailed analysis in Section 7.3 of Ref. [13] never appeared in any journals, and hence it is of interest to our review article.
The action that is quadratic in metric perturbations involves a second-order elliptic operator A with eigenvalues λ n , for which one can define a spectral ζ -function
ζ A ( s ) = Tr L 2 A s = n ( λ n ) s .
Eventually, as was shown by Schleich [12], the prefactor of the semiclassical wave function of Equation (1), with γ having the form (3), can be expressed as
P ( a ) = 1 det f f 4 π l p 2 μ 2 = D a ζ ( 0 ) ,
where f f is the Laplacian acting on transverse-traceless metric perturbations, and ζ ( 0 ) is the value at s = 0 of the analytic continuation of the spectral ζ -function (4) (for further details, see now the Appendix A on the one-loop approximation). Thus, within a functional-integral framework, the wave function of the universe may fulfill the DeWitt boundary condition if and only if ζ ( 0 ) is positive.
The linearized magnetic curvature for gravity is defined from the Weyl tensor C and from the normal n to the boundary according to the rule (with summation over repeated tensor indices)
B i j 1 2 ε j μ k l C k l i ν n μ n ν ,
and it can only vanish on S 3 if [13]
n = 3 d q n d τ ( a ) ϵ j k l G i l k ( n ) G i k l ( n ) = 0 .
The only condition on the modes that ensures the validity of (6) is
d q n d τ ( a ) = 0 , n 3 .
We are now interested in evaluating ζ ( 0 ) using (7). Thus, after setting τ = t , we study the kernel of the heat equation for the operator
P n d 2 d t 2 1 t d d t ( n 2 1 ) t 2 , n 3 ,
which results from studying the Laplacian on transverse-traceless metric perturbations. On denoting by E > 0 the eigenvalues of P n , we find that its eigenfunctions regular at the origin are (up to a multiplicative constant)
u n ( t ) = t J n ( E t ) = q n ( t ) .
Thus, the boundary condition (7) implies the eigenvalue condition
J n ( E a ) + E a J ˙ n ( E a ) = 0 , n 3 .
This equation is of the general kind studied in Ref. [16]. Setting now a = 1 for simplicity, we define the function
F n ( z ) J n ( z ) + z J ˙ n ( z ) , n 3 .
Of course, the consideration of such F n ( z ) is suggested by (10). It only has real simple zeros apart from z = 0 (page 482 of Ref. [17]). The basic idea is now the following [16]. Given the zeta-function at large x
ζ ( s , x 2 ) n ( λ n + x 2 ) s ,
one has, in four dimensions (see Theorem 2 on page 6 of Ref. [18]),
Γ ( 3 ) ζ ( 3 , x 2 ) = 0 t 2 e x 2 t G ( t ) d t n = 0 B n Γ 1 + n 2 x n 2 ,
where we have used the asymptotic expansion of the heat kernel for t 0 + , i.e.,
G ( t ) n = 0 B n t n 2 2 .
Strictly speaking, since we have not proved general results on the existence of the asymptotic expansion of the heat kernel, our Formula (14) could be initially regarded as an assumption. However, existence theorems hold for the problems studied in this paper [19,20].
On the other hand, one also has the identity
Γ ( 3 ) ζ ( 3 , x 2 ) = p = 0 N p 1 2 x d d x 3 log ( ( i x ) p F p ( i x ) ) ,
where N p is the degeneracy of the problem. Thus the comparison of (13) and (15) can yield the coefficients B n and in particular ζ ( 0 ) = B 4 , provided we carefully perform a uniform Debye expansion of F p ( i x ) . It should be emphasized that this technique seems to be the most efficient. In fact, by using this algorithm, Moss [16] was able to compute ζ ( 0 ) for a real scalar field subject to Robin boundary conditions, whereas the technique of Kennedy based on charge layers on the plane tangent to S 3 failed to provide such a value [21,22]. Indeed, the eigenvalue condition (10) is of the Robin type (just set β = 1 in Equation (22) of Ref. [16]). Thus, on passing to the variable x i x and then defining α p p 2 + x 2 , C log ( 2 π ) , we can write
log ( ( i x ) p F p ( i x ) ) C p log ( p + α p ) + 1 2 log ( α p ) + α p + n = 1 r = 0 n a n r p 2 r α p n 2 r .
The coefficients a n r in (16) can be computed by comparison using the formula
n = 1 r = 0 n a n r t 2 r = m = 1 a m ( t ) ,
because the a m ( t ) values are known polynomials in t arising from uniform asymptotic expansions of Bessel functions and their first derivatives. Thus, setting β = 1 in the Formulae (29)–(31) of Ref. [16] for the a m ( t ) , we find in our case that
a 10 = 5 8 , a 11 = 7 24 ,
a 20 = 3 16 , a 21 = 1 8 , a 22 = 7 16 ,
a 30 = 17 384 , a 31 = 389 640 , a 32 = 203 128 , a 33 = 1463 1152 ,
plus infinitely many other coefficients that we do not strictly need here. We can now insert (16)–(20) into (15), apply three times the differential operator 1 2 x d d x , and finally use the contour formula for positive integer values of k [16]
p = 0 p 2 k α p 2 k m = Γ k + 1 2 Γ m 2 1 2 2 Γ k + m 2 x 1 m , k = 1 , 2 , 3 , ,
and the known properties of the Γ -function [23]. Now, writing the asymptotic expansion of the right-hand side of (15) in the form
Γ ( 3 ) ζ ( 3 , x 2 ) n = 0 b n x n 2 ,
the comparison with (15) shows that
ζ ( 0 ) = B 4 = b 4 2 = ζ I ( 0 ) + ζ I I ( 0 ) ,
since it is well-known that the asymptotic expansion, if it exists, is unique. The two contributions to ζ ( 0 ) are obtained from the following formulae:
Γ ( 3 ) ζ ( 3 , x 2 ) [ σ 1 + σ 2 ] n = 0 b n x n 2 ,
σ 1 p = 0 N p 1 2 x d d x 3 p log ( p + α p ) + 1 2 log ( α p ) + α p ,
σ 2 p = 0 N p 1 2 x d d x 3 n = 1 r = 0 n a n r p 2 r α p n 2 r .
Bearing in mind (15) and (16), we write (24)–(26) because we can apply Theorem 3 on page 7 of Ref. [18], concerning the differentiation of asymptotic expansions.
Thus, ζ I ( 0 ) (respectively, ζ I I ( 0 ) ) is half the coefficient of x 6 in the asymptotic expansion of σ 1 (respectively, σ 2 ). We first study the asymptotic expansion of σ 2 , since it is easier to perform this calculation. In our problem, the degeneracy N p is [12]
N p = 0 p = 0 , 1 , 2 , N p = 2 ( p 2 4 ) p 3 .
This is why we find
σ 2 n = 1 r = 0 n a n r r + n 2 r + n 2 + 1 r + n 2 + 2 [ ( G H ) ( r , x , n ) ] ,
where, setting A = 8 , B = 2 (cf. (27)), we have, using also (21),
G ( r , x , n ) = p = 0 ( A + B p 2 ) p 2 r α p n 2 r 6 = O ( x n 6 ) + A 2 Γ r + 1 2 Γ n 2 + 5 2 Γ r + n 2 x 5 n r + n 2 r + n 2 + 1 r + n 2 + 2 + B 2 Γ r + 3 2 Γ n 2 + 3 2 Γ r + n 2 x 3 n r + n 2 r + n 2 + 1 r + n 2 + 2 ,
H ( r , x , n ) = p = 0 2 2 ( p 2 4 ) p 2 r α p n 2 r 6 = 6 x n 2 r 6 1 + 1 x 2 n 2 r 3 8 δ r 0 x n 6 .
Thus, H ( r , x , n ) gives rise to terms in (28) that contain x k with k 7 , and it does not contribute to ζ I I ( 0 ) . This is why (28) and (29) lead to
ζ I I ( 0 ) = 1 2 [ A ( a 10 + a 11 ) B ( a 30 + a 31 + a 32 + a 33 ) ] .
The insertion of (18), (20) and (27) into (31) finally yields
ζ I I ( 0 ) = 11 3 121 360 = 1199 360 .
The calculation of (25) is more involved. By performing the three derivatives and using the identity 1 2 x d α p d x = 1 2 α p , we find
1 2 x d d x 3 log 1 p + α p = ( p + α p ) 3 α p 3 9 8 p α p 4 3 8 p 2 α p 5 .
This intermediate step is very important because it proves that by summing over all integer values of p from 0 to , we obtain a convergent series. However, to be able to perform the ζ ( 0 ) calculation, it is convenient to use the identity
( p + α p ) 3 = ( α p p ) 3 x 6 .
Upon inserting (34) into (33) and re-expressing p 2 as α p 2 x 2 , we obtain
1 2 x d d x 3 [ p log ( p + α p ) ] = p x 6 + p 2 x 6 α p 1 + p 2 2 x 4 α p 3 + 3 8 p 2 x 2 α p 5 M ( x , α p , p ) ,
which implies
σ 1 p = 0 N p M ( x , α p , p ) + σ 1 [ σ 1 + σ 1 ] ,
where
σ 1 = p = 0 N p α p 6 2 3 8 α p 5 = p = 0 ( A + B p 2 ) α p 6 2 + 3 8 α p 5 + p = 0 2 ( A + B p 2 ) α p 6 2 3 8 α p 5 .
The infinite sum on the right-hand side of (37) contributes to ζ ( 0 ) only through the following part:
p = 0 A 2 α p 6 = A 2 x 6 2 + π 2 3 ! ! 4 ! ! x 5 .
The result (38) is proved by applying the Euler–Maclaurin formula [18] to the calculation of p = 0 ( p 2 + x 2 ) 3 , and then using the Formula (3.249.1) on page 294 of Ref. [24]. In addition, the finite sum on the right-hand side of (37) contributes to ζ ( 0 ) . In fact, one finds (we have x ) that
p = 0 2 ( A + B p 2 ) α p 6 2 3 8 α p 5 = A 2 + B 2 x 6 1 3 x 2 + 6 x 4 + . . . A 2 x 6 3 8 A x 5 3 8 ( A + B ) x 5 1 5 2 x 2 + 35 8 x 4 + . . . ,
which implies that the total contribution of σ 1 to ζ ( 0 ) is given by
ζ I b ( 0 ) = 1 2 A B 2 + A 8 = 7 2 1 = 5 2 .
Thus, we have so far
ζ ( 0 ) = ζ I ( 0 ) + ζ I I ( 0 ) = ζ I a ( 0 ) + ζ I b ( 0 ) + ζ I I ( 0 ) ,
where
ζ I b ( 0 ) + ζ I I ( 0 ) = 5 2 + 1199 360 .
It now remains to compute ζ I a ( 0 ) , i.e., the contribution to ζ ( 0 ) due to σ 1 in (36). Indeed, one has
σ 1 A p = 0 M ( x , α p , p ) + B p = 0 p 2 M ( x , α p , p ) p = 0 2 ( A + B p 2 ) M ( x , α p , p ) .
Let us now denote by Σ ( a ) , Σ ( b ) and Σ ( c ) the three sums on the right-hand side of (43). Both Σ ( a ) and Σ ( b ) contain divergent parts in view of (35). These fictitious divergences may be regularized by dividing by α p 2 s and then taking the limit as s tends to zero, as shown in Ref. [16]. It might not appear a priori obvious that this technique leads to unambiguous results, since the limit s 0 is a delicate mathematical point. However, a fundamental consistency check is presented in Section 7.4 of Ref. [13] for all one-loop calculations involving only physical degrees of freedom of bosonic fields, showing that the method is correct. In performing the calculation, we must use the contour Formula (21) and also the following asymptotic expansion [16]:
p = 0 p α p 1 n x 1 n π r = 0 2 r r ! B ˜ r x r Γ r 2 + 1 2 Γ n 2 1 2 + r 2 2 Γ 1 2 + n 2 cos r π 2 ,
where B ˜ 0 = 1 , B ˜ 1 = 1 2 , B ˜ 2 = 1 6 , B ˜ 4 = 1 30 etc., are Bernoulli numbers. Thus, using the label R for the regularized quantities, we define
Σ R ( a ) A [ x 6 lim s 0 p = 0 p α p 1 ( 2 s 1 ) + x 6 lim s 0 p = 0 p 2 α p 2 ( 2 s 1 ) + x 4 2 lim s 0 p = 0 p 2 α p 2 ( 2 s + 1 ) + 3 8 x 2 lim s 0 p = 0 p 2 α p 2 ( 2 s + 3 ) ] .
In view of (44), the first limit in (45) gives the following contribution to ζ ( 0 ) :
δ 1 = A 2 B ˜ 2 π Γ 3 2 = A 24 = 1 3 ,
whereas the other limits in (45) do not contribute to ζ ( 0 ) in view of (21), because one only obtains terms proportional to x 4 .
Moreover, bearing in mind the identity
p = 0 p 3 α p 2 s = p = 0 p α p 1 ( 2 s 3 ) x 2 p = 0 p α p 1 ( 2 s 1 ) ,
we also define
Σ R ( b ) B [ x 6 lim s 0 p = 0 p 3 α p 2 s + x 6 lim s 0 p = 0 p 4 α p 4 ( 2 s 3 ) + x 4 2 lim s 0 p = 0 p 4 α p 4 ( 2 s 1 ) + 3 8 x 2 lim s 0 p = 0 p 4 α p 4 ( 2 s + 1 ) ] .
In view of (44) and (47), the first limit in (48) gives the following contribution to ζ ( 0 ) :
δ 2 = B 2 B ˜ 4 4 = B 240 = 1 120 .
Note that the second sum in (47) does not contribute to δ 2 because its only constant term contains Γ ( s + 1 ) Γ ( s ) , which tends to 0 as s 0 . The other limits in (48) do not contribute to ζ ( 0 ) in view of (21), because they only yield terms proportional to x 2 .
Last, the sum Σ ( c ) in (43) has the following asymptotic behavior as x :
Σ ( c ) ( 3 A + 9 B ) x 6 + k = 0 ( A C k + B D k ) x 7 k ,
which yields the following contribution to ζ ( 0 ) :
δ 3 = ( 3 A + 9 B ) 2 = 3 .
To sum up, we find
ζ I a ( 0 ) = δ 1 + δ 2 + δ 3 = 1 3 1 120 3 ,
Therefore, the full ζ ( 0 ) for physical degrees of freedom is given by (cf. (41) and (42))
ζ ( 0 ) = ζ I a ( 0 ) + 5 2 + 1199 360 = 112 45 .

3. First Example of Mixed Boundary Conditions on the Whole Set of Metric Perturbations and Ghost Modes

The previous example is very instructive, but of course it would be desirable to compute the effect of boundary conditions on the whole set of metric perturbations and Feynman–DeWitt–Faddeev–Popov ghost fields [25,26,27]. For this purpose, the work in Ref. [28] studied the following one-parameter family of mixed boundary conditions (with λ being a freely specifiable real parameter):
h i j τ + λ τ h i j M = 0 ,
[ h 0 i ] M = 0 ,
[ h 00 ] M = 0 ,
φ i τ + λ τ φ i M = 0 ,
φ 0 τ + ( λ + 1 ) τ φ 0 M = 0 .
With our notation, τ lies in the closed interval [ 0 , a ] ; h i j , h 0 i , h 00 are the components of metric perturbations; and φ i and φ 0 are covariant components of the ghost field of quantum gravity. One therefore deals with transverse-traceless modes, scalar modes, vector modes, decoupled scalar modes, decoupled vector modes, scalar ghost modes, vector ghost modes, and decoupled ghost modes.
A one-parameter family of full ζ ( 0 ) values is therefore obtained [28]:
ζ λ ( 0 ) = 89 90 + λ 3 ( λ 2 9 λ 3 ) .
The λ -dependent part of (59) is always positive, either for all
λ > 9 + 93 2 ,
or for all
λ 9 93 2 , 0 .
Equations (60) and (61) are sufficient conditions for the positivity of the full ζ λ ( 0 ) , and other suitable values of λ can be computed numerically.
This model is more complete than the one in Section 2, since it deals with all perturbative modes in the one-loop functional integral. However, it still suffers from a non-trivial drawback: the whole set of boundary conditions (54)–(58) is not completely invariant under infinitesimal diffeomorphisms on metric perturbations. For this reason, we resort to the boundary conditions of Section 4.

4. Completely Diff-Invariant Boundary Conditions

The boundary conditions that we study are part of a unified scheme for Maxwell, Yang–Mills, and Quantized General Relativity at one loop, i.e., [29,30]
[ π A ] B = 0 ,
[ Φ ( A ) ] B = 0 ,
[ φ ] B = 0 .
With our notation, π is a projector acting on the gauge field A , Φ is the gauge-fixing functional, and φ is the full set of ghost fields. Both Equations (62) and (63) are preserved under infinitesimal gauge transformations provided that the ghost obeys homogeneous Dirichlet conditions as in (64). For gravity, we choose Φ so as to have an operator P of Laplace type on metric perturbations in the one-loop Euclidean theory.

5. Eigenvalue Conditions for Scalar Modes

On the Euclidean four-ball, we expand metric perturbations h μ ν in terms of scalar, transverse vector, and transverse-traceless tensor harmonics on S 3 . For the vector, tensor, and ghost modes, boundary conditions reduce to Dirichlet or Robin. For scalar modes, one finds eventually the eigenvalues E = x 2 from the roots x of [31,32]
J n ( x ) ± n x J n ( x ) = 0 ,
J n ( x ) + x 2 ± n x J n ( x ) = 0 .
Note that both x and x solve the same equation.

6. Four Spectral ζ -Functions for Scalar Modes

From Equations (65) and (66), we obtain the following integral representations of the resulting ζ -functions upon exploiting the Cauchy theorem and rotation of contour:
ζ A , B ± ( s ) ( sin π s ) π n = 3 n ( 2 s 2 ) 0 d z z 2 s z log F A , B ± ( z n ) ,
where (here β + n , β n + 2 )
F A ± ( z n ) z β ± ( z n I n ( z n ) ± n I n ( z n ) ) ,
F B ± ( z n ) z β ± z n I n ( z n ) + ( z n ) 2 2 ± n I n ( z n ) ,
with I n being the modified Bessel functions of the first kind. Regularity at the origin is easily proved in the elliptic sectors, corresponding to ζ A ± ( s ) and ζ B ( s ) .

7. Regularity of ζ B + at s = 0

We now define T ( 1 + z 2 ) 1 / 2 and consider the uniform asymptotic expansion (away from T = 1 )
z β + F B + ( z n ) e n η ( T ) h ( n ) T ( 1 T 2 ) T 1 + j = 1 r j , + ( T ) n j ,
the functions r j , + being obtained from the Olver polynomials for the uniform asymptotic expansion of I n and I n . On splitting 0 1 d T = 0 μ d T + μ 1 d T with small μ , we obtain an asymptotic expansion of the l.h.s. by writing, in the first interval on the r.h.s.,
log 1 + j = 1 r j , + ( T ) n j j = 1 R j , + ( T ) n j ,
and then computing
C j ( τ ) R j , + T = ( 1 T ) j 1 a = j 1 4 j K a ( j ) T a .
The integral μ 1 d T is instead found to yield a vanishing contribution in the μ 1 limit. Remarkably, by virtue of the spectral identity
g ( j ) a = j 4 j Γ ( a + 1 ) Γ ( a j + 1 ) K a ( j ) = 0 ,
which holds j = 1 , , , we find
lim s 0 s ζ B + ( s ) = 1 6 a = 3 12 a ( a 1 ) ( a 2 ) K a ( 3 ) = 0 ,
and
ζ B + ( 0 ) = 5 4 + 1079 240 1 2 a = 2 12 ω ( a ) K a ( 3 ) + j = 1 f ( j ) g ( j ) = 296 45 ,
where
ω ( a ) 1 6 Γ ( a + 1 ) Γ ( a 2 ) [ log ( 2 ) ( 6 a 2 9 a + 1 ) 4 Γ ( a 2 ) Γ ( a + 1 ) + 2 ψ ( a + 1 ) ψ ( a 2 ) ψ ( 4 ) ] ,
f ( j ) ( 1 ) j j ! [ 1 2 2 j + ζ R ( j 2 ) ( 1 δ j , 3 ) + γ δ j , 3 ] .
The spectral cancellation (73) achieves three goals: (i) vanishing of the log 2 coefficient in Equation (75); (ii) vanishing of j = 1 f ( j ) g ( j ) in Equation (75); and (iii) regularity at the origin of ζ B + .
To cross-check our analysis, we evaluate r j , + ( T ) r j , ( T ) and hence obtain R j , + ( T ) R j , ( T ) for all j. Only j = 3 contributes to ζ B ± ( 0 ) , and we find
ζ B + ( 0 ) = ζ B ( 0 ) 1 24 l = 1 4 Γ ( l + 1 ) Γ ( l 2 ) ψ ( l + 2 ) 1 ( l + 1 ) κ 2 l + 1 ( 3 ) = 206 45 + 2 = 296 45 ,
in agreement with Equation (75), where κ 2 l + 1 ( 3 ) are the four coefficients on the right-hand side of
T ( R 3 , + R 3 , ) = ( 1 T 2 ) 4 ( 80 T 3 24 T 5 + 32 T 7 8 T 9 ) .
Within this framework, the spectral cancellation reads as
l = 1 4 Γ ( l + 1 ) Γ ( l 2 ) κ 2 l + 1 ( 3 ) = 0 ,
which is a particular case of
a = a min ( j ) a = a max ( j ) Γ ( ( a + 1 ) / 2 ) Γ ( ( a + 1 ) / 2 j ) κ a ( j ) = 0 .
Interestingly, the full ζ ( 0 ) value for pure gravity (i.e., including the contribution of tensor, vector, scalar, and ghost modes) is then found to be positive [31,32]:
ζ ( 0 ) = 142 45 ,
which suggests in light of (5) a quantum avoidance of the cosmological singularity driven by full diffeomorphism invariance of the boundary-value problem for one-loop quantum theory.

8. Open Problems

The DeWitt boundary condition lies at the very heart of deep issues in quantum gravity. As far as we can see, the main open problems are as follows.
(1)
Among the three schemes studied in our Section 2, Section 3, Section 4, Section 5, Section 6, Section 7, the latter, i.e., the choice of completely diff-invariant boundary conditions on all perturbative modes, might seem the most satisfactory, but unfortunately, the strong ellipticity of the boundary-value problem is not fulfilled in such a case [30,33,34,35,36,37]. However, our analysis shows that, in the particular case of flat Euclidean four-space bounded by a three-sphere boundary, peculiar cancellations occur, and the resulting ζ ( 0 ) value can be defined and is positive. The deeper underlying reason might be that, in order to define a spectral ζ -function, it is sufficient to find a sector of the complex plane free of eigenvalues of the leading symbol of the elliptic operator under consideration (we are grateful to Professor Gerd Grubb for correspondence about this property a long time ago). An alternative approach might consist in considering non-local boundary conditions in Euclidean quantum gravity [38,39,40], or the normalizability criterion for the wave function of the universe [41].
(2)
The outstanding work in Ref. [10] looked for solutions of the quantum constraint equations in order to check the validity of DeWitt’s proposal. However, although one can obtain under suitable assumptions a formal proof of the equivalence of canonical and functional-integral approaches [42], DeWitt himself provided an enlightening example of a sum over histories that does not solve the Wheeler–DeWitt equation [43]. This remark might therefore account for the inequivalence between our conclusions and the results in Ref. [10].
The fascinating question of whether our universe can be non-singular in a semiclassical theory of quantum gravity [44] is therefore still waiting for a fully satisfactory answer.

Funding

This research received no external funding.

Data Availability Statement

All calculations were done by hand, hence there are no data stored.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. The One-Loop Approximation

We are here interested in the approach to quantum field theory in terms of Feynman functional integrals. Hence, we study the amplitudes of going from data on a spacelike surface Σ 1 to data on a spacelike surface Σ 2 . For example, in the case of real scalar fields ϕ in a curved background M, the data are the induced three-metric h and a linear combination of ϕ and its normal derivative: a ϕ + b ϕ n . The latter reduces to homogeneous Dirichlet conditions if b = 0 , and Neumann conditions if a = 0 . Otherwise, it is a Robin boundary condition. The quantum amplitudes are functionals of these boundary data. On making a Wick rotation and using the background-field method, we may expand both the four-metric g and the field ϕ about solutions to the classical field equations as g = g 0 + g ¯ and ϕ = ϕ 0 + ϕ ¯ . However, the more general possibility remains to consider background fields that are not solutions to any field equations, or which are (approximate) solutions in the asymptotic regions. The logarithm of the one-loop functional integral Z for a scalar field (in the main body of our paper, we study pure gravity, but here we focus on scalar fields for simplicity) has an asymptotic expansion
log ( Z ) log μ [ ϕ ¯ ] e I 2 [ ϕ ¯ ] / + O ( 2 ) ,
where μ is a suitable measure of the space of scalar-field perturbations. The part I 2 [ ϕ ¯ ] of the action that is quadratic in scalar-field perturbations involves a second-order elliptic operator B . Assuming completeness of the set { φ n } of eigenfunctions of B , with eigenvalues λ n , the corresponding contribution to one-loop quantum amplitudes involves an infinite product of Gaussian integrals, i.e.,
n = n 0 μ d y n e λ n 2 y n 2 = 1 det ( 1 2 π 1 μ 2 B ) .
In order to make sense of this infinite product of eigenvalues, one can use ζ -function regularization. This is a rigorous mathematical tool that relies on the spectral theorem, according to which for any elliptic, self-adjoint, and positive-definite operator B, its complex powers B s can be defined. Hence, its spectral ζ -function is defined as in Equation (4), and the analytic continuation of the ζ -function to the whole complex-s plane takes the form
ζ B ( s ) = k = n N a k s + k m + ϕ N ( s ) , k 0 .
Thus, on using analytic continuations, ζ B ( 0 ) is actually finite, and its value gives information about scaling properties of quantum amplitudes. We can now be more precise and describe in detail some key properties. The relation
det B = e ζ ( 0 )
becomes a possible way to define the determinant of the elliptic operator B upon the analytic continuation of ζ B ( s ) . If B is a second-order operator, its eigenvalues λ n have dimension ( length ) 2 . Under conformal rescaling of the metric according to g ^ = k 2 g , one has λ ^ n = λ n / k 2 , and the new spectral ζ -function is ζ ^ ( s ) = k 2 s ζ ( s ) . This leads to
log det B ^ = log det B log ( k 2 ) ζ ( 0 ) ,
and hence the partition function scales as
log ( Z ^ ) = log ( Z ) + 1 2 log ( k 2 ) ζ ( 0 ) + log ( μ ^ / μ ) ζ ( 0 ) .
The parameter μ is the one occurring in the one-loop semiclassical evaluation of the functional integral. This formula allows for the more general case when the normalization parameter μ changes under scale transformations. One can avoid this complication by assuming that the measure in the functional integral is defined on scalar densities of weight 1 2 .
Equation (A5) can also be used to deduce that the one-loop effective action (for the scalar field) reads as
Γ ( 1 ) = 1 2 log det B ^ = 1 2 ζ ( 0 ) 1 2 ζ ( 0 ) log ( k 2 ) .
Note that the resulting one-loop out in amplitude is measure-dependent unless ζ ( 0 ) = 0 . This is why ζ ( 0 ) is frequently called the anomalous scaling factor.

References

  1. Friedmann, A. Über die krümmung des raumes. Z. Phys. 1922, 10, 377. [Google Scholar] [CrossRef]
  2. Penrose, R. Gravitational collapse and space-time singularities. Phys. Rev. Lett. 1965, 14, 57. [Google Scholar] [CrossRef] [Green Version]
  3. Hawking, S.W. The occurrence of singularities in cosmology. Proc. R. Soc. Lond. A 1966, 294, 511. [Google Scholar]
  4. Hawking, S.W. The occurrence of singularities in cosmology. II. Proc. R. Soc. Lond. A 1966, 295, 490. [Google Scholar]
  5. Geroch, R.P. Singularities in closed universes. Phys. Rev. Lett. 1966, 17, 445. [Google Scholar] [CrossRef]
  6. Hawking, S.W. The occurrence of singularities in cosmology. III. Causality and singularities. Proc. R. Soc. Lond. A 1967, 300, 187. [Google Scholar]
  7. Hawking, S.W.; Penrose, R. The singularities of gravitational collapse and cosmology. Proc. R. Soc. Lond. A 1970, 314, 529. [Google Scholar]
  8. Christodoulou, D.; Klainerman, S. The Global Nonlinear Stability of Minkowski Space; Princeton University Press: Princeton, NJ, USA, 1993. [Google Scholar]
  9. DeWitt, B.S. Quantum theory of gravity. I. The canonical theory. Phys. Rev. 1967, 160, 1113. [Google Scholar] [CrossRef] [Green Version]
  10. Matsui, H.; Mukoyama, S.; Naruko, A. DeWitt boundary condition is consistent in Horava-Lifshitz gravity. Phys. Lett. B 2022, 833, 137340. [Google Scholar] [CrossRef]
  11. Kamenshchik, A.Y. Quantum cosmology and late-time singularities. Class. Quantum Grav. 2013, 30, 173001. [Google Scholar] [CrossRef]
  12. Schleich, K. Semiclassical wave function of the universe at small three-geometries. Phys. Rev. D 1985, 32, 1889. [Google Scholar] [CrossRef] [PubMed]
  13. Esposito, G. Quantum Gravity, Quantum Cosmology and Lorentzian Geometries; Lecture Notes in Physics Monographs; Springer: Berlin/Heidelberg, Germany, 1994; Volume 12. [Google Scholar]
  14. Breitenlohner, P.; Freedman, D.Z. Stability in gauged extended supergravity. Ann. Phys. 1982, 144, 249. [Google Scholar] [CrossRef]
  15. Hawking, S.W. The boundary conditions for gauged supergravity. Phys. Lett. B 1983, 126, 175. [Google Scholar] [CrossRef]
  16. Moss, I.G. Boundary terms in the heat kernel expansion. Class. Quantum Grav. 1989, 6, 759. [Google Scholar] [CrossRef]
  17. Watson, G.N. A Treatise on the Theory of Bessel Functions; Cambridge University Press: Cambridge, UK, 1966. [Google Scholar]
  18. Wong, R. Asymptotic Approximations of Integrals; Academic Press: New York, NY, USA, 1989. [Google Scholar]
  19. Greiner, P. An asymptotic expansion for the heat equation. Arch. Rat. Mech. Anal. 1971, 41, 163. [Google Scholar] [CrossRef]
  20. Branson, T.P.; Gilkey, P.B. The asymptotics of the Laplacian on a manifold with boundary. Commun. Part. Diff. Eqs. 1990, 15, 245. [Google Scholar] [CrossRef] [Green Version]
  21. Kennedy, G. Boundary terms in the Schwinger-DeWitt expansion: Flat-space results. J. Phys. A 1978, 11, L173. [Google Scholar] [CrossRef]
  22. Kennedy, G. Some Finite-Temperature Quantum-Field-Theory Calculations in Manifolds with Boundaries. Ph.D. Thesis, University of Manchester, Manchester, UK, 1979. [Google Scholar]
  23. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions; Dover: New York, NY, USA, 1964. [Google Scholar]
  24. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series and Products; Academic Press: New York, NY, USA, 1965. [Google Scholar]
  25. Feynman, R.P. Quantum theory of gravitation. Acta Phys. Pol. 1963, 24, 697. [Google Scholar]
  26. DeWitt, B.S. Quantum theory of gravity. II. The manifestly covariant theory. Phys. Rev. 1967, 162, 1195. [Google Scholar] [CrossRef]
  27. Faddeev, L.D.; Popov, V. Feynman diagrams for the Yang-Mills field. Phys. Lett. B 1967, 25, 29. [Google Scholar] [CrossRef]
  28. Esposito, G.; Kamenshchik, A.Y. Mixed boundary conditions in Euclidean quantum gravity. Class. Quantum Grav. 1995, 12, 2715. [Google Scholar] [CrossRef]
  29. Esposito, G.; Kamenshchik, A.Y.; Pollifrone, G. Euclidean Quantum Gravity on Manifolds with Boundary; Fundamental Theories of Physics; Kluwer: Dordrecht, The Netherlands, 1997; Volume 85. [Google Scholar]
  30. Avramidi, I.G.; Esposito, G. Gauge theories on manifolds with boundary. Commun. Math. Phys. 1999, 200, 495. [Google Scholar] [CrossRef] [Green Version]
  31. Esposito, G.; Fucci, G.; Kamenshchik, A.Y.; Kirsten, K. Spectral asymptotics of Euclidean quantum gravity with diff-invariant boundary conditions. Class. Quantum Grav. 2005, 22, 957. [Google Scholar] [CrossRef] [Green Version]
  32. Esposito, G.; Fucci, G.; Kamenshchik, A.Y.; Kirsten, K. A non-singular one-loop wave function of the universe from a new eigenvalue asymptotics in quantum gravity. J. High Energy Phys. 2005, 9, 063. [Google Scholar] [CrossRef] [Green Version]
  33. Avramidi, I.G.; Esposito, G. Lack of strong ellipticity in Euclidean quantum gravity. Class. Quantum Grav. 1998, 15, 1141. [Google Scholar] [CrossRef] [Green Version]
  34. Witten, E. A note on boundary conditions in Euclidean gravity. Rev. Math. Phys. 2021, 33, 2140004. [Google Scholar] [CrossRef]
  35. Witten, E. A note on the canonical formalism for gravity. arXiv 2022, arXiv:2212.08270. [Google Scholar]
  36. Park, I.Y. Holographic quantization of gravity in a black hole background. J. Math. Phys. 2016, 57, 022305. [Google Scholar] [CrossRef] [Green Version]
  37. Park, I.Y. Foliation, jet bundle and quantization of Einstein gravity. Front. Phys. 2016, 4, 25. [Google Scholar] [CrossRef] [Green Version]
  38. Esposito, G. Non-local boundary conditions in Euclidean quantum gravity. Class. Quantum Grav. 1999, 16, 1113. [Google Scholar] [CrossRef]
  39. Esposito, G. New kernels in quantum gravity. Class. Quantum Grav. 1999, 16, 3999. [Google Scholar] [CrossRef] [Green Version]
  40. Schröder, M. On the Laplace operator with non-local boundary conditions. Rep. Math. Phys. 1989, 27, 259. [Google Scholar] [CrossRef]
  41. Barvinsky, A.O.; Kamenshchik, A.Y. One-loop quantum cosmology: The normalizability of the Hartle-Hawking wave function and the probability of inflation. Class. Quantum Grav. 1990, 7, L181. [Google Scholar] [CrossRef]
  42. Halliwell, J.J.; Hartle, J.B. Wave functions constructed from an invariant sum over histories satisfy constraints. Phys. Rev. D 1991, 43, 1170. [Google Scholar] [CrossRef] [PubMed]
  43. DeWitt, B.S. The quantum and gravity: The Wheeler-DeWitt equation. In Relativity, Particle Physics and Cosmology, Proceedings of the Richard Arnowitt Fest, College Station, TX, USA, 5–8 April 1998; World Scientific: Singapore, 1999; pp. 70–92. [Google Scholar]
  44. Hawking, S.W. A non-singular universe. Phys. Scripta T 2005, 117, 49. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Esposito, G. DeWitt Boundary Condition in One-Loop Quantum Cosmology. Universe 2023, 9, 187. https://doi.org/10.3390/universe9040187

AMA Style

Esposito G. DeWitt Boundary Condition in One-Loop Quantum Cosmology. Universe. 2023; 9(4):187. https://doi.org/10.3390/universe9040187

Chicago/Turabian Style

Esposito, Giampiero. 2023. "DeWitt Boundary Condition in One-Loop Quantum Cosmology" Universe 9, no. 4: 187. https://doi.org/10.3390/universe9040187

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop