On Thermodynamic Stability of Black Holes. Part I: Classical Stability

We revisit the classical thermodynamic stability of the standard black hole solutions by implementing the intrinsic necessary and sufficient conditions for stable global and local thermodynamic equilibrium. The criteria for such equilibria are quite generic and well-established in classical thermodynamics, but they have not been fully utilized in black hole physics. We show how weaker or incomplete conditions could lead to misleading or incorrect results for the thermodynamic stability of the system. We also stress the importance of finding all possible local heat capacities in order to fully describe the classical equilibrium picture of black holes. Finally, we thoroughly investigate the critical and phase transition curves and the limits of the classical analysis. This paper is the first in the line of intended works on thermodynamic stability of black holes in modified theories of gravity and holography.

Black hole physics is definitely one of the major topics in modern theoretical physics.This is largely due to the fact that a consistent description of their properties relies on fundamental principles from thermodynamics, statistical mechanics, quantum theory and gravitation.The interest in the area is also significantly boosted by the recent discovery of gravitational waves produced by black hole collisions [1], and further by the visual confirmation of the existence of supermassive black holes in the center of M87 and our own galaxy [2][3][4].Evidently, near-future experiments and observations are also expected to reveal surprising new results.
In this context thermodynamics of black holes [5][6][7][8][9][10][11][12][13][14] could provide an avenue for direct tests of our theories of gravity.The latter is a consequence of the generic feature of thermal systems to allow descriptions with only few observable parameters.For black holes these macro parameters reduce to the mass, entropy, charge, angular momentum and a few more variables depending on the underlying model of gravity.More importantly, the origins of black hole thermodynamic are deeply rooted in the theory of quantum gravity, which is under active investigation by various promising approaches such as string theory.On the other hand, insights on the quantum nature of gravity can only become available near the event horizon, where strong gravitational effects become relevant and thermodynamics can no longer be treated classically.These arguments make black holes one of the most challenging objects to study.
As apparent by the images of supermassive black holes the thermodynamic states of such compact objects strongly depend on their surrounding environment.For example, this could be achieved by matter accretion onto the even horizon [15][16][17][18][19], or via some radiative process such as Hawking radiation [10].In these cases, the natural question to ask is under what condition the system is in thermodynamic equilibrium with its surroundings?To our knowledge, the answer to this question has not been explored in sufficient detail even for the standard black hole solutions, although the criteria for equilibrium are quite generic and well-established in standard thermodynamics [20][21][22][23][24].There are only few known examples of using such generic criteria for black holes.For example, discussions on the necessary and sufficient conditions of thermodynamic stability for higher dimensional black holes have been presented in [25,26].A subsequent derivation of general stability criteria for black holes has been presented by A. K. Sinha [27][28][29].The latter coincide with the Sylvester criterion for positive definiteness of the Hessian of the energy potential.Applications to thermodynamic stability of quantum black holes appear in [27][28][29][30][31][32].Thermal stability of black holes with arbitrary hairs has been investigated by [31][32][33].Our work aims to supplement these studies by implementing the full classical criteria for local and global thermodynamic stability of the standard black hole solutions in General relativity.We also stress the importance of working in natural parameters for a given thermodynamic representation and raise some caution when working with non-generic or partial stability criteria.
In classical thermodynamics there are two types of equilibria -local and global1 .If a system resides in a global thermodynamic equilibrium then, by definition, it has the same temperature, the same pressure, the same chemical potentials etc, everywhere within its boundaries.In this case, one can study the global thermodynamic stability in a given representation by considering the properties of the Hessians of the corresponding thermodynamic potentials.The global thermodynamic analysis is based on two equivalent criteria: the Hessian eigenvalue method and the Sylvester criterion for positive definiteness of quadratic forms.One can use both criteria independently as sufficient conditions for global thermodynamic stability of any systems including black holes.
The system is said to be in a local thermodynamic equilibrium if it can be divided into smaller constituents, which are individually in approximate thermodynamic equilibrium.In each partial system the intensive thermodynamic state quantities assume definite constant values and do not vary too strongly from one partial system to another, i.e.only small gradients are allowed.The study of local thermodynamic stability is based on the admissible heat capacities.Specifically, for a system to be locally stable with respect to a perturbation in a set of parameters the corresponding heat capacities must be strictly positive.
It is important to note that local equilibrium does not imply a global one.On the other hand, it is natural to assume that a system in global thermodynamic equilibrium is also locally stable.This is evident by the fact that the components of the Hessians of the thermodynamic potentials can be related to the local heat capacities of the system.
The goal of this work is to present the theory of classical thermodynamic stability in details and then revisit the standard black hole solutions in the light of the necessary and sufficient condition for thermodynamic equilibrium2 .We show that all considered black hole solutions are globally unstable, but some local stability with respect to their heat capacity can be retained.The instability of asymptotically flat black holes is a persistent feature in any dimension.For example in [35] for D ≥ 5 it was shown that all asymptotically flat rotating and neutral black holes are unstable and this was extended to include the charged case in [25].
The structure of the paper is the following.In Section 2 we present the necessary and sufficient conditions for global thermodynamic stability in energy and entropy representations 3 .
In Section 3 we revisit the thermodynamic instability of the Schwarzschild black hole solution only as a didactic example.In Section 4 we study the thermodynamic stability of the Reissner-Nordström (RN) black hole.We confirm its global instability, but show that its is locally stable with respect to certain processes.In Section 5 the thermodynamic stability of Kerr solution is studied in details.We verify that Kerr is globally unstable, but it has regions of local stability for particular values of the angular momentum.We also show that for processes with constant mass the J → 0 is a regular limit to a new locally stable state, which differs from the Schwarzschild black hole.In Section 6 we verify that even the three parametric thermodynamic space of equilibrium states for the Kerr-Newman (KN) solution is not enough to support global thermodynamic stability of the system.In this case we study all admissible heat capacities and derive the various regions of local thermodynamic stability, which have not been fully investigated previously.Finally, in Section 7 we give a brief summary of our results.

