PseudoCartesian coordinates in a model of Causal Dynamical Triangulations
J. Ambjørn, Z. Drogosz, J. GizbertStudnicki, A. Görlich, J. Jurkiewicz
The Niels Bohr Institute, Copenhagen University
Blegdamsvej 17, DK2100 Copenhagen Ø, Denmark.
email: ambjorn@nbi.dk
IMAPP, Radboud University,
Nijmegen, PO Box 9010, The Netherlands
The M. Smoluchowski Institute of Physics, Jagiellonian University,
Lojasiewicza 11, PL 30348, Krakow, Poland.
email: zbigniew.drogosz@uj.edu.pl, jakub.gizbertstudnicki@uj.edu.pl, andrzej.goerlich@uj.edu.pl, jerzy.jurkiewicz@uj.edu.pl
Abstract
Causal Dynamical Triangulations is a nonperturbative quantum gravity model, defined with a lattice
cutoff. The model can be viewed as defined with a proper time but with no reference
to any threedimensional spatial background geometry. It has four phases, depending on the parameters
(the coupling constants) of the model. The particularly interesting behavior is observed in the socalled de Sitter phase, where the spatial threevolume
distribution as a function of proper time has a semiclassical
behavior which can be obtained from an effective minisuperspace action.
In the case of the threesphere spatial topology, it has been difficult to extend the
effective semiclassical description in terms of proper time and
spatial threevolume to include genuine spatial coordinates, partially because of the
background independence inherent in the model. However,
if the spatial topology is that of a threetorus, it is possible to define a number of new
observables that might serve as spatial coordinates as well as new observables related to
the winding numbers of the threedimensional torus.
The present paper outlines how to define the observables,
and how they can be used in numerical simulations of the model.
Keywords: Causal Dynamical Triangulations; Quantum Gravity
1 Introduction
The model of Causal Dynamical Triangulations (CDT) in four dimensions was originally formulated for systems where the spatial topology of the Universe was assumed to be spherical (). In most cases the (Euclidean) proper time was assumed to be periodic with some period . The spherical topology of the spatial part of the Universe gives a number of advantages, among which the main one is its relative simplicity. Unfortunately it also has disadvantages. In the original model the path integral includes a summation over all spatial simplicial geometries corresponding to abstract triangulations of . This background independence makes it almost impossible to define reference points, except in the time direction. Another problem is caused by the numerical algorithm used in the computer simulations, which made the most interesting range in the coupling constants, i.e., the one where the phase transition lines meet, practically inaccessible. These two deficiencies may be overcome using a different choice of the spatial topology, namely the toroidal topology (). In this case, the path integral defining the model will include a summation over abstract triangulations of .
CDT with a toroidal spatial topology was shown to have a phase structure similar to that observed when the
spatial topology was spherical [1]. We have a system with four distinct phases,
where the phase transition lines are approximately the same as in
the spherical case
A finite system with a toroidal geometry can be viewed as an infinite system, where a finite elementary cell is periodically repeated infinitely many times in three directions. The resulting copies of the elementary cell may be numbered by the winding numbers. The boundary of the elementary cell is not uniquely defined. However, as we will show below, it is possible to introduce additional constraints that can make such a definition unique.
In this article we will discuss the possibility of using the set of boundaries as a reference frame that allows us to define a set of pseudoCartesian coordinates in the elementary cell. The spatial volume distribution associated with these coordinates can be interpreted as a distribution of , which potentially permits measuring the effective action parametrized not only by a time variable but also by spatial variables. We will also define a set of new observables related to the topologically closed geodesic loops with nontrivial winding numbers.
2 Toroidal topology
Details of a formulation of the CDT model were discussed earlier in a number of articles (for a review see e.g. [2]). Let us mention here a few essential points. The basic idea is to consider the analogue of the Feynman path integral for the amplitude between two spatial states with initial and final threegeometries. It is assumed that trajectories satisfy a condition of causality, which means that the considered spacetimes are globally hyperbolic (permitting a global time foliation) and that the spatial topology is fixed on all leaves of the foliation. We use the Reggetype discretization of the spacetime [3]. The spacetime is constructed from elementary building blocks – fourdimensional simplices – with a unique length of spatial edges and a unique length of time edges. Each vertex in the system has a uniquely defined integer time parameter. As a consequence, we have two types of foursimplices: the simplices (with four vertices at time and one at ) and simplices (with three vertices at time and two at ). The simplices are glued together along threedimensional tetrahedral faces to form a simplicial manifold with additional regularity constraints imposed (all simplices and subsimplices with a particular set of vertex labels appear at most once). The existence of a global foliation means that each configuration can be analytically continued between the Lorentzian and Euclidean geometry. The Wick rotation can be interpreted as analytic continuation of the length of time links. In the Euclidean formulation the memory of the time orientation is preserved. In most cases studied the systems are assumed to be periodic in the (Euclidean) time, which means that we do not need to specify initial and final spatial geometric configurations.
For every configuration of the system we may define a dual lattice, where each simplex becomes a dual vertex and each face becomes a dual link, etc. On a dual lattice, from each vertex emerge exactly five links, corresponding to exactly five neighbors of a simplex on a direct lattice. For simplicity, we will assume that all links on a dual lattice have the same length. The basic observable to be used will be the length of a geodesic line between two simplices, i.e., the minimal number of steps (links on the dual lattice) necessary to connect the two simplices. The geodesic line is usually not unique, but the minimal distance between simplices is always well defined.
A spacetime trajectory appearing in the path integral is weighted by the exponential of (minus) the HilbertEinstein action, which in Euclidean formulation becomes real and can be interpreted as a probability
(1) 
For a system with a finite number of simplices the action takes a particularly simple form
(2) 
In (2) the global numbers , and denote correspondingly the number of vertices and the numbers of simplices of each of the two types in a triangulation . The dimensionless coupling constants are – related to the inverse of the gravitational constant, the cosmological constant and – a function of the ratio between the time and spatial length of edges. Typically in a configuration . Different triangulations correspond in general to different piecewise linear geometries, and the effective number of triangulations parametrized by the same set of global numbers and at (i.e., with the same action) grows in the leading order as
(3) 
As a consequence, the amplitude is defined only for and the limit is the largevolume limit where the average number of simplices becomes large and where we may try to define a continuum limit. One can define similar models in 1+1 dimensions and 2+1 dimensions. Only in 1+1 dimensions can the model be solved analytically [4]. In 2+1 dimensions there exist partial analytic results [5], but to achieve a reliable understanding of the phase structure one has to rely on Monte Carlo simulations. In 3+1 dimensions Monte Carlo simulations are the only tool available [6]. In all the simulations of a 3+1 system we use a set of 7 local moves, which preserve topology and foliation and which are ergodic in the class of triangulations we use (see [2] for details). The Monte Carlo evolution is performed at a given set of bare couplings and . The limit is studied by considering a growing sequence of systems with volumes , where the volume is forced to fluctuate around . The initial configuration with a prescribed spatial topology and the time extent fixes the topology of the studied system to be . The computer program works for arbitrary spatial topology and arbitrary , so the correct choice of the initial configuration with the desired spatial topology is a very important element of each simulation, although different choices with the same topology give equivalent results.
The choice of an initial toroidal spatial geometry was discussed in our earlier paper [7]. The initial configuration used to start a Monte Carlo simulation is not a minimal configuration (in the sense discussed in [7]), but it is a very convenient choice from the point of view of the present article. The configuration is based on a triangulation of a fourdimensional hypercube [8]. The topological hypercube is divided into 16 foursimplices, out of which 10 are of the type and 6 of the type. To build a configuration with a topology, one needs a number of hypercubes that is even and greater or equal to four. The last restriction comes from the manifold conditions: the two vertices of a link must be different. The smallest such configuration, assuming that the time period is , has simplices. If has the topology of one can use a much smaller starting configuration.
A system with spatial topology can be viewed as an infinite system, where a finite size elementary cell is periodically repeated in all spatial directions (for periodic b.c. this is also the case in the time direction). The starting configuration described above corresponds to a system composed of periodically repeated triangulated fourcubes with simplices each. The definition of an elementary cell is not unique. We may additionally request that the boundary between the neighboring cells has a minimal threevolume. This condition, for a starting configuration, does not lead to a unique choice: in fact, in every spatial direction we may choose the initial position of a boundary in four equivalent ways, which translates into equivalent spatial boundaries. This is a result of the relatively high symmetry of the initial configuration, which is invariant under time and spatial shifts by one unit. The exact shape of each boundary may be locally modified in such a way that it remains simply connected. This, in general, changes the threevolume of the boundary. To control the position of boundaries, we introduce in the coding of the geometry an additional information for each interface between neighboring simplices (i.e., for each dual link). In the fourdimensional case each simplex has five neighbors, as explained above. For each dual link we introduce the integervalued spatial vector , with entries having values in each spatial direction. Value 0 corresponds to the case where both the simplex and a particular neighbor are in the same elementary cell. Values mean that crossing the interface we enter the neighboring elementary cell in the positive or the negative direction. For the initial configuration choice we assign values to the vector field for all interfaces between neighboring simplices. In any given configuration, we may form a sum of vector fields along any loop joining the simplices on a dual lattice. For a trivial closed loop the sum must be a zero vector. For a nontrivial closed loop (closed by periodic b.c.) the sum gives topological information about the loop’s winding number.
In the numerical simulations the standard algorithm is using a finite set of local moves. Each move affects only a finite part of the configuration (i.e., of the triangulation), keeping the rest unchanged. Since HilbertEinstein action (2) does not depend on the position of the boundaries, the decision to accept (or reject) the move does not depend on a position of the boundaries. The position of the boundaries will however play an important role in the construction proposed in this paper. Therefore, before the move is performed we check if the boundaries between elementary cells pass through the affected region of the configuration. If this is not the case, the move is performed using the standard Monte Carlo algorithm. If, however, a boundary between two cells crosses the part of the triangulation where the geometry of connections between simplices is to be changed, then we modify locally the position of the boundary. The modification is done in such a way that the boundary remains connected, but is deformed so as to lay wholly outside of the region of the triangulation that is going to be affected by the move. Such a modification can always be done and does not change the action, but in general it will increase the size of the boundary. The move can now be performed in a standard way. After the move, we check if a simple local modification of the boundary can reduce its size (its threevolume). If this is the case, the modification of the position of the boundary is made. From time to time,we perform an additional check, independent of the moves, whether a local modification of any part of the boundaries can minimize their volumes. If so, then such a modification is always accepted.
All measurements described in this article were performed for the toroidal spatial topology in the de Sitter phase (the socalled C phase) at a special point in the parameter space and for systems with and a sequence of volumes . The same point in the parameter space was earlier analyzed in simulations of systems where the configurations had spherical spatial topology [6], as well as in systems with configurations of toroidal spatial topology [7].
3 PseudoCartesian coordinates
A given set of boundaries can be used to define what we will denote as the pseudoCartesian coordinates. The procedure is as follows:

