Abstract
The entanglement entropy is a fundamental quantity, which characterizes the correlations between subsystems in a larger quantummechanical system. For two subsystems separated by a surface the entanglement entropy is proportional to the area of the surface and depends on the UV cutoff, which regulates the shortdistance correlations. The geometrical nature of entanglemententropy calculation is particularly intriguing when applied to black holes when the entangling surface is the blackhole horizon. I review a variety of aspects of this calculation: the useful mathematical tools such as the geometry of spaces with conical singularities and the heat kernel method, the UV divergences in the entropy and their renormalization, the logarithmic terms in the entanglement entropy in four and six dimensions and their relation to the conformal anomalies. The focus in the review is on the systematic use of the conical singularity method. The relations to other known approaches such as ’t Hooft’s brickwall model and the Euclidean path integral in the optical metric are discussed in detail. The puzzling behavior of the entanglement entropy due to fields, which nonminimally couple to gravity, is emphasized. The holographic description of the entanglement entropy of the blackhole horizon is illustrated on the two and fourdimensional examples. Finally, I examine the possibility to interpret the BekensteinHawking entropy entirely as the entanglement entropy.
Introduction
One of the mysteries in modern physics is why black holes have an entropy. This entropy, known as the BekensteinHawking entropy, was first introduced by Bekenstein [18, 19, 20] as a rather useful analogy. Soon after that, this idea was put on a firm ground by Hawking [128] who showed that black holes thermally radiate and calculated the blackhole temperature. The main feature of the BekensteinHawking entropy is its proportionality to the area of the blackhole horizon. This property makes it rather different from the usual entropy, for example the entropy of a thermal gas in a box, which is proportional to the volume. In 1986 Bombelli, Koul, Lee and Sorkin [23] published a paper in which they considered the reduced density matrix, obtained by tracing over the degrees of freedom of a quantum field that are inside the horizon. This procedure appears to be very natural for black holes, since the black hole horizon plays the role of a causal boundary, which does not allow anyone outside the black hole to have access to the events, which take place inside the horizon. Another attempt to understand the entropy of black holes was made by ’t Hooft in 1985 [214]. His idea was to calculate the entropy of the thermal gas of Hawking particles, which propagate just outside the horizon. This calculation has uncovered two remarkable features: the entropy does turn out to be proportional to the horizon area, however, in order to regularize the density of states very close to the horizon, it was necessary to introduce the brick wall, a boundary, which is placed at a small distance from the actual horizon. This small distance plays the role of a regulator in the ’t Hooft’s calculation. Thus, the first indications that entropy may grow as area were found.
An important step in the development of these ideas was made in 1993 when a paper of Srednicki [208] appeared. In this very inspiring paper Srednicki calculated the reduced density and the corresponding entropy directly in flat spacetime by tracing over the degrees of freedom residing inside an imaginary surface. The entropy defined in this calculation has became known as the entanglement entropy. Sometimes the term geometric entropy is used as well. The entanglement entropy, as was shown by Srednicki, is proportional to the area of the entangling surface. This fact is naturally explained by observing that the entanglement entropy is nonvanishing due to the shortdistance correlations present in the system. Thus, only modes, which are located in a small region close to the surface, contribute to the entropy. By virtue of this fact, one finds that the size of this region plays the role of the UV regulator so that the entanglement entropy is a UV sensitive quantity. A surprising feature of Srednicki’s calculation is that no black hole is actually needed: the entanglement entropy of a quantum field in flat spacetime already establishes the area law. In an independent paper, Frolov and Novikov [99] applied a similar approach directly to a black hole. These results have sparked interest in the entanglement entropy. In particular, it was realized that the brickwall model of ’t Hooft studies a similar entropy and that the two entropies are in fact related. On the technical side of the problem, a very efficient method was developed to calculate the entanglement entropy. This method, first considered by Susskind [211], is based on a simple replica trick, in which one first introduces a small conical singularity at the entangling surface, evaluates the effective action of a quantum field on the background of the metric with a conical singularity and then differentiates the action with respect to the deficit angle. By means of this method one has developed a systematic calculation of the UV divergent terms in the geometric entropy of black holes, revealing the covariant structure of the divergences [33, 197, 111]. In particular, the logarithmic UV divergent terms in the entropy were found [196]. The other aspect, which was widely discussed in the literature, is whether the UV divergence in the entanglement entropy could be properly renormalized. It was suggested by Susskind and Uglum [213] that the standard renormalization of Newton’s constant makes the entropy finite, provided one considers the entanglement entropy as a quantum contribution to the BekensteinHawking entropy. However, this proposal did not answer the question of whether the BekensteinHawking entropy itself can be considered as an entropy of entanglement. It was proposed by Jacobson [141] that, in models in which Newton’s constant is induced in the spirit of Sakharov’s ideas, the BekensteinHawking entropy would also be properly induced. A concrete model to test this idea was considered in [97].
Unfortunately, in the 1990s, the study of entanglement entropy could not compete with the booming success of the string theory (based on Dbranes) calculations of blackhole entropy [209]. The second wave of interest in entanglement entropy started in 2003 with work studying the entropy in condensed matter systems and in lattice models. These studies revealed the universality of the approach based on the replica trick and the efficiency of the conformal symmetry to compute the entropy in two dimensions. Black holes again came into the focus of study in 2006 after work of Ryu and Takayanagi [189] where a holographic interpretation of the entanglement entropy was proposed. In this proposal, in the frame of the AdS/CFT correspondence, the entanglement entropy, defined on a boundary of antide Sitter, is related to the area of a certain minimal surface in the bulk of the antide Sitter spacetime. This proposal opened interesting possibilities for computing, in a purely geometrical way, the entropy and for addressing in a new setting the question of the statistical interpretation of the BekensteinHawking entropy.
The progress made in recent years and the intensity of the ongoing research indicate that entanglement entropy is a very promising direction, which, in the coming years, may lead to a breakthrough in our understanding of black holes and quantum gravity. A number of very nice reviews appeared in recent years that address the role of entanglement entropy for black holes [21, 90, 146, 54]; review the calculation of entanglement entropy in quantum field theory in flat spacetime [81, 37] and the role of the conformal symmetry [31]; and focus on the holographic aspects of the entanglement entropy [185, 11]. In the present review I build on these works and focus on the study of entanglement entropy as applied to black holes. The goal of this review is to collect a complete variety of results and present them in a systematic and selfconsistent way without neglecting either technical or principal aspects of the problem.
Entanglement Entropy in Minkowski Spacetime
Definition
Consider a pure vacuum state ψ >of a quantum system defined inside a spacelike region \({\mathcal O}\) and suppose that the degrees of freedom in the system can be considered as located inside certain subregions of \({\mathcal O}\). A simple example of this sort is a system of coupled oscillators placed in the sites of a spacelike lattice. Then, for an arbitrary imaginary surface Σ, which separates the region \({\mathcal O}\) into two complementary subregions A and B, the system in question can be represented as a union of two subsystems. The wave function of the global system is given by a linear combination of the product of quantum states of each subsystem, ψ >= ∑_{ i,a } ψ_{ ia } A > _{ i }B >_{ a }. The states A >_{ i } are formed by the degrees of freedom localized in the region A, while the states B >_{ a } are formed by those, which are defined in region B. The density matrix that corresponds to a pure quantum state ψ >
has zero entropy. By tracing over the degrees of freedom in region A we obtain a density matrix
with elements (ρ_{ B })_{ ab } = (ψψ^{†})_{ ab }. The statistical entropy, defined for this density matrix by the standard formula
is by definition the entanglement entropy associated with the surface Σ. We could have traced over the degrees of freedom located in region B and formed the density matrix (ρ_{ A })_{ ij } = (ψ^{T}ψ*)_{ ij }. It is clear that^{Footnote 1}
for any integer k. Thus, we conclude that the entropy (3) is the same for both density matrices ρ_{ A } and ρ_{ B },
This property indicates that the entanglement entropy for a system in a pure quantum state is not an extensive quantity. In particular, it does not depend on the size of each region A or B and thus is only determined by the geometry of Σ.
Shortdistance correlations
On the other hand, if the entropy (3) is nonvanishing, this shows that in the global system there exist correlations across the surface Σ between modes, which reside on different sides of the surface. In this review we shall consider the case in which the system in question is a quantum field. The shortdistance correlations that exist in this system have two important consequences:

the entanglement entropy becomes dependent on the UV cutoff ϵ, which regularizes the shortdistance (or the largemomentum) behavior of the field system

to leading order in ϵ^{−1} the entanglement entropy is proportional to the area of the surface Σ
For a free massless scalar field the 2point correlation function in d spacetime dimensions has the standard form
where \({\Omega _d} = {{\Gamma ({{d  2} \over 2})} \over {4{\pi ^{d/2}}}}\). Correspondingly, the typical behavior of the entanglement entropy in d dimensions is
where the exact prefactor depends on the regularization scheme. Although the similarity between (5) and (6) illustrates well the fieldtheoretical origin of the entanglement entropy, the exact relation between the shortdistance behavior of 2point correlation functions in the field theory and the UV divergence of the entropy is more subtle, as we shall discuss later in the paper.
Thermal entropy
Instead of a pure state one could have started with a mixed thermal state at temperature T with density matrix \({\rho _0}(A,B) = {e^{ {T^{ 1}}H(A,B)}}\), where H(A, B) is the Hamiltonian of the global system. In this case the relation (4) is no more valid and the entropy depends on the size of the total system as well as on the size of each subsystem. By rather general arguments, in the limit of large volume the reduced density matrix approaches the thermal density matrix. So that in this limit the entanglement entropy (3) reproduces the thermal entropy. For further references we give here the expression
for the thermal entropy of a massless field residing inside a spatial (d − 1)volume V_{d−1} at temperature T.
Entropy of a system of finite size at finite temperature
In a more general situation one starts with a system of finite size L in a mixed thermal state at temperature T. This system is divided by the entangling surface Σ in two subsystems of characteristic size l. Then, the entanglement entropy is a function of several parameters (if the field in question is massive then mass m should be added to the parameters on which the entropy should depend)
where ϵ is a UV cutoff. Clearly, the entanglement entropy in this general case is due to a combination of different factors: the entanglement between two subsystems and the thermal nature of the initial mixed state. In d dimensions even for simple geometries this function of 4 variables is not known explicitly. However, in two spacetime dimensions, in some particular cases, the explicit form of this function is known.
Entropy in (1+1)dimensional spacetime
The state of a quantum field in two dimensions is defined on a union of intersecting intervals A ∪ B. The 2point correlation functions behave logarithmically in the limit of coincident points. Correspondingly, the leading UV divergence of the entanglement entropy in two dimensions is logarithmic. For example, for a 2D massless conformal field theory, characterized by a central charge c, the entropy is [208, 33, 133]
where n is the number of intersections of intervals A and B, where the subsystems are defined, l_{ A } (l_{ B }) is the length of the interval A (B). The second term in Eq. 9) is a UV finite term. In some cases the conformal symmetry in two dimensions can be used to calculate not only the UV divergent term in the entanglement entropy but also the UV finite term, thus obtaining the complete answer for the entropy, as was shown by Holzhey, Larsen and Wilczek [133] (see [161, 29] for more recent developments). There are two different limiting cases when the conformal symmetry is helpful. In the first case, one considers a pure state of the conformal field theory on a circle of circumference L, the subsystem is defined on a segment of size l of the circle. In the second situation, the system is defined on an infinite line, the subsystem lives on an interval of length l of the line and the global system is in a thermal mixed state with temperature T. In Euclidean signature both geometries represent a cylinder. For a thermal state the compact direction on the cylinder corresponds to Euclidean time τ compactified to form a circle of circumference β = 1/T. In both cases the cylinder can be further conformally mapped to a plane. The invariance of the entanglement entropy under conformal transformation can be used to obtain
in the case of a pure state on a circle and
for a thermal mixed state on an infinite line. In the limit of large l the entropy (11) approaches
where the first term represents the entropy of the thermal gas (7) in a cavity of size l, while the second term represent the purely entanglement contribution (note that the intersection of A and B contains two points in this case so that n = 2). The third term is an intermediate term due to the interaction of both factors, thermality and entanglement. This example clearly shows that for a generic thermal state the entanglement entropy is due to the combination of two factors: the entanglement between two subsystems and the thermal nature of the mixed state of the global system.
The Euclidean path integral representation and the replica method
A technical method very useful for the calculation of the entanglement entropy in a field theory is the the replica trick, see [33]. Here we illustrate this method for a field theory described by a secondorder Laplacetype operator. One considers a quantum field ψ(X) in a ddimensional spacetime and chooses the Cartesian coordinates X^{μ} = (τ,x,z^{i}, i = 1, .., d−2), where τ is Euclidean time, such that the surface Σ is defined by the condition x = 0 and (z^{i}, i = 1, .., d−2) are the coordinates on Σ. In the subspace (τ,x) it will be convenient to choose the polar coordinate system τ= r sin(ϕ) and x = r cos(ϕ), where the angular coordinate ϕ varies between 0 and 2π. We note that if the field theory in question is relativistic, then the field operator is invariant under the shifts ϕ → ϕ + w, where w is an arbitrary constant.
One first defines the vacuum state of the quantum field in question by the path integral over a half of the total Euclidean spacetime defined as ρ ≤ 0 such that the quantum field satisfies the fixed boundary condition ψ(τ = 0, x, z) = ψ_{0}(x, z) on the boundary of the halfspace,
where W[ψ] is the action of the field. The surface Σ in our case is a plane and the Cartesian coordinate x is orthogonal to Σ. The codimension 2 surface Σ defined by the conditions x = 0 and τ = 0 naturally separates the hypersurface τ = 0 into two parts: x < 0 and x > 0. These are the two subregions A and B discussed in Section 2.1.
The boundary data ψ(x, z) is also separated into ψ_{−}(x, z) = ψ_{0}(x, z), x < 0 and ψ_{+} = ψ_{0}(x, z), x > 0. By tracing over ψ_{−}(x, z) one defines a reduced density matrix
where the path integral goes over fields defined on the whole Euclidean spacetime except a cut (τ = 0, x > 0). In the path integral the field ψ(X) takes the boundary value \(\psi _ + ^2\) above the cut and \(\psi _ + ^1\) below the cut. The trace of the nth power of the density matrix (14) is then given by the Euclidean path integral over fields defined on an nsheeted covering of the cut spacetime. In the polar coordinates (r, ϕ) the cut corresponds to values ϕ = 2πk, k = 1, 2, .., n. When one passes across the cut from one sheet to another, the fields are glued analytically. Geometrically this nfold space is a flat cone C_{ n } with angle deficit 2π(1 − n) at the surface Σ. Thus, we have
where Z[C_{ n }] is the Euclidean path integral over the nfold cover of the Euclidean space, i.e., over the cone C_{ n }. Assuming that in Eq. (15) one can analytically continue to noninteger values of n, one observes that
where \(\hat \rho = \rho/{\rm{Tr}}\rho\) is the renormalized matrix density. Introduce the effective action W(α) = − ln Z(α), where Z(α) = Z[C_{ α }] is the partition function of the field system in question on a Euclidean space with conical singularity at the surface Σ. In the polar coordinates (r, ϕ) the conical space C_{ α } is defined by making the coordinate ϕ periodic with period 2πα, where (1 − α) is very small. The invariance under the abelian isometry ϕ → ϕ + w helps to construct without any problem the correlation functions with the required periodicity 2πα starting from the 2πperiodic correlation functions. The analytic continuation of Trρ^{α} to α different from 1 in the relativistic case is naturally provided by the path integral Z(α) over the conical space C_{ α }. The entropy is then calculated by the replica trick
One of the advantages of this method is that we do not need to care about the normalization of the reduced density matrix and can deal with a matrix, which is not properly normalized.
Uniqueness of analytic continuation
The uniqueness of the analytic continuation of Trρ^{n} to noninteger n may not seem obvious, especially if the field system in question is not relativistic so that there is no isometry in the polar angle ϕ, which would allow us, without any trouble, to glue together pieces of the Euclidean space to form a path integral over a conical space C_{ α }. However, some arguments can be given that the analytic continuation to noninteger n is in fact unique.
Consider a renormalized density matrix \(\hat \rho = {\rho \over {{\rm{Tr}}\rho}}\). The eigenvalues of \(\hat{\rho}\) lie in the interval 0 < λ < 1. If this matrix were a finite matrix we could use the triangle inequality to show that
For infinitesize matrices the trace is usually infinite so that a regularization is needed. Suppose that ϵ is the regularization parameter and Tr_{ ϵ } is the regularized trace. Then
Thus \({\rm{Tr}}{\hat \rho ^\alpha}\) is a bounded function in the complex halfplane, Re(α) > 1. Now suppose that we know that Tr_{ ϵ }ρ^{α}_{α = n} = Z_{0}(n) for integer values of α = n, n = 1, 2, 3, … Then, in the region Re(α) > 1, we can represent Z(α) = Tr_{ ϵ }ρ^{α} in the form
where the function g(α) is analytic (for Re(α) > 1). Since by condition (17) the function Z(α) is bounded, we obtain that, in order to compensate for the growth of the sine in Eq. (18) for complex values of α, the function g(α) should satisfy the condition
By Carlson’s theorem [36] an analytic function, which is bounded in the region Re(α) > 1 and which satisfies condition (19), vanishes identically. Thus, we conclude that g(α) ≡ 0 and there is only one analytic continuation to noninteger n, namely the one given by function Z_{0}(α).
Heat kernel and the Sommerfeld formula
Consider for concreteness a quantum bosonic field described by a field operator \({\mathcal D}\) so that the partition function is \(Z = {\det ^{1/2}}{\mathcal D}\). Then, the effective action defined as
where parameter ϵ is a UV cutoff, is expressed in terms of the trace of the heat kernel \(K(s,\, X,\, X{\prime}) = < X\vert {e^{ s{\mathcal D}}}\vert X{\prime} >\). The latter is defined as a solution to the heat equation
In order to calculate the effective action W(α) we use the heat kernel method. In the context of manifolds with conical singularities this method was developed in great detail in [69, 101]. In the Lorentz invariant case the invariance under the abelian symmetry ϕ → ϕ + w plays an important role. The heat kernel (s, ϕ, ϕ′) (where we omit the coordinates other than the angle ϕ) on regular flat space then depends on the difference (ϕ − ϕ′). This function is 2π periodic with respect to (ϕ − ϕ′). The heat kernel Kα(s, ϕ, ϕ′) on a space with a conical singularity is supposed to be 2πα periodic. It is constructed from the 2π periodic quantity by applying the Sommerfeld formula [207]
That this quantity still satisfies the heat kernel equation is a consequence of the invariance under the abelian isometry ϕ → ϕ + w. The contour Γ consists of two vertical lines, going from −π+ i∞) to (−π − i∞) and from (π − i∞) to (π − + i∞) and intersecting the real axis between the poles of the cot \({\rm{cot}} {w \over {2\alpha}}:  2\pi \alpha, \, 0\) and 0, +2πα, respectively. For α = 1 the integrand in Eq. (22) is a 2πperiodic function and the contributions of these two vertical lines cancel each other. Thus, for a small angle deficit the contribution of the integral in Eq. (22) is proportional to (1 − α).
An explicit calculation
Consider an infinite (d − 2)plane in ddimensional spacetime. The calculation of the entanglement entropy for this plane can be done explicitly by means of the heat kernel method. In flat spacetime, if the operator \({\mathcal D}\) is the Laplace operator,
one can use the Fourier transform in order to solve the heat equation. In d spacetime dimensions one has
Putting z^{i} = z′^{i}, i = 1, .., d − 2 and choosing in the polar coordinate system (r, ϕ), that ϕ = ϕ′ + w we have that \({p_\mu}{(X  X{\prime})^\mu} = 2pr\sin {w \over 2}\cos \theta\), where p^{2} = p^{μ}p_{ μ } and θ is the angle between the dvectors p^{μ} and (X^{μ} − X^{′μ}). The radial momentum p and angle θ, together with the other (d−2) angles form a spherical coordinate system in the space of momenta p^{μ}. Thus, one has for the integration measure \(\int {{d^d}} p = {\Omega _{d  2}}\int\nolimits_0^\infty dp\, {p^{d  1}}\int\nolimits_0^\pi d \theta \, {\sin ^{d  2}}\theta\), where \(\Omega_{d2}={2\, \pi^{(d1)/2} \over\Gamma((d1)/2)}\) is the area of a unit radius sphere in d−1 dimensions. Performing the integration in Eq. (23) in this coordinate system we find
For the trace one finds
where A(Σ) = ∫ d^{d−2} z is the area of the surface Σ. One uses the integral \(\int\nolimits_0^\infty dx x^{1\nu} J_\nu(x)={2^{1\nu}\over \Gamma(\nu)}\) for the derivation of Eq. (25). The integral over the contour Γ in the Sommerfeld formula (22) is calculated via residues ([69, 101])
Collecting everything together one finds that in flat Minkowski spacetime
where \(V=\int d\tau d^{d1}x\) is the volume of spacetime and A(Σ) = ∫ d^{d−2}x is the area of the surface Σ. Substituting Eq. (27) into Eq. (20) we obtain that the effective action contains two terms. The one proportional to the volume V reproduces the vacuum energy in the effective action. The second term proportional to the area A(Σ) is responsible for the entropy. Applying formula (16) we obtain the entanglement entropy
of an infinite plane Σ in d spacetime dimensions. Since any surface, locally, looks like a plane, and a curved spacetime, locally, is approximated by Minkowski space, this result gives the leading contribution to the entanglement entropy of any surface Σ in flat or curved spacetime.
Entropy of massive fields
The heat kernel of a massive field described by the wave operator \({\mathcal D}=\nabla^2 +m^2\) is expressed in terms of the heat kernel of a massless field,
Thus, one finds
where the trace of the heat kernel for vanishing mass is given by Eq. (27). Therefore, the entanglement entropy of a massive field is
In particular, if d = 4, one finds that
The logarithmic term in the entropy that is due to the mass of the field appears in any even dimension d. The presence of a UV finite term proportional to the (d − 2)th power of mass is the other general feature of (30), (31).
An expression in terms of the determinant of the Laplacian on the surface
Even though the entanglement entropy is determined by the geometry of the surface Σ, in general, this can be not only its intrinsic geometry but also how the surface is embedded in the larger spacetime. The embedding is determined by the extrinsic curvature. The curvature of the larger spacetime enters through the GaussCadazzi relations. But in some particularly simple cases the entropy can be given a purely intrinsic interpretation. To see this for the case when Σ is a plane we note that the entropy (28) or (30) originates from the surface term in the trace of the heat kernel (27) (or (31)). To leading order in (1 − α), the surface term in the case of a massive scalar field is
where
can be interpreted as the trace of the heat kernel of operator −Δ(Σ) + m^{2}, where Δ(Σ) is the intrinsic Laplace operator defined on the (d − 2)plane Σ. The determinant of the operator −Δ(Σ) + m^{2} is determined by
Thus, we obtain an interesting expression for the entanglement entropy
in terms of geometric objects defined intrinsically on the surface Σ. A similar expression in the case of an ultraextreme black hole was obtained in [172] and for a generic black hole with horizon approximated by a plane was obtained in [89].
Entropy in theories with a modified propagator
In certain physicallyinteresting situations the propagator of a quantum field is different from the standard 1/(p^{2} + m^{2}) and is described by some function as 1/F(p^{2}). The quantum field in question then satisfies a modified Lorentz invariant field equation
Theories of this type naturally arise in models with extra dimensions. The deviations from the standard form of propagator may be both in the UV regime (large values of p) or in the IR regime (small values of p). If the function F(p^{2}) for large values of p grows faster than p^{2} this theory is characterized by improved UV behavior.
The calculation of the entanglement entropy performed in Section 2.9 can be generalized to include theories with operator (33). This example is instructive since, in particular, it illuminates the exact relation between the structure of 2point function (the Green’s function in the case of free fields) and the entanglement entropy [183].
In d spacetime dimensions one has
Note that we consider Euclidean theory so that p^{2} ≥ 0. The Green’s function
is a solution to the field equation with a deltalike source
and can be expressed in terms of the heat kernel as follows
Obviously, the Green’s function can be represented in terms of the Fourier transform in a manner similar to Eq. (34),
The calculation of the trace of the heat kernel for operator (33) on a space with a conical singularity goes along the same lines as in Section 2.9. This was performed in [184] and the result is
where the functions P_{ n }(s) are defined as
The entanglement entropy takes the form (we remind the reader that for simplicity we take the surface Σ to be a (d − 2)dimensional plane) [184]
It is important to note that [184, 183]

(i)
the area law in the entanglement entropy is universal and is valid for any function F(p^{2});

(ii)
the entanglement entropy is UV divergent independently of the function F(p^{2}), with the degree of divergence depending on the particular function F(p^{2});

