Counterterms for Static Lovelock Solutions

In this paper, we introduce the counterterms that remove the non-logarithmic divergences of the action in third order Lovelock gravity for static spacetimes. We do this by defining the cosmological constant in such a way that the asymptotic form of the metric have the same form in Lovelock and Einstein gravities. Thus, we employ the counterterms of Einstein gravity and show that the power law divergences of the action of Lovelock gravity for static spacetimes can be removed by suitable choice of coefficients. We find that the dependence of these coefficients on the dimension in Lovelock gravity is the same as in Einstein gravity. We also introduce the finite energy-momentum tensor and employ these counterterms to calculate the finite action and mass of static black hole solutions of third order Lovelock gravity. Next, we calculate the thermodynamic quantities and show that the entropy calculated through the use of Gibbs-Duhem relation is consistent with the obtained entropy by Wald's formula. Furthermore, we find that in contrast to Einstein gravity in which there exists no uncharged extreme black hole, third order Lovelock gravity can have these kind of black holes. Finally, we investigate the stability of static charged black holes of Lovelock gravity in canonical ensemble and find that small black holes show a phase transition between very small and small black holes, while the large ones are stable.

theory [3]. The most natural extension of general relativity with higher curvature terms and with the assumption of Einstein -that the left hand side of the field equations is the most general symmetric conserved tensor containing no more than two-derivatives of the metric -is Lovelock theory. Lovelock [4] found the most general symmetric conserved tensor satisfying this property.
The resultant tensor is nonlinear in the Riemann tensor and differs from the Einstein tensor only if the spacetime has more than 4 dimensions. Although Lovelock gravity leads to second-order field equations and it has ghost free AdS solution [5], it has been recently shown that quadratic and cubic gravities entail causality violation and there are stringent conditions on the coupling constants [6].
The problem with the total action of Einstein gravity is that it is divergent when evaluated on the solutions [7][8][9]. Due to this fact, all the other conserved quantities which is calculated through the use of this action is also divergent. One way of eliminating these divergences is through the use of background subtraction method of Brown and York [7]. In this method, the boundary surface is embedded in another (background) spacetime, and one subtracts the action evaluated on the embedded surface of the background spacetime from the total action. Such a procedure causes the resulting physical quantities to depend on the choice of reference background. Furthermore, it is not possible in general to embed the boundary surface into a background spacetime. For asymptotically AdS solutions of Einstein gravity, one may remove the non-logarithmic divergences in the action by adding a counterterm action which is a functional of the boundary curvature invariants [10,11]. Indeed, this counterterm method furnishes a means for calculating the action and conserved quantities intrinsically without reliance on any reference spacetime [12][13][14]. Although there may exist a very large number of possible invariants, only a finite number of them are non-vanishing in a given dimension on a boundary at infinity. This method has been applied to many cases such as black holes with rotation, NUT charge, various topologies, rotating black strings with zero curvature horizons and rotating higher genus black branes [15]. Although the counterterm method applies for the case of a specially infinite boundary, it was also employed for the computation of the conserved and thermodynamic quantities in the case of a finite boundary [16].
All of the works mentioned in the previous paragraph were limited to Einstein gravity. Although the counterterm of Lovelock gravity with flat horizon has been introduced [17,18], only a few works related to the counterterm method have been done for Lovelock gravity with curved horizon. This is due to the fact that even for Einstein gravity, the systematic construction that provides the form of the counterterms becomes cumbersome for high enough dimensions [10,11]. Indeed in this method, one should reconstruct the spacetime metric by solving iteratively the field equations in the Fefferman-Graham frame [19]: where g (0)ij (x) is the boundary data of an initial-value problem governed by the equations of motion. However, even for Einstein-Hilbert theory, solving the coefficients g (p) (x) in Eq. (1) as covariant functionals of g (0) (x) is only possible for low enough dimensions. Thus, it is expected that the holographic renormalization procedure would be even more complicated in Lovelock gravity because of the nonlinearity of field equations. Indeed, because of the nonlinearity of the field equations solving g (p) (x) in Eq. (1) as covariant functionals of g (0) (x) would be even more cumbersome. So, the authors in Ref. [20] presented an alternative construction of Kounterterms. Instead of adding counterterms to cancel the divergence at the boundary explained above, they circumvented the difficulties of the standard method by using Kounterterms which depend on the intrinsic and the extrinsic curvatures of the boundary. They selected the Kounterterms as the boundary terms which are regular on the asymptotic region. Indeed, the regularization process is encoded in the boundary terms already presented and there is no need to add further counterterms.
The exact rotating solutions of Lovelock gravity with curved horizon are not introduced till now. Indeed, only static solutions of Lovelock gravity with different matter fields are known [21].
So, because of the difficulties of the holographic renormalization procedure in Lovelock gravity and the nonexistance of an exact rotating solution of this theory, we limit ourselves to the case of counterterms of Lovelock gravity for static solutions. The counterterms of asymptotically AdS static solutions of Gauss-Bonnet gravity have been introduced in Ref. [22,23]. Also, the finite action and global charges of asymptotically de Sitter static solutions has been obtained in Ref. [24].
Here we like to apply the counterterm method to the case of the static solutions of the field equations of third order Lovelock gravity with curved horizon. We define the cosmological constant in such a way that the maximally symmetric AdS spacetime be the vacuum solution of Lovelock gravity. In Eq. (2) dΣ 2 k,n−1 is the metric of an (n − 1)dimensional maximally symmetric space with curvature constant (n − 1)(n − 2)k and volume V k,n−1 . Indeed, this choice of cosmological constant makes the asymptotic form of the solutions of Lovelock gravity to be exactly the same as that of Einstein gravity. Thus, we expect that the counterterm introduced for Einstein gravity in [10] may remove the power law divergences in the action of Lovelock gravity. Although the counterterms which should be added to Gauss-Bonnet gravity in order to remove the power law divergences of the action for static solutions are introduced in Ref. [22], they depend on the Gauss-Bonnet coefficient. However, because of our choice of the cosmological constant, our counterterms are the same as those of Einstein gravity and are independent of Lovelock coefficients. In order to check our counterterms, we calculate the finite action and the mass of the black hole through the use of counterterm method. Then, we use these finite quantities and the Gibbs-Duhem relation to obtain the entropy. We find that the calculated entropy of the black holes is consistent with the Wald's formula [25]. As another test of our counterterm method, we show that the mass obtained through the use of counterterm method satisfies the first law of thermodynamics. We, also, perform a stability analysis of the black hole solutions in canonical ensemble and investigate the effects of third order Lovelock term on the stability. This paper is organized as follows. In section II, we review the well-defined action of Lovelock gravity. In section III, we introduce the counterterms for third order Lovelock gravity for static spacetimes. We also, introduce the finite stress energy tensor of this theory. Section IV is devoted to the thermodynamics of the black hole solutions of the theory. We calculate the finite action, the total mass, the temperature, the charge and the electric potential. We calculate the entropy through the use of Gibbs-Duhem relation and Wald formula and find that they are consistent. We, also, investigate the first law of thermodynamics. In Sec. V, we investigate the thermal stability of the solutions in canonical ensemble. We finish our paper with some concluding remarks.

