This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Black hole entropy from the SU(2)SU(2)-invariant formulation of Type I isolated horizons

Jonathan Engle 1,2    Karim Noui 3    Alejandro Perez1    Daniele Pranzetti1 1Centre de Physique Théorique111Unité Mixte de Recherche (UMR 6207) du CNRS et des Universités Aix-Marseille I, Aix-Marseille II, et du Sud Toulon-Var; laboratoire afilié à la FRUMAM (FR 2291), Campus de Luminy, 13288 Marseille, France. 2Institut für Theoretische Physik III, Universität Erlangen-Nürnberg, Staudtstraße 7, 91058 Erlangen, Germany. 3 Laboratoire de Mathématiques et Physique Théorique222Fédération Denis Poisson Orléans-Tours, CNRS/UMR 6083, 37200 Tours, France.
(August 10, 2025)
Abstract

A detailed analysis of the spherically symmetric isolated horizon system is performed in terms of the connection formulation of general relativity. The system is shown to admit a manifestly SU(2)SU(2) invariant formulation where the (effective) horizon degrees of freedom are described by an SU(2)SU(2) Chern-Simons theory. This leads to a more transparent description of the quantum theory in the context of loop quantum gravity and modifications of the form of the horizon entropy.

pacs:
04.60.-m, 04.60.Pp, 04.20.Fy, 11.15.Yc

I Introduction

Black holes are intriguing solutions of classical general relativity describing important aspects of the physics of gravitational collapse. Their existence in our nearby universe is by now supported by a great amount of observational evidence observ . When isolated, these systems are remarkably simple for late and distant observers: once the initial very dynamical phase of collapse is passed the system is expected to settle down to a stationary situation completely described (as implied by the famous results by Carter, Israel, and Hawking wald ) by the three extensive parameters (mass MM, angular momentum JJ, electric charge QQ) of the Kerr-Newman family kerrnew .

However, the great simplicity of the final stage of an isolated gravitational collapse for late and distant observers is in sharp contrast with the very dynamical nature of the physics seen by in-falling observers which depends on all the details of the collapsing matter. Moreover, this dynamics cannot be consistently described for late times (as measured by the infalling observers) using general relativity due to the unavoidable development, within the classical framework, of unphysical pathologies of the gravitational field. Concretely, the celebrated singularity theorems of Hawking and Penrose hawking imply the breakdown of predictability of general relativity in the black hole interior. Dimensional arguments imply that quantum effects cannot be neglected near the classical singularities. Understanding of physics in this extreme regime requires a quantum theory of gravity. Black holes (BH) provide, in this precise sense, the most tantalizing theoretical evidence for the need of a more fundamental (quantum) description of the gravitational field.

Extra motivation for the quantum description of gravitational collapse comes from the physics of black holes available to observers outside the horizon. As for the interior physics, the main piece of evidence comes from the classical theory itself which implies an (at first only) apparent relationship between the properties of idealized black hole systems and those of thermodynamical systems. On the one hand, black hole horizons satisfy the very general Hawking area theorem (the so-called second law) stating that the black hole horizon area aHa_{\scriptscriptstyle H} can only increase, namely

δaH0.\delta a_{\scriptscriptstyle H}\geq 0. (1)

On the other hand, the uniqueness of the Kerr-Newman family, as the final (stationary) stage of the gravitational collapse of an isolated gravitational system, can be used to prove the first and zeroth laws: under external perturbation the initially stationary state of a black hole can change but the final stationary state will be described by another Kerr-Newman solution whose parameters readjust according to the first law

δM=κH8πGδaH+ΦHδQ+ΩHδJ,\delta M=\frac{\kappa_{\scriptscriptstyle H}}{8\pi G}\delta a_{\scriptscriptstyle H}+\Phi_{\scriptscriptstyle H}\,\delta Q+\Omega_{\scriptscriptstyle H}\,\delta J, (2)

where κH\kappa_{\scriptscriptstyle H} is the surface gravity, ΦH\Phi_{\scriptscriptstyle H} is the electrostatic potential at the horizon, and ΩH\Omega_{\scriptscriptstyle H} the angular velocity of the horizon. There is also the zeroth law stating the uniformity of the surface gravity κH\kappa_{\scriptscriptstyle H} on the event horizon of stationary black holes, and finally the third law precluding the possibility of reaching an extremal black hole (for which κH=0\kappa_{\scriptscriptstyle H}=0) by means of any physical process333The third law can only be motivated by a series of examples. Extra motivations comes from the validity of the cosmic censorship conjecture.. The validity of these classical laws motivated Bekenstein to put forward the idea that black holes may behave as thermodynamical systems with an entropy S=αa/p2S=\alpha a/\ell_{p}^{2} and a temperature kT=κH/(8πα)kT=\hbar\kappa_{\scriptscriptstyle H}/(8\pi\alpha) where α\alpha is a dimensionless constant and the dimensionality of the quantities involved require the introduction of \hbar leading in turn to the appearance of the Planck length p\ell_{p}, even though in his first paper beke Bekenstein states “that one should not regard TT as the temperature of the black hole; such identification can lead to all sorts of paradoxes, and is thus not useful”. The key point is that the need of \hbar required by the dimensional analysis involved in the argument called for the investigation of black hole systems from a quantum perspective. In fact, not long after, the semiclassical calculations of Hawking Hawking:1974sw —that studied particle creation in a quantum test field (representing quantum matter and quantum gravitational perturbations) on the space-time background of the gravitational collapse of an isolated system described for late times by a stationary black hole—showed that once black holes have settled to their stationary (classically) final states, they continue to radiate as perfect black bodies at temperature kT=κH/(2π)kT=\kappa_{\scriptscriptstyle H}\hbar/(2\pi). Thus, on the one hand, this confirmed that black holes are indeed thermal objects that radiate at a the given temperature and whose entropy is given by S=a/(4p2)S=a/(4\ell^{2}_{p}), while, on the other hand, this raised a wide range of new questions whose proper answer requires a quantum treatment of the gravitational degrees of freedom.

Among the simplest questions is the issue of the statistical origin of black hole entropy. In other words, what is the nature of the the large amount of micro-states responsible for black hole entropy. This simple question cannot be addressed using semiclassical arguments of the kind leading to Hawking radiation and requires a more fundamental description. In this way, the computation of black hole entropy from basic principles became an important test for any candidate quantum theory of gravity. In string theory it has been computed using dualities and no-normalization theorems valid for extremal black holes string . There are also calculations based on the effective description of near horizon quantum degrees of freedom in terms of effective 22-dimensional conformal theories carlip . In loop quantum gravity the first computations (valid for physical black holes) were based on general considerations and the fact that the area spectrum in the theory is discrete bhe0 . The calculation was later refined by quantizing a sector of the phase space of general relativity containing a horizon in ‘equilibrium’ with the external matter and gravitational degrees of freedom bhe1 . In all cases agreement with the Bekenstein-Hawking formula is obtained with logarithmic corrections in a/p2a/\ell^{2}_{p}.

In this work we concentrate and further develop the theory of isolated horizons in the context of loop quantum gravity. Recently, we have proposed a new computation of BH entropy in loop quantum gravity (LQG) that avoids the internal gauge-fixing used in prior works nous and makes the underlying structure more transparent. We show, in particular, that the degrees of freedom of Type I isolated horizons can be encoded (along the lines of the standard treatment) in an SU(2)SU(2) boundary connection. The results of this work clarify the relationship between the theory of isolated horizons and SU(2) Chern-Simons theory first explored in kiril-lee , and makes the relationship with the usual treatment of degrees of freedom in loop quantum gravity clear-cut. In the present work, we provide a full detailed derivation of the result of our recent work and discuss several important issues that were only briefly mentioned then.

An important point should be emphasized concerning the logarithmic corrections mentioned above. The logarithmic corrections to the Bekenstein-Hawking area formula for black hole entropy in the loop quantum gravity literature were thought to be of the (universal) form ΔS=1/2log(aH/p2)\Delta S=-1/2\log(a_{H}/\ell^{2}_{p}) amit . In majundar Kaul and Majumdar pointed out that, due to the necessary SU(2)SU(2) gauge symmetry of the isolated horizon system, the counting should be modified leading to corrections of the form ΔS=3/2log(aH/p2)\Delta S=-3/2\log(a_{H}/\ell^{2}_{p}). This suggestion is particularly interesting because it would eliminate the apparent tension with other approaches to entropy calculation. In particular their result is in complete agreement with the seemingly very general treatment (which includes the string theory calculations) proposed by Carlip carlip-log . Our analysis confirms Kaul and Majumdar’s proposal and eliminates in this way the apparent discrepancy between different approaches.

The article is organized as follows. In the following section we review the formal definition of isolated horizons. In Section III we state the main equations implied by the isolated horizon boundary conditions for fields at a spherically symmetric isolated horizon. In Section IV we prove a series of propositions that imply the main classical part of our results: we derive the form of the conserved presymplectic structure of spherically symmetric isolated horizons, and we show that degrees of freedom at the horizon are described by an SU(2)SU(2) Chern-Simons presymplectic structure. In Section VI we briefly review the derivation of the zeroth and first law of isolated horizons. In Section V we study the gauge symmetries of the Type I isolated Horizon and explicitly compute the constraint algebra. In Section VII we review the quantization of the spherically symmetric isolated horizon phase space and present the basic formulas necessary for the counting of states that leads to the entropy. We close with a discussion of our results in Section VIII. The appendix contains an analysis of Type I isolated horizons from a concrete (and intuitive) perspective that makes use of the properties of stationary spherically symmetric black holes in general relativity.

II Definition of isolated horizons

The standard definition of a BH as a spacetime region from which no information can reach idealized observers at (future null) infinity is a global definition. This notion of BH requires a complete knowledge of a spacetime geometry and is therefore not suitable for describing local physics. The physically relevant definition used, for instance, when one claims there is a black hole in the center of the galaxy, must be local. One such local definition was introduced in ack ; better ; ih_prl with the name of isolated horizons (IH). Here we present this definition according to ih_prl ; afk ; abl2002 ; abl2001 . This discussion will also serve to fix our notation. In the definition of an isolated horizon below, we allow general matter, subject only to conditions that we explicitly state.

Definition: The internal boundary Δ\Delta of a history (M,gab)({\mathfs{M}},g_{ab}) will be called an isolated horizon provided the following conditions hold:

  1. i)

    Manifold conditions: Δ\Delta is topologically S2×RS^{2}\times R, foliated by a (preferred) family of 2-spheres SS and equipped with an equivalence class [a][\ell^{a}] of transversal future pointing vector fields whose flow preserves the foliation, where a\ell^{a} is equivalent to a\ell^{\prime a} if a=ca\ell^{a}=c\ell^{\prime a} for some positive real number cc.

  2. ii)

    Dynamical conditions: All field equations hold at Δ\Delta.

  3. iii)

    Matter conditions: On Δ\Delta the stress-energy tensor TabT_{ab} of matter is such that Tabb-T^{a}{}_{b}\ell^{b} is causal and future directed.

  4. iv)

    Conditions on the metric gg determined by ee, and on its levi-Civita derivative operator \nabla: (iv.a) The expansion of a\ell^{a} within Δ\Delta is zero. This, together with the energy condition (iii) and the Raychaudhuri equation at Δ\Delta, ensures that a\ell^{a} is additionally shear-free. This in turn implies that the Levi-Civita derivative operator \nabla naturally determines a derivative operator DaD_{a} intrinsic to Δ\Delta via XaDaYb:=XaaYbX^{a}D_{a}Y^{b}:=X^{a}\nabla_{a}Y^{b}, Xa,YaX^{a},Y^{a} tangent to Δ\Delta. We then impose (iv.b) [L,D]=0[{\mathfs{L}}_{\ell},D]=0.

  5. v)

    Restriction to ‘good cuts.’ One can show furthermore that Dab=ωabD_{a}\ell^{b}=\omega_{a}\ell^{b} for some ωa\omega_{a} intrinsic to Δ\Delta. A 2-sphere cross-section SS of Δ\Delta is called a ‘good cut’ if the pull-back of ωa\omega_{a} to SS is divergence free with respect to the pull-back of gabg_{ab} to SS. As shown in abl2002 , every horizon satisfying (i)-(iv) above possesses at least one foliation into ‘good cuts’; this foliation is furthermore generically unique. We require that the fixed foliation coincide with a foliation into ‘good cuts.’

Let us discuss the physical meaning of these conditions. The first two conditions are rather weak. The third condition is satisfied by all matter fields normally used in general relativity. The fifth condition is a partial gauge fixing of diffeomorphisms in the ‘time’ direction. The main condition is therefore the fourth condition. (iv.a) requires that a\ell^{a} be expansion-free. This is equivalent to asking that the area 2-form of the 2-sphere cross-sections of Δ\Delta be constant along generators [a][\ell^{a}]. This combined with the matter condition (iii) and the Raychaudhuri equation implies that in fact the entire pull back qabq_{ab} of the metric to the horizon is Lie dragged by a\ell^{a}. Condition (iv.b) further stipulates that the derivative operator DaD_{a} be Lie dragged by a\ell^{a}. This implies, among other things, an analogue of the zeroeth law of black hole mechanics: conditions (i) and (iii) imply that a\ell^{a} is geodesic — bbaa\ell^{b}\nabla_{b}\ell_{a}\propto\ell_{a}. The proportionality constant is called the surface gravity, and condition (iv.b) ensures that it is constant on the horizon for any given a[a]\ell^{a}\in[\ell^{a}]. Furthermore, if we had not fixed [a][\ell^{a}], but only required that an [a][\ell^{a}] exist such that the isolated horizon boundary conditions hold, then condition (iv.b) would ensure that this a\ell^{a} is generically unique abl2002 . From the above discussion, one sees that the geometrical structures on Δ\Delta that are time-independent are precisely the pull-back qabq_{ab} of the metric to Δ\Delta, and the derivative operator DD. In fact, the main conditions (iv.a) and (iv.b) are equivalent to requiring Lqab=0{\mathfs{L}}_{\ell}q_{ab}=0 and [L,D]=0[{\mathfs{L}}_{\ell},D]=0. For this reason it is natural to define (qab,D)(q_{ab},D) as the horizon geometry.

Let us summarize. Isolated horizons are null surfaces, foliated by a family of marginally trapped 2-spheres such that certain geometric structures intrinsic to Δ\Delta are time independent. The presence of trapped surfaces motivates the term ‘horizon’ while the fact that they are marginally trapped — i.e., that the expansion of a\ell^{a} vanishes — accounts for the adjective ‘isolated’. The definition extracts from the definition of Killing horizon just that ‘minimum’ of conditions necessary for analogues of the laws of black hole mechanics to hold. Boundary conditions refer only to behavior of fields at Δ\Delta and the general spirit is very similar to the way one formulates boundary conditions at null infinity.

Remarks:

  1. 1.

    All the boundary conditions are satisfied by stationary black holes in the Einstein-Maxwell-dilaton theory possibly with cosmological constant. Note however that, in the non-stationary context, there still exist physically interesting black holes satisfying our conditions: one can solve for all our conditions and show that the resulting 4-metric need not be stationary on Δ\Delta lew2000 .

  2. 2.

    In the choice of boundary conditions, we have tried to strike the usual balance: On the one hand the conditions are strong enough to enable one to prove interesting results (e.g., a well-defined action principle, a Hamiltonian framework, and a realization of black hole mechanics) and, on the other hand, they are weak enough to allow a large class of examples. As we already remarked, the standard black holes in the Einstein-Maxwell-dilatonic systems satisfy these conditions. More importantly, starting with the standard stationary black holes, and using known existence theorems one can specify procedures to construct new solutions to field equations which admit isolated horizons as well as radiation at null infinity lew2000 . These examples, already show that, while the standard stationary solutions have only a finite parameter freedom, the space of solutions admitting isolated horizons is infinite dimensional. Thus, in the Hamiltonian picture, even the reduced phase-space is infinite dimensional; the conditions thus admit a very large class of examples.

  3. 3.

    Nevertheless, space-times admitting isolated horizon are very special among generic members of the full phase space of general relativity. The reason is apparent in the context of the characteristic formulation of general relativity where initial data are given on a set (pairs) of null surfaces with non trivial domain of dependence. Let us take an isolated horizon as one of the surfaces together with a transversal null surface according to the diagram shown in Figure 1. Even when the data on the isolated horizon may be infinite dimensional (for Type II and II isolated horizons, see below), in all cases no transversing radiation data is allowed by the IH boundary condition. Roughly speaking the isolated horizon boundary condition reduces to one half the number of local degrees of freedom.

  4. 4.

    Notice that the above definition is completely geometrical and does not make reference to the tetrad formulation. There is no reference to any internal gauge symmetry. In what follows we will deal with general relativity in the first order formulation which will introduce, by the choice of variables, an internal gauge group corresponding to local SL(2,)SL(2,\mathbb{C}) transformations (in the case of Ashtekar variables) or SU(2)SU(2) transformations (in the case of real Ashtekar-Barbero variables). It should be clear from the purely geometric nature of the above definition that the IH boundary condition cannot break by any means these internal symmetries.

Isolated horizon classification according to their symmetry groups

Next, let us examine symmetry groups of isolated horizons. A symmetry of (Δ,q,D,[a])(\Delta,q,D,[\ell^{a}]) is a diffeomorphism on Δ\Delta which preserves the horizon geometry (q,D)(q,D) and at most rescales elements of [a][\ell^{a}] by a positive constant. It is clear that diffeomorphisms generated by a\ell^{a} are symmetries. So, the symmetry group GΔG_{\Delta} is at least 1-dimensional. In fact, there are only three possibilities for GΔG_{\Delta}:

  1. (a)

    Type I: the isolated horizon geometry is spherical; in this case, GΔG_{\Delta} is four dimensional (SO(3)SO(3) rotations plus rescaling-translations444In a coordinate system where a=(/v)a\ell^{a}=(\partial/\partial v)^{a} the rescaling-translation corresponds to the affine map vcv+bv\to cv+b with c,bc,b\in\mathbb{R} constants. along \ell);

  2. (b)

    Type II: the isolated horizon geometry is axi-symmetric; in this case, GΔG_{\Delta} is two dimensional (rotations round symmetry axis plus rescaling-translations along \ell);

  3. (c)

    Type III: the diffeomorphisms generated by a\ell^{a} are the only symmetries; GΔG_{\Delta} is one dimensional.

Note that these symmetries refer only to the horizon geometry. The full space-time metric need not admit any isometries even in a neighborhood of the horizon. In this paper, as in the classic works bhe1 ; ack , we restrict ourselves to the Type I case. Although a revision would be necessary in light of the results of our present work, the quantization and entropy calculation in the context of Type II and Type III isolated horizons has been considered in jon .

    Refer to captionI+I\begin{array}[]{c}\includegraphics[height=113.81102pt]{IH.eps}\end{array}\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\begin{array}[]{c}\\ \\ \\ \\ \\ \,{\mathfs{I}}^{\scriptscriptstyle+}\\ \\ \\ \\ \\ {\mathfs{I}}^{\scriptscriptstyle-}\end{array}

Figure 1: The characteristic data for a (vacuum) spherically symmetric isolated horizon corresponds to Reissner-Nordstrom data on Δ\Delta, and free radiation data on the transversal null surface with suitable fall-off conditions. For each mass, charge, and radiation data in the transverse null surface there is a unique solution of Einstein-Maxwell equations locally in a portion of the past domain of dependence of the null surfaces. This defines the phase space of Type I isolated horizons in Einstein-Maxwell theory. The picture shows two Cauchy surfaces M1M_{1} and M2M_{2} “meeting” at space-like infinity i0i_{0}. A portion of I+{\mathfs{I}}^{+} and I{\mathfs{I}}^{-} are shown; however, no reference to future time-like infinity i+i^{+} is made as the isolated horizon need not to coincide with the black hole event horizon.

III Some extra details for Type I isolated horizons

In this section we first list the main equations satisfied by fields at an isolated horizon of Type I. The equations presented here can be directly derived from the IH boundary conditions implied by the definition of Type I isolated horizons given above. Most of the equations presented here can be found in ack . For completeness we prove these equations at the end of this section. As we shall see in Subsection III.2, some of the coefficients entering the form of these equations depend on the representative chosen among the equivalence class of null generators [][\ell]. Throughout this paper we shall fix an null generator []\ell\in[\ell] by the requirement that the surface gravity ω=κ\ell{\lrcorner}\omega=\kappa matches the one corresponding to the stationary black hole with the same macroscopic parameters as the Type I isolated horizon under consideration. This choice makes the first law of IH take the form of the usual first law of stationary black holes (see Section VI).

III.1 The main equations

When written in connection variables, the isolated horizon boundary condition implies the following relationship between the curvature of the Ashtekar connection A+i=Γi+iKiA^{i}_{\scriptscriptstyle+}=\Gamma^{i}+iK^{i} at the horizon and the 22-form Σi=ϵjkiejek\Sigma^{i}=\epsilon^{i}_{\ jk}e^{j}\wedge e^{k} (in the time gauge)

Fabi(A+)=2πaHΣabi,\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{F_{ab}}^{i}(A^{\scriptscriptstyle+})=-\frac{2\pi}{a_{\scriptscriptstyle H}}\,\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\Sigma_{ab}}^{i}, (3)

where aHa_{\scriptscriptstyle H} is the area of the IH, the double arrows denote the pull-back to H=ΔMH=\Delta\cap M with MM a Cauchy surface with normal τa=(a+na)/2\tau^{a}=(\ell^{a}+n^{a})/\sqrt{2} at HH, and nan^{a} null and normalized according to n=1n\cdot\ell=-1. Notice that the imaginary part of the previous equation implies that

dΓKi=0\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{d_{\Gamma}K}^{i}=0 (4)

Another important equation is

ϵjkiKjKk=2πaHΣi.\epsilon^{i}_{\ jk}\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{K^{j}}\wedge\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{K^{k}}=\frac{2\pi}{a_{\scriptscriptstyle H}}\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\Sigma}^{i}. (5)

