2+1-dimensional gravity coupled to a dust shell: quantization in terms of global phase space variables
Alexander Andrianov, Yasser Elmahalawy, Artem Starodubtsev

TL;DR
This paper analyzes a 2+1-dimensional gravity model with a dust shell, identifying global phase space variables with ADS^{2} geometry, leading to non-commutative and discrete coordinates that resolve singularities.
Contribution
It introduces a global phase space parametrization for the model using ADS^{2} geometry and Euler angles, connecting to Kuchar variables for potential 3+1 dimensional extension.
Findings
Global phase space is an ADS^{2} manifold.
Coordinate space exhibits non-commutativity and discreteness.
Singularity is resolved through this geometric and algebraic structure.
Abstract
We perform canonical analysis of a model in which gravity is coupled to a spherically symmetric dust shell in 2+1 spacetime dimensions. The result is a reduced action depending on a finite number of degrees of freedom. The emphasis is made on finding canonical variables providing the global chart for the entire phase space of the model. It turns out that all the distinct pieces of momentum space could be assembled into a single manifold which has ADS^{2}-geometry, and the global chart for it is provided by the Euler angles. This results in both non-commutativity and discreteness in coordinate space, which allows to resolve the central singularity. We also find the map between ADS^{2} momentum space obtained here and momentum space in Kuchar variables, which could be helpful in extending the present results to 3+1 dimensions.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
0
Alexander Andrianov, Yasser Elmahalawy and Artem Starodubtsev 111St.Petersburg state university, St.Petersburg, Russia, email: [email protected], [email protected], [email protected]
**2+1-dimensional gravity coupled to a dust shell: quantization in terms of global phase space variables. **
Abstract
We perform canonical analysis of a model in which gravity is coupled to a spherically symmetric dust shell in 2+1 spacetime dimensions. The result is a reduced action depending on a finite number of degrees of freedom. The emphasis is made on finding canonical variables providing the global chart for the entire phase space of the model. It turns out that all the distinct pieces of momentum space could be assembled into a single manifold which has -geometry, and the global chart for it is provided by the Euler angles. This results in both non-commutativity and discreteness in coordinate space, which allows to resolve the central singularity. We also find the map between momentum space obtained here and momentum space in Kuchar variables, which could be helpful in extending the present results to 3+1 dimensions.
1 Introduction
It has long been believed that quantum theory of gravity could regularize its singularities. The first argument for that dates back to Bronstein [1], who showed that there could be no measurable distances smaller than the Planck length. In the absence of a full theory of gravity there is still no decisive answer whether this is indeed the case or not.
One can address this question by studying reduced models for General Relativity in which all but a few degrees of freedom are removed. Of special interest are those which contain black hole solutions, which are essential for Bronstein’s argument.
The simplest model possible is a homogenous universe with a matter field. Quantization of this model has been extensively studied, and the overall conclusion is that quantum theory does’t cure the singularity in this case [2], unless we include some exotic matter[3].
The second simplest model is spherically symmetric spacetime in which matter is represented by one or more dust shells. Unlike homogenous model, it accomodates black hole like solutions, which results in a non-trivial phase space with branching of the solution to the constraints.
There is a variety of works studying such models both on classical [4, 5, 6] and quantum [7, 9, 8] level. In some of the versions of quantum theory the central singularity is removed [9, 8]. However the above results do not always agree with each other. Apart from quantization ambiguity, the other possible reason for that is a complicated phase space structure of the model. Different quantum theories could arise on different sectors of such phase space.
In such situation the common wisdom is that the wavefunction of a quantum theory has to be defined on all possible configurations, independently of whether they are classically accessible or not. In a particular way it was realized in [9] where by making use of complex coordinates different sectors of the phase space were assembled into one Riemann surface, the branching point identified with a horizon.
Another possibility is to try to find a real global chart for the phase space (if it exists). One example where it was possible to realize is 2+1 dimensional gravity coupled to a point particle [10]. The momentum of the particle turn out to live on the Lorentz group manifold, and the different branches of the solution to the constraints result from different way of intersecting this manifold by a plane.
This will be our starting point in attempting to relate the two above approaches.
In section 2 we repeat the canonical analysis [9, 6] for gravity+shell system, including Kuchar canonical transformation necessary for it, but now in 2+1 spacetime dimensions. The result is very similar to that in 3+1 dimensions, the only difference being that there is no Newtonian potential, but the branching of the constraint solution is the same.
In section 3 we extend the results of [10] and [11] from a point particle to a spherical shell, representing the later as an ensemble of infinite number of point particles. The momentum space of the shell turn out to be , which results, in particular, in non-commutativity of the coordinates. The Hamiltonian constraint is, however, different from that of a particle, which accounts for a contribution of inter-shell movement energy to the gravitational field. The relation between momenta from and canonical momenta from [9, 6] is found.
In section 4 we perform quantization in momentum representation on . Apart from non-commutativity of the coordinates it results in discreteness of the spectrum of one of them (time). The areal radius of the shell has discrete spectrum when it is timelike and continuous, but separated from zero spectrum when it is spacelike. The later provides a mechanism for resolution of the central singularity.
Finally, the implementation of the Hamiltonian constraint and finding the physical Hilbert space is discussed. We also discuss a possibility to extend the above results to 3+1 spacetime dimensions.
2 Canonical analysis of gravity+shell action and the problem of finding the global phase space variables
In this section we rederive the results of [9, 6] in space one dimension lower.
2.1 Canonical Formalism for Spherically Symmetric Spacetimes in 2+1 Gravity
The ADM apporach to (2+1)-dimensional starts with a slicing of the spacetime manifold into constant-time surfaces and an induced metric . Consider a spherically symmetric two-dimensional Riemannian space . The line element on is characterized by two functions and ,
[TABLE]
We take and to be positive and non-vanishing functions. Note that is the curvature radius and is the radial line element. The Dirac-ADM action for vacuum metric is given by
[TABLE]
which represent the standard Einstein-Hilbert action for the gravitational field and Lagrangian is
[TABLE]
where is the scalar curvature of a space metric and is the extrinsic curvature for constant time surface. The scalar curvature for line element 1 is
[TABLE]
Note that dot and prime denote differentiation with respect to t and r, respectively. The extrinsic curvature is defined by
[TABLE]
where is the full covariant derivative and shift vector . For spherically symmetric, the extrinsic curvature for line element 1 is diagonal:
[TABLE]
The ADM action in terms of Eqs. 3, 4, 7 and 8 is
[TABLE]
The canonical formalism of the action is derived by differentiating the ADM action 2.1 with respect to velocities and , we obtain
[TABLE]
where the momentum is a density, while the momentum is a scalar. The velocities and can be written in terms of and as
[TABLE]
The extrinsic curvature as a function of the canonical momenta:
[TABLE]
The ADM action 2.1 can be written in the canonical form by Legendre transformation as
[TABLE]
where the super-Hamiltonian is
[TABLE]
and supermomentum
[TABLE]
2.2 Kuchař Transformation for 2+1 Gravity
In [5], Kuchař proposed a new method to express Hamiltonian and momentum constraints into simple set of constraints in (3+1)-dimensional gravity. This method is based on canonical transformation of the old variables to new canonical set. In this section, we discuss Kuchař method in (2+1)-dimensional gravity. The general metric of a spherically symmetric spacetime [6, 9] is
[TABLE]
where and are continuous fuctions of and only. Eq. 18 represent the Arnowitt-Deser-Misner (ADM) form of the 2+1 decomposition of a space-time. The general spherically symmetric (2+1)-dimensional Schwarzchild metric, in the curvature coordinates , is
[TABLE]
where conical singularity with angle . Consider that the hypersurface be a leaf of a foliation
[TABLE]
We substitute Eq. 20 into Eq. 19 and get
[TABLE]
By comparing Eq. 2.2 with Eq. 18, we get
[TABLE]
The lapse function and shift function are given by
[TABLE]
,respectively. To calculate , we substitute Eqs. 25 and 26 into Eq. 10
[TABLE]
The Schwarzchild mass can be calculated using Eqs. 22, which is equal to
[TABLE]
We found that functions and are canonically conjugate variables and the dynamical variable can be denoted by . The new momentum is written interms of the old momentum and a dynamical variable ,
[TABLE]
where does not depend on . Finally, the transformation form is
[TABLE]
The transformation form is not valid at the Horizon. We used natural units in which . The Liouville form is
[TABLE]
Eqs, 16 and 17 can be represented by simple set of constraints:
[TABLE]
except on the horizon.
2.3 Canonical analysis for (2+1)-dimensional gravity with thin shell
In this section, we briefly describe the canonial formalism for a spherically symmetric gravitational field of the 2+1 decomposition of a space-time with thin shell. The action for this system is
[TABLE]
The first and third term in Eq. 2.3 represent the standard Einstein-Hilbert action for the gravitational field and dust thin shell action, respectively. The standard Einstein-Hilbert action is given by Eq. 2.3 and shell part of the action is
[TABLE]
where hats denote the value of variables on the shell and is the rest mass of the shell. The explicit form of the action 2.3 in the Hamiltonian form becomes
[TABLE]
where is the total mass of the combined gravity-shell system and is the momentum conjugate to which is equal to
[TABLE]
and the super-Hamiltonian of the shell is
[TABLE]
and supermomentum of the shell is
[TABLE]
The regular contribution to the constraints is the same as in vacuum and is valid in inner and outer regions of the shell. There is also a delta-functional contribution on the shell which has to be combined with the shell hamiltonian. As a result we get the shell constraints:
[TABLE]
and
[TABLE]
where the square brackets mean the jump of the field across the shell.
Now we go to the Kuchar variables, solve the constraints in the inner and outer regions, and plug the solution back into the action. It turns out that the the kinetic term in the bulk disappears due to the bulk constraints, kinetic term containing on the shell disappears due to the constraint (44), and all that remains is the contribution from boundary terms that appear after the Kuchar canonical transformation.
[TABLE]
Here is the Killing time evaluated at the shell,
[TABLE]
and is the constraint (44).
Then we express the constraint (44) in terms of the shell canonical variables and :
[TABLE]
From the resulting action we find equations of motion for
[TABLE]
which leads to another constraint
[TABLE]
Unlike 3+1 dimensional situation the constraints (47) and (49) are now first class. Substituting (48) into (47) we recover the Israel equation for the shell:
[TABLE]
Finally, squaring the two constraints (47) and (49) and adding them we find the single Hamiltonian constraint describing the dynamics of the shell:
[TABLE]
This is the constraint to be used in quantum theory.
Here we have to point out that many of the above equations contain square roots in them which are not single-valued functions. Different choices of the signs in front of square roots correspond to different sectors of the phase space of the model, which are pictured as different regions on Penrose diagrams. The same corresponds to two different points in momentum space. Besides the Killing time and the radial momentum diverge at the gravitational collapse point (horizon). Beyond this point equation (51) does not have solutions in real variables. In other words the theory is formulated in variables which do not cover the phase space globally.
One way to circumvent the problem is to introduce a complex phase space ant to assemble different patches into a Riemann surface [9].
However there is still a possibility that there exists another set of real phase space variables which would cover it globally. This possibility is studied in the next section.
3 First order formalism and anti-DeSitter the momentum space
There is an example in literature when a simple gravity+matter model was given a global chart of its phase space. This is 2+1-dimensional gravity coupled to point particles [10]. In this section we generalize these results to a spherically symmetric shell.
3.1 Action principle and phase space reduction
In this section we start with a generic (non-spherically symmetric) action for 2+1 gravity coupled to a finite number of particles.
The basic variables are the triad and the connection , where are generators of sl(2)-algebra. The action reads:
[TABLE]
where is the curvature of .
The shell will be discretized ( Consisting of particles labelled by index )
[TABLE]
where is i-th particle worldline and – a fixed element of sl(2)-algebra. The condition of spherical symmetry will be imposed later.
Gravity action is invariant with respect to gauge transformations:
[TABLE]
where is an SL(2) group element, and is an sl(2) algebra element.
The shell action is not invariant. The i-th particle term transforms as
[TABLE]
where , is a parameter along the particle worldline and dot is the derivative with respect to it.
In the last term in the r.h.s. of 55 one can recognize the standard particle action as it has the form of , where , , and given the definition of satisfies the constraint . Thus the particles degrees of freedom are represented as gauge degrees of freedom evaluated at the location of the particles.
As before, to obtain a reduced action for this model we have to solve the constraints and plug the solution back into the initial action. We choose slicing so that particle worldlines move along the time coordinate and obtain the constraints by varying action (52) with respect to and :
[TABLE]
where is the location of the i-th particle. The first constraint (56) generates the first of the transformations (54) and the second generates the second.
By using transformations (54) one can put to zero simultaniously one component of and one component of . This automatically linearizes the constraints (56). However, such a gauge choice cannot be made globally, because the model has a non-trivial moduli space, containing for example the gauge parameter evaluated at the location of one particle with respect to another. Following [12] we divide the spacial slice into regions in each of which the above gauge choice could be made. Each such region should contain no more than one particle. Around each particle we draw a circle, so that the circle are connected to a common origin, but have no common boundaries, as it is shown on Fig.1. By making cuts along the circles the manifold is divided into discs, each containing a particle, and a polygon containing no particles, but connected to infinity.
For the discs it is convenient to write down the solution in polar coordinates with the origin at the location of the particles. We choose the gauge in which the radial components of and equal zero, solve the constraints, and put the solution back into an arbitrary gauge:
[TABLE]
[TABLE]
where .
And similar for polygon, for which the gauge parameters will be denoted and .
Now this solutions have to be put back into the kinetic term of the action which reads (for i-th disk):
[TABLE]
By using the identity (notice that does not depend on time, so commutes with time derivative)
[TABLE]
we find that in the first term of (58) there is a -functional contribution which cancels the second term, plus another term which is a total derivative. Thus the action for the disk collapses to its boundary:
[TABLE]
Similar for polygon, whose boundary,however, consists of edges , and the resulting action is a sum of contributions from every edge:
[TABLE]
The next step is to assemble all the above pieces of the action together and apply the condition of continuity of metric and connection across the boundary between discs and polygon.
First we convert the covariant derivative in (57) into ordinary derivative by a gauge transformation
[TABLE]
This condition violates the periodicity, so the boundary of the disk is no longer a circle but an interval. Then disc action (59) changes to
[TABLE]
and continuity conditions for metric and connection (57) take a simple form:
[TABLE]
where and are functions only of time. In the subsequent we will put , which is possible if the positions of different particles could be related by pure rotation with no translations. This is an implementation of spherical symmetry. Substituting this into (60) and (61), and combining them one obtains
[TABLE]
The integrands are total derivatives, so the result contains contributions only from the vertices of the polygon or endpoints of discs boundaries.
[TABLE]
where is the value of at the i-th vertex of the polygon.
Introduce new variables
[TABLE]
Then taking into account that
[TABLE]
and
[TABLE]
we can rewrite the second equation in (64) as
[TABLE]
To relate ’s for different ’s we use the second of the overlap conditions (62) which in particular implies that
[TABLE]
from which using the definition (65) one can derive
[TABLE]
and
[TABLE]
where the factors in the product are ordered from right to left.
Now introduce the holonomy around the full shell, which is the product of holonomies around every particle
[TABLE]
Using an obvious identity
[TABLE]
and taking into account the relation (68) we can rewrite the kinetic term of the action (66) in a simple form
[TABLE]
The action for many particle is thus collapsed to a term depending on a single variable. This became possible only because we put in (62), making use of spherical symmetry.
3.2 Constraints
To get the complete action for the shell one has to find the constraints satisfied by the variables entering (70), and add them to the kinetic term. In the definition of
[TABLE]
choose a spherically symmetric anzatz for
[TABLE]
where (not to confuse with imaginary unit) plays the role of an angular variable, and is an independent of boost parameter. For spherical symmetry means that , where is the overall bare mass of the shell. Substituting this into (71) one obtains
[TABLE]
The product of two exponents in the last equation could be calculated by Campbell–Hausdorff formula . But the terms with commutators will be of the order of and higher, and, therefore, negligible as . Finally we obtain
[TABLE]
The conjugacy class of the above holonomy is fixed by the boundary conditions at infinity. It can be evaluated as a Wilson loop of the angular component of the connection. The later can be expressed in terms of the ADM variables as
[TABLE]
and
[TABLE]
where is the ADM mass.
On the other hand from (74) it follows that
[TABLE]
Equating two last expressions for and solving for one finds
[TABLE]
This is the Israel equation (50) from the previous section with .
Now we can write down the relations between canonical momenta used here and canonical momenta from the previous section:
[TABLE]
and
[TABLE]
One can see that is always real even when and are complex. is also always real and elliptic when and hyperbolic when . In elliptic case, , when then , and when then . And similar for the hyperbolic case: , when then , and when then . In other words, provides a real global parametrization of the momentum space of the model. This momentum space is shown on Fig.2, the four described above regions labelled as I,III,II, and VI respectively. These regions can also be found on the Penrose diagram Fig.3, where they are labelled the same way.
3.3 Poisson brackets
By a change of variables the kinetic term of the action (70) can be put in the standard canonical form
[TABLE]
where labels temporal and radial components of a coordinate (momentum). The explicit form of this transformation is as follows:
[TABLE]
They satisfy the standard commutation relations
[TABLE]
It is, however, not these variables that will subsequently be used in quantization. The reason is that because is an SL(2) group element the momenta satisfy the relation
[TABLE]
This means that the momentum space of the model is a one sheet hyperboloid or 2-dimensional anti-DeSitter space, and being coordinates of 3-dimensional space in which the ADS is embedded. Momenta do not form a global chart on . The global chart on if formed by the Euler angles which parameterize the group element as follows:
[TABLE]
These are related to ’s as:
[TABLE]
It is clear that the Euler angles together with coordinates canonically conjugate to them cannot form the standard canonical set of variables. This is because the Euler angles parameterize a curved space which cannot be mapped onto a flat space, and translations in such a space do not commute.
In fact the coordinates conjugate to the Euler angles are the original translation parameters from (70). To show this we can calculate the Poisson brackets for them using the inverse of (82):
[TABLE]
This implies
[TABLE]
and expected non-commutativity of the coordinates.
4 Quantization
In this section we shall briefly describe quantization of the model using momentum (Euler angle) representation. The exposition will basically follow [11], but take into account that there is one less degree of freedom due to spherical symmetry, and that the constraint describing dynamics is now different.
4.1 Quantum kinematics
We define the kinematical states of the model as functions of from (85)
[TABLE]
single-valued functions on the entire momentum space. From the requirement of being single-valued it immediately follows the periodicity property,
[TABLE]
which will have an important consequence on the spectra of coordinates.
Next, for defining the scalar product, we need a Lorenz-invariant measure on our momentum space. It can be inferred from the Haar measure on SL(2):
[TABLE]
and the scalar product is thus:
[TABLE]
Easiest of all is to calculate the spectrum of time coordinate, , which is canonically conjugate to , and the corresponding operator is
[TABLE]
its eigenstates are
[TABLE]
where is an integer. Thus, time operator has a discrete spectrum:
[TABLE]
Notice that it is quantized in units of the Planck length, as the above equation also contains the Newton constant, which is put to 1 in this paper.
A more interesting observable is the areal radius, , which is a lorenz-invariant quantity defining the size of the shell. Because in (88) is defined as the left-invariant derivative on the group, its square is the Beltrami–Laplace operator on our momentum space:
[TABLE]
where
[TABLE]
This operator was shown in [11] to have two series of eigenvalues. One is continuous, but separated from zero, corresponds to positive,i.e. spacelike,
[TABLE]
where is a real number. The other is discrete, but containing zero, corresponds to negative, i.e. timelike,
[TABLE]
where is a non-negative integer, subject to the condition .
Physically the areal radius is spacelike outside an event horizon and timelike outside. In 2+1 dimensional gravity with no cosmological constant considered here there is no event horizons. So, one could say that only the continuous series is relevant to the present model. Thinking about the entire universe being inside an event horizon is not compatible with any sensible boundary conditions at infinity.
On the other hand one could consider a more complicated model with several concentric shells which allows for spacetimes with no boundary. In such a situation timelike areal radius is possible and the discrete series will be relevant.
However, in both cases the result for the spectrum of provide a regularization near singularity.
4.2 Physical states
Now we shall try to apply the hamiltonian constraint to obtain physical states. Unlike in the case of a particle, the group element repersenting the momentum (74 ) has a complicated dependence on the external parameter of the model, the bare mass . Because of these complication we will use some implicit expressions in this section.
The analog of the Hamiltonian constraint for a particle is (77), which in terms of the Euler angles reads
[TABLE]
It contains an extra parameter , which the wavefunction of the kinematical Hilbert space doesn’t depend of. To fix it one has to use another equation from the set (77), e.g.
[TABLE]
or
[TABLE]
with
[TABLE]
This constraints canonically commute with each other, i.e. they are first class. So, we can start with kinematical Hilbert space as space of functions of three parameters and then apply two of the above three constraints.
The solution can be written as
[TABLE]
where is an arbitrary function. As usual it is not normalizible w.r.t. the kinematical Hilbert space.
The scalar product on the physical Hilbert space can be defined in terms of functions entering (104) as
[TABLE]
It is easy to show that in the limit and it is reduced to the standard scalar product for the states of a relativistic particles in 1+1 dimensions,
[TABLE]
as expected.
From the above it follows, in particular, that unlike in 3+1 dimensional situation the spectrum of the energy-momentum of the model is fully continuous. This is not surprising: in 2+1 dimensional gravity there is no Newtonian potential, so there is no potential well to hold the shell in a bounded region.
5 Conclusion
So far we do not have the full quantum theory of the model studied in this paper. What remains to do is to describe the dynamics, in particular calculate the transition amplitudes between different eigenvalues of the areal radius of the shell.
What was shown is that the spectrum of areal radius in the case of timelike movement of the shell does not reach zero. This means that classically existing naked singularity in quantum theory is avoided. As to singularity behind a horizon, which is classically attained by spacelike movement of the shell, it belongs to a discrete spectrum of the radius and, therefore, is also regularized.
The most interesting question is whether all this can be generalized for 3+1 dimensional gravity. Some results exist on quantum kinematics of a Schwarzschild black hole in a frame of a test particle [13]. It also has such features as coordinate non-commutativity and discreteness. However it cannot be generalized to a shell dynamics the way it was done in section 3. The reason is that many body problem in 3+1 gravity is not solvable.
On the other hand almost all the work on spherical shell mentioned in the introduction was done in 3+1 spacetime dimensions. The form of the Hamiltonian constraint (although not in global coordinates) was found. The map between phase space coordinates used there and the global ones should be analogous to that found in the end of section 3.3. The only difference is the presence of Newtonian potential. On this basis the form of the Hamiltonian constraint in global coordinates could be guessed.
Acknowledgements. We would like to thank the organizers of VI International Conference ”Models in quantum field theory” dedicated to A.N.Vassiliv.
A.S. would like to thank V.A.Berezin for extensive conversations on his work, and D.Lyozin for collaboration on a closely related topic.
This work was partially supported by RFBR projects 16-02-00348 and 18-02-00264 (A.A.)
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] M. Bronstein, Phys. Zeitschr. Sow. 9 (1936)
- 2[2] Abhay Ashtekar, Parampreet Singh, Loop Quantum Cosmology: A Status Report, Class.Quant.Grav. 28 (2011) e-Print: ar Xiv:1108.0893
- 3[3] Alexander A. Andrianov, Oleg O. Novikov, Chen Lan, Quantum cosmology of multifield scalar matter: Some exact solutions Theor.Math.Phys. 184 (2015) 1224-1233, Teor.Mat.Fiz. 184 (2015) no.3, 380-391 e-Print: ar Xiv:1503.05527
- 4[4] Israel, W. Singular hypersurfaces and thin shells in general relativity, Il Nuovo Cimento B 44,1967
- 5[5] Kuchař, Karel V., Geometrodynamics of Schwarzschild black holes, Phys. Rev. D 50, 1994,
- 6[6] Hájíček, P., Spherically symmetric gravitating shell as a reparametrization invariant system, Phys. Rev. D 57, 1998
- 7[7] Louko, Jorma and Whiting, Bernard F. and Friedman, John L., Hamiltonian space-time dynamics with a spherical null dust shell, Phys. Rev.D 57, 1998
- 8[8] Hájíček, Petr and Kiefer, Claus, Embedding variables in the canonical theory of gravitating shells, Nucl. Phys.B 603, 2001
