Geometrothermodynamics for black holes and de Sitter space

A general method to extract thermodynamic quantities from solutions of the Einstein equation is developed. In 1994, Wald established that the entropy of a black hole could be identified as a Noether charge associated with a Killing vector of a global space-time (pseudo-Riemann) manifold. We reconstruct Wald’s method using geometrical language, e.g., via differential forms defined on the local space-time (Minkowski) manifold. Concurrently, the abstract thermodynamics are also reconstructed using geometrical terminology, which is parallel to general relativity. The correspondence between the thermodynamics and general relativity can be seen clearly by comparing the two expressions. This comparison requires a modification of Wald’s method. The new method is applied to Schwarzschild, Kerr, and Kerr–Newman black holes and de Sitter space. The results are consistent with previous results obtained using various independent methods. This strongly supports the validity of the area theorem for black holes.


Introduction
Some gravitational phenomena are known to have thermodynamic interpretations. It was first pointed out by Bekenstein [1] that black holes have entropy and follow the second law of thermodynamics. After the pioneering work of Bekenstein and Hawking [2], intensive discussions revealed thermodynamical aspects of black holes. Now, four laws of black hole mechanics [3] are widely accepted: The zeroth law states that B Yoshimasa Kurihara yoshimasa.kurihara@kek.jp 1 The High Energy Accelerator Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan the thermodynamic quantities, such as the temperature and surface gravity, remain constant over the horizon of a black hole [4,5]. The first law can be expressed as [3] d M H = (κc/8π G) d A H + H d J + d Q, (1) where M is the mass of a black hole, κ is the surface gravity, A H is the area of the event horizon, G is Newton's gravitational constant, H is the surface angular velocity, J is the angular momentum, is the electromagnetic potential, and Q is a charge. For (1) to be the first law of thermodynamics for a black hole, the relationship between black hole quantities, κ, A H , and the thermodynamic quantities, temperature and entropy, must be clarified. The temperature of a black hole can be equated with the Hawking temperature [2], T H =hκ/2π c, which is obtained using a semi-classical treatment of the particle radiation from the surface of a black hole. The entropy can be obtained from the area theorem [2] for black holes: where S H is the entropy of a black hole and l p = Gh/c 3 is the Planck length. Using these relations, the first term of the r.h.s. of (1) can be rewritten as κc/(8π G)d A H = T H d S H , which allows a thermodynamic interpretation of (1). The area theorem is obtained in multiple ways for various types of black holes [6][7][8][9]. Wald proposed equating the Noether charge with the entropy of black holes [7]. Wald's method is applicable to a wide variety of Einsteinian and post-Einsteinian gravitational theories [10]. The black hole entropy has been discussed from both classical and quantum approaches [11][12][13][14][15][16][17][18][19][20][21][22][23][24][25][26], and also string theoretically [27][28][29][30][31]. A nice overall review article can also be found [32]. Gour and Mayo demonstrated that the entropy of a black hole is (nearly) linear in the area of the horizon [33]. The relation between black hole entropy and the Noether charge is discussed by many authors [34][35][36][37][38][39][40][41]. In addition to those studies, the importance of boundary terms plays the essential role of the black hole entropy, that is first pointed out by Yoke [42]. Especially, the necessity of null boundary is discussed by Parattu [43,44]. Interesting discussions related the boundary terms can be found in Refs. [45][46][47][48][49][50]. In spite of these intensive studies, the area theorem for a black hole entropy has not yet been proven generally.
We propose a new method to extract entropy from the general relativistic theory. Our method is sufficiently general to discuss the thermodynamics of a wide variety of solutions of the Einstein equation with an event horizon. The Lagrangian and Hamiltonian formalisms are given in a completely covariant way in Sect. 2. In Sect. 3, the abstract thermodynamic system is rewritten using the terminology for a geometry that is parallel to the treatment for general relativity given in Sect. 2. The new method is applied to black holes and the de Sitter space in Sect. 4, and then the area theorem is examined for these examples. Finally, a summary is presented in Sect. 5. Physical constants are set to unity, e.g., 4π G =h = c = 1, in this work unless otherwise noted.
In our formalism [58], a spin form w ab and a surface form S ab are treated as fundamental phase-space variables. The spin-and surface-forms are respectively defined in a local Lorentz manifold M as w a 1 a 2 = ω a 1 μ b η ba 2 dx μ and S a 1 a 2 = where ω a 1 μ a 2 is a spin connection and E a μ is a vierbein. We use a local Lorentz metric η ab = diag(1, −1, −1, −1) and a Levi Civita tensor (symbol) 0123 = 1 on the local Lorentz manifold M. A torsion two-form T a can be defined as T a = d w e a using a local SO(1, 3) covariant derivative d w . In this work, the Fraktur letters are used for differential forms, and Greek-and Roman-indices, both run from zero to three, are used for a coordinate on the global and local space time manifold, respectively. A curvature two-form and a volume from are introduced as R a 1 a 2 = dw a 1 a 2 + w a 1 b ∧ w ba 2 and v = a 1 a 2 a 3 a 4 e a 1 ∧ e a 2 ∧ e a 3 ∧ e a 4 /4!, respectively. The Einstein-Hilbert gravitational Lagrangian is expressed using these forms as; where is the cosmological constant. The Einstein equation and torsionless condition can be obtained as a Euler-Lagrange equation by taking a variation with respect to the surface and spin forms, respectively. Next, we rewrite the classical gravitation theory in the Hamiltonian formalism in a completely covariant way. The gravitational Lagrangian form (2) can be represented as where a 1 a 2 a 3 a 4 = − a 1 a 2 a 3 a 4 is used. Here, a new operator, called the -conjugate operator, is introduced for shorter representations including Levi-Civita tensors. This is a map that maps a rankp tensor field to a rank-(n − p) tensor field defined on an n-dimensional space-time manifold such that where a is an arbitrary form of a rankp tensor and is a completely anti-symmetric tensor in n-dimensions. Note that a = a when a is anti-symmetric with respect to all indices. The surface form can be expressed using the -conjugate operator as S ab = e a ∧ e b . Therefore, the -conjugate of the surface form, which appears frequently in later parts of this report, is S ab = S ab = e a ∧ e b . In the four-dimensional manifold, the identities are satisfied for any rank-1 tensor a a 1 , rank-2 tensors b b 1 b 2 and B b 1 b 2 , and rank-3 tensor c c 1 c 2 c 3 . Note that the -conjugate is not the same as the Hodge-star operator, whose definition includes the metric tensor. One has to avoid to use the Hodge-star operator to construct a gravitational theory because the metric tensor can be obtained only after solving the Einstein equation.
For the Lagrangian density of the gravity, it is impossible to specify which terms are kinetic and potential energies. This is reasonable because the gravitational energy cannot be defined locally in general relativity. Moreover, even on a local Minkowski manifold, the discrimination between kinetic and potential energies is frame dependent. For example, no gravitational (potential) energy exists in the local inertial frame. Start from the Lagrangian density form (3): where Using this two-form, the Einstein equation can be expressed as Under these conditions, a covariant Hamiltonian formalism is introduced as follows. First, w is identified as the first canonical variable. Then, the second canonical variable M (the canonical momentum) is introduced as The reason why the spin form is taken as the first canonical variable will be discussed in Sect. 4. Therefore, the two fundamental forms {w, M = S} form a generalized phase space. The Hamiltonian density can be obtained as Note that this Hamiltonian density does not contain any information concerning the dynamics (the derivatives of the fields). The canonical equations can then be obtained immediately: The first equation of (8) gives the equation of motion: which simply leads to the Einstein equation, (5), and the second equation of (8), leads to the torsionless condition, as expected.