The previous equations follow from equations (3.12) and (B.7) of reference ack . Nevertheless, they also follow from the abstract definition given in the introduction. From the previous equations, only equation (5) is not explicitly proven from the definition of IH in the literature. Therefore, we give here an explicit prove at the end of this section. For concreteness, as we think it is helpful for some readers to have a concrete less abstract treatment, another derivation using directly the Schwarzschild geometry is given in Appendix A. The previous equations imply in turn that

Fabi(Aβ)=π(1β2)aHΣabi,\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{F_{ab}}^{i}(A^{\scriptscriptstyle\beta})=-\frac{\pi(1-\beta^{2})}{a_{\scriptscriptstyle H}}\,\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\Sigma_{ab}}^{i}, (6)

where Aβi=Γi+βKiA^{i}_{\scriptscriptstyle\beta}=\Gamma^{i}+\beta K^{i} is the Ashtekar-Barbero connection 555In our convention the so(3)3so(3)\to\mathbb{R}^{3} isomorphism is defined by λi=12ϵjkiλjk\lambda^{i}=-\frac{1}{2}\epsilon^{i}_{\ jk}\lambda^{jk} which implies that Fi=dAi+12ϵjkiAjAkF^{i}=dA^{i}+\frac{1}{2}\epsilon^{i}_{\ jk}A^{j}\wedge A^{k} and dAλi=dλi+ϵjkiAjλkd_{A}\lambda^{i}=d\lambda^{i}+\epsilon^{i}_{\ jk}A^{j}\wedge\lambda^{k}. .

III.2 Proof of the main equations

In this subsection we use the definition of isolated horizons provided in the previous section to prove some of the equations stated above. We will often work in a special gauge where the tetrad (eI)(e^{I}) is such that e1e^{1} is normal to HH and e2e^{2} and e3e^{3} are tangent to HH. This choice is only made for convenience, as the equations presented in the previous section are all gauge covariant, their validity in one frame implies their validity in all frames.

Lemma 1: In the gauge where the tetrad is chosen so that a=21/2(e0a+e1a)\ell^{a}=2^{-1/2}(e^{a}_{0}+e^{a}_{1}) (which can be completed to a null tetrad na=21/2(e0ae1a)n^{a}=2^{-1/2}(e^{a}_{0}-e^{a}_{1}), and ma=21/2(e2a+ie3a)m^{a}=2^{-1/2}(e^{a}_{2}+ie^{a}_{3})), the shear-free and vanishing expansion (condition (iv.aiv.a) in the definition of IH) imply

ω21=ω20andω31=ω30.\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\omega}^{21}=\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\omega}^{20}\ \ \ {\rm and}\ \ \ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\omega}^{31}=\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\omega}^{30}. (7)

Proof: The expansion ρ\rho and shear σ\sigma of the null congruence of generators \ell of the horizon is given by

ρ=mam¯bab,σ=mambab.\rho=m^{a}\bar{m}^{b}\nabla_{a}\ell_{b},\ \ \ \ \ \ \sigma=m^{a}m^{b}\nabla_{a}\ell_{b}. (8)

This implies

0\displaystyle 0 =\displaystyle= ρ=122ma(e2bie3b)a(eb1ieb0)=\displaystyle\rho=\frac{1}{2\sqrt{2}}m^{a}(e^{b}_{2}-ie^{b}_{3})\nabla_{a}(e_{b}^{1}-ie_{b}^{0})= (9)
=\displaystyle= 122ma((ωa21ωa20)ı(ωa31ωa30)),\displaystyle\frac{1}{2\sqrt{2}}m^{a}((\omega^{21}_{a}-\omega^{20}_{a})-\imath(\omega^{31}_{a}-\omega^{30}_{a})), (10)

where we have used the definition of the spin connection ωaIJ=eIbaebJ\omega_{a}^{IJ}=e^{Ib}\nabla_{a}e_{b}^{J}. Similarly we have

0\displaystyle 0 =\displaystyle= σ=122ma(e2b+ie3b)a(eb1ieb0)=\displaystyle\sigma=\frac{1}{2\sqrt{2}}m^{a}(e^{b}_{2}+ie^{b}_{3})\nabla_{a}(e_{b}^{1}-ie_{b}^{0})= (11)
=\displaystyle= 122ma((ωa21ωa20)+ı(ωa31ωa30)).\displaystyle\frac{1}{2\sqrt{2}}m^{a}((\omega^{21}_{a}-\omega^{20}_{a})+\imath(\omega^{31}_{a}-\omega^{30}_{a})). (12)

As e2ae_{2}^{a} and e3ae_{3}^{a} form a non degenerate frame for H=ΔMH=\Delta\cap M, and from the definition of pull-back, the previous two equations imply the statement of our lemma. \square

The previous lemma has an immediate consequence on the form of equation (4) for the component i=1i=1 in the frame of the previous lemma. More precisely it says that dK1=0\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{dK}^{1}=0. The good-cut condition (vv) in the definition implies then that

K1=0.\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{K}^{1}=0. (13)

Another important consequence of the previous lemma is equation (3), also derived in ack . We give here for completeness and self consistency a sketch of its derivation. This equation follows from identity

Fab(A+)i=14RabcdΣcd+i,F_{ab}{{}^{i}}(A^{\scriptscriptstyle+})=-\frac{1}{4}R_{ab}^{\ \ cd}\Sigma^{\scriptscriptstyle+i}_{cd}, (14)

where RabcdR_{abcd} is the Riemann tensor and Σ+i=ϵjkiejek+i2e0ei\Sigma^{\scriptscriptstyle+i}=\epsilon^{i}_{\ jk}e^{j}\wedge e^{k}+i2e^{0}\wedge e^{i}, which can be derived using Cartan’s structure equations. A simple algebraic calculation using the null tetrad formalism (see for instance chandra page 43) with the null tetrad of Lemma 1, and the definitions Ψ2=Cabcdambm¯cnd\Psi_{2}=C_{abcd}\ell^{a}m^{b}\bar{m}^{c}n^{d} and Φ11=Rab(anb+mam¯b)/4\Phi_{11}=R_{ab}(\ell^{a}n^{b}+m^{a}\bar{m}^{b})/4, where RabR_{ab} is the Ricci tensor and CabcdC_{abcd} the Weyl tensor, yields

Fabi=(Ψ2Φ11R24)Σi,\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{F_{ab}}^{i}=(\Psi_{2}-\Phi_{11}-\frac{R}{24})\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\Sigma}^{i}, (15)

where Σi=Re[Σ+i]=ϵjkiejek\Sigma^{i}={\rm Re}[\Sigma^{{\scriptscriptstyle+}i}]=\epsilon^{i}_{\ jk}e^{j}\wedge e^{k}. An important point here is that the previous expression is valid for any two sphere S2S^{2} embedded in spacetime in an adapted null tetrad where a\ell^{a} and nan^{a} are normal to S2S^{2}. However, in the special case where S2=HS^{2}=H (where H=ΔMH=\Delta\cap M with Δ\Delta a Type I isolated horizon) it follows from spherical symmetry that (Ψ2Φ11R24)=C(\Psi_{2}-\Phi_{11}-\frac{R}{24})=C with CC a constant on the horizon HH. Moreover, in the gauge defined in the statement of Lemma 1, the only non vanishing component of the previous equation is the i=1i=1 component for which (using Lemma 1) we get

dA+1=Cϵ,dA_{\scriptscriptstyle+}^{1}=C\epsilon, (16)

with ϵ\epsilon the area element of HH. Integrating the previous equation on HH one can completely determine the constant CC, namely

C=(Ψ2Φ11R24)=2πaH,C=(\Psi_{2}-\Phi_{11}-\frac{R}{24})=-\frac{2\pi}{a_{\scriptscriptstyle H}}, (17)

from where equation (3) immediately follows.

Lemma 2: For Type I isolated horizons

KjKkϵijk=c02πaHΣi,\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{K^{j}}\wedge\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{K^{k}}\epsilon_{ijk}=c_{0}\frac{2\pi}{a_{\scriptscriptstyle H}}\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\Sigma}^{i}\,, (18)

for some constant c0c_{0}. One can choose a representative from the equivalence class [][\ell] of null normals to the isolated horizon in order to fix c0=1c_{0}=1 by making use of the translation symmetry of IH along \ell. By studying the stationary spherically symmetric back hole solutions one can show that this corresponds to the choice where the surface gravity of the IH matches the stationary surface gravity (see Appendix A).

Proof: In order to simplify the notation all free indices associated to forms that appear in this proof are pulled back to HH (this allows us to drop the double arrows from equations). In the frame of lemma 1, where e1e^{1} is normal to HH, the only non trivial component of the equation we want to prove is the i=1i=1 component, namely:

KAKBϵAB=c02πaHΣ1,{K^{A}}\wedge{K^{B}}\epsilon_{AB}=c_{0}\frac{2\pi}{a_{\scriptscriptstyle H}}{\Sigma}^{1}, (19)

where A,B=2,3A,B=2,3 and ϵAB=ϵ1AB\epsilon^{AB}=\epsilon^{1AB}. Now, in that gauge, we have that KA=cBAeBK^{A}=c^{A}_{\ B}e^{B} for some matrix of coefficients cBAc^{A}_{\ B}. Notice that the left hand side of the previous equation equals det(c)eAeBϵAB\det(c)e^{A}\wedge e^{B}\epsilon_{AB}. We first prove that det(c)\det(c) is time independent, i.e. that (detc)=0\ell(\det{c})=0. We need to use the isolated horizon boundary condition

[,Db]va=0vaT(Δ)[\mathscr{L}_{\ell},D_{b}]v^{a}=0\ \ \ v^{a}\in T(\Delta) (20)

where DaD_{a} is the derivative operator determined on the horizon by the Levi-Civita derivative operator a\nabla_{a}. One important property of the commutator of two derivative operators is that it also satisfy the Leibnitz rule (it is itself a new derivative operator). Therefore, using the fact that the null vector nan^{a} is normalized so that n=1\ell\cdot n=-1 we get

0=[,Db]ana=na[,Db]a+a[,Db]naa[,Db]na=0,0=[\mathscr{L}_{\ell},D_{b}]\ell^{a}n_{a}=n_{a}[\mathscr{L}_{\ell},D_{b}]\ell^{a}+\ell^{a}[\mathscr{L}_{\ell},D_{b}]n_{a}\ \ \ \Rightarrow\ \ \ \ell^{a}[\mathscr{L}_{\ell},D_{b}]n_{a}=0, (21)

where we have also used that aT(Δ)\ell^{a}\in T(\Delta). Evaluating the equation on the right hand side explicitly, and using the fact that n=dn+d(n)=0\mathscr{L}_{\ell}n=\ell{\lrcorner}dn+d(\ell{\lrcorner}n)=0666 Here we used that dn=0dn=0 which comes from the restriction to ‘good cuts’ in definition of Section II. More precisely, if one introduces a coordinate vv on Δ\Delta such that aav=1\ell^{a}\partial_{a}v=1 and v=0v=0 on some leaf of the foliation, then it follows—from the fact that \ell is a symmetry of the horizon geometry (q,D)(q,D), and the fact that the horizon geometry uniquely determines the foliation into ‘good cuts’—that vv will be constant on all the leaves of the foliation. As nn must be normal to the leaves one has n=dvn=-dv, whence dn=0dn=0. we get

0\displaystyle 0 =\displaystyle= a[,Db]na=a(Dbna)=12a(Db[ea1+ea0])\displaystyle\ell^{a}[\mathscr{L}_{\ell},D_{b}]n_{a}=\ell^{a}\mathscr{L}_{\ell}(D_{b}n_{a})=-\frac{1}{\sqrt{2}}\ell^{a}\mathscr{L}_{\ell}(D_{b}[e^{1}_{a}+e^{0}_{a}])
=\displaystyle= 12a(ωbμ1eaμ+ωbμ0eaμ)=12a(ωb10[ea0+ea1])+12a(ωbA1eaA+ωbA0eaA)\displaystyle\frac{1}{\sqrt{2}}\ell^{a}\mathscr{L}_{\ell}(\omega^{1}_{b\ \mu}e^{\mu}_{a}+\omega^{0}_{b\ \mu}e^{\mu}_{a})=-\frac{1}{\sqrt{2}}\ell^{a}\mathscr{L}_{\ell}(\omega_{b}^{10}[e^{0}_{a}+e^{1}_{a}])+\frac{1}{\sqrt{2}}\ell^{a}\mathscr{L}_{\ell}(\omega^{1}_{b\ A}e^{A}_{a}+\omega^{0}_{b\ A}e^{A}_{a})
=\displaystyle= a(ωb10)na,\displaystyle\ell^{a}\mathscr{L}_{\ell}(\omega_{b}^{10})n_{a},

where in the second line we have used the fact that Dbeaν=ωbμνeaμD_{b}e_{a}^{\nu}=-\omega^{\nu}_{b\ \mu}e_{a}^{\mu} plus the fact that as Lqab=0{\mathfs{L}}_{\ell}q_{ab}=0 the Lie derivative LeA=αϵABeB{\mathfs{L}}_{\ell}e^{A}=\alpha\epsilon^{AB}e_{B} for some α\alpha (moreover, one can even fix α=0\alpha=0 if one wanted to by means of internal gauge transformations). Then it follows that

LK1=0,{\mathfs{L}}_{\ell}K^{1}=0, (22)

a condition that is also valid for the so called weakly isolated horizons better . A similar argument as the one given in equation (21)—but now replacing a\ell^{a} by eBaT(Δ)e^{a}_{B}\in T(\Delta) for B=2,3B=2,3—leads to

0\displaystyle 0 =\displaystyle= eBa[,Db]na=eBa(Dbna)=12eBa(Db[ea1+ea0])\displaystyle e_{B}^{a}[\mathscr{L}_{\ell},D_{b}]n_{a}=e_{B}^{a}\mathscr{L}_{\ell}(D_{b}n_{a})=-\frac{1}{\sqrt{2}}e_{B}^{a}\mathscr{L}_{\ell}(D_{b}[e^{1}_{a}+e^{0}_{a}])
=\displaystyle= 12eBa(ωbμ1eaμ+ωbμ0eaμ)=12eBa(ωb10[ea0+ea1])+12eBa(ωbA1eaA+ωbA0eaA)\displaystyle\frac{1}{\sqrt{2}}e_{B}^{a}\mathscr{L}_{\ell}(\omega^{1}_{b\ \mu}e^{\mu}_{a}+\omega^{0}_{b\ \mu}e^{\mu}_{a})=-\frac{1}{\sqrt{2}}e_{B}^{a}\mathscr{L}_{\ell}(\omega_{b}^{10}[e^{0}_{a}+e^{1}_{a}])+\frac{1}{\sqrt{2}}e_{B}^{a}\mathscr{L}_{\ell}(\omega^{1}_{b\ A}e^{A}_{a}+\omega^{0}_{b\ A}e^{A}_{a})
=\displaystyle= 2eBa(ωbA0eaA)=2[(ωb0B)+αϵBAωb0A],\displaystyle{\sqrt{2}}e_{B}^{a}\mathscr{L}_{\ell}(\omega^{0}_{b\ A}e^{A}_{a})={\sqrt{2}}[\mathscr{L}_{\ell}(\omega_{b}^{0B})+\alpha\epsilon^{BA}\omega_{b}^{0A}],

where, in addition to previously used identities, we have made use of lemma 1, eq. (7). The previous equations imply that the left hand side of equation (19) is Lie dragged along the vector field \ell, and since Σi\Sigma^{i} is also Lie dragged (in this gauge), all this implies that

L(det(c))=(det(c))=0.{\mathfs{L}}_{\ell}(\det(c))=\ell(\det(c))=0. (23)

Now we must use the rest of the symmetry group of Type I isolated horizons. If we denote by jiT(H)j_{i}\in T(H) (i=1,2,3i=1,2,3) the three Killing vectors corresponding to the SO(3)SO(3) symmetry group of Type I isolated horizons. Spherical symmetry of the horizon geometry (q,D)(q,D) implies

Ljiq=0and[Lji,Db]va=0vaT(Δ),{\mathfs{L}}_{j_{i}}q=0\ \ \ {\rm and}\ \ \ [{\mathfs{L}}_{j_{i}},D_{b}]v^{a}=0\ \ \ \forall v^{a}\in T(\Delta), (24)

which,through similar manipulations as the one used above, lead to

ji(detc)=0j_{i}(\det{c})=0 (25)

which completes the prove that detc\det{c} is constant on Δ\Delta. We can now introduce the dimensionless constant 2πc0:=aHdet(c)2\pi c_{0}:=a_{\scriptscriptstyle H}\det(c). Finally one can fix c0=1c_{0}=1 by choosing the appropriate null generator from the equivalence class [][\ell]. \square

IV The conserved presymplectic structure

In this section we show in detail how the IH boundary condition implies the appearance of an SU(2)SU(2) Chern-Simons boundary term in the symplectic structure describing the dynamics of Type I isolated horizons. This result is key for the quantization of the system described in Section VII.

IV.1 The action principle

The conserved pre-symplectic structure in terms of Ashtekar variables can be easily obtained in the covariant phase space formalism. The action principle of general relativity in self dual variables containing an inner boundary satisfying the IH boundary condition (for asymptotically flat spacetimes) takes the form

S[e,A+]=iκMΣi+(e)Fi(A+)+iκτΣi+(e)A+iS[e,A_{\scriptscriptstyle+}]=-\frac{i}{\kappa}\int_{{\mathfs{M}}}\Sigma^{\scriptscriptstyle+}_{i}(e)\wedge F^{i}(A_{\scriptscriptstyle+})+\frac{i}{\kappa}\int_{\tau_{\infty}}\Sigma^{\scriptscriptstyle+}_{i}(e)\wedge A^{i}_{\scriptscriptstyle+} (26)

where Σi+(e)=ϵjkiejek+i2e0ei\Sigma^{\scriptscriptstyle+}_{i}(e)=\epsilon^{i}_{\ jk}e^{j}\wedge e^{k}+i2e^{0}\wedge e^{i} and A+iA^{i}_{\scriptscriptstyle+} is the self-dual connection, and a boundary contribution at a suitable time cylinder τ\tau_{\infty} at infinity is required for the differentiability of the action. No boundary term is necessary if one allows variations that fix an isolated horizon geometry up to diffeomorphisms and Lorentz transformations. This is a very general property and we shall prove it in the next section as we need a little bit of notation that is introduced there.

First variation of the action yields

δS[e,A+]=iκMδΣi+(e)Fi(A+)dA+Σi+δA+i+d(Σi+δA+i)+iκτδ(Σi+(e)A+i),\delta S[e,A_{\scriptscriptstyle+}]=\frac{-i}{\kappa}\int_{{\mathfs{M}}}\delta\Sigma^{\scriptscriptstyle+}_{i}(e)\wedge F^{i}(A_{\scriptscriptstyle+})-d_{A_{\scriptscriptstyle+}}\Sigma^{\scriptscriptstyle+}_{i}\wedge\delta A^{i}_{\scriptscriptstyle+}+d(\Sigma^{\scriptscriptstyle+}_{i}\wedge\delta A^{i}_{\scriptscriptstyle+})+\frac{i}{\kappa}\int_{\tau_{\infty}}\delta(\Sigma^{\scriptscriptstyle+}_{i}(e)\wedge A^{i}_{\scriptscriptstyle+}), (27)

from which the self dual version of Einstein’s equations follow

ϵijkejFi(A+)+ie0Fk(A+)=0\displaystyle\epsilon_{ijk}e^{j}\wedge F^{i}(A_{\scriptscriptstyle+})+ie^{0}\wedge F_{k}(A_{\scriptscriptstyle+})=0
eiFi(A+)=0\displaystyle e_{i}\wedge F^{i}(A_{\scriptscriptstyle+})=0
dA+Σi+=0,\displaystyle d_{A_{\scriptscriptstyle+}}\Sigma^{\scriptscriptstyle+}_{i}=0, (28)

as the boundary terms in the variation of the action cancel.

IV.2 The classical results in a nutshell

In the following subsections a series of technical results are explicitly proven. Here we give an account of these results. The reader who is not interested in the explicit proofs can jump directly to Section V after reading the present subsection. In this work we study general relativity on a spacetime manifold with an internal boundary satisfying the isolated boundary condition corresponding to Type I isolated horizons, and asymptotic flatness at infinity. The phase space of such system is denoted Γ\Gamma and is defined by an infinite dimensional manifold where points pΓp\in\Gamma are given by solutions to Einstein’s equations satisfying the Type I IH boundary condition. Explicitly a point pΓp\in\Gamma can be parametrized by a pair p=(Σ+,A+)p=({\Sigma^{\scriptscriptstyle+}},{A_{\scriptscriptstyle+}}) satisfying the field equations (28) and the requirements of Definition II. In particular fields at the boundary satisfy Einstein’s equations and the constraints given in Section III. Let Tp(Γ){\rm T_{p}}(\Gamma) denote the space of variations δ=(δΣ+,δA+)\delta=(\delta\Sigma^{\scriptscriptstyle+},\delta A_{\scriptscriptstyle+}) at pp (in symbols δTp(Γ)\delta\in{\rm T_{p}}(\Gamma)). A very important point is that the IH boundary conditions severely restrict the form of field variations at the horizon. Thus we have that variations δ=(δΣ+,δA+)Tp(Γ)\delta=(\delta\Sigma^{\scriptscriptstyle+},\delta A_{\scriptscriptstyle+})\in{\rm T_{p}}(\Gamma) are such that for the pull back of fields on the horizon they correspond to linear combinations of SL(2,)SL(2,\mathbb{C}) internal gauge transformations and diffeomorphisms preserving the preferred foliation of Δ\Delta. In equations, for α:Δsl(2,C)\alpha:\Delta\rightarrow sl(2,C) and v:ΔT(H)v:\Delta\rightarrow{\rm T}(H) we have that

δΣ+\displaystyle\delta\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\Sigma^{\scriptscriptstyle+}} =\displaystyle= δαΣ++δvΣ+\displaystyle\delta_{\alpha}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\Sigma^{\scriptscriptstyle+}}+\delta_{v}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\Sigma^{\scriptscriptstyle+}}
δA+\displaystyle\delta\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{A_{\scriptscriptstyle+}} =\displaystyle= δαA++δvA+\displaystyle\delta_{\alpha}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{A_{\scriptscriptstyle+}}+\delta_{v}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{A_{\scriptscriptstyle+}} (29)