We start from a boundary
^{5} , say, orthogonal to the direction , and we mark all simplices adjacent to the boundary in the direction we define to be positive as having the coordinate . 
We move (on the dual lattice) in the positive direction to the second layer of simplices at a unit distance from the first layer. We mark these simplices as having a coordinate .

We continue the same process until all simplices in the system are marked, and a maximal extent in the direction is reached.

The same procedure can be started from the layer of simplices adjacent to a boundary in the negative direction. The simplices in that layer are marked as having a coordinate .

We continue the process as before.

This way each simplex is assigned the values of the coordinates and .

Exactly the same method can be used to define the coordinates , , and .
For a (twodimensional) visualization of the procedure please refer to Figure 1.
One should note that although each threedimensional boundary is simply connected, this does not need to be the case for threedimensional surfaces separating and layers. The set of numbers represents the distance of a simplex to a boundary, i.e., the number of links on the dual lattice, in negative and positive directions.
The first question one may ask is if the definition of a boundary between elementary cells, as implemented by the minimization procedure, is unique. For the simple initial configuration described above this is definitely not the case. As already explained, in each spatial direction we may choose the initial position of a boundary between the copies of the elementary cell in four equivalent ways, giving rise to equivalent (minimal) boundaries. Since the updates of geometry do not “feel” the position of the boundaries, we can run a parallel simulation, where the evolution of geometry is exactly the same, but the initial choice of a position of a boundary is different. In this way, during the thermalization process, each simplex can be assigned its coordinates using the two different positions of a boundary. Comparing the two coordinates, say and , we define the difference between the two alternative definitions by measuring the ratio of the number of simplices for which to the total number of simplices. This quantity is initially zero. A value means that two definitions become equivalent, or in other words, that the two boundaries fully overlap. Analogous quantities can be defined in and directions. In the figure 2 we present the evolution of the parameter for three directions in a system with 160k. The initial position of boundaries was chosen as the most distant permitted by the symmetry of the starting configuration. Similar experiments were repeated for other initial choices of boundaries and other system sizes. In all cases the result was the same, indicating that the algorithm used is powerful enough to find a global minimum of the boundary size in all directions. It is a nontrivial result, because naïvely one may expect the existence of a complicated landscape of local minima. This also means that we do not have to worry about the initial choice of the position of boundaries.
Since the considered systems are periodic in (Euclidean) time with a period , each simplex can be assigned a coordinate . The layer with is chosen as the set of simplices with four vertices at . We follow the same prescription as before, for instance is composed of a set of simplices at a unit distance from the in the positive time direction. We continue until we reach the layer formed from simplices. We could define coordinates similarly as before, but for the time direction for each system we have a trivial relation , resulting from the trivial periodicity by construction in the time direction.
In the spatial directions a similar relation is not satisfied and the distribution of values for, say, and is nontrivial. Below, in the figure 3, we present the probability distribution to find a simplex with given values of coordinates and in a system with 160k. The distribution is constructed by averaging over 800 statistically independent configurations and summed over all simplices in the system.
Distributions of this type contain interesting information about geometric properties of the system in spatial directions. We should remember that simplices with coordinates and lie on a minimal boundary between elementary cells. Such a boundary separates regions where gravitational fluctuations produce volume concentrations. Qualitatively, in the small and regions the volume of the system is suppressed. A special role is played by the distribution of the quantity (and similar quantities and ). Although for the initial configuration the shape of the elementary cell is symmetric with respect to the exchange of directions , and , during the thermalization process the shape of the elementary cell may be deformed. In fact, we expect that the distributions , and should overlap, up to a finite shift in . This is indeed the case. In the figure 4 we show the distributions , and for systems with and 160k. In both cases a complete overlap is achieved by applying shifts of order one.
Comparing the two plots we realize that the shape of distributions seems to be the same, up to a scaling, depending on the total volume. In the figure 5 we compare the two volumes (80k and 160k), applying a scaling factor to the distribution for a larger system. This agrees with the expected scaling, if we assume that the Hausdorff dimension in de Sitter phase is .
Another important information about the shape of the volume distribution in spatial directions is hidden in the distribution of the difference between and . For a toroidal topology one expects this distribution to be approximately constant, as was already observed in [7] for a similar distribution in time . While the periodicity in time was exact, the period in the spatial directions, say the direction, will vary, depending on . As a consequence, we break the spatial translational symmetry and force the system to have the smallest volume at the small limit. In the central range of the plot, where we expect the distributions to be approximately flat, constant in the infinite volume limit. Below, in the figure 6, we show the amplified dependence of the volume distribution as a function of the rescaled variable for a range of values of for a system with 160k.
We see that the distributions are indeed approximately flat in the central range. The large volume limit can be extrapolated by comparing distributions of for systems with volumes 80k and 160k. We see that for a larger system the distribution gets wider (Figure 7).
4 The fractal structure of equal layers
As explained above, the constant layer is by construction connected for . We will study the geometric structure of layers with higher values of . To do this we define the concept of a slice. The idea is very similar to that used in 2d Euclidean Dynamical Triangulations. The construction is based on following the front of a diffusion wave on the dual lattice. We start at a randomly chosen simplex with spatial coordinate and constrain the diffusion process to take place in the subset of simplices with spatial coordinate , i.e., in the bounded part of the given triangulation lying in between the set of simplices with coordinate and the boundary of the elementary cell. In each diffusion step we mark simplices we meet on the way. The marked simplices form a new limiting layer for the next step.
If the geometry of the manifold were trivial, in such a process we could reach all simplices in the bounded region of the configuration described above. The fractal nature of geometry means that, in general, only a part of the bounded region is reached. Completing the process we find a set of simplices with the same coordinate as the initial simplex. The set of these simplices forms a slice. All simplices belonging to the slice can be joined by at least one path lying in the bounded region. Repeating the same process starting at a different simplex with a coordinate , we either end with the same slice or produce a different slice, disconnected from the previous one. The set of simplices reached from the set belonging to a common slice has a treelike structure, where the slice plays a role of the root.
The procedure is performed repeatedly, assigning in this way a slice index to all simplices of the configuration:

The first slice consists of all simplices with .

In the layer with we find all simplices that can be connected by a path that does not cross the elementary cell boundary and never goes below .

We find all separate slices for .

We repeat the same procedure for . In this case, we restrict the class of paths not to go below .

We continue until all slices are found.
For a (twodimensional) visualization of the procedure please refer to Figure 1.
By construction, each slice (except the first one) has exactly one parent slice, and it may have a larger number of children (or no children at all). Following the sequence of connected slices with an increasing value of , we may either reach the boundary or end in a blind alley. In the first case we say that the sequence belongs to the trunk, in the second case it belongs to a branch. The structure looks like a tree, with the trunk connecting the lower and upper boundaries. The trunk may split into several outgrowths. The branches emerge from the trunk slices and, by definition, never reach the boundary of the elementary cell. We would like to interpret the branches as the effect of quantum fluctuations and the trunk as a semiclassical background.
Below we illustrate the properties of the trunkbranch construction for one particular wellthermalized configuration in the de Sitter phase. The configuration was obtained for a system with , , and the total number of simplices . At each value of the coordinate we split the volume into the part belonging to the trunk and the part belonging to branches. In the figure 8 we show the fraction (percent) of volume in the trunk, relative to the total volume at a coordinate , and the fraction of volume of the largest slice, relative to the trunk volume at the same coordinate . In the same plot (in red) we show the total volume distribution as a function of . The distribution is normalized to reach at its maximum. The distribution in the plot (as compared to that in the figure 3) can be viewed as a projection of volume on the axis. In the plot all distributions are averaged over directions , , , , and .
As can be seen, for a configuration in the de Sitter phase, the trunk contains the biggest part of volume and the largest slice captures almost all trunk volume up to the value of where the volume is maximal (recall that by construction the volume as a function of is small near the boundaries).
In the next plot, figure 9, we show the average number of branches at a distance and the average volume per branch. As can be seen, the number of branch slices grows together with volume, but the average volume per branch slice is relatively small. The maximal length of a branch, defined as a number of steps between the consecutive slices along the branch necessary to reach the trunk slice from a branch slice, is six, which means that branches are rather short.
The tree structure representing the branching of slices with the increasing value of a coordinate is presented in the figure 10. For clarity the plot contains only the slices belonging to the trunk. The volume of slices is shown as the horizontal length of the line.
5 Loops with a nontrivial winding number
In the former section we analyzed the fractal structure of an elementary cell, parametrized by what we denoted as
pseudoCartesian coordinates. Using the same elementary cell and its periodic extensions to neighboring cells,
we will now introduce a set of new topological observables, which measure the geodesic distance between a simplex
and its copy in the neighboring cells in different spatial directions.
Like before, the measurement is performed by following paths defined by the front of
a diffusion wave, starting from a given simplex, but
in this case we consider the torus as an infinite system where simplices
reached after diffusion steps are labeled by the (unique) index in the elementary cell and also by the number of times a path crossed the
boundaries. If one finds a simplex with the same cell index and nontrivial number of crossings, say,
, it means that there exists at least one topologically closed geodesic line characterized by these winding numbers, namely the shortest such path
We expect the distributions of lengths of all loops connecting the neighboring cells in , and directions to be approximately the same, up to a possible small shift in the length . This is indeed what is observed. The figure 11 shows the length distributions of the , and loops superimposed using a shift of order smaller than 1. The match is not perfect, but we should remember that we are using just a single configuration. The need of making small finite shifts reflects the fact that during the thermalization process the shape of the elementary cell is deformed from the regular symmetric shape it had in the initial hypercube.
Comparing the loop distance in two different directions (Figure 11), we find that not only is the shape of the distributions similar. In fact loop distances for all simplices are highly correlated, which can be seen in the figure 12, which presents the correlation between the loop distance in the and directions.
For all simplices in the configuration we assigned a set of numbers measuring distances to their copies in the neighboring cells. These numbers contain interesting information about the internal geometry of the manifold.