Geometrothermodynamic formalism
In this section, before discussing the thermodynamic aspects of the gravitation, the thermodynamic system is rewritten using geometrical terminologies. The formalization is based on that of Ref. [59]; however, it is rewritten according to the conventions of this work. Here, we confine our attention to the formal relationships between the thermodynamic variables and ignore their physical meanings. A characteristic function of the thermodynamics can be introduced as Connotation of each variables is that S is the entropy, T is the temperature, p is the pressure, V is the volume, N is the number of molecules, and μ is the chemical potential. These variables can be categorized into two types: extensive and intensive variables. According to the natural connotation given above, ξ a = (S, V, N ) are categorized as extensive variables and ζ a = (T, p, μ) are intensive variables, where a = 0, 1, 2. We assume the existence of a smooth three-dimensional base-manifold M TD and ξ and ζ are vectors belonging to the tangent and cotangent bundles on the manifold, respectively. Even though the thermodynamics can be constructed starting from either ξ a or ζ a , here an extensive set is considered to be the first canonical variables (the coordinate space). Therefore, the intensive variables are considered to be functions of the extensive variables. According to Ref. [59], S is considered to be the order parameter of the system and T = ∂ f TD /∂ S is considered as the "Hamiltonian" (with a narrower meaning compared with the discussions in previous sections that induce a "time" evolution in the system). The intensive variables are not independent of each other; however, T is assumed to be independent of S. Therefore, we can write the intensive variables as The characteristic one-form is introduced as where the metric tensor is chosen to be η ab = diag(−1, 1, 1) and the one-form is defined as x a = dξ a . The characteristic one-form can be thought of as an internal energy w TD = −dU . This metric tensor suggests that the system has a SO(1, 2) symmetry and is consistent with an interpretation that the variable S plays the role of "time" (i.e., an order parameter) and "T" plays the role of an energy (i.e., the Hamiltonian). The characteristic two-form (the Lagrangian-density form) can also be written as where z a = dζ a . The Lagrangian form is a closed form because dL TD = 0 from the definition. Therefore, a base manifold M TD is a symplectic manifold with a symplectic potential of O TD = w TD . The intensive variables are recognized as the canonical momentum: The action integral can be expressed as where the integration region (2) is an appropriate two-dimensional manifold embedded in the base manifold M and ∂ (1) is its boundary. The equations of motion that represent the Maxwell relations of thermodynamics can be obtained in an explicitly coordinate-dependent manner. The variational operation of the characteristic one-form with respect to ξ 0 = S is Therefore, the equations of motion are where the surface integration vanishes without applying any boundary condition because ∂(∂ (1) ) = 0 (a boundary of a boundary does not exist). Other sets of Maxwell relations can be obtained using variations with respect to V or N . The above equations represent the canonical equation of motion with the following interpretations: S ↔ t (time) and T ↔ H (Hamiltonian). These equations can also be written as the Euler-Lagrange equation of motion in a coordinate-independent manner as δ x I TD = 0 ⇒ z a = 0, which is, of course, equivalent to the canonical equations. The surface term O δ = ζ a δx a can be considered as another type of the symplectic potential.
Next, consider an infinitesimal translation of the characteristic one-form along a vector field υ a on the base manifold M TD such that where £ υ is the Lie derivative along the vector υ a and O υ = ζ a L υ x a can be considered as yet another type of the symplectic potential. Here dx a = dz a = ι υ ξ a = ι υ ζ a = 0 are used. The Noether current according to this translation symmetry can be obtained using O υ such that It is confirmed that the Noether current is conserved via the explicit calculation of due to (20), dL TD = 0, and the equation of motion z a = 0. Therefore, the Noether current can be represented using a closed form as J υ = dQ TD + C, where Q TD is the Noether charge and C is a term that is zero for the solution of the equation of motion.
We can see that where the second term is zero for the solution of the equation of motion. Therefore, the Noether charge can be extracted as follows: The symplectic potential and the Noether current/charge are summarized in Table 1.
For the case of the adiabatic free-expansion of an isolated gas, in which the temperature remains constant, the meaning given above is the actual physical case. The vector υ a = (1, 0, 0) is the "Killing vector" in this context because υ a z a = −dT = 0. Therefore, the Noether charge can be singled out as In this case, the Noether charge is not simply the entropy; instead, it is the temperature times the entropy.

