Thermodynamics, stability and Hawking-Page transition of Kerr black holes from R\'enyi statistics

Thermodynamics of rotating black holes described by the R\'enyi formula as equilibrium and zeroth law compatible entropy function is investigated. We show that similarly to the standard Boltzmann approach, isolated Kerr black holes are stable with respect to axisymmetric perturbations in the R\'enyi model. On the other hand, when the black holes are surrounded by a bath of thermal radiation, slowly rotating black holes can also be in stable equilibrium with the heat bath at a fixed temperature, in contrast to the Boltzmann description. For the question of possible phase transitions in the system, we show that a Hawking-Page transition and a first order small black hole/large black hole transition occur, analogous to the picture of rotating black holes in AdS space. These results confirm the similarity between the R\'enyi-asymptotically flat and Boltzmann-AdS approaches to black hole thermodynamics in the rotating case as well. We derive the relations between the thermodynamic parameters based on this correspondence.

where the system is divided into n different parts. This property is called extensivity. The two properties: additivity and extensivity are not equivalent. An additive quantity is extensive, but extensive quantities can be nonadditive too [12,13]. Black holes are very peculiar creatures in this respect because they cannot be described as the union of some constituent subsystems which are endowed with their own thermodynamics, and therefore black holes are nonextensive objects.
Looking from a different perspective, phenomenological thermodynamics of macroscopic objects has a well understood theory from statistical physics where the macroscopic properties of a given body (described by the thermodynamic parameters e.g. total energy, entropy, temperature, etc.) can be uniquely obtained from the microscopic description of the system. Standard statistical descriptions, on the other hand, usually assume that long-range type interactions are negligible, i.e. that the (linear) size of the system in question is much lager than the range of the relevant interaction between the elements of the system. Under these conditions the standard local notions of mass, energy and other extensive quantities are well defined, and by applying the additive (and therefore extensive) Boltzamann-Gibbs formula: S BG = − p i ln p i , for defining the system's entropy function, the classical thermodynamic description is recovered in the macroscopic limit.
In the presence of strong gravitational fields however, and in particular when black holes are considered, the assumption of negligible long-range type interactions can not be hold, and consequently the usual definition of mass and other extensive quantities is not possible locally. Nonlocality is indeed a fundamental feature of general relativity, and corresponding nonextesive thermodynamic phenomena have been known in cosmology and gravitation theory for a long time (see e.g. [14][15][16][17][18][19][20][21][22][23] and references therein). In fact, even as early as 1902, Gibbs already pointed out in his statistical mechanics book [24], that systems with divergent partition function lie outside the validity of Boltzamann-Gibbs theory. He explicitly mentions gravitation as an example (see e.g. [25] for more details). Therefore, the standard Boltzamann-Gibbs statistics may not be the best possible choice for defining the entropy function in strongly gravitating systems, and other statistical approaches, which could also take into account the long-range type property of the relevant interaction (i.e. gravitation) and the nonextensive nature of the problem, are also relevant and important to study.
The nonextensive nature of the Bekenstein-Hawking entropy of black hole event horizons has been noticed [26] very early on after the thermodynamic theory of black holes had been formulated [6][7][8][9], and the corresponding thermodynamic and stability problem has been investigated several times with various approaches (see e.g. [14,25,[27][28][29][30][31][32][33][34] and references therein). The general theory of nonadditive thermodynamics has also advanced significantly in the past few decades (see e.g. [12,35] and references therein), and it has been shown, that by relaxing the additivity requirement in the axiomatic approach to the entropy definition (given by Shannon [36] and Khinchin [37]) to the weaker composability requirement, new possible functional forms of the entropy may arise [38]. As a consequence, there exist certain parametric extensions of the Boltzamann-Gibbs statistical entropy formula, which seem to be more appropriate to describe systems with long-range type interactions. One such statistical entropy definition has been proposed by Tsallis [39] as: where p i are the probabilities of the microscopic states of the system, and q ∈ R is the so-called nonextensivity parameter. In the limit of q → 1, S T reproduces the standard Boltzamann-Gibbs result, however, in the case when q = 0, the Tsallis entropy is not additive, and the parameter can be attributed to measure the effects of non-localities in the system. The q-parameter is usually constant in different physical situations, and its explicit value is part of the problem to be solved.
The Tsallis statistics to the black hole problem has been investigated with various approaches (see e.g. [25] and references therein), however it has been a long-standing problem in nonextensive thermodynamics that nonadditive entropy composition rules (in general) can not be compatible with the most natural requirement of thermal equilibrium in the system [10]. They usually don't satisfy the zeroth law of thermodynamics, which requires the existence of a well defined, unique, empirical temperature in thermal equilibrium which is constant all over the system. For resolving these issues, Biró and Ván developed a method [10], called the "formal logarithm approach", which maps the original, nonadditive entropy composition rule of a given system to an additive one by a simple transformation. This procedure results a new, but also well defined entropy function for the system, which in turn, also satisfies both the equilibrium and the zeroth law compatibility requirements of thermodynamics. In case of the nonextensive Tsallis statistics, this new entropy turns out to be the well known Rényi formula [41,42], defined as which had been proposed earlier by the Hungarian mathematician Alfréd Rényi in 1959 [41]. Recently, motivated by the nonextensive and nonlocal nature of black hole thermodynamics, we proposed and studied an alternative approach to the black hole entropy problem [40]. In this model, in order to satisfy both the equilibrium the zeroth law compatibility, instead of the Tsallis description, we considered its formal logarithm, the Rényi statistics (1.2) to describe the thermodynamic entropy of black hole event horizons. The explicit details of this approach is presented in the next section, and by applying the Rényi model to Schwarzschild black holes [40], we found that the temperaturehorizon radius relation of the black hole has the same form as the one obtained from a black hole in anti-de Sitter space by using the original Boltzamann-Gibbs statistics. In both cases the temperature has a minimum. By using a semi-classical estimate on the horizon radius at this minimum, we obtained a Bekenstein bound [43] for the q-parameter value in the Rényi entropy of micro black holes (q ≥ 1 + 2/π 2 ), which was surprisingly close to other q-parameter fits from very distant and unrelated physical phenomena, e.g. cosmic ray spectra [44,45], and power-law distribution of quarks coalescing to hadrons in high energy accelerator experiments [46].
Besides the statistical approach, another fundamental problem of applying standard thermodynamic methods to black holes arising from the question of stability. In ordinary thermodynamics of extensive systems, the local thermodynamic stability (defined as the Hessian of the entropy has no positive eigenvalues) is linked to the dynamical stability of the system. This stability criteria, however, strongly relies on the additivity of the entropy function, which is a property that clearly does not hold for black holes. The simplest example of this discrepancy is the Schwarzschild black hole which is known to be perturbatively stable but has a negative specific heat (positive Hessian). Black hole phase transitions are also strongly related to the stability properties of the system (in particular the Hawking-Page transition), and since the standard methods are not reliable in nonextensive thermodynamics, one has to be very careful when considering stability and phase transitions in strongly gravitating systems.
Avoiding the complications arising from the Hessian approach to the stability problem of black holes, an alternative technique was proposed in a series of paper by Kaburaki et al. [47][48][49]. In these works the so-called "Poincaré turning point method of stability" [50] has been applied to the problem, which is a topological approach and does not depend on the additivity of the entropy function. More recently this method has also been used to study critical phenomena of higher dimensional black holes and black rings [51] and to determine the conditions of stability for equilibrium configurations of charged black holes surrounded by quintessence [52]. In section 4 we present an overview of this method.
By considering the Rényi model in the black hole problem, we also investigated the thermodynamic stability question of Schwarzschild black holes [53]. First we considered the question of pure, isolated black holes in the microcanonical approach, and showed that these configurations are stable against spherically symmetric perturbations, just like in the Boltzmann picture. However, in considering the case when the black holes are surrounded by a heat bath in the canonical treatment, we found that -in contrast to the Boltzmann approach -Schwarzschild black holes can be in stable equilibrium with thermal radiation at a fixed temperature. This results a stability change at a certain value of the mass-energy parameter of the black hole which belongs to the minimum temperature solution. Black holes with smaller masses are unstable in this model, however larger black holes become stable. These findings are essentially identical to the ones obtained by Hawking and Page in AdS space within the standard Boltzmann entropy description [1]. According to this similarity, we also analyzed the question of a possible phase transition in the canonical picture and found that a Hawking-Page black hole phase transition occurs in a very similar fashion as in AdS space in the Boltzamann statistics. We showed that the corresponding critical temperature depends only on the q-parameter of the Rényi formula, just like it depends only on the curvature parameter in AdS space. For the stability analysis we considered both the Poincaré and the Hessian methods. The latter one could also be applied since the Rényi entropy is additive for composition (see the next section), and therefore the standard stability analysis is also reliable in this case. Both approaches confirmed the same stability results.
These findings might have some relevant consequences in black hole physics. In particular, if an effective physical model could be constructed on how to compute the q-parameter value for the Rényi entropy of black holes (or other strongly gravitating systems) in order to parametrize the non-local effects of the gravitational field, the Rényi statistics can provide a well behaving and additive entropy description of the system which is also compatible with the requirements of equilibrium and the zeroth law of thermodynamics. Similar considerations have been applied recently to describe the relative information entropy measure inside compact domains of an inhomogeneous universe [54], where an explicit geometric model has been proposed to compute the q-parameter of the Rényi entropy in order to measure the effects of the gravitational entanglement problem.
In the case of black hole thermodynamics, we showed that large, asymptotically flat, Schwarzschild black holes can be in stable equilibrium with a thermal heath bath in the Rényi picture, and a Hawking-Page phase transition can occur in the system. This result offers a possible explanation for the problem of cosmic black hole nucleation in the early universe, and this mechanism might be the origin of large or super massive black holes that can be found in most galaxy centers. Many other interesting consequences can be deduced from the Rényi approach, but in this work we aim to achieve a more modest goal. In the present paper, by extending our previous investigations in the problem, we study the thermodynamic, stability and phase transition properties of Kerr black holes within the Rényi model and analyze whether similar results can be obtained to what we have found in the Schwarzschild case. In this analysis the turning point method is applied to the stability problem in both the microcanonical and canonical ensembles, and we will show that stability changes appear in the latter case, which suggests that a Hawking-Page transition and a first order small black hole/large black hole phase transition occur in the system, similar to the one observed for charged and rotating black holes in AdS space [55][56][57][58][59]. This result provides a correspondence between the Kerr-Rényi and the Kerr-AdS-Boltzmann pictures, analogous to the one we reported in the Schwarzschild problem [40,53].
The plan of the paper is as follows. In Sec. 2 we discuss the foundations and motivation of the Rényi approach arising from nonextensive thermodynamics to the black hole problem. In Sec. 3 we introduce the Kerr solution and calculate its thermodynamic quantities within the Rényi statistics. In Sec. 4 we investigate the thermodynamic stability problem of Kerr black holes in the Rényi model by the Poincaré turning point method both in the microcanonical and canonical treatments. We also discuss the question of possible phase transitions in this section. In Sec. 5 the thermodynamic stability problem of Kerr-AdS black holes in the standard Boltzmann case is also presented by the turning point method, and the correspondence between the Kerr-Rényi and the Kerr-AdS-Boltzmann approaches is discussed. In Sec. 6 we summarize our results and draw our conclusions. Throughout this paper we use units such as c = G = = k B = 1.

Rényi approach from nonadditive thermodynamics
By replacing the additivity axiom to the weaker composability in the Shannon-Khinchin axiomatic definition of the entropy function, new type of entropy expressions arise. The composability axiom asserts, roughly speaking, that the entropy S 12 of a compound system consisting of two independent systems should be computable only in terms of the individual entropies S 1 and S 2 . This means that there is a function f (x, y) such that for any independent systems. This property is of fundamental importance, since it implies that an entropic function is properly defined on macroscopic states of a given system, and it can be computed without having any information on the underlying microscopic dynamics. Composability is a key feature to ensure that the entropy function is physically meaningful. In a recent paper [60], based on the concept of composability alone, Abe derived the most general functional form of those nonadditive entropy composition rules that are compatible with homogeneous equilibrium. Assuming that f (S 1 , S 2 ) is a C 2 class symmetric function, Abe showed that the most general, equilibrium compatible composition rule takes the form where H λ is a differentiable function of S and λ ∈ R is a constant parameter. Later on, this result has been extended to non-homogeneous systems as well [10], where not only the entropy, but the energy function is also considered to be nonadditive. The simplest and perhaps the most well-known nonadditive entropy composition rule can be obtained from (2.2) by setting H λ (S) to be the identity function, i.e. H λ (S) = S. In this case Abe's equation becomes which results the familiar Tsallis composition rule with λ = 1 − q [39], and the corresponding entropy definition is given in (1.2). The nonextensive Tsallis statistics is widely investigated in many research fields from natural to social sciences, an updated bibliography on the topic can be found in [61]. This approach has also been studied in the problem of black hole thermodynamics (see e.g. [25] and references therein), and our starting point in considering a more general entropy definition for black holes than the one based on the Boltzamann-Gibbs statistics is also the Tsallis formula. Generalized, nonadditive entropy definitions has been investigated in various problems from high energy physics [13] to DNA analysis [62], and it has been a longstanding problem that the zeroth law of thermodynamics (i.e. the existence of a well defined temperature function in thermal equilibrium) cannot be compatible with nonadditive entropy composition rules. A possible resolution to this problem has been proposed recently by Biró and Ván in [10], where they developed a formulation to determine the most general functional form of those nonadditive entropy composition rules that are compatible with the zeroth law of thermodynamics. They found that the general form is additive for the formal logarithms of the original quantities, which in turn, also satisfy the familiar relations of standard thermodynamics. In particular, for homogeneous systems, they showed that the most general, zeroth law compatible entropy function takes the form which is additive for composition, i.e., 5) and the corresponding zeroth law compatible temperature function can be obtained as where E is the energy of the system. In the case of the Tsallis statistics, it is easy to show that by taking the formal logarithm (2.4) of the Tsallis entropy (1.2), i.e.

L(S
the Rényi expression (1.3) is reproduced, which, unlike the Tsallis formula, is additive for composition.
In the limit of q → 1 (λ → 0), both the Tsallis-and the Rényi entropies recovers the standard Boltzmann-Gibbs description. According to these results, in the present paper, in order to describe the non-Boltzamannian nature of Kerr black holes, we consider the Tsallis statistics as the simplest, nonadditive, parametric but equilibrium compatible extension of the Boltzamann-Gibbs theory which also satisfies Abe's formula. On the other hand, in order to satisfy the zeroth law of thermodynamics, we follow the formal logarithm method of Biró and Ván, and rather the Tsallis description, we consider the Rényi entropy for the thermodynamics of the problem. Since the Rényi definition is additive, it satisfies all laws of thermodynamics, and compared to the Boltzmann picture, it has the advantage of having a free parameter which can be accounted to describe the effects of nonlocality in our approach. The thermodynamics of Schwarzschild black holes in this model has been studied in [40], and the corresponding stability problem has been investigated in [53].

Kerr black holes
The spacetime metric that describes the geometry of a rotating black hole is given by the Kerr solution Here, M is the mass-energy parameter of the black hole and a is its rotation parameter. The thermodynamic quantities of a Kerr black hole can be expressed in terms of its horizon radius r + = M + √ M 2 − a 2 , which is defined by taking ∆ = 0. The Hawking temperature of the black hole horizon is the Bekenstein-Hawking entropy is the angular momentum of the black hole is the angular velocity of the horizon is Ω = a r 2 + + a 2 , (3.6) and the mass-energy parameter can also be written as The heat capacity at constant angular velocity is given by (3.8) and the heat capacity at constant angular momentum is (3.9) C Ω and C J can be written in simpler forms if we normalize them by r 2 + , i.e.
C Ω r 2 and were we also introduced the normalized rotation parameter h [63] as The r + horizon radius exists only for |a| ≤ M , which corresponds to 0 ≤ h ≤ 1. The h = 0 value describes a Schwarzschild black hole, while the limiting value h = 1 belongs to the extreme Kerr black hole case. The heat capacities as functions of h are plotted on Fig. 1. It can be seen that C Ω is statistics. The physical meaning of the λ parameter is connected to the nonextensive and nonlocal nature of the problem, and the Rényi entropy of a general black hole in this picture is given by The zeroth law compatible Rényi temperature is then defined as For the case of a Kerr black hole, the Rényi entropy and the corresponding temperature take the forms S R = 1 λ ln(1 + πλ(r 2 + + a 2 )) (3.14) and .
The heat capacities can be obtained as The heat capacity at constant angular velocity is then , (3.17) while the heat capacity at constant angular momentum takes the form . C ΩR and C JR can also be written in the simpler, normalized forms as before, i.e. and where we also introduced the parameter k = λr 2 + . The heat capacities change their sign depending on the parameter values as plotted on Fig. 2.
For fixed angular momentum, both the Rényi and the Boltzmann entropies of a Kerr black hole are monotonically increasing functions of the mass parameter (Fig. 3). An important difference however, is that while the standard Boltzmann entropy is asymptotically convex (being proportional to M 2 as approaching the static Schwarzschild solution in the large M limit), the Rényi entropy is asymptotically concave, since it increases only logarithmically.
On Fig. 4, we also plotted the temperature-energy relations for fixed angular momentum J = J 0 . As it is well known, there is a maximum temperature in the case of the standard Boltzmann approach. In the smaller mass (low entropy) region of the T H (M ) curve, the heat capacity C J is positive, while at larger masses (high entropy) region it is negative. The two regions correspond to the phases of h > h c and h < h c on Fig. 1. In the Rényi approach, the T R (M ) curve has different behavior depending on the actual parameter values of the angular momentum and λ. It has a local maximum and a local minimum when λJ 0 is smaller than some critical value. In this case three black holes with different masses can coexist at a given temperature. The C JR heat capacity is negative in the phase between the two local extrema and positive otherwise (see also Fig. 7 in Sec. 4). Above the critical value of λJ 0 , the local extrema disappear and the Rényi temperature becomes a monotonically increasing function

Stability analysis
Kaburaki et al. [47] have shown that isolated Kerr black holes are stable against axisymmetric perturbations in the standard picture of black hole thermodynamics. They also investigated the effects of thermal radiation around Kerr black holes from the standpoint of a canonical ensemble. They concluded that slowly rotating black holes (just like Schwarzschild black holes) are unstable in a heat bath, but rapidly rotating holes are less unstable and may even be stable.
As we pointed out, the standard thermodynamic stability criteria based on the Hessian analysis fails when it is applied to black holes because of the nonadditivity of the entropy function, and also of the impossibility to define extensive thermodynamic quantities (see e.g. [51]). These problems can be mainly avoided by the Poincaré turning point method [50], which can divide a linear series of equilibrium into subspaces of unstable and less unstable states of equilibria solely from the properties of the equilibrium sequence without solving any eigenvalue equation. The method has been widely applied to problems in astrophysical and gravitating systems, in particular for the study of the thermodynamic stability of black holes in four and also in higher dimensions [47-49, 51, 52]. A well-organized review on the method has been given by Arcioni and Lozano-Tellechea in [51], which we will mainly follow in the next brief description. Formal proofs of the results can be found in the previous works of Katz [64,65] and Sorkin [66], and a more intuitive interpretation has also been given in [67].
Suppose M is the space of possible configurations of a system in which a given configuration is specified by a point X. The set of independent thermodynamic variables µ i specifies a given ensemble. for all displacements dX. The set of equilibria is a submanifold M eq of the configuration space. Points in M eq can be labeled by the corresponding values of µ i which are called control parameters. Usually all that we know is the explicit expression of the equilibrium Massieu function which is the integral of Eq. (4.1). The explicit expressions for the conjugate variables at equilibrium are calculated as The entropy maximum postulate is a statement about the behaviour of the entropy function along off-equilibrium curves in M, not along equilibrium sequences M eq . Therefore we need an expression for an extended Massieu functionẐ =Ẑ(X ρ , µ i ) where X ρ denote a set of off-equilibrium variables. The equillibrium configurations occur at the values X ρ eq = X ρ eq (µ i ) which are solutions of A stable equilibrium takes place at the point X ρ eq if and only if it is a local maximum ofẐ at fixed µ i . Therefore stable solutions have a matrix with a negative spectrum of eigenvalues. A change of stability takes place when one of these eigenvalues become zero and changes sign. It has been shown that we can obtain information about the sign of an eigenvalue near a point where it vanishes without computing the spectrum of the matrix (4.5). What we have to do is to plot a conjugate variable β a (µ a ) = ∂Zeq ∂µ a against a control parameter µ a along an equilibrium sequence for some fixed a. When a change of stability occurs, called a turning point, the plot of the stability curve β a (µ a ) has a vertical tangent. Figure 5 shows an example of stability curves with a turning point A. The branch with negative slope near this point is always unstable, since one can prove that at least one eigenvalue of the matrix (4.5) is positive. The branch with positive slope near A is more stable. If the spectrum of zero modes is nondegenerate, only one of the eigenvalues changes its sign at the turning point. Therefore the positive slope branch has one unstable mode less than the negative slope branch. Note that what we can see with this method is the existence of instabilities because there can be positive eigenvalues which do not change sign at the turning point. If there is a point in the equilibrium sequence which is fully stable, all equilibria in the sequence are fully stable until the first turning point is reached.

Pure black holes
The case of a black hole isolated from its surroundings can be described in the microcanonical emsemble. The Massieu function is the Rényi entropy (3.14) and the control parameters are M and J. The conjugate variables are β and −α, where β is the inverse of the Rényi temperature (3.15) and α = Ω TR . To study the stability of the black hole we need to plot the stability curves β(M ) at constant J and −α(J) at constant M .
When J = J 0 is a constant, the variables β/ √ J 0 and M/ √ J 0 can be written as functions of the normalized parameter h with constant λJ 0 . On Fig. 6 we plotted the stability curves β(M ) for different values of λJ 0 . For comparison, the stability curve of the standard Kerr black hole (λ = 0) is also plotted on the figure. It can be seen that no curves of the Kerr-Rényi case have a vertical tangent (similar to the standard Kerr result analyzed by Kaburaki, Okamoto and Katz in [47]), and therefore there is no stability change for any M . Like in the λ = 0 case, the λ > 0 curves also diverge asymptotically as we consider the extreme black hole limit h → 1. In the large M limit the black holes approach the static Schwarzschild solution (h = 0). It can also be seen that the λ > 0 stability curves are similar to the Schwarzschild-Rényi case in the large M region (see Fig. 3 of [53]). No vertical tangent occurs in either case. By rotating the figure clockwise with π 2 , the stability curves of the canonical treatment can be obtained, i.e. −M (β) at fixed J. In this case, the λJ0 = 0.01 (green) curve has two vertical tangents denoting the loss and the recovery of stability. In this scenario, up to three black holes with different mass-energy parameters can coexist at a given temperature.
For the standard case, it has been shown that isolated Kerr holes are thermodynamically stable with respect to axisymmetric perturbations. Also, the isolated Schwarzschild black holes have been found to be stable against spherically symmetric perturbations in the Rényi approach. Based on these results, since no turning point occurs on the stability curves in between these two extrema, we can conclude that isolated Kerr black holes are thermodynamically stable against axisymmetric perturbations in the Rényi approach as well.
By looking at Fig. 6, we can also see that there are two points where the tangent of the stability curves with smaller λ (or small angular momentum J 0 ) becomes horizontal. These correspond to the points where the heat capacity at constant J changes its sign through an infinite discontinuity, similar to the Davies point [26] of the standard Kerr black hole case. On Fig. 7 we plotted the lines of constant J in the normalized parameter space of (h, k). Here the λJ 0 = 0.01 line crosses the line where C JR diverges. Along this line the heat capacity C JR changes its sign two times on the way from the (h = 0) Schwarzschild limit to (h = 1) extremality.
When M = M 0 is a constant, α and J/M 2 0 become functions of h with a constant λM 2 0 . On Fig. 8, we plotted the −α(J) stability curves for different values of λM 2 0 . There is no vertical tangent and hence no stability change occurs at any point for any λ. The M = const. lines in the parameter space of (h, k) are plotted on Fig. 9.

Black holes in a heat bath
Let us now consider the black hole in the canonical approach. The canonical ensemble describes the system of a black hole in equilibrium with an infinite reservoir of thermal radiation at constant The stability curves of −M (β) at constant J are simply the π 2 clockwise rotated versions of Fig.  6. We can see that there is a vertical tangent along the stability curve in the standard Boltzmann treatment (λ = 0). The heat capacity C J diverges at this turning point where h = h c . The 0 < h < h c branch of this curve is less stable than the h > h c branch. As Kaburaki, Okamato and Katz have shown [47], one can conclude from this result that since Schwarzschild black holes (h = 0) are unstable in an infinite bath, so are the slowly rotating holes until the h c turning point is reached. Rapidly rotating Kerr black holes, on the other hand, can become stable if the slowly rotating (unstable) holes have only one negative eigenmode, which changes sign at h c .
For the parametrized Rényi case, the behavior of the stability curves changes depending on the value of λJ 0 . We can see that there are two turning points on the λJ 0 = 0.01 stability curve when we rotate the plots of Fig. 6 with π/2 clockwise. These turning points disappear when the value of λJ 0 is larger than a critical value. As a consequence, stability change occurs only when the parameter λ and/or J 0 is sufficiently small. In this case, there are three phases of black holes; small, intermediate unstable and large black holes. The stability property of a rapidly rotating, small black hole is the same as of a slowly rotating, large black hole, which is expected to be stable from continuity requirements to the static solution in the Rényi approach [53].
The stability curves of −α(J) at constant β are plotted on Fig. 10, while the lines of constant β in the parameter space of (h, k) are plotted on Fig. 11. One can see again that the stability curves have no vertical tangent when the λ parameter is sufficiently large, similar to the case of β(M ) with constant J. Vertical tangents to the stability curves appear however when λ and/or β 0 are smaller than some critical value. When a curve has two vertical tangents, there are two turning points where the heat capacity C JR diverges and changes its sign. The black hole changes its stability there in the order from a less unstable state to unstable state and back to a less unstable state again, where the less unstable states may even be stable but not guaranteed.
When λβ 2 0 is less than some critical value, there are two branches of the stability curves −α(J) at constant β 0 . This result is consistent with the fact that in the Schwarzschild-Rényi case there are two black holes with the same β, as it can be seen on Fig. 1 in [53]. The lower branch of the λβ 2 0 = 9 curve terminates at the h = 0 small, static black hole limit. The behavior of this branch is similar to the curve of the λ = 0 Kerr-Boltzmann case. According to these results, we can conclude that small, static or slowly rotating black holes in the Rényi approach are unstable in a heat bath, but fast rotation can stabilize them in a similar way as it is done in the Kerr-Boltzmann case, which has been shown by Kaburaki, Okamoto and Katz [47]. The upper branch of the λβ 2 0 = 9 stability curve belongs to larger mass black holes and terminates at the large, static black hole limit when h → 0. There is no vertical tangent in this branch so the corresponding rotating black holes have the same stability property as the large, static black holes in a heat bath, i.e. they are stable.

Phase Transitions
We have shown in our earlier paper [53] that a Hawking-Page phase transition can be observed for static black holes in the Rényi approach. Previously, it had also been shown [57][58][59] that Kerr-AdS black holes exhibit a first order small black hole/large black hole (SBH/LBH) phase transition in the canonical ensemble. In this subsection we will study the question of possible phase transitions of Kerr black holes in the Rényi model.
The behavior of the free energy function F = M − T R S R for λJ 0 = 0.01, 0.02 and 0.1 at constant J is displayed on Fig. 12. We can see that all curves cross the horizontal axis. Small black holes with lower temperature possess positive free energy, while larger black holes with higher temperature possess negative free energy. One can, therefore, expect a Hawking-Page transition between the thermal gas phase with angular momentum, and the large black hole state, which would be locally stable according to our analysis in Sec. 4.2. It is generally assumed that F ≈ 0 for a thermal gas, so the phase transition occurs around the temperature where the free energy of the black hole becomes zero. An SBH/LBH phase transition can also be observed for Kerr black holes in the Rényi model when we enlarge the λJ 0 = 0.01 curve of Fig. 12 on Fig 13. The swallowtail behavior of the free energy function is a typical sign of a first order transition between the SBH and LBH phases. There are three branches on the picture: small, lower temperature holes; large, higher temperature holes; and also intermediate, unstable black holes. There is a coexistence point of small and large black holes where the SBH/LBH transition occurs. The mass and entropy functions are discontinuous at this point which indicates that the phase transition is a first order kind. By increasing λ, the swallowtail behavior disappears, as it can be seen on the λJ 0 = 0.02 curve on Fig. 12. This suggests the existence of a critical point where the phase transition becomes second order.

Kerr-AdS black holes
In order to compare the obtained stability results of the Kerr-Rényi model to the Kerr-AdS-Boltzamann case in the Poincaré approach (analogous to the Schwarzschild problem), in this section we present the Poincaré stability analysis of the Kerr-AdS-Boltzamann case as well. The Kerr-AdS black hole metric is described by where ∆ r = (r 2 + a 2 ) 1 + r 2 l 2 − 2mr, ∆ θ = 1 − a 2 cos 2 θ l 2 , The thermodynamic quantities are written in terms of a, l and the horizon radius r + , which is obtained by solving ∆ = 0. The Hawking temperature of the horizon is given by while the Bekenstein-Hawking entropy of the black hole is S = π a 2 + r 2 3) The angular momentum of a Kerr-AdS black hole is the angular velocity of the horizon is Ω = a l 2 r 2 + + l 2 r 2 + + a 2 , (5.5) and the mass-energy parameter of the black hole can be re-expressed as The heat capacity at constant angular velocity can be computed as and the heat capacity at constant angular momentum takes the form where X = −l 4 r 4 + + 3l 2 r 6 + + a 6 l 2 + r 2 + + a 4 3l 4 + 13l 2 r 2 + + 6r 4 + +a 2 6l 4 r 2 + + 23l 2 r 4 + + 9r 6 + .
Here we introduced the normalized parameters p ≡ |a| r + , and s ≡ l r + . (5.9) The heat capacities C Ω and C J change their signs depending on the values of p and s. The parameter space of (p, s) can be divided into 4 regions depending on the signs of C Ω and C J as shown on Fig. 14.
Similarly to the Kerr-Rényi case in Sec. 4, the thermodynamic stability problem of Kerr-AdS black holes can also be analyzed by the Poincaré turning point method. The stability curves of the two systems are qualitatively similar. For the study of the Kerr-AdS black hole problem we will use the following normalized variables,β AdS black holes occurs only when the temperature is higher than a certain critical value. Black holes with slightly higher temperature than the critical one have an unstable branch between two, more stable branches. There is another critical temperature above which a cusp appears on the stability curve at (J, holes are unstable in the canonical ensemble.  Figure 19. Curves of the conjugate variable −α(J) at fixed β for Kerr-AdS black holes in the canonical ensemble. The curves ofβ = 3 (brown),β = 3.63 (green),β = 3.7 (red) andβ = 4 (blue) are plotted. Thẽ β = 3.63 (green) curve has two turning points, while theβ = 3 (brown) curve has two branches and the lower branch has a turning point.
As it can be clearly seen from the analysis above, the thermodynamic properties of the Kerr-Rényi and the Kerr-AdS-Boltzmann models are very similar. In the static case, we have obtained a simple relation between the entropy parameter λ and the AdS curvature parameter l for black holes with identical horizon temperatures [40]. By assuming the same condition for stationary black holes augmented with the assumption of identical horizon angular velocity, we can derive analogous relations between the (h, k) and (p, s) parameters for rotating black holes by solving the following equationŝ where we normalized the quantities by the horizon radius r + aŝ T = T r + ,Ω = Ωr + . (5.12) As a result, a quantitative analogy between the Kerr-Rényi and Kerr-AdS-Boltzmann pictures of black hole thermodynamics can be given by the parameter equations where These equations provide a very interesting correspondence between the two approaches.

Summary and Conclusions
In this paper we investigated the thermodynamic and stability properties of Kerr black holes described by the parametric, equilibrium-and zeroth law compatible Rényi entropy function. The corresponding problem of static Schwarzschild black holes has been analyzed in [40,53], where interesting similarities have been found to the picture of standard black hole thermodynamics in asymptotically AdS space. In particular, a stability change and a Hawking-Page transition have been identified, which motivated us to extend our investigations to the present (3+1)-dimensional, rotating problem as well.
The novel results of this work are the following. We derived the temperature and heat capacities of a Kerr black hole in the Rényi approach, and found that the global maximum of the temperatureenergy curve at a fixed angular momentum in the standard description becomes only a local maximum in the Rényi model. In the thermodynamic stability analysis we investigated both the microcanonical and the canonical ensembles. We have plotted the stability curves of the Boltzmann-Gibbs and Rényi entropy models, and showed that no stability change occurs for isolated black holes in either case. From this result, we concluded that, similarly to the standard Boltzmann case, isolated Kerr black holes are thermodynamically stable with respect to axisymmetric perturbations in the Rényi approach.
In case when the black holes are surrounded by a bath of thermal radiation in the canonical picture, we found that, in contrast to the standard Boltzmann case, slowly rotating Kerr black holes can be in stable equilibrium with thermal radiation at a fixed temperature if the number of negative eigenmodes of the stability matrix is one. We showed that fast rotating black holes have similar stability properties to slowly rotating ones, and there may also exist intermediate size, unstable black holes. We also analyzed the question of possible phase transitions in the canonical picture, and found that, in addition to a Hawking-Page transition, a first order small black hole/large black hole phase transition occurs in a very similar fashion as in AdS space. These findings indicate that there is a similarity between the Kerr-Rényi and Kerr-AdS-Boltzamann models, analogous to the one that we found in the static case. Based on this result we also investigated the Poincaré stability curves of Kerr-AdS black holes in the standard Boltzmann picture, and confirmed this similarity by obtaining simple algebraic relations between the parameters of the two approaches with identical surface temperature and angular velocity.
The above results may be relevant in many aspects of black hole physics. Our main motivation in the first place was to consider a statistical model to the nonextensive and nonlocal nature of black hole thermodynamics, where we do not assume a priori that the classical, additive Boltzamann statistics can describe this strongly gravitating system. The Rényi form of the black hole entropy includes a parameter λ, which seems to be a good candidate to incorporate the effects of the long-range type behavior of the gravitational field, while also being additive and satisfying both the equilibrium compatibility and the zeroth law's requirements. A specific model on how to compute the λ parameter value for the black hole problem is yet to be developed, but a similar approach has been considered in [54] to describe the mutual information between spatially separated, compact domains of an inhomogeneous universe that are entangled via the gravitational field equations. In that work, as an effective model, the λ parameter of the Tsallis/Rényi relative entropy has been defined in a geometric way in order to describe the causal connection between the domain and its surroundings during the cosmic evolution. Since black holes are essentially the final states of cosmic structure formation, one can expect that the two directions might be connected somehow in the nonlinear regime of matter collapse.
As a different direction, it is also interesting to mention that by considering the Boltzamann picture in the standard description, the Bekenstein-Hawking entropy has a nontrivial nonadditive property which also satisfies Abe's formula. In the case of Schwarzschild black holes this nonadditivity reads as H λ (S) = √ S and for Kerr black holes H λ (S) = S √ S−a 2 π with λ = 0. The corresponding thermodynamic and stability problems (by also applying the formal logarithm method) has been studied in [68] and [69], respectively.
In the present parametric approach however, the most important result is the confirmation of a stability change and the Hawking-Page transition of Kerr black holes in the Rényi model. As we discussed in the introduction, this phenomena has many interesting connections with other open problems in theoretical physics, e.g. the cosmic nucleation of matter into black holes in the early universe, or due to the similarity to the AdS-Boltzmann problem, it may also be connected to the AdS/CFT correspondence and related phenomena. A further motivation arises from a different possible interpretation of the parametric Rényi picture originating from finite size reservoir effects in the canonical ensemble. In a recent paper [70], Biró showed that from the requirement of zero mutual information between a finite subsystem and a finite reservoir in thermodynamic equilibrium, the Tsallis-and Rényi entropy formulas arise very naturally. Although we haven't worked out the details of this approach yet, it provides a nice possible interpretation of our findings as placing a black hole into a finite heat bath in the canonical approach instead of an infinite reservoir (which is an idealistic model), and require zero mutual information between the black hole and the reservoir in thermal equilibrium. In this situation the system is dominated by the bath, and Biró showed that the entropy parameter in this case is proportional to the heat capacity of the bath as λ = 1/C 0 , where instead of the classical infinite approximation, the heat capacity of the bath is a large but finite constant C 0 . This approach has been investigated e.g. for the case when a quark-gluon plasma system is connected to a finite heat bath in [71].
In conclusion, several interesting consequences can be deduced from the Rényi approach to black hole thermodynamics which is motivated by various physical considerations. Parametric corrections to the black hole entropy problem also arise from quantum considerations, e.g. from string theory, loop quantum gravity or other semi-classical theories (see e.g. [72] and references therein), and we expect that other parametric situations are also possible which might be connected to the parametric Rényi description.