In each direction we find a finite number of loops with the minimal length . The minimal loops follow the deepest valleys being surrounded by mountains which we can view as large volume fluctuations. In some sense, the minimal loops are dual to the minimal boundaries discussed in previous sections.

The simplices located on the mountains (i.e., inside the large volume fluctuations) are characterized by longer loop distances to the neighboring cells. Given a simplex, we can compare its loop distance in each direction with the loop distances of its five neighbors. For all simplices with the same loop distance , we find that the neighboring simplices are characterized by loop distances in the range , and .

Only for the few shortest loops can we follow a simple path connecting the simplex to its copy that would go through a sequence of neighboring simplices characterized by the same value of the loop distance .

For longer loops all paths to a copy pass through fragments of loops with smaller loop distances.

We deduce that the difference between the loop distance for a particular simplex and the minimal distance indicates how far a given simplex is from the minimal loop.
Finally, we will compare the distribution of loop distances for simplices in the fundamental cell in, say, the direction to the distribution of the quantity discussed earlier. We recall that this quantity is the sum of the distances from a simplex to the two cell boundaries orthogonal to the direction. We find that the two distributions are highly correlated. An example of this correlation can be seen on the contour plot shown in the figure 13.
6 Loop distances in other directions
The universality of the loop distance distribution can be extended to include loops with higher nontrivial winding numbers. The interpretation is that the bloblike structure (valleys and mountains) reflects the fractal structure of the fluctuations of geometry in the model and therefore can be expected to be independent of the direction in which we measure the loop distance distribution. For different directions we expect different minimal loop lengths, and as a consequence of comparing distributions with different argument range we may have to adjust this range to obtain a proper overlap when we compare the distributions. In Figure 14 we show the comparison of distributions for loop distances in directions (blue), (red), (green) and (orange). The last one is seen only as a fragment of a distribution, as it requires a shift of 15 units, and is not contained fully in the region we analyzed. The universality of the distributions, independent of directions, strongly suggests that the distribution of volume fluctuations is spherically symmetric. At the same time, the difference in the minimal loop distance in different directions most likely reflects that the shape of the elementary cell of our thermalized configuration is nontrivial.
The thermalized configuration discussed here has the fourdimensional topology of . In the time direction there is a foliation, which is absent in spatial directions. Although the periodicity in the time direction is trivial, this does not necessarily mean that the distribution of the loop distances in the time direction will be a Dirac delta function. In fact, due to the large spatial volume fluctuations, it has a shape very similar to that in spatial directions. Below, in the figure 15, we show the distribution of the loop distances in the time direction () (blue) as compared to a shifted and rescaled distribution in the direction (orange). The scaling factor used was 1.19.
This result indicates that although we have a foliation in time, the information encoded in the loop distance in spatial and time directions is very similar. The average length of the loop distance in the time direction is expected to depend on the time period .
7 Conclusion
CDT is a nonperturbative quantum gravity model. It is a model with a buildin cutoff,
the length of the links in the triangulations. In order to prove that it is an interesting model of
quantum gravity one has to be able to remove the cutoff in a controlled way and show that the resulting
theory satisfies the criteria suitable for a theory of quantum gravity
The use of toroidal topology also allowed us to use loops with nontrivial winding numbers to extract information about the fractal structure of our spatial geometries. This has been done successfully in twodimensional Euclidean quantum gravity [11], and also here in the fourdimensional case it provides us with a lot of information. In some sense, the shortest noncontractible loops can be considered as dual to the minimal cell boundaries and, like in the twodimensional case, the shortest loops lie in narrow “valleys”, where the surrounding “hills” look like (topologically) spherical outgrowths. Again, further computer simulations should allow one to determine the fractal structure related to “valleys” and “hills”, much in the spirit of twodimensional Euclidean quantum gravity.
Acknowledgments
JA wishes to acknowledge support from the Danish Research Council grant “Quantum Geometry”. JGS wishes to acknowledge support of the grant UMO2016/23/D/ST2/00289 from the National Science Centre, Poland. ZD acknowledges a support of the PhD project “Kartezjusz”. AG acknowledges support by the National Science Centre, Poland under grant 2015/17/D/ST2/03479.
Footnotes
 Some differences may be attributed to relatively large finite size effects.
 The detailed analysis of this range will be the subject of an independent publication.
 It should be emphasized that the existence of a quadruple point was only a conjecture in the spherical case since it was located in the region of coupling constant space that was inaccessible to computer simulations.
 For a piecewise linear fourdimensional geometry, defined by a fourdimensional triangulation, Regge [3] provided a simple geometric expression for the HilbertEinstein action, expressed in terms of link lengths and deficit angles of triangles in the triangulation. In our case, this expression simplifies further since our link lengths take only two values, corresponding to spacelike and timelike links. The end result is that the action can be expressed entirely in terms of the number of vertices , the number of foursimplices and , and the ratio between the length of spacelike and timelike links. See [2] for details.
 The boundaries are labeled and , based on their positions in the initial configuration, which is related to the hypercube.
 The “shortest” is meant in a graphtheoretical sense: we follow a link path (on the dual lattice) and the length is the number of links. This is not necessarily the shortest path if we consider our triangulation as a piecewise linear manifold. However, we do not expect this difference to be important when we discuss generic fractal properties, like the Hausdorff dimension of a generic triangulation in the limit of infinitely large triangulations. By an abuse of language we still call the shortest link path a geodesic.
 A discussion of how the CDT cutoff can in principle be removed while physics is kept constant can be found in [9].