The symplectic formalism of general relativity
A general formalization of the thermodynamics of space-time is formalized in this section. The method proposed in Wald [7] and Iyer and Wald [8,9] are carefully rewritten using our terminology.

General formalism
Let us start from the Lagrangian density of (4). This Lagrangian form has two candidates for the first canonical variables. According to the geometrothermodynamic formulation introduced in the previous section, we would like to start from the intensive variables. The two-dimensional surface form S ab can be easily recognized as an extensive variable. Conversely, the spin form w ab may be intensive because it is defined locally, independent of the total volume under consideration. Therefore, we identify the spin form as the first canonical variable, as discussed in Sect. 2. Taking a variation of the Lagrangian density with respect to the first canonical variable w ab , one gets where E a 1 a 2 = 2T a 1 ∧ e a 2 and O δ = δw a 1 a 2 ∧ S a 1 a 2 . The second term, dO δ , is the surface term. When a surface term is set to zero due to the boundary condition, Hamilton's principle leads to the equation of motion E ab = 0, which corresponds to the torsionless condition T a = 0. Here, we introduce a three-form which is called the "symplectic potential". Here, we assume the existence of an inverse of the spin form. Compared with the geometrothermodynamics given in the previous section, it is identified as the symplectic potential. Using this symplectic potential, the surface term can be expressed as Next, consider the Lie derivative of the Lagrangian form with respect to a vector field X defined on the global manifold. The Lie derivative along the vectors ξ μ ∈ X of the Lagrangian density becomes The Lie derivative and the external derivative are commutable with each other, £ ξ dw a 1 a 2 = d£ ξ w a 1 a 2 , by definition. After simple manipulations, we finally arrive at the expression from (29) without assuming the equations of motion: where If solutions to the equations of motion are assumed, the Lie derivative of the Lagrangian density vanishes up to the total derivative of the function O ξ . This symmetry may ensure the existence of a conserved current such as which is known as the Noether current. Here, ι ξ is a contraction between the form and the vector. Therefore, the second term of (32) is a three-form as well as the first term. The conservation of the Noether current can be confirmed via where dL G = 0 and two equations of motion are used. Next, let us introduce a generating function, H ξ , and its density form, H ξ , such that H ξ = H ξ , where the integration region is the appropriate three-dimensional manifold embedded in the global space-time manifold. H ξ can be understood as a generalization of the Hamiltonian introduced in the previous section, which induces a flow in the system along the general vector fields on the local manifold. For example, when a time-like vector field can be defined globally in the local manifold and is chosen as the "time coordinate", the generating function density form H t is just the standard Hamiltonian with its narrow meaning, which results in the "time evolution" of the system.
When the vector field X is defined globally on the entire global manifold, the generating function with respect to the vector ξ μ ∈ X can be defined in the follow manner. Take the variation of the Noether current with respect to the spin form to be a 1 a 2 ). Here, the equations of motion are used again. The three-form o is called the "symplectic current". The generating function density form is introduced using the symplectic current so that δ w H ξ = o. Therefore, the Noether current can be written as Conversely, because the Noether current is the closed form for the solution of the equations of motion as shown in (33), it can be expressed as due to the Poincaré lemma, at least locally. Here, the three-form C is eliminated using solutions of the equations of motion. The two-form Q ξ is the Noether charge of the gravitational field. One can actually calculate the Noether current (32) expressed in the shape of (35) such that Therefore, one can write the Noether charge as The integration of the equation (34) can be expressed as where ∂ (2) is the boundary of (3) . From (35), (38), and the equations of motion, the generating function can be obtained as follows: The symplectic potential and Noether current/charge in general relativity are summarized in Table 2. Based on the above geometrical formalization, one can now single the "entropy" out from the general relativity system by comparing the results in Sects. 3 and 4. In general relativity (G R), the first canonical variable is taken to be w ab , which may correspond to ξ a = (S, V, N ) in thermodynamics (T D). The symplectic structure can be discussed in parallel for T D and G R. For example, using the surface term O TD ↔ O δ , the conserved Noether current can be introduced as J υ ↔ J ξ . The T D analysis suggests that the entropy can be extracted from the Noether charge for Q TD ↔ Q ξ .
Before closing this section, we would like to emphasize the origin to ensure the conservation of the Noether current/charge. In above discussions, any specific properties of the vector fields are not used. Though the Killing vector-fields will be takes as the vector filed X in a next section, the conserved current/charge can be obtained for any vector filed. The continuous symmetry to induce the Noethe's theorem in this case is the variational operation with respect to the first canonical form w ab , which caused at most a term of total derivative as shown in (27). Some examples of the Noether charge are summarized in "Appendix". Moreover concrete representations for the Lie derivative £ ξ and contraction ι ξ are not important. Algebraic rules for differential forms play a essential role.