where the arrows denote pull-back to Δ\Delta, and the infinitesimal SL(2,C)SL(2,C) transformations are explicitly given by

δαΣ+=[α,Σ+],δαA+=dA+α,\delta_{\alpha}\Sigma^{\scriptscriptstyle+}=[\alpha,\Sigma^{\scriptscriptstyle+}],\ \ \ \delta_{\alpha}A_{\scriptscriptstyle+}=-d_{A_{\scriptscriptstyle+}}\alpha, (30)

while the diffeomorphisms tangent to H take the following form

δvΣi+=LvΣi+=vdA+Σi+ =0 (Gauss)+dA+(vΣ+)i[vA+,Σ+]iδvA+i=LvA+i=vF+i+dA+(vA+)i,{\delta_{v}\Sigma^{\scriptscriptstyle+}_{i}={\mathfs{L}}_{v}\Sigma^{\scriptscriptstyle+}_{i}}=\underbrace{v{\lrcorner}d_{A_{\scriptscriptstyle+}}\Sigma^{\scriptscriptstyle+}_{i}}_{{\mbox{ \tiny\bf$=0$ (Gauss)}}}+d_{A_{\scriptscriptstyle+}}(v{\lrcorner}\Sigma^{\scriptscriptstyle+})_{i}{-[v{\lrcorner}A_{\scriptscriptstyle+},\Sigma^{\scriptscriptstyle+}]_{i}}\ \ \ \ {\delta_{v}A_{\scriptscriptstyle+}^{i}={\mathfs{L}}_{v}A_{\scriptscriptstyle+}^{i}}=v{\lrcorner}F_{\scriptscriptstyle+}^{i}{+d_{A_{\scriptscriptstyle+}}(v{\lrcorner}A_{\scriptscriptstyle+})^{i}}, (31)

where (vω)b1bp1vaωab1bp1(v{\lrcorner}\omega)_{b_{1}\cdots b_{p-1}}\equiv v^{a}\omega_{ab_{1}\cdots b_{p-1}} for any pp-form ωb1bp\omega_{b_{1}\cdots b_{p}}, and the first term in the expresion of the Lie derivative of Σi+\Sigma^{\scriptscriptstyle+}_{i} can be dropped due to the Gauss constraint dAΣi+=0d_{A}\Sigma^{\scriptscriptstyle+}_{i}=0.

So far we have defined the covariant phase space as an infinite dimensional manifold. For it to become a phase space it is necessary to provide it with a presymplectic structure. As the field equations, the presymplectic structure can be obtained from the first variation of the action (27). In particular a symplectic potential density for gravity can be directly read off from the total differential term in (27) cov . The symplectic potential density is therefore

θ(δ)=iκΣi+δA+iδTpΓ.\theta(\delta)=\frac{-i}{\kappa}\Sigma^{\scriptscriptstyle+}_{i}\wedge\delta A^{i}_{\scriptscriptstyle+}\ \ \ \ \forall\ \ \delta\in T_{p}\Gamma. (32)

and the symplectic current takes the form

J(δ1,δ2)=2iκδ[1Σi+δ2]A+iδ1,δ2TpΓ.J(\delta_{1},\delta_{2})=-\frac{2i}{\kappa}\delta_{[1}\Sigma^{\scriptscriptstyle+}_{i}\wedge\delta_{2]}A_{\scriptscriptstyle+}^{i}\ \ \ \ \forall\ \ \delta_{1},\delta_{2}\in T_{p}\Gamma. (33)

Einstein’s equations imply dJ=0dJ=0. Therefore, applying Stokes theorem to the four dimensional (shaded) region in Fig. 1 bounded by M1M_{1} in the past, M2M_{2} in the future, a timelike cylinder at spacial infinity on the right, and the isolated horizon Δ\Delta on the left we obtain

M1δ[1Σi+δ2]A+iM2δ[1Σi+δ2]A+i+Δδ[1Σi+δ2]A+i=0.\displaystyle\int_{M_{1}}\!\!\delta_{[1}\Sigma^{\scriptscriptstyle+}_{i}\wedge\delta_{2]}A_{\scriptscriptstyle+}^{i}-\int_{M_{2}}\!\!\delta_{[1}\Sigma^{\scriptscriptstyle+}_{i}\wedge\delta_{2]}A_{\scriptscriptstyle+}^{i}+{\int_{\Delta}\!\!\delta_{[1}\Sigma^{\scriptscriptstyle+}_{i}\wedge\delta_{2]}A_{\scriptscriptstyle+}^{i}}=0. (34)

Now it turns out that the horizon integral in this expression is a pure boundary contribution: the symplectic flux across the horizon can be expressed as a sum of two terms corresponding to the two-spheres H1=ΔM1H_{1}=\Delta\cap M_{1} and H2=ΔM2H_{2}=\Delta\cap M_{2}. Explicitly (see Proposition 1 proven below), the symplectic flux across the horizon Δ\Delta factorizes into two contributions on Δ\partial\Delta given by SU(2)SU(2) Chern-Simons presymplectic terms according to

Δ2δ[1Σi+δ2]A+i=aH2π[H2H1]δ[1A+iδ2]A+i.\displaystyle{\int_{\Delta}2\delta_{[1}\Sigma^{\scriptscriptstyle+}_{i}\wedge\delta_{2]}A_{\scriptscriptstyle+}^{i}=\frac{a_{\scriptscriptstyle H}}{2\pi}\left[\int_{H_{2}}-\int_{H_{1}}\right]\delta_{[1}A_{{\scriptscriptstyle+}i}\wedge\delta_{2]}A^{i}_{\scriptscriptstyle+}}. (35)

Thus

M12δ[1Σi+δ2]A+iaH2πH1δ[1A+iδ2]A+i=M22δ[1Σi+δ2]A+iaH2πH2δ[1A+iδ2]A+i\displaystyle\int_{M_{1}}2\delta_{[1}\Sigma^{\scriptscriptstyle+}_{i}\wedge\delta_{2]}A_{\scriptscriptstyle+}^{i}-\frac{a_{\scriptscriptstyle H}}{2\pi}\int_{H_{1}}\delta_{[1}A_{{\scriptscriptstyle+}i}\wedge\delta_{2]}A^{i}_{\scriptscriptstyle+}=\int_{M_{2}}2\delta_{[1}\Sigma^{\scriptscriptstyle+}_{i}\wedge\delta_{2]}A_{\scriptscriptstyle+}^{i}-\frac{a_{\scriptscriptstyle H}}{2\pi}\int_{H_{2}}\delta_{[1}A_{{\scriptscriptstyle+}i}\wedge\delta_{2]}A^{i}_{\scriptscriptstyle+} (36)

which implies that the following presymplectic structure is conserved

iκΩM(δ1,δ2)=M2δ[1Σi+δ2]A+iaH2πHδ[1A+iδ2]A+i,{{i\kappa\,\Omega_{M}(\delta_{1},\delta_{2})=\int_{M}2\delta_{[1}\Sigma^{\scriptscriptstyle+}_{i}\wedge\delta_{2]}A_{\scriptscriptstyle+}^{i}-\frac{a_{\scriptscriptstyle H}}{2\pi}\int_{H}\delta_{[1}A_{{\scriptscriptstyle+}i}\wedge\delta_{2]}A^{i}_{\scriptscriptstyle+}}}, (37)

or in other words is independent of MM. The presence of the boundary term in the presymplectic structure might seem at first sight peculiar; however, we will prove in the following section that the previous symplectic structure can be written as

κΩM(δ1,δ2)=M2δ[1Σiδ2]Ki,{\kappa\,\Omega_{M}(\delta_{1},\delta_{2})=\int_{M}2\delta_{[1}\Sigma_{i}\wedge\delta_{2]}K^{i}}, (38)

where we are using the fact that, in the time gauge where e0e^{0} is normal to the space slicing, Σ+i=Re[Σ+i]=Σi\Sigma^{\scriptscriptstyle+i}={\rm Re}[\Sigma^{\scriptscriptstyle+i}]=\Sigma^{i} when pulled back on MM. The previous equation is nothing else but the familiar presymplectic structure of general relativity in terms of the Palatini ΣK\Sigma-K variables. In essence the boundary term arises when connection variables are used in the parametrization of the gravitational degrees of freedom.

Finally, as shown in Section IV.4, the key result for the quantization of Type I IH phase space: the presymplectic structure in Ashtekar Barbero variables takes the form

κβΩM(δ1,δ2)=M2δ[1Σiδ2]AiaHπ(1β2)Hδ1Aiδ2Ai.\displaystyle\!\!\!\!\!\!\!\!\!\!\!\!\!\kappa\beta\Omega_{M}(\delta_{1},\delta_{2})=\int_{M}\!\!\!2\delta_{[1}\Sigma^{i}\wedge\delta_{2]}A_{i}-{\frac{a_{\scriptscriptstyle H}}{\pi({1-\beta^{2}})}}\int_{H}\!\!\!\delta_{1}A_{i}\wedge\delta_{2}A^{i}. (39)

The above equation is the main result of the classical analysis of this paper. It shows that the conserved presymplectic structure of Type I isolated horizons aquires a boundary term given by an SU(2)SU(2) Chern-Simons presymplectic structure when the unconstrained phase space is parametrized in terms of Ashtekar-Barbero variables. In the following subsection we prove this equation.

IV.2.1 On the absence of boundary term on the internal boundary

Before getting involved with the construction of the conserved presymplectic structure let us come back to the issue of the differentiability of the action principle. In the isolated horizon literature it is argued that the IH boundary condition guaranties the differentiability of the action principle without the need of the addition of any boundary term (see afk ). As we show here, this property is satisfied by more general kind of boundary conditions. As mentioned above, the allowed variations are such that the IH geometry is fixed up to diffeomorphisms of Δ\Delta and gauge transformations. This enough for the boundary term arising in the first variation of the action (26) to vanish. The boundary term arising on Δ\Delta upon first variation of the action is

B(δ)=iκΔΣiδAiB(\delta)=-\frac{i}{\kappa}\int_{\Delta}\Sigma_{i}\wedge\delta A^{i} (40)

First let us show that B(δα)=0B(\delta_{\alpha})=0 for δα\delta_{\alpha} as given in (30). We get

B(δα)=iκΔΣidAαi=iκΔ(dAΣi)αiΔΣiαi=0,B(\delta_{\alpha})=\frac{i}{\kappa}\int_{\Delta}\Sigma_{i}\wedge d_{A}\alpha^{i}=-\frac{i}{\kappa}\int_{\Delta}(d_{A}\Sigma_{i})\alpha^{i}-\int_{\partial\Delta}\Sigma_{i}\alpha^{i}=0, (41)

where we integrated by parts in the first identity, the first term in the second identity vanishes due to Eisntein’s equations while the second term vanishes due to the fact that fields are held fixed at the initial and final surfaces M1M_{1} and M2M_{2} and so α=0\alpha=0 when evaluated at Δ\partial\Delta. Similarly we can prove that B(δv)=0B(\delta_{v})=0 for δv\delta_{v} as given in (31) with (this is the only difference) vT(Δ)v\in T(\Delta). We get

B(δv)=iκΔΣi(vFi(A)+dA(vAi))=\displaystyle B(\delta_{v})=-\frac{i}{\kappa}\int_{\Delta}\Sigma_{i}\wedge(v{\lrcorner}F^{i}(A)+d_{A}(v{\lrcorner}A^{i}))=
=iκΔΣi(vFi(A))+iκΔdAΣi(vAi)+ΔΣi(vAi)=0,\displaystyle=-\frac{i}{\kappa}\int_{\Delta}\Sigma_{i}\wedge(v{\lrcorner}F^{i}(A))+\frac{i}{\kappa}\int_{\Delta}d_{A}\Sigma_{i}(v{\lrcorner}A^{i})+\int_{\partial\Delta}\Sigma_{i}(v{\lrcorner}A^{i})=0, (42)

where in the last line the first and second terms vanish due to Einstein’s equations, and the last term vanishes because variations are such that the vector field vv vanishes at Δ\partial\Delta. Notice that we have not made use of the IH boundary condition.

IV.3 The presymplectic structure in self-dual variables

In this section we prove a series of propositions implying that the presymplectic structure of Type I isolated horizons is given by equation (37). In addition, we will prove that the symplectic structure is real and takes the simple form (38) in terms of Palatini variables. Proposition 1: The symplectic flux across a Type I isolated horizon Δ\Delta factorizes into boundary contributions at H1=ΔM1H_{1}=\Delta\cap M_{1} and H2=ΔM2H_{2}=\Delta\cap M_{2} according to

Δδ[1Σi+δ2]A+i=aH2π[H2H1]δ[1A+iδ2]A+i.\displaystyle{\int_{\Delta}\delta_{[1}\Sigma^{\scriptscriptstyle+}_{i}\wedge\delta_{2]}A_{\scriptscriptstyle+}^{i}=\frac{a_{\scriptscriptstyle H}}{2\pi}\left[\int_{H_{2}}-\int_{H_{1}}\right]\delta_{[1}A_{{\scriptscriptstyle+}i}\wedge\delta_{2]}A^{i}_{\scriptscriptstyle+}}. (43)

Proof: On Δ\Delta all variations are linear combinations of SL(2,C)SL(2,C) transformations and tangent diffeos as stated in (29), (30), and (31).

δ=δα+δv\delta=\delta_{\alpha}+\delta_{v}

for α:Δsl(2,C)\alpha:\Delta\rightarrow sl(2,C) and v:ΔT(H)v:\Delta\rightarrow{\rm T}(H). Lets start with SL(2,)SL(2,\mathbb{C}) transformations. Using (30) we get

iκΩΔ(δα,δ)\displaystyle i\kappa\,\Omega_{\Delta}(\delta_{\alpha},\delta) =\displaystyle= Δ[α,Σ]iδA+i+δΣidA(α)i=Δαiδ(dAΣi)+d(δΣiαi)=\displaystyle\int_{\Delta}[\alpha,\Sigma]^{i}\wedge\delta A^{i}_{\scriptscriptstyle+}+\delta\Sigma_{i}\wedge d_{A}(\alpha)^{i}=\int_{\Delta}-{\alpha_{i}\delta(d_{A}\Sigma^{i})}+d(\delta\Sigma_{i}\alpha^{i})= (44)
=\displaystyle= ΔδΣiαi=aH2πΔδF+iαi=aH2πΔdA(δA+i)αi=aH2πΔδA+idAαi=\displaystyle\int_{\partial\Delta}\delta\Sigma^{i}\alpha_{i}=-\frac{a_{\scriptscriptstyle H}}{2\pi}\int_{\partial\Delta}\delta F^{i}_{\scriptscriptstyle+}\alpha_{i}=-\frac{a_{\scriptscriptstyle H}}{2\pi}\int_{\partial\Delta}d_{A}(\delta A^{i}_{\scriptscriptstyle+})\alpha_{i}=-\frac{a_{\scriptscriptstyle H}}{2\pi}\int_{\partial\Delta}\delta A^{i}_{\scriptscriptstyle+}\wedge d_{A}\alpha_{i}=
=\displaystyle= aH2πΔδA+iδαAi,\displaystyle\frac{a_{\scriptscriptstyle H}}{2\pi}\int_{\partial\Delta}\delta A^{i}_{\scriptscriptstyle+}\wedge\delta_{\alpha}A_{i},

where in the first line we used the equations of motion dAΣi=0d_{A}\Sigma^{i}=0 and in the second line we used the IH boundary condition (3). We have therefore shown that

iκΩΔ(δα,δ)=aH2πΔδαA+iδA+i{i\kappa\,\Omega_{\Delta}(\delta_{\alpha},\delta)=-\frac{a_{\scriptscriptstyle H}}{2\pi}\int_{\partial\Delta}\delta_{\alpha}A_{{\scriptscriptstyle+}i}\wedge\delta A^{i}_{\scriptscriptstyle+}}

Similarly, for diffeomorphisms we first notice that (31) implies that

δv=δv+δα(A,v),\delta_{v}=\delta^{\star}_{v}+\delta_{\alpha(A,v)},

where α(A,v)=vA+\alpha(A,v)=v{\lrcorner}A_{\scriptscriptstyle+} and the explicit form of δv\delta^{\star}_{v} is defined as

δvΣi=dA(vΣ)i,δvA+i=vF+i.\delta^{\star}_{v}\Sigma_{i}=d_{A}(v{\lrcorner}\Sigma)^{i},\ \ \ \ \delta^{\star}_{v}A_{\scriptscriptstyle+}^{i}=v{\lrcorner}F_{\scriptscriptstyle+}^{i}.

We have that

iκΩΔ(δv,δ)\displaystyle i\kappa\Omega_{\Delta}(\delta^{\star}_{v},\delta) =\displaystyle= ΔdA(vΣ)iδA+iδΣi(vF+i)=\displaystyle\int_{\Delta}d_{A}(v{\lrcorner}\Sigma)^{i}\wedge\delta A^{i}_{\scriptscriptstyle+}-\delta\Sigma_{i}\wedge(v{\lrcorner}F_{\scriptscriptstyle+}^{i})= (45)
=\displaystyle= Δd((vΣ)iδA+i)+(vΣ)idA(δA+i)δΣi(vF+i)=\displaystyle\int_{\Delta}d((v{\lrcorner}\Sigma)^{i}\wedge\delta A^{i}_{\scriptscriptstyle+})+(v{\lrcorner}\Sigma)^{i}\wedge d_{A}(\delta A^{i}_{\scriptscriptstyle+})-\delta\Sigma_{i}\wedge(v{\lrcorner}F_{\scriptscriptstyle+}^{i})=
=\displaystyle= Δd((vΣ)iδA+i)+(vΣ)iδF+iδΣi(vF+i)=Δd((vΣ)iδA+i)+δ(Σi[vFi(A+)])\displaystyle\int_{\Delta}d((v{\lrcorner}\Sigma)^{i}\wedge\delta A^{i}_{\scriptscriptstyle+}){+(v{\lrcorner}\Sigma)^{i}\wedge\delta F^{i}_{\scriptscriptstyle+}-\delta\Sigma_{i}\wedge(v{\lrcorner}F_{\scriptscriptstyle+}^{i})}=\int_{\Delta}d((v{\lrcorner}\Sigma)^{i}\wedge\delta A^{i}_{\scriptscriptstyle+})+\delta(\Sigma_{i}[v{\lrcorner}F^{i}(A_{\scriptscriptstyle+})])
=\displaystyle= Δ(vΣ+)iδA+i=aH2πΔδvA+iδA+i,\displaystyle\int_{\partial\Delta}(v{\lrcorner}\Sigma_{\scriptscriptstyle+})^{i}\wedge\delta A^{i}_{\scriptscriptstyle+}=-\frac{a_{\scriptscriptstyle H}}{2\pi}\int_{\partial\Delta}\delta^{\star}_{v}A_{\scriptscriptstyle+}^{i}\wedge\delta A^{i}_{\scriptscriptstyle+},

where in the third line we used the vector constraint Σi[vFi(A+)]=0\Sigma_{i}[v{\lrcorner}F^{i}(A_{\scriptscriptstyle+})]=0, while in last line we have used the equations of motion and equation (3). Notice that the calculation leading to equation (44) is also valid for a field dependent α\alpha such as α(A,v)\alpha(A,v). This plus the linearity of the presymplectic structure lead to

iκΩΔ(δv,δ)=aH2πΔδvA+iδA+i{{i\kappa\,\Omega_{\Delta}(\delta_{v},\delta)=-\frac{a_{\scriptscriptstyle H}}{2\pi}\int_{\partial\Delta}\delta_{v}A_{{\scriptscriptstyle+}i}\wedge\delta A^{i}_{\scriptscriptstyle+}}} (46)

and concludes the proof of our proposition \square.

The previous proposition implies that the presymplectic structure (37) is indeed conserved by evolution in Γ\Gamma. Now we are ready to state the next important proposition.

Proposition 2: The presymplectic form ΩM(δ1,δ2)\Omega_{M}(\delta_{1},\delta_{2}) given by

iκΩM(δ1,δ2)=M2δ[1Σi+δ2]A+iaH2πHδ[1A+iδ2]A+i{{i\kappa\,\Omega_{M}(\delta_{1},\delta_{2})=\int_{M}2\delta_{[1}\Sigma^{\scriptscriptstyle+}_{i}\wedge\delta_{2]}A_{\scriptscriptstyle+}^{i}-\frac{a_{\scriptscriptstyle H}}{2\pi}\int_{H}\delta_{[1}A_{{\scriptscriptstyle+}i}\wedge\delta_{2]}A^{i}_{\scriptscriptstyle+}}}

is independent of MM and real. Moreover, the symplectic structure can be described entirely in terms of variables KIm(A+)K\equiv Im(A_{\scriptscriptstyle+}) and Σ\Sigma taking the familliar form

κΩM(δ1,δ2)=M2δ[1Σiδ2]Ki,{\kappa\,\Omega_{M}(\delta_{1},\delta_{2})=\int_{M}2\delta_{[1}\Sigma_{i}\wedge\delta_{2]}K^{i}}, (47)

which is manifestly real and has no boundary contribution.

Proof: The independence of the sysmplectic structure on MM follows directly from Proposition 2 and the argument presented at the end of the previous section. Now let us analyze the reality of the presymplectic structure. The symplectic potential for Ω\Omega written in terms of self dual variables is

iκΘ(δ)=MΣiδA+iaH4πHA+iδA+i{{i\,\kappa\ \Theta(\delta)=\int_{M}\Sigma_{i}\wedge\delta A_{\scriptscriptstyle+}^{i}-\frac{a_{\scriptscriptstyle H}}{4\pi}\int_{H}A_{{\scriptscriptstyle+}i}\wedge\delta A^{i}_{\scriptscriptstyle+}}} (48)

Using A+i=Γi+iKiA_{+}^{i}=\Gamma^{i}+iK^{i} we get