Description of thermodynamic stability
In this section we present the necessary and sufficient conditions for classical global and local thermodynamic stability in energy and entropy representations.

Energy and entropy representations
The thermodynamic representation of a given system is defined by the choice of thermodynamic potential used to describe the properties of the system and the constraints it is subjected.Hence, the energy representation is used when the preferable thermodynamic potential is the internal energy E of the system.In this case, one naturally imposes constraints on the entropy and other extensive variables of the system.Consequent application of Legendre transformation along one or several natural parameters of the internal energy leads to other energy derived thermodynamic representations, called free energies, which fully describe the properties of the system on their own.There exist other representations, which cannot be derived from the energy potential via Legendre transformation.Such representation is the entropy representation, where the entropy derived potentials are called Massieu-Planck potentials or free entropies.At the end the choice of a representation depends on the initial constraints imposed on the system.In energy representation one defines the set of all extensive4 ⃗ E = (E 1 , E 2 , ..., E n ) T and all intensive ⃗ I = (I 1 , I 2 , ..., I n ) T parameters, which describe the possible macro states of the system.In these terms the first law of thermodynamics in equilibrium is written by 5 : where E is the (internal) energy of the system.The form of the first law is specifically chosen to represent I a as generalized thermodynamic forces and E a as generalized thermodynamic coordinates by analogy of classical mechanics.When it is possible to express the energy potential E as a function of its natural extensive variables ⃗ E, one finds the so called fundamental relation: This is a particularly important relation due to the fact that it can be used to directly extract the relevant thermodynamic properties of the system via the equations of state: Here the parameters in the subscript are kept fixed except for E a .The set of equations (2.1)-(2.3)define the mathematical form of the energy representation for a system in equilibrium.However, it is not always practical to work with the energy, due to the fact that different constraints on the system may require different control parameters.In this case, it may be useful to transform to another representation without loosing any relevant thermodynamic information of the system.This can be archived by the well-known Legendre transformation.Performing Legendre transformation L along one or several natural parameters of the energy we can obtain all of the standard free energy potentials.In this case, the one-parameter Legendre family of energy derived potentials Φ a , is given by 6 : ) ) . . . (2.6) Consequently, the two-parameter Legendre family of energy derived potentials Φ ab , a ̸ = b, is7 : . . . (2.9) This approach can be generalized to any number of extensive variables.At the end, the Legendre transformation along all the extensive variables leads to the trivial (or null) potential Φ 1,2,...,n = 0, which is due to the Euler homogeneity relation I a E a . (2.10) However, the potential Φ 1,2,...,n may not be trivial if the energy of the system is a quasihomogeneous function of degree r and type (r 1 , ..., r n ), i.e.
where under dilatations by a scale factor τ > 0 one has the generalized Euler relation: (2.12) In black hole thermodynamics this is related to the so called Smarr relation8 .
The situation is similar if one works in the entropy representation.In this case, choosing the entropy S as a thermodynamic potential depending on its natural extensive parameters S = S(S 1 , S 2 , ..., S n ), one can write the first law of thermodynamics in the form where the intensive variables ⃗ λ = (λ 1 , λ 2 , ..., λ n ) T are the thermodynamically conjugate parameters of ⃗ S. The equations of state follow naturally by 9   λ a = ∂S( ⃗ S) ∂S a S 1 ,..., Ŝa ,...,S n . (2.14) Legendre transformation of the entropy potential along one or several of its natural parameters is used to obtain the entropy derived family of potentials.The latter are known by several names: Massieu-Planck potentials, free entropies or free information potentials.It is important to note that entropy is not a Legendre transformation of the energy and thus entropy representation and its derived potentials are generally different from the energy representation related potentials.In fact, different potentials correspond to different constraints to which the system may be subjected 10 .The thermodynamic properties of the system can be fully described once the fundamental relation in the chosen representation has been established.
Let us clarify this point with a simple example.Assume that we want to study Kerr black hole with first law in energy representation given by dM = T dS + ΩdJ. (2.15) In this case, the natural parameters of the mass are the entropy S and the angular momentum J. Hence the equilibrium manifold is defined by the embedding of the fundamental relation M = M (S, J), which is a two-dimensional surface in R 3 .This representation is useful to study the stability and the critical phenomena of the system with respect to S and J.This means that it is not a good idea to look for the critical properties of the temperature in this representation, because T is not a natural variable of the mass 11 .In such cases one looks for an appropriate thermodynamic potential, whose natural variable is T .For example, one can Legendre transform the mass to the Helmholtz free energy (canonical ensemble), It is now evident that the natural space is F = F (T, J) and one can use F to study the properties of the system in terms of the temperature.Similarly, one can refer to the Gibbs free energy, etc.
We are now ready to describe the generic conditions for thermodynamic stability.