Application to a Schwarzschild black hole
Let us calculate the Noether charge for the Schwarzschild black hole solution, as our first example of this method. The Schwarzschild solution is with a coordinate of x a = (t, r, θ, φ), where f 2 Schw (r ) = 1 − 2M/r and M is the mass of a black hole. The vierbein form can be extracted as e a = f Schw dt, f −1 Schw dr, r dθ, r sin θ dφ .
From the above solutions and the torsionless condition, a unique solution of the spin form can be obtained [51]: where the lower half is omitted because it is obvious due to the antisymmetry of the spin form. The Schwarzschild space-time has the following Killing vectors: The second Killing vector is not independent of the third and fourth vectors because ξ 2 = ξ 3,4 with (θ = ±π/4, φ = 0). Let us take a closer look at the first Killing vector ξ a 1 , which is simply denoted as ξ a hereafter. In the region r > 2M, outside of the event horizon of a Schwarzschild black hole, the zeroth component of the vierbein form e 0 is globally a time-like vector. Therefore, t can be interpreted as the time coordinate in this region. It is expected that the thermodynamic quantities are related to this Killing vector, because, for example, the entropy can only be well defined for the thermal-equilibrium state, which is static with respect to t. A dual vector of the Killing vector is given by and therefore the normalization of the Killing vector is given as and is null on the event horizon of a black hole. A contraction of the spin form with respect to the Killing vector can be calculated as Therefore, the Noether charge can be obtained as Here, we put the factor k E = (4π G/c 4 ) back in front of the Lagrangian in the action integral and write the physical parameters G, c andh explicitly. By comparing this result with the thermodynamic result given in the previous section (26), this quantity can be described using the thermodynamic variables as where S Schw and T Schw are the entropy and temperature of a Schwarzschild black hole, respectively. A purely classical treatment of the gravity cannot resolve S Schw and T Schw unambiguously because the temperature of a black hole may appear as a quantum effect. The semi-classical analysis of a black hole temperature by Hawking [2] shows the temperature as We use this temperature to extract an expression of the entropy. We note that this relationship between the temperature and the surface gravity is common to all event horizons. Subsequently, the black hole entropy S Schw can be obtained as where A Schw = 4π(2MG/c 2 ) 2 is the area of the event horizon of a black hole. Here, the diameter of the event horizon of the hole r Schw = 2G M/c 2 is used. A Hawking temperature of T Schw =hc 3 /(8π G M) can be obtained from (46) and the surface gravity κ = c 4 /(4G M). This entropy is equivalent to the celebrated Bekenstein-Hawking entropy [1,2] of a Schwarzschild black hole.
The differences between Wald's original method [7] and the method given here is as follows: in Wald's method, the Killing vector was normalized such that the Killing field had a unit surface gravity. This normalization results in a factor of 4G M/c 4 in the Noether charge. (Wald explained in [7] that this normalization makes the Noether charge local.) In addition, Wald's definition of the entropy has a factor of 2π(c/h) 1 in front of the surface integration of the Noether charge. Therefore, Wald's entropy has an additional factor of 8π G M/hc 3 compared with our definition, which corresponds to the inverse temperature of a black hole. This is why the integration of the Noether charge simply gives the entropy in [7]. From a purely mathematical point of view, the two methods are equivalent to each other.

