Resolution of curvature singularities from quantum mechanical and loop perspective

We analyze the persistence of curvature singularities when analyzed using quantum theory. First, quantum test particles obeying the Klein-Gordon and Chandrasekhar-Dirac equation are used to probe the classical timelike naked singularity. We show that the classical singularity is felt even by our quantum probes. Next, we use loop quantization to resolve singularity hidden beneath the horizon. The singularity is resolved in this case.


Introduction
One of the important predictions of the Einstein's theory of general relativity is the formation of spacetime singularities. In classical general relativity, singularities are defined as points in which the evolution of timelike or null geodesics is not defined after a finite proper time. According to the classification of the classical singularities devised by Ellis and Schmidt [1], scalar curvature singularities are the most strongest one in the sense that the spacetime posses incomplete geodesics ending in them and all the physical quantities such as the gravitational field (scalars formed from curvature tensor), energy density and tidal forces diverge at the singular point.
But such divergence of physical quantities signify the breakdown of predictive power of classical general relativity. If these singularities are covered by horizon (as supposed by Cosmic Censorship Conjecture) then at least the physically most relevant region of spacetime is a e-mail: tayabeh.tahamtan@emu.edu.tr b e-mail: ota@matfyz.cz under control. Naked singularities (those not covered by horizon), on the other hand, provide an observer with causal access to the region of diverging quantities and should be avoided. However, even singularities covered by the horizon can be accessed by an infalling observer and, more importantly, we would like to have a theory that lacks divergences, at least effectively.
The natural direction for resolving the problem of singularities in classical theory is investigating their persistence in quantum picture. Although we do not have a final quantum theory of gravity we still have several tools for analyzing quantum singularities. The first approach relies on examining properties of quantum particle wave functions on the background represented by the studied geometry. This is a frequently used technique based on well understood properties of operators on a Hilbert space. To move further, one might proceed to using quantum fields and possibly even the backreaction of background geometry using semiclassical Einstein equations with suitably regularized stress energy tensor. Finally, one can apply quantization of the geometry itself. The last approach is in principle the most precise but relies on the selected quantization method and we have no generally accepted one in case of gravity.
We will apply two of the above mentioned approaches for analysis of singularity in case of the general metric of global monopole [16], which is determined by two parameters -one characterizing the "Schwarzschild-type mass" and the other one the deficit of solid angle. The singularity is generally covered by single horizon but the class of metrics also contains, as a special case, a naked singularity which is analyzed from quantum mechanical point of view using the technique of Horowitz and Marolf [17] (who continued the pioneering work of Wald [18]). This method for analyzing timelike singularities is based on investigation of self-adjoint extensions of the evolution operator associated with the given wave equation. If it is unique the spacetime is deemed quantum mechanically non-singular. The analysis is carried out for relativistic quantum particle wave equations on a fixed background. Specifically, we review the previous results for Klein-Gordon equation and show the calculation using Newman-Penrose formalism for the Dirac equation, both in the case of pure global monopole with naked singularity for which the method was developed.
But as already mentioned, the most reliable method when trying to investigate the possible removal of the singularities from geometry is quantum gravity. Here we have selected loop quantization method inspired by [19,20,21], where the spacetime beneath the horizon (in the non-naked subclass) is isometric to the Kantowski-Sachs cosmology. Then one can apply methods from Loop Quantum Cosmology (LQC), that are based on loop quantization on the restricted configuration space. In this way, the results for resolution of initial cosmological singularity are translated to statements about the singularity at the origin r = 0.