Classical criteria for global thermodynamic stability
We say that a thermodynamic system is in equilibrium with its surroundings if the state quantities do not spontaneously change over considerably long period of time.According to the , where b = 2, 3, ..., n. 10 Thermodynamic potentials are naturally used to describe the ability of a system to perform some kind of work under given constraints.These constraints are usually the constancy of some state variables like pressure, volume, temperature, entropy, etc.Under such conditions the decrease in thermodynamic potential from one state to another is equal to the amount of work that is produced when a reversible process carries out the transition, and hence is the upper bound to the amount of work produced by any other process, [38].
11 It is well known that one may loose information of the system if not working in natural variables.
laws of thermodynamics [20][21][22][23][24] the necessary, but not sufficient, conditions for thermodynamic equilibrium between the system and its surroundings can be established by the equalities of the corresponding intensive parameters, I a = I * a , of the system I a and the reservoir I * a .These parameters may include temperature, pressure, chemical potentials etc.The conditions can easily be derived by the restriction on the first variation of the internal energy of the system during a virtual process: (2.17) The space of possible states of equilibrium (compatible with constraints and initial conditions) is called the space of virtual states.Due to the first law in equilibrium one has (2.3), thus the necessary conditions for equilibrium become One can reach to the same conclusion in the entropy representation by δ (1) S( ⃗ S) = 0. On the other hand, the sufficient conditions for global thermodynamic equilibrium, and thus global thermodynamic stability, can be derived by the sign of the second variation of the energy or the entropy consistent with the second law of thermodynamics.Considering the energy as a potential the second variation should be strictly positive due to the fact that in equilibrium the energy of the system assumes its minimum.Here H E is the symmetric n × n Hessian matrix of the energy given by (2.20) The inequality δ (2) E > 0 defines H E as a positive definite quadratic form.This means that for global equilibrium it is sufficient that all eigenvalues12 ε a > 0, a = 1, ..., n, of the Hessian of the energy be strictly positive.
In the entropy representation the second variation should be strictly negative due to the fact that in equilibrium the entropy of the system settles at its maximum.The inequality δ (2) S < 0 defines δ (2) S as a negative definite quadratic form.
For establishing a global equilibrium it is sufficient that all eigenvalues s a < 0, a = 1, ..., n, of the Hessian of the entropy be strictly negative.An alternative set of sufficient conditions for global13 thermodynamic stability is given by the Sylvester criterion for positive/negative definiteness of the Hessians.In energy representation the energy defines a global convex function, thus the Hessian of the energy is positive definite quadratic form.Therefore, all of the principal minors ∆ k > 0 of the Hessian of the energy must be strictly positive.In entropy representation this criterion has alternating signs (−1) k ∆ k > 0 due to the fact that entropy is globally concave function 14 .An example is shown in Appendix A.