κΘ(δ)=MΣiδKii(MΣiδΓiaH4πHA+iδA+i){{\kappa\ \Theta(\delta)=\int_{M}\Sigma_{i}\wedge\delta K^{i}-i\left(\int_{M}\Sigma_{i}\wedge\delta\Gamma^{i}-\frac{a_{\scriptscriptstyle H}}{4\pi}\int_{H}A_{{\scriptscriptstyle+}i}\wedge\delta A^{i}_{\scriptscriptstyle+}\right)}} (49)

Using a well known property of the spin connection thiemann , and denoting by θ0(δ)\theta_{0}(\delta) the term in parenthesis in the previous equation, we have

Θ0(δ)MΣiδΓiaH4πHA+iδA+i=HeiδeiaH4πHA+iδA+i\displaystyle\Theta_{0}(\delta)\equiv\int_{M}\Sigma_{i}\wedge\delta\Gamma^{i}-\frac{a_{\scriptscriptstyle H}}{4\pi}\int_{H}A_{{\scriptscriptstyle+}i}\wedge\delta A^{i}_{\scriptscriptstyle+}=\int_{H}-e_{i}\wedge\delta e^{i}-\frac{a_{\scriptscriptstyle H}}{4\pi}\int_{H}A_{{\scriptscriptstyle+}i}\wedge\delta A^{i}_{\scriptscriptstyle+}

The proposition follows from that fact that Θ0(δ)\Theta_{0}(\delta) vanishes as proven in the following Lemma \square.

Lemma 3: The phase space one-form Θ0(δ)\Theta_{0}(\delta) defined by

Θ0(δ)HeiδeiaH4πHA+iδA+i\Theta_{0}(\delta)\equiv\int_{H}-e_{i}\wedge\delta e^{i}-\frac{a_{\scriptscriptstyle H}}{4\pi}\int_{H}A_{{\scriptscriptstyle+}i}\wedge\delta A^{i}_{\scriptscriptstyle+} (50)

is closed.

Proof: From the definition of the phase space Γ\Gamma given in Section IV.2, in particular from Eqs. (29) we know that

δe\displaystyle\delta\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{e} =\displaystyle= δαe+δve\displaystyle\delta_{\alpha}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{e}+\delta_{v}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{e}
δA+\displaystyle\delta\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{A_{\scriptscriptstyle+}} =\displaystyle= δαA++δvA+.\displaystyle\delta_{\alpha}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{A_{\scriptscriptstyle+}}+\delta_{v}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{A_{\scriptscriptstyle+}}. (51)

Let us denote by

𝔡Θ0(δ1,δ2)=δ1(Θ0(δ2))δ2(Θ0(δ1)){\mathfrak{d}}\Theta_{0}(\delta_{1},\delta_{2})=\delta_{1}(\Theta_{0}(\delta_{2}))-\delta_{2}(\Theta_{0}(\delta_{1}))

the exterior derivative of Θ0\Theta_{0}. For infinitesimal SL(2,C)SL(2,C) transformations we have

δαe=[α,e],δαA=dAα,,\delta_{\alpha}e=[\alpha,e],\ \ \ \delta_{\alpha}A=-d_{A}\alpha,, (52)

from which it follows

𝔡Θ0(δ,δα)\displaystyle{\mathfrak{d}}\Theta_{0}(\delta,\delta_{\alpha}) =\displaystyle= H2δei[α,e]i+aH2πδA+idAαi=H2δei[α,e]i+aH2πδFi(A+)αi=\displaystyle\int_{H}-2\delta e_{i}\wedge[\alpha,e]^{i}+\frac{a_{\scriptscriptstyle H}}{2\pi}\delta A_{{\scriptscriptstyle+}i}\wedge d_{A}\alpha^{i}=\int_{H}-2\delta e_{i}\wedge[\alpha,e]^{i}+\frac{a_{\scriptscriptstyle H}}{2\pi}\delta F^{i}(A_{{\scriptscriptstyle+}})\alpha_{i}= (53)
=\displaystyle= Hδ(ejek)αiϵijk+aH2πδFi(A+)αi=Hδ[Σi+aH2πFi(A+)]αi=0,\displaystyle\int_{H}\delta(e^{j}\wedge e^{k})\alpha^{i}\epsilon_{ijk}+\frac{a_{\scriptscriptstyle H}}{2\pi}\delta F^{i}(A_{{\scriptscriptstyle+}})\alpha_{i}=\int_{H}\delta[\Sigma^{i}+\frac{a_{\scriptscriptstyle H}}{2\pi}F^{i}(A_{{\scriptscriptstyle+}})]\alpha_{i}=0,

where in the first line we have integrated by parts, and in the second line we used the IH boundary condition. The action of diffeomorphisms tangent to H on the connection and triad take the following form

δvei=Lvei=d(vei)+vdeiδvA+i=LvA+i=vFi(A+)+dA+(vA+i).{\delta_{v}e^{i}={\mathfs{L}}_{v}e^{i}=d(v{\lrcorner}e^{i})+v{\lrcorner}de^{i}}\ \ \ \ {\delta_{v}A_{\scriptscriptstyle+}^{i}={\mathfs{L}}_{v}A_{\scriptscriptstyle+}^{i}=v{\lrcorner}F^{i}(A_{\scriptscriptstyle+})+d_{A_{\scriptscriptstyle+}}(v{\lrcorner}A_{\scriptscriptstyle+}^{i})}. (54)

Now we have

𝔡Θ0(δ,δv)\displaystyle{\mathfrak{d}}\Theta_{0}(\delta,\delta_{v}) =\displaystyle= H2δeiLveiaH2πδA+iLvA+i=\displaystyle\int_{H}-2\delta e_{i}\wedge{\mathfs{L}}_{v}e^{i}-\frac{a_{\scriptscriptstyle H}}{2\pi}\delta A_{{\scriptscriptstyle+}i}\wedge{\mathfs{L}}_{v}A_{\scriptscriptstyle+}^{i}= (55)
=\displaystyle= H2δeid(vei)2δeivdeiaH2π[δA+ivFi(A+)+δA+idA+(vA+i)]=\displaystyle\int_{H}-2\delta e_{i}\wedge d(v{\lrcorner}e^{i})-2\delta e_{i}\wedge v{\lrcorner}de^{i}-\frac{a_{\scriptscriptstyle H}}{2\pi}[\delta A_{{\scriptscriptstyle+}i}\wedge v{\lrcorner}F^{i}(A_{\scriptscriptstyle+})+\delta A_{{\scriptscriptstyle+}i}\wedge d_{A_{\scriptscriptstyle+}}(v{\lrcorner}A_{\scriptscriptstyle+}^{i})]=
=\displaystyle= H2δeid(vei)2vδeideiaH2π[δ(vA+i)Fi(A+)+δFi(A+)vA+i]\displaystyle\int_{H}-2\delta e_{i}\wedge d(v{\lrcorner}e^{i})-2v{\lrcorner}\delta e_{i}\wedge de^{i}-\frac{a_{\scriptscriptstyle H}}{2\pi}[\delta(v{\lrcorner}A_{{\scriptscriptstyle+}i})\wedge F^{i}(A_{\scriptscriptstyle+})+\delta F_{i}(A_{{\scriptscriptstyle+}})\wedge v{\lrcorner}A_{\scriptscriptstyle+}^{i}]
=\displaystyle= H2δ(dei)vei2vδeideiaH2πδ[vA+iFi(A+)]=\displaystyle\int_{H}-2\delta(de_{i})\wedge v{\lrcorner}e^{i}-2v{\lrcorner}\delta e_{i}\wedge de^{i}-\frac{a_{\scriptscriptstyle H}}{2\pi}\delta[v{\lrcorner}A_{{\scriptscriptstyle+}i}\wedge F^{i}(A_{\scriptscriptstyle+})]=
=\displaystyle= Hδ[vΓiΣiv(Γi+iKi)Σi]=0,\displaystyle\int_{H}\delta[v{\lrcorner}\Gamma^{i}\wedge\Sigma_{i}-v{\lrcorner}(\Gamma^{i}+iK^{i})\wedge\Sigma_{i}]=0,

where in addition to integrating by parts and using that H=0\partial H=0, we have used the identity A(vB)(vA)B=0A\wedge(v{\lrcorner}B)-(v{\lrcorner}A)\wedge B=0 for AA a 11-form and BB a 22-form in a two dimensional manifold, and Cartan’s structure equation dei+ϵijkΓjek=0de^{i}+\epsilon_{ijk}\Gamma^{j}e^{k}=0. In the last line we used eq. (3), and eq. (5)—which implies that KiΣi=0K^{i}\Sigma_{i}=0 \square.

IV.4 Presymplectic structure in Ashtekar-Barbero variables

In the previous section (Proposition 2) we have shown how the presymplectic structure

ΩM(δ1,δ2)=1κM[δ1Σiδ2Kiδ2Σiδ1Ki]\Omega_{M}(\delta_{1},\delta_{2})=\frac{1}{\kappa}\int_{M}[\delta_{1}\Sigma^{i}\wedge\delta_{2}K_{i}-\delta_{2}\Sigma^{i}\wedge\delta_{1}K_{i}] (56)

is indeed preserved in the presence of an IH. More precisely in the shaded space-time region in Fig. 1 one has

ΩM2(δ1,δ2)=ΩM1(δ1,δ2).\Omega_{M_{2}}(\delta_{1},\delta_{2})=\Omega_{M_{1}}(\delta_{1},\delta_{2}). (57)

That is, the symplectic flux across the isolated horizon Δ\Delta vanishes due to the isolated horizon boundary condition ack . We will show now, how the very same presymplectic structure takes the form (39) when written in terms of the Ashtekar-Barbero connection variables. For this we need to prove the following lemma:

Lemma 4: The phase space one-form Θ0β(δ)\Theta^{\beta}_{0}(\delta) defined by

βΘ0β(δ)HeiδeiaH2π(1β2)HAiδAi\beta\Theta^{\beta}_{0}(\delta)\equiv\int_{H}-e_{i}\wedge\delta e^{i}-\frac{a_{\scriptscriptstyle H}}{2\pi(1-\beta^{2})}\int_{H}A_{i}\wedge\delta A^{i} (58)

is closed.

Proof: from the definition of the phase space (Section IV.2) we have

δe\displaystyle\delta\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{e} =\displaystyle= δαe+δve\displaystyle\delta_{\alpha}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{e}+\delta_{v}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{e}
δA+\displaystyle\delta\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{A_{\scriptscriptstyle+}} =\displaystyle= δαA++δvA+.\displaystyle\delta_{\alpha}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{A_{\scriptscriptstyle+}}+\delta_{v}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{A_{\scriptscriptstyle+}}. (59)

Let us first study Θ0β(δα)\Theta^{\beta}_{0}(\delta_{\alpha}) where the infinitesimal SU(2)SU(2) transformation are explicitly given by

δαe=[α,e],δαA=dAα,.\delta_{\alpha}e=[\alpha,e],\ \ \ \delta_{\alpha}A=-d_{A}\alpha,. (60)

We have

β𝔡Θ0β(δ,δα)\displaystyle\beta{\mathfrak{d}}\Theta^{\beta}_{0}(\delta,\delta_{\alpha}) =\displaystyle= H2δei[α,e]i+aHπ(1β2)δAidAαi=H2δei[α,e]i+aHπ(1β2)δFi(A)αi=\displaystyle\int_{H}-2\delta e_{i}\wedge[\alpha,e]^{i}+\frac{a_{\scriptscriptstyle H}}{\pi(1-\beta^{2})}\delta A_{i}\wedge d_{A}\alpha^{i}=\int_{H}-2\delta e_{i}\wedge[\alpha,e]^{i}+\frac{a_{\scriptscriptstyle H}}{\pi(1-\beta^{2})}\delta F^{i}(A)\alpha_{i}= (61)
=\displaystyle= Hδ[Σi+aHπ(1β2)Fi(A)]αi=0,\displaystyle\int_{H}\delta[\Sigma^{i}+\frac{a_{\scriptscriptstyle H}}{\pi(1-\beta^{2})}F^{i}(A)]\alpha_{i}=0,

where in the first line we have integrated by parts, and in the second line we used the IH boundary condition. The proof that the presymplectic potential vanishes for δv\delta_{v} mimics exactly the corresponding part of the proof of Lemma 3 \square.

The next step is to write the symplectic structure in terms of Ashtekar-Barbero connection variables necessary for the quantization program of LQG. When there is no boundary the SU(2)SU(2) connection

Aai=Γi+βKaiA_{a}^{i}=\Gamma^{i}+{\beta}K_{a}^{i} (62)

is canonically conjugate to ϵabcβ1Σbci/4\epsilon^{abc}\beta^{-1}\Sigma_{bc}^{i}/4 where β\beta is the so-called Immirzi parameter. In the presence of a boundary the situation is more subtle: the symplectic structure acquires a boundary term.

Proposition 3: In terms of Ashtekar-Barbero variables the presymplectic structure of the spherically symmetric isolated horizon takes the form

κβΩM=M2δ[1Σiδ2]AiaHπ(1β2)Hδ1Aiδ2Ai,\displaystyle\!\!\!\!\!\!\!\!\!\!\!\!\!\kappa\beta\Omega_{M}=\int_{M}\!\!\!2\delta_{[1}\Sigma^{i}\wedge\delta_{2]}A_{i}-{\frac{a_{\scriptscriptstyle H}}{\pi({1-\beta^{2}})}}\int_{H}\!\!\!\delta_{1}A_{i}\wedge\delta_{2}A^{i}, (63)

where κ=32πG\kappa={32\pi G}.

Proof: The result follows from the variation of the presymplectic potential

κβΘ(δ)=MΣi(βδKi)+βΘ0β(δ)=\displaystyle\kappa\beta\Theta(\delta)=\int_{M}\Sigma_{i}\wedge(\beta\delta K^{i})+\beta\Theta_{0}^{\beta}(\delta)=
=MΣiδ(Γi+βKi)aH2π(1β2)HAiδAi,\displaystyle=\int_{M}\Sigma_{i}\wedge\delta(\Gamma^{i}+\beta K^{i})-{\frac{a_{\scriptscriptstyle H}}{2\pi({1-\beta^{2}})}}\int_{H}\!\!\!A_{i}\wedge\delta A^{i},

which is simply the presymplectic potential leading to the conserved presymplectic structure in ΣK\Sigma-K variables (in equation (47)) to which we have added a term proportional to Θ0β(δ)\Theta^{\beta}_{0}(\delta); a closed term which does not affect the presymplectic structure according to Lemma 4 \square. Remark: Notice that one could have introduced a new connection A¯i=Γi+β¯Ki\bar{A}^{i}=\Gamma^{i}+\bar{\beta}K^{i} with a new parameter β¯\bar{\beta} independent of the Immirzi parameter. The statement of the previous lemma would have remained true if on the right hand side of equation (58) one would have replaced β\beta by β¯\bar{\beta} and AiA^{i} by A¯i\bar{A}^{i}. Consequently, the presymplectic structure can also be parametrized in terms of the analog of equation (63) with a boundary term where AiA^{i} is replaced by A¯i\bar{A}^{i} and β\beta on the prefactor of the boundary term is replaced by β¯\bar{\beta}. This implies that the description of the boundary term in terms of Chern-Simons theory allow for the introduction of a new independent parameter β¯\bar{\beta} which has the effect of modifying the Chern-Simons level. This ambiguity in the description of the boundary degrees of freedom has however no effect in the value of the entropy.

IV.5 A side remark on the triad as the boundary degrees of freedom

Here we show that one can write the presymplectic structure

ΩM(δ1,δ2)=1κM[δ1Σiδ2Kiδ2Σiδ1Ki]\Omega_{M}(\delta_{1},\delta_{2})=\frac{1}{\kappa}\int_{M}[\delta_{1}\Sigma^{i}\wedge\delta_{2}K_{i}-\delta_{2}\Sigma^{i}\wedge\delta_{1}K_{i}] (64)

in a way such that a surface term depending only on the pull back of the triad field while the bulk term coincides with the one obtained in the previous section in terms of real connection variables. In order to do this we rewrite the symplectic potential as follows:

κβΦ(δ)\displaystyle\kappa\beta\Phi(\delta) =\displaystyle= MΣiδ(βKi)\displaystyle\int_{M}\Sigma_{i}\wedge\delta(\beta K^{i}) (65)
=\displaystyle= MΣiδ(βKi+Γi)HΣiδΓi\displaystyle\int_{M}\Sigma_{i}\wedge\delta(\beta K^{i}+\Gamma^{i})-\int_{H}\Sigma_{i}\wedge\delta\Gamma^{i}
=\displaystyle= MΣiδAi+Heiδei.\displaystyle\int_{M}\Sigma_{i}\wedge\delta A^{i}+\int_{H}e_{i}\wedge\delta e^{i}.

As a result the symplectic structure becomes noni (and independently wi )

ΩM(δ1,δ2)=1κβM[δ1Σiδ2Aiδ2Σiδ1Ai]+2βκHδ1eiδ2ei.\Omega_{M}(\delta_{1},\delta_{2})=\frac{1}{\kappa\beta}\int_{M}[\delta_{1}\Sigma^{i}\wedge\delta_{2}A_{i}-\delta_{2}\Sigma^{i}\wedge\delta_{1}A_{i}]+\frac{2}{\beta\kappa}\int_{H}\delta_{1}e_{i}\wedge\delta_{2}e^{i}. (66)

The previous equation shows that the boundary degrees of freedom could be described in terms of the pull back of the triad on the horizon. One could try to quantize the IH system in this formulation in order to address the question of black hole entropy calculation. Such project would be certainly interesting. However, the treatment is clearly not immediate as it would require the background independent quantization of the triad field on the boundary for which the usual available techniques do not seem to naturally apply. Nevertheless, the previous equations provides an interesting insight already at the classical level, as the boundary symplectic structure, written in this way, has a remarkable implication for geometric quantities of interest in the first order formulation. To see this let us take SHS\subset H and α:Hsu(2)\alpha:H\to su(2) so that we can introduce the fluxes Σ(S,α)\Sigma(S,\alpha) according to

Σ(S,α)=SHTr[αΣ],\Sigma(S,\alpha)=\int_{S\subset H}{\rm Tr}[\alpha\Sigma], (67)

where Tr[αΣ]=ϵijkαiejek{\rm Tr}[\alpha\Sigma]=\epsilon_{ijk}\alpha^{i}e^{j}\wedge e^{k}. Now (66) implies the Poisson bracket {eai(x),ebj(y)}=ϵabδijδ(x,y)\{e^{i}_{a}(x),e^{j}_{b}(y)\}=\epsilon_{ab}\delta^{ij}\delta(x,y) from which the following remarkable equation follows:

{Σ(S,α),Σ(S,β)}=Σ(SS,[α,β]).\{\Sigma(S,\alpha),\Sigma(S^{\prime},\beta)\}=\Sigma(S\cap S^{\prime},[\alpha,\beta]). (68)

The Poisson brackets among surface fluxes is non vanishing and reproduces the su(2)su(2) Lie algebra! This is an interesting property that we find entirely in terms of classical considerations using smooth field configurations. However, compatibility with the bulk fields seems to single out the treatments of kinematical degrees of freedom in terms of the so called holonomy-flux algebra of classical observables for which flux variables satisfy the exact analog of (68) as described in zapata . This fact strengthens even further the relevance of the uniqueness theorems lost , as they assume the use of the holonomy-flux algebra as the starting point for quantization.

V Gauge symmetries

In this section we rederive the form of the presymplectic symplectic structure written in Ashtekar-Barbero variables by means of gauge symmetry argument. The idea is to first study the gauge symmetries of the presymplectic structure when written in Palatini variables, as in Equation (38). We will show that, due to the nature of variations at the horizon, the boundary term in Equation (39) is completely fixed by the requirement that the gauge symmetry content is unchanged when the presymplectic structure is parametrized by Ashtekar-Barbero variables. This argument is completely equivalent to the content of the previous section and was used in nous as a shortcut construction of the presymplectic structure for Type I isolated horizons in terms of real connection variables. Another important result of this section is the computation of the classical constraint algebra in Subsection V.1 which are essential for clarifying the quantization strategy implemented in Section VII

The gauge symmetry content of the phase space Γ\Gamma is implied by the following proposition.

Proposition 4: Phase space tangent vectors δα,δvTpΓ\delta_{\alpha},\delta_{v}\in T_{p}\Gamma of the form

δαΣ=[α,Σ],δαK=[α,K];\displaystyle\delta_{\alpha}\Sigma=[\alpha,\Sigma],\ \ \delta_{\alpha}K=[\alpha,K];
δvΣ=LvΣ=vdΣ+d(vΣ),δvK=LvK=vdK+d(vK)\displaystyle\delta_{v}\Sigma={\mathfs{L}}_{v}\Sigma=v{\lrcorner}d\Sigma+d(v{\lrcorner}\Sigma),\ \ \delta_{v}K={\mathfs{L}}_{v}K=v{\lrcorner}dK+d(v{\lrcorner}K) (69)

for α:M𝔰𝔲(2)\alpha:M\rightarrow\mathfrak{su}(2) and vVect(M)v\in\mathrm{Vect}(M) tangent to the horizon, are degenerate directions of ΩM\Omega_{M}.

Proof: The proof follows from manipulations very similar in spirit to the ones used for proving the previous propositions. We start with the SU(2)SU(2) transformations δα\delta_{\alpha}, and we get

κΩM(δα,δ)=M[α,Σ]iδKiδΣi[α,K]i=Mδ(ϵijkαjΣkKi))=0,\displaystyle\kappa\Omega_{M}(\delta_{\alpha},\delta)=\int_{M}[\alpha,\Sigma]_{i}\wedge\delta K^{i}-\delta\Sigma_{i}\wedge[\alpha,K]^{i}=\int_{M}\delta({\epsilon_{ijk}\alpha^{j}\Sigma^{k}\wedge K^{i})})=0, (70)

where we used the Gauss constraint ϵijkΣkKi=0\epsilon_{ijk}\Sigma^{k}\wedge K^{i}=0. In order to treat the case of the infinitesimal diffeomorphims tangent to the horizon HH it will be convenient to first write the form of the vector constraint VaV_{a} in terms of ΣK\Sigma-K variables tate . We have