II. ACTION AND FIELD EQUATIONS
The bulk action of Lovelock gravity in n + 1 dimensions may be written as [4] with [x] denoting the integer part of x, α p 's (p ≥ 2) are Lovelock coefficients, is the Euler density of a 2p-dimensional manifold, δ is the general asymmetric kronecker delta and The action (3) is written in such a way that the maximally symmetric AdS spacetime (2) is the vacuum solution of action (3). In this notation, the independent coupling constants are L and all the Lovelock coefficients.
From a geometric point of view the Lagrangian of the action (3) in 2[n/2] + 1 and 2[n/2] + 2 dimensions is the most general Lagrangian that yields second order field equations, as in the case of Einstein-Hilbert action which is the most general Lagrangian producing second order field equations in three and four dimensions. In the rest of the paper, we work in a unit system with α 1 =α 1 = 1 and the dimensionless Lovelock coefficientsα p defined as With the definition (5), the cosmological constant for AdS spacetime is In this paper, we consider the third order Lovelock gravity in the presence of electromagnetic field. Thus, the action of matter field is where A µ is the electromagnetic potential. The first term in Lovelock Lagrangian is the Einstein-Hilbert term R, the second term is the Gauss-Bonnet Lagrangian the third term is and the cosmological constant is As in the case of Einstein-Hilbert action, the action (3) does not have a well-defined variational principle, since one encounters a total derivative that produces a surface integral involving the derivatives of δg µν normal to the boundary ∂M. These normal derivatives of δg µν can be canceled by the variation of the surface action [17,26] In Eq. (8) γ ab and K ab = −γ µ a ∇ µ n b are the induced metric and extrinsic curvature of the boundary ∂M, respectively. The explicit form of the first three terms of Eq. (8) are [17]: where J and P are the traces of respectively. In Eq. (9)Ĝ (1) ab is the Einstein tensor,R abcd (γ) is the intrinsic curvature andĜ (2) ab is the Gauss-Bonnet tensor of the metric γ ab given aŝ

