The approaches I will describe below have yielded a variety of interesting results on the discrete simulation methods themselves, as well as on the geometric properties of typical configurations, the role of gauge invariance, the phase structure, and the inclusion of matter, to name just a few. So far, they have not been successful in providing convincing evidence for the existence of a nontrivial fourdimensional quantum theory of gravity (neither, of course, have other methods). There are considerable technical difficulties in performing analytical computations and setting up sufficiently big numerical simulations. It is a theoretical challenge to come up with more realistic models that may lead to a truly interacting gravitational theory.
I have divided the discrete approaches into three categories, which I will discuss in the order of their chronological appearance: the gaugetheoretic formulations using connection variables, based on firstorder descriptions of Einstein gravity and two metric formulations; the quantum Regge calculus program; and the more recent method of dynamical triangulations, which use simplicial instead of hypercubic lattice discretizations. They mainly deal with pure gravity, with the possible inclusion of a cosmological constant term and higherorder derivative terms, and in some instances, mattercoupling.
This review covers about 200 papers, the earliest of which appeared in 1979 (not counting the 1964 article by Leutwyler [
145]
, who used a lattice approximation to perform a gravitational sum over histories à la Feynman). Much of the motivation for the discrete investigations comes of course from other sources, which for the most part I do not cite explicitly. The reader is referred to the individual articles for a more complete bibliography of the background material.
2 Gaugetheoretic discretizations of gravity
2.1 Lagrangian treatment. Introduction
This area of research was inspired by the success of nonperturbative lattice methods in treating nonabelian gauge theories [
173]
. To apply some of their techniques, gravity has to be brought into a gaugetheoretic, firstorder form, with the puregravity Lagrangian
$$\begin{array}{c}S[A,e]={\int}_{M}e\wedge e\wedge R\left[A\right],\end{array}$$ 
(1)

where the
$so(3,1)$
valued spinconnection
${A}_{\mu}^{IJ}$
(with curvature
$R$
) and the vierbein
${e}_{\mu}^{I}$
are considered as independent variables. The important feature of ( 1 ) is that
${A}_{\mu}^{IJ}$
is a gauge potential, and that the action – in addition to its diffeomorphisminvariance – is invariant under local frame rotations.
Variation with respect to
${e}_{\mu}^{I}$
leads to the metricity condition,
$$\begin{array}{c}{\mathcal{D}}_{[\mu}{e}_{\nu ]}^{I}\equiv {\partial}_{[\mu}{e}_{\nu ]}^{I}+{A}_{[\mu J}^{I}{e}_{\nu ]}^{J}=0,\end{array}$$ 
(2)

which can be solved to yield the unique torsionfree spin connection
$A=A\left[e\right]$
compatible with the
${e}_{\mu I}$
. There are some obvious differences with usual gauge theories: i) the action ( 1 ) is linear instead of quadratic in the curvature twoform
${R}_{\mu \nu}^{IJ}$
of
$A$
, and ii) it contains additional fields
${e}_{\mu}^{I}$
.
Substituting the solution to (
2 ) into the action, one obtains
$S\left[e\right]={\int}_{M}{d}^{4}x(dete)R$
, where
$R$
denotes the fourdimensional curvature scalar. This expression coincides with the usual Einstein action
$\int {d}^{4}x\sqrt{detg}R$
only for
$dete>0$
.
Most of the lattice gauge formulations I will discuss below share some common features. The lattice geometry is hypercubic, defining a natural global coordinate system for labelling the lattice sites and edges. The gauge group
$G$
is
$SO(3,1)$
or its “Euclideanized” form
$SO\left(4\right)$
, or a larger group containing it as a subgroup or via a contraction limit. Local curvature terms are represented by (the traces of )
$G$
valued Wilson holonomies
${U}_{\sqcap \bigsqcup}$
around lattice plaquettes. The vierbeins are either considered as additional fields or identified with part of the connection variables. The symmetry group of the lattice Lagrangian is a subgroup of the gauge group
$G$
, and does not contain any translation generators that appear when
$G$
is the Poincaré group.
When discretizing conformal gravity (where
$G=SO(5,1)$
) or higherderivative gravity in firstorder form, the metricity condition on the connection has to be imposed by hand. This leads to technical complications in the evaluation of the functional integral.
The diffeomorphism invariance of the continuum theory is broken on the lattice; only the local gauge invariances can be preserved exactly. The reparametrization invariance reemerges only at the linearized level, i.e. when considering small perturbations about flat space.
2.2 Smolin's lattice model
The first gaugetheoretic model for lattice gravity is due to Smolin [
186]
, based on the continuum formulation of MacDowell and Mansouri [
157]
, with de Sitter gauge group
$O(3,2)$
or
$O(4,1)$
, and a Lagrangian of
${R}^{2}$
type,
$$\begin{array}{c}S=\int {d}^{4}x{\epsilon}^{\mu \nu \rho \sigma}{\epsilon}_{IJKL}{\stackrel{~}{R}}_{\mu \nu}^{IJ}{\stackrel{~}{R}}_{\rho \sigma}^{KL},\end{array}$$ 
(3)

where the components of
${\stackrel{~}{R}}_{\mu \nu}$
are related to those of the usual curvature tensor by
$$\begin{array}{c}{\stackrel{~}{R}}_{\mu \nu}^{IJ}={R}_{\mu \nu}^{IJ}\pm \frac{1}{{l}^{2}}({e}_{\mu}^{I}{e}_{\nu}^{J}{e}_{\mu}^{J}{e}_{\nu}^{I}).\end{array}$$ 
(4)

Although the underlying gauge potentials
${A}_{\mu}$
are
$o(3,2)$
or
$o(4,1)$
valued, the action is only invariant under the 6dimensional subgroup of Lorentz transformations. The theory contains a dimensionful parameter
$l$
. The gauge potentials associated with the internal 5direction are identified with the frame fields
${e}_{\mu}^{I}$
, and the action can be decomposed into the usual Einsteinterm ( 1 ) plus a cosmological constant term with
$\lambda \sim \frac{1}{{l}^{2}}$
and a topological
$R\wedge R$
term. Smolin analyzed its lattice discretization, and found both a weakand a strongcoupling phase, with respect to the dimensionless coupling constant
$g\sim \frac{\sqrt{G}}{l}$
. He performed a weakcoupling expansion about flat space and rederived the usual propagator. In the strongcoupling regime he found massive excitations and a confining property for spinors.
2.3 Numerical implementation
Caracciolo and Pelissetto [
63,
64,
65,
66,
67]
performed a numerical investigation of the phase structure of Smolin's model. Using the compact group
$SO\left(5\right)$
and its associated Haar measure, their findings confirmed the twophase structure: a strongcoupling phase with a confining property and presence of exponential clustering, and a weakcoupling phase dominated by a class of topological configurations, with vanishing vierbein. However, their Monte Carlo data (on
${4}^{4}$
and
${8}^{4}$
lattices with periodic boundary conditions) indicated strongly that the transition was firstorder, even if the measure was generalized by a factor of
$dete{}^{p}$
,
$p\in [0,150]$
[
67]
.
2.4 Other gauge formulations
Das et al [
84]
latticediscretized an
$Sp\left(4\right)$
invariant Lagrangian due to West [
194]
, leading to a functional form
$S\left[U\right]={\sum}_{n}\sqrt{Tr\left(\right({U}_{n}\frac{1}{4}Tr{U}_{n}{)}^{2})}$
, where
${U}_{n}$
denotes a sum of doubleplaquette holonomies based at the vertex
$n$
. Details of the functional integration were not spelled out. The squareroot form of the Lagrangian does not make it very amenable to numerical investigations (see also the comments in [
139]
).
Mannion and Taylor [
159]
suggested a discretization of the
$(\int e\wedge e\wedge R)$
Lagrangian without cosmological term, with the vierbeins
$e$
treated as extra variables, and the compactified Lorentz group
$SO\left(4\right)$
.
Kaku [
138]
has proposed a lattice version of conformal gravity (see also [
137,
139]
), based on the group
$O(4,2)$
, the main hope being that unitarity could be demonstrated nonperturbatively. As usual, a metricity constraint on the connection has to be imposed by hand.
A lattice formulation of higherderivative gravity, containing fourthorder
${R}^{2}$
terms was given by Tomboulis [
189]
. The continuum theory is renormalizable and asymptotically free, but has problems with unitarity. The motivation for this work was again the hope of realizing unitarity in a lattice setting. The squareroot form of the Lagrangian is similar to that of Das et al. It is
$O\left(4\right)$
invariant and supposedly satisfies reflection positivity. Again, the form of the Lagrangian and the measure (containing a
$\delta $
function of the notorsion constraint) is rather complicated and has not been used for a further nonperturbative analysis.
Kondo [
142]
employed the same framework as Mannion and Taylor, but introduced an explicit symmetrization of the Lagrangian. He claimed that the cluster expansion goes through just as in lattice YangMills theory, leading to a positive mass gap.
2.5 Proving reflection positivity
Various properties of gauge gravity models were analyzed by Menotti and Pelissetto in a series of papers in the mideighties. They first studied a discrete
$O\left(4\right)$
version of the Lagrangian ( 1 ) [
165]
, including the
$O\left(4\right)$
Haar measure and a general local
$O\left(4\right)$
invariant measure for the vierbein fields, and showed that reflectionpositivity holds only for a restricted class of functions.
Furthermore, expanding about flat space, and after appropriate gaugefixing, they discovered a doubling phenomenon similar to that found for chiral fermions in lattice gauge theory, a behaviour that also persists for different gaugefixings. One finds the same mode doubling also for a flatbackground expansion of conformal gravity [
167]
. In the same paper, they gave a unified treatment of Poincaré, de Sitter and conformal gravity, and showed that reflection positivity for
$O\left(4\right)$
gravity (as well as for the two other gauge groups) holds exactly and for general functions only provided a signature factor sign
$(dete)$
is included in the Lagrangian. To ensure the convergence of the functional integration, one has to introduce a damping factor for the vierbeins in the measure, both for Poincaré and conformal gravity [
167]
. An extension of the results of [
167]
to supergravity with the superPoincaré group is also possible [
62]
. One can prove reflection positivity and finds a matching gravitino doubling in the perturbative expansion.
2.6 The measure
One can derive nontrivial conditions on the lattice measure by imposing the SlavnovTaylor identities in a perturbative lattice calculation. Starting from the general form of the measure for Poincaré gravity,
$$\begin{array}{c}{\prod}_{n\mu \nu \rho \sigma}det\left({e}_{\mu}\right(n),{e}_{\nu}(n),{e}_{\rho}(n),{e}_{\sigma}(n\left){)}^{N/16}{\prod}_{n\mu}f\right({e}_{\mu}^{2}\left(n\right)\left){\prod}_{n,\mu >0}d{e}_{\mu}\right(n\left){\prod}_{n,\mu >0}d{U}_{L}\right(n,n+\mu ),\end{array}$$ 
(5)

Menotti and Pelissetto [
166]
performed a oneloop calculation and found that the SlavnovTaylor identity can be satisfied for a particular choice of
$N$
,
${f}_{1},{f}_{2}$
(momentum expansion coefficients of
$f$
), and of the cosmological constant
$\lambda $
. The solution still depends on a real parameter
$\xi $
(related to a residual noninvariance under rotations). This makes it difficult to draw any immediate conclusions on the structure of the full, nonperturbative measure.
2.7 Assorted topics
Caselle et al [
70]
proposed a lattice action that is genuinely Poincaréinvariant, at the price of introducing additional lattice “coordinate variables”. They also suggested a compact
$O\left(5\right)$
formulation which reduces to the Poincaré form in the limit as the length of some preferred
$O\left(5\right)$
vector is taken to infinity, as well as a superversion involving the graded Poincaré group. The same authors in [
69]
put forward an argument for why doubling should appear in general gravity plus matter systems.
Reisenberger [
175]
has recently suggested a gaugetheoretic path integral based on the Plebanski action for Euclidean gravity. He discretizes the theory on a simplicial or hypercubic lattice with groupand algebravalued fields. A metricity constraint needs to be imposed on the basic spin1 fields, which it turns out is difficult to treat exactly.
2.8 Summary
The gaugetheoretic Lagrangian lattice approaches are afflicted by a number of technical difficulties.
Reflection positivity can be shown for some of the models, but generally requires the inclusion of a factor sign
$(dete)$
in the Lagrangian. Obtaining qualitative nonperturbative information about the phase structure requires a nontrivial measure input. The complicated functional form of the Lagrangian and the metricity condition that has to be imposed via a Lagrange multiplier and the corresponding functional integration do not make conformal and higherderivative theories attractive candidates for numerical simulations. The compact version of Smolin's de Sitter gravity is still the simplest model, but its numerical investigation did not yield interesting results.
2.9 Hamiltonian treatment. Introduction
Relatively little work has been done on discretized Hamiltonian formulations of gravity. This can in part be understood from the fact that the numerical methods available for lattice gauge field theories rely mostly on the Euclidean pathintegral description. Unfortunately the relation between the Lagrangian and Hamiltonian quantizations for generally covariant theories without a fixed background is far from clear.
In the usual metric formulation, the complicated nonpolynomial form of the Hamiltonian constraint has been a longstanding problem. In this framework, neither the functional form of the quantum representations nor the nature of the quantization problems suggest that a discrete approach might yield any advantages. This situation has improved with the introduction of new Hamiltonian gaugetheoretic variables by Ashtekar [
22,
23]
.
By a Hamiltonian lattice approach one usually means a formulation in which the time variable is left continuous, and only the spatial 3slices are discretized. In continuum gravity, the 3+1 decomposition leads to the (nonLie) Dirac algebra of the threedimensional diffeomorphism generators and the Hamiltonian constraint, associated with the deformation of threesurfaces imbedded in fourspace. One usually requires this algebra to be realized in the quantum theory, without anomalous terms, for a set of selfadjoint quantum constraint operators, for some factorordering.
Since a discretization of spacetime breaks the diffeomorphism invariance, there is no reason to expect the Dirac algebra to be preserved in any discrete approach, even classically. This raises the question of whether and in what form part of the diffeomorphism symmetry can still be realized at the discrete level. Using gaugetheoretic variables, one can maintain the exact local gauge invariance with respect to the internal degrees of freedom, but there is no analogous procedure for treating the coordinate invariance.
2.10 Hamiltonian lattice gravity
There is a gaugetheoretic Hamiltonian version of gravity defined on a cubic lattice, which in many aspects resembles the Lagrangian gauge formulations described earlier. It also is virtually the only discrete Hamiltonian formulation in which some progress has been achieved in the quantization (see also [
155]
for a recent review).
Renteln and Smolin [
180]
were the first to set up a continuoustime lattice discretization along the lines of Hamiltonian lattice gauge theory. Their basic configuration variables are the link holonomies
$U\left(l\right)$
of the spatial Ashtekar connection along the edges. The lattice analogues of the canonically conjugate pairs
$\left({A}_{a}^{i}\right(x),{E}_{i}^{a}(x\left)\right)$
are the link variables
$\left(U\right(l{{)}_{A}}^{B},{p}_{i}\left(l\right))$
, with Poisson brackets
$$\begin{array}{ccc}& & \left\{{{U}_{A}}^{B}\right(n,\hat{a}),{{U}_{C}}^{D}(m,\hat{b}\left)\right\}=0,\end{array}$$ 
(6)

