in our previous paper @xcite , we extended the notion of vassiliev invariant and other invariants associated with chern
simons theory to the context of spin networks @xcite .
we showed that these invariants are loop differentiable in the sense of distributions and we defined an infinitesimal generator of diffeomorphisms in terms of the loop derivative .
the generator correctly annihilates diffeomorphism invariant states . in this paper
we will introduce a hamiltonian constraint based on the loop derivative .
we will limit the discussion to invariants of trivalent vertices and will only concentrate on the `` euclidean '' part of the hamiltonian constraint in the sense of barbero @xcite and thiemann @xcite . limiting ourselves to trivalent intersections
is clearly unphysical since the volume operator identically vanishes on such states . however , the essence of most calculations is already present with trivalent intersections and the calculational difficulty is significantly lower than with four or higher valent intersections .
we will come back to this issue in the discussion section where we will highlight which results are quite plausibly going to survive when one goes to higher valent intersections .
we will then extend the space of invariants through the construction of `` habitats '' such that the resulting states are not diffeomorphism invariant . based on these states
we will check the consistency of the quantum commutator algebra of diffeomorphisms and hamiltonian .
we will show that there are no anomalies , therefore constituting a consistent theory of canonical quantum gravity .
we will start in the next section by introducing the hamiltonian constraint of quantum gravity in terms of the loop derivative .
this derivation will be `` generic '' in the sense that we will motivate it through formal manipulations of the loop transform and therefore will not make specific assumptions about the quantum state in the connection representation on which the operator is acting upon . without being more precise about the space of states upon which it acts , we can not guarantee that the operator is well defined .
we will show that on several spaces related to the vassiliev invariants the operator is indeed well defined .
the results obtained will be later applied in specific `` habitats '' in which we will show that the operators exist and we can compute correctly the commutators .
this we will discuss in section iii .
the need to introduce new ( non - diffeomorphism - invariant ) `` habitats '' is based on the fact that the spaces we have considered in the companion paper are all diffeomorphism invariant and therefore one can not explore in a non - trivial way commutators involving the diffeomorphisms .
we will introduce an explicit example of such a habitat in detail , the `` functions with marked vertices '' . in the end
we will see that we will be able to recover the classical poisson algebra at the level of quantum commutators , but there will be subtle issues involved in the definition of the right hand side of the commutator of two hamiltonians , which we will discuss in detail .
the last section will be devoted to discussing the implications of the level of consistency achieved by the theory .
we wish to introduce a quantum version of the hamiltonian constraint of canonical general relativity .
following thiemann @xcite , one can construct this operator introducing the real version of the ashtekar variables first discussed by barbero @xcite . in this formulation ,
the canonical pair consists of a set of ( densitized ) triads @xmath0 and as conjugate variables real - valued @xmath1 connections @xmath2 .
the hamiltonian constraint of canonical , real , lorentzian general relativity in this framework consists of two terms , @xmath3 where @xmath4    in these expressions , @xmath5 is the newton constant and in ( [ k ] ) @xmath6 means that the smearing function @xmath7 is unity .
this remarkable formulation implies that the hamiltonian constraint is composed by two terms . the first one , @xmath8 , coincides with the hamiltonian constraint one would obtain in a canonical formulation of general relativity on a manifold of euclidean signature .
the lorentzian theory is therefore attained through the addition of an extra term .
both terms are given as expressions involving the poisson bracket of a connection with the volume of the space @xmath9 and with the function @xmath10 .
the latter can be in turn obtained as a poisson bracket of the euclidean part of the hamiltonian constraint with the volume .
in this paper we will only concentrate on realizing at a quantum mechanical level the euclidean part of the hamiltonian .
if the reader wishes , we are considering canonical quantum gravity in the euclidean context .
however , it should be noticed that the additional term in the hamiltonian is obtained by successive commutators of the operator we will consider , with the volume operator .
this suggests that one could generalize our construction to the lorentzian case rather straightforwardly , but this has not been studied in detail .
the work of thiemann also shows that similar structures arise if one couples the theory to matter .
we will in this paper only concentrate on the vacuum case but again one does not see a priori impediments to extend our treatment to general relativity coupled to matter .    to implement the hamiltonian constraint we will partially use the same procedure proposed by thiemann .
we will not describe it again in full detail here , we just list some of the salient features .
one starts by writing the classical expression ( [ hame ] ) as @xmath11 f_\epsilon(y , x)\ ] ] where @xmath12 is a regularization of the dirac delta function , i.e. @xmath13 . strictly speaking
, this expression is only gauge invariant in the limit , so to preserve gauge invariance in the regularization procedure we will `` join '' the @xmath14 with the connection using an infinitesimal piece of holonomy @xmath15 along an arbitrary path @xmath16 going from @xmath17 to @xmath18 , which tends to the null path in the limit @xmath19 . to unclutter the notation we will therefore denote @xmath20 , @xmath21 f_\epsilon(y , x ) \label{uno}\ ] ]    in order to prepare the classical expression for quantization , we triangulate the space - like hypersurface @xmath22 in terms of elementary tetrahedra @xmath23 . the triangulation has the following properties : we select a finite set of distinct points of @xmath22 , denoted as @xmath24 . at each of these points
we choose three independent directions @xmath25 and construct the eight tetrahedra with vertex @xmath26 and edges @xmath27 , with @xmath28 .
the eight tetrahedra saturating @xmath26 define a closed region @xmath29 of length scale @xmath30 .
the remaining region @xmath31 is triangulated by arbitrary tetrahedra @xmath32 .
the motivation for this peculiar discretization is that when we promote the classical expression to a quantum operator , we will adapt the triangulation to the spin network of the state in question by choosing the points @xmath26 to coincide with the vertices of the spin network . in terms of this triangulation of space
we can write equation ( [ uno ] ) in the form , @xmath33 f_\diamond(y)\,,\label{hamt}\ ] ] where the regions @xmath34 indicate either a box @xmath29 or a tetrahedron @xmath32 , and @xmath35 is any point interior to @xmath34 .
conventionally we will choose for the boxes @xmath36 , and for the tetrahedra @xmath37 will represent one of its vertices . to discretize the integral
we have introduced the volume of each region @xmath38 , @xmath39 where @xmath40 if @xmath41 , and @xmath42 if @xmath43 . in this last case ,
the @xmath44 s represent the edges of the tetrahedra @xmath32 adjacent to @xmath37 .
we also have adapted the regularization of the dirac delta function to the tetrahedral decomposition by defining , @xmath45 where @xmath46 is one if @xmath47 and zero otherwise .
we now replace in ( [ hamt ] ) the @xmath48 using formula ( [ vol ] ) and we represent @xmath49 using a holonomy in the fundamental representation along the edge @xmath50 of the triangulation , @xmath51f_\diamond(y)\ , , \label{classham}\ ] ] and we retraced the path described by the holonomy to preserve gauge invariance .
it should be noticed that this retracing does not contribute to the expression at leading order in @xmath30 ( the scale parameter of the regions @xmath34 ) .
we will now proceed to study the quantization of the last expression . for
that we need to `` adapt '' the triangulation we introduced to a spin network state @xmath52 @xcite .
let @xmath53 be the set of vertices of the spin network and let @xmath54 be the triples of non - coplanar edges incident at @xmath55 ( @xmath56 , being @xmath57 the number of such triples ) .
to adapt the triangulation to the spin network state we perform two operations : first we identify the points @xmath26 of the boxes @xmath29 with the vertices @xmath55 and second , given a triple of edges @xmath54 incident at @xmath55 , we orientate the three unit vectors @xmath58 along the tangents of the edges at @xmath55 .
then we can write , @xmath59 to simplify the notation we identify from now on @xmath60 . as in the previous paper , we assume we are acting on a state given by a loop transform , @xmath61 w_a[s]\ , , \label{states}\ ] ] and we realize the hamiltonian operator over the spin network wavefunctions by promoting the classical expression ( [ classham ] ) as an operator acting on the wilson net appearing in the loop transform , very much as we did in the companion paper when we discussed the diffeomorphism operator , @xmath62 \right)f_{\box_v}(y ) w_a(s)\,.\nonumber\end{aligned}\ ] ] a first observation is that the volume operator has non - vanishing contributions only at the vertices of the spin net , so we replace the sum in the hamiltonian over all @xmath34 s by a sum over all the vertices of the spin net .
the second observation is that , for an n - valent vertex , we have taken the average over all the non - coplanar triples of edges associated with this vertex .
finally notice that only one of the terms in the commutator contributes , that with the volume operator on the left .
the one with the volume operator on the right vanishes since it is proportional to the trace of the lie algebra element @xmath63 .
let us now consider a generic @xmath64valent vertex and study the action of the operator on the wilson net .
we assume that the holonomies are all outgoing from @xmath26 .
schematically , we are trying to represent the action of @xmath65 , with @xmath66 a generic edge of spin @xmath67 incident at @xmath26 .
an @xmath64valent vertex is characterized by an intertwiner that can be represented through a vector @xmath68 of @xmath69 of irreducible representations .
the factor @xmath70 acts by adding a line of spin @xmath71 ingoing from the vertex @xmath26 and parallel to the line @xmath66 , @xmath72 notice that the insertion of the holonomy leaves a wilson net that is not gauge invariant anymore , we represent this by keeping the group indices @xmath73 in the diagram .
if one was acting on an @xmath64valent vertex characterized by an intertwiner @xmath68 , this leaves an @xmath74valent vertex characterized by intertwiners @xmath75 ; @xmath67 being in the diagram the spin of the line connecting the original vertex to the point ( infinitesimally nearby ) where the holonomy was inserted . notice that @xmath67 coincides with the spin of the original line @xmath66 .
we denote that the line is infinitesimally close by the dashed circle .
we now act with the volume operator . for that we first have to reduce the product of holonomies associated with @xmath66 to a superposition of irreducible representations .
this procedure defines a new intertwiner for the vertex @xmath26 over which the action of the volume operator is well - defined , @xmath76 the action of the volume operator can be represented by a matrix that rearranges the intertwining of the vertex .
this matrix elements would depend also on the spins of the external edges of the vertex ( notice that @xmath10 is the color of one of this external edges ) .
schematically , @xmath77 where in the last step we reduce the line of spin @xmath71 of the holonomy using recoupling theory .
one is left with a double insertion at the vertex , which we can rearrange via recoupling , since everything is happening infinitesimally close to the vertex .
we write the result in the following way , @xmath78 the above expression has now to be contracted with @xmath14 at the open strands labeled with @xmath79 and @xmath80 .
one can generate @xmath14 by introducing an extension of the idea of loop derivative for spin networks .
the reason why one does not simply recover the ordinary loop derivative is that the edge before and after the point of the insertion of @xmath14 are in general in different representations @xmath67 and @xmath81 .
the ordinary loop derivative inserts an @xmath14 without changing the representation of the original line .
this definition of loop derivative appeared as natural in the context of loops , however in terms of spin networks one expects in general to have a situation like the one we have here .
one can use the same set of constructions we did for the ordinary loop derivative of the invariant @xmath82 to compute the action of this generalized loop derivative is proportional to the inverse of the coupling constant of chern - simons theory . for
a definition of the invariant @xmath82 , see the companion paper , section ii.c . ] .
however , we will not pursue the exploration of this generalized derivative in this paper .
instead we will concentrate on the case of trivalent intersections .
if the spin networks considered have trivalent intersections , when we repeat the construction we did above , the representations of the edges before and after the insertion of @xmath14 are the same .
the operator we then obtain to represent @xmath14 corresponds in this case to the ordinary loop derivative .
we have now to consider the action of the volume operator over a four - valent vertex with only three real edges @xcite , @xmath83 with @xmath84 where @xmath85 , and @xmath86 . using this result , the matrix elements of equation ( [ nval ] )
are reduced to the following simple expression for trivalent vertices ( @xmath87 in this case ) , @xmath88 where @xmath89 , and @xmath90 and @xmath91 are the other spins of the edges incident at the trivalent vertex .
therefore the spin of the link going from the intersection to the point of insertion of @xmath14 is unchanged through the action of the volume operator . by inserting the @xmath14
we obtain the final action of the hamiltonian constraint on trivalent vertices ,    @xmath92    the superscript on the loop derivative denotes at which of the incident edges at the intersection it adds the extra strand @xmath16 .
notice that since @xmath93 is summed by the einstein convention , the derivative ends up acting on all edges ( we have abused of the summation convention allowing the index @xmath93 to be repeated three times in the last expression ) .
the vector @xmath94 labels the spins of the tree edges of the spin network incident at @xmath26 ( @xmath95 is the spin of the edge @xmath96 ) , and the group factor @xmath97 is given by the expression , @xmath98 where the indices @xmath99 take cyclic values in the set @xmath100 .
we have therefore a general expression for the hamiltonian constraint acting on any loop - differentiable wavefunction with support on trivalent spin networks , that bears a relationship with a wavefunction in the connection representation @xmath101 given by ( [ states ] ) .    at this point
it is worthwhile comparing this action of the hamiltonian constraint we just introduced with the original proposal of a hamiltonian constraint in the loop representation ( doubly densitized ) in terms of loop derivative as introduced in @xcite .
such constraint was written in @xcite as , @xmath102 the hamiltonian acted on functions of loops and the operator @xmath103 had the action of re - routing one of the lobes of the partition of the loop determined by the points @xmath104 .
we see two fundamental differences with the operator we just introduced .
the first one is the re - routing operator @xmath103 .
this operator arose in the loop representation to account for the fact that states based on loops in the fundamental representation do not diagonalize operators like the volume . in the spin network context
the operator is replaced by the group factors we discussed .
more important is the difference concerning density weights . in the hamiltonian in terms of loops
we see that one obtains a doubly densitized quantity by considering the product of a dirac delta times the regulator .
to understand properly this difference , it is worthwhile considering on which spaces of functions are these operators meant to operate upon . in the hamiltonian in terms of loops one
had in mind that the loop derivative was acting on functions of loops such that the result was a smooth function .
an example of such functions would be holonomies built with a smooth connection .
the resulting expression would then consist of a smooth tensor , @xmath14 , a dirac delta integrated on three dimensional space and , in the limit , a two dimensional dirac delta given by @xmath105 integrated along the two one dimensional integrals .
the result is finite but is regularization dependent , since the two dimensional dirac delta has an inverse power of the three dimensional volume in it that is not compensated .
if we now turn our attention to the hamiltonian we introduced in this paper , and consider its action on a state of the form @xmath82 , we will see a different behavior . given the action ( companion paper : formulas 71 and 44 ) of the loop derivative on @xmath82 we get , @xmath106 where the extra edge added by the loop derivative starts and ends in the edge @xmath107 and as before , @xmath99 take cyclic values in the set @xmath100 .
the reason for this is that if it ended in any other edge incident on the vertex one would ( in the limit @xmath108 ) get zero for the action since one would have an integral repeated over one of the edges contracted with @xmath109 .
the action on @xmath110 of the loop derivative that appears in the hamiltonian constraint , for trivalent vertices reduces therefore to a chord diagram . if we evaluate the chord diagram using recoupling identities , we finally get , @xmath111 at this point it is worthwhile continuing the comparison with the doubly densitized case .
if we evaluate the integrals involved in the above expression we see that we have six one dimensional integrals and two three dimensional dirac delta functions . in the limit
we therefore have a finite result , that is defined independent of the background structures introduced by the regulator @xmath112 .
one first performs the integral in @xmath18 , which fixes @xmath113 and then one is left with three one dimensional integrals with three step functions ( two explicit , one present in the @xmath114 ) .
the result of the integrals is @xmath115 , which cancels the denominator of the @xmath116 ( @xmath117 ) and one is left with a finite result , that in the limit @xmath108 gives , @xmath118 where @xmath119 is the factor group associated with the vertex @xmath26 .
we have introduced a quantum version of thiemann s singly - densitized classical hamiltonian in terms of the loop derivative that has the remarkable following properties :    \a ) it has a significant resemblance to the original doubly - densitized hamiltonian proposed in the loop representation @xcite , which in turn closed the constraint algebra formally @xcite .
we will discuss the implications of this in the next section .
\b ) the hamiltonian introduced is finite on the space of vassiliev invariants we introduced in the previous paper .
this happens due to the detailed form in which the loop derivative acts on these states .
\c ) the hamiltonian may be well defined on other spaces of states .
for instance , on the space of functions of spin nets obtained by considering the wilson nets along smooth holonomies ( the original kind of function that was first considered as loop differentiable ) , it is straightforward that the hamiltonian identically vanishes ( essentially , since the loop derivative is smooth , the integrals that before gained contributions due to the distributional character of the loop derivative now vanish ) .
\d ) the last two points show that something surprising is happening here , in the sense that _ we did not construct the hamiltonian in an ad - hoc way to obtain properties b ) and c)_. that is , we introduced a discretization of the singly - densitized hamiltonian and it naturally turned out to be well defined on the vassiliev invariants and to vanish on holonomies of smooth connections .
this might be pointing to a certain `` naturalness '' of the vassiliev invariant space in the context of quantum gravity . in spite of the a priori huge ambiguity in the regularization of the singly - densitized hamiltonian , it is difficult to imagine a regularization involving the loop derivative that would not be finite on the vassiliev invariants .
in addition , we have not introduced in our construction any ad - hoc `` renormalization '' of the operator to obtain finiteness ( as one had to do in the doubly - densitized case ) .
we will now proceed to discuss the issue of the constraint algebra .
the operators we introduced are written purely in terms of the loop derivative and integrals along edges of the spin net .
one can discuss the successive action of such operators without making specific reference to a space of states on which one is acting upon .
such a calculation would be `` formal '' in the sense that terms arising in it could fail to be well defined when acting on particular spaces of states .
this is the underlying philosophy of the kind of ( unregulated ) formal calculations some of us undertook in reference @xcite .
such an approach can be pursued with the regulated operators we introduce in this paper .
the fact that the operators are regulated and that they might act on the space of diffeomorphism invariant states will imply certain departures from the specifics of the calculations of reference @xcite .
if one pursues calculations at this level of generality , one can not expect to reproduce the classical poisson algebra at the level of quantum commutators without making some assumptions about the distributional behavior of the various quantities involved in the calculations .
one will recover terms that correspond to the classical ones but generically there will be additional terms .
when one particularizes to a certain `` habitat '' of wavefunctions , these extra terms in many cases will vanish .
however , it is not excluded that on certain habitats pathologies could appear . it would be an interesting exercise to pursue the calculation in general and then particularize to various `` habitats '' .
this , however , is a complex and lengthy calculation in terms of loop calculus . we will not attempt a calculation at this level of generality here . in this paper
we will concentrate on computing the constraint algebra on a series of habitats on which we recover the classical constraint algebra at the level of quantum commutators .
these habitats will include as particular cases spaces of diffeomorphism invariant functions on which the appropriate algebra will be recovered .
a point to be addressed when studying the constraint algebra is that , given the action of the hamiltonian constraint we have introduced , if one starts with a function of spin networks , when one acts with a hamiltonian constraint one ends with a more general object .
this can be explicitly seen for instance in equation ( [ nu ] ) .
what we see is that the result of the action of a hamiltonian on a function of spin networks is a function of spin networks times a `` vertex function '' , that is , a function dependent on the position of the vertex where the hamiltonian acts . if one wishes to compute the constraint algebra , the action of the second operator in the calculation of the algebra therefore takes place on such a space of functions .
one can view these functions as simply more general functions of spin networks , since after all the information of the position of the vertices comes with a given spin network .
this is possible , but one has to exhibit that dependence -of the vertices on the spin networks- in an explicit enough way that operators like the loop derivative have an appropriate action on such a dependence .
one is initially tempted to say that the loop derivative simply ignores the extra factors , since they do not appear to contain an explicit loop dependence .
but the position of the vertices depends on the edges of the spin network ( more precisely , a vertex is defined by the intersection point of at least three edges of the net ) , and therefore the vertex function should be affected by the loop derivative operator . to operate with the loop derivative over this kind of functions ,
we need to make more explicit the implicit functional dependence @xmath120 $ ] of the vertex function .
we limit the analysis to the case of trivalent vertices . given three arbitrary open paths @xmath121 ( @xmath122 ) with a common origin @xmath123 , we define the quantity , @xmath124 where @xmath125 denotes generically the three curves @xmath126 , @xmath127 is @xmath128 if @xmath129 and zero otherwise , and @xmath130 is a scalar function defined on the manifold .
it is assumed that the three paths are oriented and we denote by the overbar the path in the reversed orientation .
we also assume that , starting at @xmath123 , the three curves overlap in a finite segment and then they separate following disjoint paths , as it is shown in the figure .
we call @xmath26 the bifurcation point of the three curves .
it is immediate to show that , @xmath131 and , @xmath132 the above expressions can be used to calculate the loop derivative of a vertex function . given a vertex of the net with edges @xmath133 , we identify the portion of the curve @xmath121 going form @xmath26 to the end point with the portion of the edge @xmath96 from the vertex to some intermediate point .
this automatically guarantees that the point @xmath26 defined through the family of curves @xmath121 coincides with the chosen vertex of the graph . if the end points of the edges @xmath134 lie outside of the ball of radius @xmath30 centered at @xmath26 , then one has the identity @xmath135
.    we can now proceed to compute explicitly the loop derivative . for any function @xmath136 we have ( @xcite chapter 1 ) , @xmath137z } f(z).\label{loopder}\ ] ] using this result it is straightforward to calculate the loop derivative of @xmath138 , @xmath139z } m(z ) \theta_{\epsilon}(z - y)\right .
\nonumber\\&&\left.+ ( m(z)+m(y ) ) \delta(\epsilon-|z - y| ) \frac{\delta^d_{[b}(y - z)_{a]}}{|z - y|}\right\}.\label{derr}\end{aligned}\ ] ] the loop derivative of @xmath140 is evaluated using the above result and the leibniz rule .    now that we know how to compute the loop derivative of a vertex function , we can consider the construction of a `` habitat '' of functions of spin networks with marked points of the definite following kind : @xmath141 where @xmath142 is a group factor depending on the spins of the edges incident at @xmath26 .
this type of dependence on the vertex is precisely the one that appears in the action of the hamiltonian constraint ( [ hame2 ] ) . in the last identity
, we express the relationship between @xmath143 and @xmath144 through an operator @xmath145 .
it is clear that , in the limit @xmath19 , one gets as a result a function of the spin network times a function of the vertices and incoming spins , @xmath146 to define the action of the constraints on these kinds of functions we will study the action of the constraints before taking the limit @xmath19 in these expressions in a suitable way .
let us start by considering the diffeomorphism constraint , @xmath147\,.\end{aligned}\ ] ] in the second term the diffeomorphism acts on @xmath144 leaving all the factors unaffected .
this term is immediate to evaluate ( and in the case of diffeomorphism invariant @xmath144 s vanishes ) , so we will not discuss it further . in order to analyze the first contribution
it is convenient to make the following intermediate calculation , @xmath148z } m(z ) \theta_{\epsilon}(z - y)\right . \nonumber\\ & & \left.+(m(z)+m(y ) ) \delta(\epsilon-|z - y|)\frac{\delta^d_{[b}(y - z)_{a]}}{|z - y|}\right\}+ \mbox{\footnotesize{cyclic permutations in $ e_1,e_2,e_3$}},\nonumber\end{aligned}\ ] ] where we have used ( [ derr ] ) . in the evaluation of this quantity , the result depends on the order in which the limits are taken .
if one takes the limit @xmath19 first , one finds that the diffeomorphism annihilates @xmath149 .
this means that @xmath150 . on the other hand ,
if we take the limit @xmath151 first one gets the result , @xmath152 r_{c}(1,\epsilon , e ) + o(\epsilon),\ ] ] and , @xmath153 with this choice we see that the action of the diffeomorphism is what one would have expected geometrically : it lie - drags the dependence of @xmath91 on the position of the vertex .
this action immediately ensures that the action of the diffeomorphism constraint closes the appropriate commutator algebra on the space of functions we are considering , provided @xmath144 is either invariant ( like in the case of the vassiliev invariants ) or is such that the diffeomorphism constraint has the natural geometric action . in the general case ,
@xmath154    based on these calculations for the diffeomorphism , let us now explore the action of the hamiltonian constraint on the habitat we have been discussing . from ( [ hgen ] ) and ( [ derr ] ) we get ( @xmath155 is now the scale parameter of the triangulation adapted to the spin network ) , @xmath156z } m(z ) \theta_{\epsilon}(z - y)+ ( m(z)+m(y ) ) \delta(\epsilon-|z - y| ) \frac{\delta^d_{[b}(y - z)_{a]}}{|z - y|}\right\ } + \mbox{\footnotesize{cyclic permut . in $ e_1,e_2,e_3$}},\label{hr }
\nonumber\end{aligned}\ ] ] where we have written @xmath157 in order to make explicit the dependence of the triangulation with @xmath155 .
notice that ,    @xmath158    and , @xmath159 irrespective of the order in which the limits are taken .
this means that , @xmath160 and taking into account ( [ psim ] ) , @xmath161 } h(n)o(m,\omega ) \psi ( s ) = o(m,\omega ) h(n ) \psi(s).\ ] ] this is an important result .
it will imply that the hamiltonian constraints will commute with each other .
we will discuss this in detail in the next section .      in the classical theory , the poisson bracket of two smeared hamiltonian constraints
is given by , @xmath162 with the one - form @xmath163 and @xmath164 is the ( unsmeared ) diffeomorphism constraint .
that is , the right hand side of this expression is a diffeomorphism generated by the vector obtained by contracting the one - form @xmath165 with a double - contravariant metric .
if one wishes this expression to have a quantum counterpart , it implies realizing the product of the doubly contravariant metric times a diffeomorphism as a quantum operator . in order to do this ,
we need to re - express it in terms of variables that are suitable for the ashtekar formulation .
this is done via the following classical identities , as discussed by thiemann @xcite , @xmath166 and we recall the relationship of the covariant ( undensitized ) triads with the more usual contravariant densitized triads , @xmath167 .
we can now use the key identity , due to thiemann , @xmath168 to obtain , @xmath169    this expression can not be easily promoted in a direct way to an operator due to the denominators involving the determinant of the metric .
we can achieve this noting that the poisson bracket of the connection and the volume depends only on local information of the volume operator in the surrounding of the point @xmath17 .
we can therefore replace in the above expression the volume of the whole manifold @xmath9 , by a `` localized '' volume ( as discussed by thiemann @xcite ) @xmath170 , given by , @xmath171 where @xmath172 is an infinitesimal closed region around the point @xmath17 and @xmath173 is the characteristic function associated with this region .
the following identity holds , @xmath174 and in addition , @xmath175 where @xmath176 is the euclidean volume of the infinitesimal region @xmath177 given by a fiducial flat metric ( we will see that the final result does not depend on the fiducial metric ) .
the above expressions allow us to absorb the denominators involving the determinant of the metric into poisson brackets . using this results we write ( [ metc ] ) in the form , @xmath178 we now discretize the integral introducing the triangulation @xmath179 of space defined in section ii , @xmath180 where we have identified @xmath181 to simplify the notation . using ( [ vol ] )
we replace the @xmath182 and @xmath183 in terms of the volume of the elementary regions and the edges @xmath184 of the tetrahedra , and we join four of the @xmath184 s with the @xmath185 s to construct holonomies along the edges of the triangulation , @xmath186 where we have defined , @xmath187 { \rm tr}[\tau^k h(u_r)\{h^{-1}(u_r),\sqrt{v({\cal{r}}_v)}\ } ] \times \nonumber\\ & & { \rm tr}[\tau^l h(u_t)\{h^{-1}(u_t),\sqrt{v({\cal{r}}_v)}\ } ] { \rm tr}[\tau^m h(u_n)\{h^{-1}(u_n),\sqrt{v({\cal{r}}_v)}\}]\ , .
\label{q}\end{aligned}\ ] ] we are now ready to promote the above expression as an operator acting over the spin network wavefunctions , which we assume to be the loop transform of a state in the connection representation . as in the case of the hamiltonian constraint ,
we adapt the triangulation to the graph of the spin network choosing the points @xmath24 of the boxes @xmath188 coincident with the vertices of the spin network , and we choose the regions @xmath189 equal to that defined by the triangulation around @xmath26 but with a length scale @xmath155 . with this prescriptions we get from ( [ metc2 ] ) , @xmath190 where we have used ( [ tangent ] ) and @xmath191 with @xmath30 the length scale of the regions @xmath29 .
the action of the operator @xmath192 on the spin network wavefunctions is calculated as usual through its action on the wilson net appearing in the loop transform .
it is important to notice that this operator looses all dependence with the scale parameters involved in the triangulation and the localized volume .
this is due to the fact that , acting on the wilson net , the holonomies @xmath193 that appear in the quantum version of ( [ q ] ) generate ( in the limit @xmath194 ) finite recoupling coefficients and that the action of the volume operator is local ( it depends only on the vertex included in the region @xmath195 .    with this result at hand it is straightforward to quantize the right hand side of the poisson bracket of two smeared hamiltonian constraints . from the classical expression we immediately write , @xmath196 f_{\delta}(x , y ) \delta_{cb}^{(e)}(\pi_y^x ) \psi(s)\ , , \nonumber\end{aligned}\ ] ] where we have used the regularized expression of the unsmeared diffeomorphism operator ( see the companion paper , section ii.b ) .
introducing now the regularized metric operator we get , @xmath197 \,f_{\delta}(x , y)\times\nonumber\\ & & \int_{e_p } dw^a \theta_{\epsilon}(w;v ) \int_{e_s } dt^b\,\theta_{\epsilon}(t;v ) { q}^{ps}(v)\delta_{cb}^{(e)}(\pi_y^x)\ , \psi(s)\,.\end{aligned}\ ] ] at this point it is worthwhile mentioning that this expression has a regularization ambiguity , given by the two regularization parameters @xmath30 and @xmath155 .
the latter was introduced , as we discussed above , in the definition of the localized volume operator .
this ambiguity , associated with the fact that the localized volume and the ordinary one are the same on these spaces of functions was first noticed by lewandowski and discussed in @xcite . if one chooses @xmath198 one notices quickly that the expression for the right - hand - side vanishes .
the powers @xmath199 cancel and one is left with two integrals along @xmath200 , @xmath201 which are of order @xmath30 each .
the rest of the expression is simply the action of a diffeomorphism . assuming that the latter is finite , the expression therefore vanishes .
if we trace back the origin of this cancellation , one notices that the expression ( [ metq ] ) for the contravariant metric operator vanishes identically over spin network states of any kind , in particular , the vassiliev invariants we are considering in this paper .
this property is a consequence of the following two peculiarities of the quantization of @xmath202 in terms of spin network wavefunctions : the sum over the triangulation reduces to a sum over the vertices of the spin net ( which includes a finite number of terms ) , and in each term the result of @xmath203 is finite . notice
that in the classical result ( [ metc ] ) one has an infinite sum of @xmath204 , each of which tends to zero in the limit @xmath205 .
this limit gives in general a nonzero result .
but in ( [ metq ] ) we have instead a finite sum of terms @xmath203 , which are independent of @xmath30 and @xmath155 .
this fact alters the power counting of the factors @xmath30 and @xmath155 in such a way that all the terms tends naturally ( i.e. choosing @xmath198 ) to zero as the triangulation shrinks to a point .    that the right - hand - side may vanish was already observed in @xcite .
could one `` tune '' the limits in @xmath30 and @xmath155 so this quantity is nonvanishing ? one could , but the result would be dependent on the background structures introduced to regularize .
if one made the expression non - vanishing it would be proportional to the normalized tangent vectors at the intersection , which are background dependent .
this is not surprising : there is no naturally defined second order symmetric contravariant tensor defined in a manifold without metric .
this result is quite strong .
it implies that the commutator of two hamiltonians will have to vanish if one wishes to have consistency .
we anticipated that this would happen in the previous subsection and we will now see in detail how it happens .
as a final remark , we notice that if one wished to define the doubly covariant metric , it is straightforward to compute it using the identities introduced by thiemann and one finds that regularized with the same procedures we followed up to now , it diverges . should one worry about a theory of quantum gravity where the doubly contravariant metric vanishes and the doubly covariant metric diverges ?
we will return to this in the discussion section .
the loop derivative was originally introduced in the context of yang
mills theories , where the natural functions to act upon were holonomies of smooth connections . these functions are not the most natural ones to consider in the context of diffeomorphism invariant theories like general relativity , but our constraints are well defined on them , so we can consider them to be a `` habitat '' where to test the constraint algebra , at least as a mathematical exercise . on these kinds of functions
the diffeomorphism constraints that we have introduced here are known to close the appropriate algebra in the limit in which regulators are removed ( see @xcite , although the regularization is slightly different than the one we use in this paper , it is immediate to check that the same calculations go through ) .
we therefore will not repeat the calculation here .
the hamiltonian constraint vanishes identically on this habitat .
starting from equation ( [ hgen ] ) , if the loop derivative is a smooth function , one is left with two one - dimensional integrals along the edges of a cell of a finite function . in the limit in which the triangulation is shrank , the result vanishes .
the commutator of a diffeomorphism with a hamiltonian therefore immediately reproduces the classical result .    the right - hand side of the commutator of two hamiltonians that we introduced in the previous subsection also vanishes , here one has three one dimensional integrals along the edges of a triangulation of a quantity that goes as @xmath206 .
the final result therefore goes as @xmath207 .
therefore , as long as one chooses @xmath155 shrinking to zero faster than @xmath208 , the right hand side vanishes .    on this space
therefore , we reproduce the classical poisson algebra at the level of the commutators
. the diffeomorphisms close among themselves and the hamiltonian vanishes identically .
could one simply claim that this is the `` right '' habitat to do quantum gravity ?
the answer is that this is unlikely .
although all states are solutions of the hamiltonian constraint , this space does not contain any solution to the diffeomorphism constraint .
solutions to the diffeomorphism constraint in terms of holonomies can only be constructed as infinite superpositions , functional integrals or `` group averaging '' , and in these cases one includes connections that are not smooth .
an interesting point is that in this habitat , since one does not have solutions of the diffeomorphism constraint but has solutions to the hamiltonian constraint , the only way that one could achieve consistency in the algebra is if the right hand side of the commutator of two hamiltonians vanishes , which is the case .
another point to consider is that if one examines the expression of the commutator of two hamiltonians , although both members vanish in the limit in which one shrinks the triangulation given the smoothness of the loop derivative away from the limit the calculation is problematic .
for instance , it is not clear what is the action of the volume operator on a state given by the loop derivative of a smooth holonomy .
one might consider introducing a definition of the volume on this kind of space , but this has yet to be done in detail .
remarkably , this space of invariants also leads to a reproduction of the classical poisson algebra at a trivial level only ( as usual , we limit our discussion to trivalent intersections ) .
this is based on the fact proved in the appendix of the previous paper , that these invariants are annihilated by the loop derivative that appears in both the diffeomorphism and the hamiltonian constraints and of the expression of the rhs .
therefore the ambient isotopic ( framing - independent ) vassiliev invariants for spin networks with trivalent intersections are annihilated by all the constraints of quantum gravity and consistently , by the right hand side of the commutator of two hamiltonian constraints .
it should be pointed out that the space is in no way trivial : as we discussed in the appendix of the previous paper , the annihilation is a detailed property of the space of vassiliev invariants , related to the decomposition of the invariants in framing independent and framing dependent components and the detailed structure of chord diagrams appearing in the coefficients .
hints that the vassiliev invariants for trivalent intersections were annihilated by the hamiltonian constraint were found for the lower invariants in terms of loops @xcite,@xcite , in the lattice @xcite .
one could mention as a more trivial example of a space with similar property the states based on diffeomorphism invariants of spin networks with no vertices , which are also trivially annihilated by all the constraints we consider .
this is an extension of results also first suggested in the context of loops @xcite .
it should be remarked that this property of the loop derivative is not true for higher valence intersections .
future extensions of the hamiltonian constraint to higher valence intersections could be tested in this habitat for non - trivial consistency .      by this space
we mean the invariants that appear in the power series expansion of the expectation value of the wilson net in a chern  simons theory in terms of the inverse coupling constant @xmath209 .
in particular , the whole series is an example of such an invariant .
each coefficient of the series and certain portions of them are also examples . in general ,
these invariants are sums and products of the independent invariants that appear at each order , both framing - independent and framing - dependent ( see previous paper ) .
these states have the important property that the loop derivative that appears in the expressions of the constraints ( which is evaluated on a path @xmath16 of infinitesimal length ) can be rewritten simply as , @xmath210 where @xmath211 is the edge of the spin net containing the point @xmath17 and where the invariant @xmath212 is another vassiliev invariant of one order less than @xmath213 .
what happens is that because of the infinitesimal length of @xmath16 one can rearrange the action of the loop derivative in terms of the original spin network @xmath52 using recoupling identities , at the expense of some additional group factors , which we reabsorb notationally in @xmath212 .
we will see that we do not need the details of the relation of @xmath213 to @xmath212 for proving the consistency . as usual , all our discussion is limited to trivalent intersections .
the relation above implies that for the vassiliev invariants , one can write , @xmath214 this can be seen by considering equation ( [ nu ] ) and recalling that @xmath215 .
in particular , we see that the action of the hamiltonian on these states can be written ( up to a constant factor which we will omit for simplicity ) as , @xmath216      on this space , the diffeomorphism constraint vanishes identically , as we discussed in the previous paper .
so the algebra of diffeomorphisms is trivially satisfied .
nontrivial commutators to be realized will be those of a diffeomorphism with a hamiltonian and that of two hamiltonians .      in the case of a diffeomorphism with a hamiltonian ,
since the wavefunctions are diffeomorphism invariant , of the two terms of the commutator one is left with the one in which the diffeomorphism acts at the left .
since the action of the hamiltonian on @xmath213 is not diffeomorphism invariant , one should recover the proper action of the diffeomorphism through such term . here
the derivations of section [ marked ] become useful .
we start from the expression ( [ honvas ] ) , and then use ( [ diffeoono ] ) to get , @xmath217 and from here the correct commutator follows immediately .
this implies that the quantum hamiltonian transforms covariantly in a correct way .
this calculation is one of the main differences of our construction with respect to the one of thiemann @xcite , since in that context one does not consider an infinitesimal generator of diffeomorphisms .
as we discussed in section [ marked ] , the action of a hamiltonian on a function with marked points as one gets after acting with a hamiltonian on a vassiliev invariant was given by expression ( [ [ ho ] ] ) , which we can combine with the action of the hamiltonian ( [ honvas ] ) to get , @xmath218 and if one performs the calculation in the reverse order , one obtains the @xmath103 operators in the reverse order .
however , these operators are multiplicative , so they commute .
therefore one has the result , @xmath219 v_n(s)=0,\ ] ] which is the expected commutation relation on states that are invariant under diffeomorphisms .
the states we have already considered , where one has a vassiliev invariant times a scalar function that depends on the position of the vertices of the spin network , are an attractive habitat where the diffeomorphism constraint does not vanish .      as discussed in section ( [ marked ] ) the diffeomorphisms on these states reduce to lie dragging of the scalar functions of the vertices ( see formula ( [ diffeoono ] ) ) .
since they correspond to a natural geometric action , the consistency of the algebra of diffeomorphisms is immediate .
we start from the action of a hamiltonian on these kinds of states , @xmath220 where we have used the fact that @xmath221 and @xmath222 commute and that the action of the hamiltonian on these states produces a vassiliev invariant of order lowered by one unit times a vertex function .
we now act with a diffeomorphism , @xmath223 and now use the fact that leibniz rule applies to the action of the diffeomorphism constraint ( stemming from the fact that it also applies to the action of the loop derivative ) , @xmath224 and we now reconstruct the hamiltonian in the first and second terms , @xmath225 and using the fact that @xmath222 and @xmath221 commute , we get , @xmath226 which we can rewrite as , @xmath227 from where we get the correct commutation relation , @xmath228 o(m,\omega ) v_n   =   h(l^a \partial_a n ) o(m,\omega ) v_{n}.\ ] ]
this calculation proceeds along the same lines as the one in the previous subsection , we essentially rewrite the action of the hamiltonian in terms of the @xmath221 operator , and note that the @xmath221 operators commute , @xmath229 and we therefore get that @xmath230=0 $ ] , which as we discussed before , is consistent with the representation in these spaces of functions of the right hand side of the quantum commutator .
we have presented a canonical quantization of the constraints of canonical general relativity .
we represented the diffeomorphism and hamiltonian constraints using two novel ingredients : the loop derivative to represent the field tensor @xmath14 and the use of spaces related to the generalization of the vassiliev invariants to spin networks ( restricted to trivalent intersections ) as wavefunctions . in terms of the latter we constructed several `` habitats '' , including spaces of functions that are not invariant under diffeomorphisms and we checked that one obtained a consistent algebra of quantum commutators of the constraints . consistency in this sense
implies that the quantization of the canonical poisson identities between the classical constraints is implemented correctly in the quantization .
we have observed that this consistency is achieved at the price of having a vanishing right hand - side for the commutator of two hamiltonians and that we can trace back this fact to the vanishing nature of the doubly covariant metric operator in these spaces of functions .
we should point out common elements and differences with the quantization presented by thiemann @xcite . in thiemann
s case the hamiltonian was implemented on the space of diffeomorphism invariant cylindrical functions of spin networks . on these spaces
one does not have a well defined notion of the field tensor @xmath14 and the functions are not loop - differentiable . in our case
, the availability of the loop derivatives allows to have at hand a `` differential calculus '' that allows for several novel constructions .
examples of them are the calculations performed in computing the constraint algebra , and the possibility of finding novel solutions to the hamiltonian constraint , based on the behavior of the vassiliev invariants under loop differentiation .
another difference involves the implementation of an infinitesimal generator of diffeomorphisms . in thiemann s original construction
one worked directly in terms of diffeomorphism invariant states and therefore one did not have an infinitesimal generator .
the construction was extended by lewandowski and marolf @xcite to `` habitats '' that are dependent on diffeomorphisms . with that extension
, thiemann s work achieves the same level of consistency that we have in this paper , in the sense that one can check non - trivial commutators of diffeomorphisms and diffeomorphisms with hamiltonians .
one still has the feature of a vanishing right hand side of the commutator of two hamiltonians @xcite . as we discussed in this paper
, this feature appears as inescapable in the context of wavefunctions only dependent on spin networks as the ones we considered here , since as we mentioned one does not have enough structures to construct a naturally defined symmetric metric tensor .
other differences arise insofar as the space of wavefunctions considered . in the case of thiemann
, one had spaces of functions with well defined inner products , which allow to discuss normalizability and study the spaces of solutions with a level of rigor that is not available currently in our approach , since we do not have an inner product on the habitats we are considering .
it might not be impossible to find a suitable inner product with the same techniques that led to the construction of the measures on the spaces of cylindrical functions , but it has not been achieved yet .
concerning solutions of the constraints , one has some available both in our approach and in thiemann s , that appear as quite distinct in their features . in thiemann s
approach , the solutions are obtained using group averaging techniques .
this leads to structures like `` tassels '' and others @xcite in which accumulation of lines and vertices take place . in our approach
, one has a number of solutions of the hamiltonian constraint ( the framing independent vassiliev invariants for trivalent intersections ) that do not depend on such structures . on the other hand , they appear as `` trivial '' solutions in the sense that they may not exist if one considers intersections of higher valences . in our approach
one may also construct solutions to the hamiltonian constraint with a cosmological constant , following similar ideas that led to the construction of states in terms of loops @xcite .
although we have not pursued this in detail yet , it appears quite plausible that these types of solutions exist , given the structure of the extra term due to the cosmological constant in the hamiltonian constraint in terms of the spin network approach @xcite .
thiemann s approach has also been studied in @xmath231 dimensions @xcite , and appears to lead to a satisfactory quantization , provided one chooses in an ad - hoc way an inner product that rules out certain infinite dimensional set of solutions . in a forthcoming paper
we will discuss the quantization of @xmath231 dimensional gravity using an approach that has elements in common with the one we pursue here , in particular the requirement of loop differentiability of the states .
we will see that this requirement limits us ( at least for low valence intersections ) to the correct solution space in a natural way .
the reader might find unsatisfactory that the right hand side of the commutator of two hamiltonians vanishes .
even more unsatisfactory may appear the fact that this is due to the vanishing of the doubly covariant metric in this approach to quantum gravity . in this subsection
we will address this and other issues .
the first question that one may raise is if this is not just a pathology that is introduced by our limitation to trivalent intersections for reasons of calculational convenience .
after all , one knows that this subspace in in some sense degenerate since the volume operator vanishes identically .
unfortunately , it appears unlikely that extending our results to higher valent intersections ( to which we do not see any technical obstruction , apart from greater calculational complexity ) will change things insofar as the commutator of two hamiltonians .
the action of the hamiltonian constraint on higher valent intersections is more complicated largely because the volume operator and the action of the loop derivative can not be reduced to a simple group - dependent prefactor times the original state , but will in general involve a linear combination of states , with a non - trivial group - dependent matrix of coefficients
. however , the dependence of the hamiltonian on the smearing function will remain as a multiplicative one .
therefore it appears that the commutator of two hamiltonians will again vanish .
moreover , the reasons we gave for the vanishing of the right hand side of the commutator are independent of the valences of the intersections of the spin networks .
the expression for the doubly - covariant metric will be more complicated , but will still include a prefactor with the same dependence on the regulators as in the trivalent case and it will vanish .
one therefore expects consistency at the same level we achieved for the trivalent case .
can one consider satisfactory a theory with a vanishing doubly contravariant metric ( and a divergent doubly covariant metric ) ?
the answer to this will only appear in a definitive way when one constructs physical predictions from the theory , which at the moment are lacking .
one can get partial indications that the pathology might not be as severe as it appears at first sight from the fact that one can define reasonable quantities like areas , volumes and lengths in this framework in spite of having the ill defined metric operators .
an attractive feature of the vassiliev states is that other diffeomorphism invariant operators can be defined as well . for instance , if one considers the integral @xmath232 , given that the chern
simons states are such that the densitized triads are proportional to the magnetic vector constructed from the curvature , one finds that quantum representation of this integral is identical to the volume operator .
an important observation concerning how satisfatory is a theory with a vanishing metric operator like the one we propose is that in @xmath231-dimensional gravity similar pathologies appear , yet the correct physical theory is recovered . if one pursues a quantization similar to thiemann s dimensional gravity along the same lines as the one proposed in this paper , and the observations we make in this subsection apply to it as well .
we will also expand the discussion on the physical observables and their connection with the kinematical metric we present in this subsection . ] in @xmath231=dimensions ( this is discussed in detail in @xcite ) one notices that the hamiltonians also commute .
moreover , the geometrical arguments leading to the vanishing of the @xmath202 operator still hold in this context : there is no quantity one can build out of loop states that will yield a doubly contravariant symmetric tensor .
yet , it was shown in @xcite that one can recover the correct physical theory . how can these apparently contradictory elements be reconciled ? it has to do with the nature of the kinematical calculations we are performing . to check the `` off shell '' contraint algrebra one
is required to operate with a kinematial space of wavefunctions that are not annihilated by the constraints .
when one recovers `` physics '' one should really do it with states that are annihilated by the constraint .
this can be carried out in detail in the @xmath231 dimensional case .
suppose one wishes to ask questions about the `` metric '' of the theory .
the first observation is that on physical states one can not define a metric tensor since the latter states are diffeomorphism invariant and the metric tensor is not .
one could gauge fix a metric tensor , say by asking questions about the value of its components in a fixed coordinate system .
these questions can be answered : the components of the metric in a fixed coordinate system can be related to the values of the invariant operators @xmath233 and @xmath234 in such a coordinate system ( and therefore in general , since the latter are coordinate invariant ) .
one therefore has well defined operators associated with the metric that can be evaluated in the physical space of states .
but such operators have little to do with the @xmath202 we introduced in the kinematical space .
it is not obvious at all that the two operators will be related in any way since the former are non - vanishing whereas the latter vanishes identically .
there is a certain disquieting element in the last observation , since it seems to imply that one should be careful before drawing conclusions from calculations at a kinematical level .
this in particular , implies all calculations involving the constraint algebra . when we set out to work on the current paper , our expectation was that reproducing the contraint algebra would be a strong test for our quantization procedure , that perhaps would rule out all but a few of the possible theories .
what we have learned is that this is probably not the correct view on the issue : one can obtain consistency in many ways at a kinematical level ( in particular with @xmath235 ) .
one can construct many consistent theories of quantum gravity . and still , what the @xmath231 dimensional example shows is that the `` physics '' of all these theories is really deeply buried in the states that solve the constraints . in a sense this is good ,
since the theory we are proposing here differs significantly at this level from the one proposed by thiemann and therefore further adds to our options for trying to match experimental results .
the observation that kinematical calculations may have quite nontrivial connections with `` real physics '' poses problems for computations that try to obtain heuristic physical insights from looking at properties of the kinematical states ( the `` weave '' approach to the semiclassical theory ) .
an example of this are our own results on gamma - ray - burst light dispersion @xcite . from the @xmath231 example
we see that quantum fluctuations in the `` physical metric '' ( the ones one would expect to influence the propagation of matter fields ) may not have a direct connection with fluctuations of the kinematical metric . yet it is the latter that are used in the concrete computations of @xcite .
concerning the latter point , it is worthwhile mentioning explicitly that the theory presented in this paper can be coupled to matter and have well defined expressions for the hamiltonians of matter fields , which also involve the metric in a non - trivial way ( the coupling of the theory to matter can , for instance , be achieved using the same construction as thiemann @xcite ; or one could consider alternative settings in which the loop derivative is also used in the representation of certain matter fields ) .
it is somewhat unfortunate that the remarkable results on black hole entropy that have been achieved within this context @xcite do not involve in a detailed way the action of the hamiltonian constraint in the `` bulk '' of the spacetime to be used as guideline for constructing the constraints , and therefore can not distinguish at the moment our proposal from thiemann s .
the vanishing of the doubly contravariant metric tensor , apart from appearing as a `` robust '' feature based only on elementary notions of covariance of the elements involved in constructing the spin network states , is the way in which our approach bypasses the `` hermiticity problem '' of the canonical quantization . as we mentioned , if one were to demand that the diffeomorphism and hamiltonian constraints and the doubly covariant metric be hermitian operators ( since they correspond to real classical quantities )
, the commutator of two hamiltonians ( which is hermitian ) can not be simply equal to the product of a metric times a diffeomorphism ( which is not hermitian ) .
one could fix the hermiticity of the right hand side by `` symmetrizing '' the operator , but then one could face an anomaly problem , since the diffeomorphism constraint would not act to the right of the product .
this difficulty is bypassed ( in the context of non - diffeomorphism invariant states ) if the metric vanishes .
it should be noticed that it is not obvious that one should promote the constraints to hermitian operators , for instance , see @xcite for counterexamples in the context of the quantization of non - unimodular gauge groups .
related with the issue of the lack of guidelines to construct the theory is the problem of uniqueness .
we already have a manifest non - uniqueness in that the theory presented here and the theory introduced by thiemann appear as both satisfactory yet distinct .
the ambiguity is worse than this . even if one stays within one general approach , for instance considering the space of vassiliev invariants as `` arena '' for quantization ,
there are many ways in which one could implement the hamiltonian constraint .
we have not analyzed them in any detail , but some salient features of them are worthwhile mentioning .
the ambiguities arise in the various limits involved in constructing the hamiltonian .
one of them is associated with the specific action of the loop derivative .
the loop derivative is dependent on a path . in this paper
we have taken such a path to coincide with one of the lines of the spin network .
we then acted with the loop derivative and collapsed its action using recoupling identities .
the final action of the hamiltonian was therefore `` ultra - local '' in the sense that it returned a wavefunction with the same vertex structure times a vertex dependent prefactor ( for valences higher than three one gets a linear combination with different spin weights for the intertwiners yet the topology of the vertex is unchanged ) .
this may raise concerns that `` super - selections '' could appear in the sense that it could be easy to construct operators that commute with the hamiltonian .
one could define a different hamiltonian in a straightforward way , by assuming a different topology for the path associated with the loop derivative .
for instance , one could assume that the path starts at the vertex , advances along one of the edges of the spin network and then crosses towards another edge of the spin network along one of the edges of the tetrahedra introduced in the discretization .
at the end of the path , the loop derivative acts .
one then is left with an action of the hamiltonian resembling the one introduced by thiemann : the constraint produces a vertex dependent prefactor , but it also alters the structure of the spin network by `` adding a line '' at the vertex .
a different proposal would be to have a path that starts at the vertex , advances along one of the edges of the spin network , crosses towards another edge along the tetrahedron and then continues to cross towards another edge of the spin network , as shown in figure [ fifi ] .    in the first of these proposals ( as in thiemann s ) , the hamiltonian adds vertices called `` exceptional '' in the sense that they are planar vertices .
the hamiltonian vanishes identically on such vertices , because the volume vanishes on planar vertices , even if they are four - valent . in the second proposal ,
the intermediate vertex that is created is non - planar .
this might raise hopes that the commutator of two hamiltonians may be less trivial , since the second hamiltonian acts at the newly created vertex , but more careful analyses seem to indicate that the hamiltonians still commute @xcite .
if it they were not to commute one would still be faced with how to reproduce a similar vertex structure with the right hand side operator in a natural way . summarizing this portion of the discussion
: it appears that there is a non - trivial , possibly infinite , amount of ambiguity in the definition of the theory , that is not significantly constrained by the imposition of the correct quantum commutator algebra .
related to the latter point is the fact that due to the vanishing of the metric , the consistency check provided by the commutator of two hamiltonians is less detailed than if one had to prove the equality of non - vanishing operators .
we already see that we have two distinct and apparently consistent , quantizations of the constraints ( with several possible variants of each ) .
it will require further study to determine if one quantization is `` better '' in the sense of reproducing expected results than the other .
one possibility would be to consider the commutation of the hamiltonian with various operators , and study if inconsistencies appear .
this procedure , however , might be limited by the fact that most commutators have non - trivial right hand sides which will inevitably involve quite a bit of ambiguity at the time of their quantization .
another point that might be raised is that these approaches appear confined to four dimensions and to the einstein - hilbert action , and therefore one may have little hopes of making contact with other approaches , as those based on string theories .
these approaches not only may present avenues to understand quantum gravity but they also have the attractive feature of unifying gravity with other interactions , a goal some may consider desirable in itself .
it should be pointed out however , that progress is being made @xcite towards describing @xmath236 dimensional general relativity in terms of connections .
although constructions like the ones we discussed here have not been pursued in detail in this context , they appear as plausible .
to conclude , we have at the moment canonical quantizations of general relativity ( possibly coupled to matter ) that appear as mathematically consistent at the kinematical level at which they have been studied .
this level of consistency had never been achieved before in other approaches .
further exploration of the consequences of the quantization will be needed to determine if any of them are physically satisfactory theories of the quantum gravitational field . in particular the exploration of
the space of states that solve the constraints and how physical quantities can be evaluated on them appears as a natural next step in the quantization program .
we wish to thank abhay ashtekar , laurent freidel , john klauder , karel kucha and thomas thiemann for comments and discussions .
this work was supported in part by the national science foundation under grants nsf - phy-9423950 , nsf - int-9811610 , nsf - phy-9407194 , research funds of the pennsylvania state university , the eberly family research fund at psu .
jp acknowledges support of the alfred p. sloan and john simon guggenheim foundations .
we acknowledge support of pedeciba ( uruguay ) .
rg and jp wish to thank the institute for theoretical physics of the university of california at santa barbara and cdb , rg and jg the center for gravitational physics and geometry at penn state for hospitality during the completion of this work .    c. di bartolo , r. gambini , j. griego , j. pullin `` canonical quantum gravity in the vassiliev invariants arena : i. kinematical structure '' , companion paper . c. rovelli , l. smolin , nucl .
phys . * b442 * , 593 ( 1995 ) .
j. barbero , phys . rev .
* d51 * , 5507 ( 1995 ) .
t. thiemann , class .
grav . * 15 * , 839 ( 1998 ) .
r. borissov , r. de pietri , c. rovelli , class .
* 14 * , 2793 ( 1997 ) .
t. thiemann , j. math . phys . * 39 * , 3347 ( 1988 ) .
r. gambini , phys . lett . *
b255 * , 180 ( 1991 ) .
r. gambini , a. garat , j. pullin , int . j. mod
d4 * , 589 ( 1995 ) .
r. gambini , j. pullin `` loops , knots , gauge theories and quantum gravity '' cambridge university press ( 1996 ) .
t. thiemann , class .
grav . * 15 * , 1207 ( 1998 ) .
r. gambini , j. lewandowski , d. marolf , j. pullin , int .
* d7 * 97 ( 1997 ) . b. brgmann , r. gambini , j. pullin , phys .
lett . * 68 * 431 ( 1992 ) .
j. griego , phys .
* d53 * 6966 ( 1996 ) .
r. gambini , j. pullin , phys . rev . *
d54 * 5935 ( 1996 ) . c. rovelli , l. smolin , phys .
lett . * 61 * 1155 ( 1988 ) .
j. lewandowski , d. marolf , int . j. mod
d7 * , 299 ( 1998 ) .
j. lewandowski , t. thiemann , class .
* 16 * 2299 ( 1999 ) .
b. brgmann , r. gambini , j. pullin , nucl . phys . *
b385 * , 587 ( 1992 ) .
r. gambini , j. griego , j. pullin , nucl . phys . *
b534 * , 675 - 696 ( 1998 ) .
r. gambini , j. pullin , phys . rev . *
d59 * , 124021 ( 1999 )
. t. thiemann , class .
grav . * 15 * , 1249 ( 1998 ) .
d. marolf , d. giulini , class .
* 16 * 2479 ( 1999 ) .
see for instance a. ashtekar , j. baez , a. corichi , k. krasnov phys .
lett . * 80 * , 904 ( 1998 ) , and references therein .
l. freidel , k. krasnov , r. puzio , hep - th/9901069 .