III. COUNTERTERM METHOD FOR STATIC SOLUTIONS OF THIRD ORDER LOVELOCK GRAVITY
It is well known that the action I bulk +I sur is not finite for asymptotically AdS solutions. Inspired by AdS/CFT correspondence, one needs to add counterterms to the gravity action in order to get a finite action. These counterterms are made from the curvature invariants of the boundary metric with the coefficients of the higher curvature terms chosen so that power law divergences in the bulk where The coefficients A p , B p , C p and D p in Eq. (12) should depend on L. These coefficients for Einstein gravity (α p = 1,α 2 =α 3 = 0) are [10] One may note that the coefficients A p , B p , C p and D p are independent of Lovelock coefficients. We apply the counterterms (12) to various static solutions of Gauss-Bonnet and third order Lovelock gravity with different topolgy and find that these coefficients are and for Gauss-Bonnet and third order Lovelock gravity, respectively. The reason that we use exactly the counterterms of Einstein gravity is as follows. First, the boundary at r = const. for static solutions is a constant curvature hypersurface and therefore no six-derivative term will be appeared in the counterterms. In other words, all the terms of a specific order of counterterms [for example R 2 and R ab R ab in C 2 (R abR ab − nR 2 /[4(n − 1)])] for static solutions are proportional to r −4 . Also, R abcd R abcd is proportional to r −4 . Therefore, in order to remove the divergences of the action which are proportional to r −4 , any combination of R 2 , R ab R ab and R abcd R abcd can be used. So, we just use exactly the counterterms of Einstein gravity. Second, A p , B p , C p , D p , for p = 2 and 3 are proportional to those of Einstein counterterm independent of the dimensions. That is, the dimensional-dependence of these coefficients are the same as those in Einstein gravity. Third, as we will see in the next section, the entropy calculated through the use of Gibbs-Duhem relation and the mass and action calculated by our counterterms is consistent with the entropy obtained by use of Wald's formula. Fourth, the mass calculated by our counterterms satisfies the first law of thermodynamics.
While the total action I bulk + I sur + I ct is appropriate in grand-canonical ensemble where δA µ is zero at the boundary, the appropriate action in the canonical ensemble where the electric charge is fixed is [27] Thus both in canonical and grand-canonical ensemble, the variation of total action about the solutions of the field equations is So, the energy-momentum tensor can be written as: The explicit expressions of T To compute the conserved charges of the spacetime, we choose a spacelike hypersurface B in ∂M with metric σ ij , and write the boundary metric in ADM form: where the coordinates ϕ i are the angular variables parameterizing the hypersurface of constant r around the origin, and N and V i are the lapse and shift functions, respectively. When there is a Killing vector field ς on the boundary, then the quasilocal conserved quantities associated with the stress tensors of Eq. (18) can be written as where σ is the determinant of the metric σ ij , and n a is the timelike unit normal vector to the boundary B. In the context of counterterm method, the limit in which the boundary B becomes infinite (B ∞ ) is taken, and the counterterm prescription ensures that the action and conserved charges are finite. No embedding of the surface B into a reference of spacetime is required and the quantities which are computed are intrinsic to the spacetimes.