$$\begin{array}{ccc}& & \left\{{p}_{i}\right(n,\hat{a}),{{U}_{A}}^{C}(m,\hat{b}\left)\right\}=\frac{1}{2}{\delta}_{nm}{\delta}_{\hat{a}\hat{b}}{{\tau}_{iA}}^{B}{{U}_{B}}^{C}(n,\hat{a}),\end{array}$$ 
(7)

$$\begin{array}{ccc}& & \left\{{p}_{i}\right(n,\hat{a}),{p}_{j}(m,\hat{b}\left)\right\}={\delta}_{nm}{\delta}_{\hat{a}\hat{b}}{\epsilon}_{ijk}{p}_{k}(n,\hat{a}),\end{array}$$ 
(8)

with the
$SU\left(2\right)$
generators satisfying
$[{\tau}_{i},{\tau}_{j}]=2{\epsilon}_{ijk}{\tau}_{k}$
. Lattice links
$l=(n,\hat{a})$
are labelled by a vertex
$n$
and a lattice direction
$\hat{a}$
. These relations go over to the usual continuum brackets in the limit as the lattice spacing
$a$
is taken to zero. In this scheme, they wrote down discrete analogues of the seven polynomial firstclass constraints, and also attempted to interpret the action of the discretized diffeomorphism and Hamiltonian constraints in terms of their geometric action on lattice Wilson loop states.
2.11 The measure
There is a natural measure for the quantum theory, given by the product over all lattice edges of the Haar measures
$dg$
. However, since the Ashtekar connections
$A$
are complexvalued, the gauge group is the noncompact group
$SO(3,C)=SL(2,C)$
, and the gaugeinvariant Wilson loop functions are not squareintegrable. For the alternative formulation in terms of real SU(2)variables (see below), these problems are not present. An alternative heat kernel measure
$d\nu $
for holomorphic
$SL(2,C)$
holonomies on the lattice was used in [
146,
99]
.
2.12 The constraint algebra
This line of research was continued by Renteln [
179]
, who proved that for a particular factorordering (all momenta to the left), the subalgebra of the discretized diffeomorphism constraints (smeared by lapse functions
$N$
),
$$\begin{array}{c}{\sum}_{n}{\mathcal{V}}^{latt}[N,n):={\sum}_{n}{N}^{latt}(n,\hat{a})Tr\left({U}_{\sqcap {\bigsqcup}_{ab}}{\tau}_{i}\right){p}_{i}(n,\hat{b}),\end{array}$$ 
(9)

closes in the limit as the lattice spacing is taken to zero. This calculation was later extended to a variety of different symmetrizations for the lattice operator and to an arbitrary factorordering of the form
$\alpha Tr\left(\hat{U}\tau \right)\hat{p}+(1\alpha )\hat{p}Tr\left(\hat{U}\tau \right)$
, with
$0\le \alpha \le 1$
[
156]
. Again, one does not find any quantum anomalies. It would be highly desirable to extend this result to commutators involving also the discretized Hamiltonian constraint and to find the explicit functional form of the anomalies, if there were any.
2.13 Solutions to the WheelerDeWitt equation
Part of solving the canonical quantum theory is to determine the states annihilated by the Hamiltonian constraint
$\hat{H}$
. It was shown by Loll [
146]
that solutions exist in the RentelnSmolin formulation, where
$\hat{H}\left(n\right)={\sum}_{a,b}{\epsilon}^{ijk}Tr\left(\hat{U}\right(n,\sqcap {\bigsqcup}_{\hat{a}\hat{b}}\left){\tau}_{k}\right){\hat{p}}_{i}(n,a){\hat{p}}_{j}(n,b)$
. They are given by multiple, nonintersecting Polyakov loops (the lattice is assumed to have compact topology
${T}^{3}$
). Such solutions are trivial in the sense that they correspond to quantum states “without volume”.
The difficulties one encounters when trying to find other solutions is illustrated by the explicit calculations for the
$1\times 1\times 1$
lattice in [
146]
.
The search for solutions was continued by Ezawa [
99]
(see also [
100]
for an extensive review), who used a symmetrized form of the Hamiltonian. His solutions depend on multiple plaquette loops
$({U}_{\sqcap \bigsqcup}{)}^{k}$
, where a single lattice plaquette
$\sqcap \bigsqcup $
is traversed by the loop
$k$
times. The solutions are less trivial than those formed from Polyakov loops, since they involve kinks, but they are still annihilated by the volume operator.
A somewhat different strategy was followed by Fort et al [
102]
, who constructed a Hamiltonian lattice regularization for the calculation of certain knot invariants. They defined lattice constraint operators in terms of their geometric action on lattice Wilson loop states, and reproduced some of the formal continuum solutions to the polynomial Hamiltonian constraint of complex Ashtekar gravity on simple loop geometries.
2.14 The role of diffeomorphisms
As in other discrete formulations, the spatial diffeomorphism group cannot be realized exactly on the lattice, and the only obvious symmetries of a cubic lattice are discrete rotations and overall translations. The commutator computation described above indicates that one may be able recover the diffeomorphism invariance in a suitable continuum limit. Corichi and Zapata [
83]
have suggested the presence of a residual diffeomorphism symmetry in the lattice theory, under which, for example, all nonintersecting Wilson loop lattice states would be identified.
One can try to interpret the lattice theory as a manifestly diffeomorphisminvariant construction, with the lattice representing an entire diffeomorphism equivalence class of lattices embedded in the continuum [
146]
. In order to make this interpretation consistent, one should modify the functional form of either the Hamiltonian or the measure, in such a way that the commutator of two lattice Hamiltonians vanishes, as
$a\to 0$
.
2.15 The volume operator
An important quantity in Hamiltonian lattice quantum gravity is the volume operator, the quantum analogue of the classical volume function
$\int {d}^{3}x\sqrt{detg}$
. The continuum dreibein determinant
$detE\left(x\right)$
(with
$detE(x\left)\right=detg$
), has a natural lattice analogue, given by
$$\begin{array}{c}detp\left(n\right):=\frac{1}{3!}{\epsilon}_{abc}{\epsilon}^{ijk}{p}_{i}(n,\hat{a}){p}_{j}(n,\hat{b}){p}_{k}(n,\hat{c}).\end{array}$$ 
(10)

Apart from characterizing geometric properties of lattice quantum states, it is needed in the construction of the quantum Hamiltonian of the real connection approach.
Loll [
147]
showed that a lattice Wilson loop state has to have intersections of valence at least 4 in order not to be annihilated by the volume operator
${\sum}_{n}\sqrt{\widehat{detp\left(n\right)}}$
. This result is independent of the choice of gauge group (
$SU\left(2\right)$
or
$SL(2,C)$
). The volume operator has discrete eigenvalues, and part of its nonvanishing spectrum, for the simplest case of fourvalent intersections, was first calculated by Loll [
149]
. These spectral calculations were confirmed by De Pietri and Rovelli [
94]
, who derived a formula for matrix elements of the volume operator on intersections of general valence (as did Thiemann [
188]
).
This still leaves questions about the spectrum itself unanswered, since the eigenspaces of
$\widehat{detp\left(n\right)}$
grow rapidly, and diagonalization of the matrix representations becomes a technical problem. Nevertheless, one can achieve a better understanding of some general spectral properties of the lattice volume operator, using symmetry properties. It was observed in [
149]
that all nonvanishing eigenvalues of
$\widehat{detp\left(n\right)}$
come in pairs of opposite sign. Loll subsequently proved that this is always the case [
151]
. A related observation concerns the need for imposing an operator condition
$\widehat{detp}>0$
on physical states in nonperturbative quantum gravity [
151]
, a condition which distinguishes its state space from that of a gauge theory already at a kinematical level.
The symmetry group of the cubic threedimensional lattice is the discrete octagonal group
$\mathcal{O}$
, leaving the classical local volume function
$detp\left(n\right)$
invariant. Consequently, one can find a set of operators that commute among themselves and with the action of the volume operator, and simplify its spectral analysis by decomposing the Hilbert space into the corresponding irreducible representations [
154]
. This method is most powerful when applied to states which are themselves maximally symmetric under the action of
$\mathcal{O}$
, in which case it leads to a dramatic reduction of the dimension of the eigenspaces of
$\widehat{detp\left(n\right)}$
.
In addition to the volume operator, one may define geometric lattice operators measuring areas and lengths [
150,
152]
. They are based on (nonunique) discretizations of the continuum spatial integrals of the square root of the determinant of the metric induced on subspaces of dimension 2 and 1. For the case of the length operator, operatorordering problems arise in the quantization.
2.16 The real dynamics
To avoid problems with the noncompactness of the gauge group and the formulation of suitable “quantum reality conditions”, Barbero [
25]
advocated to use a real su(2)connection formulation for Lorentzian continuum gravity. This can be achieved, at the price of having to deal with a more complicated Hamiltonian constraint. Loll [
148,
153]
translated the real connection formulation to the lattice and studied some of the differences that arise in comparison with the complex approach.
Adding for generality a cosmological constant term, this leads to a lattice Hamiltonian
$$\begin{array}{c}{H}^{latt}={\sum}_{n}{N}^{latt}\left(n\right)\left[{\mathcal{\mathscr{H}}}_{kin}^{latt}+{\mathcal{\mathscr{H}}}_{pot}^{latt}+\sqrt{\lambda}G\sqrt{detp\left(n\right)}\right],\end{array}$$ 
(11)

where schematically
${\mathcal{\mathscr{H}}}_{kin}^{latt}=(detp{)}^{1/2}Tr(U\tau ){p}^{2}$
,
${\mathcal{\mathscr{H}}}_{pot}^{latt}=(Tr(\tau Up\tau U)p{)}^{2}{p}^{6}$
$(detp{)}^{5/2}$
.
This regularized Hamiltonian is welldefined on states with
$\widehat{detp}\ne 0$
, but its functional form is not simple. The negative powers of the determinant of the metric can be defined in terms of the spectral resolution of
$\widehat{detp}$
. The type of representation and regularization enables one to handle this nonpolynomiality.
2.17 Summary
Some progress has been achieved in Hamiltonian lattice gravity, using discrete analogues of the Ashtekar variables. The quantization program is still openended, and no physically nontrivial solutions to the WheelerDeWitt equation are known. The fact that one can study and evaluate geometric operators provides useful characterization of quantum states.
An analysis of the spectrum of the volume operator is crucial for handling the nonpolynomial terms in the quantum Hamiltonian constraint. It will be necessary to find a suitable truncation or approximation to simplify further the spectral analysis of the WheelerDeWitt operator. A suitable quantum analogue of the continuum limit
$a\to 0$
has not yet been established, and in this regard the Hamiltonian ansatz does not go beyond the results obtained in the Lagrangian formulations described earlier.
Why should one bother with a Hamiltonian quantization at all? Typical quantities one wants to study in a discrete pathintegral approach to gravity are transition amplitudes between threegeometries on different spatial slices. This is not complete without a specification of the corresponding quantum states, which are in principle elements of Hilbert spaces of discrete threegeometries of the type described above.
3 Quantum Regge Calculus
3.1 Path integral for Regge calculus
A pathintegral quantization of 4d Regge calculus was first considered in the early eighties [
181,
104,
81]
. This approach goes back to Regge [
174]
, who proposed approximating Einstein's continuum theory by a simplicial discretization of the metric spacetime manifold and the gravitational action.
Its local building blocks are foursimplices
$\sigma $
. The metric tensor associated with each simplex is expressed as a function of the squared edge lengths
${l}^{2}$
of
$\sigma $
, which are the dynamical variables of this model. For introductory material on classical Regge calculus and simplicial manifolds, see [
168,
187,
104,
124,
106,
197,
190,
8]
; various quantum aspects are reviewed in [
106,
107,
112,
110,
116,
133,
40,
160,
190,
198]
.
One may regard a Regge geometry as a special case of a
continuum Riemannian manifold, a socalled piecewise flat manifold, with a flat metric in the interior of its 4simplices
$\sigma $
, and singular curvature assignments to its twosimplices
$b$
(the bones or hinges).
The Einstein action with cosmological term in the Regge approach is given by
$$\begin{array}{c}{S}^{Regge}\left[{l}_{i}^{2}\right]={\sum}_{bonesb}{V}_{b}(\lambda k\frac{{A}_{b}{\delta}_{b}}{{V}_{b}})\sim \int {d}^{4}x\sqrt{g}(\lambda \frac{k}{2}R),\end{array}$$ 
(12)

where
$k=\frac{1}{8\pi G}$
,
${A}_{b}$
is the area of a triangular bone,
${\delta}_{b}=2\pi {\sum}_{\sigma \supset b}\theta (\sigma ,b)$
the deficit angle there, and
${V}_{b}$
a local fourvolume element.
$\theta (\sigma ,b)$
is the angle between the two 3simplices of
$\sigma $
(
$\pi $
minus the angle between their inward normals) intersecting in
$b$
. Hartle and Sorkin [
128]
have generalized ( 12 ) to the case of manifolds with boundary. The boundary contribution to the action is given by
$$\begin{array}{c}{\sum}_{b\subset boundary}{A}_{b}{\psi}_{b}\sim \int {d}^{3}x\sqrt{{g}^{\left(3\right)}}K,\end{array}$$ 
(13)