Application to a Kerr-Newman black hole
For the second example, we will treat the general solution of a Kerr-Newman black hole. The line element of a black hole having a mass M, an angular momentum J = α M, and an electric charge q can be expressed using the Boyer-Lindquist metric [60] as [61,62] where ρ = √ r 2 + α 2 cos 2 θ and = r 2 + α 2 − 2Mr + q 2 . The vierbein form can be read out from the line-element easily. The spin form can be obtained by solving the torsionless condition with the above vierbein form: The event horizon is given as one of the solutions of the equation = 0 = (r −r + )(r − r − ), which are given as r ± = M ± M 2 − α 2 − q 2 . The horizon corresponds to the larger solution r + . A Kerr-Newman black hole has two Killing vectors of ξ t = ∂ t , and ξ φ = ∂ φ . The Killing vector for the thermodynamic analysis is chosen to be ξ = β t ξ t + β φ ξ φ , where β t and β φ are the appropriate normalization factors. We note that a linear combination of Killing vectors is also a Killing vector. Two cases of Killing vectors are considered in this work: where H = α/(r 2 + + α 2 ) is the surface angular velocity. A physical constant c (the speed of light) is written explicitly in (49) for future convenience. In the other parts, the physical constants are still set to unity. Details of the treatment of the physical constants are given later in this section. Even though each case corresponds to different thermodynamic processes, the calculations of the Noether charge are the same. Therefore, in the following, the β's are not immediately specified.
The Noether charge with respect to this Killing vector has two non-zero surface forms of where the suffix "a" can be t or φ. Compared with the case of a Schwarzschild black hole, an additional factor of 1/2 is placed in front of (ι ξ w ab ) to avoid double counting owing to the two surface forms that contribute to the results. The other surface forms are zero at θ = π/2 and therefore do not contribute to the final results. These two terms have the following integrands: where the integrand f a (θ ; β t , β φ ) is given as; and Following the θ and φ integrations, the Noether charge is obtained as where Q t,φ (β t , β φ ) can be obtained by integrating a function f t,ϕ (θ ; β t , β φ ). The area of the event horizon of a Kerr-Newman black hole can be obtained as Now, we can discuss the thermodynamics of a Kerr-Newman black hole.