IV. THERMODYNAMICS OF ADS CHARGED BLACK HOLES
The field equation of third order Lovelock gravity for the static metric in the presence of electromagnetic field may be written as where ψ(r) ≡ [f (r) − k]L 2 /r 2 , q is the charge of the black hole and m is the integration constant which is related to the mass of the solution. The electromagnetic potential for the above metric is The mass parameter m may be written in terms of horizon radius and q by using the fact that ψ(r + ) = −kL 2 /r 2 + , and therefore (m = 0 and q = 0). The three solutions of the cubic equation (21) are where All of the above three roots could be real in the appropriate range ofα 2 andα 3 . The second and third solutions are real provided u 3 > v 2 . Here, we will consider only the first solution ψ 1 (r) which is real provided u 3 < v 2 at any r. Now, we investigate the thermodynamics of the black hole solutions. The temperature of the event horizon may be calculated through the use of analytic continuation of the metric. One obtains where η = 1 + 2α 2 kL 2 r 2 The charge of the black holes per unit volume can be calculated by integrating the flux of the electric field as The electric potential Φ, measured at infinity with respect to the horizon is defined as where χ = ∂/∂t is the null generator of the horizon. One obtains use of the counterterm method introduced in the last section. It is a matter of straightforward calculations to show that the total action is finite. Since we are interested in the stability of the solutions in canonical ensemble, we calculate the finite Euclidian action per unit volume V k,n−1 in this ensemble. One obtains where β is the Euclidean time period (the inverse of temperature), η is given in Eq. (27), ξ is and I 0 is One should note that I 0 appears only in odd dimensions (even n).
For our static solution, there is a Killing vector field ς = ∂/∂t on the boundary, and therefore the quasilocal conserved quantity associated with the stress tensors of Eq. (18) is the mass of the black hole. Using the counterterm method introduced in the last section, it is a matter of calculation to obtain the mass of black hole as where M 0 = I 0 /β is the Casimir energy for the vacuum AdS metric per volume V k,n−1 , which is nonzero only in odd dimensions. In order to have positive energy M − M 0 , one should restrict the range of Lovelock coefficients in terms of Q, L and r + as That is, the third order Lovelock coefficient has a lower limit in terms ofα 2 , q, L and r + .
Using the mass and action per unit volume V k,n−1 calculated through the use of counterterm method, one may calculate the entropy per unit volume through the use of Gibbs-Duhem relation by the use of Gibbs-Duhem relation is consistent with the calculation through the use of Wald formula [25,28] where the integration is done on the (n − 1)-dimensional spacelike hypersurface of Killing horizon, g µν is the induced metric on it,g is the determinant ofg µν andL 2 is the 2nd order Lovelock Lagrangian ofg µν . Also, one may note that the thermodynamic quantities calculated in this section satisfy the first law of thermodynamics

