arxiv: v1 [gr-qc] 19 Mar 2015

Transcription

1 Curvatures and discrete Gauss-Codazzi equation in +1-dimensional loop quantum gravity Seramika Ariwahoedi 1,, Jusak Sali Kosasih, Carlo Rovelli 1,, Freddy P. Zen 1 Aix Marseille Université, CNRS, CPT, UMR 7, 188 Marseille, France. Université de Toulon, CNRS, CPT, UMR 7, 8957 La Garde, France. Institut Teknologi Bandung, Bandung 401, West Java, Indonesia. We derive the Gauss-Codazzi equation in the holonomy and plane-angle representations and we use the result to write a Gauss-Codazzi equation for a discrete +1-dimensional manifold, triangulated by isosceles tetrahedra. This allows us to write operators acting on spin network states in +1- dimensional loop quantum gravity, representing the -dimensional intrinsic, -dimensional intrinsic, and -dimensional extrinsic curvatures. arxiv: v1 [gr-qc] 19 Mar 015 I. INTRODUCTION Three distinct notions of curvature are used in general relativity: the intrinsic curvature of the spacetime manifold, M, the intrinsic curvature of the spacial hypersurface Σ embedded in M which is utilised in the canonical framework, and the extrinsic curvature of Σ [1]. The three are related by the Gauss-Codazzi equation. On a discrete geometry, the definition of extrinsic curvature is not entirely clear []. Even less so in loop quantum gravity LQG, where the phase-space variables are derived from the first order formalism [], which generically does not yield a direct interpretation as spacetime geometry. But a definition of these curvatures is important in LQG, because we expect a discrete twisted geometry to emerge from the theory in an appropriate semi-classical limit [4]. In n-dimensional discrete geometry, the manifold is formed by n-simplices [, 5]. The intrinsics curvature sits on the n -simplices, called hinges, and can be defined by the angle of rotation: a vector parallel-transported around the hinge gets rotated by this angle. The rotation is in the hyperplane dual to the hinge. In the canonical formulation of general relativity it is convenient to use the ADM formalism [6] or a generalisation [7]. The phase space variables are defined on an n 1 surface: the initial time, or, more generally, boundary surface Σ. At the core of the ADM formalism is the Gauss-Codazzi equation, relating the intrinsic curvature of the 4D spacetime with the intrinsic and extrinsic curvatures of Σ. In the context of a discrete geometry, Σ is defined as an n 1-dimensional simplicial manifold formed by n 1-simplices. To develop a formalism analogous to the ADM one, we need an equation relating the intrinsic and extrinsic curvatures of Σ and M, a discrete version of the Gauss-Codazzi equation. In this paper, we explore a definition of extrinsic curvature for the discrete geometry which allows us to write an equation relating the extrinsic and intrinsics curvatures, which we refer as the Gauss-Codazzi equation for discrete geometry. To get some insight into this problem, we first review the standard Gauss-Codazzi equation, both in second and first order formalism. In Section II we derive the holonomy matrix representation and plane-angle representation of the Gauss-Codazzi equation. We consider the discrete version of the Gauss-Codazzi equation in Section III, where we define all the curvatures. In Section IV, we move to the loop quantum gravity picture and study the operators corresponding to these curvatures into operators. II. GAUSS-CODAZZI EQUATION A. Standard Gauss-Codazzi equation The Gauss-Codazzi equation relates the Riemann intrinsic curvatures of a manifold and its submanifold with the extrinsic curvature of the submanifold. We will briefly review the continuous +1-dimensional Gauss-Codazzi equation in this section, but the formula will be valid in n + 1-dimension.