vV=dKivΣi+vKidΣi0\displaystyle v{\lrcorner}V=dK^{i}\wedge v{\lrcorner}\Sigma_{i}+v{\lrcorner}K^{i}\ d\Sigma_{i}\approx 0 (71)

variations of the previous equation yields

vδV\displaystyle v{\lrcorner}\delta V =\displaystyle= d(δK)ivΣi+dKivδΣi+vδKidΣi+vKid(δΣ)i=\displaystyle d(\delta K)^{i}\wedge v{\lrcorner}\Sigma_{i}+dK^{i}\wedge v{\lrcorner}\delta\Sigma_{i}+v{\lrcorner}\delta K^{i}\ d\Sigma_{i}+v{\lrcorner}K^{i}\ d(\delta\Sigma)_{i}= (72)
=\displaystyle= vΣid(δK)iδΣivdKi+vdΣiδKi+d(δΣ)ivKi=0,\displaystyle v{\lrcorner}\Sigma_{i}\wedge d(\delta K)^{i}-\delta\Sigma_{i}\wedge v{\lrcorner}dK^{i}+v{\lrcorner}d\Sigma_{i}\wedge\delta K^{i}+d(\delta\Sigma)_{i}v{\lrcorner}K^{i}=0,

where in the second line we have put all the KK’s to the right, and modified the second and third terms using the identities A(vB)+(vA)B=0A\wedge(v{\lrcorner}B)+(v{\lrcorner}A)\wedge B=0 that is valid for any two 22-forms AA and BB on a 33-manifold, and A(vB)(vA)B=0A\wedge(v{\lrcorner}B)-(v{\lrcorner}A)\wedge B=0 for a 11-form AA and a 33-form BB on a 33-manifold respectively. We are now ready to show that δv\delta_{v} is a null direction of ΩM\Omega_{M}. Explicitly:

κΩM(δv,δ)=M(vdΣ+d(vΣ))iδKiδΣi(vdK+d(vK))i=\displaystyle\kappa\Omega_{M}(\delta_{v},\delta)=\int_{M}(v{\lrcorner}d\Sigma+d(v{\lrcorner}\Sigma))_{i}\wedge\delta K^{i}-\delta\Sigma_{i}\wedge(v{\lrcorner}dK+d(v{\lrcorner}K))^{i}=
=Mv𝑑ΣiδKi+d(vΣ)iδKiδΣivdKiδΣid(vK)i=\displaystyle=\int_{M}v{\lrcorner}d\Sigma_{i}\wedge\delta K^{i}+d(v{\lrcorner}\Sigma)_{i}\wedge\delta K^{i}-\delta\Sigma^{i}\wedge v{\lrcorner}dK_{i}-\delta\Sigma_{i}\wedge d(v{\lrcorner}K)^{i}=
=MvdΣiδKi+vΣid(δK)iδΣivdKi+d(δΣ)ivKivδV=0+\displaystyle=\int_{M}\underbrace{v{\lrcorner}d\Sigma_{i}\wedge\delta K^{i}+v{\lrcorner}\Sigma_{i}\wedge d(\delta K)^{i}-\delta\Sigma^{i}\wedge v{\lrcorner}dK_{i}+d(\delta\Sigma)^{i}\wedge v{\lrcorner}K_{i}}_{v{\lrcorner}\delta V=0}+
+MvΣiδKiδΣivKi=Mδ(vΣiKi)=0,\displaystyle+\int_{\partial M}v{\lrcorner}\Sigma_{i}\wedge\delta K^{i}-\delta\Sigma_{i}\wedge v{\lrcorner}K^{i}=\int_{\partial M}\delta(v{\lrcorner}\Sigma_{i}\wedge K^{i})=0, (73)

where in the last line we have used the identity vAB+AvB=0v{\lrcorner}A\wedge B+A\wedge v{\lrcorner}B=0 valid for an arbitrary 22-form AA and arbitrary 11-form BB on a 22-manifold, the fact that vv is tangent to HH, and the IH boundary condition Eq. (5) implying ΣiKi=0\Sigma_{i}\wedge K^{i}=0 when pulled back on HH \square.

The previous proposition shows that the IH boundary condition breaks neither the symmetry under diffeomorphisms preserving HH nor the SU(2)SU(2) internal gauge symmetry introduced by the use of triad variables.

The gauge invariances of the IH system have been made explicit in the ΣK\Sigma-K parametrization of the presymplectic structure. However, due to the results of Propositions 2 and 3, these can also be made explicit in the parametrization of the presymplectic structure using either self-dual connection variables or real Ashtekar-Barbero variables. It is in fact possible to uniquely determine the horizon contributions to the presymplectic structure in connection variables entirely in terms of the requirement the transformations (69) be gauge invariances of the standard bulk presymplectic contribution plus a suitable boundary term. More precisely, the requirement of SU(2)SU(2) local invariance becomes

0=κβΩM(δα,δ)=MδαΣiδAiδΣiδαAi+κβΩHδTp(Γ),0=\kappa\beta\Omega_{M}(\delta_{\alpha},\delta)=\int_{\scriptscriptstyle M}\!\!\!\delta_{\alpha}\Sigma_{i}\wedge\delta A^{i}\!-\!\delta\Sigma_{i}\wedge\delta_{\alpha}A^{i}\!+\kappa\beta\Omega_{H}\ \ \ \ \forall\ \ \ \delta\in{\rm T_{p}}(\Gamma), (74)

for an (in principle) unknown horizon contribution to the presymplectic structure ΩH\Omega_{H}. This gives some information about the nature of the boundary term, namely

κβΩH=MδαΣiδAiδΣiδαAi=M[α,Σ]iδAi+δΣidAαi\displaystyle-\kappa\beta\Omega_{\scriptscriptstyle H}\!=\!\int_{\scriptscriptstyle M}\!\!\!\delta_{\alpha}\Sigma_{i}\wedge\delta A^{i}\!-\!\delta\Sigma_{i}\wedge\delta_{\alpha}A^{i}\!=\!\int_{\scriptscriptstyle M}\!\!\![\alpha,\Sigma]_{i}\wedge\delta A^{i}\!+\!\delta\Sigma_{i}\wedge d_{A}\alpha^{i}
=Md(αiδΣi)αiδ(dAΣi)=aHπ(1β2)HαiδFi(A)\displaystyle\!=\!\int_{\scriptscriptstyle M}\!\!\!d(\alpha_{i}\delta\Sigma^{i})\!-\!\alpha_{i}\delta(d_{A}\Sigma^{i})\!=\!-\frac{a_{\scriptscriptstyle H}}{\pi(1-\beta^{2})}\int_{\scriptscriptstyle H}\alpha_{i}\delta F^{i}(A)
=aHπ(1β2)HδαAiδAi.\displaystyle\!=\!\frac{a_{\scriptscriptstyle H}}{\pi(1-\beta^{2})}\int_{\scriptscriptstyle H}\delta_{\alpha}A_{i}\wedge\delta A^{i}.

where we used the Gauss law δ(dAΣ)=0\delta(d_{A}\Sigma)=0, condition (6), and that boundary terms at infinity vanish. A similar calculation for diffeomorphisms tangent to the horizon gives an equivalent result. This together with the nature of variations at the horizon (see Eqs. (29)) provides an independent way of establishing the results of Proposition 3. This alternative derivation of the conserved presymplectic structure was used in nous .

V.1 On the first class nature of the IH constraints

The previous equation above can be written as

κβΩ(δα,δ)=Mαiδ(dAΣi)Hαi[aHπ(1β2)δFi+δΣi],\kappa\beta\Omega(\delta_{\alpha},\delta)=-\int_{M}\alpha_{i}\delta(d_{A}\Sigma^{i})-\int_{H}\alpha_{i}\left[\frac{a_{\scriptscriptstyle H}}{\pi(1-\beta^{2})}\delta F^{i}+\delta\Sigma^{i}\right], (75)

or equivalently

Ω(δα,δ)+δG[α,A,Σ]=0,\Omega(\delta_{\alpha},\delta)+\delta G[\alpha,A,\Sigma]=0, (76)

where

G[α,A,Σ]=Mαi(dAΣi/(κβ))+Hαi[aHπκβ(1β2)Fi+1κβΣi].G[\alpha,A,\Sigma]=\int_{M}\alpha_{i}(d_{A}\Sigma^{i}/(\kappa\beta))+\int_{H}\alpha_{i}\left[\frac{a_{\scriptscriptstyle H}}{\pi\kappa\beta(1-\beta^{2})}F^{i}+\frac{1}{\kappa\beta}{\Sigma^{i}}\right]. (77)

In the canonical framework Equation 76 implies that local SU(2)SU(2) transformations δα\delta_{\alpha} are Hamiltonian vector fields generated by the “Hamiltonian” G[α,A,Σ]G[\alpha,A,\Sigma]. It follows immediately from the definition of Poisson brackets that the Poisson algebra of G[α,A,Σ]G[\alpha,A,\Sigma] closes. More precisely, one has

{G[α,A,Σ],G[β,A,Σ]}=Ω(δα,δβ)=δβG(α,A,Σ)\{G[\alpha,A,\Sigma],G[\beta,A,\Sigma]\}=\Omega(\delta_{\alpha},\delta_{\beta})=\delta_{\beta}G(\alpha,A,\Sigma) (78)

from where we get

{G[α,A,Σ],G[β,A,Σ]}=G([α,β],A,Σ).\{G[\alpha,A,\Sigma],G[\beta,A,\Sigma]\}=G([\alpha,\beta],A,\Sigma). (79)

Not surprisingly we get the SU(2)SU(2) Lie algebra a local SU(2)SU(2) transformations. In the previous section we showed that these local transformations are indeed gauge transformations. This implies, in the canonical picture, that canonical variables must satisfy the constraints

G(α,A,Σ)0α:HMsu(2).G(\alpha,A,\Sigma)\approx 0\ \ \ \forall\ \ \ \alpha:H\cup M\to su(2). (80)

Now let us look at diffeomorphisms. A calculation based on the analog of equation (74) for an infinitesimal diffeormorphism preserving HH yields

Ω(δv,δ)+δV[v,A,Σ]=0,\Omega(\delta_{v},\delta)+\delta V[v,A,\Sigma]=0, (81)

where

V[v,A,Σ]=M1κβ[ΣivFivAidAΣi]HvAi[aHπκβ(1β2)Fi+1κβΣi].V[v,A,\Sigma]=\int_{M}\frac{1}{\kappa\beta}\left[\Sigma_{i}\wedge v{\lrcorner}F^{i}-v{\lrcorner}A_{i}d_{A}\Sigma^{i}\right]-\int_{H}v{\lrcorner}A_{i}\left[\frac{a_{\scriptscriptstyle H}}{\pi\kappa\beta(1-\beta^{2})}F^{i}+\frac{1}{\kappa\beta}{\Sigma^{i}}\right]. (82)

Finally, a simple calculation as the one leading to (79), leads to the following first-class constraint algebra

{G[α,A,Σ],G[β,A,Σ]}=G([α,β],A,Σ)\displaystyle\{G[\alpha,A,\Sigma],G[\beta,A,\Sigma]\}=G([\alpha,\beta],A,\Sigma)
{G[α,A,Σ],V[v,A,Σ]}=G(Lvα,A,Σ)\displaystyle\{G[\alpha,A,\Sigma],V[v,A,\Sigma]\}=G({\mathfs{L}}_{v}\alpha,A,\Sigma)
{V[v,A,Σ],V[w,A,Σ]}=V([v,w],A,Σ),\displaystyle\{V[v,A,\Sigma],V[w,A,\Sigma]\}=V([v,w],A,\Sigma), (83)

where we have ignored the Poisson brackets involving the scalar constraint777Recall that the smearing of the scalar constraints must vanish on HH and hence the full constraint algebra including the scalar constraint will remain first class.. Using α\alpha and vv with support only on the horizon HH we can now conclude that the IH boundary condition is first class which justifies the Dirac implementation that will be carried our in the quantum theory.

VI The zeroth and first laws of BH mechanics for (spherical) isolated horizons

The definition given in Section II implies authomaticaly the zeroth law of BH mechanics as κH\kappa_{\scriptscriptstyle H} is constant on Δ\Delta. In turn, the first law cannot be tested unless a definition of energy of the IH is given. Due to the fully dynamical nature of the gravitation field in the neighbourhood of the horizon this might seem problematic. Of course one can in addition postulate an energy formula for the IH in order to satisfy de facto the first law. Fortunately, there is a more elegant way. This consists in requiring the time evolution along vector fields taT(M)t^{a}\in{\rm T}({\mathfs{M}}) which are time translations at infinity and proportional to the null generators \ell at the horizon to correspond to a Hamiltonian time evolution afk . More precisely, denote by δt:ΓT(Γ)\delta_{t}:\Gamma\rightarrow T(\Gamma) the phase space tangent vector field associated to an infinitesimal time evolution along the vector field tat^{a} (which we allow to depend on the phase space point). Then δt\delta_{t} is Hamiltonian if there exists a functional HH such that

δH=ΩM(δ,δt)\delta H=\Omega_{M}(\delta,\delta_{t}) (84)

This requirement fixes a family of good energy formula and translates into the first law of isolated horizons

δEH=κHκδaH+ΦHδQH+other work terms,\delta E_{\scriptscriptstyle H}=\frac{\kappa_{\scriptscriptstyle H}}{\kappa}\delta a_{\scriptscriptstyle H}+\Phi_{\scriptscriptstyle H}\delta Q_{\scriptscriptstyle H}+\mbox{other work terms}, (85)

where we have put the explicit expression of the electromagnetic work term where ΦH\Phi_{\scriptscriptstyle H} is the electromagnetic potential (constant due to the IH boundary condition) and QHQ_{\scriptscriptstyle H} is the electric charge. The above equation implies that κH\kappa_{\scriptscriptstyle H} and ΦH\Phi_{\scriptscriptstyle H} to be functions of the IH area aHa_{H} and charge QHQ_{\scriptscriptstyle H} alone. A unique energy formula is singled out if we require κH\kappa_{\scriptscriptstyle H} to coincide with the surface gravity of Type I stationary BHs, i.e., those in the Reissner-Nordstrom family:

κH=(M2Q2)2M[M+(M2Q2)]Q2.\kappa_{\scriptscriptstyle H}=\frac{\sqrt{(M^{2}-Q^{2})}}{2M[M+\sqrt{(M^{2}-Q^{2})}]-Q^{2}}. (86)

Here we can explicitly prove the above statement in terms of our variables. We shall make here the simplifying assumption that there are no matter fields, i.e. , we work in the vacuum case. The explicit form of δt\delta_{t} is given by

δtΣ\displaystyle\delta_{t}\Sigma =\displaystyle= LtΣ=tdΣ+d(tΣ)\displaystyle{\mathfs{L}}_{t}\Sigma=t{\lrcorner}d\Sigma+d(t{\lrcorner}\Sigma)
δtK\displaystyle\delta_{t}K =\displaystyle= LtK=tdK+d(tK).\displaystyle{\mathfs{L}}_{t}K=t{\lrcorner}dK+d(t{\lrcorner}K). (87)

We can now explicitly write the main condition, namely

16πGΩM(δt,δ)=M(tdΣ+d(tΣ))iδKiδΣi(tdK+d(tK))i=\displaystyle 16\pi G\ \Omega_{M}(\delta_{t},\delta)=\int_{M}(t{\lrcorner}d\Sigma+d(t{\lrcorner}\Sigma))_{i}\wedge\delta K^{i}-\delta\Sigma_{i}\wedge(t{\lrcorner}dK+d(t{\lrcorner}K))^{i}=
=Mt𝑑ΣiδKi+d(tΣ)iδKiδΣitdKiδΣid(tK)i=\displaystyle=\int_{M}t{\lrcorner}d\Sigma_{i}\wedge\delta K^{i}+d(t{\lrcorner}\Sigma)_{i}\wedge\delta K^{i}-\delta\Sigma^{i}\wedge t{\lrcorner}dK_{i}-\delta\Sigma_{i}\wedge d(t{\lrcorner}K)^{i}=
=MΣiδKiδΣiKi=MδΣiKi=\displaystyle=\int_{\partial M}\ell{\lrcorner}\Sigma_{i}\wedge\delta K^{i}-\delta\Sigma_{i}\wedge\ell{\lrcorner}K^{i}=-\int_{\partial M}\delta\Sigma_{i}\wedge\ell{\lrcorner}K^{i}=
=2κHδaH+δEADM,\displaystyle=2\kappa_{\scriptscriptstyle H}\ \delta a_{\scriptscriptstyle H}+\delta E_{\scriptscriptstyle ADM}, (88)

Where we have used the same kind of manipulations used in equation (73) paying special attention to the fact that the relevant vector field t=t=\ell is (at the horizon) no longer tangent to the horizon cross-section, and the fact that the first term in the third line vanishes due to the IH boundary condition 888This follows from equations (5), (175), and (178) implying that ΣiδKi=eαe3δ(e22πaH)+eαe2δ(e32πaH)=\displaystyle\ell{\lrcorner}\Sigma_{i}\wedge\delta K^{i}=-e^{\alpha}e^{3}\wedge\delta\left({e^{2}}{\sqrt{\frac{2\pi}{a_{\scriptscriptstyle H}}}}\right)+e^{\alpha}e^{2}\wedge\delta\left({e^{3}}{\sqrt{\frac{2\pi}{a_{\scriptscriptstyle H}}}}\right)= =eα(2πaHδ(e2e3)+2e2e3δ(2πaH)).\displaystyle=e^{\alpha}\left(\sqrt{\frac{2\pi}{a_{\scriptscriptstyle H}}}\delta\left(e^{2}\wedge e^{3}\right)+2e^{2}\wedge e^{3}\delta\left(\sqrt{\frac{2\pi}{a_{\scriptscriptstyle H}}}\right)\right). Integrating the previous expression on the horizon gives zero..

The condition δHt=ΩM(δt,δ)\delta H_{t}=\Omega_{M}(\delta_{t},\delta) is solved by Ht=EADMEHH_{t}=E_{\scriptscriptstyle ADM}-E_{\scriptscriptstyle H} with

δEH=κHκδaH.\delta E_{\scriptscriptstyle H}=\frac{\kappa_{\scriptscriptstyle H}}{\kappa}\delta a_{\scriptscriptstyle H}. (89)

Demanding time evolution to be Hamiltonian singles out a notion of isolated horizon energy which automatically satisfies, by this requirement, the first law of black hole mechanics (now extended from the static or locally static context to the isolated horizon context). The general treatment and derivation of the first law can be found in afk ; abl2001 .

VII Quantization

The form of the symplectic structure motivates one to handle the quantization of the bulk and horizon degrees of freedom (d.o.f.) separately. We first discuss the bulk quantization. As in standard LQG [8] one first considers (bulk) Hilbert spaces HγB{\mathfs{H}}^{B}_{\gamma} defined on a graph γM\gamma\subset M and then takes the projective limit containing the Hilbert spaces for arbitrary graphs. Along these lines let us first consider HγB{\mathfs{H}}^{\scriptscriptstyle B}_{\gamma} for a fixed graph γM\gamma\subset M with end points on HH, denoted γH\gamma\cap H. The quantum operator associated with Σ\Sigma in (6) is

ϵabΣ^abi(x)=16πGβpγHδ(x,xp)J^i(p)\epsilon^{ab}\hat{\Sigma}^{i}_{ab}(x)=16\pi G\beta\sum_{p\in\gamma\cap H}\delta(x,x_{p})\hat{J}^{i}(p) (90)

where [J^i(p),J^j(p)]=ϵkijJ^k(p)[\hat{J}^{i}(p),\hat{J}^{j}(p)]=\epsilon^{ij}_{\ \ k}\hat{J}^{k}(p) at each pγHp\in\gamma\cap H. Consider a basis of HγB{\mathfs{H}}^{{\scriptscriptstyle B}}_{\gamma} of eigenstates of both JpJpJ_{p}\cdot J_{p} as well as Jp3J_{p}^{3} for all pγHp\in\gamma\cap H with eigenvalues 2jp(jp+1)\hbar^{2}j_{p}(j_{p}+1) and mp\hbar m_{p} respectively. These states are spin network states, here denoted |{jp,mp}1n;|\{j_{p},m_{p}\}_{\scriptscriptstyle 1}^{\scriptscriptstyle n};{\scriptstyle\cdots}\rangle, where jpj_{p} and mpm_{p} are the spins and magnetic numbers labeling the nn edges puncturing the horizon at points xpx_{p} (other labels are left implicit). They are also eigenstates of the horizon area operator a^H\hat{a}_{\scriptscriptstyle H}

a^H|{jp,mp}1n;=8πβp2p=1njp(jp+1)|{jp,mp}1n;.\hat{a}_{\scriptscriptstyle H}|\{j_{p},m_{p}\}_{\scriptscriptstyle 1}^{\scriptscriptstyle n};{\scriptstyle\cdots}\rangle=8\pi\beta\ell_{p}^{2}\,\sum_{p=1}^{n}\sqrt{j_{p}(j_{p}+1)}|\{j_{p},m_{p}\}_{\scriptscriptstyle 1}^{\scriptscriptstyle n};{\scriptstyle\cdots}\rangle.

Now substituting the expression (90) into the quantum version of (6) we get

aHπ(1β2)ϵabF^abi=16πGβpγHδ(x,xp)J^i(p)-\frac{a_{\scriptscriptstyle H}}{\pi(1-\beta^{2})}\epsilon^{ab}\hat{F}^{i}_{ab}=16\pi G\beta\sum_{p\in\gamma\cap H}\delta(x,x_{p})\hat{J}^{i}(p) (91)

As we will show now, the previous equation tells us that the surface Hilbert space, HγHH{\mathfs{H}}^{{\scriptscriptstyle H}}_{\gamma\cap H} that we are looking for is precisely the one corresponding to (the well studied) CS theory in the presence of particles. More precisely, consider the SU(2)SU(2) Chern-Simons action