Heat capacities and local thermodynamic stability
One of the major effects of heat transfer is temperature change defined by dQ = CdT = T dS, where the extensive quantity C is called the total heat capacity of the system.One can also write and since dQ depends on the nature of the process 15 , so does C. Hence, for different processes one has different heat capacities.The general definition of a heat capacity C x 1 ,x 2 ,...,x n−1 , at fixed set of thermodynamic parameters (x 1 , x 2 , ..., x n−1 ), is given by the derivative of the entropy in a certain space of variables (y 1 , y 2 , ..., y n ), namely 16 [39]: The Nambu brackets { }, used in the formula above, generalize the Poisson brackets for three or more independent variables (see Appendix B).Furthermore, the set of constant parameters x 1 , x 2 , ..., x n−1 could be a mix of all kinds of intensive and extensive variables.Additionally, all the relevant state quantities become functions of the independent parameters y 1 , y 2 , ..., y n .In this case we say that (y 1 , ..., y n ) define the coordinates of our space of equilibrium states.
Local thermodynamic equilibrium can be defined by quasi equilibrium between different parts of the system, where sufficiently small gradients of the parameters are still allowed.The identification of local stability with the positivity of certain heat capacity is related to the components of the Hessian, where imposing the generic conditions for stability always require C > 0. This is most evident for simple systems (see for example [20][21][22][23][24]).Therefore, one can insist that the classical condition for local thermodynamic stability, with respect to some fixed parameters (x 1 , x 2 , ..., x n−1 ), is 17C x 1 ,x 2 ,...,x n−1 > 0. (2.24) Heat capacities are also important for identifying critical and phase structures in the system.Specifically, if a given heat capacity diverges 18 or changes sign this would signal the presence of a phase transition and the breakdown of the equilibrium thermodynamic description.
In the forthcoming sections we are going to revisit the thermodynamic stability of the standard black hole systems from general relativity in the light of the strict classical criteria presented above.

Thermodynamic instability of the Schwarzschild solution
The simplest space-time solution of general theory of relativity is the Schwarzschild solution which describes a static spherically symmetric black hole.It is purely didactic to study its thermodynamic properties, which are defined only by three parameters: the mass M , the entropy S and the Hawking temperature 19 T on the event horizon of the black hole.In the energy representation the mass of the black hole is identified by the energy of the system, hence in equilibrium the first law of thermodynamics is simply written by20 dM = T dS.
In terms of the energy natural parameter S one has Assuming M, S, T > 0 the Schwarzschild thermodynamics can be presented in a more convenient form where we have introduced the following notations: The global thermodynamic instability of the Schwarzschild black hole solution follows directly from the Hessian of the mass, which has only one element It is evident that the sign contradicts the general thermodynamic stability criteria (A.4)-(A.8).
In order to analyze the local thermodynamic stability one looks at the heat capacity of the system Its negative sign defines the local thermodynamic instability of the Schwarzschild black hole, so it can radiate.One can think of this as a process of evaporation where the temperature T rises on the expense of the decreasing mass M of the black hole 21 .At the end the system will evaporate explosively unless quantum effects are taken into account.The Schwarzschild black hole can be stabilized by placing it in a cavity with a heat bath at a finite distance from the horizon [41] or in the presence of negative cosmological constant [42] of sufficient magnitude.Quantum corrections also lead to thermodynamically stable Schwarzschild configurations [43][44][45][46].Next we revisit the classical thermodynamic stability of the Reissner-Nordström and the Kerr black hole solutions.

Thermodynamics in Reissner-Nordström spacetime
Reissner-Nordström (RN) solution is the charged generalization of the Schwarzschild black hole 22 : Here M > 0 is the mass and Q ∈ R is the charge of the black hole.The event horizon is located at r + = M + M 2 − Q 2 and its existence (r + > 0) assumes the condition M 2 > Q 2 .Note that the extremal case M = |Q| leads to T = 0, which is in contradiction to the third law of thermodynamics 23 .In equilibrium the first law yields where S is the entropy, T is the Hawking temperature, and Φ is the electric potential of the black hole.The relation (4.2) defines the thermodynamics of the RN black hole in (S, Q) space: After introducing the set of new parameters, one can write the first law in the form dm = τ ds + ϕdq, where In this representation the existence condition The latter is satisfied above the blue parabola in (s, q) space (Fig. 1).Note that the curve s = q 2 is forbidden by the third law of thermodynamics.The condition for existence (4.8) together with (4.5) and (4.7) also lead to (4.9) Figure 1: The existence of the black hole is defined above the blue parabola s = q 2 .The orange parabola s = 3q 2 represents the Davies curve for the heat capacity C q .One has local thermodynamic stability with respect to fixed charge q = const in the region between the blue and the orange curves.
For fixed mass m = const RN is locally stable everywhere above the blue curve.No locally stable regions exist for RN with respect to constant electric potential ϕ.