References
 J. Ambjorn, J. GizbertStudnicki, A. Görlich, J. Jurkiewicz, D. Németh, The phase structure of Causal Dynamical Triangulations with toroidal spatial topology, JHEP 1806 (2018) 111, arXiv:1802.10434[hpth];
 J. Ambjorn, A. Görlich, J. Jurkiewicz, R.Loll, Nonperturbative Quantum Gravity, Phys.Rept. 519 (2012) 127210, arXiv:1203.3591[hepth];
 T. Regge, General relativity without coordinates, Nuovo Cim. 19 (1961) 558.

J. Ambjorn and R. Loll,
Nonperturbative Lorentzian quantum gravity, causality and topology change,
Nucl. Phys. B 536 (1998) 407, arXiv:
[hepth/9805108].
J. Ambjorn, R. Loll, W. Westra and S. Zohren, Putting a cap on causality violations in CDT, JHEP 0712 (2007) 017, arXiv: [arXiv:0709.2784 [grqc]].
J. Ambjorn, R. Loll, Y. Watabiki, W. Westra and S. Zohren, A String Field Theory based on Causal Dynamical Triangulations, JHEP 0805 (2008) 032, arXiv: [arXiv:0802.0719 [hepth]]. 
J. Ambjorn, J. Jurkiewicz, R. Loll and G. Vernizzi,
Lorentzian 3D gravity with wormholes via matrix models,
JHEP 0109 (2001) 022, arXiv:
[hepth/0106082];
3D Lorentzian quantum gravity from the asymmetric ABAB matrix model,,
Acta Phys. Polon. B 34 (2003) 4667, arXiv:
[hepth/0311072].
J. Ambjorn, J. Jurkiewicz and R. Loll, Renormalization of 3d quantum gravity from matrix models, Phys. Lett. B 581 (2004) 255, arXiv: [hepth/0307263]. 
J. Ambjorn, J. Jurkiewicz and R. Loll,
Emergence of a 4D world from causal quantum gravity,
Phys. Rev. Lett. 93 (2004) 131301, arXiv:
[hepth/0404156];
Reconstructing the universe,
Phys. Rev. D 72 (2005) 064014, arXiv:
[hepth/0505154].
J. Ambjorn, A. Gorlich, J. Jurkiewicz and R. Loll, Planckian Birth of the Quantum de Sitter Universe, Phys. Rev. Lett. 100 (2008) 091304, arXiv: [arXiv:0712.2485 [hepth]]; Phys. Rev. D 78 (2008) 063544, arXiv: [arXiv:0807.4481 [hepth]]. 
J. Ambjorn, Z. Drogosz, J. GizbertStudnicki, A. Görlich, J. Jurkiewicz, D. Németh,
Impact of topology in causal dynamical triangulations quantum gravity,
Phys.Rev. D94 (2016) 044010, arXiv:1604.08786[hepth];
J. Ambjorn, J. GizbertStudnicki, A. Görlich, K. Grosvenor, J. Jurkiewicz, Fourdimensional CDT with toroidal topology, Nucl.Phys. B922 (2017) 226246, arXiv:1705.07653[hepth];  P.S. Mara. Triangulations for the Cube, Journal of Combinatorial Theory (A)20:170177, 1976;
 J. Ambjorn, A. GÃ¶rlich, J. Jurkiewicz, A. Kreienbuehl and R. Loll, Renormalization Group Flow in CDT, Class. Quant. Grav. 31 (2014) 165003, arXiv:1405.4585 [hepth].
 H. Kawai, N. Kawamoto, T. Mogami and Y. Watabiki, Transfer matrix formalism for twodimensional quantum gravity and fractal structures of spacetime, Phys. Lett. B 306 (1993) 19, arXiv: [hepth/9302133].

J. AmbjÃ¸rn and T. Budd,
The toroidal Hausdorff dimension of 2d Euclidean quantum gravity,
Phys. Lett. B 724 (2013) 328, arXiv:
[arXiv:1305.3674 [hepth]];
Semiclassical Dynamical Triangulations,,
Phys. Lett. B 718 (2012) 200, arXiv: 1209.6031 [hepth]].
J. Ambjorn, J. Barkley and T. G. Budd, Roaming moduli space using dynamical triangulations, Nucl. Phys. B 858 (2012) 267, arXiv: 1110.4649 [hepth]].
J. Ambjorn, J. Barkley, T. Budd and R. Loll, Baby Universes Revisited, Phys. Lett. B 706 (2011) 86, arXiv: 1110.3998 [hepth].