where
${\psi}_{b}$
is the angle between the normals of the two threesimplices meeting at
$b$
. The Euclidean path integral on a finite simplicial complex of fixed connectivity takes the form
$$\begin{array}{c}Z(k,\lambda )=\int \mathcal{D}l{e}^{{S}^{Regge}\left[{l}_{i}^{2}\right]}\end{array}$$ 
(14)

with
$\int \mathcal{D}l$
representing the discrete analogue of the sum over all metrics. A crucial input in ( 14 ) is the choice of an appropriate measure
$\mathcal{D}l$
. In general, a cutoff is required for both short and long edge lengths to make the functional integral convergent [
104,
81,
18]
. One is interested in the behaviour of expectation values of local observables as the simplicial complex becomes large, and the existence of critical points and longrange correlations, in a scaling limit and as the cutoffs are removed.
A first implementation of these ideas was given by Roček and Williams [
181,
183]
. They obtained a simplicial lattice geometry by subdividing each unit cell of a hypercubic lattice into simplices.
Their main result was to rederive the continuum free propagator (see also [
101]
for related results) in the limit of weak perturbations about flat space. This calculation can be repeated for Lorentzian signature [
196]
.
Some nonperturbative aspects of the path integral were investigated in [
183]
(see also [
182]
). In this work, discrete analogues of spacetime diffeomorphisms are defined as the local link length transformations which leave the action invariant, and go over to translations in the flat case. It is argued that an approximate invariance should exist in 4d. One may define analogues of local conformal transformations on a simplicial complex by multiplication with a positive scale factor at each vertex, but the global group property is incompatible with the existence of the generalized triangle inequalities.
3.2 Higherderivative terms
Simplicial analogues of higherderivative terms were introduced in [
117,
119]
. In the continuum, with an appropriate choice of coupling constants, their inclusion makes the path integral less illbehaved. The simplest higherderivative term in Regge calculus is given by
$$\begin{array}{c}{\sum}_{b}\frac{{A}_{b}^{2}{\delta}_{b}^{2}}{{V}_{b}}\sim \int {d}^{4}x\sqrt{g}{R}_{\mu \nu \rho \sigma}{R}^{\mu \nu \rho \sigma},\end{array}$$ 
(15)

with
${V}_{b}$
denoting the local Voronoi fourvolume at
$b$
. The fact that ( 15 ) should not be identified with
$\int {d}^{4}x\sqrt{g}{R}^{2}$
is less surprising in light of the classical result [
81,
82]
, that Regge's expression
${A}_{b}{\delta}_{b}$
for the scalar curvature (for
$d>2$
) converges to its continuum counterpart not pointwise, but only after integration, i.e. “in the sense of measures” (see also [
103,
101]
, where similar convergence properties were studied by using an imbedding into a sufficiently large vector space
$I{R}^{N}$
).
More complicated highercurvature terms can in principle be constructed, using a simplicial analogue of the Riemann tensor (see, for example, [
174,
106,
119,
55]
), but have up to now not been used in numerical simulations. A related proposal by Ambjørn et al [
19]
is to include terms in the action that depend on higher powers of the deficit angle
${\delta}_{b}$
, as well as terms containing powers of the solid angle
${\delta}_{v}$
at a vertex. The introduction of local vierbeins and parallel transporters is also necessary if one considers fermion coupling [
104,
176]
.
3.3 First simulations
The first numerical studies of the Regge action were undertaken by Berg [
38,
39]
and Hamber and Williams [
118,
119]
. Berg performed a MonteCarlo simulation of the purecurvature action for hypercubic
${2}^{4}$
and
${3}^{4}$
lattices with simplicial subdivision (see [
40]
for a description of the method).
He used the scaleinvariant measure
$\mathcal{D}l={\prod}_{i}\frac{d{l}_{i}}{{l}_{i}}$
. To avoid the divergence that results from rescaling the link lengths, he kept the total volume constant by performing an overall length rescaling of all links after each move. This amounts to fixing a typical length scale
${l}_{0}:=({v}_{0}{)}^{\frac{1}{4}}$
, where
${v}_{0}$
is the expectation value of the 4simplex volume.
For
$k=0$
, he found a negative average curvature
$\langle R\rangle $
, and some evidence for a canonical scaling behaviour of lengths, areas and volumes. For
$k=\pm 0.3$
, he obtained a negative (positive) average deficit angle
$\langle \delta \rangle $
and a positive (negative)
$\langle R\rangle $
. A more detailed analysis for
$0\le k\le 0.1$
led Berg [
39]
to conclude that there exists a critical value
${k}^{c}$
(presumably a firstorder transition [
40]
), below which
$\langle R\rangle $
is convergent, whereas above it diverges. (Myers [
169]
has conjectured that it may be possible to perform a similar analysis for MonteCarlo data for the Lorentzian action.) By contrast, Hamber and Williams [
118]
simulated the higherderivative action
$$\begin{array}{c}S\left[{l}_{i}^{2}\right]={\sum}_{b}[\lambda {V}_{b}k{\delta}_{b}{A}_{b}+a{A}_{b}^{2}{\delta}_{b}^{2}/{V}_{b}]\end{array}$$ 
(16)

on
${2}^{4}$
and
${4}^{4}$
lattices, using a timediscretized form of the Langevin evolution equation (see also [
106]
). For technical reasons, one uses barycentric instead of Voronoi volumes. They employed the scaleinvariant measure
$\mathcal{D}l={\prod}_{i}\frac{d{l}_{i}}{{l}_{i}}{F}_{\epsilon}\left(l\right)$
, where
${F}_{\epsilon}$
enforces an ultraviolet cutoff
$\epsilon $
. They investigated the average curvature
$\mathcal{\mathcal{R}}\sim \langle R\rangle $
and squared curvature
${\mathcal{\mathcal{R}}}^{2}$
(scaled by powers of
$\langle {l}^{2}\rangle $
to make them dimensionless), as well as
$\langle {\delta}_{b}^{2}\rangle $
and
$\langle {V}_{b}\rangle /\langle {l}^{2}\rangle $
. For
$\lambda =k=a=0$
, one finds a negative
$\mathcal{\mathcal{R}}$
and a large
${\mathcal{\mathcal{R}}}^{2}$
, indicating a rough geometry. For small
$a$
, one observes a sudden decrease in
${\mathcal{\mathcal{R}}}^{2}$
, as well as a jump from large positive to small negative values of
$\mathcal{\mathcal{R}}$
as
$\lambda $
is increased. For large
$a$
,
$\mathcal{\mathcal{R}}$
is small and negative, and the geometry appears to be smooth. Like Berg, they advocated a fundamentallength scenario, where the dynamically determined average link length provides an effective UVcutoff.
3.4 The phase structure
Further evidence for a transition between a region of rough and smooth geometry comes from Monte Carlo simulations by Hamber [
108,
111]
on
${4}^{4}$
and
${8}^{4}$
lattices, this time with the lattice
$\mathcal{D}l={\prod}_{i}{l}_{i}d{l}_{i}$
(see also [
114]
for a summary of results, and [
109]
for more details on the method).
There is a value
${k}_{c}$
at which the average curvature vanishes. For
$k>{k}_{c}$
, the curvature becomes large and the simplices degenerate into configurations with very small volumes. He performed a simultaneous fit for
${k}_{c}$
,
${A}_{\mathcal{\mathcal{R}}}$
and
$\delta $
in the scaling relation
$$\begin{array}{c}\mathcal{\mathcal{R}}\left(k\right)\stackrel{k\to {k}_{c}}{\sim}{A}_{\mathcal{\mathcal{R}}}({k}_{c}k{)}^{\delta}.\end{array}$$ 
(17)

This leads to a scaling exponent
$\delta \approx 0.60$
, with only a weak dependence on
$a$
. There are even points with
$a=0$
that lie in the welldefined, smooth phase. Hamber also investigated the curvature and volume susceptibilities
${\chi}_{\mathcal{\mathcal{R}}}$
and
${\chi}_{V}$
. At a continuous phase transition,
${\chi}_{\mathcal{\mathcal{R}}}$
should diverge, reflecting longrange correlations of a massless graviton excitation. The data obtained are not incompatible with such a scenario, but the extrapolation to the transition point
${k}_{c}$
seems somewhat ambiguous. On the other hand, one does not expect
${\chi}_{V}$
to diverge at
${k}_{c}$
, which is corroborated by the simulations.
3.5 Influence of the measure
For the pure Einstein action, ( 12 ) with
$\lambda =0$
, some differences between the
${M}_{1}=dl/l$
and the
${M}_{2}=ldl$
measure were investigated on a
${4}^{4}$
lattice by Beirl et al [
29,
30]
. A constantvolume constraint was used for simulations with the
${M}_{1}$
measure, and a cutoff
${l}^{2}<const$
for
${M}_{2}$
. A study of the
$k$
dependence of bulk geometric quantities agreed with previous simulations, wherever applicable. For small
$k$
, both measures lead essentially to identical results. For the
${M}_{1}$
measure only, one finds that
$\langle R\rangle /\langle V\rangle $
and
$\langle A\rangle $
exhibit a jump at some
${k}_{0}$
, due to the formation of spikes (isolated long link lengths, with the areas staying small). ¿From this, and the study of edge length distributions, one concludes that the DeWitttype measure
${M}_{2}$
is generally better behaved.
3.6 Evidence for a secondorder transition?
The same
$ldl$
measure was used by Hamber in an extension of a previous simulation of the higherderivative action ( 16 ) [
113]
, for lattice sizes of up to
${16}^{4}$
. At the transition point (
$k={k}_{c}\simeq 0.244$
,
$\lambda =1$
,
$a=0.005$
), the distributions of edge lengths, volumes and curvatures are smooth and Gaussianlike, and the average curvature vanishes. The location of
${k}_{c}\left(a\right)$
from fits ( 17 ) to the average curvature coincide with those from [
111]
, leading to
$\delta \simeq 0.626$
. The data at
$a=0$
do not seem to match this interpretation. This leads to the tentative conclusion that only for sufficiently large
$a$
the observed transition is of second order. It is in general “difficult to entirely exclude the presence of a weak firstorder transition, if it has a very small latent heat”. For
$a=0.005$
, one finds some evidence for a decrease in the fractal dimension as
$k$
grows.
3.7 Avoiding collapse
Starting from the observation that Regge configurations exist with
${R}_{tot}\to \infty $
and
${V}_{tot}\to 0$
, Beirl et al [
31,
32]
investigated the influence of a cutoff
$f$
on the fatness of a simplex
$\sigma $
, defined by
$$\begin{array}{c}{\phi}_{\sigma}:=\frac{{V}_{\sigma}^{2}}{{max}_{l\in \sigma}{l}^{8}}\ge f>0.\end{array}$$ 
(18)

Such a uniform shrinking of simplices is known to be necessary in order for piecewise flat manifolds to approach their continuum counterparts [
82]
. On the
${4}^{4}$
lattice, with the
$ldl$
measure and
$\lambda =0$
, they determined
$\mathcal{\mathcal{R}}\left(k\right)$
for decreasing
$f$
. For small
$k$
, the choice of
$f$
seems to have only little influence, but towards the transition point
${k}_{c}$
,
$\langle {\phi}_{\sigma}\rangle $
simultaneously decreases. Next they investigated a variety of measures of the form
${l}^{2\sigma 1}dl$
(see also [
35]
). For
$\sigma \le 1$
, only a mild
$\sigma $
dependence of
$\mathcal{\mathcal{R}}$
is observed as
$k\to {k}_{c}$
. However, for
$\sigma =1.5$
, there are significant differences for the entire range of
$k$
, and some evidence that the geometry degenerates.
3.8 Twopoint functions
To understand the nature of the possible excitations at the phase transition, one needs to study correlation functions in the vicinity of
${k}_{c}$
, which is difficult numerically. Some data are available on the connected correlation functions of the curvatures and the volumes at fixed geodesic distance
$d$
,
${G}_{R}\left(d\right)$
and
${G}_{V}\left(d\right)$
, for lattice sizes
$\le {16}^{4}$
[
115]
, using a scalar field propagator to determine
$d$
. Both correlators were measured at various
$k$
values, leading to similar results for both
$a=0$
and
$a=0.005$
. The data, taken for
$d\le 16$
(
$\simeq $
7 lattice spacings), can be fitted to decaying exponentials.
Some further data (for
$a=0$
) were reported by the Vienna group [
37,
34]
. These authors simply used the lattice distance
$n$
instead of the true geodesic distance
$d$
. In [
37]
, the measure was taken to be of the form
${l}^{2\sigma 1}dl$
. They looked at
${G}_{V}\left(n\right)$
on
${3}^{3}\times 8$
and
${4}^{3}\times 16$
lattices, for
$\lambda =\sigma =1$
and
$\lambda =\sigma =0.1$
, and found a fast decay for all investigated values of
$k$
, and
$n\le 8$
.
3.9 Nonhypercubic lattices
Simulations on lattices with irregular link geometry (still with
${T}^{4}$
topology) have been performed by Beirl et al [
31,
36]
. They were obtained by adding a few vertices of low coordination number to otherwise regular lattices. One finds that the average curvature
$\langle R\rangle $
increases from negative to positive values, even for
$k<{k}_{c}$
, as a result of the formation of spikes. In [
36]
, the averages
$\langle {l}^{2}\rangle $
were monitored separately at the regular and the inserted vertices.
The dependence of
$\mathcal{\mathcal{R}}$
on
$k$
is rather interesting: One observes two “critical” points, a smaller one
${k}_{1}$
, where the extra vertices develop spikes, and a second one
${k}_{2}$
where the remaining vertices follow.
$\mathcal{\mathcal{R}}$
undergoes a small jump at
${k}_{1}$
, and a larger one at
${k}_{2}$
. There is also a transition point to a phase with collapsed simplices at large negative
$k$
, with a jump to large negative
$\mathcal{\mathcal{R}}$
. Additional transition points at negative
$k$
were also found in simulations of the “compactified” Regge action
$S={\sum}_{b}(k{A}_{b}sin{\delta}_{b}+\lambda {V}_{b})$
(this action was discussed in [
104]
; see also [
71,
140]
) and a
$\mathsf{Z}{\mathsf{Z}}_{2}$
version of Regge gravity [
33]
. Correlation functions at those points were computed in [
34]
for short distances, but no evidence for longrange correlations was found.
The same authors studied the inclusion of the higherderivative term (
15 ) in [
36]
. On the regular lattice, their findings for
$\mathcal{\mathcal{R}}(k,a)$
confirmed those by Hamber, apart from the fact that they found stable expectation values even for positive
$\mathcal{\mathcal{R}}$
. Inserting irregular vertices pushes
$\mathcal{\mathcal{R}}$
to larger values and leads again to the appearance of an additional transition point.
3.10 Coupling to SU(2)gauge fields
Berg and collaborators [
42,
43,
44,
27,
28,
41]
coupled the purecurvature action geometrically to the Wilson action for SU(2)gauge fields via dimensionless weight factors
${W}_{b}$
,
$$\begin{array}{c}S={S}^{Regge}\frac{\beta}{2}{\sum}_{b}{W}_{b}Re[Tr(1{U}_{b}\left)\right],{W}_{b}=const\frac{{V}_{b}}{({A}_{b}{)}^{2}},\end{array}$$ 
(19)