Global thermodynamic instability of Reissner-Nordström black hole
The global thermodynamic stability of the RN black hole solution can be determined by the properties of the Hessian of the mass in (s, q) space: According to the general theory the sufficient conditions for having a stable global equilibrium in the mass-energy representation insists on λ 1 > 0 and λ 2 > 0. The two eigenvalues of H are A closer inspection of their signs shows that λ 1 > 0 and λ 2 < 0 for all admissible values of s and q in the region of existence (4.8) of the black hole.This indicates that the Reissner-Nordström black hole cannot be globally stable from thermodynamic standpoint.An additional check of the Sylvester criterion for the positive definiteness of H also confirms our this.The conditions for global thermodynamic stability require: In this case, the first condition H ss > 0 is satisfied below the orange parabola s < 3q 2 (Fig. 1).
The second condition H qq > 0 is always true.The final condition can not be true in the region of existence (4.8), which is evident from the determinant of the Hessian: Hence, global thermodynamic stability criteria (A.4)-(A.8)cannot be simultaneously satisfied and the RN solution is globally unstable from thermodynamic point of view.
A note of caution is advised here when using the Sylvester criterion.It would be misleading to consider only the weak convexity conditions of the Hessian of the mass along s and q: which are satisfied in the region q 2 < s < 3q 2 .This would falsely indicate that RN is globally stable for q 2 < s < 3q 2 .Furthermore, it is not recommended to weaken the strict positive definiteness of the Hessian by positive semi-definiteness.The latter may incorrectly indicate that the system is stable on some of the critical or phase transition curves.

Local thermodynamic stability of Reissner-Nordström black hole
The local thermodynamic stability of the RN black hole is determined by the admissible heat capacities of the solution in (s, q) space.By definition, for a fixed parameter x, one has where the Nambu brackets are given by simple Poison brackets (Jacobians): Therefore, the relevant heat capacities of the RN black hole in (s, q) space are 24 : C m = τ ∂s ∂τ m = τ {s, m} s,q {τ, m} s,q = s (s − q 2 ) q 2 , (4.17) C q = τ ∂s ∂τ q = τ {s, q} s,q {τ, q} s,q = 2s(s − q 2 ) 3q 2 − s .
For a processes with constant mass one notes that C m is always positive in the region of existence (4.8).Thus RN is locally stable above the blue parabola (Fig. 1).Using the relation (4.5) we can express C m as a function of only one variable s or q and the constant parameter m c , i.e.
For a processes with constant electric potential the heat capacity C ϕ is always negative, which leads to unstable black hole for this kind of processes.Using (4.7) we can find C ϕ as a function of the electric charge q and the constant parameter ϕ c Finally, when the charge is fixed, the region of local stability C q > 0 is q 2 < s < 3q 2 .It is located between the blue and the orange parabolas (Fig. 1).Above the orange parabola the RN black hole is unstable for processes with constant charge.

Critical curves and phase transitions of Reissner-Nordström black hole
The critical curves, also known as Davies curves, are defined by the divergences of the heat capacities (C → ±∞) or by the curves where a change of sign occurs (C = 0).In the (s, q) space the RN black hole has the following critical curves: The heat capacities C m and C q change sign on the blue parabola s = q 2 (Fig. 1), where the temperature of the black hole is zero.The latter corresponds to the extremal case, which cannot be reached for finite number of fluctuations due to the third law of thermodynamics.If we consider a process with constant mass, the line q = 0 is a Davies curve.On this line the RN black hole admits a phase transition, where the classical equilibrium description breaks down.
For processes with constant electric potential (4.7) the limit q → 0 leads to s → 0, which is forbidden by the third law of thermodynamics.
For processes with constant electric charge the orange parabola s = 3q 2 is a Davies curve.Here RN phases from locally stable (below the orange curve) to locally unstable (above the orange curve) thermodynamic state.
5 Thermodynamic stability of Kerr solution

Thermodynamics in Kerr spacetime
The Kerr solution describes a rotating uncharged axially symmetric black hole with metric25 (5.1) Here we use the standard notations: 2) The existence of the event horizon, r + = M + √ M 2 − a 2 , leads to M > |J|.We do not include the extremal case M = |J| (a = M ) due to the violation of the third law of thermodynamics.In the energy representation the first law is written by where the parameters of the solution in (S, J) space take the form: To simplify the expressions we introduce new parameters: where dm = τ ds + ωdj and ) Insisting on s, τ > 0 the existence of the Kerr black hole in in (s, j) space is determined by s > |j|, (5.9) which is satisfied above the blue line depicted on Fig. 2. The condition for existence (5.9) together with (5.6) and (5.8) also leads to (5.10)

Global thermodynamic instability of Kerr black hole
The Hessian of the mass in (s, j) space is given by (5.11) The two eigenvalues of the Hessian (5.12) differ by signs: λ 1 > 0 and λ 2 < 0, for all admissible values of s and j.This is sufficient to indicate that Kerr black hole is globally unstable solution.The same conclusion can be drawn from the Sylvester criterion, where the determinant of the Hessian is always negative, (5.13)