(iii)
in the coincidence limit, X = Y, the Green’s function (38)
$$G(X,X) = {2 \over {\Gamma ({d \over 2})}}\,{1 \over {{{(4\pi)}^{{d \over 2}}}}}\int\nolimits_0^\infty d p{p^{d  1}}G({p^2})$$(42)may take a finite value if G(p^{2}) = 1/F(p^{2}) is decaying faster than 1/p^{d}. However, even for this function F(p^{2}), the entanglement entropy is UV divergent.
As an example, consider a function, which grows for large values of p as F(p^{2}) ∼ p^{2k}. The 2point correlation function in this theory behaves as
and for k > d/2 it is regular in the coincidence limit. On the other hand, the entanglement entropy scales as
and remains divergent for any positive value of k. Comparison of Eqs. (43) and (44) shows that only for k = 1 (the standard form of the wave operator and the propagator) the shortdistance behavior of the 2point function is similar to the UV divergence of the entanglement entropy.
Entanglement entropy in nonLorentz invariant theories
NonLorentz invariant theories are characterized by a modified dispersion relation, ω^{2} + F(p^{2}) = 0, between the energy ω and the 3momentum p. These theories can be described by a wave operator of the following type
where \(\Delta_x=\sum\nolimits_i^{d1}\partial_i^2\) is the spatial Laplace operator. Clearly, the symmetry with respect to the Lorentz boosts is broken in operator (45) if F(q) ≠ q.
As in the Lorentz invariant case to compute the entanglement entropy associated with a surface Σ we choose (d − 1) spatial coordinates {x^{i}, i = 1,.., d − 1} = {x, z^{a}, a = 1, .., d − 2}, where x is the coordinate orthogonal to the surface Σ and z^{a} are the coordinates on the surface Σ. Then, after going to Euclidean time τ= it, we switch to the polar coordinates, τ= r sin(ϕ), x = cos(ϕ). In the Lorentz invariant case the conical space, which is needed for calculation of the entanglement entropy, is obtained by making the angular coordinate ϕ periodic with period 2πα by applying the Sommerfeld formula (22) to the heat kernel. If Lorentz invariance is broken, as it is for the operator (45), there are certain difficulties in applying the method of the conical singularity when one computes the entanglement entropy. The difficulties come from the fact that the wave operator \({\mathcal D}\), if written in terms of the polar coordinates r and ϕ, becomes an explicit function of the angular coordinate ϕ. As a result of this, the operator \({\mathcal D}\) is not invariant under shifts of ϕ to arbitrary ϕ + w. Only shifts with w = 2πn, where n is an integer are allowed. Thus, in this case one cannot apply the Sommerfeld formula since it explicitly uses the symmetry of the differential operator under shifts of angle ϕ. On the other hand, a conical space with angle deficit 2π(1 − n) is exactly what we need to compute Trρ^{n} for the reduced density matrix. In [184], by using some scaling arguments it was shown that the trace of the heat kernel \(K(s) = {e^{ s{\mathcal D}}}\) on a conical space with 2πn periodicity, is
where nTrK_{n=1}(s) is the bulk contribution. By the arguments presented in Section 2.7 there is a unique analytic extension of this formula to noninteger n. A simple comparison with the surface term in the heat kernel of the Lorentz invariant operator, which was obtained in Section 2.12, shows that the surface terms of the two kernels are identical. Thus, we conclude that the entanglement entropy is given by the same formula
where P_{ n }(s) is defined in Eq. (40), as in the Lorentz invariant case (41). A similar property of the entanglement entropy was observed for a nonrelativistic theory described by the Schrödinger operator [205] (see also [59] for a holographic derivation). For polynomial operators, F(q) ∼ q^{k}, some scaling arguments can be used [205] to get the form of the entropy that follows from Eq. (47).
In the rest of the review we shall mostly focus on the study of Lorentz invariant theories, with field operator quadratic in derivatives, of the Laplace type, \({\mathcal D} =  ({\nabla ^2} + X)\).
Arbitrary surface in curved spacetime: general structure of UV divergences
The definition of the entanglement entropy and the procedure for its calculation generalize to curved spacetime. The surface Σ can then be any smooth closed codimension two surface^{Footnote 2}, which divides the space into two subregions. In Section 3 we will consider in detail the case where this surface is a blackhole horizon. Before proceeding to the blackhole case we would like to specify the general structure of UV divergent terms in the entanglement entropy. In ddimensional curved spacetime, entanglement entropy is presented in the form of a Laurent series with respect to the UV cutoff ϵ (for d = 4 see [204])
where s_{d−2} is proportional to the area of the surface Σ. All other terms in the expansion (48) can be presented as integrals over Σ of local quantities constructed in terms of the Riemann curvature of the spacetime and the extrinsic curvature of the surface Σ. Of course, the intrinsic curvature of the surface Σ can be expressed in terms of \({\mathcal R}\) and k using the GaussCodazzi equations. Since nothing should depend on the direction of vectors normal to Σ, the integrands in expansion (48) should be even powers of extrinsic curvature. The general form of the s_{d−2−2n} term can be symbolically presented in the form
where \({\mathcal R}\) stands for components of the Riemann tensor and their projections onto the subspace orthogonal to Σ and k labels the components of the extrinsic curvature. Thus, since the integrands are even in derivatives, only terms ϵ^{d−2n−2}, n = 0, 1, 2, … appear in Eq. (48). If d is even, then there also may appear a logarithmic term s_{0}. The term s(g) in Eq. (48) is a UV finite term, which may also depend on the geometry of the surface Σ, as well as on the geometry of the spacetime itself.
Entanglement Entropy of NonDegenerate Killing Horizons
The geometric setting of blackhole spacetimes
The notion of entanglement entropy is naturally applicable to a black hole. In fact, probably the only way to separate a system into two subsystems is to place one of them inside a blackhole horizon. The important feature that, in fact, defines the black hole is the existence of a horizon. Many useful definitions of a horizon are known. In the present paper we shall consider only the case of the eternal black holes for which different definitions of the horizon coincide. The corresponding spacetime then admits a maximal analytic extension, which we shall use in our construction. The simplest example is the Schwarzschild black hole, the maximal extension of which is demonstrated on the wellknown Penrose diagram. The horizon of the Schwarzschild black hole is an example of a Killing horizon. The spacetime in this case possesses a global Killing vector, ξ_{ t } = ∂_{ t }, which generates the time translations. The Killing horizon is defined as a null hypersurface on which the Killing vector ξ_{ t } is null, \(\xi_t^2=0\). The null surface in the maximal extension of an eternal black hole consists of two parts: the future horizon and the past horizon. The two intersect on a compact surface of codimension two, Σ, called the bifurcation surface. In the maximally extended spacetime a hypersurface \({\mathcal H}_t\) of constant time t is a Cauchy surface. The bifurcation surface Σ naturally splits the Cauchy surface into two parts, \({{\mathcal H}_}\) and \({{\mathcal H}_+}\), respectively inside and outside the black hole. For asymptoticallyflat spacetime, such as the Schwarzschild metric, the hypersurface \({{\mathcal H}_t}\) has the topology of a wormhole. (In the case of the Schwarzschild metric it is called the EinsteinRosen bridge.) The surface Σ is the surface of minimal area in \({{\mathcal H}_t}\). In fact the bifurcation surface Σ is a minimal surface not only in the (d − 1)dimensional Euclidean space \({{\mathcal H}_t}\), but also in the ddimensional spacetime. As a consequence, as we show below, the components of the extrinsic curvature defined for two vectors normal to Σ, vanish on Σ.
The spacetime in question admits a Euclidean version by analytic continuation t → iτ. It is a feature of regular metrics with a Killing horizon that the direction of Euclidean time τ is compact with period 2πβ_{ H }, which is determined by the condition of regularity, i.e., the absence of a conical singularity. In a vicinity of the bifurcation surface Σ, the spacetime then is a product of a compact surface Σ and a twodimensional disk, the time coordinate τ playing the role of the angular coordinate on the disk. The latter can be made more precise by introducing a new angular variable \(\phi=\beta^{1}_H\tau\), which varies from 0 to 2π. In this paper we consider the spacetime with Euclidean metric of the general type
The radial coordinate ρ is such that the surface Σ is defined by the condition ρ = 0. Near this point the functions g(ρ) and λ_{ ij } (ρ, θ) can be expanded as
where \(\gamma_{ij}^{(0)}(\theta)\) is the metric on the bifurcation surface Σ equipped with coordinates {θ^{i}, i = 1, .., d − 2}. This metric describes what is called a nondegenerate horizon. The Hawking temperature of the horizon is finite in this case and equal to T_{ H } = 1/(2 π β_{ H }).
It is important to note that the metric (50) does not have to satisfy any field equations. The entanglement entropy can be defined for any metric, which possesses a Killingtype horizon. In this sense the entanglement entropy is an offshell quantity. It is useful to keep this in mind when one compares the entanglement entropy with some other approaches in which an entropy is assigned to a black hole horizon. Even though the metric (50) with (51) does not have to satisfy the Einstein equations we shall still call the complete space described by the Euclidean metric (50) the Euclidean black hole instanton and will denote it by E.
Extrinsic curvature of horizon, horizon as a minimal surface
The horizon surface Σ defined by the condition ρ = 0 in the metric (50) is a codimension 2 surface. It has two normal vectors: a spacelike vector n^{1} with the only nonvanishing component \(n^1_{\rho}=1\) and a timelike vector n^{2} with the nonvanishing component \(n^2_\phi=1/\rho\). With respect to each normal vector one defines an extrinsic curvature, \(k_{ij}^a =  \gamma _i^{\;l}\gamma _j^{\;p}{\nabla _l}n_p^a\), a = 1, 2. The extrinsic curvature \(k_{ij}^2\) identically vanishes. It is a consequence of the fact that n^{2} is a Killing vector, which generates time translations. Indeed, the extrinsic curvature can be also written as a Lie derivative, \({k_{\mu \nu}} =  {1 \over 2}{{\mathcal L}_n}{g_{\mu \nu}}\), so that it vanishes if n is a Killing vector. The extrinsic curvature associated to the vector n^{1},
is vanishing when restricted to the surface defined by the condition ρ = 0. It is due to the fact that the term linear in ρ is absent in the ρexpansion for γ_{ ij } (ρ, θ) in the metric (50). This is required by the regularity of the metric (50): in the presence of such a term the Ricci scalar would be singular at the horizon, R ∼ 1/ρ.
The vanishing of the extrinsic curvature of the horizon indicates that the horizon is necessarily a minimal surface. It has the minimal area considered as a surface in ddimensional spacetime. On the other hand, in the Lorentzian signature, the horizon Σ has the minimal area if considered on the hypersurface of constant time t, \({{\mathcal H}_t}\), thus, the latter has the topology of a wormhole.
The wave function of a black hole
Although the entanglement entropy can be defined for any codimension two surface, when the surface is a horizon particular care is required. In order to apply the general prescription outlined in Section 2.1, we first of all need to specify the corresponding wave function. Here we will follow the prescription proposed by Barvinsky, Frolov and Zelnikov [15]. This prescription is a natural generalization of the one in flat spacetime discussed in Section 3.8. On the other hand, it is similar to the “noboundary” wave function of the universe introduced in [127]. We define the wave function of a black hole by the Euclidean path integral over field configurations on the halfperiod Euclidean instanton defined by the metric (50) with angular coordinate ϕ changing in the interval from 0 to π. This halfperiod instanton has Cauchy surface \({\mathcal H}\) (on which we can choose coordinates x = (ρ,θ)) as a boundary where we specify the boundary conditions in the path integral,
where \(W[\psi ] = {1 \over 2}\;\int \psi \hat {{\mathcal D}}\psi\) is the action of the quantum field ψ. The functions ψ_{−}(x) and ψ_{+}(x) are the boundary values defined on the part of the hypersurface \({\mathcal H}\), which is respectively inside \({\mathcal H}_\) and outside \({{\mathcal H}_+}\) the horizon Σ. As was shown in [15], the wave function (53) corresponds to the HartleHawking vacuum state [126].
Reduced density matrix and entropy
f \(\rho(\psi_+^1,\psi_+^2)\) defined by tracing over ψ_{−}modes is given by the Euclidean path integral over field configurations on the complete instanton (0 < ϕ < 2π) with a cut along the axis ϕ = 0 where the field ψ(X) in the path integral takes the values \(\psi _ + ^1(x)\) and \(\psi _ + ^2(x)\) below and above the cut respectively. The trace Trπ is obtained by equating the fields across the cut and doing the unrestricted Euclidean path integral on the complete Euclidean instanton E. Analogously, Trπ^{n} is given by the path integral over field configurations defined on the nfold cover E_{ n } of the complete instanton. This space is described by the metric (50) where angular coordinate ϕ is periodic with period 2πn. It has a conical singularity on the surface Σ so that in a small vicinity of Σ the total space E_{ n } is a direct product of Σ and a twodimensional cone \({{\mathcal C}_n}\) with angle deficit δ = 2π(1 − n). Due to the abelian isometry generated by the Killing vector ∂_{ ϕ } this construction can be analytically continued to arbitrary (noninteger) n → α. So that one can define a partition function
by the path integral over field configurations over E_{ α }, the αfold cover of the instanton E. For a bosonic field described by the field operator \(\hat {\mathcal D}\) one has that \(Z(\alpha)=\det^{1/2}\hat{\mathcal D}\). Defining the effective action as W(α) = − ln Z(α), the entanglement entropy is still given by formula (16), i.e., by differentiating the effective action with respect to the angle deficit. Clearly, only the term linear in (1 − α) contributes to the entropy. Thus, the problem reduces to the calculation of this term in the effective action.
The role of the rotational symmetry
We emphasize that the presence of the rotational symmetry with respect to the Killing vector ∂_{ ϕ }, which generates rotations in the 2plane orthogonal to the entangling surface Σ, plays an important role in our construction. Indeed, without such a symmetry it would be impossible to interpret Trπ^{α} for an arbitrary α as a partition function in some gravitational background. In general, two points are important for this interpretation:

i)
that the spacetime possesses, at least locally near the entangling surface, a rotational symmetry so that, after the identification ϕ → ϕ + 2πα, we get a welldefined spacetime E_{ α }, with no more than just a conical singularity; this holds automatically if the surface in question is a Killing horizon;