The General Metric for Global monopole
It is well known that different types of non-standard topological objects may have been formed during initial Universe evolution, such as domain walls, cosmic strings and monopoles [16,22]. The basic idea is that these topological defects have formed as a result of a breakdown of local or global gauge symmetries. The simplest model that gives rise to global monopole is described by the Lagrangian where φ a is a triplet of scalar fields, a = 1, 2, 3. The model has a global O (3) symmetry, which is spontaneously broken to U (1). The field configuration describing the monopole is where x a x a = r 2 . We assume that underlying geometry is general static spherically symmetric described by the line element with the usual relation between the spherical coordinates, r, θ, φ and the Cartesian coordinates x a . The Lagrangian for the above given field configuration simplifies in the following way and the diagonal energy momentum tensor is given by these components The general solution of the Einstein equations with this T ν µ is where M is a constant of integration. The metric describes a black hole of mass M , carrying a global monopole charge characterized by η. Such a black hole can be formed if a global monopole is swallowed by an ordinary black hole [16]. The Kretschmann scalar which indicates the formation of curvature singularity is given by It is obvious that r = 0 is a typical central curvature singularity (scalar curvature singularity according to above mentioned classification) and the dominant contribution comes from term corresponding to black hole mass M . If M > 0 the singularity is evidently spacelike and covered by a single horizon.

Global monopole and its singularity
If we assume that the mass term is negligible on the astrophysical scale or vanishing, we will have For simplicity we choose α 2 = 1 − 8πGη 2 and by rescaling r and t variables, we can rewrite the monopole metric as If we calculate the Kretschmann scalar, still there is a weaker singularity at r = 0. From the metric (7) one can immediately see that the singularity is timelike. This time, because our simplified metric does not have the horizon the singularity is naked.

Naked Singularity
As mentioned in the Introduction naked singularity poses a serious problems and its resolution would be desirable. In this section, the occurrence of naked singularities in global monopole will be analyzed from quantum mechanical point of view. In probing the singularity, quantum test particles obeying the Klein-Gordon and Dirac equations are used. The reason for using two different types of fields is to clarify whether the classical singularity is sensitive to spin of the fields.
According to Horowitz and Marolf (HM) [17], the singular character of the spacetime is defined as the ambiguity in the evolution of the wave functions. That is to say, the singular character is determined based on the number of self-adjoint extensions of the evolution operator to the entire Hilbert space. If the extension is unique, it is said that the spacetime is quantum mechanically regular. The brief review of the method follows: Consider a static spacetime (M, g µν ) with a timelike Killing vector field ξ µ . Let t denote the Killing parameter and Σ denote a static slice. The Klein-Gordon equation in this space is This equation can be written in the form in which f = −ξ µ ξ µ and D i is the spatial covariant derivative on Σ. We assume that the Hilbert space H = L 2 (Σ, µ) is the space of square integrable functions on Σ with appropriate measure µ. Initially the operator A is defined on smooth functions with compact support C ∞ 0 (Σ). Since the operator A is real, positive and symmetric its self-adjoint extensions always exist. If it has a unique extension A E , then A is called essentially self-adjoint [23,24,25]. Accordingly, the Klein-Gordon equation for a free particle satisfies with the solution If A is not essentially self-adjoint, the future time evolution of the wave function (12) is ambiguous. Then, HM criterion defines the spacetime as quantum mechanically singular. However, if there is only a single selfadjoint extension, the operator A is said to be essentially self-adjoint and the quantum evolution described by equation (12) is uniquely determined by the initial conditions. According to the HM criterion, this spacetime is said to be quantum mechanically non-singular.
In order to determine the number of self-adjoint extensions, the concept of deficiency indices is used. The deficiency subspaces N ± are defined by ( see Ref. [26] for a detailed mathematical background), that belong to the Hilbert space H. If there are no square integrable solutions ( i.e. n + = n − = 0), the operator A possesses a unique self-adjoint extension and it is essentially self-adjoint. Consequently, a sufficient condition for the operator A to be essentially self-adjoint is to find only solutions satisfying Eq. (14) that do not belong to the Hilbert space.

Klein -Gordon Fields
The Klein-Gordon equation for a massless scalar particle is given by, For the metric (8), the Klein-Gordon equation becomes, In analogy with the equation (10), the spatial operator A is and the equation to be solved is (A * ± i) ψ = 0. Using separation of variables, ψ = R (r) Y m l (θ, ϕ), we get the radial portion of equation (14) as, The square integrability of the above solution is checked by calculating the squared norm of the above solution in which the function space on each t = constant hypersurface Σ is defined as H =L 2 (Σ, µ) where µ is the measure given by the spatial metric volume element.
We easily recover the results showed in [8]: The spacetime of global monopole remains singular in the view of relativistic quantum mechanics: the future of a given initial wave packet obeying the Klein-Gordon equation is not generally well determined, similarly to the future of a classical particle which reaches the classical singularity at r = 0.