Local thermodynamic stability of Kerr black hole
The admissible heat capacities of the Kerr black hole in (s, j) space are (5.15) (5.16) For processes with constant mass C m is positive in the region of existence (5.9), hence Kerr is locally stable above the blue line shown in Fig. 2. Using the relation (5.6), together with the condition for existence (5.9), we find C m as a function of one variable s or q and the constant mass m c : (5.17) For processes with constant angular velocity the heat capacity C ω is always negative in the region of existence, thus the black hole is locally unstable for this kind of processes.Using (5.8) we can find C ω as a function of s and the constant parameter ω c : Finally, local thermodynamic stability for processes with fixed angular momentum requires (5. 19) This is the region between the blue and the orange lines (Fig. 2), where Kerr is locally stable with respect to C j .Above the orange line the black hole is locally unstable for such processes.The heat capacity C j (s) is a function of one variable s, because j c is a constant.

Critical curves and phase transitions of Kerr black hole
The set of critical curves in (s, j) space are defined by (5.20) The heat capacities C m , C ω and C j change sign on the blue line s = |j| (Fig. 2), which corresponds to the extremal case.For process with constant mass, on the line j = 0 (the angular velocity is ω = 0) the Kerr black hole reduces to a new locally stable state with respect to C m .The latter is not a Schwarzschild black hole, since the heat capacity C m differs from (3.7).It is a regular limit and not a phase transition.
For processes with constant angular velocity (5.8), the limit j → 0 leads to s → 0, which is forbidden by the third law.
Finally, for fixed angular momentum, the (orange) line s = |j| 3 + 2 √ 3 is a Davies curve.It indicates a phase transition from locally stable to unstable thermodynamic state of the Kerr black hole with respect to C j .
6 Thermodynamic stability of Kerr-Newman solution

Thermodynamics in Kerr-Newman spacetime
The Kerr-Newman (KN) spacetime is the charged version of the Kerr solution with line element where we have the notations: In addition to the gravitational field, the KN black hole is surrounded by a stationary electromagnetic field which is completely determined by the charge Q and the parameter a.The standard form of the KN thermodynamics is presented in Appendix C.Here we prefer to work with new set of parameters defined by The latter renders the KN thermodynamics in the form dm = τ ds + ωdj + ϕdq with 2 s 3 j 2 + (q 2 + s) 2 , (6.5) ω = ∂m ∂j s,q = j s j 2 + (q 2 + s) 2 , (6.6) The region of existence of the Kerr-Newman black hole in (s, j, q) space is now given by s > j 2 + q 4 .(6.8) This condition is satisfied above the blue surface Fig. 3.

Global thermodynamic instability of Kerr-Newman black hole
The global thermodynamic stability is determined by the components of the Hessian: where the explicit expressions are given in Appendix D. The eigenvalues of the Hessian satisfy a cumbersome cubic equation, which makes them difficult for an analytical treatment 26 .Fortunately, we can use the conditions imposed by the Sylvester criterion.In mass-energy representation global thermodynamic stability insists on positive definiteness of the Hessian of the mass.The latter suggests that the first level principal minors of the Hessian should satisfy together with the conditions for the second level principal minors: and the determinant of the Hessian itself: In this case, it only suffices to calculate the determinant of the Hessian and show that it is always negative in the region of existence (6.8): Therefore the RN black hole is globally unstable from thermodynamic standpoint.
The heat capacity with respect to fixed angular velocity and electric potential C ω,ϕ = τ ∂s ∂τ ω,ϕ = τ {s, ω, ϕ} s,j,q {τ, ω, ϕ} s,j,q = −2s q 2 + s 3 s 2 − (j 2 + q 4 ) is always negative in the region of existence.Hence the black hole is locally unstable for a processes with constant angular velocity and electric potential.The heat capacity with fixed angular velocity and charge C ω,q = τ ∂s ∂τ ω,q = τ {s, ω, q} s,j,q {τ, ω, q} s,j,q = 2s q 2 + s 2 s 2 − (j 2 + q 4 ) j 2 + (q 2 + s) 2 R(s, j, q) (6.26) is positive between the blue and the red surfaces (Fig. 3), which is the region of local stability.Above these surfaces the KN black hole is unstable thermodynamically for such processes.The region of existence of the KN black hole is above the blue surface s = j 2 + q 4 .The red surface (6.16) represents the Davies surface for C ω,q .One has local thermodynamic stability with respect to C ω,q between the two surfaces.
For a process with constant angular momentum and electric potential, the heat capacity C j,ϕ = τ ∂s ∂τ j,ϕ = τ {s, j, ϕ} s,j,q {τ, j, ϕ} s,j,q = 2s s 2 − (j 2 + q 4 ) j 2 (3q 2 + s) + (q 2 + s) 3 O(s, j, q) (6.27) is positive between the blue and the orange surfaces (Fig. 4), hence the black hole is locally stable in this region.Above the orange surfaces the black hole is unstable.Finally, we consider a process with constant angular momentum and electric charge.The corresponding heat capacity C j,q = τ ∂s ∂τ j,q = τ {s, j, q} s,j,q {τ, j, q} s,j,q = 2s s 2 − (j 2 + q 4 ) j 2 + (q 2 + s) 2 P (s, j, q) (6.28) is positive between the blue and the purple surfaces (Fig. 5), hence the KN black hole is locally stable against such fluctuations.Above the purple surface the black hole is unstable and can radiate.
Figure 5: The purple surface (6.18) represents the Davies surface for C j,q .One has local thermodynamic stability with respect to C j,q between the blue and the purple surfaces.Above the purple surface the black hole is unstable.