ii)
and that the field operator is invariant under the “rotations”, ϕ → ϕ + w; this is automatic if the field operator is a covariant operator.
In particular, point ii) allows us to use the Sommerfeld formula (more precisely its generalization to a curved spacetime) in order to define the Green’s function or the heat kernel on the space E_{ α }. As is shown in [184] (see also discussion in Section 2.13) in the case of the nonLorentz invariant field operators in flat Minkowski spacetime, the lack of the symmetry ii) makes the whole “conical space” approach rather obscure. On the other hand, in the absence of rotational symmetry i) there may appear terms in the entropy that are “missing” in the naively applied conical space approach: the extrinsic curvature contributions [204] or even some curvature terms [134].
In what follows we consider the entanglement entropy of the Killing horizons and deal with the covariant operators so that we do not have to worry about i) or ii).
Thermality of the reduced density matrix of a Killing horizon
The quantum state defined by Eq. (53) is the HartleHawking vacuum [126]. The Green’s function in this state is defined by analytic continuation from the Euclidean Green’s function. The periodicity t → t + iβ_{ H } is thus inherent in this state. This periodicity indicates that the correlation functions computed in this state are in fact thermal correlation functions when continued to the Lorentzian section. This fact generalizes to an arbitrary interacting quantum field as shown in [121]. On the other hand, being globally defined, the HartleHawking state is a pure state, which involves correlations between modes localized on different sides of the horizon. However, this state is described by a thermal density matrix if reduced to modes defined on one side of the horizon as was shown by Israel [138]. That the reduced density matrix obtained by tracing over modes inside the horizon is thermal can be formally seen by using angular quantization. Introducing the Euclidean Hamiltonian H_{ E }, which is the generator of rotations with respect to the angular coordinate ϕ defined above, one finds that \(\rho(\psi_+^1,\psi_+^2)=< \psi_+^1\vert e^{2\pi H_E}\vert \psi^2_+>\), i.e., the density matrix is thermal with respect to the Hamiltonian H_{ E } with inverse temperature 2π. This formal proof in Minkowski space was outlined in [152]. The appropriate Euclidean Hamiltonian is then the Rindler Hamiltonian, which generates Lorentz boosts in a direction orthogonal to the surface Σ. In [142] the proof was generalized to the case of generic static spacetimes with bifurcate Killing horizons admitting a regular Euclidean section.
Useful mathematical tools
Curvature of space with a conical singularity
Consider a space Eα, which is an αfold covering of a smooth manifold E along the Killing vector ∂_{ φ }, generating an abelian isometry. Let surface Σ be a stationary point of this isometry so that near Σ the space E_{ α } looks like a direct product, \(\Sigma \times {{\mathcal C}_\alpha}\), of the surface Σ and a twodimensional cone \({{\mathcal C}_\alpha}\) with angle deficit δ = 2π(1 − α). Outside the singular surface Σ the space E_{ α } has the same geometry as a smooth manifold E. In particular, their curvature tensors coincide. However, the conical singularity at the surface Σ produces a singular (deltafunction like) contribution to the curvatures. This was first demonstrated by Sokolov and Starobinsky [195] in the twodimensional case by using topological arguments. These arguments were generalized to higher dimensions in [7]. One way to extract the singular contribution is to use some regularization procedure, replacing the singular space E_{ α } by a sequence of regular manifolds \({\tilde E_\alpha}\). This procedure was developed by Fursaev and Solodukhin in [111]. In the limit \({\tilde E_\alpha} \rightarrow {E_\alpha}\) one obtains the following results [111]:
where δ_{Σ} is the deltafunction, \(\int\nolimits_{\mathcal M} {f{\delta _\Sigma} =} \int\nolimits_\Sigma f; \,{n^k} = n_k^\mu {\partial _\mu}\,,\,k = 1, \, 2\) are two orthonormal vectors orthogonal to the surface \(\Sigma, \, \,({n_\mu}{n_\nu}) = \sum\nolimits_{k = 1}^2 {n_\mu ^k} n_\nu ^k\) and the quantities \(\bar{R}^{\mu\nu}_{\ \ \alpha\beta},\, \bar{R}^{\mu}_{\ \ \nu}\) and \(\bar{R}\) are computed in the regular points E_{ α }/Σ by the standard method.
These formulas can be used to define the integral expressions^{Footnote 3} [111]
where \(\bar{R}_{ii}=\bar{R}_{\mu\nu}n_i^{\mu}n_i^{\nu}\) and \(\bar{R}_{ijij}=\bar{R}_{\mu\nu\lambda\rho}n^{\mu}_in^{\lambda}_i n^{\nu}_j n^{\rho}_j\). We use a shorthand notation for the surface integral \(\int\nolimits_\Sigma\equiv\int\nolimits_\Sigma \sqrt{\gamma}d^{d2}\theta\).
The terms proportional to α in Eqs. (56)–(59) are defined on the regular space E. The terms O((1 − α)^{2}) in Eqs. (57)–(59) are something like a square of the δfunction. They are not welldefined and depend on the way the singular limit \({\tilde E_\beta} \rightarrow {\varepsilon _\beta}\) is taken. However, those terms are not important in the calculation of the entropy since they are of higher order in (1 − α). However, there are certain invariants, polynomial in the Riemann tensor, in which the terms O((1 − α) ^{2}) do not appear at all. Thus, these invariants are well defined on the manifolds with conical singularity. Below we consider two examples of such invariants [111].
Topological Euler number. The topological Euler number of a 2pdimensional smooth manifold ℰ is given by the integral^{Footnote 4}
Suppose that E_{ α } has several singular surfaces (of dimension 2(p − 1)) Σ_{ i }, each with conical deficit 2π(1 − α_{ i }), then the Euler characteristic of this manifold is [111]
A special case is when E_{ α } possesses a continuous abelian isometry. The singular surfaces Σ_{ i } are the fixed points of this isometry so that all surfaces have the same angle deficit α_{ i } = α. The Euler number in this case is [111]
An interesting consequence of this formula is worth mentioning. Since the introduction of a conical singularity can be considered as the limit of certain smooth deformation, under which the topological number does not change, one has χ[E_{ α }] = χ[E_{α= 1}]. Then one obtains an interesting formula reducing the number χ of a manifold E to that of the fixed points set of its abelian isometry [111]
A simple check shows that Eq. (63) gives the correct result for the Euler number of the sphere \(S_\alpha ^d\). Indeed, the fixed points of 2sphere \(S_\alpha ^2\) are its “north” and “south” poles. Each of these points has χ =1 and one gets from Eq. (63): χ[S^{2}] = 1 + 1 = 2. On the other hand, the singular surface of \(S_\alpha ^d\,(d \geq 3)\) is S^{d−2} and from Eq. (63) the known identity χ[S^{d}] = χ[S^{d−2}] follows. Note that Eq. (63) is valid for spaces with continuous abelian isometry and it may be violated for an orbifold with conical singularities.
Lovelock gravitational action. The general Lovelock gravitational action is introduced on a ddimensional Riemannian manifold as the following polynomial [166]
where \(\delta _{[ \ldots ]}^{[ \ldots ]}\) is the totally antisymmetrized product of the Kronecker symbols and k_{ d } is (d − 2)/2 (or (d − 1)/2) for even (odd) dimension d. If the dimension of spacetime is 2p, the action W_{ p } reduces to the Euler number (60) and is thus topological. In other dimensions the action (64) is not topological, although it has some nice properties, which make it interesting. In particular, the field equations, which follow from Eq. (64), are quadratic in derivatives even though the action itself is polynomial in curvature.
On a conical manifold \({{\mathcal M}_\alpha}\), the Lovelock action is the sum of volume and surface parts [111]
where the first term is the action computed at the regular points. As in the case of the topological Euler number, all terms quadratic in (1 − α) mutually cancel in Eq. (65). The surface term in Eq. (65) takes the form of the Lovelock action on the singular surface Σ. It should be stressed that integrals W_{ p } [Σ] are defined completely in terms of the intrinsic Riemann curvature \({R^{ij}}_{kn}\) of Σ
and W_{0} ≡ ∫_{Σ}. Eq. (65) allows us to compute the entropy in the Lovelock gravity by applying the replica formula. In [145] this entropy was derived in the Hamiltonian approach, whereas arguments based on the dimensional continuation of the Euler characteristics have been used for its derivation in [7].
The heat kernel expansion on a space with a conical singularity
The useful tool to compute the effective action on a space with a conical singularity is the heat kernel method already discussed in Section 2.8. In Section 2.9 we have shown how, in flat space, using the Sommerfeld formula (22), to compute the contribution to the heat kernel due to the singular surface Σ. This calculation can be generalized to an arbitrary curved space E_{ α } that possesses, at least locally, an abelian isometry with a fixed point. To be more specific we consider a scalar field operator \({\mathcal D} =  ({\nabla ^2} + X)\), where X is some scalar function. Then, the trace of the heat kernel \(K = {e^{ s{\mathcal D}}}\) has the following small s expansion
where the coefficients in the expansion decompose into bulk (regular) and surface (singular) parts
The regular coefficients are the same as for a smooth space. The first few coefficients are
The coefficients due to the singular surface Σ (the stationary point of the isometry) are
The form of the regular coefficients (69) in the heat kernel expansion has been well studied in physics and mathematics literature (for a review see [219]). The surface coefficient \(a_1^{\Sigma}\) in Eq. (70) was calculated by the mathematicians McKean and Singer [174] (see also [42]). In physics literature this term has appeared in the work of Dowker [69]. (In the context of cosmic strings one has focused more on the Green’s function rather on the heat kernel [3, 100].) The coefficient \(a_2^{\Sigma}\) was first obtained by Fursaev [101] although in some special cases it was known before in works of Donnelly [64, 65].
It should be noted that due to the fact that the surface Σ is a fixed point of the abelian isometry, all components of the extrinsic curvature of the surface Σ vanish. This explains why the extrinsic curvature does not appear in the surface terms (70) in the heat kernel expansion.
General formula for entropy in the replica method, relation to the Wald entropy
As a consequence of the expressions (55) for the curvature of space with a conical singularity that were presented in Section 3.7.1 one obtains a general expression for the entropy. Consider a Euclidean general covariant action
which describes the gravitational field coupled to some matter fields φ_{ A }. In the replica trick we first introduce a conical singularity at the horizon surface Σ with a small angle deficit δ = 2π(1 − α) so that the Riemann curvature obtains a deltalike surface contribution (55) and the gravitational action (71) becomes a function of α. Then applying the replica formula
we get
for the entropy associated to Σ, where tensor Q_{ αβμν } is defined as a variation of action (71) with respect to the Riemann tensor,
If action (71) is local and it does not contain covariant derivatives of the Riemann tensor, then the tensor \(Q^{\mu\nu}_{\ \ \alpha\beta}\) is a partial derivative of the Lagrangian,
Now, as was observed by Myers and Sinha [181] (see also [4]), one can reexpress
where \({\epsilon^{\alpha \beta}} = n_1^\alpha n_2^\beta  n_2^\alpha n_1^\beta\) is the twodimensional volume form in the space transverse to the horizon surface Σ. Then, for a local action (71) polynomial in the Riemann curvature, the entropy (72) takes the form
which is exactly the Wald entropy [221, 144]. It should be noted that Wald’s Noether charge method is an onshell method so that the metric in the expression for the Wald entropy is supposed to satisfy the field equations. On the other hand, the conical singularity method is an offshell method valid for any metric that describes a blackhole horizon. The relation between the onshell and the offshell descriptions will be discussed in Section 4.1.
UV divergences of entanglement entropy for a scalar field
For a bosonic field described by a field operator \({\mathcal D}\) the partition function is \(Z(\alpha) = {\det ^{ 1/2}}{\mathcal D}\). The corresponding effective action W(α) = − ln Z(α) on a space with a conical singularity, E_{ α }, is expressed in terms of the heat kernel \({K_{{E_\alpha}}}(s)\) in a standard way
The entanglement entropy is computed using the replica trick as
Using the small expansion one can, in principle, compute all UV divergent terms in the entropy. However, the surface terms are known only for the first few terms in the expansion (67). This allows us to derive an explicit form for the UV divergent terms in the entropy.
In two dimensions the horizon is just a point and the entanglement entropy diverges logarithmically [33, 152, 71, 85, 196]
In three dimensions the horizon is a circle and the entropy
is linearly divergent.
The leading UV divergence in d dimensions can be computed directly by using the form of the coefficient \(a_1^{\Sigma}\) (70) in the heat kernel expansion [33]
It is identical to expression (28) for the entanglement entropy in flat Minkowski spacetime. This has a simple explanation. To leading order the spacetime near the blackhole horizon is approximated by the flat Rindler metric. Thus, the leading UV divergent term in the entropy is the entanglement entropy of the Rindler horizon. The curvature corrections then show up in the subleading UV divergent terms and in the UV finite terms.
The fourdimensional case is the most interesting since in this dimension there appears a logarithmic subleading term in the entropy. For a scalar field described by a field operator −(∇^{2} + X) the UV divergent terms in the entanglement entropy of a generic 4dimensional black hole read [197]
We note that for a massive scalar field X = −m^{2}.
Of special interest is the case of the 4d conformal scalar field. In this case \(X =  {1 \over 6}R\) and the entropy (82) takes the form
The logarithmic term in Eq. (83) is invariant under the simultaneous conformal transformations of bulk metric g_{ μν } → e^{2σ}g^{μν} and the metric on the surface Σ, γ_{ ij } → e^{2σ}γ_{ ij }. This is a general feature of the logarithmic term in the entanglement entropy of a conformallyinvariant field.
Let us consider some particular examples.
The ReissnerNordström black hole
A black hole of particular interest is the charged black hole described by the ReissnerNordström metric,
This metric has a vanishing Ricci scalar, \(\bar {R} = 0\). It has inner and out horizons, r_{−} and r_{+} respectively, defined by
where is the mass of the black hole and q is the electric charge of the black hole. The two vectors normal to the horizon are characterized by the nonvanishing components \(n_1^\tau = {g^{ 1/2}}(r),\,n_2^r = \sqrt {g(r)}\). The projections of the Ricci and Riemann tensors on the subspace orthogonal to Σ are
Since R = 0 for the ReissnerNordström metric, the entanglement entropy of a massless, minimally coupled, scalar field (X = 0) and of a conformallycoupled scalar field \(X =  {1 \over 6}R\) coincide [197],
where \(A(\Sigma)=4\pi r_+^2\) and \(s({{{r_ }} \over {{r_ +}}})\) represents the UV finite term. Since s is dimensionless it may depend only on the ratio \({{{r_ }} \over {{r_ +}}}\) of the parameters, which characterize the geometry of the black hole.
If the black hole geometry is characterized by just one dimensionful parameter, the UV finite term in Eq. (87) becomes an irrelevant constant. Let us consider two cases when this happens.
The Schwarzschild black hole. In this case r_{−} = 0 (q = 0) and r_{+} = 2m so that the entropy, found by Solodukhin [196], is
Historically, this was the first time when the subleading logarithmic term in entanglement entropy was computed. The leading term in this entropy is the same as in the Rindler space, when the actual blackhole spacetime is approximated by flat Rindler spacetime. This approximation is sometimes argued to be valid in the limit of infinite mass M. However, we see that, even in this limit, there always exists the logarithmic subleading term in the entropy of the black hole that was absent in the case of the Rindler horizon. The reason for this difference is purely topological. The Euler number of the blackhole spacetime is nonzero while it vanishes for the Rindler spacetime; the Euler number of the blackhole horizon (a sphere) is 2, while it is zero for the Rindler horizon (a plane).
The extreme charged black hole. The extreme geometry is obtained in the limit r_{−} → r_{+} (q = m). The entropy of the extreme black hole is found to take the form [197]
Notice that we have omitted the irrelevant constants (0) and (1) in Eq. (88) and (89) respectively.
The dilatonic charged black hole
The metric of a dilatonic black hole, which has mass m, electric charge q and magnetic charge P takes the form [120]:
with the metric functions
where D is the dilaton charge, \(D={P^2q^2 \over 2m}\). The outer and the inner horizons are defined by
The entanglement entropy is defined for the outer horizon at r = r_{+}. The Ricci scalar of metric (90)
vanishes at the outer horizon, r = r_{+}. Therefore, the entanglement entropy associated with the outer horizon is the same for a minimal scalar field (X = 0) and for a conformallycoupled scalar field (\(X =  {1 \over 6}R\)),
where \({A_\Sigma} = 4\pi (r_ + ^2  {D^2})\) is the area of the outer horizon.
It is instructive to consider the black hole with only electric charge (the magnetic charge P = 0 in this case). This geometry is characterized by two parameters: m and q. In this case one finds
so that expression (93) takes the form
In the extremal limit, 2m^{2} = q^{2}, the area of the outer horizon vanishes, A_{Σ} = 0, and the whole blackhole entropy is determined only by the logarithmicallydivergent term^{Footnote 5} (using a different brickwall method a similar conclusion was reached in [114])
In this respect the extreme dilatonic black hole is similar to a twodimensional black hole. Notice that Eq. (95) is positive as it should be since the entanglement entropy is, by definition, a positive quantity.
The calculation of the entanglement entropy of a static black hole is discussed in the following papers [102, 94, 110, 61, 104, 82, 28, 227, 48, 47,46, 135, 137, 176, 196, 197, 117, 118, 114, 115, 116].
Entanglement Entropy of the KerrNewman black hole
The geometry of the rotating black hole is more subtle than that of a static black hole: near the horizon the rotating spacetime is no longer a product of a horizon sphere S_{2} and a twodimensional disk. The other difficulty with applying the technique of the heat kernel to this case is that the Euclidean version of the geometry requires the rotation parameter to be complex. Nevertheless with some care these difficulties can be overcome and the entanglement entropy of a rotating black hole can be computed along the same lines as for a static black hole [170]. In this section we briefly review the results of Mann and Solodukhin [170].
Euclidean geometry of KerrNewman black hole
First we describe the Euclidean geometry in the nearhorizon limit of the KerrNewmann black hole. The Euclidean KerrNewman metric can be written in the form
where the Euclidean time is t = ιτ and the rotation and charge parameters have also been transformed \(a = \imath \hat{a},\, q = \imath \hat {q}\), so that the metric (96) is purely real. Here \(\hat \Delta (r) = (r  {\hat r_ +})(r  {\hat r_ })\), where \({\hat r_ \pm} = m \pm \sqrt {{m^2} + {{\hat a}^2} + {{\hat q}^2}}\), the quantities ω and ῶ take the form
with \({\hat \rho ^2} = {r^2}  {\hat a^2}{\cos ^2}\theta\). This spacetime has a pair of orthogonal Killing vectors
which are the respective analogs of the vectors d_{ τ } and d_{ ϕ } in the (Euclidean) Schwarzschild case. The horizon surface Σ defined by \(r=\hat{r}_+\) is the stationary surface of the Killing vector K. Near this surface the metric (96) is approximately
where \(\hat{\rho}^2_+=\hat{r}^2_+ \hat{a}^2\cos^2\theta\) and
is the metric on the horizon surface Σ up to O(x^{2}), where variable x is defined by the relation \((r\hat{r}_+)={\gamma x^2 \over 4}\) and \(\gamma = 2\sqrt {{m^2} + {{\hat a}^2} + {{\hat q}^2}}\). The angle coordinate \(\psi=\phi+{\hat{a} \over (\hat{r}^2_+\hat{a}^2)} \tau\) and is well defined on Σ. The metric \(ds_{{C_2}}^2\) is that of a twodimensional disk C_{2}
attached to Σ at a point (θ, ψ), where χ = τ − â sin^{2} θ ϕ is an angle coordinate on C_{ 2 }.
Regularity of the metric near the horizon implies the identifications ψ ↔ ψ + 2π and \(\chi \leftrightarrow \chi+4\pi\gamma^{1}\hat{\rho}^2_+\). For this latter condition to hold, independently of θ on the horizon, it is also necessary to identify (τ, ϕ) with (τ+2πβ_{ H }, ϕ − 2πΩβ_{ H }), where \(\Omega={\hat{a} \over(\hat{r}^2_+ \hat{a}^2)}\) is the (complex) angular velocity and \({\beta _H} = (\hat r_ + ^2  {\hat a^2})/\sqrt {{m^2} + {{\hat a}^2} + {{\hat q}^2}}\). The identified points have the same coordinate ψ.
Therefore, near Σ we have the following description of the Euclidean KerrNewman geometry: attached to every point (θ,ψ) of the horizon is a twodimensional disk C_{2} with coordinates (x, χ). The periodic identification of points on C_{2} holds independently for different points on the horizon Σ, even though χ is not a global coordinate. As in the static case, there is an abelian isometry generated by the Killing vector K, whose fixed set is Σ. Locally we have K = ∂_{ χ }. The periodicity is in the direction of the vector K and the resulting Euclidean space E is a regular manifold.
Now consider closing the trajectory of K with an arbitrary period β ≠ β_{ H }. This implies the identification (τ + 2πβ, ϕ − 2πΩβ), and the metric on C_{2} becomes
where \(\chi = \beta \hat \rho _ + ^2{(\hat r_ + ^2  {\hat a^2})^{ 1}}\bar \chi\) is a new angular coordinate, with period 2π. This is the metric of a two dimensional cone with angular deficit \(\delta = 2\pi (1  \alpha),\;\alpha \equiv {\beta \over {{\beta _H}}}\). With this new identification the metric (96) now describes the Euclidean conical space E_{ α } with singular surface Σ.
The difference of the KerrNewman metric from the static case considered above is that the Euclidean space near the bifurcation surface is not a direct product of the surface Σ and twodimensional cone C_{2,α}. Instead, it is a nontrivial foliation of C_{2,α} over Σ. However, this foliation shares certain common features with the static case. Namely, the invariants constructed from quadratic combinations of extrinsic curvature of Σ vanish identically.
Extrinsic curvature of the horizon
In the case of a static black hole we have argued that the presence of an abelian isometry with horizon being the stationary point of the isometry guarantees that the extrinsic curvature identically vanishes on the horizon. In fact this is also true in the case of a rotating black hole. The role of the abelian isometry generated by the Killing vector K is less evident in this case. That is why, in this subsection, following the analysis of [170], we explicitly evaluate the extrinsic curvature for the KerrNewman black hole and demonstrate that quadratic invariants, that can be constructed with the help of the extrinsic curvature, vanish on the horizon.
With respect to the Euclidean metric (96) we may define a pair of orthonormal vectors \(\{{n_a} = n_a^\mu {\partial _\mu},\quad a = 1,2\}\):
Covariantly, these are
The vectors n^{1} and n^{2} are normal to the horizon surface Σ (defined as r = r_{+}, Δ(r = r_{+}) =0), which is a twodimensional surface with induced metric \({\gamma _{\mu \nu}} = {g_{\mu \nu}}  n_\mu ^1n_\nu ^1  n_\mu ^2n_\nu ^2\). With respect to the normal vectors n^{a}, a = 1, 2 one defines the extrinsic curvatures of the surface \(\Sigma:\kappa _{\mu \nu}^a =  \gamma _\mu ^\alpha \gamma _\nu ^\beta {\nabla _\alpha}n_\beta ^a\). The exact expression for the components of extrinsic curvature is given in [170]. The trace of the extrinsic curvature, \({\kappa ^a} = \kappa _{\mu \nu}^a{g^{\mu \nu}}\),
vanishes when restricted to the horizon surface Σ defined by condition \(\hat{\Delta}(r=\hat{r}_+)=0\). Moreover, the quadratic combinations
vanish on the horizon Σ. Consequently, we have \(\kappa^a_{\mu\nu} \kappa^{a\mu\nu}=0\) on the horizon.
Entropy
Applying the conical singularity method to calculate the entanglement entropy of a rotating black hole we have to verify that i) the curvature singularity at the horizon of a stationary black hole behaves in the same way as in the static case and ii) there are no extra surface terms in the heat kernel expansion for the rotating black hole. The first point was explicitly checked in [170]: the curvature formulas (56)–(59) are still valid in the stationary case. Regarding the second point, it was shown by Dowker [70] that for a generic metric with conical singularity at some surface Σ the only modification of the surface terms in the heat kernel expansion (70) are due to the extrinsic curvature of Σ. For example, the surface coefficient \(a_2^\Sigma\) may be modified by integrals over Σ of terms κ^{a}κ_{ a } and \(\kappa _{\mu \nu}^a{\kappa ^{a\mu \nu}}\). Since, as was shown in Section 3.10.2, these terms identically vanish for the KerrNewman metric there is no modification of the surface terms in this case. Thus, the expression for the entropy (82) remains unchanged in the case of a rotating black hole. The Ricci scalar for the KerrNewmann metric is zero, R = 0. The integrals of the projections of Ricci and Riemann tensors over horizon surface are
The analytic continuation of these expressions back to real values of the parameters a and q requires the substitution
With these identities the quantum entropy of the KerrNewman black hole reads [170]
where \(A(\Sigma) = 4\pi (r_ + ^2 + {a^2})\) is the area of the horizon Σ. In the limit a → 0 this expression reduces to that of the ReissnerNordström black hole (87). An interesting and still somewhat puzzling feature of this result is that, in the case of the Kerr black hole, described by the KerrNewman metric with vanishing electric charge (q = 0), the logarithmic term in the entropy does not depend on the rotation parameter a and is the same as in the case of the Schwarzschild black hole. In particular for the extreme Kerr black hole (q = 0, m = a) one has
The entropy of the Kerr black hole in the brickwall model was calculated in [45] and a result different from Eq. (109) was found. However, the subsequent study in [92] has confirmed Eq. (109).
Entanglement entropy as oneloop quantum correction
A natural point of view on the entanglement entropy of a black hole is that this entropy, as was suggested by Callan and Wilczek [33], is the first quantum correction to the BekensteinHawking entropy^{Footnote 6}. Indeed, the BekensteinHawking entropy S_{BH} can be considered as classical, or treelevel, entropy. If we restore the presence of the Planck constant ħ the BekensteinHawking entropy S_{BH} is proportional to 1/ħ while the entanglement entropy S_{ent} is an ħ^{0} quantity. The total entropy of a black hole is then the sum
where all particles that exist in nature contribute to the entanglement entropy S_{ent}.
The statement on the renormalization of the entropy
As explained in Sections 3.9 and 3.10, the entanglement entropy is a UV divergent quantity. The other wellknown quantity, which possesses UV divergences is the effective action. The standard way to handle the UV divergences in the action is to absorb them into a redefinition of the couplings, which appear in the gravitational action. In four dimensions the gravitational action should also include the terms quadratic in the Riemann curvature. The renormalization procedure is well studied and is described in textbooks (see for instance [22]). The idea now is that exactly the same procedure renormalizes the UV divergences in the entropy. In order to demonstrate this statement, consider a minimallycoupled scalar field. For simplicity suppose that the mass of the field vanishes. The bare (treelevel) gravitational action in four dimensions is the sum of the EinsteinHilbert term and all possible combinations quadratic in the Riemann curvature,
where G_{ B }, Λ_{ B }, c_{1,B}, c_{2, B}, c_{3,B} are the bare coupling constants in the gravitational action.
The UV divergences of the gravitational action are computed by the heat kernel method using the small s expansion (69). For a minimal massless field (X = 0 in the scalar field equation) one finds
These divergences are removed by standard renormalization of the gravitational couplings in the bare gravitational action
where G_{ren} and c_{ i,ren } are the renormalized couplings expressed in terms of the bare ones and the UV parameter ϵ;
The treelevel entropy can be obtained by means of the same replica trick, considered in Sections 3.8, 3.9 and 3.10, upon introduction of the conical singularity with a small angle deficit 2π(1 − α), S(G_{ B }, c_{ i,B }) = (α∂_{ α } − 1)W_{gr}(α). The conical singularity at the horizon Σ manifests itself in that a part of the Riemann tensor for such a manifold E_{ α } behaves as a distribution having support on the surface Σ. Using formulas (56)–(59) one finds for the treelevel entropy
Thus, the BekensteinHawking entropy \(S={1 \over 4G}A(\Sigma)\) is modified due to the presence of R^{2}terms in the action (112). It should be noted that Eq. (116) exactly coincides with the entropy computed by the Noether charge method of Wald [221, 144] (the relation between Wald’s method and the method of conical singularity is discussed in [140]).
The UV divergent part of the entanglement entropy of a black hole has already been calculated, see Eq. (82). For a minimal massless scalar, one has
The main point now is that the sum of the UV divergent part (117) of the entanglement entropy and the treelevel entropy (116)
takes again the treelevel form (116) if expressed in terms of the renormalized coupling constants G_{ren}, c_{i,ren} defined in Eq. (115). Thus, the UV divergences in entanglement entropy can be handled by the standard renormalization of the gravitational couplings, so that no separate renormalization procedure for the entropy is required.
It should be noted that the proof of the renormalization statement is based on a nice property of the heat kernel coefficients a_{ n } (68) on space with conical singularity. Namely, up to (1 − α)^{2} terms the exact coefficient \({a_n} = a_n^{{\rm{reg}}} + a_n^\Sigma\) on conical space E_{ α } is equal to the regular volume coefficient \(a_n^{{\rm{reg}}}\) expressed in terms of the complete curvature, regular part plus a deltalike contribution, using relations (55)
The terms quadratic in R^{sing} are not well defined. However, these terms are proportional to (1 − α)^{2} and do not affect the entropy calculation. Thus, neglecting terms of order (1 − α)^{2} in the calculation of entropy, the renormalization of entropy (118) directly follows from the renormalization of the effective action (114).
That the leading 1/ϵ^{2} divergence in the entropy can be handled by the standard renormalization of Newton’s constant G has been suggested by Susskind and Uglum [213] and by Jacobson [141]. That one also has to renormalize the higher curvature couplings in the gravitational action in order to remove all divergences in the entropy of the Schwarzschild black hole was suggested by Solodukhin [196]. For a generic static black hole the renormalization statement was proven by Fursaev and Solodukhin in [112]. In a different approach based on ’t Hooft’s “brickwall model” the renormalization was verified for the ReissnerNordström black hole by Demers, Lafrance and Myers [62]. For the rotating black hole described by the KerrNewman metric the renormalization of the entropy was demonstrated by Mann and Solodukhin [170]. The nonequilibrium aspect (as defining the rate in a semiclassical decay of hot flat space by black hole nucleation) of the black hole entropy and the renormalization was discussed by Barbon and Emparan [12].
Renormalization in theories with a modified propagator
Let us comment briefly on the behavior of the entropy in theories described by a wave operator \({\mathcal D} = F( {\nabla ^2})\), which is a function of the standard Laplace operator ∇^{2}. In flat space this was analyzed in Section 2.12. As is shown in [184] there is a precise relation between the small s expansion of the heat kernel of operator F(−∇^{2}) and that of the Laplace operator −∇^{2}. The latter heat kernel has the standard decomposition
The heat kernel of operator F(−∇^{2}) then has the decomposition [184]
where
In even dimension d the term \({{\mathcal T}_{d/2}}(s) = 1\). This decomposition is valid both for regular manifolds and manifolds with a conical singularity. If a conical singularity is present, the coefficients a_{ n } have the standard decomposition into regular \(a_n^{{\rm{reg}}}\) and surface \(a_n^\Sigma\) parts as in Eq. (68). The surface term for n = 1 is just the area of the surface Σ, while the surface terms with n ≥ 2 contain surface integrals of (n − 1)th power of the Riemann curvature. Thus, Eq. (121) is a decomposition in powers of the curvature of the spacetime.
The functions P_{ n } are defined in Eq. (40). In particular, if F(q) = q^{k} (k > 0) one finds that
The terms with n ≤ d/2 in decomposition (121) produce the UV divergent terms in the effective action and entropy. The term n = d/2 gives rise to the logarithmic UV divergence. In d dimensions the area term in the entropy is the same as in flat spacetime (see Eq.(41)). In four dimensions (d = 4) the UV divergent terms in the entropy are
We note that an additional contribution to the logarithmic term may come from the first term in Eq. (124) (for instance, this is so for the Laplace operator modified by the mass term, F(q) = q + m^{2}).
In the theory with operator F(−∇^{2}) Newton’s constant is renormalized as [184]
while the higher curvature couplings c_{ i }, i = 1, 2, 3 in the effective action are renormalized in the same way as in Eq. (115). The renormalization of G and {c_{ i }} then makes both the effective action and the entropy finite in the exact same way as in the case of the Laplace operator −∇^{2}. Thus, the renormalization statement generalizes to the theories with modified wave operator F(−∇^{2}).
Area law: generalization to higher spin fields
In this section, we will focus only on the leading UV divergent term, proportional to the area of the horizon. The proportionality of the entanglement entropy to the area is known as the “area law”. As we have discussed already for the case of a scalar field, this term in the entanglement entropy of a black hole is the same as in flat spacetime. In flat Minkowski spacetime, for a field of spin s, massive or massless, including the gauge fields, the calculation of entanglement entropy effectively reduces to the scalar field calculation, provided the number of scalar fields is equal to the number of physical degrees of freedom of the spins field in question. The contribution of fermions comes with the weight 1/2. Thus, we can immediately write down the general expression for the entanglement entropy of a quantum field of spin s in d dimensions,
where \({{\mathcal D}_s}(d)\) is (with weight 1/2 for fermionic fields) the number of physical (onshell) degrees of freedom of a particle of spin s in d dimensions. For gauge fields this assumes gauge fixing. In particular one has
for Dirac fermions,
for the gauge vector fields (including the contribution of ghosts), where dim G is the dimension of the gauge group,
for the RaritaSchwinger particles of spin 3/2 with gauge symmetry (gravitinos), and
for massless spin2 particles (gravitons).
The result for Dirac fermions was first obtained by Larsen and Wilczek [162, 163] and later in a paper by Kabat [150]. The contribution of the gauge fields to the entropy was derived by Kabat [150]. The entropy of the RaritaSchwinger spin3/2 particle and of a massless graviton was analyzed by Fursaev and Miele^{Footnote 7} [110].
Renormalization of entropy due to fields of different spin
The effective action of a field of spin s can be written as
The secondorder covariant operators acting on the spins field can be represented in the following general form
where the matrices X_{(s)} depend on the chosen representation of the quantum field and are linear in the Riemann tensor. Here are some examples [43, 44]
where \(\gamma _{AB}^\alpha\) are gammamatrices. The coefficient \(a_1^{{\rm{reg}}}\) in the small s expansion (67)–(69) of the heat kernel of operator (132) has the general form
where D_{ s } (d) is the dimension of the representation of spin s,
D_{ s }(d) can be interpreted as the number of offshell degrees of freedom. Let us consider some particular cases.
Dirac fermions (s = 1/2). The partition function for Dirac fermions is Z_{1/2} = det^{1/2} Δ^{(1/2)}. In this case \({\rm{Tr}}{X_{(1/2)}} = {1 \over 4}\,{2^{[d/2]}}\,R\) R and hence
where \({{\mathcal D}_{1/2}}\) was introduced in Eq. (127). We note that the negative sign in Eq. (136) in combination with the negative sign for fermions in the effective action (131) gives the total positive contribution to Newton’s constant. The renormalization of Newton’s constant due to Dirac fermions is
Comparison of this equation with the UV divergence of entropy (126) for spin1/2 shows that the leading UV divergence in the entropy of spin1/2 field is handled by the renormalization of Newton’s constant in the same manner as it was for a scalar field.
The RaritaSchwinger field (s = 3/2). The partition function, including gauge fixing and the FaddeevPopov ghost contribution, in this case, is
so that the appropriate heat kernel coefficient is
where \({{\mathcal D}_{3/2}}\) is introduced in Eq. (129). The renormalization of Newton’s constant
then, similarly to the case of Dirac fermions, automatically renormalizes the entanglement entropy (126).
However, this property, does not hold for all fields. The main role in the mismatch between the UV divergences in the entanglement entropy and in Newton’s constant is played by the nonminimal coupling terms X^{(s)}, which appear in the field operators (132).
The puzzle of nonminimal coupling
The simplest case to consider is that of a
Nonminimally coupled scalar field. In this case, one has TrX^{(0)} = ξR, where ξ is the parameter of nonminimal coupling. The renormalization of Newton’s constant
is modified due to the presence of the nonminimal coupling ξ in the scalar field operator. At the same time, the entropy calculation on a Ricci flat background (\(\bar {R} = 0\)) is not affected by the nonminimal coupling since the field operator for this background is identical to the minimal one. This simple reasoning shows that the area law in the case of a nonminimally coupled scalar field is the same as in the case of the minimal scalar field (81). Clearly, there is a mismatch between the renormalization of Newton’s constant and the renormalization of the entanglement entropy. One concludes that, in the presence of nonminimal coupling, when the Riemann tensor appears explicitly in the action of the quantum field, the UV divergence of the entanglement entropy cannot be handled by the standard renormalization of Newton’s constant. The mismatch in the entropy is
It is an important fact that there is no known way to give a statistical meaning to this entropy. Moreover, (142) does not have a definite sign and may become negative if ξ is positive. In some respects, this term is similar to the classical BekensteinHawking entropy: both entropies, at least in the framework of the conventional field theory, do not have a welldefined statistical meaning. There is a hope that in string theory the terms similar to Eq. (142) may acquire a better meaning. However, this question is still open.
We should note that on a space with a conical singularity one can consider the Ricci scalar in the nonminimal scalar operator as the complete curvature including the δlike singular term as in Eq. (55). Then, the differential operator −(∇^{2} +ξR) contains a deltalike potential concentrated on the horizon surface Σ. The presence of this potential modifies the surface terms in the heat kernel in such a way that [198]
where \(a_1^\Sigma (\xi = 0)\) is the surface term \(a_1^\Sigma\) (70) without the nonminimal coupling and the term O(1 − α)^{2} is ill defined (something like δ^{2}(0)). However, it does not affect the entropy calculation. If we now apply the replica trick and calculate the entropy corresponding to the theory with the heat kernel with the surface term (143) we get that [198, 163, 16]
This divergence takes a form consistent with the UV divergence of Newton’constant (141). However, we cannot interpret this entropy as a contribution to the entanglement entropy since the presence of the deltalike potential in the Euclidean field operator is not motivated from the point of view of the original Lorentzian theory, for which the entanglement entropy is calculated. Moreover, Eq. (144) is not positive if ξ > 1/6, while the entanglement entropy is supposed to be a positive quantity.
Similar features are shared by other nonminimallycoupled fields.
Abelian vector field. After gauge fixing, the partition function of an abelian gauge field is
where \(\Delta _ + ^{(1)}\) is operator defined in Eq. (132) with sign + in the matrix \(X_{\mu \nu}^{(1)}\) (133). For the effective action W_{eff} = − ln Z we find,
where
and \({{\mathcal D}_1}(d) = d  2\) is the number of onshell degrees of freedom of the abelian vector field. The renormalization of Newton’s constant is
Comparison with Eq. (126) shows that there is again a mismatch between the UV divergences in the entropy and in Newton’s constant. This mismatch stems from the nonminimal term \(X_{\mu \nu}^{(1)}\) in the Laplace type field operator for the vector field.
Massless graviton. The partition function of a massless graviton in d dimensions, after gauge fixing and adding the FaddeevPopov ghost contribution, is
where \(\Delta _  ^{(1)}\) is operator defined in Eq. (132) with sign — in the matrix \(X_{\mu \nu}^{(1)}\) (133). The operator Δ^{(2)} governs the dynamics of the tensor perturbations, which satisfy the condition \({\nabla ^\mu}({h_{\mu \nu}}  {1 \over 2}{g_{\mu \nu}}h) = 0\). The operator Δ^{(0)} is due to the contribution of the conformal mode, while the determinant of operator \(\Delta _  ^{(1)}\) is due to the FaddeevPopov ghosts. Hence, one has in this case that
where \({{\mathcal D}_2}(d) = {{d(d  3)} \over 2}\) is the number of onshell degrees of freedom of a massless spin2 particle. The renormalization of Newton’s constant is
Again, we observe the mismatch between the UV divergent terms in the entropy (126) and in Newton’s constant, this time due to the graviton.
To summarize, the UV divergences in the entanglement entropy of minimallycoupled scalars and fermions are properly renormalized by the redefinition of Newton’s constant. It happens that each minimallycoupled field (no matter bosonic or fermionic) contributes positively to Newton’s constant and positively to the entropy of the black hole. The contributions to both quantities come proportionally, which allows the simultaneous renormalization of both quantities. The mismatch between the UV divergences in the entropy and in Newton’s constant appears for gauge bosons: the abelian (and nonabelian) vector fields and gravitons. The source of the mismatch are those nonminimal terms X^{(s)} in the field operator, which contribute negatively to Newton’s constant and do not make any contribution to the entanglement entropy of the black hole. At the moment of writing of this review the appropriate treatment of the entropy of nonminimallycoupled fields is not yet available.
Comments on the entropy of interacting fields
So far we have considered free fields in a fixed gravitational (black hole) background. The interaction can be included by adding a potential term ∫ d^{d}xV(ψ) to the classical action, here ψ = {ψ^{i}, i = 1, .., N} is the set of fields in question. In the oneloop approximation, one splits ψ = ψ_{ c } + ψ_{ q }, where ψ_{ c } is the classical background field and ψ_{ q } is the quantum field. The integration over ψ_{ q } then reduces to the calculation of the functional determinant of operator \({\mathcal D} + {M^2}({\psi _c})\), where \(M_{ij}^2 = \partial _{ij}^2V({x_c})\). The fields ψ_{ c } representing the classical background are, in general, functions on the curved spacetime. In some cases these fields are constants that minimize the potential V(ψ). The matrix M_{ ij }(x) plays the role of an xdependent mass matrix. In the approximation, when one can neglect the derivatives of the matrix M, the heat kernel of operator \({\mathcal D} + {M^2}\) is presented as the product \({\rm{Tr}}{e^{ s{\mathcal D}}}\cdot{\rm{Tr}}{e^{ s{M^2}}}\). Using the already calculated trace of the heat kernel \({\rm{Tr}}{e^{ s{\mathcal D}}}\) on space with a conical singularity one obtains at oneloop the entanglement entropy of the interacting fields. In d dimensions one obtains [205] (see also [175] for a related discussion)
where we used that \(\int\nolimits_{\epsilon^2}^\infty s^{d/2}e^{M^2s}=\Gamma(1{d\over 2},M^2\epsilon^2)\). Using the asymptotic behavior Γ(−α, x) = α^{−1}x^{−α} +.., we find that the leading UV divergence of the entropy (152) is again (multiplied by N) (81) and is thus not affected by the presence of the interaction in the action. However, the interaction shows up in the subleading UV divergent and the UV finite terms. For instance, in four dimensions on a flat background we find
We see that the leading UV divergent term proportional to the area is not modified by the presence of the selfinteraction. The mass matrix M^{2} is a function of the background field ψ_{ c }. Thus, the result (153) indicates that at treelevel the entropy should contain terms additional to those of the standard area law, which depend on the value of field ψ at the horizon. In order to illustrate this point consider a ψ^{4} model of a single field (a twodimensional model of this type was considered in [151], in four dimensions the role of selfinteraction was discussed in [46])
where we included the term with the nonminimal coupling. In fact, if we had not included this term, it would have been generated by the quantum corrections due to the selfinteraction of the field ψ^{4}. This is a wellknown fact, established in [182]. The renormalized nonminimal coupling in the model (154) is
where we omit the terms of higher order in λ. Splitting the field ψ into classical and quantum parts in Eq. (154) we find \({M^2} = \xi R + \lambda \psi _c^2\). Suppose for simplicity that the background metric is flat. Then to leading order in λ the entanglement entropy (omitting the UV finite terms) is
This entropy should be considered as a quantum correction to the treelevel entropy
which follows from the action W_{gr} + W[ψ]. We see that the logarithmic divergences in Eqs. (156) and (155) agree if ξ = 0. On the other hand, the renormalization of ξ (155) does not make the total entropy S_{tree} + S_{div} completely UV finite. This is yet another manifestation of the puzzling behavior of nonminimal coupling.
Other Related Methods
Euclidean path integral and thermodynamic entropy
In 1977, Gibbons and Hawking [119] developed a method based on the Euclidean path integral for studying the thermodynamics of black holes. In this method one obtains what may be called a thermodynamic entropy. One deals with metrics, which satisfy the gravitational field equations, and thus avoid the appearance of metrics with conical singularities. The entanglement entropy, on the other hand, has a welldefined statistical meaning. In ordinary systems, the thermodynamic entropy and the statistical (microscopic) entropies coincide. For black holes the exact relation between the two entropies can be seen from the following reasoning^{Footnote 8} [198].
Consider a gravitationally coupled system (gravity plus quantum matter fields) at some arbitrary temperature T = (β)^{−1}. A standard way to describe a thermal state of a field system is to use an Euclidean path integral over all fields in question defined on manifold with periodicity 2πβ along the timelike Killing vector. Suppose that it is a priori known that the system includes a black hole. Thus, there exists a surface Σ (horizon), which is a fixed point of the isometry generated by the killing vector. This imposes an extra condition on the possible class of metrics in the path integral. The other condition to be imposed on metrics in the path integral is the asymptotic behavior at infinity: provided the mass M and the electric charge Q of the gravitational configuration are fixed, one has to specify the falloff of the metrics for large values of r. Thus, the Euclidean path integral is
where the integral is taken over βperiodic fields ψ(τ, x^{i}) = ψ(τ + β, x) and over metrics, which satisfy the following conditions:

i)
g_{ μν } possesses an abelian isometry with respect to the Killing vector ∂_{ τ };

ii)
there exists a surface Σ (horizon) where the Killing vector ∂_{ τ } becomes null;