where
${U}_{b}$
denotes the SU(2)holonomy around
$b$
and
$\beta $
is proportional to the inverse square coupling constant,
$\beta =\frac{4}{{g}^{2}}$
. One motivation was to understand whether in the welldefined puregravity region, one can choose the elementary particle masses to be
$\ll {m}_{Pl}$
as
$\beta \to \infty $
, as one might expect for a realistic gravity+matter system. This seems a rather distant hope, since in the simulations performed so far, the ratio
${m}_{Pl}/{m}_{hadron}$
is of order unity.
Initial computations were performed on a
$2\times {4}^{3}$
lattice with the scaleinvariant measure, and at
$k=0.01$
[
42,
44]
, and extended to larger
$k$
values in [
43]
. For
$k\le 0.04$
, one finds some evidence for a (firstorder?) transition; the region of
$\beta $
where the transition occurs does not change much with
$k$
. Beirl et al [
27,
28]
extended this analysis by measuring the static potential
$V$
of a quarkantiquark pair on lattices of size
$4\times {6}^{3}$
and
$4\times {8}^{3}$
. With and without gravity, one finds both a confined and a deconfined phase; in the presence of gravity, the transition occurs at a smaller
$\beta $
value. More recently, Berg et al [
41]
have gathered further data on the location and stability of the welldefined phase in the
$(k,\beta )$
plane, and extracted a string tension for various
$\beta $
values.
3.11 Coupling to scalar matter
The coupling of Regge gravity to scalar matter was investigated by Hamber and Williams [
120,
116]
for
$\lambda =1$
,
$a=0.005$
, and various
$k$
values. They considered a single, massive scalar field, with an action contribution
$$\begin{array}{c}{S}_{\phi}=\frac{1}{2}{\sum}_{<ij>}{V}_{ij}{\left(\frac{{\phi}_{i}{\phi}_{j}}{{l}_{ij}}\right)}^{2}+\frac{{m}^{2}}{2}{\sum}_{i}{V}_{i}{\phi}_{i}^{2},\end{array}$$ 
(20)

where
${V}_{i}$
,
${V}_{ij}$
denote the 4volumes associated with the vertices and edges. Most data were obtained on the
${4}^{4}$
lattice, with the measure
${V}_{l}^{(1/30)}ldl$
. The inclusion of the scalar field leads to a slight decrease in the average volume and edge length, as well as a very slight decrease in
$\mathcal{\mathcal{R}}$
.
They again performed a fit for
$\mathcal{\mathcal{R}}$
according to ( 17 ) and found that the critical exponent and the location of
${k}_{c}$
were almost unchanged.
3.12 Recovering the Newtonian potential
Hamber and Williams [
121]
used an ansatz coming from a weakfield calculation for the amplitude of two gravitationally interacting particles separated by a fixed geodesic distance. They used the
$ldl$
measure on a
${16}^{4}$
lattice, with
$\lambda =1$
and
$a=0$
. The potential
$V\left(d\right)$
is extracted from the connected Wilsonline to Wilsonline correlator. ¿From data in
$d\in [0,5]$
, and for sufficiently large
$k$
, one finds an attractive potential
$V\left(d\right)<0$
. They considered fits to various shapes of the potential
$V$
, in order to extract an effective Newton's constant. – They also suggested an effective action for
$\mathcal{\mathcal{R}}$
depending on
$k,\lambda $
and
$a$
that shares some of the features observed in the simulations.
3.13 Gauge invariance in Regge calculus?
A view that has been expressed frequently is that away from configurations with special symmetries, different edge length assignments correspond to inequivalent geometries, and in this sense Regge calculus possesses no gauge invariance ([
187,
81,
124]
; for a dissenting opinion, see [
184]
).
The diffeomorphism invariance can be recovered in a weakfield perturbation about flat space, as was shown by Roček and Williams [
181,
183]
, and there is some evidence for the existence of analogous zeromodes in perturbations of regular, nonflat tesselations, at least in 2d [
122]
. Inspired by the perturbative analysis around flat space, Hamber and Williams [
122]
argue that a similar gauge invariance should persist even if one perturbs around an arbitrary nonflat background.
They propose as a possible definition for such gauge transformations local variations
$\delta {l}_{i}^{2}$
of the link lengths that leave both the local volume and the local curvature terms invariant.
One may hope that in the nonperturbative Regge regime no gaugefixing is necessary, since the contributions from zeromodes cancel out in the pathintegral representation for operator averages [
124,
122]
. Menotti and Peirano [
162,
164,
163,
161]
, following a strategy suggested by Jevicki and Ninomiya [
134]
, have argued vigorously that the functional integral should contain a nontrivial FaddeevPopov determinant. Their starting point is somewhat different from that adopted in the pathintegral simulations (see also [
199]
). They treat piecewise flat spaces as special cases of differentiable manifolds (with singular metric), with the action of the full diffeomorphism group still welldefined. To arrive at a concrete representation for the FaddeevPopov term which could be used in simulations seems at present out of reach.
Recently, Hamber and Williams [
122,
123]
have argued that the
$ldl$
lattice measure is the essentially unique local lattice measure over squared edge lengths (this is a special case of the oneparameter family of local measures of the form
${\prod}_{\sigma}\left[V\right(\sigma ){]}^{\nu}{\prod}_{l}ldl$
; see also [
24]
for a related derivation). It does of course require a term with positive cosmological term in the action in order to suppress long edge lengths.
3.14 Assorted topics
Hartle [
124,
125]
has suggested computing the wave functional of the universe in a simplicial approximation, and evaluating the discrete path integral semiclassically near stationary points of the Regge action. He investigated numerically the extrema of the action ( 12 ) on small simplicial manifolds with topologies
${S}^{4}$
,
$C{P}^{2}$
and
${S}^{2}\times {S}^{2}$
[
126]
. The properties of a HartleHawking wave functional for a small complex with an
${S}^{3}$
boundary were studied in [
127]
.
Fröhlich [
104]
has advocated the need for a proof of reflection positivity of the Regge path integral, which one may expect to play a role in proving the unitarity of the theory. This can be formulated as a condition on the pathintegral measure (including the action) under the gluing of two simplicial fourmanifolds along a threedimensional boundary.
Other authors have suggested associating gaugetheoretic instead of metric variables with the building blocks of a simplicial complex, for the case of the Poincaré group [
71]
, the Lorentz group [
140]
, and for Ashtekar gravity with gauge group
$SU(2,C)$
[
131,
132]
, and reformulating the quantum theory in terms of them.
Hamiltonian 3+1 versions of Regge calculus have been studied classically (see [
190]
for a review), but attempts to quantize them have not progressed very far. One meets problems with the definition of the constraints and the (non)closure of their Poisson algebra. A recent proposal for constructing a canonical quantum theory is due to Mäkelä [
158]
, who constructed a simplicial version of the WheelerDeWitt equation, based on the use of area instead of length variables (which however are known to be overcomplete). In a similar vein, Khatsymovsky [
141]
has suggested that the operators measuring spatial areas ought to have a discrete spectrum.
3.15 Summary
Quantum Regge calculus is based on the wellexplored classical discretization of the Einstein action due to Regge. Its weakfield limit around flat space agrees with the continuum result.
Numerical simulations of the Euclidean path integral indicate the existence of a welldefined phase with small (negative) average curvature for sufficiently small
$k$
and sufficiently large
$\lambda $
, even in the absence of higherorder curvature terms. Hamber [
111,
113]
has found some evidence for a secondorder phase transition in the presence of a small higherorder derivative coupling, with a vanishing average curvature at the transition point. These findings have not been confirmed by other groups. The recent controversy in the dynamical triangulations approach teaches us to treat this issue with some caution.
Almost all simulations have been done on hypercubic, subdivided lattices with
${T}^{4}$
topology, which may introduce a systematic bias in the results. There is evidence that the choice of measure plays a role in the appearance and suppression of singular geometries, socalled spikes. The study of irregular lattices suggests a direct link of the transition points with the appearance of such spikes.
This feature is reminiscent of the appearance of singular structures in dynamical triangulations.
The coupling of a single scalar or SU(2)gauge field seems to have little influence on the phase structure of the gravitational sector.
4 Dynamical triangulations
4.1 Introduction
This quantization approach has received a lot of attention since the early nineties [
2,
1,
11]
, inspired by analogous studies in twodimensional gravity, where dynamical triangulation methods have been a valuable tool in complementing analytical results (see, for example, [
85,
5]
). I will here exclusively concentrate on the 4d results. Other overview material is contained in [
6,
136,
195,
4,
3,
17,
59,
73,
144,
135,
8,
9,
54]
.
Dynamical triangulations are a variant of quantum Regge calculus, where the dynamical variables are not the edge lengths of a given simplicial complex, but its connectivity. A precursor is Weingarten's [
193]
prescription for computing transition amplitudes between threegeometries, by summing over all interpolating fourgeometries, built from equilateral 4d hypercubes living on an imbedding
$p$
dimensional hypercubic lattice with lattice spacing
$a$
. Evaluating the Einstein action on such a configuration amounts to a simple counting of hypercubes of dimension 2 and 4, c.f.
(
21 ).
To avoid a potential overcounting in the usual Regge calculus, Römer and Zähringer [
184]
proposed a gaugefixing procedure for Regge geometries. They argued for an essentially unique association of Riemannian manifolds and equilateral triangulations that in a certain sense are best approximations to the continuum manifolds. The resulting “rigid Regge calculus” is essentially the same structure that nowadays goes by the name of “dynamical triangulations”. In this ansatz one studies the statistical mechanical ensemble of triangulated fourmanifolds with fixed edge lengths, weighted by the Euclideanized Regge action, with a cosmological constant term, and optionally higherderivative contributions. Each configuration represents a discrete geometry, i.e. the discrete version of a Riemannian fourmetric modulo diffeomorphisms. At least for fixed total volume, the state sum converges for appropriate values of the bare coupling constants, if one restricts the topology (usually to that of a sphere
${S}^{4}$
).
4.2 Path integral for dynamical triangulations
Denoting by
$\mathcal{T}$
the set of all triangulations of the foursphere, the partition function for the model is given by
$$\begin{array}{c}Z({\kappa}_{2},{\kappa}_{4})={\sum}_{T\in \mathcal{T}}\frac{1}{C\left(T\right)}{e}^{S\left[T\right]},S\left[T\right]={\kappa}_{2}{N}_{2}\left(T\right)+{\kappa}_{4}{N}_{4}\left(T\right),\end{array}$$ 
(21)

where
${N}_{2}$
and
${N}_{4}$
denote the numbers of 2and 4simplices contained in the simplicial manifold
$T$
, and
$C\left(T\right)$
is the order of the automorphism group of
$T$
. One may think of ( 21 ) as a grand canonical ensemble, with chemical potential
${\kappa}_{4}$
. It is related to the canonical ensemble with fixed volume,
$Z({\kappa}_{2},{N}_{4})$
, by a Legendre transform
$$\begin{array}{c}Z({\kappa}_{2},{\kappa}_{4})={\sum}_{{N}_{4}}{e}^{{\kappa}_{4}{N}_{4}}Z({\kappa}_{2},{N}_{4})\equiv {\sum}_{{N}_{4}}{e}^{{\kappa}_{4}{N}_{4}}{\sum}_{T\in \mathcal{T}\left({N}_{4}\right)}{e}^{{\kappa}_{2}{N}_{2}\left(T\right)}.\end{array}$$ 
(22)