Dirac Fields
The Newman-Penrose formalism will be used here to analyze massless Dirac particle propagating in the space of global monopole. The signature of the metric (8) is changed to −2 in order to use the Dirac equation in Newman-Penrose formalism. Thus, the metric is given by, The Chandrasekhar-Dirac (CD) [10] equations in Newman-Penrose formalism are given by where F 1 , F 2 , G 1 and G 2 are the components of the wave function, ǫ, ρ, π, α, µ, γ, β and τ are the spin coefficients to be found and the "bar" denotes complex conjugation. The null tetrad vectors for the metric (19) are defined by The directional derivatives in the Dirac equation are defined by D = l a ∂ a , ∇ = n a ∂ a and δ = m a ∂ a . We define operators in the following way The nonzero spin coefficients are, Substituting nonzero spin coefficients and the definitions of the operators given above into the CD equations leads to For the solution of the CD equations, we assume separable solution in the form of Here {f 1 , f 2 , g 1 , g 2 } and {Y 1 , Y 2 , Y 3 , Y 4 } are functions of r and θ respectively, m is the azimuthal quantum number and k is the frequency of the Dirac spinor, which is assumed to be positive and real. By substituting (25) in (24) we will see that with these assumptions f 1 (r) = g 2 (r) and f 2 (r) = g 1 (r) , Dirac equation reduces to two equations. The radial part of the Dirac equations become where λ comes from separation of variables. We further assume that then equation (28) transforms into, where λ ′ = λ α , so we will have d dr In order to write the above equation in a more compact form we combine the solutions in the following way, After doing some calculations we end up with a pair of one-dimensional Schrödinger-like wave equations with effective potentials, In analogy with the equation (10), the spatial operator A for the massless case is The solutions of the above equations are expressible using Bessel functions of the first and second kind in the following way Using the asymptotic formulas for Bessel functions when r → ∞ (Y (κ, z) ≈ z −1/2 sin(z−κπ/2−π/4) and J(κ, z) ≈ z −1/2 cos(z − κπ/2 − π/4)) and noting the complex argument in both solutions one can find a combination of constants C 1 , C 2 or C ′ 1 , C ′ 2 which is square integrable near infinity. (But, it is also possible to choose the constants differently so that both solutions are not square integrable!).
When r → 0 the approximate expressions for Bessel functions (Y (κ, z) ≈ z −κ for κ = 0, Y (0, z) ≈ ln(z/2) and J(κ, z) ≈ z κ ) imply that for C 2 = 0 and C ′ 2 = 0 we have square integrable solution near zero. (Here again if we suppose C 1 = 0 and C ′ 1 = 0, for κ ≥ 3/2, the solutions are not square integrable!. One could restrict an analysis to only certain wave modes and purposely choose the modes to be quantum regular).
But since we have a solution of equations valid on the whole domain (not just asymptotic forms of equations) we can match the behaviour at zero and infinity. Based on the results we can have solution square integrable over the whole domain and therefore our deficiency indices are nonzero. The operator is not essentially self-adjoint and the spacetime is quantum mechanically singular.

Quantum Gravity
Now we are going to investigate the singularity of general global monopole using techniques from loop quantization in the manner of [20]. Consider equation (2), for r < 2GM 1−8πGη 2 . This metric describes spacetime inside the horizon of a black hole. The coordinate r is timelike and the coordinate t is spatial there; for convenience we rename them as r ≡ T and t ≡ r with T ∈ [0, 2GM 1−8πGη 2 ] and r ∈ [−∞, +∞] and the metric becomes we eliminate the coefficient of dT 2 by defining a new temporal variable τ via Accordingly, the metric becomes We introduce two functions a 2 (τ ) ≡ 2Gm T − α 2 and b 2 (τ ) ≡ T 2 (τ ) and redefine τ ≡ t. The metric becomes this metric describes a homogeneous, anisotropic Kantowski-Sachs cosmological mmodel with spatial section having topology R × S 2 . From this observation comes the motivation to use LQC approach. In our case a (t) is a function of b (t).