iii)
asymptotic falloff of metric g_{ μν } at large values of radial coordinate r is fixed by the mass M and electric charge Q of the configuration.
Since the inverse temperature β and mass M in the path integral are two independent parameters, the path integral (158) is mostly over metrics, which have a conical singularity at the surface Σ. The integration in Eq. (158) can be done in two steps. First, one computes the integral over matter fields ψ on the background of a metric, which satisfies conditions i), ii) and iii). The result of this integration is the quantity (15) used in the computation of the entanglement entropy,
Semiclassically, the functional integration over metrics in Eq. (158) can be performed in a saddlepoint approximation,
where metric g_{ μν }(β) is a solution to the saddlepoint equation
with the inverse temperature β kept fixed. The solution of this equation is a regular (without conical singularities) metric g_{ μν }(β). This is the onshell metric, which incororates the quantum corrections due to the vacuum polarization by the matter fields. It can also be called an equilibrium configuration, which corresponds to the fixed temperature β^{−1}. In the saddle point approximation there is a constraint relating the charges at infinity M and Q and the inverse temperature β: β = β(M, Q).
The thermodynamic entropy is defined by the total response of the free energy F = −β^{−1} ln Z(β) to a small change in temperature,
and involves, in particular, the derivative of the equilibrium configuration gμν(β) with respect to β
For an equilibrium configuration, satisfying Eq. (161), the second term in Eq. (163) vanishes and thus the total derivative with respect to β coincides with a partial derivative.
Thus, in order to compute the thermodynamic entropy, one may proceed in two steps. First, for a generic metric, which satisfies the conditions i), ii) and iii) compute the offshell entropy using the replica method, i.e., by introducing a small conical singularity at the horizon. This computation is done by taking a partial derivative with respect to β. Second, consider this offshell entropy for an equilibrium configuration, which solves Eq. (161). Since for the classical gravitational action (112) one finds (β∂_{ β } − 1)W_{gr}[β,g] = S(G_{ B }, c_{ iB }) (116) and for the quantum effective action one obtains the entanglement entropy (β∂_{ β } − 1)W[β, g] = S_{ent}, the relation between the entanglement entropy and thermodynamic entropy is given by
Therefore, the entanglement entropy constitutes only a (quantum) part of the thermodynamic entropy of the black hole. The thermodynamic entropy is defined for equilibrium configurations satisfying the quantum corrected Einstein equations (161). Thus, these configurations are not classical solutions to the Einstein equations but incorporate the quantum (oneloop) corrections. These configurations are regular metrics without conical singularities. The UV divergences in the free energy for these configurations are renormalized in a standard way and thus for the thermodynamic entropy the renormalization statement discussed above holds automatically^{Footnote 9}.
In flat spacetime the quantum (oneloop) thermodynamic and statistical entropies coincide as was shown by Allen [2] due to the fact that the corresponding partition functions differ by terms proportional to β. In the presence of black holes the exact relation between the two entropies has been a subject of some debate (see, for example, [86, 199]). However, the analysis made in [104] shows that in the presence of black hole the Euclidean and statistical free energies coincide, provided an appropriate method of regularization is used to regularize both quantities.
’t Hooft’s brickwall model
In 1985 ’t Hooft [214] proposed a model, which was one of the first successful demonstrations that an entropy that scales as an area can be associated, in a rather natural way, to a blackhole horizon. The idea of ’t Hooft’s calculation was to consider a thermal gas of Hawking particles propagating just outside the blackhole horizon. The entropy in the canonical description of the system is calculated by means of the WentzelKramersBrillouin (WKB) approximation. Provided the temperature of the gas is equal to the Hawking temperature, the result of this calculation is unambiguous. However, there is an important subtlety: the density of states of a Hawking particle becomes infinite as one gets closer to the horizon. The reason for this is simple. Close to the horizon all particles effectively propagate in the optical metric. The later is conformally related to the blackhole metric
as follows
where \(d\omega^2_{d2}\) is the metric of the (d − 2)unit sphere. In the optical metric, the nearhorizon region, where the metric function in Eq. (165) can be approximated as
occupies an infinite volume. Clearly, the infinite volume contains an infinite number of states. In order to regularize this infinity ’t Hooft introduced a brick wall, an imaginary boundary at some small distance ϵ from the actual horizon. Then, the regularized optical volume is divergent when ϵ is taken to zero
where \(A(\Sigma) = r_ + ^{d  2}{\Omega _{d  2}}\) is the area of the horizon and \(\epsilon\sim\sqrt{\beta_H(r_\epsilonr_+)}\) is the invariant distance between the brick wall (r = r_{ ϵ }) and the actual horizon (r = r_{+}). The entropy of a gas of massless particles at temperature \(T = \beta _H^{ 1}\) confined in volume V_{opt} in d spacetime dimensions
in the optical metric, is proportional to the horizon area. We should note that the universal behavior of the regularized optical volume (167) in the limit of small ϵ and its proportionality to the horizon area in this limit was important in establishing the result (168).
WKB approximation, PauliVillars fields
In the original calculation by ’t Hooft one considers a minimallycoupled scalar field, which satisfies the KleinGordon equation
on the background of a blackhole metric (165) and imposes a brickwall boundary condition
Consider, for simplicity, the fourdimensional case. Expanding the scalar field in spherical coordinates φ = e^{−iωt}Y_{ l,m }(θ,ϕ)f(r) one finds that Eq. (169) becomes
One uses the WKB approximation in order to find a solution to this equation. In this approximation one represents f(r) = ρ(r)e^{iS(r)}, where ρ(r) is a slowlyvarying function of r while S(r) is a rapidlyvarying phase. One neglects derivatives of ρ(r) and the second derivative of S(r), one obtains the radial function in the form
valid in the region where k^{2}(r) ≥ 0. The latter condition defines a maximal radius r_{ω, l}, which is a solution to the equation k^{2}(r_{ ω,l }) = 0. For a fixed value of the energy E, by increasing the mass m of the particle or the angular momentum l, the radius r_{ ω,l } approaches r_{+} so that the characteristic region where the solution (172) is valid is in fact the near horizon region. One imposes an extra Dirichlet condition φ = 0 at r = r_{ ω,l } so that the oneparticle spectrum becomes discrete
where n is an integer. This relation is used to count the number of oneparticle states that correspond to fixed values of energy ω and angular momentum l,
Calculating the total number of states, which have the same energy E, one has to sum over l. This sum can be approximated by an integral
where r_{ ω } is determined by the condition that k(r, ω) = 0.
In the near horizon region one approximates the metric function g(r) in Eq. (166) by the first two terms in the expansion in powers of (r − r_{+}),
where β_{ H } is the inverse Hawking temperature and r_{+} is the horizon radius. Constant C is related to the curvature of spacetime near the horizon. The radial position of the brick wall is \({r_{\epsilon}} = {r_ +} + {{\pi {{\epsilon}^2}} \over {{\beta _H}}}\), where ϵ is the geodesic distance between the brick wall and the horizon. Focusing only on the brickwall divergent terms, one obtains for the number of states (175)
In a thermal ensemble of scalar particles at fixed temperature T = β^{−1}, each state in the oneparticle spectrum can be occupied by any integer number of quanta. One gets for the free energy
or, integrating by parts,
Substituting Eq. (177) here and using the integrals
one calculates the divergent terms in the free energy
and, using equation \(S = {\beta ^2}{\partial _\beta}F{\vert_{\beta = {\beta _H}}}\), the entropy
Due to the relations
this expression for the entropy can be rewritten in a completely geometric form
The leading term proportional to the area was first calculated in the seminal paper of ’t Hooft [214]. The area law in the brickwall model was also studied in [173, 213, 9, 10, 178].
It is an important observation made by Demers, Lafrance and Myers in [62] that the brickwall divergences are in fact the UV divergences. This can be seen in the PauliVillars regularization as was first done in [62]. Applying the PauliVillars regularization scheme for the fourdimensional scalar field theory studied here, one introduces five regulator fields {φ_{ i }, i = 1, …, 5} of different statistics and masses {m_{ i }, i = 1, …, 5} dependent on the UV cutoff μ [62]. Together with the original scalar φ_{0} = φ (m_{0} = m) these fields satisfy two constraints
where Δ_{ i } = +1 for the commuting fields, and Δ_{ i } = −1 for the anticommuting fields. Not deriving the exact expressions for m_{ i }, we just quote here the following asymptotic behavior
(where b_{1} and b_{2} are some constants), valid in the limit μ → ∞. The total free energy is the sum of all contributions, from the original scalar field and the regulators
It is clear that, due to the constraints (185), all brickwall divergences (with respect to the parameter ϵ) in the free energy (187) and in the entropy cancel. On the other hand, both the free energy and the entropy become divergent if the PauliVillars regulator μ is taken to infinity, thus confirming their identification as UV divergences. For the free energy one finds
and for the entropy one has
Several remarks are in order.

1)
Comparing the entropy calculated in the brickwall model, (184) or (189), with the entanglement entropy (82) we see that the structure of the UV divergent terms in two entropies is similar. The logarithmic terms in Eqs. (82) and (189) (or (184)) are identical if the blackhole metric has vanishing Ricci scalar, R = 0. This is the case, for example, for the ReissnerNordström black hole considered in [62]. However, the logarithmic terms in the two calculations are different if the Ricci scalar is nonzero. This discrepancy appears to arise due to certain limitations of the WKB approximation. In the exact solutions, known explicitly, for example, for a scalar field in a constant curvature spacetime, the mass m always appears in combination \({m^2}  {1 \over 6}R\). However, this is not seen in the WKB approximation (172). In fact, if one makes this substitution everywhere in the above brickwall calculation, the Ricci scalar would appear in the brickwall entropy in a manner, which agrees with the entanglement calculation (82). Moreover, an alternative calculation [104] of the density of states, which does not make use of the WKB approximation, results in an expression for the entropy, which agrees with Eq. (82).

2)
The similarity between the two entropies suggests that the UV divergences in the brickwall entropy (189) can be renormalized by the renormalization of the couplings in the gravitational action in the same way as for the entanglement entropy. That this indeed works was demonstrated in [62].

3)
For a nonminimallycoupled field we have the same problem as in the case of the entanglement entropy. For metrics with R = 0 the nonminimal coupling does not show up in the scalar field equation and does not change the density of states, the free energy and entropy. On the other hand, the nonminimal coupling affects the renormalization of Newton’s constant, even if the background metric is Ricci flat. An attempt was made in [200] to modify the Dirichlet boundary condition at the brick wall and replace it by a more sophisticated condition, which would depend on the value of the nonminimal coupling ξ, so that the resulting entropy would have the UV divergences consistent with the renormalization of Newton’s constant. However, this attempt cannot be considered as successful since it does not reproduce the expected behavior of the entropy for large positive values of ξ.
The calculation of the brickwall entropy for a rotating black hole is more complicated due to the presence of the superradiance modes in the spectrum. This issue is considered in [45, 92, 155, 41, 149, 224, 132, 148, 153].
Euclidean path integral approach in terms of optical metric
Field equation in optical metric. Consider a slightly more general equation than Eq. (169), by including a nonminimal coupling,
in the background of the blackhole metric g_{ μν }, which takes the form (165). The optical metric is conformally related to the blackhole metric, ḡ_{ μν } = e^{2σ}g_{ μν }, where e^{2σ} = 1/g_{ tt }  (in the metric (165) we have that g_{ tt } = g(r)). Eq. (190) can be rewritten entirely in terms of the optical metric ḡ_{ μν } (166) as follows
where \({\varphi _{{\rm{opt}}}} = {e^{{{(d  2)} \over 4}\sigma}}\varphi, \, {\Delta _{{\rm{opt}}}}\), Δ_{opt} is the Laplace operator for spatial part \(\gamma _{ij}^{opt} = {\bar g_{ij}}\) of the optical metric, while the scalar curvature R and the covariant derivative ∇ are defined with respect to the original metric g_{ μν } = e^{−2σ}ḡ_{ μν }. We notice that, since e^{−2σ} = g(r), the effective potential \({\mathcal V}\) in Eq. (191) vanishes at the horizon. This is a general feature of wave equations in the blackhole background: the fields become, effectively, massless in the near horizon region. The frequency ω, which appears in Eq. (171) in the brickwall calculation, is thus an eigenvalue of the operator Ĥ^{2},
The canonical free energy and Euclidean path integral. The canonical free energy (178)
where ω are eigenvalues of the spatial operator Ĥ^{2} and n(ω) is the degeneracy of the energy level ω, can be represented in terms of the Euclidean path integral for a field theory with wave operator (\(\partial _\tau ^2 + {\hat H^2}\)), provided that the Euclidean time τ is a circle with period β. (This property was first clearly formulated by Allen [2].) It order to see this in a rather elementary way, we first notice that
The sum in this expression can be rewritten as a difference of two sums
Each of these sums should be understood in terms of the zetafunction regularization. In particular, using the properties of the Riemann ζfunction, we find
Collecting together Eqs. (194), (195) and (196) one obtains that
The second term in Eq. (197) can be expressed in terms of the Euclidean path integral. This can be seen as follows. In the Euclidean formulation one first makes a Wick rotation of time t → −iτ. The effective action W_{opt} then is defined by means of the Euclidean path integral
At finite temperature one closes the Euclidean time by identifying τ and τ + β. Then, the eigenvalues of operator ∂_{ τ } then are \(i{{2\pi} \over \beta}k\), where k = 0, ±1, ±2, … The effective action can be expressed in terms of the logarithm of the determinant
Comparing this with Eq. (197) and defining the vacuum energy as
one arrives at an expression for the free energy (197)
Evaluation of the effective action in the optical metric. The effective action W_{opt} (198) can be calculated using the heat kernel method. One notes that the heat kernel of the operator \( \partial _\tau ^2 + {\hat H^2}\) takes the form of a product of two heat kernels, for commuting operators \( \partial _\tau ^2\) and Ĥ^{2}. The heat kernel for operator \(\partial _\tau ^2\) is computed explicitly, provided the periodicity condition, τ → τ + β, is imposed. One finds for the trace
So that the effective action takes the form
The operator \({\hat H^2} =  {\Delta _{{\rm{opt}}}} + {\mathcal V}\) (191) is defined for the (d − 1)dimensional metric \(\gamma _{ij}^{opt}\), the spatial part of the optical metric (166). The trace of the heat kernel of operator \({\hat H^2} =  {\Delta _{{\rm{opt}}}} + {\mathcal V}\) (191) can be represented as a series expansion in powers of s,
where the integration is taken over the spatial part of the optical metric (166) and R_{opt} is the Ricci scalar of (d − 1)metric \(\gamma _{ij}^{{\rm{opt}}}\).
For n ≠ 0 the integration over the proper time s in Eq. (203) is regularized for small s due to the thermal exponential factor, so that the UV regulator ϵ can be removed. Interchanging the sum and the integral one obtains
Only the term with n = 0 in Eq. (203) contains the UV divergences. This term in the effective action is proportional to the inverse temperature β and thus it does not make any contribution to the entropy. On the other hand, the n = 0 term gives the free energy at zero temperature. Thus, only the zero temperature contribution to the free energy is UV divergent as was shown by Dowker and Kennedy [73]. In fact the usual way to renormalize the free energy is to subtract the n = 0 term in Eq. (203) or, equivalently, to subtract the zero temperature free energy, F^{R} = F(T) − F(T = 0). With this regularization and using Eq. (205) one obtains a sort of high temperature expansion of the effective action. In d dimensions one finds for the regularized action [74, 75]
where, for d = 4, the term O(β^{4−d}) also contains a logarithmic term ln β. In four dimensions the first two terms in Eq. (206) are the only terms in the effective action, which are divergent when the integration in the optical metric is taken up to the horizon. The first term in Eq. (206) and the respective term in the free energy and entropy is the contribution of a thermal gas in (d − 1) spatial volume V_{opt} at temperature T = β^{−1} in flat spacetime. The other terms in Eq. (206) are curvature corrections to the flat spacetime result as discussed by Dowker and Schofield [74, 75].
Let us focus on the fourdimensional case. Defining the regularized free energy \({F^R} = {\beta ^{ 1}}W_{{\rm{opt}}}^R\) and entropy S_{opt} = β^{2}∂_{ β } F^{R} one finds (provided one imposes the condition β = β_{ H }, after taking the derivative with respect to β)
where we omit terms, which are finite when the volume integration is extended to the horizon. The important observation now is that
where e^{−2σ} = g(r) and R is the scalar curvature of the original blackhole metric. We recall that the latter is conformally related to the optical metric, \(g_{\mu \nu}^{{\rm{opt}}} = {e^{2\sigma}}{g_{\mu \nu}}\), so the relation between the scalar curvature in two spacetimes is
Using this relation and the form of the potential term \({\mathcal V}\) one arrives at Eq. (208).
Introducing the cutoff ϵ as before, \({r_{\epsilon}} = {r_ +} + {{\pi {{\epsilon}^2}} \over {{\beta _H}}}\), one finds for the volume in the optical metric
and
where \({A_ +} = 4\pi r_ + ^2\) is the horizon area and R(r_{+}) is the value of the scalar curvature at the horizon.
Putting everything together, one finds that the entropy in the optical metric is
Comparing this result with the entanglement entropy (82) computed earlier, we find complete agreement. Notice that ϵ here is in fact the IR regulator, the brickwall cutoff, which regularizes the integration in the radial direction in the optical metric. As in’t Hooft’s original calculation, this divergence can be transformed into a UV divergence by introducing the PauliVillars regulator fields of characteristic mass μ The UV divergences of the entropy when μ is taken to infinity are then of the same type as in Eq. (212). This was analyzed in [104].
The Euclidean path integral approach in the optical metric was considered in [15, 58, 57, 9, 10, 177, 136].
We remind the reader that the entanglement entropy (82) is obtained by using the Euclidean path integral in the original blackhole metric with a conical singularity at the horizon (for β ≠ β_{ H }). The blackhole metric and the optical metric are related by a conformal transformation. This transformation is singular at the horizon and in fact produces a topology change: there appears a new boundary at r = r_{+} in the optical metric, which was a tip of the cone in the original blackhole metric. Because of this singular behavior of the conformal transformation, the exact relation between the two Euclidean path integrals is more subtle than for a regular conformal transformation. That the UV divergences in the entropy calculated in these two approaches coincide suggests that the equivalence between the two approaches might extend to the UV finite terms. Although in arbitrary dimension this equivalence may be difficult to prove, the analysis in two dimensions [199] shows that the entanglement entropy and the brickwall entropy are indeed equivalent. This is, of course, consistent with the formal proof outlined in Section 3.6.
Some Particular Cases
In this section we shall consider some particular examples in which the entanglement entropy, including the UV finite terms, can be calculated explicitly.
Entropy of a 2D black hole
In two dimensions the conformal symmetry plays a special role. This has many manifestations. In particular, the conformal symmetry can be used in order to completely reproduce, for a conformal field theory (CFT), the UV finite part of the corresponding gravitational effective action. This is done by integration of the conformal anomaly. For regular twodimensional spacetimes, the result is the wellknown nonlocal Polyakov action. In the presence of a conical singularity the derivation is essentially the same, although one has to take into account the contribution of the singularity. Consider a twodimensional CFT characterized by a central charge c. For a regular twodimensional manifold, the Polyakov action can be written in the form
where the field equation for the field ψ is ∇^{2}ψ = R. On a manifold M^{α} with a conical singularity with angle deficit δ = 2π(1 − α) the Polyakov action is modified by the contribution from the singularity at the horizon Σ (which is just a point in two dimensions) so that [196, 112]
where ψ_{ h } is the value of the field ψ on the horizon. Applying the replica method to the Polyakov action (214) one obtains that the corresponding contribution to the entanglement entropy from the UV finite term in the effective action is
This result agrees with a derivation of Myers [179] who used the Noether charge method of Wald [221] in order to calculate the entropy. The easiest way to compute the function ψ is to use the conformal gauge g_{ μν } = e^{2σ}δ_{ μν } in which ψ = 2σ. Together with the UV divergent part, the complete entanglement entropy in two dimensions is
where Λ is an IR cutoff.
Let the blackhole geometry be described by a 2D metric
where the metric function f(x) has a simple zero at x = x_{+}. Assume that this black hole is placed inside a box of finite size L so that x_{+} ≤ x ≤ L. In order to get a regular space, one closes the Euclidean time τ with period \({\beta _H},\, {\beta _H} = {{4\pi} \over {f{\prime}({x_ +})}}\). It is easy to see that Eq. (217) is conformal to the flat disk of radius \({z_0}\, (\ln z = {{2\pi} \over {{\beta _H}}}\int\nolimits_L^x {{{dx} \over {f(x)}}})\):
where \(\tilde \tau = {{2\pi \tau} \over {{\beta _H}}}\,(0 \leq \tilde \tau \leq 2\pi),\;0 \leq z \leq 1\). So that the entanglement entropy of the 2D black hole takes the form [199, 98]
where we omit the irrelevant term that is a function of (Λ, z_{0}) but not of the parameters of the black hole and have retained dependence on the UV regulator ϵ.
As was shown in [199], the entanglement entropy (219) is identical to the entropy of the thermal atmosphere of quantum excitations outside the horizon in the “brickwall” approach of ’t Hooft [214].
The black hole resides inside a finitesized box and L is the coordinate of the boundary of the box. The coordinate invariant size of the subsystem complimentary to the black hole is \({L_{{\rm{inv}}}} = \int\nolimits_{{x_ +}}^L d x/\sqrt {f(x)}\). Two limiting cases are of interest. In the first, the size of the system L_{inv} is taken to infinity. Then, assuming that the blackhole spacetime is asymptotically flat, we obtain that the entanglement entropy (219) approaches the entropy of the thermal gas,
This calculation illustrates an important feature of the entanglement entropy of a black hole placed in a box of volume V. Namely, the entanglement entropy contains a contribution of the thermal gas that, in the limit of large volume in dimension d, takes the form (7). This is consistent with the thermal nature of the reduced density matrix obtained from the HartleHawking state by tracing over modes inside the horizon.
The other interesting case is when L_{inv} is small. In this case, we find the universal behavior
The universality of this formula lies in the fact that it does not depend on any characteristics of the black hole (mass, temperature) other than the value of the curvature R(x_{+}) at the horizon.
Consider two particular examples.
2D de Sitter spacetime is characterized by the metric function \(f(x) = 1  {{{x^2}} \over {{l^2}}}\) and the Hawking temperature T_{ H } = 1/2πl. In this spacetime the size of the box is bounded from above, L_{inv} ≤ πl. The corresponding entanglement entropy
is a periodic function of L_{inv}.
The string inspired black hole [222, 169] is described by the metric function f(x) = 1 − e^{−λx}. It described an asymptoticallyflat spacetime. The Hawking temperature is \({T_H} = {\lambda \over {4\pi}}\). The entanglement entropy in this case is
The entropy in these two examples resembles the entanglement entropy in flat spacetime at zero temperature (11) and at a finite temperature (12) respectively.
Entropy of 3D BanadosTeitelboimZanelli (BTZ) black hole
BTZ blackhole geometry
The blackhole solution in threedimensional gravity with negative cosmological constant was first obtained in [6] (see also [5] for global analysis of the solution). We start with the blackhole metric written in a form that makes it similar to the fourdimensional Kerr metric. Since we are interested in its thermodynamic behavior, we write the metric in the Euclidean form:
where the metric functions f(r) and N(r) read
and we use the notation
Obviously one has that r_{+}r_{−} = jl. The coordinate ϕ in Eq. (224) is assumed to be periodic with period 2π.
In order to transform the metric (224) to a Lorentzian signature we need to make the analytic transformation τ → it, j → −ij so that
where \(r_ + ^L\) and \(r_  ^L\) are the values in the Lorentzian spacetime. These are the respective radii of the outer and inner horizons of the Lorentzian black hole in (2 + 1) dimensions. Therefore, we must always apply the transformation (227) after carrying out all calculations in the Euclidean geometry in order to obtain the result for the Lorentzian black hole. The Lorentzian version of the metric (224) describes a black hole with mass m and angular momentum J = 2j. The outer horizon is located at r = r_{+}; the respective inverse Hawking temperature is
In the (τ, r) sector of the metric (224) there is no conical singularity at the horizon if the Euclidean time τ is periodic with period β_{ H }. The horizon Σ is a onedimensional space with metric \(ds_\Sigma ^2 = {l^2}d{\psi ^2}\), where \(\psi = {{{r_ +}} \over l}\phi  {{\vert{r_ }\vert} \over {{l^2}}}\tau\) is a natural coordinate on the horizon.
The BTZ space is obtained from the threedimensional maximallysymmetric hyperbolic space H_{3} (sometimes called the global Euclidean antide Sitter space) by making certain identifications. In order to see this one may use the coordinate transformation
In new coordinates (ρ, θ, ψ) the BTZ metric takes the form
which is the metric on the hyperbolic space H_{3}. In this metric the BTZ geometry is defined by identifications

i)
θ → θ + 2π