The metric information is encoded in the connectivity of the simplicial decomposition, since the individual 4simplices are assumed equilateral, with the edge length
$a$
set to 1.
To understand the simple form of the action
$S$
, recall that the curvature term in Regge calculus (c.f.( 12 ) is represented by
${\sum}_{b}2{\delta}_{b}{A}_{b}$
, which for fixed edge length is proportional to
$({c}_{4}{N}_{2}10{N}_{4})$
.
The constant
${c}_{4}=2\pi /arccos\frac{1}{4}=4.767$
is determined from the condition that a triangulation of flat space should have average vanishing curvature [
2,
11]
. (Because the foursimplices
$\sigma $
are equilateral, zero curvature can only be achieved upon averaging. This explains the absence of a conventional perturbation theory around flat space.) The cosmological term is represented by
$\lambda {N}_{4}V\left(\sigma \right)\sim \lambda {N}_{4}$
.
It is sometimes convenient to reexpress
${N}_{2}$
as a function of
${N}_{0}$
, using
${N}_{2}=2{N}_{0}+2{N}_{4}4$
, valid for the
${S}^{4}$
topology. The corresponding partition function is
$Z({\kappa}_{0},{\kappa}_{4})$
(where
${\kappa}_{0}=2{\kappa}_{2}$
).
4.3 Existence of an exponential bound?
As a consequence of identities and inequalities satisfied by the
${N}_{i}$
([
8]
contains a detailed discussion), the action ( 21 ) is bounded above and below for fixed volume
${N}_{4}$
. If the number of configurations for fixed
${N}_{4}$
, is exponentially bounded as
${N}_{4}\to \infty $
, that is,
$Z({\kappa}_{2},{N}_{4})$
grows at most as
$Z({\kappa}_{2},{N}_{4})\sim {e}^{const{N}_{4}}$
, there is a “critical line”
${\kappa}_{4}={\kappa}_{4}^{c}\left({\kappa}_{2}\right)$
in the
$({\kappa}_{2},{\kappa}_{4})$
plane, where for fixed
${\kappa}_{2}$
,
$Z({\kappa}_{2},{\kappa}_{4})$
converges for
${\kappa}_{4}>{\kappa}_{4}^{c}\left({\kappa}_{2}\right)$
. True critical behaviour may be found by approaching suitable points on this line from the region above the line, where
$Z$
is well defined.
Doubts on the existence of an exponential bound were raised by Catterall et al [
74]
, who considered the behaviour of
$\Omega $
in
$Z({\kappa}_{0},{\kappa}_{4})={\sum}_{{N}_{4}}{e}^{{\kappa}_{4}{N}_{4}}\Omega ({\kappa}_{0},{N}_{4})$
. Their data (taken for
${N}_{4}\le 32k$
) were consistent with a leading factorial behaviour
$\Omega \sim ({N}_{4}!{)}^{\delta}$
. The same scenario was favoured by de Bakker and Smit [
89]
, who performed further investigations of
${\kappa}_{4}^{c}$
. Subsequently, Ambjørn and Jurkiewicz [
12]
and Brügmann and Marinari [
60]
added further data points at
${N}_{4}=64k$
and
${N}_{4}=128k$
respectively. Their numerical results, as well as those by Catterall et al [
78]
, who employed an alternative method for measuring
$\Omega $
, favour the existence of an exponential bound, although they cannot claim to be conclusive.
There have also been theoretical arguments for the existence of an exponential bound, based on the proofs of such bounds for the counting of minimal geodesic ball coverings of Riemannian spaces of bounded geometry [
68,
26]
, and the counting of discrete curvature assignments to unordered sets of bones [
8]
.
4.4 Performing the state sum
The partition function is evaluated numerically with the help of a MonteCarlo algorithm (see [
57,
50,
72]
for details). There is a set of five topologypreserving moves which change a triangulation locally, and which are ergodic in the grand canonical ensemble [
171,
172,
105]
(see also [
8]
for a discussion).
No ergodic finite set of topology
and volumepreserving moves exists for generic fourdimensional manifolds. This prevents one from using the canonical (i.e. volumepreserving) ensemble. If
${S}^{4}$
is not algorithmically recognizable in the class of all piecewise linear (or smooth) 4d manifolds, the numerical simulations may miss out a substantial part of the state space because of the absence of “computational ergodicity” [
170]
.
An attempt was made by Ambjørn and Jurkiewicz [
13]
to link the nonrecognizability to the presence of large
${N}_{4}$
barriers, which should manifest themselves as an obstacle to cooling down an large initial random triangulation to the minimal
${S}^{4}$
configuration. No such barriers were found for system sizes
$\le 64k$
, but unfortunately they were also absent for an analogous simulation (for
${N}_{5}\le 32k$
) performed by de Bakker [
86]
for
${S}^{5}$
, which is not recognizable.
Since the local moves alter the volume, one works in practice with a “quasicanonical” ensemble, i.e. one uses the grand canonical ensemble
$Z({\kappa}_{2},{\kappa}_{4})$
, but adds a potential term to the action so that the only relevant contributions come from states in an interval
$[V\Delta {N}_{4},V+\Delta {N}_{4}]$
around the target volume
$V$
. There have been several crosschecks which have found no dependence of the results on the width and shape of the potential term [
79,
51]
, but the lattice sizes and fluctuations may still be too small to detect a potential failure of ergodicity, c.f. [
53]
.
To improve the efficiency of the algorithm, Ambjørn and Jurkiewicz [
14]
used additional global (topologypreserving) “baby universe surgery” moves, by cutting and gluing pieces of the simplicial complex . In the branched polymer phase, one can estimate the entropy exponent
$\gamma $
, assuming a behaviour of the form
$Z({\kappa}_{2},{N}_{4})={N}_{4}^{\gamma \left({\kappa}_{2}\right)3}{e}^{{\kappa}_{4}^{c}\left({\kappa}_{2}\right){N}_{4}}\times (1+O(1/{N}_{4}\left)\right)$
, by counting baby universes of various sizes. At the transition point, one finds
$\gamma \left({\kappa}_{2}^{c}\right)\approx 0$
[
10,
14,
97]
. More recently, Egawa et al [
98]
have reported a value of
$\gamma \approx 0.26$
.
4.5 The phase structure
Already the first simulations by Agishtein and Migdal [
2,
1]
and Ambjørn and Jurkiewicz [
11]
of dynamically triangulated 4d gravity exhibited a clear twophase structure. After tuning to the infinitevolume or critical line
${\kappa}_{4}^{c}\left({\kappa}_{2}\right)$
, one identifies two regions,
${\kappa}_{2}<{\kappa}_{2}^{c}$
and
${\kappa}_{2}>{\kappa}_{2}^{c}$
. The critical value
${\kappa}_{2}^{c}$
depends on the volume
${N}_{4}$
, and it was conjectured that it may move out to
$+\infty $
as
${N}_{4}\to \infty $
[
75,
89]
, but it was later shown to converge to a finite value [
14]
. One can characterize the region with small
${\kappa}_{2}<{\kappa}_{2}^{c}$
as the hot, crumpled, or condensed phase. It has small negative or positive curvature, large (possibly infinite) Hausdorff dimension
${d}_{H}$
and a high connectivity. By contrast, for
${\kappa}_{2}>{\kappa}_{2}^{c}$
one is in the cold, extended, elongated, or fluid phase. It has large positive curvature, with an effective treelike branchedpolymer geometry, and
${d}_{H}\approx 2$
.
The location of the critical point on the infinitevolume line may be estimated from the peak in the curvature susceptibility
$\chi ({\kappa}_{2},{N}_{4})=(\langle {N}_{2}^{2}\rangle \langle {N}_{2}{\rangle}^{2})/{N}_{4}$
[
2,
1,
11,
192,
129,
90]
, or the node susceptibility
$\chi ({\kappa}_{0},{N}_{4})=(\langle {N}_{0}^{2}\rangle \langle {N}_{0}{\rangle}^{2})/{N}_{4}$
[
75,
76]
, as well as higher cumulants of
${N}_{0}$
[
48]
. In [
14]
, it was suggested that one may alternatively estimate
${\kappa}_{2}^{c}\left({N}_{4}\right)$
by looking at the behaviour of the entropy exponent
$\gamma $
, approaching
${\kappa}_{2}^{c}$
from the elongated phase. More recently, Catterall et al [
77]
have used the fluctuations
${\chi}_{0}=(\langle {\omega}_{0}^{2}\rangle \langle {\omega}_{0}{\rangle}^{2})/{N}_{4}$
of the local volume
${\omega}_{0}$
around singular vertices as an order parameter.
4.6 Evidence for a secondorder transition?
It sometimes seems to be assumed that if one were to find a continuum theory at a secondorder phase transition, it would have flat Minkowski space as its ground state (in spite of the
${S}^{4}$
topology), and gravitonic spin2 excitations. An alternative scenario with a constantcurvature spheremetric has been put forward by de Bakker and Smit [
90,
92]
.
Agishtein and Migdal [
2]
initially reported a hysteresis in the average curvature
$\langle R\rangle \left({\kappa}_{2}\right)$
, indicating a firstorder transition. However, subsequent authors [
11,
191,
192,
75,
76]
found numerical data not incompatible with the existence of a secondorder transition, and also Agishtein and Migdal [
1]
retracted their original claim as a result of a closer examination of the fixed point region.
One tries to discriminate between a firstand second(or higher)order transition by looking at the Binder parameter [
191,
192,
7]
, or scaling exponents
$\alpha $
governing the scaling behaviour
$\sim {\kappa}_{2}^{c}{\kappa}_{2}{}^{\alpha}$
of suitable observables [
1,
192]
, or the peak height of susceptibilities as a function of the volume
${N}_{4}$
[
16,
75,
76,
14,
48,
87]
. Other scaling relations are discussed in [
90,
91,
96,
95]
. However, since critical parameters are hard to measure, and it is difficult to estimate finitesize effects, none of the data can claim to be conclusive.
Some doubts were cast on the conjectured continuous nature of the phase transition by Bialas et al [
48]
, who found an unexpected twopeak structure in the distribution of nodes near the fixed point. This was strengthened further by data taken at 64k by de Bakker [
87]
, with an even more pronounced double peak (see also [
51]
). Most likely previous simulations were simply too small to detect the true nature of the phase transition. Both Bialas et al and de Bakker observed that the finite size scaling exponents extracted from the node susceptibility
${\chi}_{0}$
grow with volume and may well reach the value 1 expected for a firstorder transition as
$V\to \infty $
.
The origin of this behaviour was further elucidated by Bialas et al [
49,
47]
and Bialas and Burda [
46]
, who found a simple meanfield model that reproduces qualitatively the phase structure of 4d dynamically triangulated quantum gravity. With an appropriate choice of local weights, this model has a condensed and a fluid phase, with a firstorder transition in between. A similar behaviour was found by Catterall et al [
77]
, who made a related meanfield ansatz, with the local weights depending on the local entropies around the vertices.
4.7 Influence of the measure
The inclusion of a local measure term corresponding to
$(detg{)}^{\frac{\mu}{2}}$
was studied by Brügmann and Marinari [
56,
57,
58]
. This amounts to adding to the action a term of the form
$\mu {\sum}_{v}lno\left(v\right)/5$
, where
$o\left(v\right)$
is the number of 4simplices containing a vertex
$v$
. For
$\mu =5,1,1,5$
, this seems to lead merely to a shift of the critical line. This setting has been revived recently by Renken [
178,
177]
in the context of a renormalization group analysis.
4.8 Higherderivative terms
Higherorder derivative terms were considered by Ambjørn et al [
16]
(see also [
143]
), who added a term of the form
$$\begin{array}{c}h{c}_{4}^{2}{\sum}_{b}o\left(b\right){\left(\frac{{c}_{4}o\left(b\right)}{o\left(b\right)}\right)}^{2}\end{array}$$ 
(23)

to the action, where again
${c}_{4}=4.767$
and
$o\left(b\right)$
denotes the number of 4simplices sharing a bone
$b$
. For
$h\in [0,20]$
, and with volumes up to
$32k$
, no major qualitative changes of the geometrical observables were found. The inclusion of the higherderivative term also does not improve the behaviour of the average curvature
$\langle R\rangle $
, which continues to be positive at the critical point, whereas from a naïve comparison with the continuum theory one would expect it to scale to zero. (This is also incompatible with the prediction of Antoniadis et al [
20]
, should dynamical triangulations possess an infrared stable fixed point.) De Bakker and Smit [
89]
have argued that this may not be a reason of concern, since one expects the volume and curvature terms to mix under renormalization.
4.9 Coupling to matter fields
Ambjørn et al [
7]
considered the influence of both Ising spins and Gaussian scalar fields on bulk geometric quantities. The phase structure remained essentially unchanged, and no improvement in the scaling behaviour of
$\langle R\rangle $
was found. Taken together with their results on higherderivative gravity [
16]
, one finds a universal linear dependence of the cosmological constant
${\kappa}_{4}^{c}\left({\kappa}_{2}\right)$
, with slope
$\sim 2.5$
.
Coupling to
$\mathsf{Z}{\mathsf{Z}}_{2}$
spin variables
$s\left(l\right)$
located on edges
$l$
was considered by Ambjørn et al [
15]
, who added a Wilson loop term
${S}_{W}[s;T]=\beta {\sum}_{b\in T}o\left(b\right){\prod}_{l\in b}s\left(l\right)$
to the action. The matter sector behaves largely as expected when
${\kappa}_{2}$
is varied between the crumpled and the elongated gravity phase. However, in the common critical region of both sectors, where a priori one might have expected interesting effects, the critical behaviour seems to agree with that of the pure gravity system.
More recently, Bilke et al [
52]
have reported a nontrivial backreaction of matter on geometry, when considering coupling to several noncompact
$U\left(1\right)$
gauge fields. Their study was in part motivated by a continuum analysis of the dynamics of the conformal factor of Antoniadis et al [
21]
. Including three fields
$U\left(1\right)$
fields seems to lead to a total suppression of the branched polymer phase, which is replaced by a new weakcoupling phase with negative susceptibility exponent
$\gamma $
and a fractal dimension
$\approx 4$
. These are clearly interesting results, but should be treated with some caution because of the small lattice sizes involved (
${N}_{4}\le 16k$
).
4.10 Nonspherical lattices
Alternative topologies for the underlying simplicial complex were considered by Bilke et al [
51,
53]
. They showed that for topologies
${S}^{1}\times {S}^{3}$
and
${T}^{4}$
, the free energy agrees to leading order with that of
${S}^{4}$
. Neither the critical value
${\kappa}_{4}^{c}({\kappa}_{2},{N}_{4})$
, nor the leading contributions to the action density
$\langle {N}_{2}\rangle /{N}_{4}$
are changed. However, the simulations become more involved, since the size of the minimal configurations is increased.
4.11 Singular configurations
The singular nature of the geometry in the phase below the critical value
${\kappa}_{2}^{c}$
can be quantified by the distribution
$\rho \left(n\right)$
of the vertex order
$o\left(v\right)$
,
$$\begin{array}{c}\rho \left(n\right)=\frac{1}{{N}_{0}}<{\sum}_{v}{\delta}_{o\left(v\right),n}>,\end{array}$$ 
(24)

first considered by Hotta et al [
129]
. It has a continuum part and a separate peak at high vertex order. One can suppress this effect by adding a term
$\sim {\sum}_{v}\left(o\right(v)5{)}^{2}$
to the action (see also [
19]
for a discussion of terms of a similar nature), but this leads to a simultaneous disappearance of the phase transition. Hotta et al [
130]
have checked that for a variety of initial configurations the singular structure is a generic feature of the model.
Catterall et al [
79,
80]
observed that the pair of singular vertices form the end points of a singular link. They also offered a possible explanation for the formation of these singular structures:
Simplices of sufficiently low dimension can maximize their local entropy by acquiring large local volumes (see also [
49]
for a meanfield argument). Catterall et al [
77]
found two pseudocritical points,
${\kappa}_{0}^{\left(1\right)}$
and
${\kappa}_{0}^{\left(2\right)}$
, associated with the creation of singular vertices and links, which seem to merge into a single critical point
${\kappa}_{0}^{c}$
as
$V\to \infty $
. One concludes that the observed phase transition in the 4d dynamical triangulations model is driven by the appearance and disappearance of singular geometries.
4.12 Renormalization group
There have been attempts to apply renormalization group techniques, assuming that the transition is indeed continuous. Burda et al [
61]
and later Bialas et al [
48]
considered the “cutting of the last generation of minimalneck baby universes”. This step can only be performed once, which severely limits the power of the procedure. Renken [
178,
177]
has applied a different blocking move, involving node deletion, and studied the RG flow using the volume and the vertex order as observables.
4.13 Exploring geometric properties
Since the effective geometry in both phases is rather singular, different ways of measuring length may lead to inequivalent definitions of “dimension”. A common (local) notion is derived from the volume of a geodesic ball with radius
$r$
. Usually the radius is measured in terms of the geodesic distance
${d}_{1}$
(the minimal number of links), or dual geodesic distance
${d}_{4}$
(the minimal number of links of the dual graph). Alternatively, one may consider the number
$n\left(r\right)$
of 4simplices in spherical shells of thickness 1 at distance
$r$
, and define a fractal dimension by
$\langle n\left(r\right)\rangle \sim {r}^{{d}_{F}1}$
[
90,
14]
.
To extract a global length scale, one may use the averages
$\langle {d}_{1}\rangle $
,
$\langle {d}_{4}\rangle $
, or consider the average “radius of the universe”
$\langle r\left(T\right)\rangle $
[
11]
to obtain a cosmological Hausdorff dimension
${d}_{CH}$
, or the “average intrinsic linear extent”
$L=\langle {V}^{2}{\sum}_{i,j}{d}_{4}(i,j,T)\rangle $
(see, for example, [
75,
76]
). Away from the phase transition, the fractal dimensions associated with these geometric construction are more or less equivalent and give
${d}_{F}=\infty $
in the crumpled phase and
${d}_{F}\approx 2$
in the elongated phase. It is difficult to measure the dimension close to the transition point.
It will not be straightforward to interpret the behaviour of observables (defined in analogy with the continuum theory), since in most of the phase space the geometry of the simplicial complex is far from approximating a metric 4manifold. In search of a semiclassical interpretation for geometric observables, an alternative notion of local curvature for a simplicial manifold was suggested by de Bakker and Smit [
90]
, based on a continuum expansion of the volume of a geodesic ball. Assuming furthermore that independent of
$n$
,
$n$
volumes of balls with radius
$r$
behave like regions on
${S}^{n}\subset I{R}^{n+1}$
, they extracted scaling relations for various geometric quantities for an intermediate range for
$r$
. This line of thought was pursued further in [
185]
.
Close to the phase transition, one may investigate the behaviour of test particles (ignoring backreactions on the geometry). Comparing the mass extracted from the oneparticle propagator with the energy of the combined system obtained from the twoparticle propagator [
88,
92,
93]
, one does indeed find evidence for gravitational binding.
4.14 Twopoint functions
It is possible to define twopoint correlation functions on random geometries [
16,
14]
, which are to be thought of as the discrete analogues of formal continuum correlators
$$\begin{array}{c}{G}_{\mathcal{O}}\left(r\right)={\int}_{\frac{Riem{S}^{4}}{Diff{S}^{4}}}D{g}_{\mu \nu}{e}^{S}\int {d}^{4}x\int {d}^{4}{x}^{\prime}\sqrt{detg\left(x\right)}\sqrt{detg\left({x}^{\prime}\right)}\mathcal{O}\left(x\right)\mathcal{O}\left({x}^{\prime}\right)\delta \left({d}_{g}\right(x,{x}^{\prime})r),\end{array}$$ 
(25)

for local observables
$\mathcal{O}\left(x\right)$
, with
${d}_{g}$
denoting the geodesic distance with respect to the metric
${g}_{\mu \nu}$
. There is an ambiguity in defining the connected part of the correlator ( 25 ), as was pointed out by de Bakker and Smit [
91,
92]
. Contrary to expectations, after subtraction of the square of the curvature expectation value, the resulting quantity
$\langle RR\rangle \left(r\right)\langle R{\rangle}^{2}$
does not scale to zero with large distances. They therefore proposed an alternative definition of the connected twopoint function, by subtracting the square of a “curvaturetonothing” correlator
$\langle R\rangle \left(r\right)$
. This definition was compared in more detail by Bialas [
45]
with a more conventional notion, as, for example, the one used in [
48]
. For the case of curvature correlators, their behaviour differs significantly, especially at short distances.
4.15 Summary
In the dynamical triangulations approach, one studies the properties of a statistical ensemble of simplicial fourgeometries à la Regge of fixed edge lengths. By summing over such discrete configurations according to ( 21 ), one has implicitly assumed that this leads to a uniform sampling of the space of smooth Riemannian manifolds. There is no obvious weakfield limit, but this is no obstacle in principle to the pathintegral construction. Numerical simulations indicate the existence of a welldefined phase for sufficiently small
${\kappa}_{2}$
(inverse Newton's constant) and a sufficiently large cosmological constant. For small
${\kappa}_{2}$
, one finds a “crumpled” phase, with small average curvature and a large Hausdorff dimension, and for large
${\kappa}_{2}$
an elongated, effectively twodimensional polymer phase. At present, the consensus seems to be that the corresponding phase transition is of first order, with a finite average curvature at the transition point.
Almost all simulations have been done on simplicial manifolds with
${S}^{4}$
topology. Neither the inclusion of factors of
$(detg{)}^{\rho}$
in the measure nor the addition of higherorder curvature terms to the action seem to have a substantial influence on the phase structure. Also matter coupling to spinorial and scalar fields does not seem to lead to a change of universality class, although the inclusion of several gauge fields may have a more drastic effect. The study of singular structures (vertices of high coordination number) has led to a qualitative understanding of the phase structure of the model.
5 Conclusions and Outlook
In this review, I have collected a variety of results on discrete fourdimensional models of quantum gravity, mainly coming from Euclidean pathintegral approaches. Numerical simulations have yielded information on the phase structure of these models, the behaviour of twopoint functions and a number of other properties of their partition functions. All of the pathintegral models have some qualitative features in common. They need a (sufficiently large, positive) cosmological constant
$\lambda $
to be well defined. For sufficiently small values of Newton's constant
$G$
, one finds a phase of collapsed geometry, with effective dimension
$<4$
. In the gaugetheoretic model, the metric is degenerate; in Regge calculus, one finds spiky configurations; and in dynamical triangulations, the ensemble behaves like that of a branched polymer.
In all cases, one observes a transition on the boundary of this phase, but so far no convincing evidence of longrange correlations has been found in its vicinity. Within the accuracy of the numerical simulations, this main conclusion is not altered by the inclusion of determinantal factors
$(detg{)}^{\rho}$
in the measure, the inclusion of higherorder derivative terms, or the addition of matter fields. Why does this happen? Each of the approaches can claim that its state space represents, at least roughly, an approximation to the space of smooth Riemannian metrics or geometries. This leaves only the pathintegral measure as a possible culprit. The measures used up to now were the simplest ones compatible with considerations of locality and gaugeinvariance. It seems premature to blame the absence of diffeomorphism invariance (whose status in the gaugetheoretic formulation and the Regge calculus program remains unclear), since the explicitly diffeomorphisminvariant dynamical triangulations approach suffers from similar problems. Further analytical insights are needed to understand which modifications of the measure would make these models more interacting.
There are a number of loop holes which could change the picture just presented. It is possible that adding enough matter of the correct type could have a nontrivial effect, or that Regge calculus with the inclusion of higherorder curvature terms does indeed possess a secondorder phase transition. Since we have very little experience with universality properties of 4d generally covariant theories, it is not
a priori clear whether the choice of measure and the initial restrictions on the lattice geometry can affect the final results.
One may of course take the attitude that something is fundamentally wrong with trying to construct a theory of quantum gravity via a statistical field theory approach, and that a different starting point is needed, an obvious candidate being a nonperturbative theory of superstrings, or of more general extended objects. In any case, these different approaches need not be mutually exclusive, and one may therefore take the results of the discrete approaches presented here as an indication that other attempts of constructing quantum gravity nonperturbatively may run into similar difficulties.
A further unresolved problem is the “analytic continuation” of the pathintegral results to Lorentzian signature. The Hamiltonian ansatz circumvents this problem, and some progress has been made in the Hamiltonian gaugetheoretic discrete approach. Although the kinematical structure is in place and some information on the constraint algebra has been obtained, the physical state space has not yet been identified. Its results are therefore not sufficiently complete to admit comparison with the pathintegral simulations. For the simplicial formulations, only little is known about their canonical counterparts. One would hope that future research will throw further light on these issues.
Acknowledgement. I am indebted to R. Williams, P. Menotti, H. Hamber, W. Beirl and J. Ambjørn for comments and criticism, to E. Schlenk for help with the references, and to R. Helling for assistance in reformatting them.
References

Agishtein, M.E., and Migdal, A.A., “Critical behavior of dynamically triangulated quantum gravity in four dimensions”, Nucl. Phys., 385, 395–412, (1992). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Agishtein, M.E., and Migdal, A.A., “Simulations of fourdimensional simplicial quantum gravity as dynamical triangulation”, Mod. Phys. Lett. A, 7, 1039–1061, (1992).

Ambjørn, J., “Quantum gravity represented as dynamical triangulations”, Class. Quantum Grav., 12, 2079–2134, (1995).

Ambjørn, J., “Recent progress in the theory of random surfaces and simplicial quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 42, 3–16, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Ambjørn, J., “Quantization of geometry”, in David, F., Ginsparg, P., and ZinnJustin, J., eds., Fluctuating Geometries in Statistical Mechanics and Field Theory, Proceedings of the Les Houches Summer School, Session LXII, 2 August 9 September 1994, 77–193, (Elsevier, Amsterdam, Netherlands; New York, U.S.A., 1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Ambjørn, J., Burda, Z., Jurkiewicz, J., and Kristjansen, C.F., “Quantum gravity represented as dynamical triangulations”, Acta Phys. Pol. B, 23, 991–1030, (1992).

Ambjørn, J., Burda, Z., Jurkiewicz, J., and Kristjansen, C.F., “Fourdimensional dynamically triangulated gravity coupled to matter”, Phys. Rev. D, 48, 3695–3703, (1993). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Ambjørn, J., Carfora, M., and Marzuoli, A., The Geometry of Dynamical Triangulations, vol. m50 of Lecture Notes in Physics, (Springer, Berlin, Germany; New York, U.S.A., 1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Ambjørn, J., Durhuus, B., and Jonsson, T., Quantum geometry, (Cambridge University Press, Cambridge, 1997).

Ambjørn, J., Jain, S., Jurkiewicz, J., and Kristjansen, C.F., “Observing 4d baby universes in quantum gravity”, Phys. Lett. B, 305, 208–213, (1993). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Ambjørn, J., and Jurkiewicz, J., “Fourdimensional simplicial quantum gravity”, Phys. Lett. B, 278, 42–50, (1992).

Ambjørn, J., and Jurkiewicz, J., “On the exponential bound in fourdimensional simplicial gravity”, Phys. Lett. B, 335, 355–358, (1994). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Ambjørn, J., and Jurkiewicz, J., “Computational ergodicity of
${S}^{4}$
”, Phys. Lett. B, 345, 435–440, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Ambjørn, J., and Jurkiewicz, J., “Scaling in four dimensional quantum gravity”, Nucl. Phys. B, 451, 643–676, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Ambjørn, J., Jurkiewicz, J., Bilke, S., Burda, Z., and Petersson, B., “Z(2) gauge matter coupled to 4d simplicial quantum gravity”, Mod. Phys. Lett. A, 9, 2527–2542, (1994).

Ambjørn, J., Jurkiewicz, J., and Kristjansen, C.F., “Quantum gravity, dynamical triangulation and higher derivative regularization”, Nucl. Phys. B, 393, 601–632, (1993). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Ambjørn, J., Jurkiewicz, J., and Watabiki, Y., “Dynamical triangulations, a gateway to quantum gravity?”, J. Math. Phys., 36, 6299–6339, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Ambjørn, J., Nielsen, J.L., Rolf, J., and Savvidy, G.K., “Spikes in quantum Regge calculus”, Class. Quantum Grav., 14, 3225–3241, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Ambjørn, J., Savvidy, G.K., and Savvidy, K.G., “Alternative actions for quantum gravity and the intrinsic rigidity of the spacetime”, Nucl. Phys. B, 486, 390–412, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Antoniadis, I., Mazur, P.O., and Mottola, E., “Scaling behavior of quantum fourgeometries”, Phys. Lett. B, 323, 284–291, (1994). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Antoniadis, I., Mazur, P.O., and Mottola, E., “Criticality and scaling in 4d quantum gravity”, Phys. Lett. B, 394, 49–56, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Ashtekar, A., “New variables for classical and quantum gravity”, Phys. Rev. Lett., 57(18), 2244–2247, (1986).

Ashtekar, A., “New Hamiltonian formulation of general relativity”, Phys. Rev. D, 36(6), 1587–1602, (1987).

Bander, M., “Functional measure for lattice gravity”, Phys. Rev. Lett., 57, 1825–1827, (1986).

Barbero G, J.F., “Real Ashtekar variables for Lorentzian signature spacetimes”, Phys. Rev. D, 51, 5507–5510, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Bartocci, C., Bruzzo, U., Carfora, M., and Marzuoli, A., “Entropy of random coverings and 4D quantum gravity”, J. Geom. Phys., 18, 247–294, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Beirl, W., Berg, B.A., Krishnan, B., Markum, H., and Riedler, J., “Static quark potentials in quantum gravity”, Phys. Lett. B, 348, 355–359, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Beirl, W., Berg, B.A., Krishnan, B., Markum, H., and Riedler, J., “SU(2) potentials in quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 42, 707–709, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Beirl, W., Gerstenmayer, E., and Markum, H., “Exploration of simplicial quantum gravity in four dimensions”, Nucl. Phys. B (Proc. Suppl.), 26, 575–577, (1992).

Beirl, W., Gerstenmayer, E., and Markum, H., “Influence of the measure on simplicial quantum gravity in fourdimensions”, Phys. Rev. Lett., 69, 713–716, (1992). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Beirl, W., Gerstenmayer, E., Markum, H., and Riedler, J., “Gravitational action versus entropy on simplicial lattices in fourdimensions”, Nucl. Phys. B (Proc. Suppl.), 30, 764–767, (1993). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Beirl, W., Gerstenmayer, E., Markum, H., and Riedler, J., “The well defined phase of simplicial quantum gravity in fourdimensions”, Phys. Rev. D, 49, 5231–5239, (1994). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Beirl, W., Hauke, A., Homolka, P., Krishnan, B., Kröger, H., Markum, H., and Riedler, J., “The phase structure of pure Regge gravity”, Nucl. Phys. B (Proc. Suppl.), 47, 625–628, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Beirl, W., Hauke, A., Homolka, P., Markum, H., and Riedler, J., “Correlation functions in lattice formulations of quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 53, 735–738, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Beirl, W., Markum, H., and Riedler, J., “On the measure of simplicial quantum gravity in fourdimensions”, in Carr, J., and Perrottet, M., eds., High Energy Physics (HEP 93), Proceedings of the International Europhysics Conference, Marseille, France, 22–28 July 1993, 214–216, (Editions Frontières, GifsurYvette, France, 1994). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Beirl, W., Markum, H., and Riedler, J., “Regge gravity on general triangulations”, Phys. Lett. B, 341, 12–18, (1994). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Beirl, W., Markum, H., and Riedler, J., “Two point functions of fourdimensional simplicial quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 34, 736–738, (1994). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Berg, B.A., “Exploratory numerical study of discrete quantum gravity”, Phys. Rev. Lett., 55, 904–907, (1985).

Berg, B.A., “Entropy versus energy on a fluctuating fourdimensional Regge skeleton”, Phys. Lett. B, 176, 39–44, (1986).

Berg, B.A., “Quantum gravity motivated computer simulations”, in Mitter, H., and Widder, F., eds., Particle Physics and Astrophysics: Current Viewpoints, Proceedings of the XXVII Internationale Universitätswochen für Kernphysik, Schladming, Austria, February 1988, 223–259, (Springer, Berlin, Germany; New York, U.S.A., 1989).

Berg, B.A., Beirl, W., Krishnan, B., Markum, H., and Riedler, J., “Phase diagram of Regge quantum gravity coupled to SU(2) gauge theory”, Phys. Rev. D, 54, 7421–7425, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Berg, B.A., and Krishnan, B., “Asymptotic freedom and ReggeEinstein quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 30, 768–770, (1993).

Berg, B.A., and Krishnan, B., “Phase structure of SU(2) lattice gauge theory with quantum gravity”, Phys. Lett. B, 318, 59–62, (1993). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Berg, B.A., Krishnan, B., and Katoot, M., “Asymptotic freedom and Euclidean quantum gravity”, Nucl. Phys. B, 404, 359–384, (1993). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Bialas, P., “Correlations in fluctuating geometries”, Nucl. Phys. B (Proc. Suppl.), 53, 739–742, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Bialas, P., and Burda, Z., “Collapse of 4d random geometries”, Phys. Lett. B, 416, 281–285, (1998). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Bialas, P., Burda, Z., and Johnston, D.A., “Balls in boxes and quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 63, 763–765, (1998). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Bialas, P., Burda, Z., Krzywicki, A., and Petersson, B., “Focusing on the fixed point of 4d simplicial gravity”, Nucl. Phys. B, 472, 293–308, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Bialas, P., Burda, Z., Petersson, B., and Tabaczek, J., “Appearance of mother universe and singular vertices in random geometries”, Nucl. Phys. B, 495, 463–476, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Bilke, S., Burda, Z., and Jurkiewicz, J., “Simplicial quantum gravity on a computer”, Computer Phys. Commun., 85, 278–292, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Bilke, S., Burda, Z., Krzywicki, A., and Petersson, B., “Phase transition and topology in 4d simplicial gravity”, Nucl. Phys. B (Proc. Suppl.), 53, 743–745, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Bilke, S., Burda, Z., Krzywicki, A., Petersson, B., Tabaczek, J., and Thorleifsson, G., “4d simplicial quantum gravity interacting with gauge matter fields”, Phys. Lett. B, 418, 266–272, (1998). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Bilke, S., Burda, Z., and Petersson, B., “Topology in 4d simplicial quantum gravity”, Phys. Lett. B, 395, 4–9, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Bowick, M., “Random surfaces and lattice gravity”, Nucl. Phys. B (Proc. Suppl.), 63, 77–88, (1998). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Brewin, L., “The Riemann and extrinsic curvature tensors in the Regge calculus”, Class. Quantum Grav., 5, 1193–1203, (1988).

Brügmann, B., “Measure of fourdimensional simplicial quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 30, 760–763, (1993).

Brügmann, B., “Nonuniform measure in fourdimensional simplicial quantum gravity”, Phys. Rev. D, 47, 3330–3338, (1993). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Brügmann, B., and Marinari, E., “4D simplicial quantum gravity with a nontrivial measure”, Phys. Rev. Lett., 70, 1908–1911, (1993). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Brügmann, B., and Marinari, E., “Monte Carlo simulations of 4d simplicial quantum gravity”, J. Math. Phys., 36, 6340–6352, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Brügmann, B., and Marinari, E., “More on the exponential bound of fourdimensional simplicial quantum gravity”, Phys. Lett. B, 349, 35–41, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Burda, Z., Kownacki, J.P., and Krzywicki, A., “Towards a nonperturbative renormalization of Euclidean quantum gravity”, Phys. Lett. B, 356, 466–471, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Caracciolo, S., Menotti, P., and Pelissetto, A., “Lattice supergravity and gravitongravitino doubling”, Nucl. Phys. B, 296, 868–876, (1988).

Caracciolo, S., and Pelissetto, A., “Analysis of the critical behavior or the de Sitter quantum gravity on a hypercubic lattice”, Phys. Lett. B, 193, 237–240, (1987).

Caracciolo, S., and Pelissetto, A., “Nonperturbative lattice gravity”, Nucl. Phys. B (Proc. Suppl.), 4, 78–82, (1987).

Caracciolo, S., and Pelissetto, A., “A numerical investigation about quantum measure in lattice gravity”, Phys. Lett. B, 207, 468–470, (1988).

Caracciolo, S., and Pelissetto, A., “Phases and topological structures of de Sitter lattice gravity”, Nucl. Phys. B, 299, 693–718, (1988).

Caracciolo, S., and Pelissetto, A., “From lattice gauge theory towards gravity”, in Damgaard, P.H., Hüffel, H., and Rosenblum, A., eds., Probabilistic Methods in Quantum Field Theory and Quantum Gravity, Proceedings of the NATO Advanced Research Workshop, held August 21–27, 1989, in Cargèse, France, vol. 224 of NATO ASI Series B, 37–54, (Plenum Press, New York, U.S.A., 1990).

Carfora, M., and Marzuoli, A., “Entropy estimates for simplicial quantum gravity”, J. Geom. Phys., 16, 99–119, (1995).

Caselle, M., D'Adda, A., and Magnea, L., “Doubling of all matter fields coupled with gravity on a lattice”, Phys. Lett. B, 192, 411–414, (1987).

Caselle, M., D'Adda, A., and Magnea, L., “Lattice gravity and supergravity as spontaneously broken gauge theories of the (super)Poincaré group”, Phys. Lett. B, 192, 406–410, (1987).

Caselle, M., D'Adda, A., and Magnea, L., “Regge Calculus as a local theory of the Poincaré group”, Phys. Lett. B, 232, 457–461, (1989).

Catterall, S.M., “Simulations of dynamically triangulated gravity”, Computer Phys. Commun., 87, 409–415, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Catterall, S.M., “Lattice quantum gravity: review and recent developments”, Nucl. Phys. B (Proc. Suppl.), 47, 59–70, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Catterall, S.M., Kogut, J.B., and Renken, R.L., “Is there an exponential bound in fourdimensional simplicial gravity?”, Phys. Rev. Lett., 72, 4062–4065, (1994). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Catterall, S.M., Kogut, J.B., and Renken, R.L., “Phase structure of fourdimensional simplicial quantum gravity”, Phys. Lett. B, 328, 277–283, (1994). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Catterall, S.M., Kogut, J.B., and Renken, R.L., “Simulations of fourdimensional simplicial quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 34, 733–735, (1994).

Catterall, S.M., Kogut, J.B., and Renken, R.L., “Singular structure in 4d simplicial gravity”, Phys. Lett. B, 416, 274–280, (1998). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Catterall, S.M., Kogut, J.B., Renken, R.L., and Thorleifsson, G., “Baby universes in 4d dynamical triangulation”, Phys. Lett. B, 366, 72–76, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Catterall, S.M., Thorleifsson, G., Kogut, J.B., and Renken, R.L., “Singular vertices and the triangulation space of the d sphere”, Nucl. Phys. B, 468, 263–276, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Catterall, S.M., Thorleifsson, G., Kogut, J.B., and Renken, R.L., “Simplicial gravity in dimension greater than two”, Nucl. Phys. B (Proc. Suppl.), 53, 756–759, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Cheeger, J., Müller, W., and Schrader, R., “Lattice gravity or Riemannian structure on piecewise linear spaces”, in Breitenlohner, P., and Dürr, H.P., eds., Unified Theories of Elementary Particles: Critical Assessment and Prospects, Proceedings of the Heisenberg Symposium, held in Munich, July 16–21, 1981, vol. 160 of Lecture Notes in Physics, 176–188, (Springer, Berlin, Germany; New York, U.S.A., 1982).

Cheeger, J., Müller, W., and Schrader, R., “On the curvature of piecewise flat spaces”, Commun. Math. Phys., 92, 405–454, (1984).

Corichi, A., and Zapata, J.A., “On diffeomorphism invariance for lattice theories”, Nucl. Phys. B, 493, 475–490, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Das, A., Kaku, M., and Townsend, P.K., “Lattice formulation of general relativity”, Phys. Lett. B, 81, 11–14, (1979).

David, F., “Simplicial quantum gravity and random lattices”, in Julia, B., and ZinnJustin, J., eds., Gravitation and Quantizations, Proceedings of the Les Houches Summer School, Session LVII, 5 July 1 August 1992, 679–750, (Elsevier, Amsterdam, Netherlands, New York, U.S.A., 1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

de Bakker, B.V., “Absence of barriers in dynamical triangulation”, Phys. Lett. B, 348, 35–38, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

de Bakker, B.V., “Further evidence that the transition of 4d dynamical triangulation is first order”, Phys. Lett. B, 389, 238–242, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

de Bakker, B.V., and Smit, J., “Euclidean gravity attracts”, Nucl. Phys. B (Proc. Suppl.), 34, 739–741, (1994). Related online version (cited on 19 November 1998):
.
☻ open access ✓

de Bakker, B.V., and Smit, J., “Volume dependence of the phase boundary in 4d dynamical riangulation”, Phys. Lett. B, 334, 304–308, (1994). Related online version (cited on 19 November 1998):
.
☻ open access ✓

de Bakker, B.V., and Smit, J., “Curvature and scaling in 4d dynamical triangulation”, Nucl. Phys. B, 439, 239–258, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

de Bakker, B.V., and Smit, J., “Two point functions in 4d dynamical triangulation”, Nucl. Phys. B, 454, 343–356, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

de Bakker, B.V., and Smit, J., “Correlations and binding in 4d dynamical triangulation”, Nucl. Phys. B (Proc. Suppl.), 47, 613–616, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

de Bakker, B.V., and Smit, J., “Gravitational binding in 4d dynamical triangulation”, Nucl. Phys. B, 484, 476–494, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

De Pietri, R., and Rovelli, C., “Geometry eigenvalues and the scalar product from recoupling theory in loop quantum gravity”, Phys. Rev. D, 54, 2664–2690, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Egawa, H.S., Hotta, T., Izubuchi, T., Tsuda, N., and Yukawa, T., “Scaling behaviour in 4d simplicial quantum gravity”, Prog. Theor. Phys., 97, 539–552, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Egawa, H.S., Hotta, T., Izubuchi, T., Tsuda, N., and Yukawa, T., “Scaling structures in fourdimensional simplicial gravity”, Nucl. Phys. B (Proc. Suppl.), 53, 760–762, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Egawa, H.S., Tsuda, N., and Yukawa, T., “Common structures in 2d, 3d and 4d simplicial quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 63, 736–738, (1998). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Egawa, H.S., Tsuda, N., and Yukawa, T., “Common structures in simplicial quantum gravity”, Phys. Lett. B, 459, 97–104, (1998). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Ezawa, K., “Multiplaquette solutions for discretized Ashtekar gravity”, Mod. Phys. Lett. A, 11, 2921–2932, (1996).

Ezawa, K., “Nonperturbative solutions for canonical quantum gravity: An overview”, Phys. Rep., 286, 271–348, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Feinberg, G., Friedberg, R., Lee, T.D., and Ren, H.C., “Lattice gravity near the continuum limit”, Nucl. Phys. B, 245, 343–368, (1984).

Fort, H., Gambini, R., and Griego, J., “Lattice knot theory and quantum gravity in the loop representation”, Phys. Rev. D, 56, 2127–2143, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Friedberg, R., and Lee, T.D., “Derivation of Regge's action from Einstein's theory of general relativity”, Nucl. Phys. B, 242, 145–166, (1984).

Fröhlich, J., “Regge calculus and discretized gravitational functional integrals”, in NonPerturbative Quantum Field Theory: Mathematical Aspects and Applications. Selected Papers of Jürg Fröhlich, vol. 15 of Advanced Series in Mathematical Physics, 523–545, (World Scientific, Singapore; River Edge, U.S.A., 1992).

Gross, M., and Varsted, S., “Elementary moves and ergodicity in ddimensional simplicial quantum gravity”, Nucl. Phys. B, 378, 367–380, (1992).

Hamber, H.W., “Simplicial quantum gravity”, in Osterwalder, K., and Stora, R., eds., Critical Phenomena, Random Systems, Gauge Theories, Proceedings of the Les Houches Summer School, Session XLIII, 1 August 7 September 1984), 375–439, (Elsevier, Amsterdam, Netherlands; New York, U.S.A., 1986).

Hamber, H.W., “Simplicial quantum gravity from twodimensions to fourdimensions”, in Damgaard, P.H., Hüffel, H., and Rosenblum, A., eds., Probabilistic Methods in Quantum Field Theory and Quantum Gravity, Proceedings of the NATO Advanced Research Workshop, held August 21–27, 1989, in Cargèse, France, vol. 224 of NATO ASI Series B, 77–87, (Plenum Press, New York, U.S.A., 1990).

Hamber, H.W., “Critical behavior in simplicial quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 20, 728–732, (1991).

Hamber, H.W., “Simulations of discrete quantized gravity”, Int. J. Supercomput. Appl., 5, 84–97, (1991).

Hamber, H.W., “Fluctuations and correlations in simplicial quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 26, 581–583, (1992).

Hamber, H.W., “Phases of 4d simplicial quantum gravity”, Phys. Rev. D, 45, 507–512, (1992).

Hamber, H.W., “Phases of simplicial quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 25, 150–175, (1992).

Hamber, H.W., “Phases of simplicial quantum gravity in four dimensions: estimates for the critical exponents”, Nucl. Phys. B, 400, 347–389, (1993).

Hamber, H.W., “Smooth and rough phases of quantized gravity”, Nucl. Phys. B (Proc. Suppl.), 30, 751–755, (1993).

Hamber, H.W., “Invariant correlations in simplicial gravity”, Phys. Rev. D, 50, 3932–3941, (1994). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Hamber, H.W., “Scalar fields coupled to fourdimensional lattice gravity”, in Carr, J., and Perrottet, M., eds., High Energy Physics (HEP 93), Proceedings of the International Europhysics Conference, Marseille, France, 22–28 July 1993, 211–213, (Editions Frontières, GifsurYvette, France, 1994). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Hamber, H.W., and Williams, R.M., “Higher derivative quantum gravity on a simplicial lattice”, Nucl. Phys. B, 248, 392–415, (1984).

Hamber, H.W., and Williams, R.M., “Nonperturbative simplicial quantum gravity”, Phys. Lett. B, 157, 368–374, (1985).

Hamber, H.W., and Williams, R.M., “Simplicial quantum gravity with higher derivative terms: formalism and numerical results in fourdimensions”, Nucl. Phys. B, 269, 712–734, (1986).

Hamber, H.W., and Williams, R.M., “Simplicial gravity coupled to scalar matter”, Nucl. Phys. B, 415, 463–496, (1994). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Hamber, H.W., and Williams, R.M., “Newtonian potential in quantum Regge gravity”, Nucl. Phys. B, 435, 361–398, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Hamber, H.W., and Williams, R.M., “Gauge invariance in simplicial gravity”, Nucl. Phys. B, 487, 345–408, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Hamber, H.W., and Williams, R.M., “On the measure in simplicial gravity”, Phys. Rev. D, 59, 064014–1–8, (1997).

Hartle, J.B., “Simplicial minisuperspace. I. General discussion”, J. Math. Phys., 26, 804–814, (1985).

Hartle, J.B., “Numerical quantum gravity”, in Sato, H., and Nakamura, T., eds., Gravitational Collapse and General Relativity, Proceedings of Yamada Conference XIV, Kyoto International Conference Hall, Japan, 7–11 April 1986, 329–338, (World Scientific, Singapore, 1986).

Hartle, J.B., “Simplicial minisuperspace. II. Some classical solutions on simple triangulations”, J. Math. Phys., 27, 287–295, (1986).

Hartle, J.B., “Simplicial minisuperspace. III. Integration contours in a fivesimplex model”, J. Math. Phys., 30, 452–460, (1989).

Hartle, J.B., and Sorkin, R.D., “Boundary terms in the action for the Regge calculus”, Gen. Relativ. Gravit., 13, 541–549, (1981).

Hotta, T., Izubuchi, T., and Nishimura, J., “Singular vertices in the strong coupling phase of fourdimensional simplicial gravity”, Prog. Theor. Phys., 94, 263–270, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Hotta, T., Izubuchi, T., and Nishimura, J., “Singular vertices in the strong coupling phase of fourdimensional simplicial gravity”, Nucl. Phys. B (Proc. Suppl.), 47, 609–612, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Immirzi, G., “Quantizing Regge Calculus”, Class. Quantum Grav., 13, 2385–2394, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Immirzi, G., “Quantum gravity and Regge calculus”, Nucl. Phys. B (Proc. Suppl.), 57, 65–72, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Isham, C.J., “Quantum gravity”, in MacCallum, M.A.H., ed., General Relativity and Gravitation, Proceedings of the 11th International Conference on General Relativity and Gravitation, Stockholm, July 6–12, 1986, 99–127, (Cambridge University Press, Cambridge, U.K., New York, U.S.A., 1987).

Jevicki, A., and Ninomiya, M., “Functional formulation of Regge gravity”, Phys. Rev. D, 33, 1634–1637, (1986).

Johnston, D.A., “Gravity and random surfaces on the lattice: a review”, Nucl. Phys. B (Proc. Suppl.), 53, 43–55, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Jurkiewicz, J., “Simplicial gravity and random surfaces”, Nucl. Phys. B (Proc. Suppl.), 30, 108–121, (1993).

Kaku, M., “Lattices, gravity and supergravity”, in Hawking, S.W., and Roček, M., eds., Superspace and Supergravity, Proceedings of the Nuffield workshop, Cambridge, June 16 July 12, 1980, 503–515, (Cambridge University Press, Cambridge, U.K.; New York, U.S.A., 1981).

Kaku, M., “Strongcoupling approach to the quantization of conformal gravity”, Phys. Rev. D, 27, 2819–2834, (1983).

Kaku, M., “Generally covariant lattices, the random calculus, and the strong coupling approach to the renormalization of gravity”, in Batalin, I.A., Isham, C.J., and Vilkovisky, G.A., eds., Quantum Field Theory and Quantum Statistics. Essays in Honour of the Sixtieth Birthday of E.S. Fradkin, Vol. 2: Models of Field Theory, 141–163, (A. Hilger, Bristol, U.K., 1987).

Kawamoto, N., and Nielsen, H.B., “Lattice gauge gravity with fermions”, Phys. Rev. D, 43, 1150–1156, (1991).

Khatsymovsky, V., “On the quantization of Regge links”, Phys. Lett. B, 323, 292–295, (1994).

Kondo, K., “Euclidean quantum gravity on a flat lattice”, Prog. Theor. Phys., 72, 841–852, (1984).

Kristjansen, C.F., “Higher derivative regularization in 4d quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 30, 756–759, (1993).

Krzywicki, A., “Perspectives in lattice gravity”, Acta Phys. Pol. B, 27, 827–838, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Leutwyler, H., “Gravitational field: equivalence of Feynman quantization and canonical quantization”, Phys. Rev., 134, 1155–1182, (1964).

Loll, R., “Nonperturbative solutions for lattice quantum gravity”, Nucl. Phys. B, 444, 619–639, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Loll, R., “The volume operator in discretized quantum gravity”, Phys. Rev. Lett., 75, 3048–3051, (1995). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Loll, R., “A real alternative to quantum gravity in loop space”, Phys. Rev. D, 54, 5381–5384, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Loll, R., “Spectrum of the volume operator in quantum gravity”, Nucl. Phys. B, 460(1), 143–154, (1996). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Loll, R., “Further results on geometric operators in quantum gravity”, Class. Quantum Grav., 14, 1725–1741, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Loll, R., “Imposing det E
$>$
0 in discrete quantum gravity”, Phys. Lett. B, 399, 227–232, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Loll, R., “Latticing quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 57, 255–258, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Loll, R., “Quantizing canonical gravity in the real domain”, in Dremin, I.M., and Semikhatov, A.M., eds., Second International A. D. Sakharov Conference on Physics, Proceedings of the conference held in Moscow, Russia 20–24 May 1996, 280–283, (World Scientific, Singapore; River Edge, U.S.A., 1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Loll, R., “Simplifying the spectral analysis of the volume operator”, Nucl. Phys. B, 500, 405–420, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Loll, R., “Still on the way to quantizing gravity”, in Bassan, M. et al., ed., Proceedings of the 12th Italian Conference on General Relativity and Gravitational Physics, Rome, Italy, September 23–27, 1996, 193–206, (World Scientific, Singapore; River Edge, U.S.A., 1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Loll, R., “On the diffeomorphism commutators of lattice quantum gravity”, Class. Quantum Grav., 15, 799–809, (1998). Related online version (cited on 19 November 1998):
.
☻ open access ✓

MacDowell, S.W., and Mansouri, R., “Unified geometric theory of gravity and supergravity”, Phys. Rev. Lett., 38, 739–742, (1977).

Mäkelä, J., “Phase space coordinates and the Hamiltonian constraint of Regge calculus”, Phys. Rev. D, 49, 2882–2896, (1994).

Mannion, C.L.T., and Taylor, J.G., “General relativity on a flat lattice”, Phys. Lett. B, 100, 261–266, (1981).

Menotti, P., “Nonperturbative quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 17, 29–38, (1990).

Menotti, P., “The role of diffeomorphisms in the integration over a finite dimensional space of geometries”, Nucl. Phys. B (Proc. Suppl.), 63, 760–762, (1998). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Menotti, P., and Peirano, P.P., “Diffeomorphism invariant measure for finite dimensional geometries”, Nucl. Phys. B, 488, 719–734, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Menotti, P., and Peirano, P.P., “Functional integration for Regge gravity”, Nucl. Phys. B (Proc. Suppl.), 57, 82–90, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Menotti, P., and Peirano, P.P., “Functional integration on Regge geometries”, Nucl. Phys. B (Proc. Suppl.), 53, 780–782, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Menotti, P., and Pelissetto, A., “Reflection positivity and graviton doubling in Euclidean lattice gravity”, Ann. Phys. (N.Y.), 170(2), 287–309, (1986).

Menotti, P., and Pelissetto, A., “Gauge invariance and functional integration measure in lattice gravity”, Nucl. Phys. B, 288, 813–831, (1987).

Menotti, P., and Pelissetto, A., “Poincaré, de Sitter, and conformal gravity on the lattice”, Phys. Rev. D, 35, 1194–1204, (1987).

Misner, C.W., Thorne, K.S., and Wheeler, J.A., Gravitation, (W.H. Freeman, San Francisco, U.S.A., 1973).

Myers, E., “Unbounded action and `density of states”', Class. Quantum Grav., 9, 405–411, (1992).

Nabutovsky, A., and BenAv, R., “Noncomputability arising in dynamical triangulation model of fourdimensional quantum gravity”, Commun. Math. Phys., 157, 93–98, (1993).

Pachner, U., “Konstruktionsmethoden und das kombinatorische Homöomorphieproblem für Triangulationen kompakter semilinearer Mannigfaltigkeiten”, Abh. Math. Sem. Univ. Hamburg, 57, 69–85, (1986).

Pachner, U., “P.L. homeomorphic manifolds are equivalent by elementary shellings”, Europ. J. Combinatorics, 12, 129–145, (1991).

Rebbi, C., ed., Lattice gauge theories and Monte Carlo simulations, (World Scientific, Singapore, 1983).

Regge, T., “General Relativity without coordinates”, Nuovo Cimento A, 19, 558–571, (1961).

Reisenberger, M.P., “A lattice worldsheet sum for 4d Euclidean general relativity”, (1997). URL (cited on 19 November 1998):
.
☻ open access ✓

Ren, H.C., “Matter fields in lattice gravity”, Nucl. Phys. B, 301, 661–684, (1988).

Renken, R.L., “The renormalization group and dynamical triangulations”, Nucl. Phys. B (Proc. Suppl.), 53, 783–785, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Renken, R.L., “A renormalization group for dynamical triangulations in arbitrary dimensions”, Nucl. Phys. B, 485, 503–516, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Renteln, P., “Some results of SU(2) spinorial lattice gravity”, Class. Quantum Grav., 7, 493–502, (1990).

Renteln, P., and Smolin, L., “A lattice approach to spinorial quantum gravity”, Class. Quantum Grav., 6, 275–294, (1989).

Roček, M., and Williams, R.M., “Quantum Regge calculus”, Phys. Lett. B, 104, 31–37, (1981).

Roček, M., and Williams, R.M., “Introduction to quantum Regge calculus”, in Isham, C., and Duff, M., eds., Quantum Structure of Space and Time, Proceedings of the Nuffield Workshop, Imperial College, London, 3–21 August, 1981, 105–116, (Cambridge University Press, Cambridge, U.K.; New York, U.S.A., 1982).

Roček, M., and Williams, R.M., “The quantization of Regge calculus”, Z. Phys. C, 21, 371–381, (1984).

Römer, H., and Zähringer, M., “Functional integration and the diffeomorphism group in Euclidean lattice quantum gravity”, Class. Quantum Grav., 3, 897–910, (1986).

Smit, J., “Remarks on the quantum gravity interpretation of 4D dynamical triangulation”, Nucl. Phys. B (Proc. Suppl.), 53, 786–788, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Smolin, L., “Quantum gravity on a lattice”, Nucl. Phys. B, 148, 333–372, (1979).

Sorkin, R.D., “Time evolution problem in Regge calculus”, Phys. Rev. D, 12, 385–396, (1975).

Thiemann, T., “Closed formula for the matrix elements of the volume operator in canonical quantum gravity”, J. Math. Phys., 39, 3347–3371, (1998). Related online version (cited on 19 November 1998):
.
☻ open access ✓

Tomboulis, E.T., “Unitarity in higherderivative quantum gravity”, Phys. Rev. Lett., 52, 1173–1176, (1984).

Tuckey, P.A., and Williams, R.M., “Regge calculus: a brief review and bibliography”, Class. Quantum Grav., 9, 1409–1422, (1992).

Varsted, S., “Threeand fourdimensional simplicial quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 26, 578–580, (1992).

Varsted, S., “Fourdimensional quantum gravity by dynamical triangulations”, Nucl. Phys. B, 412, 406–414, (1994).

Weingarten, D., “Euclidean quantum gravity on a lattice”, Nucl. Phys. B, 210, 229–245, (1982).

West, P.C., “A geometric gravity Lagrangian”, Phys. Lett. B, 76, 569–570, (1978).

Wheater, J.F., “Random surfaces and lattice quantum gravity”, Nucl. Phys. B (Proc. Suppl.), 34, 15–28, (1994).

Williams, R.M., “Quantum Regge calculus in the Lorentzian domain and its Hamiltonian formulation”, Class. Quantum Grav., 3, 853–869, (1986).

Williams, R.M., “Discrete quantum gravity: the Regge calculus approach”, Int. J. Mod. Phys. B, 6, 2097–2108, (1992).

Williams, R.M., “Simplicial quantum gravity in four dimensions and invariants of discretized manifolds”, J. Math. Phys., 36, 6276–6287, (1995).

Williams, R.M., “Recent progress in Regge calculus”, Nucl. Phys. B (Proc. Suppl.), 57, 73–81, (1997). Related online version (cited on 19 November 1998):
.
☻ open access ✓
Note: The reference version of this article is published by Living Reviews in Relativity