SCS[A]=aH32π2Gβ(1β2)ΔAidAi+13Ai[A,A]i,\displaystyle S_{\scriptscriptstyle CS}[{A}]\;=\;\frac{-a_{\scriptscriptstyle H}}{32\pi^{2}G\beta(1-\beta^{2})}\int_{\Delta}{A_{i}}\wedge d{A^{i}}+\frac{1}{3}{A_{i}}\wedge[{A},{A}]^{i},

to which we couple a collection of particles by adding the following source term:

Sint[A,Λ1Λn]=p=1nλpcptr[τ3(Λp1dΛp+Λp1AΛp)],\displaystyle S_{\scriptscriptstyle\rm int}[{A},\Lambda_{1}{\scriptstyle\cdots}\Lambda_{n}]\;=\sum_{p=1}^{n}\lambda_{p}\int_{c_{p}}{\rm tr}[\tau_{3}(\Lambda_{p}^{-1}{d\Lambda_{p}}+\Lambda_{p}^{-1}{A}\Lambda_{p})],

where cpΔc_{p}\subset\Delta are the particle world-lines, λp\lambda_{p} coupling constants, and ΛpSU(2)\Lambda_{p}\in SU(2) are group valued d.o.f. of the particles. The particle d.o.f. being added will turn out to correspond precisely to the d.o.f. associated to the bulk J^(p)i\hat{J}(p)^{i} appearing in (90). The horizon and bulk will then be coupled by identifying these d.o.f. The gauge symmetries of the full action are

AgAg1+gdg1,Λpg(xp)Λp,\displaystyle{A}\to gAg^{-1}+gdg^{-1}\!\!,\ \ \ \Lambda_{p}\to g(x_{p})\Lambda_{p}, (92)

where gC(Δ,SU(2))g\in C^{\infty}({\Delta},SU(2)), and

ΛpΛpexp(ϕτ3)\Lambda_{p}\to\Lambda_{p}\exp({\phi\tau^{3}}) (93)

where ϕC(cp,[0,2π])\phi\in C^{\infty}(c_{p},[0,2\pi]).

In order to perform the canonical analysis we assume that Λp(r)=exp(rpατα)\Lambda_{p}(r)=\exp(-r_{p}^{\alpha}\tau_{\alpha}) (α=1,2,3\alpha=1,2,3). Under the left action of the group we have

exp(κeτe)Λp(r)=Λp(f(r,κ))\exp(-\kappa^{e}\tau_{e})\Lambda_{p}(r)=\Lambda_{p}(f(r,\kappa)) (94)

for a function f(r,κ)f(r,\kappa) whose explicit form will not play any role in what follows. The infinitesimal version of the previous action is

τeΛp(r)=Λprαfακe-\tau_{e}\Lambda_{p}(r)=\frac{\partial\Lambda_{p}}{\partial r^{\alpha}}\frac{\partial f^{\alpha}}{\partial\kappa^{e}} (95)

If we define the (spin) momentum SpiS^{i}_{p} as

Spi=παrfακi,S_{p}^{i}=-\pi^{r}_{\alpha}\frac{\partial f^{\alpha}}{\partial\kappa^{i}}, (96)

where παr\pi^{r}_{\alpha} are the conjugate momenta of rαr^{\alpha} then it is easy to recover the following simple Poisson brackets

{Spα,Λp}=ταΛpδpp\displaystyle\{S_{p}^{\alpha},\Lambda_{p^{\prime}}\}=-\tau^{\alpha}\Lambda_{p}\ \delta_{pp^{\prime}}
{Spα,Spβ}=ϵγαβSpγδpp,\displaystyle\{S_{p}^{\alpha},S_{p^{\prime}}^{\beta}\}=\epsilon^{\alpha\beta}_{\ \ \gamma}S_{p}^{\gamma}\ \delta_{pp^{\prime}}, (97)

where the last equation follows from the Jacobi identity. Explicit computation shows that Spi=λpTr[τiΛpτ3Λp1]S_{p}^{i}=\lambda_{p}{\rm Tr}[\tau^{i}\Lambda_{p}\tau_{3}\Lambda_{p}^{-1}]. Therefore; we have three primary constraints per particle

Ψi(Sp,Λp)SpiλpTr[τiΛpτ3Λp1] 0,\displaystyle\Psi^{i}(S_{p},\Lambda_{p})\;\equiv\;S_{p}^{i}-\lambda_{p}{\rm Tr}[\tau^{i}\Lambda_{p}\tau_{3}\Lambda_{p}^{-1}]\;\approx\;0\;, (98)

The primary Hamiltonian is simply given by

H({Sp},{Λp})=pηipΨi(Sp,Λp)H(\{S_{p}\},\{\Lambda_{p}\})=\sum_{p}\eta^{p}_{i}\ \Psi^{i}(S_{p},\Lambda_{p})

the requirement that the constraints be preserved by the time evolution reads

{Ψi(Sp,Λp),H}ϵijkTr[τiΛpτ3Λp1]ηpj\displaystyle\{\Psi_{i}(S_{p},\Lambda_{p}),H\}\approx-\epsilon_{ij}^{\ \,k}{\rm Tr}[\tau_{i}\Lambda_{p}\tau_{3}\Lambda_{p}^{-1}]\eta_{p}^{j} (99)

and the constraint algebra is

{Ψi(Sp,Λp),Ψj(Sp,Λp)}ϵijk(Ψk(Sp,Λp)λpTr[τiΛpτ3Λp1])δpp.\displaystyle\{\Psi_{i}(S_{p},\Lambda_{p}),\Psi_{j}(S_{p^{\prime}},\Lambda_{p^{\prime}})\}\approx\epsilon_{ij}^{\ \,k}(\Psi_{k}(S_{p},\Lambda_{p})-\lambda_{p}{\rm Tr}[\tau_{i}\Lambda_{p}\tau_{3}\Lambda_{p}^{-1}])\ \delta_{pp^{\prime}}.

If we write ηp=ηp+η||p\eta^{p}=\eta^{p}_{\bot}+\eta^{p}_{\scriptscriptstyle||}, where ηp\eta^{p}_{\bot} is the component normal to Λpτ3Λp1\Lambda_{p}\tau_{3}\Lambda_{p}^{-1} while η||p\eta^{p}_{\scriptscriptstyle||} is the parallel one, equations (99) completely fix the Lagrange multipliers ηp\eta^{p}_{\bot}. This means that, per particle, two (out of three) constraints are second class. The fact that η||p\eta^{p}_{\scriptscriptstyle||} remains unfixed by the equations of motion implies the presence of first class contraints which are in fact given by

SpSpλp2 0.\displaystyle S_{p}\cdot S_{p}-\lambda_{p}^{2}\;\approx\;0\;. (100)

This constraint generates rotations Λexpϕτ3Λ\Lambda\rightarrow\exp{\phi\tau_{3}}\Lambda which conserve the quantity Tr[τiΛτ3Λ1]{\rm Tr}[\tau_{i}\Lambda\tau_{3}\Lambda^{-1}]. Now, the fact that there are secons class constraints implies that in order to quantize the theory one has to either work with Dirac brackets, solve the constraints classically before quantizing, or parametrize the phase space in terms of Dirac observables. In this case the third option turns out to be immediate. The reason is that the SpS_{p} turn out to be Dirac observables of the particle system as far as the constraints (98)(\ref{yyyy}) is concerned, namely

{Spi,Ψj(Sp,Λp)}=ϵijkΨk(Sp,Λp)δpp0.\{S^{i}_{p},\Psi^{j}(S_{p^{\prime}},\Lambda_{p^{\prime}})\}=\epsilon^{ijk}\Psi_{k}(S_{p},\Lambda_{p})\ \delta_{pp^{\prime}}\approx 0. (101)

This implies that the Poisson bracket relations (97) remain unchanged when one replaces the brackets {,}\{,\} by Dirac brackets {,}D\{,\}_{\scriptscriptstyle D}. Due to this fact and for notational simplicity we shall keep using the standard Poisson bracket notation.

In summary, the phase space of each particle is T(SU(2))T^{\star}(SU(2)), where the momenta conjugate to Λp\Lambda_{p} are given by SpiS^{i}_{p} satisfying the Poisson bracket relations

{Spi,Λp}=τiΛpδppand{Spi,Spj}=ϵkijSpkδpp.\displaystyle\{S^{i}_{p},\Lambda_{p^{\prime}}\}=-\tau^{i}\Lambda_{p}\ \delta_{pp^{\prime}}\ \ \ {\rm and}\ \ \ \{S_{p}^{i},S_{p^{\prime}}^{j}\}=\epsilon^{ij}_{\ \ k}S^{k}_{p}\ \delta_{pp^{\prime}}. (102)

Explicit computation shows that Spi+λptr[τiΛpτ3Λp1]=0S^{i}_{p}+\lambda_{p}\ {\rm tr}[\tau^{i}\Lambda_{p}\tau^{3}\Lambda_{p}^{-1}]=0 are primary constraints (two of which are second class). In the Hamiltonian framework we use Δ=H×\Delta=H\times\mathbb{R}, and the symmetries (92) and (93) arise from (and imply) the following set of first class constraints on HH:

aHπ(1β2)ϵabFab(x)\displaystyle-\frac{a_{\scriptscriptstyle H}}{\pi(1-\beta^{2})}\epsilon^{ab}{F}_{ab}(x) =\displaystyle= 16πGβp=1nδ(x,xp)Sp,\displaystyle{16\pi G\beta}\sum_{p=1}^{n}\delta(x,x_{p})S_{p}, (103)
SpSpλp2\displaystyle S_{p}\cdot S_{p}-\lambda_{p}^{2} =\displaystyle= 0.\displaystyle 0. (104)

The first constraint tells us that the level of the Chern-Simons theory is 999 If we use the connection A¯i\bar{A}^{i} introduced in the remark below (63) then the level takes the form kaH/(4πp2β(1β¯2))k\equiv a_{\scriptscriptstyle H}/(4\pi\ell_{p}^{2}\beta(1-\bar{\beta}^{2})).

kaH/(4πp2β(1β2)),k\equiv a_{\scriptscriptstyle H}/(4\pi\ell_{p}^{2}\beta(1-\beta^{2})), (105)

and that the curvature of the Chern-Simons connection vanishes everywhere on HH except at the position of the defects where we find conical singularities of strength proportional to the defects’ momenta.

The theory is topological which means in our case that non trivial d.o.f. are only present at punctures. Note that due to (102) and (104) the λp\lambda_{p} are quantized according to λp=sp(sp+1)\lambda_{p}=\sqrt{s_{p}(s_{p}+1)} where sps_{p} is a half integer labelling a unitary irreducible representation of SU(2)SU(2).

From now on we denote HCS(s1sn){\mathfs{H}}^{\scriptscriptstyle CS}({s_{1}{\scriptscriptstyle\cdots}s_{n}}) the Hilbert space of the CS theory associated with a fixed choice of spins sps_{p} at each puncture pγHp\in\gamma\cap H. This will be a proper subspace of the ‘kinematical’ Hilbert space HkinCS(s1sn):=s1sn{\mathfs{H}}^{\scriptscriptstyle CS}_{kin}({s_{1}{\scriptscriptstyle\cdots}s_{n}}):=s_{1}\otimes\cdots\otimes s_{n}. In particular there is an important global constraint that follows from (103) and the fact that the holonomy around a contractible loop that goes around all particles is trivial. It implies

HCS(s1sn)Inv(s1sn).{\mathfs{H}}^{\scriptscriptstyle CS}(s_{1}{\scriptstyle\cdots}s_{n})\subset{\rm Inv}(s_{1}\otimes{\scriptstyle\cdots}\otimes s_{n}). (106)

Moreover, the above containment becomes an equality in the limit kaH/(4πp2β(1β2))k\equiv a_{\scriptscriptstyle H}/(4\pi\ell_{p}^{2}\beta(1-\beta^{2}))\to\infty witten , i.e. in the large BH limit. In that limit we see that the constraint (103) has the simple effect of projecting the particle kinematical states in s1sns_{1}\otimes{\scriptstyle\cdots}\otimes s_{n} into the SU(2)SU(2) singlet.

To make contact with the bulk theory, we first note that the bulk Hilbert space HγB{\mathfs{H}}^{B}_{\gamma} can be written

HγB={jp}pγHH{jp}(pγHjp){\mathfs{H}}_{\gamma}^{B}=\underset{\{j_{p}\}_{p\in\gamma\cap H}}{\bigoplus}{\mathfs{H}}_{\{j_{p}\}}\otimes\left(\underset{p\in\gamma\cap H}{\otimes}{j_{p}}\right) (107)

for certain spaces H{jp}{\mathfs{H}}_{\{j_{p}\}}, and where, for each pp, the generators J^(p)i\hat{J}(p)^{i} appearing in (90) act on the representation space jp{j_{p}}. If we now identify (pjp)\left(\otimes_{p}{j_{p}}\right) with the ‘kinematical’ Chern-Simons Hilbert space HkinCS(j1jn){\mathfs{H}}^{\scriptscriptstyle CS}_{kin}({j_{1}{\scriptscriptstyle\cdots}j_{n}}), the Ji(p)J^{i}(p) operators in (91) are identified with the Si(p)S^{i}(p) of (103). The constraints of the CS theory then restrict HkinCS{\mathfs{H}}^{CS}_{kin} to HCS{\mathfs{H}}^{CS} yielding

Hγ={jp}pγHH{jp}HCS(j1jn),{\mathfs{H}}_{\gamma}=\underset{\{j_{p}\}_{p\in\gamma\cap H}}{\bigoplus}{\mathfs{H}}_{\{j_{p}\}}\otimes{\mathfs{H}}^{\scriptscriptstyle CS}({j_{1}{\scriptscriptstyle\cdots}j_{n}}), (108)

as the full kinematical Hilbert space for fixed γ\gamma.

So far we have dealt with a fixed graph. The Hilbert space satisfying the quantum version of (6) is obtained as the projective limit of the spaces Hγ{\mathfs{H}}_{\gamma}. The Gauss and diffeomorphism constraints are imposed in the same way as in ack ; bhe1 . The IH boundary condition implies that lapse must be zero at the horizon so that the Hamiltonian constraint is only imposed in the bulk.

VII.1 State counting

The entropy of the IH is computed by the formula S=tr(ρIHlogρIH)S={\rm tr}(\rho_{\scriptscriptstyle IH}\log\rho_{\scriptscriptstyle IH}) where the density matrix ρIH\rho_{\scriptscriptstyle IH} is obtained by tracing over the bulk d.o.f., while restricting to horizon states that are compatible with the macroscopic area parameter aHa_{\scriptscriptstyle H}. Assuming that there exist at least one solution of the bulk constraints for every state in the CS theory, the entropy is given by S=log(N)S=\log({\mathfs{N}}) where N{\mathfs{N}} is the number of horizon states compatible with the given macroscopic horizon area aHa_{\scriptscriptstyle H}. After a moment of reflection one sees that

N=n;(j)1ndim[HCS(j1jn)],{{\mathfs{N}}}=\sum_{n;(j)_{\scriptscriptstyle 1}^{\scriptscriptstyle n}}{\rm dim}[{{\mathfs{H}}}^{\scriptscriptstyle CS}(j_{1}{\scriptstyle\cdots}j_{n})], (109)

where the labels j1jpj_{1}\cdots j_{p} of the punctures are constrained by the condition

aHϵ8πβp2p=1njp(jp+1)aH+ϵ.a_{\scriptscriptstyle H}-\epsilon\leq 8\pi\beta\ell_{p}^{2}\,\sum_{p=1}^{n}\sqrt{j_{p}(j_{p}+1)}\leq a_{\scriptscriptstyle H}+\epsilon. (110)

Similar formulae, with a different kk value, were first used in majundar .

Notice that due to (110) we can compute the entropy for aH>>βp2a_{\scriptscriptstyle H}>>\beta\ell_{p}^{2} (not necessarily infinite). The reason is that the representation theory of Uq(SU(2))U_{q}(SU(2))—describing HCS{{\mathfs{H}}}_{\scriptscriptstyle CS} for finite kk—implies

dim[HCS(j1jn)]=dim[Inv(pjp)],{\rm dim}[{{\mathfs{H}}}_{\scriptscriptstyle CS}(j_{1}{\scriptstyle\cdots}j_{n})]={\rm dim}[{\rm Inv}(\otimes_{p}j_{p})], (111)

as long as all the jpj_{p} as well as the interwining internal spins are less than k/2=aH/(8πβ(1β2)p2)k/2=a_{\scriptscriptstyle H}/(8\pi\beta(1-\beta^{2})\ell_{p}^{2}). But for Immirzi parameter in the range |γ|2|\gamma|\leq\sqrt{2} this is precisely granted by (110) according to the Lemma below. All this simplifies the entropy formula considerably. The previous dimension corresponds to the number of independent states one has if one models the black hole by a single SU(2)SU(2) intertwiner!

Lemma 5: The Hilbert spaces HCS(j1jn){{\mathfs{H}}}^{\scriptscriptstyle CS}(j_{1}{\scriptstyle\cdots}j_{n}) of Chern-Simons theory with level kk selected by the restriction

p=1njp(jp+1)k2\sum_{p=1}^{n}\sqrt{j_{p}(j_{p}+1)}\leq\frac{k}{2} (112)

are isomorphic to Inv[(j1jn)]{\rm Inv}[(j_{1}{\scriptstyle\cdots}j_{n})].

Proof: The Chern-Simons Hilbert space CS(j1jn)\mathscr{H}^{CS}(j_{1}\cdots j_{n}) will be isomorphic to Inv[(j1jn)]{\rm Inv}[(j_{1}\cdots j_{n})] if for instance all element of a given basis (of intertwiners) of Inv[(j1jn)]{\rm Inv}[(j_{1}{\scriptstyle\cdots}j_{n})] if (voir for instance babaez )

jpk2p=1,,nj_{p}\leq\frac{k}{2}\ \ \ \forall\ \ p=1,\cdots,n (113)

and

ιak2a=1,,n3.\iota_{a}\leq\frac{k}{2}\ \ \ \forall\ \ a=1,\cdots,n-3. (114)

Equation (113) is immediately implied by (112) as the latter implies

pjpk2.\sum_{p}j_{p}\leq\frac{k}{2}. (115)

The condition (114) requires a more precise analysis. Notice the fact that, being intertwining spins, the ιa\iota_{a} satisfy the following set of nested restrictions which imply the result:

0ι1Min[j1+j2,j3+ι2]j1+j2k2\displaystyle 0\leq\iota_{1}\leq{\rm Min}[j_{1}+j_{2},j_{3}+\iota_{2}]\leq j_{1}+j_{2}\leq\frac{k}{2}
0ι2Min[j3+ι1,j4+ι3]j1+j2+j3k2\displaystyle 0\leq\iota_{2}\leq{\rm Min}[j_{3}+\iota_{1},j_{4}+\iota_{3}]\leq j_{1}+j_{2}+j_{3}\leq\frac{k}{2}
\displaystyle\ \ \ \ \ \ \ \ \cdots
0ιn4Min[jn3+ιn5,jn2+ιn3]p=1n3jpk2\displaystyle 0\leq\iota_{n-4}\leq{\rm Min}[j_{n-3}+\iota_{n-5},j_{n-2}+\iota_{n-3}]\leq\sum_{p=1}^{n-3}j_{p}\leq\frac{k}{2}
0ιn3Min[jn2+ιn4,jn+jn1]p=1n2jpk2,\displaystyle 0\leq\iota_{n-3}\leq{\rm Min}[j_{n-2}+\iota_{n-4},j_{n}+j_{n-1}]\leq\sum_{p=1}^{n-2}j_{p}\leq\frac{k}{2},

where in each line we have used (115) \square.

Remark: An interesting point can be made here as a further developement of the remark below equation (63). Notice that if we had worked with the connection A¯i=Γi\bar{A}^{i}=\Gamma^{i} as our boundary field degree of freedom—corresponding to the choice β¯=0\bar{\beta}=0 in the notation of the remark below (63)—then the boundary Chern-Simons level would be k=aH/(4πp2β)k=a_{\scriptscriptstyle H}/(4\pi\ell_{p}^{2}\beta) (see Footnote 9). This implies that the condition (113) imposed on representations labelling the punctures would take the simple form

jpjmaxaH8πp2β,j_{p}\leq j_{\scriptscriptstyle max}\equiv\frac{a_{\scriptscriptstyle H}}{8\pi\ell_{p}^{2}\beta}, (116)

or equivalently

jjmaxs.t.amax(1)=8πp2βjmax(jmax+1)8πp2βjmax=aHj\leq j_{\scriptscriptstyle max}\ \ \ \ s.t.\ \ \ \ a^{\scriptscriptstyle(1)}_{\scriptscriptstyle max}=8\pi\ell_{p}^{2}\beta\sqrt{j_{\scriptscriptstyle max}(j_{\scriptscriptstyle max}+1)}\approx 8\pi\ell_{p}^{2}\beta j_{max}=a_{\scriptscriptstyle H} (117)

where amax(1)a^{\scriptscriptstyle(1)}_{\scriptscriptstyle max} is the maximum single-puncture eigenvalue allowed. Our effective treatment depends on a classical input: the macroscopic area. One would perhaps hope that this effective treatment would only allow for states where the microscopic area is close to aHa_{\scriptscriptstyle H}, unfortunately such a strong requirement is not satisfied as the allowed eigenvalues can be very far away from aHa_{\scriptscriptstyle H}. However, the effective theory at least forbids quantum states where individual area quanta are larger than aHa_{\scriptscriptstyle H}. This is a nice interplay between the classical input and the associated effective quantum description. Of course this interplay is still qualitatively valid for the case in which one works with the Ashtekar-Barbero connection on the boundary (i.e., β=β¯\beta=\bar{\beta}).

VIII Conclusion