V. STABILITY IN THE CANONICAL ENSEMBLE
Now, we study the thermal stability of the black hole solutions of third order Lovelock gravity in canonical ensemble. First, we investigate the conditions of having black hole solution in Lovelock gravity. The solution given by Eq. (23) presents a black hole solution provided f (r) has at least one real positive root. This occurs if r + ≥ r ext , or q ≤ q ext where r ext and q ext satisfy the following equation: Also, the temperature of a physical black hole should be positive. One may note that the temperature changes its sign at the root of η = 0 (r crit ). The radius r crit depends on L and Lovelock coefficients, while r ext depends on q too. In Fig. 1, the vertical line is η = 0 line. This figure shows that one may have only small (η < 0) and large black holes, and there is no medium black hole solution (r crit < r + < r ext ). This feature does not occur for black holes of Einstein gravity or Lovelock gravity with positive Lovelock coefficients, since η = 0. Thus, for the case of Lovelock gravity with negative Lovelock coefficient(s) one may divide the black holes into two classes with negative and positive η.
In the canonical ensemble Q is fixed and therefore the black hole solutions are stable provided (∂ 2 M/∂S 2 ) Q > 0 in the range that T is positive. Using the expressions for mass and entropy given in Eq. (32) and (34), one may calculate (∂ 2 M/∂S 2 ) Q : where  Figure 3 shows that there is no large black holes for large negativeα 3 , while the large black holes for small negativeα 3 are not stable.  Figure 4 shows that there is no small black hole solution for positiveα 3 , while the small black hole with negativeα 3 are stable. Figure 5 shows that there is no negative-η (small) uncharged black hole with event horizon, while the large black holes (η > 0) show a phase transition between large and very large ones. One may see in

VI. CLOSING REMARKS
The concepts of action plays a central role in gravitation theories, but the sum of the bulk and the surface terms diverges. In this paper, we introduced the counterterms that remove the non-logarithmic divergences of static solutions of third order Lovelock gravity. We did this by defining the cosmological constant in such a way that the AdS metric (2) is the vacuum solution of Lovelock gravity. Indeed, the cosmological constant (4) makes the asymptotic form of the solutions of Lovelock gravity to be exactly the same as that of Einstein gravity. Thus, we employed the counterterms introduced for Einstein gravity in [10] and found that the power law divergences of static solutions in the action of Lovelock gravity can be removed by suitable choice of coefficients.
We found that the counterterms are independent of Lovelock coefficients and the dimensionally dependent of them is the same as those of Einstein gravity. The main difference of our work with that of Ref. [16] is that the counterterms are exactly the same as those of Einstein gravity and do not depend on Lovelock coefficients. The only job which remains is that one needs to calculate the coefficients A p , B p , C p and D p . For example if one wants to remove the divergences of the action due to the Gauss-Bonnet term, one should calculate A 2 , B 2 , C 2 and D 2 that remove the divergences of √ γL 2 , without regarding the other Lovelock terms. This enables one to generalize this method to other higher curvature theories of gravity easily, including fourth order Lovelock gravity or f (R) gravity. We also introduced the finite energy-momentum tensor in third order Lovelock gravity.
In addition, we employed these counterterms to calculate the finite action and mass of the static black hole solutions of third order Lovelock gravity. Calculating the temperature, the electric charge and electric potential and using the calculated finite action and mass, we showed that the entropy calculated through the use of Gibbs-Duhem relation is consistent with the calculated entropy by Wald's formula. We, also, showed that the conserved and thermodynamic quantities satisfy the first law of thermodynamics. Finally, we investigated the stability of charged black holes of Lovelock gravity in canonical ensemble. We found that the black holes with respect to the sign of η given in Eq. (27)  The first term of Eq. (18) which is the energy-momentum tensor of the surface term is T (sur) ab = 1 8π {K ab − Kγ ab + 2α 2 [3J ab − Jγ ab − 2Ĝ (a c K b)c + 2R ab K − K abR + 2K cdĜ cd γ ab − 2K cdR acbd ] +3α 3 [5P ab − P γ ab + 2KĜ (2) ab + L 2 (K ab + Kγ ab ) + 4JR ab − 24J (a cd K cd + 6J cdR cd − JR + 2KK cd K efR cedf − 4K c e K cd K f hR df eh )]}