ii)
\(\theta \rightarrow \theta + 2\pi {{\vert{r_ }\vert} \over l},\quad \psi \rightarrow \psi + 2\pi {{{r_ +}} \over l}\). The outer horizon r = r_{+} in the coordinate system (ρ, θ, ψ) is located at ρ = 0 and ψ is the angular coordinate on the horizon. Notice that the geodesic distance σ between two points with coordinates (ρ, ψ, θ) and (ρ, ψ′, θ′) is
$${\sinh ^2}{\sigma \over {2l}} = {\cosh ^2}\rho \,{\sinh ^2}{{\psi  \psi{\prime}} \over 2} + {\sinh ^2}\rho \,{\sin ^2}{{\theta  \theta{\prime}} \over 2}\,.$$(231)
Heat kernel on regular BTZ geometry
Consider a scalar field with the operator \({\mathcal D} =  ({\nabla ^2} + \xi/{l^2})\). The maximallysymmetric constantcurvature space is a nice example of a curved space in which the heat equation \(({\partial _s} + {\mathcal D})K(x,{x{\prime}},s) = 0\) has a simple, exact, solution. The heat kernel in this case is a function of the geodesic distance σ between two points X and x′. On the global space H_{3} one finds
where μ =1 − ξ. The regular BTZ geometry is defined by identifications i and ii defined above. As is seen from Eq. (231) the geodesic distance and the heat kernel (2320), expressed in coordinates (ρ,ψ,θ), are automatically invariant under identification i. Thus, it remains to maintain identification ii. This is done by summing over images
Using the path integral representation of the heat kernel we would say that the n = 0 term in Eq. (233) is due to the direct way of connecting points x and x′ in the path integral. On the other hand, the n ≠ 0 terms are due to uncontractible winding paths that go n times around the circle.
Heat kernel on conical BTZ geometry
The conical BTZ geometry, which is relevant to the entanglement entropy calculation, is obtained from global hyperbolic space H_{3} by replacing identification i as follows
and not changing identification ii. For α ≠ 1 this Euclidean space has a conical singularity at the horizon (ρ = 0). The heat kernel on the conical BTZ geometry is constructed via the heat kernel (233) on the regular BTZ space by means of the Sommerfeld formula (22)
where K_{BTZ} is the heat kernel (233). The contour Γ is defined in Eq. (22).
For the trace of the heat kernel (234) one finds [171] after computing by residues the contour integral
where \(\bar s = s/{l^2},\,\, {\gamma _n} = {A_ }n/l\), γ_{ n } = A_{−}n/l and Δψ_{ n } = A_{+} = 2πr_{+} and A_{−}.=2πr_{−}). Notice that we have already made the analytical continuation to the values of r_{+} and r_{−} in the Lorentzian geometry.
The entropy
When the trace of the heat kernel on the conical geometry is known one may compute the entanglement entropy by using the replica trick. Then, the entropy is the sum of UV divergent and UV finite parts [171]
where the UV divergent part is
This divergence is renormalized by the standard renormalization of Newton’s constant
in the threedimensional gravitational action.
The UV finite part in the entropy is
where A_{±} = 2πr_{±}/l.
After the renormalization of Newton’s constant, the complete entropy of the BTZ black hole, S_{BTZ} = S_{BH} + S_{ent}, is a rather complicated function of the area of inner and outer horizons. Approximating in Eq. (239) the infinite sum by an integral one finds [171]
where k = A_{−}/A_{+}. The second term on the righthand side of Eq. (240) can be considered to be the oneloop quantum (UVfinite) correction to the classical entropy of a black hole.
For large enough \(\bar{A}_+\equiv {A_+\over l}>>1\) 1 the integral in Eq. (240) goes to zero exponentially and we have the classical BekensteinHawking formula for entropy. On the other hand, for small Ā_{+}, the integral in Eq. (240) behaves logarithmically so that one has [171]
This logarithmic behavior for small values of A_{+} (provided the ratio k = A_{−}/A_{+} is fixed) is universal and independent of the constant ξ (or μ) in the field operator and the area of the inner horizon (A_{−}) of the black hole. Hence, the rotation parameter J enters Eq. (241) only via the area A_{+} of the outer horizon.
The other interesting feature of the entropy (240) is that it always develops a minimum, which is a solution to the equation
This black hole of minimal entropy may be interesting in the context of the final stage of the Hawking evaporation in three dimensions. As follows from the analysis of Mann and Solodukhin [171], the minimum of the entropy occurs for a black hole whose horizon area is of the Planck length, A_{+} ∼ l_{PL} (in threee dimensions l_{PL} ∼ G_{ren}).
Entropy of ddimensional extreme black holes
The extremal black holes play a special role in gravitational theory. These black holes are characterized by vanishing Hawking temperature T_{ H }, which means that in the metric (165) the nearhorizon expansion in the metric function g(r) starts with the quadratic term (r − r_{+})^{2}. Topologically, the true extremal geometry is different from the nonextremal one. Near the horizon the nonextremal static geometry looks like a product of a twodimensional disk (in the plane (r, τ)) and a (d − 2)dimensional sphere. Then, the horizon is the center in the polar coordinate system on the disk. Contrary to this, an extremal geometry in the nearhorizon limit is a product of a twodimensional cylinder and a (d − 2)dimensional sphere. Thus, the horizon in the extremal case is just another boundary rather than a regular inner point, as in the nonextremal geometry. However, one may consider a certain limiting procedure in which one approaches the extremal case staying all the time in the class of nonextremal geometries. This limiting procedure is what we shall call the “extremal limit”. A concrete procedure of this type was suggested by Zaslavsky [226]. One considers a sequence of nonextreme black holes in a cavity at r = r_{ B } and finds that there exists a set of data (r_{+}, r_{ B }, r_{−}) such that the limit \({{{r_ +}} \over {{r_ }}} \rightarrow 1,\quad {{{r_B}} \over {{r_ +}}} \rightarrow 1\) is well defined. Even if one may have started with a rather general nonextremal metric, the limiting geometry is characterized by very few parameters. In this sense, one may talk about “universality” of the extremal limit. In fact, in the most interesting (and tractable) case the limiting geometry is the product of twodimensional hyperbolic space H_{2} with the (d − 2)dimensional sphere. Since the limiting geometry belongs to the nonextreme class, its classical entropy is proportional to the horizon area in accord with the BekensteinHawking formula. Then, the entanglement entropy of the limiting geometry is a oneloop quantum correction to the classical result. The universality we have just mentioned suggests that this correction possesses a universal behavior in the extreme limit and, since the limiting geometry is rather simple, the limiting entropy can be found explicitly. The latter was indeed shown by Mann and Solodukhin in [172].
Universal extremal limit
Consider a static sphericallysymmetric metric in the following form
where \(d\omega _{d  2}^2\) is the metric on the (d − 2)dimensional unit sphere, describing a nonextreme hole with an outer horizon located at r = r_{+}. However, the analysis can be made for a more general metric, in which \(g_{\tau\tau}\neq g^{1}_{rr}\), the limiting geometry is the simplest in the case we consider in Eq. (243). The function g(r) in Eq. (243) can be expanded as follows
It is convenient to consider the geodesic distance l = ∫g^{−1/2}dr as a radial coordinate. Retaining the first two terms in Eq. (244), we find, for r > r_{+}, that
In order to avoid the appearance of a conical singularity at r = r_{+}, the Euclidean time τ in Eq. (243) must be compactified with period 4π/a, which goes to infinity in the extreme limit a → 0. However, rescaling τ → ϕ = τa/2 yields a new variable ϕ having period 2π. Then, taking into account Eq. (245), one finds for the metric (243)
where we have introduced the variable x = lb^{1/2}. To obtain the extremal limit one just takes a → 0. The limiting geometry
is that of the direct product of a 2dimensional space and a (d − 2)sphere and is characterized by a pair of dimensional parameters b^{−1/2} and r_{+}. The parameter r_{+} sets the radius of the (d − 2)dimensional sphere, while the parameter b^{−1/2} is the curvature radius for the (x, ϕ) 2space. Clearly, this twodimensional space is the negative constant curvature space H_{2}. This is the universality we mentioned above: although the nonextreme geometry is in general described by an infinite number of parameters associated with the determining function g(r), the geometry in the extreme limit depends only on two parameters b and r_{+}. Note that the coordinate r is inadequate for describing the extremal limit (247) since the coordinate transformation (245) is singular when a → 0. The limiting metric (247) is characterized by a finite temperature, determined by the 2π periodicity in angular coordinate ϕ.
The limiting geometry (247) is that of a direct product H_{2}× S_{2} of 2d hyperbolic space H_{2} with radius l = b^{−1/2} and a 2D sphere S_{2} with radius l_{1} = r_{+}. It is worth noting that the limiting geometry (247) precisely merges near the horizon with the geometry of the original metric (243) in the sense that all the curvature tensors for both metrics coincide. This is in contrast with, say, the situation in which the Rindler metric is considered to approximate the geometry of a nonextreme black hole: the curvatures of both spaces do not merge in general.
For a special type of extremal black holes l= l_{1}, the limiting geometry is characterized by just one dimensionful parameter. This is the case for the ReissnerNordström black hole in four dimensions. The limiting extreme geometry in this case is the wellknown BertottiRobinson space characterized by just one parameter r_{+}. This space has remarkable properties in the context of supergravity theory that are not the subject of the present review.
Entanglement entropy in the extremal limit
Consider now a scalar field propagating on the background of the limiting geometry (247) and described by the operator
where R_{(d)} is the Ricci scalar. For the metric (247) characterized by two dimensionful parameters l and l_{1}, one has that \({R_{(d)}} =  2/{l^2} + (d  2)(d  3)/l_1^2\). For a ddimensional conformallycoupled scalar field we have \(\xi = {{d  2} \over {4(d  1)}}\). In this case
The calculation of the respective entanglement entropy goes along the same lines as before. First, one allows the coordinate ϕ, which plays the role of the Euclidean time, to have period 2πα. For α ≠ 1 the metric (247) then describes the space \({E^\alpha} = H_2^\alpha \times {S_{d  2}}\) where \(H_2^\alpha\) is the hyperbolic space coinciding with H_{2} everywhere except the point x = 0, where it has a conical singularity with an angular deficit δ = 2π(1 − α). The heat kernel of the Laplace operator ∇^{2} on E^{α} is given by the product
where \({K_{H_2^\alpha}}\) and \(H_2^\alpha\) are the heat kernels, of the Laplace operator on \(H_2^\alpha\) and S_{d−2} respectively. The effective action reads
where ϵ is a UV cutoff. On spaces with constant curvature the heat kernel function is known explicitly [34]. In particular, on a 2D space H_{2} of negative constant curvature, the heat kernel has the following integral representation:
where \(\bar{s}=sl^{2}\). In Eq. (251) σ is the geodesic distance between the points on H_{2}. Between two points (x, ϕ) and (x, ϕ + ∇ϕ) the geodesic distance is given by \({\sinh ^2}{\sigma \over 2} = {\sinh ^2}x{\sin ^2}{{\Delta \phi} \over 2}\). The heat kernel on the conical hyperbolic space \(H_2^\alpha\) can be obtained from (251) by applying the Sommerfeld formula (22). Skipping the technical details, available in [172], let us just quote the result for the trace
where \(\bar s = s/{l^2}\).
Let us denote \(\Theta_{d2}(s)={\rm Tr} K_{S_{d2}}(s)\) the trace of the heat kernel of the Laplace operator −∇^{2} on a (d − 2)dimensional sphere of unit radius. The entanglement entropy in the extremal limit then takes the form [172, 206]
The function k_{ H }(s) has the following smalls expansion
The trace of the heat kernel on a sphere is known in some implicit form. However, for our purposes a representation in a form of an expansion is more useful,
where \({\Omega _{d  2}} = {{2{\pi ^{(d  1)/2}}} \over {\Gamma ((d  1)/2)}}\) is the area of a unit radius sphere S_{d−2}. The first few coefficients in this expansion can be calculated using the results collected in [219],
We shall consider some particular cases. d = 4. The entanglement entropy in the extreme limit is
where \(s({l_1\over l})\) is the UV finite part of the entropy. For minimal coupling (ξ = 0) this result was obtained in [172]. The first term in Eq. (257) is proportional to the horizon area \(A = 4\pi l_1^2\) while the second term is a logarithmic correction to the area law. For conformal coupling ξ= 1/6, the logarithmic term is
d = 5. The entropy is
To simplify the expressions in higher dimensions we consider only the case of the conformal coupling \(\xi={d2\over 4(d1)}\).
The entropy takes the form: d = 6.
d = 7.
d = 8.
Two examples of the extreme geometry are of particular interest.
Entanglement entropy of the round sphere in Minkowski spacetime. Consider a sphere of radius R in flat Minkowski spacetime. One can choose a spherical coordinate system (τ, r, θ^{i}) so that the surface Σ is defined as τ = 0 and r = R, and variables θ^{i}, i = 1,.., d − 2 are the angular coordinates on Σ. The dmetric reads
where \(d\omega^2_{d2}\) is a metric on (d − 2) sphere of unit radius. Metric (264) is conformal to the metric
which describes the product of twodimensional hyperbolic space H_{2} with coordinates (τ, r) and the sphere S_{d−2}. Note that both spaces, H_{2} and S_{d−2}, have the same radius R. Metric (265) describes the spacetime, which appears in the extremal limit of a ddimensional static black hole. In the hyperbolic space H_{2} we can choose a polar coordinate system (ρ, ϕ) with its center at point r = R,
(for small ρ one has that r = R + ρ cos ϕ and τ = ρ sin ϕ as in the polar system in flat spacetime) so that the metric takes the form
In this coordinate system the surface Σ is defined by the condition ρ = 0. In the entanglement entropy of a conformallycoupled scalar field the logarithmic term s_{0} is conformally invariant. Therefore, it is the same [206] for the entropy of a round sphere of radius R in Minkowski spacetime and in the extreme limiting geometry (267). In various dimensions s_{0} can be obtained from the results (258)–(262) by setting l = l_{1} = R. One finds \(s_0={1\over 90}\) in d = 4, \({s_0} =  {1 \over {756}}\) in d = 6 and \({s_0} = {{{\rm{23}}} \over {{\rm{113400}}}}\) in d = 8. For d > 4 the logarithmic term in the entropy of a round sphere has been calculated by Casini and Huerta [38] directly in Minkowski spacetime. They have obtained s_{0} in all even dimensions up to d = 14. Subsequently, Dowker [72] has extended this result to d =16 and d = 18. In arbitrary dimension d the logarithmic term s_{0} can be expressed in terms of Bernulli numbers as is shown in [38] and [72]
Entanglement entropy of the extreme ReissnerNordström black hole in d dimensions. As was shown by Myers and Perry [180] a generalization of the ReissnerNordström solution to higher dimension d > 4 is given by Eq. (243) with
In the extreme limit r_{+} → r_{−}. Expanding Eq. (268) near the horizon one finds, in this limit, that \(b = {(d  3)^2}/r_ + ^2\). Thus, this extreme geometry is characterized by values of radii l = r_{+}/(d − 3) and l_{1} = r_{+}. In dimension d = 4 we have l = l_{1} as in the case considered above. In dimension d > 4 the two radii are different, l ≠ l_{1}. For the conformal coupling the values of logarithmic term s_{0} in various dimensions are presented in Table 1.
Logarithmic Term in the Entropy of Generic Conformal Field Theory
As we have already seen, in even dimensions, there typically appears a logarithmic term in the entanglement entropy. This term is universal in the sense that it does not depend on the scheme, which is used to regularize the UV divergences. In conformal field theories the logarithmic terms in the entropy are closely related to the conformal anomaly. In this section we discuss in detail this aspect and formulate precisely the relation between entanglement entropy and conformal anomalies.
Consider a conformal field theory in d spacetime dimensions. As we have discussed throughout this review the most efficient way to calculate the corresponding entanglement entropy for a nonextremal black hole is to introduce a small angle deficit δ = 2π(1 − α) at the horizon surface Σ, compute the effective action W(E^{α}) on a manifold E^{α} with a singular surface and then apply the replica formula S = (α∂_{ α } − 1)_{α=1}W(E^{α}) and obtain from it the entanglement entropy. In d dimensions the effective action has the general structure
The logarithmic term in this expansion appears only if dimension d is even. Thus, only even d will be considered in this section. The terms a_{ d }, a_{d−2}, … representing the power UV divergences, are not universal, while the term a_{d/2} is universal and is determined by the integrated conformal anomaly. The term w(g^{(α)}) is the UV finite part of the effective action. Under a global rescaling of the metric on E^{α}, g^{(a)} − λ^{2}g^{(α)}, one has
An important property of the expansion (269) for a quantum field theory, which classically is conformally invariant, is that the logarithmic term a_{d/2} is conformally invariant (see [219] and references therein),
On a manifold with a conical singularity at the surface Σ, the coefficients a_{d−2n} have a bulk part and a surface part. To first order in (1 − α) one finds that
For η = d/2 the coefficients \(a_{d/2}^{\rm{bulk}} (E)\) and a_{d/2}(Σ) are the integrated bulk and the surface conformal anomalies respectively. The bulk and surface terms are independently invariant under conformal transformation g → e^{ω}g,
Applying the replica formula one obtains the entanglement entropy
for a regular manifold each term a_{ d }−2_{ n } is the integral of a polynomial of degree η in the Riemann curvature. Respectively, the surface terms \(a_{d  2n}^\Sigma\) are the integral over the singular surface Σ of polynomial of degree η − 1 in the Riemann curvature and its projections onto the subspace orthogonal to surface Σ. The concrete structure of the polynomials depends on the dimension d.
Logarithmic terms in 4dimensional conformal field theory
In four dimensions the bulk conformal anomaly is a combination of two terms, the topological Euler term and the square of the Weyl tensor,
These are, respectively, the conformal anomalies of type A and B. In a theory with n_{ s } particles of spin s, one finds [76] (the contributions of fields of spin 3/2 and 2 can be obtained from Table 2 on p. 180 of the book of Birrell and Davies [22])
The surface contribution to the conformal anomaly can be calculated directly by, for example, the heat kernel method, as in [101]. Although straightforward, the direct computation is technically involved. However, one has a short cut: there is a precise balance, observed in [196] and [112], between the bulk and surface anomalies; this balance is such that, to first order in (1 − α), one can take \({a_2}({E^\alpha}) = a_2^{{\rm{bulk}}}({E^\alpha}) + O{(1  \alpha)^2}\) and use for the Riemann tensor of E^{α} the representation as a sum of regular and singular (proportional to a deltafunction concentrated on surface Σ) parts. The precise expressions are given in [112, 111]. However, this representation is obtained under the assumption that the surface Σ is a stationary point of an abelian isometry and thus has vanishing extrinsic curvature. Under this assumption, one finds that [112, 111] (see also [188])
where \({R_{ijij}} = {R_{\alpha \beta \mu \nu}}n_i^\alpha n_j^\beta n_i^\mu n_j^\nu, \quad {R_{ii}} = {R_{\alpha \beta}}n_i^\alpha n_i^\beta\).
Each surface term in Eq. (276) is invariant under a subclass of conformal transformations, g → e^{− 2ω}g, such that the normal derivatives of ω vanish on surface Σ. Moreover, the surface term due to the bulk Euler number is a topological invariant: using the GaussCodazzi equation
where R_{Σ} is the intrinsic Ricci scalar of the surface and \(k^i_{\alpha\beta}\) is the extrinsic curvature, and in the assumption of vanishing extrinsic curvature the \(a_A^\Sigma\) term, as shown in [111], is proportional to the Euler number of the 2D surface Σ,
where R_{Σ} is the intrinsic curvature of Σ.
For completeness we note that this result can be generalized to an arbitrary codimension 2 surface in 4dimensional spacetime. Then, the conformal transformation is generalized to any function with nonvanishing normal derivative at Σ. The terms with the normal derivatives of ω in the conformal transformation of \(a_2^\Sigma\) can be canceled by adding the quadratic combinations of extrinsic curvature, Trk^{2} and k_{ a }k_{ a }. The analysis presented by Solodukhin [204] (this analysis is based on an earlier consideration by Dowker [70]) results in the following expressions
This is the most general form of the logarithmic term in the entanglement entropy in four spacetime dimensions.
Thus, as follows from Eq. (274), the logarithmic term in the entanglement entropy of a black hole in four dimensions is
For conformal fields of various spin the values of A and B are presented in Eq. (275).
Consider some particular examples. Extreme static geometry. For an extreme geometry, which has the structure of the product H_{2} × S_{2} and characterized by two dimensionful parameters l (radius of H_{2}) and l_{1} (radius of S_{2}), the logarithmic term in the entropy
is determined by both the anomalies of type A and B. In the case of the extreme ReissnerNordström black hole one has l = l_{1} and the logarithmic term (281) is determined only by the anomaly of type A. For a conformal scalar field one has that A = B/3 = 1/90π^{2} and this equation reduces to Eq. (258). As we already discussed, the geometry H_{2} × S_{2} for l = l_{1} is conformal to flat 4dimensional space. Thus, the Weyl tensor vanishes in this case as does its projection to the subspace orthogonal to horizon S_{2}. That is why the type B anomaly does not contribute in this case to the logarithmic term.
The Schwarzschild black hole. In this case, the background is Ricci flat and the logarithmic term is determined by the difference of A and B,
The same is true for any Ricci flat metric. For a conformal scalar field, Eq. (282) reduces to Eq. (88). For a scalar field the relation of the logarithmic term in the entropy and the conformal anomaly was discussed by Fursaev [103]. The logarithmic term vanishes if A = B. In this case the Riemann tensor does not appear in the conformal anomaly (274) so that the anomaly vanishes if the metric is Ricci flat. In particular, the relation A = B can be found from the \({\mathcal N}=4\) superconformal gauge theory, dual to supergravity on AdS_{5}, according to the AdS/CFT correspondence of Maldacena [167]. The conformal anomaly in this theory was calculated in [131]. Nonextreme and extreme Kerr black hole. For a Kerr black hole (characterized by mass m and rotation a) the logarithmic term does not depend on the parameter of the rotation and it takes the same form
as in the case of the Schwarschild metric. In the extreme limit a = m the logarithmic term takes the same value as (283).
Logarithmic terms in 6dimensional conformal field theory
In six dimensions, omitting the total derivative terms, the conformal anomaly is a combination of four different conformal invariants [17]
where I_{1}, I_{2} and I_{3} are cubic in the Weyl tensor
and E_{6} is the Euler density (60)
As was shown in [17] in a free conformal field theory with n_{0} scalars, n_{1/2} Dirac fermions and n_{ B } 2form fields, one has that^{Footnote 10}
Applying the formulas (55) to I_{1}, I_{2} and I_{3} and using the relation (62) for the Euler number, one finds for the logarithmic term in the entanglement entropy of 4dimensional surface Σ in a 6dimensional conformal field theory
where
is the Euler number of the surface Σ, and
where tensors with Latin indices are obtained by contraction with components of normal vectors \(n_\alpha ^i,\,\,i = 1,\, 2\). Note that in Eq. (295) we used for brevity the notation \({\nabla ^2}{C^{ijij}} \equiv n_\alpha ^in_\beta ^jn_\mu ^in_\nu ^j{\nabla ^2}{C^{\alpha \beta \mu \nu}}\). Eqs. (293), (294), and (295) agree with the results obtained in [134].
Let us consider some examples.
6dimensional Schwarzschild black hole. The 6dimensional generalization of the Schwarzschild solution is [180]
where \(d\omega _{{S_4}}^2\) is a metric of unit 4sphere. The area of horizon is \(A_+={8\pi^2\over 3}r_+^4\). The Euler number of the horizon χ[S_{4}] = 2. This metric is Ricci flat so that only the Riemann tensor contributes to the Weyl tensor. The logarithmic term in this case is
It is interesting to note that this term vanishes in the case of the interacting (2, 0) conformal theory, which is dual to supergravity on AdS_{7}. Indeed, in this case one has [131, 17]
so that \(s^{(2,0)}_0=0\). This is as expected. The Riemann tensor does not appear in the conformal anomaly of the strongly interacting (2, 0) theory so that the anomaly identically vanishes if the spacetime is Ricci flat. This property is not valid in the case of the free (2, 0) tensor multiplet [17] so that the logarithmic term of the free multiplet is nonvanishing.
Conformallyflat extreme geometry. In conformallyflat spacetime the Weyl tensor C_{ αβμν } = 0 so that terms (293), (294) and (295) identically vanish. The logarithmic term (291) then is determined by the anomaly of type A only. In particular, this is the case for the extreme geometry H_{2} × S_{4} with equal radii l = l_{1} of two components. One has
for this extreme geometry. This geometry is conformal to flat spacetime and the logarithmic term (299) is the same as for the entanglement entropy of a round sphere in flat 6dimensional spacetime. This generalizes the result discussed in Section 5.3.2 for the entropy of a round sphere due to a conformal scalar field. The result (299), as is shown in [39, 181], generalizes to a spherical entangling surface in a conformallyflat spacetime of any even dimension.
Why might logarithmic terms in the entropy be interesting?
By a logarithmic term we mean both the logarithmically UVdivergent term in the entropy and the UVfinite term, which scales logarithmically with respect to the characteristic size of the black hole. As we have seen, these terms are identical. However, after the UV divergences in the entropy have been renormalized, it is the UV finite term, which scales logarithmically, that will interest us here.

i)
First of all, the logarithmic terms are universal and do not depend on the way the entropy was calculated and on the scheme in which the UV divergences have been regularized. This is in contrast to the power UV divergences in the entropy that depend both on the calculation procedure and on the regularization scheme.

ii)
The logarithmic terms are related to the conformal anomalies. As the conformal anomalies play an important role in the modern theoretical models, any new manifestation of the anomalies merits our special attention. This may be even more important in view of ideas that the conformal symmetry may play a more fundamental role in nature than is usually thought. As is advocated by ’t Hooft in a number of recent papers [217, 218, 216] a crucial ingredient for understanding Hawking radiation and entropy is to realize that gravity itself is a spontaneouslybroken conformal theory.

iii)
For a large class of extremal blackhole solutions, which arise in supergravity theories considered as low energy approximations of string theory, there exists a microscopic calculation of the entropy. This calculation requires a certain amount of unbroken supersymmetry, so that the black holes in question are the Bogomol’nyiPrasadSommerfeld (BPS) type solutions and use the conformal field theory tools, such as the Cardy formula. The Cardy formula predicts certain logarithmic corrections to the entropy (these corrections are discussed, in particular, in [35] and [201]). One may worry whether exactly the same corrections are reproduced in the macroscopic, field theoretical, computation of the entropy. This aspect was studied recently in [8] for black holes in \({\mathcal N}=4\) supergravity and at least some partial (for the entropy due to matter multiplet of the supergravity) agreement with the microscopic calculation has indeed been observed.

iv)
Speaking about the already renormalized entropy of black holes and taking into account the backreaction of the quantum matter on the geometry, the blackhole entropy can be represented as a series expansion in powers of Newton’s constant, \(G/r_g^2\) (the quantity \(1/r_g^2\), where r_{ g } = 2GM is the size of a black hole, is the scale of the curvature at the horizon; thus, the ratio \(G/r^2_g\) measures the strength of gravitational selfinteraction at the horizon) or, equivalently, in powers of \(M^2_{\rm{PL}}/M^2\). In particular, for the Schwarzschild black hole of mass M in four spacetime dimensions, one finds
$$S = 4\pi {{{M^2}} \over {M_{{\rm{PL}}}^2}} + \sigma \ln M + O\left({{{M_{{\rm{PL}}}^2} \over {{M^2}}}} \right)\,.$$(300)The logarithmic term is the only correction to the classical BekensteinHawking entropy that is growing with mass M.