We have shown that the spherically symmetric isolated horizon (or Type I isoated horizon) can be described as a dynamical system by a pre-symplectic form ΩM\Omega_{M} that, when written in the (connection) variables suitable for quantization, acquires a horizon contribution corresponding to an SU(2) Chern-Simons theory. There are different ways to prove this important statement. In nous we first observed that SU(2)SU(2) gauge transformations and diffeomorphism preserving HH are not broken by the IH boundary condition. Moreover, infinitesimal diffeomorphisms tangent to HH and SU(2)SU(2) local transformations continue to be degenerate directions of ΩM\Omega_{M} on shell. This by itself is then sufficient for deriving the boundary term that arises when writing the symplectic structure in terms of Ashtekar-Barbero connection variables. Here we have reviewed this construction in Section V. A result that was not explicitly presented in nous is the precise form of the constraint algebra found in Subsection V.1. There we see in a precise way how the canonical gauge symmetry structure of our system is precisely that of an SU(2)SU(2) Chern-Simons theory: in particular, at the boundary, infinitesimal diffeomorphisms, preserving HH, form a subalgebra of SU(2)SU(2) gauge algebra, as in the topological theory.

A different, more direct approach is based on a subtle fact about the canonical transformation that takes us from the Palatini (Σabi,Kai)(\Sigma^{i}_{ab},K^{i}_{a}) phase space parametrization to the Ashtekar-Barbero (Σabi,Aai)(\Sigma^{i}_{ab},A^{i}_{a}) connection formulation, in the presence of an internal boundary. In the case of Type I isolated horizons, the term to be added to the symplectic potential producing the above transformation gives rise to a boundary contribution that eventually leads to a boundary Chern-Simons term in the presymplectic structure. This is the content of Section IV.4. The boundary Chern-Simons term appears due to the use of connection variables which in turn are the ones in terms of which the quantization program of loop quantum gravity is applicable.

Finally, at a fundamental level, what actually fixes the surface term in the symplectic structure is the requirement that it be conserved in time. The above mentioned proofs show that the various expressions for the symplectic structure using different variables are in fact one and the same symplectic structure. That this symplectic structure is preserved in time was proven in Section IV.

There is a certain freedom in the choice of boundary variables leading to different parametrizations of the boundary degrees of freedom. The most direct description would appear, at first sight, to be the one defined simply in terms of the triad field (pulled back on H) along the lines exhibited in Section IV.5. Such parametrization is however less preferable from the point of view of quantization as one is confronted to the background independent quantization of form fields for which the usual techniques are not directly applicable. In contrast, the parametrization of the boundary degrees of freedom in terms of a connection directly leads to a description in terms of SU(2)SU(2) Chern-Simons theory which, being a well studied topological field theory, drastically simplifies the problem of quantization. However, such description comes with the freedom of the introduction of an extra dimensionless parameter β¯\bar{\beta} (as pointed out in the Remark below equation (63)). Such appearance of extra parameters is very much related to what happens in the general context of the canonical formulation of gravity in terms of connections (see Appendix in rezendeyo ). Therefore, this observation is by no means a new feature proper of IHs. The existence of this extra parameter has a direct influence on the value of the Chern-Simons level; however, the value of the entropy is independent of this extra parameter future2 .

Note that no d.o.f. is available at the horizon in the classical theory as the IH boundary condition completely fixes the geometry at Δ\Delta (the IH condition allows a single (characteristic) initial data once aHa_{\scriptscriptstyle H} is fixed (see fig. 1)). Nevertheless, non trivial d.o.f. arise as would be gauge d.o.f. upon quantization. These are described by SU(2) Chern-Simons theory coupled to (an arbitrary number of) defects through a dimensionless parameter proportional to 4π(1β2)aj/aH4\pi(1-\beta^{2})a_{j}/a_{\scriptscriptstyle H}, where aj=8πp2j(j+1)a_{j}=8\pi\ell_{p}^{2}\sqrt{j({j+1})} is the basic quantum of area carried by the defect. These would be gauge excitations are entirely responsible for the entropy in this approach 101010More insight on the nature of these degrees of freedom could be gained by studying simpler models. In liu a theory with no local degrees of freedom has been introduced. The attractive feature of this model is that it admits an (unconstrained) phase parametrization in terms of the same field content as gravity. Moreover, one can argue that it contains the minimal structure to serve as a toy model to study some generic features of the Type I isolated horizon quantization..

We obtain a remarkably simple formula for the horizon entropy: the number of states of the horizon is simply given in terms of the (well studied) dimension of the Hilbert spaces of Chern-Simons theory with punctures labeled by spins. In the large aHa_{\scriptscriptstyle H} limit the latter is simply equal to the dimension of the singlet component in the tensor product of the representations carried by punctures. In this limit the black hole density matrix ρIH\rho_{\scriptscriptstyle IH} is the identity on Inv(pjp){\rm Inv}(\otimes_{p}j_{p}) for admissible jpj_{p}. Similar counting formulae have been proposed in the literature models ; majundar . Our derivation from first principles clarifies these previous proposals.

Remarkably, the counting of states necessary to compute the entropy of the above Type I isolated horizons can be exactly done barberos using the novel counting techniques introduced in barba . It turns out to be SBH=β0aH/(4βp2)S_{\scriptscriptstyle BH}={\beta_{0}}a_{\scriptscriptstyle H}/({4\beta\ell^{2}_{p}}), where β0=0.274067\beta_{0}=0.274067.... However, the subleading corrections turn out to have the form ΔS=32logaH\Delta S=-\frac{3}{2}\log a_{\scriptscriptstyle H} (instead of the ΔS=12logaH\Delta S=-\frac{1}{2}\log a_{\scriptscriptstyle H} that follows the classic treatment bhe1 ; amit ) matching other approaches majundar . This is due to the full SU(2)SU(2) nature of the IH quantum constraints imposed here. We must mention that the proposal of Majumdar et al. majundar is most closely related to our result. Their intuition was particularly insightful as it yielded a universal form of logarithmic corrections in agreement with those found in different quantum gravity formulations carlip-log . Our work clarifies the relevance of their proposal.

We have concentrated in this work on Type I isolated horizons. The natural question that follows from this analysis is whether we can generalize the SU(2)SU(2) invariant treatment in order to include distortion. The classical formulation and quantization of Type II isolated horizons in the U(1)U(1) (gauge fixed) treatment has been studied in jon . Work in progress future1 shows that, in the SU(2)SU(2) invariant formulation, it is possible to include distortion in a simple way as long as the isolated horizon is non rotating (i.e., when Im[Ψ2]=0{\rm Im}[\Psi_{2}]=0). The rotating case is more subtle but we believe that there are no insurmountable obstacles to its SU(2)SU(2) invariant treatment (this will be studied elsewhere).

IX Acknowledgements

We thank A. Ashtekar, C. Beetle, E. Bianchi, K. Krasnov, M. Montesinos, M. Reisenberger, and C. Rovelli for discussions. This work was supported in part by the Agence Nationale de la Recherche; grant ANR-06-BLAN-0050. J.E. was supported by NSF grant OISE 0601844 and the Alexander von Humboldt Foundation of Germany, and thanks Florida Atlantic University for hospitality during his visit there. A.P. was supported by l’Institut Universitaire de France. D.P. was supported by Marie Curie EU-NCG network.

Appendix A Type I Isolated Horizons: Horizon geometry from the Reissner-Nordstrom family

The spherically symmetric isolated horizons or Type I isolated horizons are easy to visualise in terms of the characteristic formulation of general relativity with initial data given on null surfaces lewa . This observation is very useful if one is looking for a concrete visualisation of the horizon geometry and properties of the matter fields at the horizon. In this appendix we chose to derive the main properties of Type I isolated horizons by studying their geometry in the context of Einstein-Maxwell theory (which is general enough for the most relevant applications of the formalism). An additional motivation for the explicit approach presented here is its complementarity with more abstract discussions available in the literature ack ; better ; ih_prl . In the context of Einstein-Maxwell theory, spacetimes with a Type I IH are solutions to Einstein-Maxwell equations where Reissner-Nordstrom horizon data are given on a null surface Δ=S2×\Delta=S^{2}\times\mathbb{R} and suitable free radiation is given at the transversal null surface for both geometric as well as electromagnetic degrees of freedom. This allows to derive the main equations of IH directly from the Reissner-Nordstrom geometry as far as we are careful enough only to use the information that is intrinsic to the IH geometry.

A.1 The Reissner-Nordstrom solution in Kruskal-like coordinates

The Reissner-Nordstrom metric can be written in Kruskal-like coordinates chandra as

ds2=Ω2(x,t)(dt2+dx2)+r2(dθ2+sin(θ)dϕ2)ds^{2}=\Omega^{2}(x,t)(-dt^{2}+dx^{2})+r^{2}(d\theta^{2}+\sin(\theta)d\phi^{2}) (118)

where

Ω(x,t)=(rr)1+b2earar,\Omega(x,t)=\frac{(r-r_{-})^{\frac{1+b}{2}}e^{-ar}}{ar}, (119)

with a=(r+r)/(2r+2)a=(r_{+}-r_{-})/(2r_{+}^{2}), b=r2/r+2b=r_{-}^{2}/r_{+}^{2}, and the function r(x,t)r(x,t) is determined by the following implicit equation:

F(r)=x2t2;withF(r)=(rr+)e2ar(rr)b.F(r)=x^{2}-t^{2};\ \ \ {\rm with}\ \ \ F(r)=\frac{(r-r_{+})e^{2ar}}{(r-r_{-})^{b}}. (120)

The previous Kruskal-like coordinates are valid for the external region rr+r\geq r_{+}. The metric is smooth at the horizon r=r+r=r_{+} which in the new coordinates corresponds to the null surface x=tx=t. An important identity is:

dr|Δ=2xF(dxdt),dr|_{\Delta}=\frac{2x}{F^{\prime}}\ (dx-dt), (121)

where |Δ|_{\Delta} denotes that the equality holds at the horizon Δ\Delta for which x=tx=t. Here we are interested in the first order formalism. Thus we are interested in an associated tetrad eμIe_{\mu}^{I} with I=0,1,2,3I=0,1,2,3. It is immediate to verify that a possible such tetrad is given by

e0=Ω(x,t)dte1=Ω(x,t)dxe2=rdθe3=rsin(θ)dϕ\displaystyle\begin{array}[]{lll}e^{0}=\Omega(x,t)dt\\ e^{1}=\Omega(x,t)dx\end{array}\ \ \ \ \begin{array}[]{lll}e^{2}=rd\theta\\ e^{3}=r\sin(\theta)d\phi\end{array} (126)

We now want to compute the components of the spin connection ωaIJ\omega_{a}^{IJ} at the horizon. Therefore, we will use Cartan’s first structure equations de+ωe=0de+\omega\wedge e=0 at Δ\Delta. The solution is (all details are given in Section A.3)

ω01|Δ=2xΩFΩ(dtdx)ω02|Δ=2xFΩdθω03|Δ=2xFΩsin(θ)dϕω12|Δ=2xFΩdθω13|Δ=2xFΩsin(θ)dϕω23|Δ=cos(θ)dϕ.\displaystyle\begin{array}[]{lll}\omega^{01}|_{\scriptscriptstyle\Delta}=\frac{2x\Omega^{\prime}}{F^{\prime}\Omega}\ (dt-dx)\\ \omega^{02}|_{\scriptscriptstyle\Delta}=-\frac{2x}{F^{\prime}\Omega}\ d\theta\\ \omega^{03}|_{\scriptscriptstyle\Delta}=-\frac{2x}{F^{\prime}\Omega}\sin(\theta)\ d\phi\end{array}\ \ \ \ \begin{array}[]{lll}\omega^{12}|_{\scriptscriptstyle\Delta}=-\frac{2x}{F^{\prime}\Omega}\ d\theta\\ \omega^{13}|_{\scriptscriptstyle\Delta}=-\frac{2x}{F^{\prime}\Omega}\sin(\theta)\ d\phi\\ \omega^{23}|_{\scriptscriptstyle\Delta}=-\cos(\theta)\ d\phi.\end{array} (133)

At this stage we consider a Lorentz transformation of the form

ΛJI=[cs 0 0sc 0 00 0 1 00 0 0 1],\Lambda^{I}_{J}=\left[\begin{array}[]{ccc}c\ \ s\ \ 0\ \ 0\\ s\ \ c\ \ 0\ \ 0\\ 0\ \ 0\ \ 1\ \ 0\\ 0\ \ 0\ \ 0\ \ 1\end{array}\right], (134)

where c=cosh(α(x))c=\cosh(\alpha(x)) and s=sinh(α(x))s=\sinh(\alpha(x)). It is immediate to see that under such transformation the connection above transforms to

ω~01=α(x)dxω~02=λ(x)dθω~03=λ(x)sin(θ)dϕω~12=λ(x)dθω~13=λ(x)sin(θ)dϕω~23=cos(θ)dϕ\displaystyle\begin{array}[]{lll}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\tilde{\omega}^{01}}=-\alpha^{\prime}(x)\ dx\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\tilde{\omega}^{02}}=-\lambda(x)\ d\theta\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\tilde{\omega}^{03}}=-\lambda(x)\sin(\theta)\ d\phi\end{array}\ \ \ \ \begin{array}[]{lll}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\tilde{\omega}^{12}}=-\lambda(x)\ d\theta\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\tilde{\omega}^{13}}=-\lambda(x)\sin(\theta)\ d\phi\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\tilde{\omega}^{23}}=-cos(\theta)\ d\phi\end{array} (141)

where the arrows below the components denote the pull-back of the one forms to Δ\Delta , and λ(x)=2xFΩexp(α(x))\lambda(x)=\frac{2x}{F^{\prime}\Omega}\exp(\alpha(x)). We can obviously chose this Lorentz transformation in order for λ(x)=λ0\lambda(x)=\lambda_{0} with λ0\lambda_{0} an arbitrary constant. We have

λ0=2xFΩexp(α0(x))\lambda_{0}=\frac{2x}{F^{\prime}\Omega}\exp(\alpha_{0}(x)) (142)

This can be made compatible with the time gauge by changing the spacetime foliation just at the intersection with the horizon Δ\Delta so that e~0=(Λe)0\tilde{e}^{0}=(\Lambda\cdot e)^{0} is the new normal111111 Recently, a similar analysis as the one presented here—and also in beigin —has been done maju . In that reference the authors derive a result which is compatible with the above equations in the singular vanishing extrinsic curvature slicing λ0=0\lambda_{0}=0. Such (null) slicing is however inconsistent with the canonical formulation that is necessary for the LQG quantization of the bulk degrees of freedom. . Now we are ready to write the quantities we were looking for

K1=0K2=λ0dθK3=λ0sin(θ)dϕΓ3=λ0dθΓ2=λ0sin(θ)dϕΓ1=cos(θ)dϕ\displaystyle\begin{array}[]{lll}\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{K^{1}}=0\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{K^{2}}=-\lambda_{0}\ d\theta\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{K^{3}}=-\lambda_{0}\sin(\theta)\ d\phi\end{array}\ \ \ \ \begin{array}[]{lll}\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\Gamma}^{3}=\lambda_{0}\ d\theta\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\Gamma^{2}}=-\lambda_{0}\sin(\theta)\ d\phi\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\Gamma^{1}}=\cos(\theta)\ d\phi\end{array} (149)

where Γi=12ϵijkωjk\Gamma^{i}=-\frac{1}{2}\epsilon^{ijk}\omega_{jk} and Ki=ω0iK^{i}=\omega^{0i}. The self dual connection A+iΓi+iKiA^{i}_{\scriptscriptstyle+}\equiv\Gamma^{i}+iK^{i} and the Ashtekar-Barbero connection become

A+3=λ0(isin(θ)dϕ+dθ)A+2=λ0(sin(θ)dϕidθ)A+1=cos(θ)dϕAβ3=λ0(βsin(θ)dϕ+dθ)Aβ2=λ0(sin(θ)dϕβdθ)Aβ1=cos(θ)dϕ\displaystyle\begin{array}[]{lll}\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{A_{\scriptscriptstyle+}^{3}}=\lambda_{0}\ (-i\ \sin(\theta)d\phi+d\theta)\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{A_{\scriptscriptstyle+}^{2}}=\lambda_{0}\ (-\sin(\theta)\ d\phi-i\ d\theta)\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{A_{\scriptscriptstyle+}^{1}}=\cos(\theta)\ d\phi\end{array}\ \ \ \ \begin{array}[]{lll}\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{A_{\scriptscriptstyle\beta}^{3}}=\lambda_{0}\ (-\beta\ \sin(\theta)d\phi+d\theta)\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{A_{\scriptscriptstyle\beta}^{2}}=\lambda_{0}\ (-\sin(\theta)\ d\phi-\beta\ d\theta)\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{A_{\scriptscriptstyle\beta}^{1}}=\cos(\theta)\ d\phi\end{array} (156)

The curvature of the self-dual and Ashtekar-Barbero connections is (when pulled back to the cross sections HH)

F+3=0F+2=0F+1=sin(θ)dθdϕFβ3=0Fβ2=0Fβ1=(1λ02[1+β2])sin(θ)dθdϕ\displaystyle\begin{array}[]{lll}\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{F^{3}_{\scriptscriptstyle+}}=0\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{F^{2}_{\scriptscriptstyle+}}=0\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{F^{1}_{\scriptscriptstyle+}}=-\sin(\theta)\ d\theta\wedge d\phi\end{array}\ \ \ \ \begin{array}[]{lll}\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{F^{3}_{\scriptscriptstyle\beta}}=0\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{F^{2}_{\scriptscriptstyle\beta}}=0\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{F^{1}_{\scriptscriptstyle\beta}}=-(1-\lambda^{2}_{0}[1+\beta^{2}])\sin(\theta)\ d\theta\wedge d\phi\end{array} (163)

Using that aH=4πr2a_{\scriptscriptstyle H}=4\pi r^{2} we can write the previous equations as

F+i=2πaHΣi\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{F_{\scriptscriptstyle+}^{i}}=-\frac{2\pi}{a_{\scriptscriptstyle H}}\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\Sigma^{i}} (164)

and

Fβi=(1λ02(1+β2))F+i=2π(1λ02(1+β2))aHΣi.\hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{F^{i}_{\scriptscriptstyle\beta}}=(1-\lambda^{2}_{0}(1+\beta^{2}))\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{F_{\scriptscriptstyle+}^{i}}=-\frac{2\pi(1-\lambda^{2}_{0}(1+\beta^{2}))}{a_{\scriptscriptstyle H}}\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\Leftarrow$}\hss}{\Sigma^{i}}. (165)

In the following subsection we will show that λ0=1/2\lambda_{0}=-1/\sqrt{2} defines the frame where the IH surface gravity matches the stationary black hole one. With this value of λ0\lambda_{0}, the previous two equations and equation (149) imply eqs. (3), (5), and (6) respectively. For completeness we write the componets of ΣIJ\Sigma^{IJ}

Σ01=0Σ02=rΩexp(α)dxdθΣ03=rΩexp(α)sin(θ)dxdϕΣ12=rΩexp(α)dxdθΣ13=rΩexp(α)sin(θ)dxdϕΣ23=r2sin(θ)dθdϕΣ+3=rΩexp(α)dxdθ+irΩexp(α)sin(θ)dxdϕΣ+2=exp(α)Ωrsin(θ)dxdϕ+irΩexp(α)dxdθΣ+1=r2sin(θ)dθdϕ\displaystyle\begin{array}[]{lll}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\Sigma^{01}}=0\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\Sigma^{02}}=r\Omega\exp(\alpha)\ dx\wedge d\theta\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\Sigma^{03}}=r\Omega\exp(\alpha)\sin(\theta)\ dx\wedge d\phi\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\Sigma^{12}}=r\Omega\exp(\alpha)\ dx\wedge d\theta\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\Sigma^{13}}=r\Omega\exp(\alpha)\sin(\theta)\ dx\wedge d\phi\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\Sigma^{23}}=r^{2}\sin(\theta)\ d\theta\wedge d\phi\end{array}\ \ \ \ \ \begin{array}[]{lll}\hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\Sigma^{3}_{\scriptscriptstyle+}}=r\Omega\exp(\alpha)\ dx\wedge d\theta+ir\Omega\exp(\alpha)\sin(\theta)\ dx\wedge d\phi\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\Sigma^{2}_{\scriptscriptstyle+}}=-\exp(\alpha)\Omega r\sin(\theta)\ dx\wedge d\phi+ir\Omega\exp(\alpha)\ dx\wedge d\theta\\ \hbox to0.0pt{\lower 6.45831pt\hbox{$\leftarrow$}\hss}{\Sigma^{1}_{\scriptscriptstyle+}}=r^{2}\sin(\theta)\ d\theta\wedge d\phi\end{array} (175)

where on the right we have written the corresponding self-dual components.

A.2 Surface gravity and the value of λ0\lambda_{0}

For stationary black holes, the surface gravity κH\kappa_{\scriptscriptstyle H} is defined by the equation

aab=κHb\ell^{a}\nabla_{a}\ell^{b}=\kappa_{\scriptscriptstyle H}\,\ell^{b} (176)

where a\ell^{a} is the Killing vector field tangent to the horizon. For isolated horizons there is no unique notion of a\ell^{a}. We shall define a\ell_{a} in terms of the tetrad in the usual way with a(ea1ea0)/2\ell_{a}\equiv(e^{1}_{a}-e^{0}_{a})/\sqrt{2} 121212The future pointing null generators of the horizon a\ell^{a} are such that a(/x)a+(/t)a.\ell^{a}\propto(\partial/\partial x)^{a}+(\partial/\partial t)^{a}. This implies that adxadta\ell_{a}\propto dx_{a}-dt_{a} from wich we get a=(ea1ea0)/2\ell_{a}=(e_{a}^{1}-e_{a}^{0})/\sqrt{2} and na=(ea1+ea0)/2n_{a}=-(e_{a}^{1}+e_{a}^{0})/\sqrt{2} so that n=1n\cdot\ell=-1.. However, this definition still allows the freedom associated to the Lorentz transformations (134) which send aexp(α(x))a\ell^{a}\rightarrow\exp(-\alpha(x))\ell^{a}. We can fix this freedom by demanding the surface gravity to match that of a Reissner-Nordstrom black hole with mass MM and charge QQ for which

κH=(M2Q2)2M[M+(M2Q2)]Q2.\kappa_{\scriptscriptstyle H}=\frac{\sqrt{(M^{2}-Q^{2})}}{2M[M+\sqrt{(M^{2}-Q^{2})}]-Q^{2}}. (177)