2 1. Second order formulation Consider a -dimensional manifold M and a -dimensional surface Σ embedded in M. In the second order formulation of general relativity [8], the Riemann intrinsic curvature of M is a tensor: 1 R = R α µνβ α dx β dx µ dx ν, α, β, µ, ν = 0, 1, which can be thought as a map that rotates a vector v T p M parallel-transported along an infinitesimal square loop defined by unit vectors { µ, ν. Next, we define the proection tensor as: q µν = g µν n, n n µ n ν, with n µ is the normal to hypersurface Σ M. n, n = ±1, depends on the signature of the metric. The proected Riemann curvature of R on Σ is defined as: R Σ = Rα µνβ α dx β dx µ dx ν = qα α qµ µ qν ν q β β Rµ α ν β α dx β dx µ dx ν. Taking the proected part R Σ from the full part R, we have relation as follow: R = R Σ + S, decomposition formula, 1 S is the residual part of R. The proected part R Σ can be written as: R Σ = R + [K, K] Gauss equation, which can be written in terms of components: Rα µνβ = Rµνβ α + Kµ α K βν Kν α K βµ, Kαβ = q αq µ β ν µ n ν. R and K are, respectively, the -dimensional Riemannian curvature and extrinsic curvature of hypersurface Σ. The residual part S can be written as: Sµνβ α = n α ν K µβ µ K νβ Codazzi equation, with is the covariant derivative on the slice Σ. Writing 1 and together, we obtain: R = R + [K, K] +S R Σ.. First order formulation The gravitational field is a gauge field and can be written in a form closer to Yang-Mills theory. This way of representing gravity is known as first order formulation of general relativity [, 8]. Let manifold M be a spacetime. Let e be a local trivialization, a diffeomorphism map between the trivial vector bundle M R with the tangent bundle over M: T M = p {p T p M, and A be the connection on M R. The -dimensional intrinsic curvature of the connection is the curvature -form: F = F I µνjξ I ξ J dx µ dx ν, µ, ν, I, J = 0, 1,. which comes from the exterior covariant derivative of the connection: F = d D A. 4 { µ, dx µ and { ξ I, ξ I are local coordinate basis on M and R, respectively. In terms of components: F I µνj = µ A I νj ν A I µj + A I µka K νj A I νka K µj. 5

3 We use nice coordinates such that the time coordinate ξ 0 in fibre R is mapped by e to the time coordinate 0 in the basespace M, i.e., the map e is fixed into: e ξ 0 = e 0 = δ 0 µdx µ = dx 0. For the foliation, we use the time gauge, where the normal of the hypersurface Σ is taken to be the time direction: n = 0, e 1 n = ξ 0. Then the + 1 split is simply carried by spliting the indices as I = 0, a, and µ = 0, i, with 0 as the temporal part: F = F a0i0 + F 00ia + F ai00 F 0i0a ξ a ξ 0 dx i dx 0 + F ai0 F 0ia ξ a ξ 0 dx i dx F ai0b F a0ib ξ a ξ b dx i dx 0 + F aib ξ a ξ b dx i dx The proected part of the curvature -form on Σ is: and the residual part is: F Σ = F aib ξ a ξ b dx i dx, S = F a0i0 + F 00ia + F ai00 F 0i0a ξ a ξ 0 dx i dx 0 + F ai0 F 0ia ξ a ξ 0 dx i dx + F ai0b F a0ib ξ a ξ b dx i dx 0. Therefore, the decomposition formula is clearly: Next, we take only the proected part: which can be written in terms of components as follow: F = F Σ + S. 7 F Σ = F aib ξ a ξ b dx i dx, F aib = F aib + A a i0a 0 b A a 0A 0 ib Gauss equation. 8 The closed part of F Σ is clearly the D intrinsic curvature of connection in Σ and the rest is the extrinsic curvature part: F = F aib ξ a ξ b dx i dx = i A ab A aib + A aic A c b A ac A c ib. ξ a ξ b dx i dx, K = K a i ξ α dx i = A a i0ξ α ξ 0 dx i. The residual part S, in this special coordinates and special gauge fixing satisfies: F 0ia = D K ia D i K a Codazzi equation, and the other temporal components are zero. It must be kept in mind that this is the Gauss-Codazzi equation for a special local coordinate and special choice of gauge fixing, for a general case, they are not this simple. To conclude, we have the decomposition and Gauss-Codazzi equation for a fibre bundle of gravity: F = F + [ K, K ] +S. F Σ B. Gauss-Codazzi equation in holonomy representation 1. Holonomy around a loop In this section, we will write the Gauss-Codazzi equation in terms of holonomy. The holonomy is defined by the parallel transport of any section of a bundle, say, s 0, so that it satisfies the equation as follow: D v s 0 = s 0 + A s 0 = 0. 9

4 4 Solving 9 using recursive method [, 7, 8], we obtain the solution: s t = U γ t, D s 0, U γ t, D is holonomy of connection D along path γ t: U γ t, D = U γ P exp A, 10 with P is the path-ordered operator See [, 7] for the details of the derivation. Consider a square loop γ embedded in Σ which encloses a -dimensional area. The holonomy around the square loop can be written as a product of four holonomies, since holonomy is piecewise-linear: U γ = U γ4 U γ U γ U γ1, Taylor expanding the holonomy in 10 up to the second order [, 8, 9], we obtain: U γ = 1 + aµν F µν + O 4, with a µν is an infinitesimal square area inside loop γ. The formula: will be used to write the Gauss-Codazzi equation in terms of holonomies. a µν F µν = U γ 1 O 4 11 Contracting 6 by an infinitesimal area a µν, we obtain: a µν. First order formulation in holonomy representation F IµνJ = U IJ δ IJ O 4 = ai0 F a0i0 + F 00ia + F ai00 F 0i0a ξ a ξ 0 + ai F ai0 F 0ia ξ a ξ ai0 F ai0b F a0ib ξ a ξ b + ai F aibξ a ξ b, where we have used the result in 11, relating the holonomy with the curvature of the connection. Taking the proected part 8 contracted with a i, and using the relation between holonomy with the curvature -form, we can write the proected part in terms of -dimensional holonomy and the contracted extrinsic curvature as follow: a i F aib = ai F aib + ai K aik b K a K ib = U ab δ ab O 4 + ai K aik b K a K ib. Finally, collecting all these result together, we obtain: or simply: a µν F IµνJ ξ I ξ J {{ U IJ δ IJ O 4 ξ I ξ J F aib ξ a ξ b + ai {{ K aik b K a K ib ξ a ξ b U ab δ ab O 4 ξ a ξ b [K,K] ab ξ a ξ b Fa0i0 + F 00ia + F ai00 F 0i0a ξ a ξ 0 + F ai0b F a0ib ξ a ξ b = ai + ai0 + ai F ai0 F 0ia ξ a ξ 0 U = U + [K, K] + S + I I, 1 T where we have written the Gauss-Codazzi and the decomposition formula together, noting the residual terms as S. Remember that U is the holonomy of the connection, which describes the curvature of the connection of fibre F, not the curvature of the basespace M. S,

5 5. Second order formulation in holonomy representation In the same manner as above, for a tangent bundle T M, we obtain the relation a µν or simply: R αµνβ dx α dx β = + ai R kil dx k dx l {{ {{ H αβ δ αβ O 4 dx α dx β H kl δ kl O 4 dx k dx l + + ai K kik l K k K il dx k dx l [K,K] kl dx k dx l a i0 R 0i0 + R 00i + R i00 R 0i0 + ai R ik0 R 0ik + ai0 R i0k R 0ik dx dx k dx dx 0 S, H = H + [K, K] + S + I I, 1 T with H is the holonomy around any loop γ embedded in Σ M, coming from the Riemann tensor R. C. Gauss-Codazzi equation in plane-angle representation Rotations can be represented in two ways, using the holonomy matrix representation, or using the plane-angle representation. For a continuous theory, the holonomy representation provides a simpler way to do calculation concerning rotation. In a discrete theory, where we would like to get rid of coordinates, the plane-angle representation is a natural way to represent rotations and curvatures. The plane-angle representation is only a representation of the rotation group using its Lie algebra for the plane of rotation and one real parameter group times the norm of the algebra for the angle of rotation. There exist a biective map sending the holonomy to the plane-angle representation, known as the exponential map [10]. In this section, we rewrite the Gauss-Codazzi equation using the plane-angle representation. We only do the calculation for the second order formalism, the first order formalism version can be obtained in a similar way. Firstly, let us define the variables; we have two equations: the decomposition formula and the Gauss-Codazzi equation written compactly in 1. Since we are working in +1-dimension, it is natural to use matrix group SO to represents the -dimensional intrinsic curvature, and the subgroup SO to represents the -dimensional intrinsic curvature. But for simplicity of the calculation, we use the unitary group SU the double cover of SO instead of SO. Having information of H SU is equivalent with having information of the plane and angle of the - dimensional rotation, which are denoted, respectively, by {J su, δθ R. The same way goes with H SO, it contains same amount of information with the plane-angle of the -dimensional rotation { so, δφ R. Using the exponential map, we can write: H = exp δφ = I cos δφ + sin δφ, 14 with = [ ] 0 1, is the generator of SO. Doing the same way to the element of SU, we obtain: H exp Jδθ = I cos δθ + J sin δθ, 16 with J su is an element of the Lie algebra of SU, satisfying: σ a are the basis of su, namely, the generator of SU: satisfying the algebra structure relation as follow: J = J a σ a, J a J a = 1. σ a = iσ a, [ σ a, σ b ] = ε c ab σ c, 17

6 6 FIG. 1. Suppose we have a -dimensional curved manifold M discretized by four tetrahedron in the figure above in flat case it is known as 1-4 Pachner move. Then we take an embedded slice Σ as the surface of one tetrahedron the dark blue surface discretized by three triangle. Embedded on Σ, we take a loop γ circling a point of a tetrahedron. Attached to γ, are the SU and SO holonomy, which are related to the D and D intrinsic curvature, respectively. with σ a are the Pauli matrices: σ x = [ ] [ ] [ ] i 1 0, σ 1 0 y =, σ i 0 z = 0 1 The 1 factor in 16 comes out from the normalization occuring when we write SO using SU representation, i.e., to reduce the factor in relation 17 Taking the trace of 1 gives the relation between -dimensional and -dimensional rotation angle: tr H cos δθ = tr {{ H + trt, 18 cos δφ using the fact that the elements of the algebra are skew-symmetric. The other relation we need to have the full information contained by 1 is the planes of rotation. Since the embedded surface is -dimensional, the plane of rotation for δφ is trivial, which is 15. The plane of rotation J for δθ can be obtained by solving complex matrix linear equation 16: H I cos δθ J = sin δθ. 19 In conclusion, the contracted Gauss-Codazzi equation concerning the relation of curvatures of a manifold and its submanifold in +1 dimension can be written using holonomy representation 1 or using the plane-angle representation, i.e., the relation between rotation angles by 18 and the condition for the planes of rotation, by 15 and 19. III. DISCRETE +1 GEOMETRY A. Geometrical setting Using the Gauss-Codazzi equation in terms of holonomy and plane-angle representation obtained in the previous section, we can write the Gauss-Codazzi equation for discrete +1 geometry; this is the obective of this section. Let a portion of curved -dimensional manifold M be discretized by flat tetrahedra, we call this discretized manifold M. See FIG Angle relation In this section we derive the Gauss-Codazzi equation purely from the angles between simplices, i.e, the relation between geometrical obects, without reference to coordinates this is, indeed, the reason Tullio Regge developed his Regge calculus and discrete general relativity, in his famous paper titled General relativity without coordinates [5]. There are two kinds of relevant angles in a flat tetrahedron: the angle at a vertex between two segments of the

7 7 a b a-c b-a c c-b a ca ab bc abac b c {U5,J5 {U1,J1 {U,J {U,J {U4,J4 {U6,J6 FIG.. Given angles φ ab, φ bc, φ ca at point p of a tetrahedron, we could obtain the dihedral angle θ abac. In fact, θ abac is only φ bc proected on the plane normal to segment a. a^b a ab b FIG.. Suppose Σ is a D surface discretized by three triangles. Point p is the hinge where the deficit angle δφ is located. ' δφ is the D intrinsic curvature of Σ at point p. tetrahedron meeting at a vertex, which we denote φ; and the dihedral angle at a segment, namely the angle between two triangles, which we denote θ. These two types of angles are related by the angle formula of the tetrahedron : cos θ abac = cos φ bc cos φ ab cos φ ac. 0 sin φ ab sin φ ac See FIG... Isosceles tetrahedron For simplicity of the derivation in the next section, we consider isosceles tetrahedron. An isosceles tetrahedron is build by three isosceles triangles and one equilateral triangle. In this case, all angles around one point p of the tetrahedron are equal, say φ i = φ, and the dihedral angle relation 0 becomes: cos θ = cos φ 1 + cos φ. This implies that the dihedral angles between two isosceles triangles are equal as well, say θ i = θ. Another geometrical property of a tetrahedron which is useful for the derivation in this section is the volume. The volume of unit isosceles parallelepiped is: vol = n 1 n n = 1 1 cos φ 1 + cos φ = 1 sin φ 4 cos φ 1. 1 B. Curvatures in discrete geometry In n-dimensional discrete geometry, the notion of curvature is represented using the plane-angle representation. The intrinsic curvature is represented by the angle of rotation, which in this case is called as deficit angle, and this deficit angle is located on the hinge, the n form duals to the plane of rotation. The geometrical interpretation of deficit angle is illustrated by the FIG.. To obtain the -dimensional and -dimensional intrinsic curvatures in discrete geometry picture, firstly we need to define a discretized surface embedded in M, say Σ. The simplest way is to take the surface of the tetrahedron as a

8 8 n p n1 n s1 s s FIG. 4. A portion of a discretized curved manifold M with an embedded discretized hypersurface Σ. It is clear that the loop γ circles point p, where the D deficit angle δφ is located. δφ represents D intrinsic curvature since it is located on a point. The loop also crosses segments which are s 1, s, s. Each segments contains its own D deficit angle δθ, which represent D intrinsic curvature since it is located on a line. Using the exponential map, we can obtain holonomies on each line corresponding to the deficit angles. From these holonomies, we can obtain the D total holonomy along loop γ, which is the product of holonomies on the segments. e e1 slice see FIG. 1. Having Σ embedded on M, we consider a loop γ on Σ. We consider the holonomy related to the n D and D curvatures definede by this loop. The D curvature is defined by the D holonomy H SO SU n1 around loop γ, which describes the intrinsic curvature n of M n n, while the D curvature is the D holonomy H SO, describing the intrinsic curvature of slice Σ, still along the same loop. Taking trace of the holonomies, we obtain the angles of rotation, i.e., the amount of a components of a vector parallel to the plane of rotation get rotated when parallel transported along loop γ. In discrete geometry, this rotation angle is the same as the deficit angle on the hinges. 1. D intrinsic n1curvature SO Using FIG. 4 as the simplest case, we take loop γ circling point p on the surface Σ, therefore the D curvature is the deficit angle on p: R = δφ = π i φ i, or in terms of the cosine function: cos δφ = cos φ i cos φ i sin φ sin φ k.,i,k,<k For the isosceles tetrahedron case, - become: R = δφ = π φ, cos δφ = cos φ cos φ sin φ = cos φ 4 cos φ. We could obtain the holonomy by the exponential map 14, using the SO algebra 15.. D intrinsic curvature SU We could write D curvature by the holonomy of SU using 16, but the -dimensional holonomy along loop γ circling point p is a product of three holonomies on the three hinges, because the loop crosses three hinges. See FIG. 4. H = H i 4

9 Using 16, we obtain: H = I cos δθ i + J i sin δθ i = I cos δθ i +,i,k,<k J i sin δθ i cos δθ cos δθ k +,i<,i, k J i sin δθ i J sin δθ cos δθ k + J i sin δθ i, where each hinge segment i is dual to plane J i = J i ˆ σ and is attached by a deficit angle δθ i. Taking trace of the equation above gives: tr H cos δθ or simply: = tri + tr J i + J 1 J J cos δθ cos δθ i + = But Ji J = cos φ i, and,i,k,<k sin δθ i, cos δθ i trj i sin δθ i cos δθ cos δθ k + 0,i<,i, k tr J i J,i<,i, k J i J sin δθ i sin δθ cos δθ k Ji J sin δθ i sin δθ cos δθ k + J1 J J sin δθ i, J1 J J is ust the volume of the parallelepiped defined by 9 { J1, J, J. Remember that J = J, or 1-form J is the Hodge dual of the -form J in -dimensional vector space. The direction of J is normal to the plane discribed by J. Therefore, the D intrinsic curvature is: 1 or cos δθ = cos δθ i,i<,i, k R = δθ = arccos cos δθ i,i<,i, k cos φ i sin δθ i sin δθ cos δθ k + vol cos φ i sin δθ i sin δθ cos δθ k + vol sin δθ i, 5 sin δθ i, with vol is the volume form described in 1. For isosceles tetrahedron case, we have φ i = φ and δθ i = δθ, so the holonomy is: H = I cos δθ + cos δθ sin δθ J i ˆ σ + cos δθ δθ sin Ji ˆ σ J ˆ σ + sin δθ Ji ˆ σ, i,=1,i Therefore, the D intrinsic curvature for isosceles tetrahedra case is: R = δθ = arccos cos δθ cos φ cos δθ sin δθ + vol sin δθ.. D extrinsic curvature Equation 18 is a linear relation between the trace of D holonomy with the D holonomy and its residual part T. Inserting and 5 to 18, we could obtain trt, but this is not the extrinsic curvature, since T = 1 5 has a similar form with, except the cos φ i and the term which contains the volume form.

10 10 {U5,J5 {U1,J1 {U {U,J {U4,J4 FIG. 5. From the dihedral angle relation, we could obtain the internal dihedral angle θ on a segment. α is the external dihedral angle, while θ is the angle between normals of the two triangles. S + [K, K] + I I. Next, we would like to obtain an angle describing the extrinsic curvature K. We do not expect the relation between this angle with R and R to be linear, since linearity is held on the holonomy representation of Gauss-Codazzi 1, while the relation between their rotation angle R and R does not need to be linear. Consider the isosceles tetrahedron case. On each segment of the tetrahedron, lies the D deficit angle δθ. This deficit angle is defined as: δθ = π θ + α, where θ is the internal dihedral angle, and α is the external dihedral angle coming from the dihedral angles of other tetrahedra remember that the full discretized + 1 picture is the 1-4 Pachner moves. Let s introduce the quantity θ = α π. 6 For the case where M is flat, δθ = 0. This causes θ + α = π, and using definition 6, we obtain: θ + θ = π. In this flat case, it is clear that θ is the angle between the normals of two triangles, see FIG. 5. Therefore, θ is in accordance with the definition of extrinsic curvature in, where K is defined as the covariant derivative of the normal n to the hypersurface Σ. Because of this reason, we define θ as the D extrinsic curvature, since in a general curved case, it will inherit the curvature of the D manifold. C. Discrete Gauss-Codazzi equation in +1 dimension A simpler but equivalent way to obtain the D instrinsic curvature is to use the purely geometrical picture of the discretized manifold. See FIG. 6. From Figure 6, we obtain the proected D intrinsic curvature of each segment on the artificial segment dual to the total plane of rotation as: Using cos δθ γ = sin χ cos δθ. 7 1 cos χ cos δθ J i = sin χ ˆk + cos χˆri, we obtain: Ji, J = cos φ = sin χ + cos χ cos π, or: cos χ = 1 cos φ, 8 sin χ = cos φ, 9

11 11 ' FIG. 6. The D deficit angle on each segment can be thought as a piece of space inserted along the hinge. The size of this piece of space is described by the deficit angle δθ, which is located on the plane normal to the hinge. The total deficit angle δθ obtained by adding three rotation in 4, is located in the plane normal to an artificial hinge in the direction ˆk, since we use isosceles tetrahedron. Therefore, proecting δθ to this plane, we obtain δθ γ. since we use isosceles tetrahedron. All the curvatures: D intrinsic, D intrinsic, and D extrinsic curvatures are contained in 7, through: δθ γ = δθ = R, 0 φ = π δφ = π R 1 cos φ δθ = π θ + θ = π arccos + K. 1 + cos φ Using the half-angle formula and inserting 8- to 7, we obtain the relation between angles of the discrete Gauss-Codazzi equation for a special discretization using isosceles tetrahedra: cos π R 1 + cos R 1 cos K + arccos = 1 + cos π R cos R 1 cos π R Together with their planes of rotation: in 15 for R and ˆk σ for R since we use isosceles tetrahedron, is the Gauss-Codazzi equation in plane-angle representation, which contains same amount of information provided by the linear Gauss-Codazzi equation in holonomy representation in 1. IV. SPIN NETWORK STATES AND D LOOP QUANTUM GRAVITY A. -complex and its slice In this section, we apply the results obtained above to loop quantum gravity LQG. We will briefly review the quantization of gravity in first order formalism using the loop representation. From the real space representation, we move to the Hodge dual space [7]. The dual space is introduced in analogy to what is usually done in higher dimension. In LQG, the fundamental geometrical obect is the -complex. A -complex is a -dimensional geometrical obect dual to the discrete manifold M. Let us call the -complex as M. We can take a slice on a -complex, which is a 1-complex called as Σ, dual to the discrete hypersurface Σ of M. A 1-complex is a graph consisting oriented lines called links l and points called nodes n. The variables attached to the graphs are the variables coming from

12 ' 1 FIG. 7. The figure above is a -complex dual to the geometric figure in FIG. 1. Suppose we have a -dimensional curved manifold M discretized by four tetrahedron as the first figure above in flat case it is known as 1-4 Pachner move. Then we take an embedded slice Σ as the surface of one tetrahedron the dark blue surface discretized by three triangle. Going to the dual space picture, M is dual to a -complex M called as the D bubble [7] its geometry is defined by black lines and black points, while Σ is dual to a 1-complex Σ called as D bubble defined by purple lines and point. Clearly, the D bubble Σ is embedded on the D bubble M. The bubble graph define complete curvature of the portion of space. first order formalism: the holonomy U l SU and algebra J l su on each links l of the graph. A graph with elements of SU su is called a spin network. In the canonical quantization, we promote {U l, J l to operators {Ûl Ĵl,. The space SU su is the phase-space of LQG, and the pair {Ûl, Ĵl are the pair conugate to each other, satisfying a quantized-algebra structure as follow [, 7, 11]: ] [Ûl, Ûl = 0 ] [Ûl, Ĵ l i = δ ll σ i Û l ] i [Ĵ l, Ĵ l = δ ll ε i kĵ l k. On each node, the gauge invariant condition must be satisfied: Ĝ ψ = l n Ĵ l ψ = 0. [ This condition selects an invariant Hilbert space K = L SU l /SU n] ψ of a closed triangle formed by {Ĵ1, Ĵ, Ĵ [, 7, 11]. As an example for our case, see FIG. 7. B. Curvature operators Given a D bubble graph with a phase space variable on each link, we can construct the operators of D, D intrinsics and D extrinsic curvature from these variables. We construct these operators using the basic phase-space operators, for example, the angle operator from the phase-space variables [1]: Ĵ + Ĵk Ĵi ˆφ i = arccos, i,, k n. Ĵ Ĵk From this angle operator, we obtain the D intrinsic curvature along loop γ as: ˆR = π ˆφ i, i γ

13 1 {U5,J5 {U1,J1 {U,J {U,J {U4,J4 {U6,J6 FIG. 8. Taking only the slice on the dual space representation, we obtain the spin network Σ. γ is the loop chosen. Inside γ is where the point p located in the real space representation. which is clearly the deficit angle on a point in the real space representation. See FIG. 8. The D intrinsic curvature can be obtained from the holonomy around a loop. Remember that the direct geometrical interpretation of spacetime can only be obtained from second order formalism, while in LQG, the fundamental variables comes from first order formalism: the holonomy U comes from the curvature -form F, instead of the Riemannian curvature R see Section 1. The relation between F and R, subected to the torsionless condition is: R = e F, 4 with e is the local trivialization map between them. Contracting 4 with an infinitesimal area inside the loop carried by the derivation in Section 1 gives: since e is a diffeomorphism and must have inverse. The D intrinsic curvature is R trh, with H is the holonomy coming from the contraction of Riemanian curvature R with an infinitesimal area. But since e is a diffeomorphism and trace of the holonomy is invariant under diffeomorphism and gauge transformation, we obtain trh = tru. Therefore, we can write the D intrinsic curvature operator as: The extrinsic curvature operator can beobtained directly from the discrete Gauss-Codazzi equation : ˆK = arccos H = e U, U = e 1 H, ' trû1ûû ˆR = arccos. 5 1 cos ˆR 1 + cos ˆR 1 cos π ˆR arccos cos π ˆR 1 + cos π ˆR. C. The semi-classical limit In this section, we show that the linearity of the relation between tr H and tr H in equation 18 is recovered from the spin network calculation in the semi-classical limit, by dropping the Maslov phase. Let us take a simple spin network illustrated in FIG. 9. For simplicity in our calculation, we use a closed graph as an example, but the result is valid for any spin network graph. The vector and covector state of the spin network Σ in the algebra representation are: ψ = i m6mn i m5mn1 i m4m1n i n6n5n4 1 m 1 n m 6 n 6, 6 m 1...m 6 n 1..n 6 m ψ = i 6 m n i m 5 m n 1 i m 4 m 1 n i n 6 n 5 n 4 1 m 1n m 6n 6, 7 m 1...m 6 n 1..n 6 with i is the intertwinner attached on the node. See [7, 11, 1] for a detail explanation about the spin network state. Changing the representation basis using: U mn = D mn U,

14 14 J4 J J J J5 J6 6 a l FIG. 9. A simple spin network Σ the graph in purple, which will give a geometrical interpretation of a fuzzy quantum tetrahedron in the semi-classical limit. This spin network can be obtained from the spin network in FIG. 8, by identifying their open links together. In D theory, the nodes purple points are dual to triangles, and the links purple lines are dual to segments. with D mn U is the matrix representation of U SU in + 1-dimension, we can write the vector and covector states in the group representation basis: ψ = m 1...m 6 n 1..n 6 i m6mn i m5mn1 i m4m1n i n6n5n4 du 1..dU 6 6 D i m in i U i U 1..U 6, m ψ = i 6 m n i m 5 m n 1 i m 4 m 1 n i n 6 n 5 n 4 m 1...m 6 n 4..n 6 du 1..dU 6 6 D i m U i U 1..U 6. i n i We choose a specific closed path on spin network Σ, say, loop γ, defined by links {1 See FIG. 9. The holonomy U = U 1 U U are attached along γ, therefore we have a well-defined operator related to the D intrinsic curvature operator 5 on the loop: tr H = Since the matrix representation satisfies: tru cos ˆR = tr U 1 U U = m,n δ mn D mn U 1 U U. 8, m ˆD U 1 U U, n = D mn U 1 U U = k,l D mk U 1 D kl U D ln U, 9 we can write tr H = δ mn D mk U 1 D kl U D ln U. 40 k,l,m,n The geometrical interpretation of the spin- in 8-40 is the measure on the artificial hinge dual to the plane of rotation inside loop γ. See FIG. 6 and Section III C. Acting tr H with the spin network states gives: ψ tr H ψ = i m6mn i m6m n i m lm i n ln i m5mn1 i m5m n 1 i m km i n 1kn 1.. m 6,l,m,n,m,n... m 4,m,m 1,n,m 1,n m 5,k,m,n 1,m,n 5 i m4m1n i m4m 1 n i n mn i m 1mm 1 i n6n5n4 i n6n5n4, n 4..n 6 but the first three-terms in the parantheses are only the Wigner 6-symbols: { { { ψ tr 6 H ψ = i n6n5n4 i n6n5n4. 41 n 4..n 6 {6 i

15 15 J4 J J J J5 J6 6 a l 6 a 6 { 6 FIG. 10. The 6-symbol can be represented using a graph left, where it can be translated into its dual right which is a tetrahedron. Node a on the graph is the triangle a on the tetrahedron. θ l is the dihedral angle on. We are ready to take the semi-classical limit by setting all the i s to be large. Using the result proved by Roberts [14, 15], for large spin, the 6-symbol can be approximated as follow: { 6 1 {6 = 1 cos 4 5 i θ i π, 6 1π V 4 where the six i constructs a tetrahedron, θ i are the internal dihedral angle on segment i, and V is the volume of the tetrahedron. The π 4 term is known as the Maslov phase. Using this result to our case, we obtain: { 6 1 cos 6 θ 6 + θ + θ + θ + θ + θ l π 1π V 4 1 cos 6 θ 6 + θ + θ + θ l π. 4 1π V 4 where the last step we use the fact that i = i, noting that i comes from the covector state 7. This fact also cause θ i = θ i. See FIG. 10. Remember that is the measure on the artificial hinge inside loop γ, in other words, the spin number is the length of the artificial hinge. Along this artificial hinge, there is a point p where the D deficit angle is located. Taking i, we have: See FIG. 10. Using this approximation, 4 becomes: { 6 θ 6 0, θ π, θ π, θ l φ = φ 6, Doing the same way for the other two 6-symbols, we obtain: { { By writing the cosine function using the Euler formula: 1 1π V1 6,,, cos + π + φ 6 π 4 1 1π V 5,, 1, cos + 1 π + φ 5 π 4., 1 1π V 4, 1,, cos 1 + π + φ 4 π. 4 cos x = 1 exp ix + exp ix, we can write the products of the three 6-symbols as: {6 i π 1π V 1 V V cos... + cos π + φ 4 φ 5 + φ 6 π cos 1 π + φ 4 + φ 5 φ 6 π π + φ 4 + φ 5 + φ 6 π 4 + cos π + φ 4 + φ 5 + φ 6 π 4. 4

16 From this point, we take the sum of to be an integer/natural number we discard the half-integer possibility and choose the spin number = 1, then we could write: {6 i π 1π V 1 V V cos δφ π + O, 44 4 where we denoted the remaining terms by O and using φ 4 + φ 5 + φ 6 = π δφ see FIG. 9. Inserting 44 to 41, we obtain: ˆR π cos cos ˆR. 4 The linearity relation 18 can be recovered by dropping the Maslov phase: tr H tr H. 16 An important point is we have a freedom to choose the spin number, since this spin number comes from 8, which is the dimension of the SU representation. If we set = 0, then, the volume of the 6 tetrahedra V i will be zero and relation 4 will diverge, given any i s. But since the range of spin is 1, with multiple of half, then the spin number acts as an ultraviolet cut-off which prevents relation 4 from divergence. V. CONCLUSION We have obtained the Gauss-Codazzi equation for a discrete +1-dimensional manifold, discretized by isosceles tetrahedra. We have studied the definition of extrinsic curvature in the discretized context. With definitions of the -dimensional intrinsic, -dimensional intrinsic, and -dimensional extrinsic curvature in the discrete picture, we have promoted them to curvature operators, acting on the spin network states of +1-dimensional loop quantum gravity. In the semi classical limit, we have shown that the linearity between tr H and tr H in equation 18 is recovered in the spin network calculation, by dropping the Maslov phase. There exist a natural ultraviolet cut-off which prevents the discrete Gauss-Codazzi equation from divergences in the semi-classical limit. [1] S. M. Carroll. Spacetime and geometry: An introduction to general relativity. San Francisco, CA, USA. Addison-Wesley. ISBN [] T. Regge, R. M. Williams. Discrete structures in gravity. J. Math. Phys. 41, arxiv:gr-qc/00105v1. [] P. Doná, S. Speziale. Introductory lectures to loop quantum gravity arxiv:gr-qc/ [4] J. W. Barrett and I. Naish-Guzman. The Ponzano-Regge model. Class. Quant. Grav. 6: arxiv:grqc/ [5] T. Regge. General relativity without coordinates. Nuovo Cim [6] R. Arnowitt, S. Deser, C. Misner. Dynamical Structure and Definition of Energy in General Relativity. Phys. Rev : [7] C. Rovelli, F. Vidotto. Covariant Loop Quantum Gravity: An Elementary Introduction to Quantum Gravity and Spinfoam Theory. UK. Cambridge University Press. ISBN [8] J. C. Baez, J. P. Muniain. Gauge Fields, Knots, and Gravity. Series on Knots and Everything: vol 4. World Scientific Pub Co Ltd. ISBN [9] H. W. Hamber. Quantum Gravitation: The Feynman Path Integral Approach. Springer-Verlag. ISBN [10] B. C. Hall. An Elementary Introduction to Groups and Representations. Graduate Texts in Mathematics 00. Springer-Verlag. arxiv:math-ph/00050v1. [11] C. Rovelli. Zakopane lectures on loop gravity arxiv:gr-qc/ [1] M. Seifert. Angle and Volume Studies in Quantized Space. arxiv:gr-qc/ [1] S. A. Maor. A Spin Network Primer. Am. J. Phys arxiv:gr-qc/ [14] J. Roberts. Classical 6-symbols and the tetrahedron. Geom. Topol arxiv:math-ph/ [15] G. Ponzano and T. Regge. Semiclassical Limit of Racah Coefficients. Spectroscopy and Group Theoretical Methods in Physics. Amsterdam. pp