v)
Although the logarithmic term is still negligibly small compared to the classical entropy for macroscopic black holes, it becomes important for small black holes especially at the latest stage of black hole evaporation. In particular, it manifests itself in a modification of the Hawking temperature as a function of mass M [103]. Indeed, neglecting the terms \(O({{M_{{\rm{PL}}}^2} \over 2})\) in the entropy (300) one finds
$$1/{T_H} = {{\partial S} \over {\partial M}} = 8\pi M/M_{{\rm{PL}}}^2 + \sigma {M^{ 1}}$$(301)so that the Hawking temperature T_{ H } ∼ M for small black holes. Depending on the sign of the coefficient σ in Eq. (300) there can be two different scenarios. If σ < 0, then the entropy S(M) as function of mass develops a minimum at some value of M_{min} ∼ M_{PL}. For this value of the mass the temperature (301) becomes infinite. This is the final point (at least in this approximation) of the evaporation for black holes of mass M > M_{min}. It is reached in finite time. Not worrying about exact numerical factors one has
$${{dM} \over {dt}}\sim  T_H^4{A_ +}\sim  T_H^4{M^2}$$(302)for the evaporation rate. For large black holes (M_{0} ≫ M_{min}) the evaporation time is \({t_{{\rm{BH}}}}\sim M_0^3/M_{{\rm{PL}}}^4\). This evaporation time can be obtained by solving Eq. (302) with the classical expression for the Hawking temperature \({T_H}\sim M_{{\rm{PL}}}^2/M\), i.e., without the correction as in Eq. (301). Thus, if there is no logarithmic term in Eq. (300), any black hole evaporates in finite time. If the correction term is present, it becomes important for M_{0} ∼ M_{min}. Then, assuming that M_{min} ∼ M_{PL}, one finds \(t_{\rm{BH}}\sim (M_0M_{rm{PL}})^5/M^4_{\rm{PL}}\) for the evaporation time. On the other hand, a black hole of mass M_{0} < M_{min}, if it exists, evaporates down to zero mass in infinite time. Similar behavior is valid for black holes of arbitrary mass M_{0} if σ > 0 in Eqs. (300) and (301). The evaporation rate considerably slows down for small black holes since the Hawking temperature T_{ h } ∼ M for small M. The black hole then evaporates to zero mass in infinite time. Asymptotically, for large time t, the mass of the black hole decreases as M(t) ∼ t^{−1/5}. Thus, the subPlanckian black holes if (σ < 0) and any black holes if (σ > 0) live much longer than one would have expected if one used the classical expression for the entropy and for the Hawking temperature.
A Holographic Description of the Entanglement Entropy of Black Holes
A popular trend in modern fundamental physics is to reconsider various, sometimes very well known, phenomena from the point of view of holography. Holography is a rather general statement that the physics inside a spatial region can be understood by looking at a certain theory defined on the boundary of the region. This holographic principle was first formulated by ’t Hooft [215] and later generalized by Susskind [212]. For a review of the holographic principle see [24]. A concrete realization of holography is the AdS/CFT correspondence [167, 223, 125]. According to this correspondence the theory of supegravity (more precisely string theory, the low energy regime of which is described by supegravity) in a (d + 1)dimensional antide Sitter spacetime (AdS) is equivalent to a quantum conformal field theory (CFT) defined on a ddimensional boundary of the antide Sitter. There is a precise dictionary definition of how phenomena on one side of the correspondence can be translated into phenomena on the other side. The correspondence has proven to be extremely useful, both for better understanding gravitational physics and quantum field theory. If d = 4, then the CFT on the boundary is known to be an \({\mathcal N}=4\) superconformal gauge theory. This theory is strongly coupled and in many aspects resembles the QCD. Thus utilizing the correspondence one, in particular, may gain some information on how theories of this type behave (for a review on the correspondence and its applications see [1]).
One of the aspects of the AdS/CFT correspondence is geometrical. The boundary theory provides certain boundary conditions for the gravitational theory in bulk so that one may decode the hologram: reconstruct the bulk spacetime from the boundary data. As was analyzed in [60], for this reconstruction, the boundary data one has to specify consist of the boundary metric and the vacuum expectation of the stressenergy tensor of the boundary CFT. The details are presented in [60]; see also [193].
Entanglement entropy is one of the fundamental quantities, which characterize the boundary theory. One would think that it should have an interpretation within the AdS/CFT correspondence. This interpretation was suggested in 2006 by Ryu and Takayanagi [189, 188] (for a review on this proposal see [185]). This proposal is very interesting since it allows one to compute the entanglement entropy in a purely geometrical way (see also [109]).
Holographic proposal for entanglement entropy
Let M be a (d + 1)dimensional asymptotically antide Sitter spacetime. Its conformal boundary is a ddimensional spacetime ∂M. On a slice of constant time t in ∂M one picks a closed (d − 2)dimensional surface Σ and defines the entropy of entanglement with respect to Σ. Now, on the constant t slice of the (d + 1)dimensional antide Sitter spacetime, consider a (d − 1)dimensional minimal surface γ such that its boundary in ∂M is the surface Σ, dγ = Σ. According to [189, 188], the following quantity
where A(γ) is the area of minimal surface γ and G_{d+1} is Newton’s constant in the (d+1)dimensional gravitational theory, is equal to the entanglement entropy one has calculated in the boundary conformal field theory. This holographic proposal has passed many tests and never failed. It has correctly reproduced the entropy in all cases where it is explicitly known. In particular, in two spacetime dimensions (d = 2) it correctly reproduced (10) and (11) for the entropy at finite size and at finite entropy respectively. In higher dimensions (d > 2) the area of minimal surface γ diverges when it is extended till ∂M. This is an important feature, typical of the AdS/CFT correspondence. In fact, instead of a conformal boundary ∂M, one has to consider a regularized boundary ∂M_{ ϵ } located at a small distance ϵ (measured in terms of some radial coordinate ρ). In the AdS/CFT correspondence the divergence in ϵ has the interpretation of a UV divergence in the boundary quantum field theory. Considering the regularized surface γ, which extends to ∂M_{ ϵ }, one finds that its area, to leading order in ϵ, behaves as A(γ) ∼ A(Σ)/ϵ^{d−2}. Taking this behavior of the area, one sees that the proposal (303) correctly reproduces the proportionality of the entropy to the area of surface Σ and its dependence on the UV cutoff ϵ.
However, we note that in certain situations the choice of the minimal surface Σ may not be unique. In particular, if the dividing surface Σ has several components or if the quantum field resides inside a cavity instead of being defined on an infinite space, there is more than one natural choice of the minimal surface γ. Different choices may correspond to different phases in the boundary theory [157].
Proposals for the holographic entanglement entropy of black holes
If one wants to generalize the proposal of Ryu and Takayanagi to black holes, the first step would be to find a (d + 1)dimensional asymptoticallyAdS metric, which solves the Einstein equations in bulk and whose conformal boundary describes a ddimensional black hole. To find such a metric explicitly may be a difficult task, although some exact solutions are known. First of all, it is easy to construct an AdS metric, which gives a de Sitter spacetime on the boundary. The de Sitter horizon is in many aspects, similar to a black hole horizon. Entanglement entropy associated to the de Sitter horizon [130, 139] has the same properties as the entropy of any other Killing horizon. In four dimensions (d = 3) an exact solution to the Einstein equations has been found in [84] that describes a 3D black hole on the boundary. In three dimensions (d = 2), an exact solution, which describes a generic static twodimensional black hole on the boundary, has been found in [194]. On the other hand, the results of [60] show that for any chosen metric on the boundary ∂M, one can find, at least in a small region close to the boundary, an exact solution to the Einstein equations with negative cosmological constant. Exact formulas are given in [60]. Thus, at least principally, it is not a problem to find an asymptotically AdS metric which describes a black hole on the boundary.
The next question is how to choose the minimal surface 7. Emparan’s proposal [83] consists of the following. Suppose the metric on the boundary of asymptoticallyAdS spacetime describes a black hole with a Killing horizon at surface Σ. Presumably, the horizon on the boundary ∂M is extended to the bulk. The bulk horizon is characterized by vanishing extrinsic curvature and is a minimal (d − 1)dimensional surface. Thus, one can choose the bulk horizon to be that minimal surface γ, the area of which should appear in the holographic formula (303). In this construction the Killing horizon Σ is the only boundary of the minimal surface γ. This prescription is perfectly eligible if one computes the entanglement entropy of a black hole in infinite spacetime. In [83] it was applied to the entropy of a black hole residing on the 2brane in the 4d solution of [84]. In [139] a similar prescription is used to compute entanglement entropy of the de Sitter horizon.
On the other hand, in certain situations it is interesting to consider a black hole residing inside a cavity, the “black hole in a box”. Then, as we have learned in the twodimensional case, the entanglement entropy will depend on the size L of the box so that, in the limit of large L, the entropy will have a thermal contribution proportional to the volume of the box. This contribution is additional to the pure entanglement part, which is due to the presence of the horizon Σ. In order to reproduce this dependence on the size of the box, one should use a different proposal. A relevant proposal was suggested by Solodukhin in [202].
Let a ddimensional sphericallysymmetric static black hole with horizon Σ lie on the regularized boundary (with regularization parameter ϵ) of asymptotically antide Sitter spacetime AdS_{d+1} inside a spherical cavity Σ_{ L } of radius L. Consider a minimal dsurface Γ, whose boundary is the union of Σ and Σ_{ L }. Γ has saddle points, where the radial AdS coordinate has an extremum. By spherical symmetry the saddle points form a (d − 2)surface \({\mathcal C}\) with the geometry of a sphere. Consider the subset Γ_{ h } of Γ whose boundary is the union of Σ and \({\mathcal C}\). According to the prescription of [202], the quantity
is equal to the entanglement entropy of the black hole in the boundary of AdS. In particular, it gives the expected dependence of the entropy on the UV regulator ϵ. The minimal surface Γ_{ h }“knows” about the existence of the other boundary Σ_{ L }. That is why Eq. (304) reproduces correctly the dependence of the entropy on the size of the “box”. In [202] this proposal has been verified for d =2 and d = 4.
It should be noted that as far as the UV divergent part of the entanglement entropy is concerned, the two proposals [83, 139] and [202] give the same result. This is due to the fact that the UV divergences come from that part of the minimal surface, which approaches the boundary ∂M of the AdS spacetime. In both proposals this part of the surfaces γ and Γ_{ h } is the same. Thus, the difference is in the finite terms, which are due to global properties of the minimal surface.
From the geometrical point of view, the holographic calculation of the logarithmic term in the entanglement entropy is related to the surface anomalies studied by Graham and Witten [123] (this point is discussed in [191]).
The holographic entanglement entropy of 2D black holes
In order to check the proposal (304) one needs a solution to the bulk Einstein equations that describes a black hole on the boundary of AdS. In three dimensions a solution of this type is known explicitly [194],
where the AdS radius is set to unity. At asymptotic infinity (ρ = 0) of the metric (305) one has the 2D black hole metric
where f(x) has a simple zero in x = x_{+}. The cavity Σ_{ L } is placed at x = L so that x_{+} ≤ χ ≤ L. The regularity of the metric (306) at the horizon x = x_{+} requires that 0 ≤ τ ≤ β_{ H }, = 4π/f′(x_{+}). Note that Eq. (305) is a vacuum solution of the Einstein equations for any function f(x). The regularity of the 3D metric (305) requires that the constant should be related to the Hawking temperature of the twodimensional horizon by b = f′(x_{+}). The geodesic Γ lies in the hypersurface of constant time τ. The induced metric on the hypersurface (ρ, x) has a constant curvature equal to −2 for any function f(x) and is, thus, related by a coordinate transformation to the metric
The exact relation between the two coordinate systems is
where \(z(x)={b\over 2}\int\nolimits_L^x {dx\over f(x)}\). The equation for the geodesic in the metric (307) is \(r={1\over C^2}(ww_0)^2\). The geodesic length between two points lying on the geodesic with radial coordinates r_{1} and r_{2} is
The saddle point of the geodesic is at r_{ m } = 1/C^{2}. The constant C is determined from the condition that the geodesic Γ joins points (x = x_{+}, ρ = ϵ^{2}) and (x = L, ρ = ϵ^{2}) lying on the regularized (with regularization parameter ϵ) boundary. In the limit of small e one finds that
The length of the geodesic Γ_{ h } joining point r_{+} corresponding to (x = x_{+}, ρ = ϵ) and thus the saddle point is \(\gamma ({\Gamma _h}) =  {1 \over 2}\ln {{{C^2}{r_ +}} \over 4}\).
Now, one has to take into account that, in the AdS/CFT correspondence, the value of Newton’s constant in the bulk is related to the number of quantum fields living in the boundary ∂M. In the AdS_{3}/CFT_{2} case, one has that \({1\over G_N}={2c\over 3}\), where c is the central charge of boundary CFT. Thus, one finds that the holographic entropy (304)
where b = f′(x_{+}), indeed coincides with expression (219) for the holographic entanglement entropy of the 2D black hole in conformal field theory. In particular, for large values of L the holographic formula for the entropy correctly reproduces the entropy of thermal gas \({S_{th}} = {{c\pi} \over 3}{T_H}L\) (we remember that T_{ h } = f′(x_{+})/4π). This is a consequence of the choice of the minimal surface Γ_{ h } made in the proposal (304).
Holographic entanglement entropy of higher dimensional black holes
In higher dimensions there is no known exact solution similar to Eq. (305). However, a solution in the form of a series expansion in ρ is available. In the rest of this section we focus on the case of boundary dimension 4. Then one finds [131]
where g_{(0)ij} (x) is the boundary metric and coefficient [131]
is the local covariant function of the boundary metric. Coefficient g_{(4)} is not expressed as a local function of the boundary metric and is related to the stressenergy tensor of the boundary CFT [60], which has an essentially nonlocal nature. h_{(4)} is a local, covariant, function of the boundary metric and is obtained as a variation of the integrated conformal anomaly with respect to the boundary metric [60]. Its explicit form was computed in [60]. h_{(4)} is a traceless tensor and in mathematics literature it is known as the “obstruction” tensor [122]. The explicit form of h_{(4)} or g_{(4)} is not important if one wants to compute the UV divergence terms in the entropy.
One may choose the boundary metric describing a static sphericallysymmetric black hole to take the form
The minimal surface Γ lies in the hypersurface of the constant of 5dimensional spacetime (310). The induced metric on the hypersurface takes the form
where functions \(F(r,\rho)=g^{rr}_{(0)}g_{rr}\) and R^{2}(r, ρ) = gθθ have a ρexpansion due to Eq. (310). The minimal surface Γ can be parameterized by (ρ, θ, φ). Instead of the radial coordinate r, it is convenient to introduce the coordinate \(y=\int dr/\sqrt{f}\) so that, near the horizon, one has r = r_{+} + by^{2}/4 + O(y^{4}). The coordinate y measures the invariant distance along the radial direction. By spherical symmetry, the area to be minimized is
where ρ_{ m } is the saddle point. In the vicinity of the horizon (y ≪ 1), we can neglect the dependence of the functions F(y, ρ) and R^{2}(y, ρ) on the coordinate y. The minimization of the area of the surface gives the equation
The area of the minimal surface Γ_{ h } is then given by the integral
Using Eq. (310) we find that \({\mathcal A}(\rho) = [{{r_ + ^2} \over {{\rho ^2}}} + {{{g_{(2)\theta \theta}}({r_ +})} \over \rho} +..]\). Substituting this expansion into Eq. (316) we find that the first two terms produce divergences (when ϵ goes to zero), which, according to our proposal, are to be interpreted as UV divergences of the entanglement entropy. At the black hole horizon, one has the relation \(2{R_{\theta \theta}}{_{{r_ +}}} = r_ + ^2(R  {R_{aa}})\). Putting everything together and applying proposal (304), one finds for the divergent part
where \(A(\Sigma) = 4\pi r_ + ^2\) is the horizon area.
The logarithmic term in Eq. (317) is related to the logarithmic divergence as (calculated holographically in [131])
in the quantum effective action of boundary CFT. This relation is a particular manifestation of the general formula (280) that relates the logarithmic term in the entropy to the conformal anomalies of type A and B. One notes that in the \({\mathcal N}=4\) superconformal SU(N) gauge theory one has that \(A = B = {{{N^2}} \over {{\pi ^2}}}\).
It should be noted that the UV finite terms and their dependence on the size L_{inv} of the box can be computed in the limit of small L_{inv}. This calculation is given in [202]. In particular, in any even dimension d one finds a universal term in the entropy that takes the form (up to numerical factor) \(S\sim r_ + ^{d  2}{h_{(d)\theta \theta}}({r_ +})L_{{\rm{inv}}}^2\ln {L_{{\rm{inv}}}}\) and is proportional to the value of the “obstruction tensor” on the black hole horizon. The direct calculation of such terms in the entanglement entropy on the CFT side is not yet available.
Can Entanglement Entropy Explain the BekensteinHawking Entropy of Black Holes?
Entanglement entropy of a black hole is naturally proportional to the area of the black hole horizon. This property makes it very similar to the BekensteinHawking entropy assigned to the horizon. This apparent similarity between the two entropies is the main motivation to raise the question of whether the BekensteinHawking entropy is in fact entirely the entropy of entanglement. In this section we discuss problems, which this interpretation has to face, different approaches to solve them and difficulties, which still remain unsolved.
Problems of interpretation of the BekensteinHawking entropy as entanglement entropy
Any approach, which wants to treat the BekensteinHawking entropy as an entanglement entropy, has to answer the following questions:

(i)
The entanglement entropy is a UV divergent quantity, while the BekensteinHawking entropy is a finite quantity, defined with respect to Newton’s constant, which has been measured in experiments. How can these two quantities be equal?

(ii)
The entanglement entropy is proportional to the number of different field species, which exist in nature. On the other hand, the BekensteinHawking entropy does not seem to depend on any number of fields. This problem is known as the “species puzzle”.

(iii)
We have seen that entanglement entropy due to fields, which are nonminimally coupled to gravity, the gauge bosons and gravitons, behave differently from the entropy due to minimallycoupled fields. Since the gauge bosons and gravitons are fields, which are clearly present in nature and thus should contribute to the entropy, does this contribution spoil the possibility of interpreting the black hole entropy as an entanglement entropy?
Entanglement entropy in induced gravity
One, possibly very natural, way, originally proposed by Jacobson [141], to attack these problems is to consider gravity as an induced phenomenon, in the spirit of Sakharov’s ideas [190] (for a review on a modern touch on these ideas see [220]). In this approach the gravitational field is not fundamental but arises as a mean field approximation of the underlying quantum field theory of fundamental particles (the constituents). This is based on the fact that, even if there is no gravitational action at tree level, it will appear at oneloop. The details of this mechanism will, of course, depend on the concrete model.
Model with minimallycoupled fields. To start with, let us consider a simple model in which the constituents are minimallycoupled fields: we consider N_{0} scalars and N_{1/2} Dirac fermions. The induced gravitational action in this model, to lowest order in curvature, is
where the induced Newton’s constant is
N being the number of field species in this model. The renormalization statement, which is valid for the minimallycoupled fields, guarantees that there is a precise balance between the induced Newton’s constant and the entanglement entropy, so that
i.e., the entanglement entropy of the constituents is precisely equal to the BekensteinHawking entropy, expressed in terms of the induced Newton’s constant (319). Thus, if at a fundamental level the constituents in nature were only minimal fields, the BekensteinHawking entropy, as this example shows, would be explained as the entropy of entanglement. Of course, this example ignores the fact that there are elementary particles, namely the gauge bosons, which are nonminimally coupled.
Models with nonminimally coupled fields. In the model with minimal fields the induced Newton’s constant (319) is set by the UV cutoff ϵ. If one wants to deal with the UV finite quantities, one has to add fields, which contribute negatively to Newton’s constant. Excluding nonphysical fields with wrong statistics, the only possibility is to include nonminimallycoupled fields, scalars or vectors. Models of this type have been considered by Frolov, Fursaev and Zelnikov [95, 97, 87]. One considers a multiplet of scalar fields of mass m_{ s } and nonminimal coupling ξ_{ s } and a set of massive Dirac fields with mass m_{ d }. The number of fields and their parameters are fine tuned so that the ultraviolet divergences in the cosmological constant and in Newton’s constant are canceled. The induced Newton’s constant then
is dominated by the mass of the heaviest constituents. However, as soon as we include the nonminimallycoupled fields the precise balance between Newton’s constant and the entanglement entropy is violated, so that the BekensteinHawking entropy S_{BH} = A(Σ)/4G_{ind}, defined with respect to the induced Newton’s constant, is no longer equal to the entanglement entropy. In the model considered in [95, 97, 87] (various models of a similar nature are considered in [88, 96, 93, 91, 92, 89, 106, 105]), the exact relation between two entropies is
where the quantity Q is determined by the expectation value of the nonminimallycoupled scalar fields on the horizon Σ
This quantity is UV divergent. For a single field it is similar to the quantity (142). Thus, the sharp difference between the entanglement entropy and the BekensteinHawking entropy in this model can be summarized as follows: even though the induced Newton’s constant is made UV finite, the entanglement entropy still (and, in fact, always) remains UV divergent. Thus, we conclude that, in the model of Frolov, Fursaev and Zelnikov, the entanglement entropy is clearly different from the BekensteinHawking entropy.
Entropy in braneworld scenario
An interesting example where the BekensteinHawking entropy is apparently induced in the correct way is given in [130]. This example is closely related to the AdS/CFT correspondence discussed in Section 7. In the RandallSundrum setup [187] one may consider the regularized boundary, which appeared in our discussion of Section 7, as a 3brane with Z_{2} symmetry in an antide Sitter spacetime. In the framework of the AdS/CFT correspondence this brane has a description in terms of CFT coupled to gravity at a UV cutoff [124]. If the brane is placed at the distance ρ = ϵ^{2} from the antide Sitter boundary, one obtains that there is a dynamic gravity induced on the brane, with the induced Newton’s constant
where N is the number of colors in the superconformal SU(N) YangMills theory. N^{2} in this case plays the role of the number of species. We notice that, according to the AdS/CFT dictionary, the parameter ϵ, which is an infrared cutoff on the antide Sitter side, is, in fact, a UV cutoff on the CFT side. Now consider a black hole on the 3brane. The BekensteinHawking entropy can then be represented as follows
As Hawking, Maldacena and Strominger [130] suggested, the right hand side of Eq. (325) can be interpreted as the entanglement entropy of N^{2} fields. This interpretation turns out to be the right one, if one uses the holographic entanglement entropy discussed in Section 7. Indeed, taking the leading divergent term in Eq. (317) and noting that, in a Z_{2} brane configuration, this result should be multiplied by a factor of 2, we get exactly the righthand side of Eq. (325). In [130] one considers de Sitter spacetime (so that Eq. (325) is the entropy of the de Sitter horizon in this case) on the brane since it is the simplest brane configuration one can construct in an antide Sitter spacetime. In [83] this proposal was extended to the holographic entanglement entropy of a black hole on a 2brane solution found in [84]. The twodimensional black hole is considered in [107]. Entropy of a generic black hole on the 3brane was considered in [202].
However, there are certain open questions regarding this example. First of all, we should note that the weakly coupled \({\mathcal N} = 4\; \; SU(N)\) supermultiplet contains the YangMills fields (gauge bosons), conformallycoupled scalars and the Weyl fermions [167]. Thus, it is a bit of a mystery how the entanglement entropy of these, mostly nonminimally coupled, fields (gauge bosons and scalars) has managed to become equal to the BekensteinHawking entropy, when recalling the problems with the nonminimal coupling we have discussed in Section 3.16. A part of this mystery is the fact that the holographic regularization (which corresponds to infrared cutoff on the antide Sitter side) does not have a clear analog on the boundary side. Indeed, if we take, for example, a standard heat kernel regularization, we find that the term linear in the scalar curvature R does not appear at all in the effective action produced by the weakly coupled \({\mathcal N} = 4\) superconformal gauge multiplet.
Gravity cutoff
If we compare the two examples when the BekensteinHawking entropy is correctly reproduced, the model of induced gravity with minimallycoupled constituents and the brane world model, we find that the success of the two models is strongly based on a precise relation between Newton’s constant, the number of species and the UV cutoff. This relation can be reformulated in terms of the Planck mass \({M_{{\rm{PL}}}}(1/{G_N}\sim M_{{\rm{PL}}}^2)\) and energy cutoff Λ ∼ 1/ϵ,
where the precise numerical prefactor depends on the concrete model. It is amazing to note that exactly this relation (326) was proposed by Dvali [77, 79, 78] to hold in general in a theory of quantum gravity coupled to a large number of matter species. The arguments, which were used to get this relation, although they include some thought experiments with black holes, are, in fact, unrelated to (and are thus independent of) the entropy. However, it is clear that relation (326), provided it is correct, helps to reproduce precisely the BekensteinHawking entropy as an entropy of entanglement and automatically solves the species puzzle [80].
KaluzaKlein example
One example of when relation (326) holds is the KaluzaKlein model. In this model one starts with a (4 + n)dimensional theory of gravity, which is then compactified so that one has n compact directions, forming, for example, an ntorus with characteristic size R, and 4 noncompact directions, which form our 4dimensional geometry. The higherdimensional Planck scale Λ is considered fundamental in this model and plays the role of the UV cutoff, while the 4dimensional Planck scale M_{PL} (or 4dimensional Newton’s constant) is derived,
Suppose that, in higher dimensions, there are only one particle — the massless graviton. From the fourdimensional point of view one has, additionally to a single massless graviton, a theory of the tower of spin2 massive KaluzaKlein (KK) particles. Truncating the tower at the cutoff Λ, one finds that N = (RΛ)^{n} is precisely the numebr of these KaluzaKlein species. Thus, as was noted in [77, 79], the relation (327) is a particular example of relation (326) in which N should be understood as the number of the KK species.
In the KK example, the entanglement entropy is equal to the BekensteinHawking entropy as demonstrated by Dvali and Solodukhin [80]. Now consider a large black hole with horizon size r_{ g } ≫ R. Such a large black hole fills all compact directions so that, from the higherdimensional point of view, the black hole horizon is a product of a 2sphere of radius r_{ g } and an ndimensional torus of size R. The BekensteinHawking entropy in the (4 + n)dimensional theory is
where \(4\pi r^2_g R^n\) is the area of (4 + n)dimensional horizon. From the 4dimensional point of view this horizon is a 2sphere of radius r_{ g } and the BekensteinHawking entropy in the 4dimensional theory is
We observe that these two entropies are equal so that the two pictures, the higher dimensional and 4dimensional one, are consistent. Let us now discuss the entanglement entropy. In the (4 + n)dimensional theory, there is only one field, the massless graviton. Its entropy is
where the cutoff is the higher dimensional Planck scale Λ. On the other hand, in the 4dimensional theory one computes the entanglement entropy of N KK fields
These two entropies are equal to each other so that the two ways to compute the entanglement entropy agree. Moreover, the entanglement entropy (331) and (330) exactly reproduces the BekensteinHawking entropy (328) and (329).
However, discussing this result, we should note that the massless and massive gravitons are nonminimally coupled particles. It remains to be understood how the problem of the nonminimal coupling is overcome in this example.
Other Directions of Research
In this section we briefly mention some other interesting directions of ongoing research.
Entanglement entropy in string theory
It is generally believed that the problem of the entanglement entropy of a black hole can and should be resolved in string theory. This was originally suggested by Susskind and Uglum [213]. Indeed, taking that string theory is UV finite, the corresponding entropy calculation should result in a finite quantity. More specifically, the effective action of a closed string can be decomposed into powers of string coupling g as g^{2(n−1)}, where n is the genus of the stringworld sheet. The string configurations with spherical topology, n = 0, give a 1/g^{2} contribution. In a low energy approximation this is exactly the contribution to Newton’s constant G ∼ g^{2}. Thus, one may expect that taking into account just n = 0 closed string configurations, one will correctly reproduce both the entanglement entropy and Newton’s constant. In the Euclidean formulation, the prescription of [213] is to look at the zero genus string world sheet, which intersects the Killing horizon. In the Lorentzian picture this corresponds to an open string with both ends attached to the horizon. The higher genus configurations should give some corrections to the n = 0 result. This is a very attractive idea. However, a very little progress has been made in the literature to actually calculate the entanglement entropy directly in string theory. The reason is of course the technical complexity of the problem. Some support to the idea of Susskind and Uglum was found in the work of Kabat, Shenker and Strassler [151], where the entropy in a twodimensional O(N) invariant σmodel was considered. In particular, it was found that the state counting of the entropy in the UV regime may be lost if considered in the low energy (IR) regime. This type of behavior models the situation with the classical BekensteinHawking entropy. Presumably this analysis could be generalized to the string theory σmodel considered either in optical target metric [9, 10] or in the Euclidean metric with a conical singularity at the horizon (as suggested in [33]). Possibly in the latter case the results obtained for strings on orbifolds [63] can be useful (see [50, 49, 51] for earlier works in this direction).
Another promising approach to attack the problem is to use some indirect methods based on dualities. For example, the AdS/CFT correspondence has been used in [25] to relate the entanglement entropy of a string propagating on a gravitational AdS background with a Killing horizon to the thermal entropy of field theory defined on a boundary of AdS and then, eventually, the thermal entropy to the BekensteinHawking entropy of the horizon.
An interesting approach to the entanglement entropy of extremal black holes via AdS_{2}/CFT_{1} duality is considered in [4, 192], where, in particular, one can identify the entanglement entropy and the microcanonical statistical entropy. This approach is based on the earlier work of Maldacena [168] in which the HartleHawking state is identified with an entangled state of two copies of CFT, defined on two boundaries of the maximallyextended BTZ spacetime. In the accurately taken zero temperature limit, the reduced density matrix, obtained by tracing over the states of one copy of CFT, of the extremal black hole is shown to take the form
which describes the maximallyentangled state in the two copies of the CFT_{1} living on the two boundaries of global AdS_{2}. d(N) is the dimension of the Hilbert space of CFT_{1}. The corresponding entanglement entropy S = −Trp ln ρ = ln d(N) then is precisely equal to the microcanonical entropy in the familiar counting of BPS states and thus is equal to the blackhole entropy [32].
Entanglement entropy in loop quantum gravity
Another approach to quantum gravity, sometimes considered as competing with string theory, is loop quantum gravity. In this theory one considers polymeric excitations of the gravitational field represented by the states of spin networks. A spin network is a graph, a network of points with links representing the relation between points. Each link is labeled by a halfinteger j (the label stands for SU(2) representations). To points, or vertices, of a spin network are attached a SU(2) intertwiner, a SU(2) invariant tensor between the representations attached to all the edges linked to the considered vertex. A simpler and more familiar object in particle physics is the Wilson loop. A surface Σ is represented by vertices (punctures), which divide the spin network into two parts. By tracing over states of just one part of the network, one obtains a density matrix. The entanglement entropy then reduces to a sum over intersections of the spin network with the surface Σ [56, 165, 66],
where P is the number of punctures representing Σ. This quantity should be compared to eigenvalues of the operator of area,
Both quantities scale as P for large P, which indicates that the area law is correctly reproduced. However, the exact relation between the two quantities and the classical entropy S_{BH} = A(Σ)/G is not obvious due to ambiguities present in the formalism. The Immirzi parameter γ is one of them. The question, whether the BekensteinHawking entropy is correctly reproduced in this approach, is eventually related to the continuum limit of the theory [143]. As discussed by Jacobson [143], answering this question may require a certain renormalization of Newton’s constant as well as area renormalization. Indeed, quantity (334) represents a microscopic area, which may be related to the macroscopic quantity in a nontrivial way. These issues remain open.
Entropy in noncommutative theories and in models with minimal length
One might have hoped that the UV divergence of the entanglement entropy could be cured in a natural way were the structure of spacetime modified on some fundamental level. For example, if spacetime becomes noncommutative at short distances. This idea was tested in the case of simple fuzzy spaces in [68, 67]. Although the area law has been verified, the entanglement entropy appears to be sensitive to the size of the ignored region, a phenomenon, which may be understood as a UVIR mixing typical for the noncommutative models.
A holographic calculation of the entanglement entropy in noncommutative YangMills theory was considered in [13, 14]. This calculation for a strip of width l shows that for large values of l ≫ l_{ c } compared to some characteristic length l_{ c } ∼ θλ^{1/2}/ϵ, where θ is the parameter of noncommutativity and \(\lambda = g_{Y\, M}^2N\) is the ’t Hooft coupling, then the shortdistance contribution to the entanglement entropy shows an area law of the form
while for smaller values l ∼ l_{ c } the entropy is proportional to the volume. As seen from Eq. (335) the noncommutativity does not improve the UV behavior of the entropy but leads to the renormalization of the effective number of degrees of freedom that may be interpreted as a manifestation of nonlocality of the model.
The other related idea is to consider models in which the Heisenberg uncertainty relation is modified as \(\Delta x\Delta p \geq \,\,{\hbar \over 2}(1 + {\lambda ^2}{(\Delta p)^2})\), which shows that there exists a minimal length Δx ≥ ħλ (for a review on the models of this type see [113]). In a brickwall calculation the presence of this minimal length will regularize the entropy as discussed in [27, 225, 210, 154, 156].
Transplanckian physics and entanglement entropy
One way to check whether the entanglement entropy is sensitive to the way the conventional theory is completed in the UV regime is to study the possible modifications of the standard Lorentz invariant dispersion relation ω^{2} = k^{2} at large values of momentum k (or at short distances). A typical modification is to break the Lorentz invariance as follows ω^{2} = k^{2} + f(k^{2}). This issue was studied in [147] and [40] in the context of the brickwall model. However, the conclusions made in these papers are opposite. According to [147] the entropy is still UV divergent, although the degree of divergence is modified in a way, which depends on the form of function f(k^{2}). On the other hand, [40] claims that the entropy can be made completely UV finite. In a similar claim [186] suggests that the shortdistance finiteness of the 2point correlation function should imply the UV finiteness of the entanglement entropy. The entanglement entropy in a wide class of theories characterized by modified (Lorentz invariant or not) field operators (so that the UV behavior of the modified propagator is improved compared to the standard one) was calculated in [184]. The conclusion reached in [184] (see also discussion in Sections 2.12 and 3.13 of this review) agrees with that of [147]: no matter how good the UV behavior of the propagator is, the entanglement entropy remains UV divergent. That the shortdistance regularity of correlation functions does not necessarily imply that the entanglement entropy is UV finite was pointed out in [183].
Entropy of more general states
The pure quantum state, which is the starting point in the entanglemententropy calculation, should not necessarily be a vacuum state. According to [55, 53, 52, 54], if one starts with a mixed state of the vacuum and an excited state, the entanglement entropy receives some power law corrections,
with β always less than unity, and the power law correction is due to the excited state.
Nonunitary time evolution
An interesting issue discussed in the literature is the time evolution of the entanglement entropy. It was suggested in [26] that the eigenvalues of a reduced density matrix depend on time t. This is not possible if the time evolution of the density matrix is described by a unitary operator. Thus, the time evolution should be nonunitary. In particular, the entanglement entropy should depend on time t. Similar conclusions have been made in [30] and [158, 159, 160], where, in particular, it was shown that the entanglement entropy is an increasing function of time. These observations may have interesting applications for black holes. As was proposed by Hawking [129] the evolution in time of a black hole should be nonunitary, so that a pure initial state may evolve into a mixed state. From the entanglement point of view, this behavior appears not to be in contradiction with the principles of quantum mechanics, rather it is a simple consequence of the entangled nature of the system. The irreversible loss of information due to entanglement is also seen from the evolution of the entropy under the RenormGroup (RG) flow [203, 164].
Concluding remarks
Since the inspiring work of Srednicki in 1993 we have come a long way in understanding the entanglement entropy. Of course, the main motivation for this research always was and still is the attempt to consider the BekensteinHawking entropy as an entropy of entanglement. If successful, this approach would give a universal explication for the entropy of black holes, valid for black holes of any size and mass and carrying any charges. Unfortunately, we are not yet at that point. Many important unresolved problems remain. However, after 17 years of continuous progress it would not be too surprising if we were actually not that far from the final answer. Speaking of future developments, I think that further progress could be made in two main directions: understanding entanglement entropy directly in string theory and resolving the puzzle of nonminimal coupling. Hopefully, at some notsodistant point in time, the efforts made in these directions will provide us with the still missing elements in the picture.
Notes
 1.
For finite matrices this property indicates that the two density matrices have the same eigenvalues.
 2.