Conclusion
We conducted a thorough investigation of the thermodynamic stability of standard black hole solutions within the framework of general relativity, employing established techniques from classical thermodynamics.While these methods are widely applied in the study of conventional systems, their utilization in gravitational contexts, particularly in the thermodynamics of black holes, remains relatively unexplored.Our approach offers a significantly more systematic and comprehensive analysis compared to the scattered results found in existing literature.We employ two rigorous global criteria for equilibrium: the Hessian eigenvalue method and the Sylvester criterion, which evaluates the positive definiteness of the mass-energy Hessian quadratic form.We illustrate that, with the exception of the simplest case of the Schwarzschild black hole, which is consistently thermodynamically unstable, Reissner-Nordström, Kerr, and Kerr-Newman solutions exhibit local stability in certain subregions against fluctuations concerning specific fixed parameters.However, comprehensive assessments of global stability reveal that all solutions in general relativity are globally unstable from a classical perspective.This situation does not appear to improve upon considering additional thermodynamic parameters in the Kerr-Newman black hole solution.Thus, it can be inferred that Schwarzschild, Kerr, and Reissner-Nordström black holes inherit their global thermodynamic instability from the Kerr-Newman black hole.
Our methodology offers the benefit of addressing perils often overlooked, which could lead to incomplete or incorrect conclusions regarding system stability.Firstly, we show that semidefinite criteria could lead to certain contradictions with the stability of the system and should be considered with care.In essence, they allow the system to settle in the neighborhood of some of the classically forbidden regions in the state space, as shown for the RN black hole in Section 4.2.Furthermore, we highlight potential inconsistencies arising from the consideration of partial stability conditions, as demonstrated again in Section 4.2.
Moreover, we emphasize the significance of evaluating all admissible heat capacities, as they provide crucial information about phase transition points and the black hole's responses to various perturbations.However, as demonstrated by [47], the positivity of all heat capacities does not guarantee global stability of the black hole.This prompts us to suggest that solely examining heat capacities may not suffice for studying the thermodynamic stability of black holes.It is plausible that a comprehensive assessment involving all thermodynamic response functions, such as latent heats and compressibilities is necessary.We plan to investigate this matter further in a separate study.
This method has recently been employed to investigate the AdS family of black hole solutions, wherein the cosmological constant emerges naturally [48].In the aforementioned study, the authors underscore the crucial role played by the cosmological constant in stabilizing the thermodynamics of these black holes.A logical step forward from this point is to regard the cosmological constant as a thermodynamic variable, thereby introducing an effective pressure.We will explore this framework in the third installment of this series of papers.Indeed, our methodology lends itself readily to the extension of such cases, wherein the complexity of the thermodynamic state space is further heightened.
This paper marks the initial step in a series of planned investigations into the thermodynamic stability of black holes within modified theories of gravity and holography.In addition it would be intriguing to extend the full generic criteria to explore other types of black hole solutions, including, but not limited to quantum-corrected black holes, hairy black holes, black holes in lower and higher dimensions, regular black holes, black holes with extended thermodynamics, and related systems.
Furthermore, within the framework of holography, the dual quantum field theory can naturally undergo finite temperature embedding, induced by the corresponding black hole in the bulk.A prime example is the Kerr-AdS/QGP correspondence [49][50][51][52][53], where the characteristics of strongly correlated quark gluon plasma can be explored within the supergravity approximation.It would be intriguing to examine the thermodynamic stability of such systems, offering potential insights into nonperturbative quantum effects beyond the constraints of the supergravity approximation.
Another avenue of investigation involves examining the interplay between thermodynamic and dynamic stability.Investigating fluctuation theory and its connection to thermodynamic information geometry could provide valuable insights as well.Finally, one could also study the relationship between thermodynamic stability and the holographic complexity of black holes.
A Sylvester's criterion for systems with three independent parameters.Local vs. global thermodynamic stability Following Callen [21] we define the sufficient strong intrinsic global condition for thermodynamic stability as the strict convexity/concavity of the energy/entropy: which should be valid for all admissible values of the parameters and fluctuations ∆E a = E a − Ēa or ∆S a = S a − Sa .The sufficient (differential) conditions for thermodynamic stability follow from (A.1) and (A.2) by Taylor expansion up to second order in the fluctuations ∆E a (or ∆S a ) and are given by the Sylvester criterion for positive/negative definiteness of the Hessians of the energy/entropy.This criterion is in general less restrictive than the convexity/concavity conditions above, but it is sufficient to assure strict convexity/concavity of the energy/entropy in certain intervals of the independent parameters.For this reason we call Sylvester criterion "global", because it establishes strict convexity/concavity and hence global equilibrium of the system only within these intervals.
On the other hand, due to the fact that the components of the corresponding Hessians can be related to the thermodynamic coefficients such as heat capacities and compresibilities, we can define local equilibrium only by considering positiveness of these coefficients in any given ensemble.Strictly (global) convex/concave functions can be achieved in a certain parameter region if and only if all thermodynamic coefficients are stable (positive).Therefore we make a bit artificial distinction between local and global thermodynamic stability, due to the fact that one expects the Sylvester criterion to fail when one or more of the thermodynamic coefficients are negative 27 .
In energy representation the energy defines a global convex function, thus the Hessian of the energy is positive definite quadratic form.In this case Sylvester's criterion states that all the principal minors ∆ k > 0 of the Hessian of the energy must be strictly positive.In entropy representation this criterion has alternating signs (−1) k ∆ k > 0 due to the fact that entropy is globally concave function.
Let us show what this implies for n = 3 parametric thermodynamics in the energy representation.In this case the energy E is a function of its natural parameters (E 1 , E 2 , E 3 ) and its Hessian H is the following 3 × 3 symmetric matrix According to the Sylvester criterion the fist condition for global thermodynamic stability requires the first level principal minors of Ĥ be strictly positive: One should consider the previous conditions, together with the restrictions on the determinants of second level principal minors: Here, the lower index in ∆ i indicates that the i th row and column of the Hessian have been removed.The final part of the Sylvester criterion is a condition on the determinant of the Hessian itself: (A.8) Similar conditions can be stated for n = 3 in the entropy representation.In this case the first part of the Sylvester criterion yields which just reflects the fact that entropy is a concave function along its natural parameters (S 1 , S 2 , S 3 ).The second part of the criterion requires (A.10) (A.12) Finally, the third part is (A.13) Note the differences in the signs of (A.4) and (A.9), and also between (A.8) and (A.13), which are due to the convex/concave nature of the energy and the entropy.