Indeed this choice is the one that makes the zero, and first law of IH look just as the corresponding laws of stationary black hole mechanics.

This choice is then physically motivated. In turn this will fix the value of λ0\lambda_{0} in (165). If we define na(ea0+ea1)/2n_{a}\equiv-(e_{a}^{0}+e_{a}^{1})/\sqrt{2} then we have that (176) implies

anbab=κH\displaystyle\ell^{a}n_{b}\nabla_{a}\ell^{b}=-\kappa_{\scriptscriptstyle H}
12a(eb0+eb1)a(e1be0b)=κH\displaystyle-\frac{1}{2}\ell^{a}(e^{0}_{b}+e^{1}_{b})\nabla_{a}(e^{1b}-e^{0b})=-\kappa_{\scriptscriptstyle H}
aωa01=κH.\displaystyle\ell^{a}\omega_{a}^{01}=\kappa_{\scriptscriptstyle H}. (178)

Notice that after the Lorentz transformation (134) we have

κH\displaystyle\kappa_{\scriptscriptstyle H} =\displaystyle= aωa01=αadxa=αgabadxb=exp(α)Ω2αgax(dtadxa)=exp(α)Ω2αgxx\displaystyle\ell^{a}\omega_{a}^{01}=-\alpha^{\prime}\ell^{a}dx_{a}=-\alpha^{\prime}g^{ab}\ell_{a}dx_{b}=-\exp{(-\alpha)}\frac{\Omega}{\sqrt{2}}\alpha^{\prime}g^{ax}(dt_{a}-dx_{a})=\exp{(-\alpha)}\frac{\Omega}{\sqrt{2}}\alpha^{\prime}g^{xx} (179)
=\displaystyle= (exp(α))Ω2gxx=(exp(α))12Ω.\displaystyle-(\exp{(-\alpha)})^{\prime}\frac{\Omega}{\sqrt{2}}g^{xx}=-(\exp{(-\alpha)})^{\prime}\frac{1}{\sqrt{2}\Omega}.

Now we can fix α(x)=α0(x)\alpha(x)=\alpha_{0}(x) so that κH\kappa_{\scriptscriptstyle H} takes the RN value. Recalling equation (142) and using the above equations, a simple calculation shows that this happens for

λ0=12\lambda_{0}=-\frac{1}{\sqrt{2}} (180)

which implies the desired result

Fβi=12(1β2)F+i.\mbox{$F^{i}_{\scriptscriptstyle\beta}=\frac{1}{2}(1-\beta^{2})\ F_{\scriptscriptstyle+}^{i}$}. (181)

Notice that

ab=ωa01b,\nabla_{a}\ell_{b}=\omega^{01}_{a}\ell_{b}, (182)

and that (according to (141)) we also have

dω01=0.d\omega^{01}=0. (183)

All this implies that Lω01=d(ω01)+dω01=dκH=0{\mathfs{L}}_{\ell}\omega^{01}=d(\ell{\lrcorner}\omega^{01})+\ell{\lrcorner}d\omega^{01}=d\kappa_{\scriptscriptstyle H}=0 as expected from [L,D]=0[{\mathfs{L}}_{\ell},D]=0 (general proof in Lemma 2). In other words, the \ell we have chosen by means of fixing the boost freedom exp(α(x))\ell\rightarrow\exp(-\alpha(x))\ell is a member of the equivalence class [][\ell] in Definition II.

A.3 Solving Cartan’s equation

For this we first compute dede, namely:

de0=Ω(x,t)drdt|Δ=2ΩxFdxdt\displaystyle de^{0}=\Omega^{\prime}(x,t)dr\wedge dt|_{\scriptscriptstyle\Delta}=2\Omega^{\prime}\frac{x}{F^{\prime}}\ dx\wedge dt
de1|Δ=2ΩxFdxdt\displaystyle de^{1}|_{\scriptscriptstyle\Delta}=2\Omega^{\prime}\frac{x}{F^{\prime}}\ dx\wedge dt
de2|Δ=2xF(dxdθdtdθ)\displaystyle de^{2}|_{\scriptscriptstyle\Delta}=\frac{2x}{F^{\prime}}\ (dx\wedge d\theta-dt\wedge d\theta)
de3|Δ=rcos(θ)dθdϕ+2xFsin(θ)(dxdϕdtdϕ).\displaystyle de^{3}|_{\scriptscriptstyle\Delta}=-r\cos(\theta)\ d\theta\wedge d\phi+\frac{2x}{F^{\prime}}\sin(\theta)\ (dx\wedge d\phi-dt\wedge d\phi). (184)

Now we are ready to explicitly write Cartan’s first structure equations. They are

0|Δ=2ΩxFdxdt+Ωω01dx+rω02dθ+rsin(θ)ω03dϕ\displaystyle 0|_{\scriptscriptstyle\Delta}=2\Omega^{\prime}\frac{x}{F^{\prime}}\ dx\wedge dt+\Omega\ \omega^{01}\wedge dx+r\omega^{02}\wedge d\theta+r\sin(\theta)\omega^{03}\wedge d\phi
0|Δ=2ΩxFdxdt+Ωω01dt+rω12dθ+rsin(θ)ω13dϕ\displaystyle 0|_{\scriptscriptstyle\Delta}=2\Omega^{\prime}\frac{x}{F^{\prime}}\ dx\wedge dt+\Omega\ \omega^{01}\wedge dt+r\omega^{12}\wedge d\theta+r\sin(\theta)\omega^{13}\wedge d\phi
0|Δ=2xF(dxdθdtdθ)+Ωω02dt+Ωω21dx+rsin(θ)ω23dϕ\displaystyle 0|_{\scriptscriptstyle\Delta}=\frac{2x}{F^{\prime}}\ (dx\wedge d\theta-dt\wedge d\theta)+\Omega\omega^{02}\wedge dt+\Omega\ \omega^{21}\wedge dx+r\sin(\theta)\omega^{23}\wedge d\phi
0|Δ=rcos(θ)dθdϕ+2xFsin(θ)(dxdϕdtdϕ)+\displaystyle 0|_{\scriptscriptstyle\Delta}=-r\cos(\theta)\ d\theta\wedge d\phi+\frac{2x}{F^{\prime}}\sin(\theta)\ (dx\wedge d\phi-dt\wedge d\phi)+
+Ωω03dt+Ωω31dx+rω32dθ.\displaystyle\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ +\Omega\omega^{03}\wedge dt+\Omega\ \omega^{31}\wedge dx+r\omega^{32}\wedge d\theta. (185)

Let us now study the previous equation individually. The six components of the first, equation (A.3) become

0|Δ=dxdt(2ΩxFωt01Ω)\displaystyle 0|_{\scriptscriptstyle\Delta}=dx\wedge dt\ (2\Omega^{\prime}\frac{x}{F^{\prime}}-\omega^{01}_{t}\Omega)
0|Δ=dxdθ(Ωωθ01+ωx02r)\displaystyle 0|_{\scriptscriptstyle\Delta}=dx\wedge d\theta\ (-\Omega\omega^{01}_{\theta}+\omega^{02}_{x}r)
0|Δ=dxdϕ(Ωωϕ01+ωx03rsin(θ))\displaystyle 0|_{\scriptscriptstyle\Delta}=dx\wedge d\phi\ (-\Omega\omega^{01}_{\phi}+\omega^{03}_{x}r\sin(\theta))
0|Δ=dtdθ(rωt02)\displaystyle 0|_{\scriptscriptstyle\Delta}=dt\wedge d\theta\ (r\omega^{02}_{t})
0|Δ=dtdϕ(ωt03rsin(θ))\displaystyle 0|_{\scriptscriptstyle\Delta}=dt\wedge d\phi\ (\omega^{03}_{t}r\sin(\theta))
0|Δ=dθdϕ(rωϕ02+ωθ03rsin(θ)).\displaystyle 0|_{\scriptscriptstyle\Delta}=d\theta\wedge d\phi\ (-r\omega^{02}_{\phi}+\omega^{03}_{\theta}r\sin(\theta)). (186)

The six components of the second, equation (LABEL:2), become

0|Δ=dxdt(2ΩxF+ωx01Ω)\displaystyle 0|_{\scriptscriptstyle\Delta}=dx\wedge dt\ (2\Omega^{\prime}\frac{x}{F^{\prime}}+\omega^{01}_{x}\Omega)
0|Δ=dxdθ(ωx12r)\displaystyle 0|_{\scriptscriptstyle\Delta}=dx\wedge d\theta\ (\omega^{12}_{x}r)
0|Δ=dxdϕ(ωx13rsin(θ))\displaystyle 0|_{\scriptscriptstyle\Delta}=dx\wedge d\phi\ (\omega^{13}_{x}r\sin(\theta))
0|Δ=dtdθ(Ωωθ01+rωt12)\displaystyle 0|_{\scriptscriptstyle\Delta}=dt\wedge d\theta\ (-\Omega\omega^{01}_{\theta}+r\omega^{12}_{t})
0|Δ=dtdϕ(Ωωϕ01+rsin(θ)ωt13)\displaystyle 0|_{\scriptscriptstyle\Delta}=dt\wedge d\phi\ (-\Omega\omega^{01}_{\phi}+r\sin(\theta)\omega^{13}_{t})
0|Δ=dθdϕ(rωϕ12+ωθ13rsin(θ)).\displaystyle 0|_{\scriptscriptstyle\Delta}=d\theta\wedge d\phi\ (-r\omega^{12}_{\phi}+\omega^{13}_{\theta}r\sin(\theta)). (187)

The six components of the third, equation (A.3), become

0|Δ=dxdt(ωx02Ω+ωt21Ω)\displaystyle 0|_{\scriptscriptstyle\Delta}=dx\wedge dt\ (\omega^{02}_{x}\Omega+\omega^{21}_{t}\Omega)
0|Δ=dxdθ(2xFωθ21Ω)\displaystyle 0|_{\scriptscriptstyle\Delta}=dx\wedge d\theta\ (2\frac{x}{F^{\prime}}-\omega^{21}_{\theta}\Omega)
0|Δ=dxdϕ(ωϕ21Ω+ωx23rsin(θ))\displaystyle 0|_{\scriptscriptstyle\Delta}=dx\wedge d\phi\ (-\omega^{21}_{\phi}\Omega+\omega^{23}_{x}r\sin(\theta))
0|Δ=dtdθ(2xFΩωθ02)\displaystyle 0|_{\scriptscriptstyle\Delta}=dt\wedge d\theta\ (-2\frac{x}{F^{\prime}}-\Omega\omega^{02}_{\theta})
0|Δ=dtdϕ(Ωωϕ02+rsin(θ)ωt23)\displaystyle 0|_{\scriptscriptstyle\Delta}=dt\wedge d\phi\ (-\Omega\omega^{02}_{\phi}+r\sin(\theta)\omega^{23}_{t})
0|Δ=dθdϕ(ωθ23rsin(θ)).\displaystyle 0|_{\scriptscriptstyle\Delta}=d\theta\wedge d\phi\ (\omega^{23}_{\theta}r\sin(\theta)). (188)

Finally, the six components of the fourth, equation (A.3), become

0|Δ=dxdt(ωx03Ωωt31Ω)\displaystyle 0|_{\scriptscriptstyle\Delta}=dx\wedge dt\ (\omega^{03}_{x}\Omega-\omega^{31}_{t}\Omega)
0|Δ=dxdθ(ωθ31Ω+ωx32r)\displaystyle 0|_{\scriptscriptstyle\Delta}=dx\wedge d\theta\ (-\omega^{31}_{\theta}\Omega+\omega^{32}_{x}r)
0|Δ=dxdϕ(2xFsin(θ)ωϕ31Ω)\displaystyle 0|_{\scriptscriptstyle\Delta}=dx\wedge d\phi\ (2\frac{x}{F^{\prime}}\sin(\theta)-\omega^{31}_{\phi}\Omega)
0|Δ=dtdθ(ωθ03Ω+rωt32)\displaystyle 0|_{\scriptscriptstyle\Delta}=dt\wedge d\theta\ (-\omega^{03}_{\theta}\Omega+r\omega^{32}_{t})
0|Δ=dtdϕ(2xFsin(θ)Ωωϕ03)\displaystyle 0|_{\scriptscriptstyle\Delta}=dt\wedge d\phi\ (-2\frac{x}{F^{\prime}}\sin(\theta)-\Omega\omega^{03}_{\phi})
0|Δ=dθdϕ(rcos(θ)ωϕ32r).\displaystyle 0|_{\scriptscriptstyle\Delta}=d\theta\wedge d\phi\ (-r\cos(\theta)-\omega^{32}_{\phi}r). (189)

At this point we make the following ansatz ω01θ=0,ω01ϕ=0,ω23x=0,andω23t=0\omega^{\theta}_{01}=0,\omega^{\phi}_{01}=0,\omega^{x}_{23}=0,\ {\rm and}\ \omega^{t}_{23}=0. From which we get the solution (133).

References

  • (1) M. J. Reid, “Is there a Supermassive Black Hole at the Center of the Milky Way?,” Int. J. Mod. Phys.  D 18 (2009) 889 [arXiv:0808.2624 [astro-ph]]. A. Mueller, “Experimental evidence of black holes,” PoS P2GC (2006) 017 [arXiv:astro-ph/0701228]. A. E. Broderick, A. Loeb and R. Narayan, “The Event Horizon of Sagittarius A*,” Astrophys. J.  701 (2009) 1357 [arXiv:0903.1105 [astro-ph.HE]].
  • (2) R. M. Wald, “General Relativity,” Chicago, Usa: Univ. Pr. ( 1984) 491p.
  • (3) R. P.  Kerr, “Gravitational field of a spinning mass as an example of algebraically special metrics” Phys. Rev. Lett. 11 (1963) 237-238. E. T. Newman, R. Couch, K. Chinnapared, A. Exton, A. Prakash and R. Torrence, “Metric of a Rotating, Charged Mass,” J. Math. Phys.  6 (1965) 918.
  • (4) S. W. Hawking and G. F. R. Ellis, Cambridge University Press, Cambridge, 1973
  • (5) J. D. Bekenstein, “Black holes and entropy,” Phys. Rev.  D 7 (1973) 2333.
  • (6) S. W. Hawking, “Particle Creation By Black Holes,” Commun. Math. Phys.  43 (1975) 199 [Erratum-ibid.  46 (1976) 206].
  • (7) A. Strominger and C. Vafa, “Microscopic Origin of the Bekenstein-Hawking Entropy,” Phys. Lett.  B 379 (1996) 99 [arXiv:hep-th/9601029].
  • (8) S. Carlip, “Entropy from conformal field theory at Killing horizons,” Class. Quant. Grav.  16 (1999) 3327 [arXiv:gr-qc/9906126]. S. Carlip, “Black hole entropy from conformal field theory in any dimension,” Phys. Rev. Lett.  82 (1999) 2828 [arXiv:hep-th/9812013].
  • (9) C. Rovelli, “Black hole entropy from loop quantum gravity,” Phys. Rev. Lett.  77 (1996) 3288 [arXiv:gr-qc/9603063].
  • (10) Ashtekar A, Baez B and Krasnov K. “Quantum geometry of isolated horizons and black hole entropy” Adv. Theor. Math. Phys.  4 1-94.
  • (11) J. C. Baez, “An introduction to spin foam models of BF theory and quantum gravity,” Lect. Notes Phys.  543 (2000) 25 [arXiv:gr-qc/9905087].
  • (12) J. Engle, A. Perez and K. Noui, “Black hole entropy and SU(2) Chern-Simons theory,” arXiv:0905.3168 [gr-qc].
  • (13) K. V. Krasnov, “On statistical mechanics of gravitational systems,” Gen. Rel. Grav.  30 (1998) 53 [arXiv:gr-qc/9605047]. L. Smolin, “Linking topological quantum field theory and nonperturbative quantum gravity,” J. Math. Phys.  36 (1995) 6417.
  • (14) A. Ashtekar, A. Corichi and K. Krasnov, “Isolated horizons: The classical phase space,” Adv. Theor. Math. Phys.  3 (2000) 419 [arXiv:gr-qc/9905089].
  • (15) A. Ashtekar, J. Engle and C. Van Den Broeck, “Quantum horizons and black hole entropy: Inclusion of distortion and rotation,” Class. Quant. Grav.  22 (2005) L27 [arXiv:gr-qc/0412003]. C.  Beetle, J.  Engle ”Quantization of generic isolated horizons in loop quantum gravity.” To appear.
  • (16) A. Ashtekar, C. Beetle and J. Lewandowski, “Geometry of Generic Isolated Horizons,” Class. Quant. Grav.  19 (2002) 1195 [arXiv:gr-qc/0111067].
  • (17) Ashtekar A, Beetle C, Dreyer O, Fairhurst S, Krishnan B, Lewandowski J and Wiśniewski J 2000 Generic Isolated Horizons and their applications Phys. Rev. Lett. 85 3564-3567
  • (18) J. Lewandowski, “Spacetimes Admitting Isolated Horizons,” Class. Quant. Grav.  17 (2000) L53.
  • (19) Ashtekar A, Fairhurst S and Krishnan B 2000 Isolated horizons: Hamiltonian evolution and the first law Phys.  Rev. D 62 104025
  • (20) Ashtekar A, Beetle C and Lewandowski J 2002 Geometry of generic isolated horizons Class. Quantum Grav. 19 1195-1225
  • (21) Ashtekar A, Beetle C and Lewandowski J 2001 Mechanics of rotating isolated horizons Phys. Rev.  D 64 044016
  • (22) Lewandowski J 2000 Space-times admitting isolated horizons Class.Quant.Grav. 17 L53-L59
  • (23) T. Thiemann, “Modern canonical quantum general relativity,” Cambridge, UK: Cambridge Univ. Pr. (2007) 819 p
  • (24) C. Crnkovic and E. Witten, in ‘Three hundred years of gravitation’; ed. S. Hawking, W. Israel. J. Lee and R. M. Wald, “Local symmetries and constraints,” J. Math. Phys. 31 (1990) 725. A. Ashtekar, L. Bombelli and O. Reula, in ’200 Years After Lagrange’, Ed. by M. Francaviglia, D. Holm.
  • (25) E. Witten, “Quantum field theory and the Jones polynomial,” Commun. Math. Phys.  121 (1989) 351.
  • (26) A. Ashtekar, “Lectures on nonperturbative canonical gravity,” Singapore, Singapore: World Scientific (1991) 334 p. (Advanced series in astrophysics and cosmology, 6)
  • (27) R. K. Kaul and P. Majumdar, “Quantum black hole entropy,” Phys. Lett.  B 439 (1998) 267. R. K. Kaul and P. Majumdar, “Logarithmic correction to the Bekenstein-Hawking entropy,” Phys. Rev. Lett.  84 (2000) 5255.
  • (28) E. Livine and D. Terno, “Quantum black holes: Entropy and entanglement on the horizon,” Nucl. Phs.  B 741 (2006) 131.
  • (29) A. Ghosh and P. Mitra, “A bound on the log correction to the black hole area law,” Phys. Rev.  D 71 (2005) 027502. G. Gour, “Algebraic approach to quantum black holes: Logarithmic corrections to black hole entropy,” Phys. Rev.  D 66 (2002) 104022.
  • (30) I. Agullo, J. F. Barbero G., J. Diaz-Polo, E. Fernandez-Borja and E. J. S. Villasenor, “Black hole state counting in LQG: A number theoretical approach,” Phys. Rev. Lett.  100 (2008) 211301. J. F. Barbero G. and E. J. S. Villasenor, “Generating functions for black hole entropy in Loop Quantum Gravity,” Phys. Rev.  D 77 (2008) 121502. “On the computation of black hole entropy in loop quantum gravity,” Class. Quant. Grav.  26 (2009) 035017
  • (31) I. Agullo, J. F. B. G., E. F. Borja, J. Diaz-Polo and E. J. S. Villasenor, “The combinatorics of the SU(2) black hole entropy in loop quantum gravity,” arXiv:0906.4529 [gr-qc].
  • (32) S. Carlip, “Logarithmic corrections to black hole entropy from the Cardy formula,” Class. Quant. Grav.  17 (2000) 4175 [arXiv:gr-qc/0005017].
  • (33) D. J. Rezende and A. Perez, Phys. Rev.  D 79 (2009) 064026 [arXiv:0902.3416 [gr-qc]].
  • (34) L. Liu, M. Montesinos and A. Perez, arXiv:0906.4524 [gr-qc].
  • (35) A. Perez, “Black Hole Entropy and SU(2) Chern Simons Theory”, International Loop quantum Gravity Seminar, http://relativity.phys.lsu.edu/ilqgs/perez050410.pdf May 4, 2010.
  • (36) A. Corichi and E. Wilson-Ewing, “Surface terms, Asymptotics and Thermodynamics of the Holst Action,” arXiv:1005.3298 [gr-qc].
  • (37) A. Ashtekar, A. Corichi and J. A. Zapata, “Quantum theory of geometry. III: Non-commutativity of Riemannian structures,” Class. Quant. Grav.  15 (1998) 2955 [arXiv:gr-qc/9806041].
  • (38) J. Lewandowski, A. Okolow, H Sahlmann, and Thiemann T. Uniqueness of the diffeomorphism invariant state on the quantum holonomy-flux algebra. 2004. C. Fleischhack, “Representations of the Weyl Algebra in Quantum Geometry,” Commun. Math. Phys.  285 (2009) 67 [arXiv:math-ph/0407006].
  • (39) S. Chandrasekhar, “The mathematical theory of black holes,” Oxford, UK: Clarendon (1992) 646 p.
  • (40) A. Perez, D. Pranzetti, in preparation.
  • (41) A. Perez, “SU(2)SU(2) Chern-Simons theory and black hole entropy”, Plenary talk presented at LOOPS09, Beiging, China, August 2009.
  • (42) K. Noui, J. Engle, A. Perez, ”The SU(2) Black Hole entropy revisited”, in preparation.
  • (43) R. K. Kaul and P. Majumdar, “Schwarzschild horizon dynamics and SU(2) Chern-Simons theory,” arXiv:1004.5487 [gr-qc].