If the boundary of Σ is not empty there could be extra terms in the entropy proportional to the “area” of the boundary ∂Σ as was shown in [108]. We do not consider this case here.
 3.
 4.
 5.
 6.
This statement should be taken with some care. Entanglement entropy is a small correction compared to the BekensteinHawking entropy if the UV cutoff 1/ϵ is, for example, on the order of a few GeV (energy scale of the standard model). However, the two entropies are of the same order if the cutoff is at the Planck scale. I thank ’t Hooft for his comments on this point.
 7.
Among other things the authors of [110] observe certain nonsmooth behavior of the heat kernel coefficients for the spin3/2 and spin2 fields in the limit of vanishing angle deficit.
 8.
As I have recently learned (private communication from Myers), Jacobson and Myers (unpublished) had similar ideas back in the 1990s.
 9.
This is true for minimallycoupled matter fields. In the presence of nonminimal couplings there appear extra terms in the thermodynamic entropy, which are absent in the entanglement entropy, as we discussed earlier in Section 3.16.
 10.
Note that the coefficient b_{6} of [17] is related to a_{3} as \({b_6} =  {a_3}\).
References
 [1]
Aharony, O., Gubser, S.S., Maldacena, J.M., Ooguri, H. and Oz, Y., “Large N field theories, string theory and gravity”, Phys. Rep., 323, 183–386, (2000). [DOI], [arXiv:hepth/9905111 [hepth]]. (Cited on page 68.)
 [2]
Allen, B., “Does statistical mechanics equal one loop field theory?”, Phys. Rev. D, 33, 3640, (1986). [DOI]. (Cited on pages 41 and 46.)
 [3]
Allen, B. and Ottewill, A.C., “Effects of curvature couplings for quantum fields on cosmic string spacetimes”, Phys. Rev. D, 42, 2669–2677, (1990). [DOI]. (Cited on page 24.)
 [4]
Azeyanagi, T., Nishioka, T. and Takayanagi, T., “Near extremal black hole entropy as entanglement entropy via AdS_{2}/CFT_{1}”, Phys. Rev. D, 77, 064005, (2008). [DOI], [arXiv:0710.2956 [hepth]]. (Cited on pages 24 and 78.)
 [5]
Bañados, M., Henneaux, M., Teitelboim, C. and Zanelli, J., “Geometry of the (2+1) black hole”, Phys. Rev. D, 48, 1506–1525, (1993). [DOI], [arXiv:grqc/9302012 [grqc]]. (Cited on page 52.)
 [6]
Bañados, M., Teitelboim, C. and Zanelli, J., “The black hole in threedimensional spacetime”, Phys. Rev. Lett., 69, 1849–1851, (1992). [DOI], [arXiv:hepth/9204099 [hepth]]. (Cited on page 52.)
 [7]
Banados, M., Teitelboim, C. and Zanelli, J., “Black hole entropy and the dimensional continuation of the GaussBonnet theorem”, Phys. Rev. Lett., 72, 957–960, (1994). [DOI], [arXiv:grqc/9309026 [grqc]]. (Cited on pages 21 and 23.)
 [8]
Banerjee, S., Gupta, R.K. and Sen, A., “Logarithmic Corrections to Extremal Black Hole Entropy from Quantum Entropy Function”, J. High Energy Phys., 2011(03), 147, (2011). [DOI], [arXiv:1005.3044 [hepth]]. (Cited on page 66.)
 [9]
Barbón, J.L.F., “Horizon divergences of fields and strings in black hole backgrounds”, Phys. Rev. D, 50, 2712–2718, (1994). [DOI], [arXiv:hepth/9402004]. (Cited on pages 44, 49, and 78.)
 [10]
Barbón, J.L.F., “Remarks on thermal strings outside black holes”, Phys. Lett. B, 339, 41–48, (1994). [DOI], [arXiv:hepth/9406209]. (Cited on pages 44, 49, and 78.)
 [11]
Barbón, J.L.F., “Holographic avatars of entanglement entropy”, Nucl. Phys. B (Proc. Suppl.), 192–193, 12–26, (2009). [DOI]. (Cited on page 6.)
 [12]
Barbón, J.L.F. and Emparan, R., “On quantum black hole entropy and Newton constant renormalization”, Phys. Rev. D, 52, 4527–4539, (1995). [DOI], [arXiv:hepth/9502155]. (Cited on page 33.)
 [13]
Barbón, J.L.F. and Fuertes, C.A., “Holographic entanglement entropy probes (non)locality”, J. High Energy Phys., 2008(04), 096, (2008). [DOI], [arXiv:0803.1928 [hepth]]. (Cited on page 79.)
 [14]
Barbón, J.L.F. and Fuertes, C.A., “A Note on the extensivity of the holographic entanglement entropy”, J. High Energy Phys., 2008(05), 053, (2008). [DOI], [arXiv:0801.2153 [hepth]]. (Cited on page 79.)
 [15]
Barvinsky, A.O., Frolov, V.P. and Zelnikov, A., “Wave function of a black hole and the dynamical origin of entropy”, Phys. Rev. D, 51, 1741–1763, (1995). [DOI], [arXiv:grqc/9404036]. (Cited on pages 19 and 49.)
 [16]
Barvinsky, A.D. and Solodukhin, S.N., “Nonminimal coupling, boundary terms and renormalization of the EinsteinHilbert action and black hole entropy”, Nucl. Phys. B, 479, 305–318, (1996). [DOI], [arXiv:grqc/9512047]. (Cited on page 37.)
 [17]
Bastianelli, F., Frolov, S. and Tseytlin, A.A., “Conformal anomaly of (2,0) tensor multiplet in six dimensions and AdS/CFT correspondence”, J. High Energy Phys., 2000(02), 013, (2000). [DOI], [arXiv:hepth/0001041 [hepth]]. (Cited on pages 64 and 65.)
 [18]
Bekenstein, J.D., “Black holes and the second law”, Lett. Nuovo Cimento, 4, 737–740, (1972). [DOI]. (Cited on page 5.)
 [19]
Bekenstein, J.D., “Black holes and entropy”, Phys. Rev. D, 7, 2333–2346, (1973). [DOI]. (Cited on page 5.)
 [20]
Bekenstein, J.D., “Generalized second law of thermodynamics in blackhole physics”, Phys. Rev. D, 9, 3292–3300, (1974). [DOI]. (Cited on page 5.)
 [21]
Bekenstein, J.D., “Do We Understand Black Hole Entropy?”, in Jantzen, R.T. and Mac Keiser, G., eds., The Seventh Marcel Grossmann Meeting on Recent Developments in Theoretical and Experimental General Relativity, Gravitation and Relativistic Field Theories, Proceedings of the meeting held at Stanford University, 24–30 July 1994, pp. 39–58, (World Scientific, Singapore; River Edge, NJ, 1994). [grqc/9409015]. (Cited on page 6.)
 [22]
Birrell, N.D. and Davies, P.C.W., Quantum fields in curved space, Cambridge Monographs on Mathematical Physics, (Cambridge University Press, Cambridge; New York, 1982). [Google Books]. (Cited on pages 31 and 62.)
 [23]
Bombelli, L., Koul, R.K., Lee, J. and Sorkin, R.D., “Quantum source of entropy for black holes”, Phys. Rev. D, 34, 373–383, (1986). [DOI]. (Cited on page 5.)
 [24]
Bousso, R., “The holographic principle”, Rev. Mod. Phys., 74, 825–874, (2002). [DOI], [arXiv:hepth/0203101 [hepth]]. (Cited on page 68.)
 [25]
Brustein, R., Einhorn, M.B. and Yarom, A., “Entanglement interpretation of black hole entropy in string theory”, J. High Energy Phys., 2006(01), 098, (2006). [DOI], [arXiv:hepth/0508217 [hepth]]. (Cited on page 78.)
 [26]
Brustein, R., Einhorn, M.B. and Yarom, A., “Entanglement and Nonunitary Evolution”, J. High Energy Phys., 2007(04), 086, (2007). [DOI], [arXiv:hepth/0609075 [hepth]]. (Cited on page 80.)
 [27]
Brustein, R. and Kupferman, J., “Black hole entropy divergence and the uncertainty principle”, Phys. Rev. D, 83, 124014, (2011). [DOI], [arXiv:1010.4157 [hepth]]. (Cited on page 80.)
 [28]
Bytsenko, A.A., Vanzo, L. and Zerbini, S., “Quantum correction to the entropy of the (2+1)dimensional black hole”, Phys. Rev. D, 57, 4917–4924, (1998). [DOI], [arXiv:grqc/9710106 [grqc]]. (Cited on page 28.)
 [29]
Calabrese, P. and Cardy, J.L., “Entanglement entropy and quantum field theory”, J. Stat. Mech., 2004(06), P06002, (2004). [DOI], [arXiv:hepth/0405152]. (Cited on page 9.)
 [30]
Calabrese, P. and Cardy, J.L., “Evolution of entanglement entropy in onedimensional systems”, J. Stat. Mech., 2005(04), P04010, (2005). [DOI], [arXiv:condmat/0503393 [condmat]]. (Cited on page 80.)
 [31]
Calabrese, P. and Cardy, J., “Entanglement entropy and conformal field theory”, J. Phys. A: Math. Theor., 42, 504005, (2009). [DOI], [arXiv:0905.4013 [condmat.statmech]]. (Cited on page 6.)
 [32]
Callan, C.G. and Maldacena, J.M., “Dbrane Approach to Black Hole Quantum Mechanics”, Nucl. Phys. B, 472, 591–610, (1996). [DOI], [arXiv:hepth/9602043]. (Cited on page 78.)
 [33]
Callan Jr, C.G. and Wilczek, F., “On geometric entropy”, Phys. Lett. B, 333, 55–61, (1994). [DOI], [arXiv:hepth/9401072]. (Cited on pages 5, 9, 25, 31, and 78.)
 [34]
Camporesi, R., “Harmonic analysis and propagators on homogeneous spaces”, Phys. Rep., 196, 1–134, (1990). [DOI]. (Cited on page 57.)
 [35]
Carlip, S., “Logarithmic corrections to black hole entropy from the Cardy formula”, Class. Quantum Grav., 17, 4175–4186, (2000). [DOI], [arXiv:grqc/0005017]. (Cited on page 66.)
 [36]
Carlson, F., Sur une classe de séries de Taylor, Ph.D. Thesis, (Uppsala University, Uppsala, 1914). (Cited on page 11.)
 [37]
Casini, H. and Huerta, M., “Entanglement entropy in free quantum field theory”, J. Phys. A: Math. Theor., 42, 504007, (2009). [DOI], [arXiv:0905.2562 [hepth]]. (Cited on page 6.)
 [38]
Casini, H. and Huerta, M., “Entanglement entropy for the nsphere”, Phys. Lett. B, 694, 167–171, (2010). [DOI], [arXiv:1007.1813 [hepth]]. (Cited on page 59.)
 [39]
Casini, H., Huerta, M. and Myers, R.C., “Towards a derivation of holographic entanglement entropy”, J. High Energy Phys., 2011(05), 036, (2011). [DOI], [arXiv:1102.0440 [hepth]]. (Cited on page 65.)
 [40]
Chang, D., Chu, C.S. and Lin, F.L., “Transplanckian entanglement entropy”, Phys. Lett B., 583, 192–198, (2004). [DOI], [arXiv:hepth/0306055 [hepth]]. (Cited on page 80.)
 [41]
ChangYoung, E., Lee, D. and Yoon, M., “Rotating Black Hole Entropy from Two Different Viewpoints”, Class. Quantum Grav., 26, 155011, (2009). [DOI], [arXiv:0811.3294 [hepth]]. (Cited on page 45.)
 [42]
Cheeger, J., “Spectral geometry of singular Riemannian spaces”, J. Differ. Geom., 18(4), 575–657, (1983). (Cited on page 24.)
 [43]
Christensen, S.M. and Duff, M.J., “New Gravitational Index Theorems and Supertheorems”, Nucl. Phys. B, 154, 301–342, (1979). [DOI]. (Cited on page 35.)
 [44]
Christensen, S.M. and Duff, M.J., “Quantizing Gravity with a Cosmological Constant”, Nucl. Phys. B, 170, 480–506, (1980). [DOI]. (Cited on page 35.)
 [45]
Cognola, G., “Thermodynamics of scalar fields in Kerr’s geometry”, Phys. Rev. D, 57, 6292–6296, (1998). [DOI], [arXiv:grqc/9710118 [grqc]]. (Cited on pages 31 and 45.)
 [46]
Cognola, G., Kirsten, K. and Vanzo, L., “Free and selfinteracting scalar fields in the presence of conical singularities”, Phys. Rev. D, 49, 1029–1038, (1994). [DOI], [arXiv:hepth/9308106 [hepth]]. (Cited on pages 28 and 39.)
 [47]
Cognola, G., Vanzo, L. and Zerbini, S., “One loop quantum corrections to the entropy for a fourdimensional eternal black hole”, Class. Quantum Grav., 12, 1927–1936, (1995). [DOI], [arXiv:hepth/9502006 [hepth]]. (Cited on page 28.)
 [48]
Cognola, G., Vanzo, L. and Zerbini, S., “One loop quantum corrections to the entropy for an extremal ReissnerNordström black hole”, Phys. Rev. D, 52, 4548–4553, (1995). [DOI], [arXiv:hepth/9504064 [hepth]]. (Cited on page 28.)
 [49]
Dabholkar, A., “Quantum corrections to black hole entropy in string theory”, Phys. Lett. B, 347, 222–229, (1995). [DOI], [arXiv:hepth/9409158 [hepth]]. (Cited on page 78.)
 [50]
Dabholkar, A., “Strings on a cone and black hole entropy”, Nucl. Phys. B, 439, 650–664, (1995). [DOI], [arXiv:hepth/9408098 [hepth]]. (Cited on page 78.)
 [51]
Dabholkar, A., “Tachyon condensation and black hole entropy”, Phys. Rev. Lett., 88, 091301, (2002). [DOI], [arXiv:hepth/0111004 [hepth]]. (Cited on page 78.)
 [52]
Das, S. and Shankaranarayanan, S., “How robust is the entanglement entropyarea relation?”, Phys. Rev. D, 73, 121701, (2006). [DOI], [arXiv:grqc/0511066]. (Cited on page 80.)
 [53]
Das, S., Shankaranarayanan, S. and Sur, S., “Powerlaw corrections to entanglement entropy of black holes”, Phys. Rev. D, 77, 064013, (2008). [DOI], [arXiv:0705.2070 [grqc]]. (Cited on page 80.)
 [54]
Das, S., Shankaranarayanan, S. and Sur, S., “Black hole entropy from entanglement: A review”, in Everett, M. and Pedroza, L., eds., Horizons in World Physics, Horizons in World Physics, 268, (Nova Science, Hauppauge, NY, 2009). [arXiv:0806.0402 [grqc]]. (Cited on pages 6 and 80.)
 [55]
Das, S., Shankaranarayanan, S. and Sur, S., “Entanglement and corrections to BekensteinHawking entropy”, arXiv, eprint, (2010). [arXiv:1002.1129 [grqc]]. (Cited on page 80.)
 [56]
Dasgupta, A., “Semiclassical quantisation of spacetimes with apparent horizons”, Class. Quantum Grav., 23, 635–672, (2006). [DOI], [arXiv:grqc/0505017 [grqc]]. (Cited on page 79.)
 [57]
de Alwis, S.P. and Ohta, N., “On the entropy of quantum fields in black hole backgrounds”, arXiv, eprint, (1994). [arXiv:hepth/9412027 [hepth]]. (Cited on page 49.)
 [58]
de Alwis, S.P. and Ohta, N., “Thermodynamics of quantum fields in black hole backgrounds”, Phys. Rev. D, 52, 3529–3542, (1995). [DOI], [arXiv:hepth/9504033 [hepth]]. (Cited on page 49.)
 [59]
de Boer, J., Kulaxizi, M. and Parnachev, A., “Holographic Entanglement Entropy in Lovelock Gravities”, arXiv, eprint, (2011). [arXiv:1101.5781 [hepth]]. (Cited on page 16.)
 [60]
de Haro, S., Solodukhin, S.N. and Skenderis, K., “Holographic reconstruction of spacetime and renormalization in the AdS/CFT correspondence”, Commun. Math. Phys., 217, 595–622, (2001). [DOI], [arXiv:hepth/0002230 [hepth]]. (Cited on pages 68, 69, and 71.)
 [61]
De Nardo, L., Fursaev, D.V. and Miele, G., “Heatkernel coefficients and spectra of the vector Laplacians on spherical domains with conical singularities”, Class. Quantum Grav., 14, 1059–1078, (1997). [DOI], [arXiv:hepth/9610011]. (Cited on page 28.)
 [62]
Demers, J.G., Lafrance, R. and Myers, R.C., “Black hole entropy without brick walls”, Phys. Rev. D, 52, 2245–2253, (1995). [DOI], [arXiv:grqc/9503003]. (Cited on pages 33, 44, and 45.)
 [63]
Dixon, L.J., Harvey, J.A., Vafa, C. and Witten, E., “Strings on Orbifolds”, Nucl. Phys. B, 261, 678–686, (1985). [DOI]. (Cited on page 78.)
 [64]
Donnelly, H., “Spectrum and the fixed point sets of isometries. I”, Math. Ann., 224, 161–170, (1976). [DOI]. (Cited on page 24.)
 [65]
Donnelly, H., “Asymptotic expansions for the compact quotients of properly discontinuous group actions”, Illinois J. Math., 23, 485–496, (1979). Online version (accessed 20 October 2011): http://projecteuclid.org/euclid.ijm/1256048110. (Cited on page 24.)
 [66]
Donnelly, W., “Entanglement entropy in loop quantum gravity”, Phys. Rev. D, 77, 104006, (2008). [DOI], [arXiv:0802.0880 [grqc]]. (Cited on page 79.)
 [67]
Dou, D., “Comments on the Entanglement Entropy on Fuzzy Spaces”, Mod. Phys. Lett. A, 24, 2467–2480, (2009). [DOI], [arXiv:0903.3731 [grqc]]. (Cited on page 79.)
 [68]
Dou, D. and Ydri, B., “Entanglement entropy on fuzzy spaces”, Phys. Rev. D, 74, 044014, (2006). [DOI], [arXiv:grqc/0605003 [grqc]]. (Cited on page 79.)
 [69]
Dowker, J.S., “Quantum field theory on a cone”, J. Phys. A: Math. Gen., 10, 115–124, (1977). [DOI]. (Cited on pages 11, 12, and 24.)
 [70]
Dowker, J.S., “Effective actions with fixed points”, Phys. Rev. D, 50, 6369–6373, (1994). [DOI], [arXiv:hepth/9406144 [hepth]]. (Cited on pages 30 and 63.)
 [71]
Dowker, J.S., “Remarks on geometric entropy”, Class. Quantum Grav., 11, L55–L60, (1994). [DOI], [arXiv:hepth/9401159]. (Cited on page 25.)
 [72]
Dowker, J.S., “Hyperspherical entanglement entropy”, J. Phys. A: Math. Theor., 43, 445402, (2010). [DOI], [arXiv:1007.3865 [hepth]]. (Cited on page 59.)
 [73]
Dowker, J.S. and Kennedy, G., “Finite temperature and boundary effects in static spacetimes”, J. Phys. A: Math. Gen., 11, 895–920, (1978). [DOI]. (Cited on page 48.)
 [74]
Dowker, J.S. and Schofield, J.P., “High temperature expansion of the free energy of a massive scalar field in a curved space”, Phys. Rev. D, 38, 3327, (1988). [DOI]. (Cited on page 48.)
 [75]
Dowker, J.S. and Schofield, J.P., “Chemical potentials in curved space”, Nucl. Phys. B, 327, 267, (1989). [DOI]. (Cited on page 48.)
 [76]
Duff, M.J., “Twenty years of the Weyl anomaly”, Class. Quantum Grav., 11, 1387–1404, (1994). [DOI], [arXiv:hepth/9308075]. (Cited on page 62.)
 [77]
Dvali, G., “Black holes and large species solution to the hierarchy problem”, Fortschr. Phys., 58, 528–536, (2010). [DOI], [arXiv:0706.2050 [hepth]]. (Cited on pages 76 and 77.)
 [78]
Dvali, G., “Nature of Microscopic Black Holes and Gravity in Theories with Particle Species”, Int. J. Mod. Phys. A, 25, 602–615, (2010). [DOI], [arXiv:0806.3801 [hepth]]. (Cited on page 76.)
 [79]
Dvali, G. and Redi, M., “Black Hole Bound on the Number of Species and Quantum Gravity at LHC”, Phys. Rev. D, 77, 045027, (2008). [DOI], [arXiv:0710.4344 [hepth]]. (Cited on pages 76 and 77.)
 [80]
Dvali, G. and Solodukhin, S.N., “Black Hole Entropy and Gravity Cutoff”, arXiv, eprint, (2008). [arXiv:0806.3976 [hepth]]. (Cited on page 77.)
 [81]
Eisert, J., Cramer, M. and Plenio, M.B., “Colloquium: Area laws for the entanglement entropy”, Rev. Mod. Phys., 82, 277–306, (2010). [DOI], [arXiv:0808.3773 [quantph]]. (Cited on page 6.)
 [82]
Emparan, R., “Heat kernels and thermodynamics in Rindler space”, Phys. Rev. D, 51, 5716–5719, (1995). [DOI], [arXiv:hepth/9407064 [hepth]]. (Cited on page 28.)
 [83]
Emparan, R., “Black hole entropy as entanglement entropy: a holographic derivation”, J. High Energy Phys., 06(2006), 012, (2006). [DOI], [arXiv:hepth/0603081 [hepth]]. (Cited on pages 69, 70, and 76.)
 [84]
Emparan, R., Horowitz, G.T. and Myers, R.C., “Exact description of black holes on branes”, J. High Energy Phys., 2000(01), 007, (2000). [DOI], [arXiv:hepth/9911043 [hepth]]. (Cited on pages 69 and 76.)
 [85]
Fiola, T.M., Preskill, J., Strominger, A. and Trivedi, S.P., “Black hole thermodynamics and information loss in twodimensions”, Phys. Rev. D, 50, 3987–4014, (1994). [DOI], [arXiv:hepth/9403137 [hepth]]. (Cited on page 25.)
 [86]
Frolov, V.P., “Why the entropy of a black hole is A/4?”, Phys. Rev. Lett., 74, 3319–3322, (1995). [DOI], [arXiv:grqc/9406037]. (Cited on page 41.)
 [87]
Frolov, V.P. and Fursaev, D.V., “Mechanism of generation of black hole entropy in Sakharov’s induced gravity”, Phys. Rev. D, 56, 2212–2225, (1997). [DOI], [arXiv:hepth/9703178]. (Cited on page 75.)
 [88]
Frolov, V.P. and Fursaev, D.V., “Plenty of nothing: Black hole entropy in induced gravity”, arXiv, eprint, (1997). [arXiv:hepth/9705207]. (Cited on page 75.)
 [89]
Frolov, V.P. and Fursaev, D., “Black hole entropy in induced gravity: Reduction to 2D quantum field theory on the horizon”, Phys. Rev. D, 58, 124009, (1998). [DOI], [arXiv:hepth/9806078]. (Cited on pages 14 and 75.)
 [90]
Frolov, V.P. and Fursaev, D.V., “Thermal fields, entropy, and black holes”, Class. Quantum Grav., 15, 2041–2074, (1998). [DOI], [arXiv:hepth/9802010]. (Cited on page 6.)
 [91]
Frolov, V.P. and Fursaev, D.V., “Statistical mechanics of charged black holes in induced EinsteinMaxwell gravity”, Phys. Rev. D, 61, 064010, (2000). [DOI], [arXiv:hepth/9910006]. (Cited on page 75.)
 [92]
Frolov, V.P. and Fursaev, D.V., “Statistical mechanics on axiallysymmetric spacetimes with the Killing horizon and entropy of rotating black holes in induced gravity”, Phys. Rev. D, 61, 024007, (2000). [DOI], [arXiv:grqc/9907046]. (Cited on pages 31, 45, and 75.)
 [93]
Frolov, V.P., Fursaev, D., Gegenberg, J. and Kunstatter, G., “Thermodynamics and statistical mechanics of induced Liouville gravity”, Phys. Rev. D, 60, 024016, (1999). [DOI], [arXiv:hepth/9901087]. (Cited on page 75.)
 [94]
Frolov, V.P., Fursaev, D.V. and Zelnikov, A., “Black Hole Entropy: OffShell vs OnShell”, Phys. Rev. D, 54, 2711–2731, (1996). [DOI], [arXiv:hepth/9512184]. (Cited on page 28.)
 [95]
Frolov, V.P., Fursaev, D.V. and Zelnikov, A., “Black Hole Entropy: Thermodynamics, StatisticalMechanics and Subtraction Procedure”, Phys. Lett. B, 382, 220–226, (1996). [DOI], [arXiv:hepth/9603175]. (Cited on page 75.)
 [96]
Frolov, V.P., Fursaev, D.V. and Zelnikov, A., “Black hole statistical mechanics and induced gravity”, Nucl. Phys. B (Proc. Suppl.), 57, 192–196, (1997). [DOI]. (Cited on page 75.)
 [97]
Frolov, V.P., Fursaev, D.V. and Zelnikov, A., “Statistical origin of black hole entropy in induced gravity”, Nucl. Phys. B, 486, 339–352, (1997). [DOI], [arXiv:hepth/9607104]. (Cited on pages 6 and 75.)
 [98]
Frolov, V.P., Israel, W. and Solodukhin, S.N., “On Oneloop Quantum Corrections to the Thermodynamics of Charged Black Holes”, Phys. Rev. D, 54, 2732–2745, (1996). [DOI], [arXiv:hepth/9602105]. (Cited on page 51.)
 [99]
Frolov, V.P. and Novikov, I., “Dynamical origin of the entropy of a black hole”, Phys. Rev. D, 48, 4545–4551, (1993). [DOI], [arXiv:grqc/9309001]. (Cited on page 5.)
 [100]
Frolov, V.P. and Serebryanyi, E.M., “Vacuum polarization in the gravitational field of a cosmic string”, Phys. Rev. D, 35, 3779–3782, (1987). [DOI]. (Cited on page 24.)
 [101]
Fursaev, D.V., “Spectral geometry and one loop divergences on manifolds with conical singularities”, Phys. Lett. B, 334, 53–60, (1994). [DOI], [arXiv:hepth/9405143]. (Cited on pages 11, 12, 24, and 62.)
 [102]
Fursaev, D.V., “Black hole thermodynamics and renormalization”, Mod. Phys. Lett. A, 10, 649–656, (1995). [DOI], [arXiv:hepth/9408066]. (Cited on page 28.)
 [103]
Fursaev, D.V., “Temperature and entropy of a quantum black hole and conformal anomaly”, Phys. Rev. D, 51, 5352–5355, (1995). [DOI], [arXiv:hepth/9412161]. (Cited on pages 63 and 66.)
 [104]
Fursaev, D.V., “Euclidean and canonical formulations of statistical mechanics in the presence of Killing horizons”, Nucl. Phys. B, 524, 447–468, (1998). [DOI], [arXiv:hepth/9709213]. (Cited on pages 28, 41, 45, and 49.)
 [105]
Fursaev, D.V., “Energy, Hamiltonian, Noether charge, and black holes”, Phys. Rev. D, 59, 064020, (1999). [DOI], [arXiv:hepth/9809049]. (Cited on page 75.)
 [106]
Fursaev, D.V., “Black hole entropy in induced gravity and information loss”, Nucl. Phys. B (Proc. Suppl.), 88, 277–280, (2000). [DOI], [arXiv:grqc/9911065]. (Cited on page 75.)
 [107]
Fursaev, D.V., “Black Hole Thermodynamics, Induced Gravity and Gravity in Brane Worlds”, arXiv, eprint, (2000). [arXiv:hepth/0009164]. (Cited on page 76.)
 [108]
Fursaev, D.V., “Entanglement entropy in critical phenomena and analogue models of quantum gravity”, Phys. Rev. D, 73, 124025, (2006). [DOI], [arXiv:hepth/0602134]. (Cited on page 16.)
 [109]
Fursaev, D.V., “Proof of the holographic formula for entanglement entropy”, J. High Energy Phys., 2006(09), 018, (2006). [DOI], [arXiv:hepth/0606184]. (Cited on page 68.)
 [110]
Fursaev, D.V. and Miele, G., “Cones, Spins and Heat Kernels”, Nucl. Phys. B, 484, 697–723, (1997). [DOI], [arXiv:hepth/9605153]. (Cited on pages 28 and 35.)
 [111]
Fursaev, D.V. and Solodukhin, S.N., “On the description of the Riemannian geometry in the presence of conical defects”, Phys. Rev. D, 52, 2133–2143, (1995). [DOI], [arXiv:hepth/9501127]. (Cited on pages 5, 21, 22, 23, and 62.)
 [112]