A Kerr black hole
First, a Kerr black hole without an electric charge is treated, which can be obtained by setting q → 0. For the first case of (49), the Killing vector is null at the black hole surface ("event horizon" = "Killing horizon"). When α is zero, which corresponds to non-rotating black holes, the normalization is the same as that for a Schwarzschild black hole. Therefore, in this case, the Noether charge is expected to be T KN S KN , where T KN and S KN are the temperature and the entropy of a Kerr-Newman black hole, respectively The Noether charge in this case can be obtained as where c K = αc/(r 2 K + + α 2 ) is the angular velocity of the event horizon and r K + = M + √ M 2 − α 2 is the diameter of the event horizon. The definition and meaning of M is given below. Here, physical parameters are put back into the expressions, which can be performed using the following replacements: c ξ t , and by multiplying a factor of c 4 /(4π G). In this unit, α has "length" dimensions. Therefore, one can confirm that the dimension of c K is t −1 , which is consistent with its interpretation. Because the angular velocity is constant on the event horizon (which is the zeroth law of black hole thermodynamics), this can be expressed as d(ω K )/dt = c K , where ω K is the rotation angle (which is constant) per unit time. When an integration with respect to time is performed during one cycle of a black hole's rotation, one obtains Therefore, the Noether charge can be obtained by summing two charges as Again, the temperature of the Kerr-Newman black hole, is borrowed from [63] and is set to T K = T KN (q → 0) and r + → r K + . As a result, the entropy of the Kerr black hole can be extracted as Therefore, the area theorem for the black hole entropy is valid for a Kerr black hole.
Next, let us take the second Killing vector, ξ = c∂ φ . In this case, one can get where J K = Mcα is the angular momentum of a black hole. After the t-integration, as before, the Noether charge becomes Q ξ = J K . The Noether charge for this Killing vector consists of J K and K , both of which are defined in classical physics, in contrast to the temperature, which can only be defined via the quantum effect. Because the angular momentum is correctly obtained as the Noether charge, the method proposed here is ensured to work correctly. Compared with the thermodynamics developed in Sect. 3, this may correspond to the isobaric process, which has a "Killing vector" of υ a = (0, 1, 0). The Noether charge induced by this Killing vector becomes Q TD = pV . This result supports the correspondence

A Kerr-Newman black hole
From (51), the Noether charge for a Kerr-Newman black hole can be obtained as where the physical constants are written explicitly. The charge q is defined to have a length dimension. The entropy of a Kerr-Newman black hole can be singled out using the Hawking temperature (53) as This corresponds to the area theorem for a Kerr-Newman black hole.