Classical observables
The corresponding action for gravity minimally coupled with scalar field can be written in the form by considering the metric (38), the action becomes by using the relation between a and b, we will be able to write the action in terms of a single function Now, we will compute the Hamiltonian (Hamiltonian constraint). The momentum associated to the chosen configuration variable is and therefore we obtain

Now, we calculate the Hamiltonian constraint in terms ofḃ
and immediately get the following solutioṅ which is exactly the equation (36). When the horizon radius, r h = 2GM α 2 , is much larger than the scale on which we are probing the singularity, we can write so the Hamiltonian would be The volume simplifies when using the above approximation and we obtain The canonical pair is given by b ≡ x and p b , with Poisson bracket {x, p b } = 1.
For isotropic models, only holonomies evaluated in isotropic connections A i a =cδ i a appear. Along straight lines in the direction of translation symmetries X a I = ∂/∂X I a , holonomies exp X a I A i a τ i in the fundamental representation of SU (2) have matrix elements of the form exp (iµc), where µ depends on the length of the curve used. Here, it turns out to be useful to introduce c := V 1/3 0c defined in terms of the coordinate size V 0 of the region used to define the isotropic phase space [27].
Using this motivation we introduce the following function which will be used instead of the momentum (from now on we leave out the subscript b for momentum associated with this observable) [20] where γ is a real parameter and L fixes the length scale. The parameter γ determines the separation of momentum points in the phase space. The pair (x, U γ (p)) has the following Poisson bracket algebra A straightforward calculation gives We are concerned with the quantity 1 |x| which can serve as an indicator for singularity presence because classically it diverges for |x| → 0 thus producing singularity. From this moment we choose n = 1/3

Quantization
We will use the basis of Hillbert space introduced in [19,20], which is formed by eigenstates ofx. This implies the existence of a self-adjoint operatorx, acting on the basis states according tô Next, we want to promote the classical momentum function U γ = e (8πG iγ L p) to operator. We can do so by defining the action ofÛ γ on the basis states with the help of the definition equation (53) and using commutation relation based on the Poisson bracket between x and U γ we obtain

Volume operator and disappearance of the singularity
In the vicinity of the singularity we assume the approximate equation (48). Then the volume operator acts in the following way on the basis stateŝ Using the equation (52) and promoting the Poisson brackets to commutators, while setting γ = 1, we find On the basis states this operator acts in the following waŷ (58) so finally we get We can see that the spectrum is bounded from above and so the singularity is resolved in the quantum theory (the theory gives finite predictions for observables related to singularity). In fact, the eigenvalue of operator 1 |x| corresponding to the state |0 which probes the classical singularity is equal to 2 πl 2 p , which is the highest eigenvalue of the spectrum. Specifically, the operator corresponding to the curvature invariant R µνρσ R µνρσ = 48M 2 G 2 r 6 + 128M πG 2 η 2 r 5 + 256G 2 π 2 η 4 r 4 is then automatically finite in quantum mechanics. Promoting it to operator and evaluating on |0 we get On the other hand, when |µ| → ∞ the eigenvalue of 1 |x| goes to zero which is natural behaviour for large |x|. Also, it is possible to show that the quantum Hamiltonian constraint gives a discrete difference equation for the coefficients of the physical states.

Conclusion
We have seen that we have not been successful in removing the naked singularity by using relativistic quantum mechanics (for both Klein-Gordon and Dirac equations). On the other hand we have shown that the curvature singularity of general global monopole is resolved when the geometry is quantized using loop techniques. Unfortunately, one cannot directly compare the results because the loop quantization relied on radial coordinate being timelike beneath the horizon which is not the case for naked singularity of pure monopole. But still, this might be an indication that the first method is not reliable for determining the fate of singularities in quantum theory and one should rather focus on quantization of the geometry itself. But even the approach using loop quantization that relied on restricted class of geometries should not be trusted completely. One should allow, e.g., for deviations from spherical symmetry to be completely sure about the fate of singularities.