Fursaev, D.V. and Solodukhin, S.N., “On oneloop renormalization of black hole entropy”, Phys. Lett. B, 365, 51–55, (1996). [DOI], [arXiv:hepth/9412020]. (Cited on pages 33, 50, and 62.)
 [113]
Garay, L.J., “Quantum gravity and minimum length”, Int. J. Mod. Phys. A, 10, 145–166, (1995). [DOI], [arXiv:grqc/9403008 [grqc]]. (Cited on page 80.)
 [114]
Ghosh, A. and Mitra, P., “Entropy in dilatonic black hole background”, Phys. Rev. Lett., 73, 2521–2523, (1994). [DOI], [arXiv:hepth/9406210 [hepth]]. (Cited on page 28.)
 [115]
Ghosh, A. and Mitra, P., “Entropy for extremal ReissnerNordstrom black holes”, Phys. Lett. B, 357, 295–299, (1995). [DOI], [arXiv:hepth/9411128 [hepth]]. (Cited on page 28.)
 [116]
Ghosh, A. and Mitra, P., “Quantum correction to the thermodynamical entropy of a black hole”, Phys. Rev. D, 56, 3568–3570, (1997). [DOI], [arXiv:grqc/9706054 [grqc]]. Erratum: 10.1103/PhysRevD.57.1316. (Cited on page 28.)
 [117]
Ghosh, K., “Entropy of scalar fields in (3+1)dimensional constant curvature black hole background”, Phys. Rev. D, 60, 104003, (1999). [DOI], [arXiv:grqc/9809063 [grqc]]. (Cited on page 28.)
 [118]
Ghosh, K., “Entropy of scalar field in (3+1)dimensional ReissnerNordstrom de Sitter black hole background”, arXiv, eprint, (2002). [arXiv:grqc/0212060 [grqc]]. (Cited on page 28.)
 [119]
Gibbons, G.W. and Hawking, S.W., “Action Integrals and Partition Functions in Quantum Gravity”, Phys. Rev. D, 15, 2752–2756, (1977). [DOI]. (Cited on page 40.)
 [120]
Gibbons, G.W. and Maeda, K., “Black holes and membranes in higherdimensional theories with dilaton fields”, Nucl. Phys. B, 298, 741–775, (1988). [DOI]. (Cited on page 27.)
 [121]
Gibbons, G.W. and Perry, M.J., “Black holes and thermal Green functions”, Proc. R. Soc. London, Ser. A, 358, 467–494, (1978). (Cited on page 20.)
 [122]
Graham, C.R. and Hirachi, V., “The Ambient Obstruction Tensor and QCurvature”, arXiv, eprint, (2004). [arXiv:math.DG/0405068]. (Cited on page 72.)
 [123]
Graham, C.R. and Witten, E., “Conformal anomaly of submanifold observables in AdS/CFT correspondence”, Nucl. Phys. B, 546, 52–64, (1999). [DOI], [arXiv:hepth/9901021]. (Cited on page 70.)
 [124]
Gubser, S.S., “AdS/CFT and gravity”, Phys. Rev. D, 63, 084017, (2001). [DOI], [arXiv:hepth/9912001 [hepth]]. (Cited on page 76.)
 [125]
Gubser, S.S., Klebanov, I.R. and Polyakov, A.M., “Gauge theory correlators from noncritical string theory”, Phys. Lett. B, 428, 105–114, (1998). [DOI], [arXiv:hepth/9802109 [hepth]]. (Cited on page 68.)
 [126]
Hartle, J.B. and Hawking, S.W., “Pathintegral derivation of black hole radiance”, Phys. Rev. D, 13, 2188–2203, (1976). [DOI]. (Cited on pages 19 and 20.)
 [127]
Hartle, J.B. and Hawking, S.W., “Wave Function of the Universe”, Phys. Rev. D, 28, 2960–2975, (1983). [DOI]. (Cited on page 19.)
 [128]
Hawking, S.W., “Particle creation by black holes”, Commun. Math. Phys., 43, 199–220, (1975). [DOI]. (Cited on page 5.)
 [129]
Hawking, S.W., “Breakdown of Predictability in Gravitational Collapse”, Phys. Rev. D, 14, 2460–2473, (1976). [DOI]. (Cited on page 80.)
 [130]
Hawking, S., Maldacena, J.M. and Strominger, A., “DeSitter entropy, quantum entanglement and ADS/CFT”, J. High Energy Phys., 2001(05), 001, (2001). [DOI], [arXiv:hepth/0002145 [hepth]]. (Cited on pages 69, 75, and 76.)
 [131]
Henningson, M. and Skenderis, K., “The Holographic Weyl anomaly”, J. High Energy Phys., 1998(07), 023, (1998). [DOI], [arXiv:hepth/9806087 [hepth]]. (Cited on pages 63, 65, 71, and 72.)
 [132]
Ho, J. and Kang, G., “Superradiance and the statistical mechanical entropy of rotating BTZ black holes”, Phys. Lett. B, 445, 27–35, (1998). [DOI], [arXiv:grqc/9806118 [grqc]]. (Cited on page 45.)
 [133]
Holzhey, C., Larsen, F. and Wilczek, F., “Geometric and renormalized entropy in conformal field theory”, Nucl. Phys. B, 424, 443–467, (1994). [DOI], [arXiv:hepth/9403108]. (Cited on page 9.)
 [134]
Hung, L.Y., Myers, R.C. and Smolkin, M., “On Holographic Entanglement Entropy and Higher Curvature Gravity”, J. High Energy Phys., 2011(4), 025, (2011). [DOI], [arXiv:1101.5813 [hepth]]. (Cited on pages 20 and 65.)
 [135]
Iellici, D. and Moreira Jr, E.S., “Ambiguity in the evaluation of effective action on the cone”, Phys. Rev. D, 60, 124015, (1999). [DOI], [arXiv:hepth/9807184 [hepth]]. (Cited on page 28.)
 [136]
Iellici, D. and Moretti, V., “Thermal partition function of photons and gravitons in a Rindler wedge”, Phys. Rev. D, 54, 7459–7469, (1996). [DOI], [arXiv:hepth/9607015 [hepth]]. (Cited on page 49.)
 [137]
Iellici, D. and Moretti, V., “Kabat’s Surface Terms in the ZetaFunction Approach”, in Bassan, M., Ferrari, V., Francaviglia, M., Fucito, F. and Modena, I., eds., General Relativity and Gravitational Physics, Proceedings of the 12th Italian Conference, Rome, Italy, September 23–27, 1996, pp. 317–321, (World Scientific, Singapore; River Edge, NJ, 1997). [arXiv:hepth/9703088]. (Cited on page 28.)
 [138]
Israel, W., “Thermofield dynamics of black holes”, Phys. Lett. A, 57, 107–110, (1976). [DOI]. (Cited on page 21.)
 [139]
Iwashita, Y., Kobayashi, T., Shiromizu, T. and Yoshino, H., “Holographic entanglement entropy of de Sitter braneworld”, Phys. Rev. D, 74, 064027, (2006). [DOI], [arXiv:hepth/0606027 [hepth]]. (Cited on pages 69 and 70.)
 [140]
Iyer, V. and Wald, R.M., “Comparison of Noether charge and Euclidean methods for computing the entropy of stationary black holes”, Phys. Rev. D, 52, 4430–4439, (1995). [DOI], [arXiv:grqc/9503052 [grqc]]. (Cited on page 32.)
 [141]
Jacobson, T., “Black Hole Entropy and Induced Gravity”, arXiv, eprint, (1994). [arXiv:grqc/9404039]. (Cited on pages 5, 33, and 74.)
 [142]
Jacobson, T., “A Note on HartleHawking vacua”, Phys. Rev. D, 50, 6031–6032, (1994). [DOI], [arXiv:grqc/9407022 [grqc]]. (Cited on page 21.)
 [143]
Jacobson, T.A., “Renormalization and black hole entropy in Loop Quantum Gravity”, Class. Quantum Grav., 24, 4875–4879, (2007). [DOI], [arXiv:0707.4026 [grqc]]. (Cited on page 79.)
 [144]
Jacobson, T., Kang, G. and Myers, R.C., “On black hole entropy”, Phys. Rev. D, 49, 6587–6598, (1994). [DOI], [arXiv:grqc/9312023]. (Cited on pages 25 and 32.)
 [145]
Jacobson, T. and Myers, R.C., “Black hole entropy and higher curvature interactions”, Phys. Rev. Lett., 70, 3684–3687, (1993). [DOI], [arXiv:hepth/9305016]. (Cited on page 23.)
 [146]
Jacobson, T. and Parentani, R., “Horizon Entropy”, Found. Phys., 33, 323–348, (2003). [DOI], [arXiv:grqc/0302099]. (Cited on page 6.)
 [147]
Jacobson, T. and Parentani, R., “Black hole entanglement entropy regularized in a freely falling frame”, Phys. Rev. D, 76, 024006, (2007). [DOI], [arXiv:hepth/0703233]. (Cited on page 80.)
 [148]
Jing, J. and Yan, M.L., “Entropies of the general nonextreme stationary axisymmetric black hole: Statistical mechanics and thermodynamics”, Phys. Rev. D, 61, 044016, (2000). [DOI], [arXiv:grqc/9907011 [grqc]]. (Cited on page 45.)
 [149]
Jing, J. and Yan, M.L., “Quantum entropy of the Kerr black hole arising from gravitational perturbation”, Phys. Rev. D, 64, 064015, (2001). [DOI], [arXiv:grqc/0104054 [grqc]]. (Cited on page 45.)
 [150]
Kabat, D.N., “Black hole entropy and entropy of entanglement”, Nucl. Phys. B, 453, 281–302, (1995). [DOI], [arXiv:hepth/9503016]. (Cited on page 35.)
 [151]
Kabat, D.N., Shenker, S.H. and Strassler, M.J., “Black hole entropy in the O(N) model”, Phys. Rev. D, 52, 7027–7036, (1995). [DOI], [arXiv:hepth/9506182]. (Cited on pages 39 and 78.)
 [152]
Kabat, D.N. and Strassler, M.J., “A Comment on entropy and area”, Phys. Lett. B, 329, 46–52, (1994). [DOI], [arXiv:hepth/9401125]. (Cited on pages 21 and 25.)
 [153]
Kenmoku, M., Ishimoto, K., Nandi, Kamal Kanti and Shigemoto, K., “Scalar field contribution to rotating black hole entropy”, Phys. Rev. D, 73, 064004, (2006). [DOI], [arXiv:grqc/0510012 [grqc]]. (Cited on page 45.)
 [154]
Kim, W., Kim, Y.W. and Park, Y.J., “Entropy of 2+1 de Sitter space with the GUP”, J. Korean Phys. Soc., 49, 1360–1364, (2006). [arXiv:grqc/0604065 [grqc]]. (Cited on page 80.)
 [155]
Kim, W., Park, Y.J., Son, E.J. and Yoon, M.S., “Statistical Entropy and Superradiance in 2+1 Dimensional Acoustic Black Holes”, J. Korean Phys. Soc., 49, 15–20, (2006). [DOI], [arXiv:grqc/0504127 [grqc]]. (Cited on page 45.)
 [156]
Kim, Y.W. and Park, Y.J., “Entropy of the Schwarzschild black hole to all orders in the Planck length”, Phys. Lett. B, 655, 172–177, (2007). [DOI], [arXiv:0707.2128 [grqc]]. (Cited on page 80.)
 [157]
Klebanov, I.R., Kutasov, D. and Murugan, A., “Entanglement as a Probe of Confinement”, Nucl. Phys. B, 796, 274–293, (2008). [DOI], [arXiv:0709.2140 [hepth]]. (Cited on page 69.)
 [158]
Koksma, J.F., Prokopec, T. and Schmidt, M.G., “Entropy and correlators in quantum field theory”, Ann. Phys. (N.Y.), 325, 1277–1303, (2010). [DOI], [arXiv:1002.0749 [hepth]]. (Cited on page 80.)
 [159]
Koksma, J.F., Prokopec, T. and Schmidt, M.G., “Decoherence and Dynamical Entropy Generation in Quantum Field Theory”, arXiv, eprint, (2011). [arXiv:1101.5323 [quantph]]. (Cited on page 80.)
 [160]
Koksma, J.F., Prokopec, T. and Schmidt, M.G., “Decoherence in an Interacting Quantum Field Theory: Thermal Case”, Phys. Rev. D, 83, 085011, (2011). [DOI], [arXiv:1102.4713 [hepth]]. (Cited on page 80.)
 [161]
Korepin, V.E., “Universality of Entropy Scaling in One Dimensional Gapless Models”, Phys. Rev. Lett., 92, 096402, (2004). [DOI]. (Cited on page 9.)
 [162]
Larsen, F. and Wilczek, F., “Geometric Entropy, Wave Functionals, and Fermions”, Ann. Phys. (N.Y.), 243, 280–298, (1995). [DOI], [arXiv:hepth/9408089]. (Cited on page 35.)
 [163]
Larsen, F. and Wilczek, F., “Renormalization of black hole entropy and of the gravitational coupling constant”, Nucl. Phys. B, 458, 249–266, (1996). [DOI], [arXiv:hepth/9506066]. (Cited on pages 35 and 37.)
 [164]
Latorre, J.I., Lütken, C.A., Rico, E. and Vidal, G., “Finegrained entanglement loss along renormalizationgroup flows”, Phys. Rev. A, 71, 034301, (2005). [DOI], [arXiv:quantph/0404120 [quantph]]. (Cited on page 80.)
 [165]
Livine, E.R. and Terno, D.R., “Quantum black holes: Entropy and entanglement on the horizon”, Nucl. Phys. B, 741, 131–161, (2006). [DOI], [arXiv:grqc/0508085]. (Cited on page 79.)
 [166]
Lovelock, D., “The Einstein tensor and its generalizations”, J. Math. Phys., 12, 498–501, (1971). [DOI]. (Cited on page 22.)
 [167]
Maldacena, J.M., “The Large Limit of Superconformal Field Theories and Supergravity”, Adv. Theor. Math. Phys., 2, 231–252, (1998). [arXiv:hepth/9711200 [hepth]]. (Cited on pages 63, 68, and 76.)
 [168]
Maldacena, J.M., “Eternal black holes in antide Sitter”, J. High Energy Phys., 2003(04), 021, (2003). [DOI], [arXiv:hepth/0106112]. (Cited on page 78.)
 [169]
Mandal, G., Sengupta, A.M. and Wadia, S.R., “Classical solutions of 2dimensional string theory”, Mod. Phys. Lett. A, 6, 1685–1692, (1991). [DOI]. (Cited on page 51.)
 [170]
Mann, R.B. and Solodukhin, S.N., “Conical geometry and quantum entropy of a charged Kerr black hole”, Phys. Rev. D, 54, 3932–3940, (1996). [DOI], [arXiv:hepth/9604118]. (Cited on pages 28, 29, 30, 31, and 33.)
 [171]
Mann, R.B. and Solodukhin, S.N., “Quantum scalar field on threedimensional (BTZ) black hole instanton: Heat kernel, effective action and thermodynamics”, Phys. Rev. D, 55, 3622–3632, (1997). [DOI], [arXiv:hepth/9609085]. (Cited on pages 54 and 55.)
 [172]
Mann, R.B. and Solodukhin, S.N., “Universality of quantum entropy for extreme black holes”, Nucl. Phys. B, 523, 293–307, (1998). [DOI], [arXiv:hepth/9709064]. (Cited on pages 14, 55, 57, and 58.)
 [173]
Mann, R.B., Tarasov, L. and Zelnikov, A., “Brick walls for black holes”, Class. Quantum Grav., 9, 1487–1494, (1992). [DOI]. (Cited on page 44.)
 [174]
McKean Jr, H.P. and Singer, I.M., “Curvature and the eigenvalues of the Laplacian”, J. Differ. Geom., 1(1–2), 43–69, (1967). (Cited on page 24.)
 [175]
Metlitski, M.A., Fuertes, C.A. and Sachdev, S., “Entanglement Entropy in the O(N) model”, Phys. Rev. B, 80, 115122, (2009). [DOI], [arXiv:0904.4477 [condmat.statmech]]. (Cited on page 39.)
 [176]
Moretti, V., “Geometric entropy and curvature coupling in conical spaces: Zeta function approach”, Class. Quantum Grav., 14, L123–L127, (1997). [DOI], [arXiv:hepth/9701099 [hepth]]. (Cited on page 28.)
 [177]
Moretti, V. and Iellici, D., “Optical approach for the thermal partition function of photons”, Phys. Rev. D, 55, 3552–3563, (1997). [DOI], [arXiv:hepth/9610180 [hepth]]. (Cited on page 49.)
 [178]
Mukohyama, S. and Israel, W., “Black holes, brick walls and the Boulware state”, Phys. Rev. D, 58, 104005, (1998). [DOI], [arXiv:grqc/9806012 [grqc]]. (Cited on page 44.)
 [179]
Myers, R.C., “Black hole entropy in twodimensions”, Phys. Rev. D, 50, 6412–6421, (1994). [DOI], [arXiv:hepth/9405162 [hepth]]. (Cited on page 50.)
 [180]
Myers, R.C. and Perry, M.J., “Black Holes in Higher Dimensional SpaceTimes”, Ann. Phys. (N.Y.), 172, 304, (1986). [DOI]. (Cited on pages 60 and 65.)
 [181]
Myers, R.C. and Sinha, A., “Holographic ctheorems in arbitrary dimensions”, J. High Energy Phys., 2011(01), 125, (2011). [DOI], [arXiv:1011.5819 [hepth]]. (Cited on pages 24 and 65.)
 [182]
Nelson, B.L. and Panangaden, P., “Scaling behavior of interacting quantum fields in curved spacetime”, Phys. Rev. D, 25, 1019–1027, (1982). [DOI]. (Cited on page 39.)
 [183]
Nesterov, D. and Solodukhin, S.N., “Shortdistance regularity of Green’s function and UV divergences in entanglement entropy”, J. High Energy Phys., 2010(09), 041, (2010). [DOI], [arXiv:1008.0777 [hepth]]. (Cited on pages 14, 15, and 80.)
 [184]
Nesterov, D. and Solodukhin, S.N., “Gravitational effective action and entanglement entropy in UV modified theories with and without Lorentz symmetry”, Nucl. Phys. B, 842, 141–171, (2011). [DOI], [arXiv:1007.1246 [hepth]]. (Cited on pages 14, 15, 16, 20, 33, 34, and 80.)
 [185]
Nishioka, T., Ryu, S. and Takayanagi, Tadashi, “Holographic entanglement entropy: an overview”, J. Phys. A: Math. Theor., 42, 504008, (2009). [DOI], [arXiv:0905.0932 [hepth]]. (Cited on pages 6 and 68.)
 [186]
Padmanabhan, T., “Finite entanglement entropy from the zeropointarea of spacetime”, Phys. Rev. D, 82, 124025, (2010). [DOI], [arXiv:1007.5066 [grqc]]. (Cited on page 80.)
 [187]
Randall, L. and Sundrum, R., “An alternative to compactification”, Phys. Rev. Lett., 83, 4690–4693, (1999). [DOI], [arXiv:hepth/9906064 [hepth]]. (Cited on page 76.)
 [188]
Ryu, S. and Takayanagi, T., “Aspects of Holographic Entanglement Entropy”, J. High Energy Phys., 2006(08), 045, (2006). [DOI], [arXiv:hepth/0605073 [hepth]]. (Cited on pages 62 and 68.)
 [189]
Ryu, S. and Takayanagi, T., “Holographic derivation of entanglement entropy from AdS/CFT”, Phys. Rev. Lett., 96, 181602, (2006). [DOI], [arXiv:hepth/0603001 [hepth]]. (Cited on pages 6 and 68.)
 [190]
Sakharov, A.D., “Vacuum quantum fluctuations in curved space and the theory of gravitation”, Sov. Phys. Dokl., 12, 1040–1041, (1968). (Cited on page 74.)
 [191]
Schwimmer, A. and Theisen, S., “Entanglement Entropy, Trace Anomalies and Holography”, Nucl. Phys. B, 801, 1–24, (2008). [DOI], [arXiv:0802.1017 [hepth]]. (Cited on page 70.)
 [192]
Sen, A., “State Operator Correspondence and Entanglement in AdS_{2}/CFT_{1}”, arXiv, eprint, (2011). [arXiv:1101.4254 [hepth]]. (Cited on page 78.)
 [193]
Skenderis, K., “Asymptotically Antide Sitter spacetimes and their stress energy tensor”, Int. J. Mod. Phys. A, 16, 740–749, (2001). [DOI], [arXiv:hepth/0010138 [hepth]]. (Cited on page 68.)
 [194]
Skenderis, K. and Solodukhin, S.N., “Quantum effective action from the AdS/CFT correspondence”, Phys. Lett. B, 472, 316–322, (2000). [DOI], [arXiv:hepth/9910023 [hepth]]. (Cited on pages 69 and 70.)
 [195]
Sokolov, D.D. and Starobinsky, A.A., “The structure of the curvature tensor at conical singularities”, Sov. Phys. Dokl., 22, 312–313, (1977). (Cited on page 21.)
 [196]
Solodukhin, S.N., “The conical singularity and quantum corrections to entropy of black hole”, Phys. Rev. D, 51, 609–617, (1995). [DOI], [arXiv:hepth/9407001]. (Cited on pages 5, 25, 27, 28, 33, 50, and 62.)
 [197]
Solodukhin, S.N., “On ‘Nongeometric’ contribution to the entropy of black hole due to quantum corrections”, Phys. Rev. D, 51, 618–621, (1995). [DOI], [arXiv:hepth/9408068]. (Cited on pages 5, 26, 27, and 28.)
 [198]
Solodukhin, S.N., “One loop renormalization of black hole entropy due to nonminimally coupled matter”, Phys. Rev. D, 52, 7046–7052, (1995). [DOI], [arXiv:hepth/9504022]. (Cited on pages 37 and 40.)
 [199]
Solodukhin, S.N., “Black hole entropy: statistical mechanics agrees thermodynamics”, Phys. Rev. D, 54, 3900–3903, (1996). [DOI], [arXiv:hepth/9601154]. (Cited on pages 41, 49, and 51.)
 [200]
Solodukhin, S.N., “Nonminimal coupling and quantum entropy of black hole”, Phys. Rev. D, 56, 4968–4974, (1997). [DOI], [arXiv:hepth/9612061]. (Cited on page 45.)
 [201]
Solodukhin, S.N., “Entropy of Schwarzschild black hole and stringblack hole correspondence”, Phys. Rev. D, 57, 2410–2414, (1998). [DOI], [arXiv:hepth/9701106]. (Cited on page 66.)
 [202]
Solodukhin, S.N., “Entanglement entropy of black holes and AdS/CFT correspondence”, Phys. Rev. Lett., 97, 201601, (2006). [DOI], [arXiv:hepth/0606205 [hepth]]. (Cited on pages 69, 70, 73, and 76.)
 [203]
Solodukhin, S.N., “Entanglement entropy and the Ricci flow”, Phys. Lett. B, 646, 268–274, (2007). [DOI], [arXiv:hepth/0609045]. (Cited on page 80.)
 [204]
Solodukhin, S.N., “Entanglement entropy, conformal invariance and extrinsic geometry”, Phys. Lett. B, 665, 305–309, (2008). [DOI], [arXiv:0802.3117 [hepth]]. (Cited on pages 16, 20, and 63.)
 [205]
Solodukhin, S.N., “Entanglement Entropy in NonRelativistic Field Theories”, J. High Energy Phys., 2010(04), 101, (2010). [DOI], [arXiv:0909.0277 [hepth]]. (Cited on pages 16 and 39.)
 [206]
Solodukhin, S.N., “Entanglement entropy of round spheres”, Phys. Lett. B, 693, 605–608, (2010). [DOI], [arXiv:1008.4314 [hepth]]. (Cited on pages 57 and 59.)
 [207]
Sommerfeld, A., “Über verzweigte Potentiale im Raum”, Proc. London Math. Soc. (1), XXVIII, 395–429, (1897). [DOI]. (Cited on page 11.)
 [208]
Srednicki, M., “Entropy and area”, Phys. Rev. Lett., 71, 666–669, (1993). [DOI], [arXiv:hepth/9303048]. (Cited on pages 5 and 9.)
 [209]
Strominger, A. and Vafa, C., “Microscopic Origin of the BekensteinHawking Entropy”, Phys. Lett. B, 379, 99–104, (1996). [DOI], [arXiv:hepth/9601029]. (Cited on page 6.)
 [210]
Sun, X. and Liu, W., “Improved black hole entropy calculation without cutoff”, Mod. Phys. Lett. A, 19, 677–680, (2004). [DOI]. (Cited on page 80.)
 [211]
Susskind, L., “Some Speculations about Black Hole Entropy in String Theory”, arXiv, eprint, (1993). [arXiv:hepth/9309145]. (Cited on page 5.)
 [212]
Susskind, L., “The world as a hologram”, J. Math. Phys., 36, 6377–6396, (1995). [DOI], [arXiv:hepth/9409089 [hepth]]. (Cited on page 68.)
 [213]
Susskind, L. and Uglum, J., “Black hole entropy in canonical quantum gravity and superstring theory”, Phys. Rev. D, 50, 2700–2711, (1994). [DOI], [arXiv:hepth/9401070]. (Cited on pages 5, 33, 44, and 78.)
 [214]
’t Hooft, G., “On the Quantum Structure of a Black Hole”, Nucl. Phys. B, 256, 727–745, (1985). [DOI]. (Cited on pages 5, 41, 44, and 51.)
 [215]
’t Hooft, G., “Dimensional Reduction in Quantum Gravity”, arXiv, eprint, (1993). [arXiv:grqc/9310026]. (Cited on page 68.)
 [216]
’t Hooft, G., “Quantum gravity without spacetime singularities or horizons”, arXiv, eprint, (2009). [arXiv:0909.3426 [grqc]]. (Cited on page 66.)
 [217]
’t Hooft, G., “The Conformal Constraint in Canonical Quantum Gravity”, arXiv, eprint, (2010). [arXiv:1011.0061 [grqc]]. (Cited on page 66.)
 [218]
’t Hooft, G., “Probing the small distance structure of canonical quantum gravity using the conformal group”, arXiv, eprint, (2010). [arXiv:1009.0669 [grqc]]. (Cited on page 66.)
 [219]
Vassilevich, D.V., “Heat kernel expansion: User’s manual”, Phys. Rep., 388, 279–360, (2003). [DOI], [arXiv:hepth/0306138 [hepth]]. (Cited on pages 24, 58, and 61.)
 [220]
Visser, M., “Sakharov’s induced gravity: A modern perspective”, Mod. Phys. Lett. A, 17, 977–992, (2002). [DOI], [arXiv:grqc/0204062]. (Cited on page 74.)
 [221]
Wald, R.M., “Black hole entropy is the Noether charge”, Phys. Rev. D, 48, R3427–R3431, (1993). [DOI], [arXiv:grqc/9307038]. (Cited on pages 25, 32, and 50.)
 [222]
Witten, E., “On string theory and black holes”, Phys. Rev. D, 44, 314–324, (1991). [DOI]. (Cited on page 51.)
 [223]
Witten, E., “Antide Sitter space and holography”, Adv. Theor. Math. Phys., 2, 253–291, (1998). [arXiv:hepth/9802150 [hepth]]. (Cited on page 68.)
 [224]
Wu, S.Q. and Yan, M.L., “Entropy of Kerrde Sitter black hole due to arbitrary spin fields”, Phys. Rev. D, 69, 044019, (2004). [DOI], [arXiv:grqc/0303076 [grqc]]. (Cited on page 45.)
 [225]
Yoon, M., Ha, J. and Kim, W., “Entropy of ReissnerNordstrom black holes with minimal length revisited”, Phys. Rev. D, 76, 047501, (2007). [DOI], [arXiv:0706.0364 [grqc]]. (Cited on page 80.)
 [226]
Zaslavsky, O.B., “Geometry of nonextreme black holes near the extreme state”, Phys. Rev. D, 56, 2188–2191, (1997). [DOI], [arXiv:grqc/9707015 [grqc]]. (Cited on page 55.)
 [227]
Zerbini, S., Cognola, G. and Vanzo, L., “Euclidean approach to the entropy for a scalar field in Rindlerlike spacetimes”, Phys. Rev. D, 54, 2699–2710, (1996). [DOI], [arXiv:hepth/9603106 [hepth]]. (Cited on page 28.)
Acknowledgments
I would like to thank my collaborators and friends A. Barvinsky, G. Dvali, G. ’t Hooft, W. Israel, V. Frolov, D. Fursaev, R. Mann, R. Myers, S. Mukhanov, D. Nesterov and A. Zelnikov for fruitful collaboration and many inspiring discussions. I thank G. ’t Hooft, R. Myers and R. Mann for reading this review and giving valuable remarks. I am indebted to my colleagues N. Mohammedi, L. Villain and, especially, to S. Nicolis for reading this review and their many comments, which helped to improve the presentation. As always, the writing of this review would not be possible without the constant support of my family.