3 8 4 4 . 1 13 5 15 6 7 20 A
Thermodynamic instability of the Schwarzschild solution Thermodynamic stability of Reissner-Nordström solution 9 Thermodynamics in Reissner-Nordström spacetime 9 4.2 Global thermodynamic instability of Reissner-Nordström black hole 11 4.3 Local thermodynamic stability of Reissner-Nordström black hole 12 4.4 Critical curves and phase transitions of Reissner-Nordström black hole Thermodynamic stability of Kerr solution 13 5.1 Thermodynamics in Kerr spacetime 13 5.2 Global thermodynamic instability of Kerr black hole 14 5.3 Local thermodynamic stability of Kerr black hole 15 5.4 Critical curves and phase transitions of Kerr black hole Thermodynamic stability of Kerr-Newman solution 16 6.1 Thermodynamics in Kerr-Newman spacetime 16 6.2Global thermodynamic instability of Kerr-Newman black hole 16 6.3Local thermodynamic stability of Kerr-Newman black hole and critical points 17 Conclusion Sylvester's criterion for systems with three independent parameters.Local vs. global thermodynamic stability 22 1 Introduction

Figure 2 :
Figure 2: The existence of the black hole is defined above the blue lines s = |j|.The orange lines s = |j| 3 + 2 √ 3 represent the Davies curves for C j .One has local thermodynamic stability with respect to fixed angular momentum j = const in the region between the blue and the orange lines.For fixed mass m = const Kerr is locally stable everywhere above the blue lines.No locally stable regions exist for Kerr with respect to constant angular velocity ω.

Figure 3 :
Figure3: The region of existence of the KN black hole is above the blue surface s = j 2 + q 4 .The red surface (6.16) represents the Davies surface for C ω,q .One has local thermodynamic stability with respect to C ω,q between the two surfaces.

Figure 4 :
Figure 4: The orange surface (6.17) indicates the Davies surface for C j,ϕ .One has local thermodynamic stability with respect to C j,ϕ between the blue and the orange surfaces.Above the orange surface the black hole is unstable.