Application to de Sitter space
The de Sitter space is one of the cosmological solutions of the Einstein equation, which has a positive cosmological constant ( > 0) without any matter or gauge fields. The line element and vierbein forms of the static coordinate solution are similar to those of the Schwarzschild metrics, (40) and (41). The result can be obtained from a simple replacement of f 2 Schw (r ) → f 2 dS (r ) = 1 − r 2 /L 2 dS , with the coordinate x a = (t, r, θ, φ), where L dS = √ 3/ is the radius of the de Sitter horizon. The physical constants are again written explicitly. We note that the cosmological constant has a dimension of (length dimension) −2 in this convention. According to this change, the spin form can be obtained as Therefore, the Noether charge with respect to the Killing vector ξ t = ∂ t /c can be calculated to be Together with the surface temperature of the horizon, such that the entropy can be obtained as where A dS = 4π L 2 dS is the surface area of the de Sitter horizon. The area theorem is also valid in this case.

Concluding remarks
The method proposed in this work, called the "geometrothermodynamic method (GTmethod)", is sufficiently general to treat a large variety of gravitational objects that have event horizons, such as black holes and de Sitter space. The Noether charge associated with an appropriate Killing vector shows a clear relationship with thermodynamical objects, such as the temperature and entropy. One advantage of the GT-method is that it is a purely classical and thermodynamic method that does not assume any (Euclidean) ensembles. Moreover, its formalism is highly geometrical and can be described in a coordinate independent way. The same formalism can also be applied to thermodynamics, and therefore the relationships between the thermodynamics and gravitational theory are apparent. The resultant entropies obtained using the GT-method are consistent with previous results and confirm the area theorem for black holes and the de Sitter space.
To understand the GT-method more fully, let us recall a method proposed in [3] and discussed in [64]. The area of the event horizon of a Kerr-Newman black hole is again shown here as A = 4π 2M M + M 2 − (J/M) 2 − q 2 − q 2 , where the subscript "KN" is omitted from the expression. Here, this relation is considered to be a state equation with three independent variables, M, J , and q, and has a variation with respect to these variables: One can see after simple manipulations that Here and correctly correspond to the surface velocity and the electric potential, respectively. A question that needs to be asked here is what is the relationship between F and the temperature/entropy. For any black hole, the Hawking temperature and the area of the horizon have the relationship of T A = (r 2 + − α 2 − q 2 )/r + . Therefore, one can easily confirm the relation F = T /4. Therefore, the first term of the r.h.s. of (57) can be written as with S = A/4. This simple analysis also verifies the area theorem and the first law of black hole thermodynamics. The implicit assumption in this discussion is that the term Fδ A corresponds to the SδT term in thermodynamics. If we accept that the entropy of a black hole is only a function of the area and that the temperature is given by the Hawking radiation, the differential equation d S(A)/d A = 1/4, with a boundary condition of S = 0 at A = 0 uniquely gives the solution S = A/4. Conversely, the GT-method does not assume any relationship between the entropy and the area. The area of the event horizon appears in non-zero components of the surface form S ab in the Noether charge. A second term on the r.h.s. of (57) is given by the same procedure with a different Killing vector. The reason why the Noether charge is proportional to the surface form is very simple. In diffeomorphic theory for a d-dimensional manifold, the Lagrangian form and the Noether current must be the d-form and the (d − 1)-form, respectively. Therefore, the Noether charge, which is obtained by integrating the Noether current, must be a (d − 2)-form. In fourdimensional space-time, the Noether charge is a two-form, which can be expanded by the bases e a 1 ∧e a 2 = S a 1 a 2 . Therefore, the GT-method provides independent evidence that the entropy can be represented using only the area. In summary, the GT-method, which can extract thermodynamic quantities from a large variety of solutions for the Einstein equation with an event horizon, is proposed in this work. This method is mathematically consistent with Wald's method. However, the relationships between general relativity and thermodynamics are clearer because abstract thermodynamics can be reconstructed using geometrical terminologies that are parallel to general relativity.
The discussions given in this work are purely thermodynamic without any statistical (Euclidean) ensembles. Therefore, the conclusion, for example, that S KN T KN of a Kerr-Newman black hole, such as (55), must be independent of the microscopic details is expected to be correct, even if the expression of the Hawking temperature will receive some corrections from the (still-unknown) quantum gravity.