---
author:
- 'Viktor G. Czinner'
- and Hideo Iguchi
title: 'Thermodynamics, stability and Hawking–Page transition of Kerr black holes from Rényi statistics'
---

Introduction
============

Gravitational phase transitions, in particular the ones connected to black hole thermodynamics, are essential constituents of many open problems in modern theoretical physics. The Hawking–Page phase transition [@Hawking:1982dh] of black holes in anti-de Sitter space is one of the most important ones due to its role in the AdS/CFT correspondence [@Maldacena:1997re; @Witten:1998qj] and also in related phenomena of confinement/deconfinement transitions at finite temperature in various gauge theories [@MMT1; @MMT2]. Because of the different background geometry, asymptotically flat black holes have different stability properties than AdS ones, and in the standard black hole thermodynamic picture [@Bekenstein:1973ur; @Bardeen:1973gs; @Hawking:1974sw; @Hawking:1976de], they mostly tend to be unstable for any large masses when surrounded by an infinite bath of thermal radiation. A Hawking–Page transition does not occur under these conditions, and a cosmic black hole nucleation is not present in asymptotically flat spacetimes. Apart from the gravity interest, the above phenomenon is interesting from a thermodynamic viewpoint as well, and for a clear understanding of the physics behind, the underlying theory of black hole thermodynamics is also necessary to be well understood. In the past 40 years, after the foundations of the standard thermodynamic theory of black holes [@Bekenstein:1973ur; @Bardeen:1973gs; @Hawking:1974sw; @Hawking:1976de], numerous achievements have been made in the field. In spite of the active research and successes however, there still are some unsettled and important issues which could not be resolved satisfactorily so far. From a classical thermodynamic perspective, one of the most interesting ones is the nonextensive nature of black holes and the corresponding problem of thermodynamic stability.

A basic group of physical quantities in classical thermodynamics is the group called extensive variables $X$ (like energy, entropy, etc.), where it is assumed that these quantities are *additive* for composition, i.e. $X_{12}=X_1+X_2$ when thermodynamic systems are joined together [@PhysRevE.83.061147]. On the other hand, it is also customary to assume that these quantities characterize the system down to the smallest scales [@Mackey], i.e. when working with finite densities: $$\varrho_X = \lim_{n\rightarrow\infty}\frac{1}{n}\sum_{i=1}^n X_i < \infty,$$ 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 [@T2; @B1]. 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}=-\sum 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. [@Landsberg1984; @cg1; @cg2; @cg3; @cg4; @cg5; @cg6; @cg7; @cg8; @cg9] and references therein). In fact, even as early as 1902, Gibbs already pointed out in his statistical mechanics book [@Gibbs], that systems with divergent partition function lie outside the validity of Boltzamann–Gibbs theory. He explicitly mentions gravitation as an example (see e.g. [@Tsallis:2012js] 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 [@Davies:1978mf] very early on after the thermodynamic theory of black holes had been formulated [@Bekenstein:1973ur; @Bardeen:1973gs; @Hawking:1974sw; @Hawking:1976de], and the corresponding thermodynamic and stability problem has been investigated several times with various approaches (see e.g. [@landsberg1980entropies; @bishop1987thermodynamics; @Landsberg1984; @pavon1986some; @1993Natur.365..103M; @Gour:2003pd; @Oppenheim:2002kx; @Pesci:2006sb; @Aranha:2008ni; @Tsallis:2012js] and references therein). The general theory of nonadditive thermodynamics has also advanced significantly in the past few decades (see e.g. [@T2; @Landsberg] 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 [@Shannon] and Khinchin [@Khinchin]) to the weaker *composability* requirement, new possible functional forms of the entropy may arise [@Tempesta]. 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 [@Tsallis:1987eu] as: $$\label{TsS}
 S_T=\frac{1}{1-q}\sum_i(p^q_i-p_i),$$ where $p_i$ are the probabilities of the microscopic states of the system, and $q\in\mathbb{R}$ is the so-called nonextensivity parameter. In the limit of $q\rightarrow 1$, $S_T$ reproduces the standard Boltzamann–Gibbs result, however, in the case when $q\neq 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. [@Tsallis:2012js] 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 [@PhysRevE.83.061147]. 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 [@PhysRevE.83.061147], 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 [@renyi1959dimension; @renyi1970probability], defined as $$\label{Srenyi}
S_R=\frac{1}{1-q}\ln\sum_ip^q_i,$$ which had been proposed earlier by the Hungarian mathematician Alfréd Rényi in 1959 [@renyi1959dimension].

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 [@Biro:2013cra]. 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 (\[TsS\]) 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 [@Biro:2013cra], we found that the temperature-horizon 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 [@Bekenstein:1980jp] for the $q$-parameter value in the Rényi entropy of micro black holes ($q \geq 1 + 2/\pi^2$), which was surprisingly close to other $q$-parameter fits from very distant and unrelated physical phenomena, e.g. cosmic ray spectra [@Bek1; @Bek2], and power-law distribution of quarks coalescing to hadrons in high energy accelerator experiments [@BU].

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. [@Kaburaki:1993ah; @Katz:1993up; @kaburaki1996critical]. In these works the so-called “Poincaré turning point method of stability” [@poincare1885equilibre] 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 [@Arcioni:2004ww] and to determine the conditions of stability for equilibrium configurations of charged black holes surrounded by quintessence [@AzregAinou:2012hy]. In section \[stab\] 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 [@Czinner:2015eyk]. 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 [@Hawking:1982dh]. 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 [@CzM], 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 [@Chamblin:1999tk; @Chamblin:1999hg; @Caldarelli:1999xj; @Tsai:2011gv; @Altamirano:2014tva]. 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 [@Biro:2013cra; @Czinner:2015eyk].

The plan of the paper is as follows. In Sec. \[sec:Renyi\] we discuss the foundations and motivation of the Rényi approach arising from nonextensive thermodynamics to the black hole problem. In Sec. \[sec:Kerr\] we introduce the Kerr solution and calculate its thermodynamic quantities within the Rényi statistics. In Sec. \[sec:Stability\] 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. \[sec:Kerr-AdS\] 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. \[sec:summary\] we summarize our results and draw our conclusions. Throughout this paper we use units such as $c=G=\hbar=k_B=1$.

Rényi approach from nonadditive thermodynamics {#sec:Renyi}
==============================================

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 $$S_{12} = f (S_1, S_2),$$ 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 [@abe2001general], 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 $$\label{eq:Abe}
H_{\lambda}(S_{12})=H_{\lambda}(S_1)+H_{\lambda}(S_2)+\lambda H_{\lambda}(S_1)H_{\lambda}(S_2),$$ where $H_{\lambda}$ is a differentiable function of $S$ and $\lambda\in\mathbb{R}$ is a constant parameter. Later on, this result has been extended to non-homogeneous systems as well [@PhysRevE.83.061147], 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 (\[eq:Abe\]) by setting $H_{\lambda}(S)$ to be the identity function, i.e. $H_{\lambda}(S)=S$. In this case Abe’s equation becomes $$\label{eq:tsallis_composition}
S_{12}=S_1+S_2+\lambda S_1 S_2,$$ which results the familiar Tsallis composition rule with $\lambda=1-q$ [@Tsallis:1987eu], and the corresponding entropy definition is given in (\[TsS\]). 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 [@Tsallis-bib]. This approach has also been studied in the problem of black hole thermodynamics (see e.g. [@Tsallis:2012js] 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 [@B1] to DNA analysis [@DNA], 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 [@PhysRevE.83.061147], 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 $$\label{eq:formallog}
L(S)=\frac{1}{\lambda}\ln[1+\lambda H_{\lambda}(S)],$$ which is additive for composition, i.e., $$L(S_{12})=L(S_{1})+L(S_{2}),$$ and the corresponding zeroth law compatible temperature function can be obtained as $$\frac{1}{T}=\frac{\partial L(S(E))}{\partial E},$$ 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 (\[eq:formallog\]) of the Tsallis entropy (\[TsS\]), i.e. $$L(S_T)=\frac{1}{1-q}\ln\left[1+(1-q)S_T\right] \equiv S_R,$$ the Rényi expression (\[Srenyi\]) is reproduced, which, unlike the Tsallis formula, is additive for composition. In the limit of $q\rightarrow 1$ ($\lambda \rightarrow 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 [@Biro:2013cra], and the corresponding stability problem has been investigated in [@Czinner:2015eyk].

Kerr black holes {#sec:Kerr}
================

The spacetime metric that describes the geometry of a rotating black hole is given by the Kerr solution $$\begin{aligned}
 ds^2 & = & -dt^2 + \frac{2 M r}{\Sigma} \left( dt -a \sin^2 \theta d \phi \right)^2 
 + \frac{\Sigma}{\Delta} dr^2 + \Sigma d \theta^2  \nonumber \\
 & &+ (r^2 +a^2) \sin \theta d \phi^2\ ,\end{aligned}$$ where $$\Sigma = r^2 + a^2 \cos^2 \theta , ~~~ \Delta = r^2 + a^2 - 2Mr .$$ 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 + \sqrt{M^2 - a^2}$, which is defined by taking $\Delta = 0$. The Hawking temperature of the black hole horizon is $$T_H = \frac{1}{2\pi}\left[ \frac{r_{+}}{r_{+}^2 + a^2} - \frac{1}{2 r_{+}} \right],$$ the Bekenstein-Hawking entropy is $$S_{BH} = \pi (r_{+}^2 + a^2),$$ the angular momentum of the black hole is $$J = \frac{a}{2 r_{+}}(r_{+}^2 + a^2),$$ the angular velocity of the horizon is $$\Omega = \frac{a}{r_{+}^2 + a^2},$$ and the mass-energy parameter can also be written as $$M = \frac{r_{+}^2 + a^2}{2 r_{+}}.$$ The heat capacity at constant angular velocity is given by $$C_{\Omega} = T_H \left( \frac{\partial S_{BH}}{\partial T_H} \right)_\Omega 
 = \frac{2 \pi r_{+}^2(a^2 - r_{+}^2)}{r_{+}^2 +a^2},$$ and the heat capacity at constant angular momentum is $$C_J = T_H \left( \frac{\partial S_{BH}}{\partial T_H} \right)_J 
 = \frac{2 \pi (r_{+}^2 -a^2)(r_{+}^2 + a^2)^2}{3 a^4 + 6 r_{+}^2 a^2 - r_{+}^4}\ .$$ $C_\Omega$ and $C_J$ can be written in simpler forms if we normalize them by $r_{+}^2$, i.e. $$\frac{C_{\Omega}}{r_{+}^2} = - \frac{2 \pi (1 - h^2)}{h^2 + 1 }$$ and $$\frac{C_{J}}{r_{+}^2} =  \frac{2 \pi (1 - h^2 )(h^2 + 1)^2}{3 h^4 + 6 h^2 - 1},$$ were we also introduced the normalized rotation parameter $h$ [[@okamoto1990thermodynamical]]{} as $$h \equiv \frac{|a|}{r_{+}} .$$ The $r_{+}$ horizon radius exists only for $|a| \le M$, which corresponds to $0 \le h \le 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. \[fig:heat\_capacity\_Kerr\]. It can be seen that $C_\Omega$ is negative for $ 0\le h<1$ and $C_J$ diverges at $h_c = \sqrt{\frac{2}{3}\sqrt{3} - 1}$, where a pole occurs. $C_J$ is negative for $h<h_c$ and positive for $h>h_c$ values. The heat capacities coincide at the limit values $h=0$ and 1.

![Plots of the heat capacities $C_J$ (red solid line) and $C_\Omega$ (blue dotted line) against $h$. $C_J$ diverges at $h_c$ where a pole occurs.[]{data-label="fig:heat_capacity_Kerr"}](figure1.eps)

The Rényi entropy function of black holes can be obtained by taking the formal logarithm of the Bekenstein-Hawking entropy, which in our non-Boltzmannian approach follows the nonadditive Tsallis statistics. The physical meaning of the $\lambda$ 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 $$S_R = \frac{1}{\lambda} \ln (1 + \lambda S_{BH}).$$ The zeroth law compatible Rényi temperature is then defined as $$T_R = \frac{1}{\partial S_R / \partial M} = T_H ( 1 + \lambda S_{BH} ).$$

For the case of a Kerr black hole, the Rényi entropy and the corresponding temperature take the forms $$\label{eq:Renyi_entropy}
 S_R = \frac{1}{\lambda} \ln ( 1 + \pi \lambda (r_{+}^2 + a^2) )$$ and $$\label{eq:Renyi_temperature}
 T_R = \frac{(1 + \pi \lambda (r_{+}^2 + a^2))(r_{+}^2 - a^2)}{4 \pi r_{+} (r_{+}^2 + a^2)}.$$ The heat capacities can be obtained as $$C_{R} = T_R \left( \frac{\partial S_R}{\partial T_R} \right) 
 = \frac{C_{BH}}{1 + \lambda(S_{BH} + C_{BH})},$$ where $C_{BH} = T_H \left( \frac{\partial S_{BH}}{\partial T_H} \right)$. The heat capacity at constant angular velocity is then $$C_{\Omega R} = \frac{2 \pi r_{+}^2(a^2 - r_{+}^2)}{r_{+}^2 + a^2 
 + \pi \lambda (a^4 - r_{+}^4 + 4 a^2 r_{+}^2)},$$ while the heat capacity at constant angular momentum takes the form $$C_{JR} = \frac{2 \pi (r_{+}^2 - a^2)(a^2 + r_{+}^2)^2}{3 a^4 +6 a^2 r_{+}^2 - r_{+}^4 
 + \pi \lambda (a^2 + r_{+}^2)(a^4 +6 r_{+}^2 a^2 +r_{+}^4)}.$$ $C_{\Omega R}$ and $C_{JR}$ can also be written in the simpler, normalized forms as before, i.e. $$\frac{C_{\Omega R}}{r_{+}^2} = \frac{2 \pi (h^2 -1)}{h^2 + 1 + \pi k (h^4+ 4 h^2 -1)},$$ and $$\frac{C_{JR}}{r_{+}^2} = - \frac{2 \pi (h^2 - 1)(h^2 + 1)^2}{3 h^4 + 6 h^2 - 1 
 + \pi k (h^2 + 1)(h^4 + 6 h^2 + 1)},$$ where we also introduced the parameter $$k = \lambda r_{+}^2.$$ The heat capacities change their sign depending on the parameter values as plotted on Fig. \[fig:heat\_capacity\].

![Phase diagram of Kerr black holes with Rényi entropy. On the dashed curve $C_{JR}$ diverges while on the dotted curve $C_{\Omega R}$ diverges. In region I, $C_{\Omega R}<0$ and $C_{JR} <0$, in region II, $C_{\Omega R}<0$ and $C_{JR} >0$, while in region III, $C_{\Omega R}>0$ and $C_{JR} >0$.[]{data-label="fig:heat_capacity"}](figure2.eps)

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. \[fig:MS\_Kerr\]). 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.]{}]{}

![Plots of the Rényi entropy as a function of the mass-energy parameter at fixed $J = J_0$ for the parameter values $\lambda J_0 = 0.05$, 0.04, 0.03, 0.02 and 0.01 starting from the bottom curve respectively. The top, bold curve belongs to the standard Boltzmann entropy of the Kerr black hole.[]{data-label="fig:MS_Kerr"}](figure3.eps)

On Fig. \[fig:MT\_Kerr\], 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. \[fig:heat\_capacity\_Kerr\]. In the Rényi approach, the $T_R(M)$ curve has different behavior depending on the actual parameter values of the angular momentum and $\lambda$. It has a local maximum and a local minimum when $\lambda 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. \[fig:hc\_J0\] in Sec. \[sec:Stability\]). Above the critical value of $\lambda J_0$, the local extrema disappear and the Rényi temperature becomes a monotonically increasing function of the energy. From these properties one can expect a stability change and a thermodynamic phase transition of Kerr black holes in the canonical ensemble. We will investigate this question in the next section.

![Plots of the Rényi temperature against the mass-energy parameter of Kerr black holes at fixed $J = J_0$ with $\lambda J_0 = 0.05$, 0.04, 0.03, 0.02 and 0.01 parameter values starting from the top respectively. The bottom, bold curve belongs to the standard Hawking temperature of the black hole in the Boltzmann model.[]{data-label="fig:MT_Kerr"}](figure4.eps)

Stability analysis {#stab}
==================

Kaburaki [*et al.*]{} [@Kaburaki:1993ah] 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.

\[sec:Stability\] 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. [@Arcioni:2004ww]). These problems can be mainly avoided by the Poincaré turning point method [@poincare1885equilibre], 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 [@Kaburaki:1993ah; @Katz:1993up; @kaburaki1996critical; @Arcioni:2004ww; @AzregAinou:2012hy]. A well-organized review on the method has been given by Arcioni and Lozano-Tellechea in [@Arcioni:2004ww], which we will mainly follow in the next brief description. Formal proofs of the results can be found in the previous works of Katz [@1978MNRAS.183..765K; @1979MNRAS.189..817K] and Sorkin [@Sorkin:1981jc], and a more intuitive interpretation has also been given in [@1994PhLA..185...21K].

Suppose $\mathcal{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 ${\mu^i}$ specifies a given ensemble. Let $Z$ be the corresponding Massieu function. Both ${\mu^i}$ and $Z$ are functions on $\mathcal{M}$. Equilibrium states occur at points in $\mathcal{M}$ which are extrema of $Z$ under displacements $dX$ for which $d \mu^i = 0$. At equilibrium the Massieu function depends only on the values of the thermodynamic variables $\mu^i$. We can define conjugate variables $\beta_i$ such that $$dS = \beta_i d \mu^i
 \label{eq:dS}$$ for all displacements $dX$. The set of equilibria is a submanifold $\mathcal{M}_{eq}$ of the configuration space. Points in $\mathcal{M}_{eq}$ can be labeled by the corresponding values of $\mu_i$ which are called control parameters. Usually all that we know is the explicit expression of the equilibrium Massieu function $$Z_{eq} = Z (\mu^i),$$ which is the integral of Eq. (\[eq:dS\]). The explicit expressions for the conjugate variables at equilibrium are calculated as $$\beta_i (\mu^i) = \frac{\partial Z_{eq}}{\partial \mu^i}.$$ The entropy maximum postulate is a statement about the behaviour of the entropy function along off-equilibrium curves in $\mathcal{M}$, not along equilibrium sequences $\mathcal{M}_{eq}$. Therefore we need an expression for an extended Massieu function $\hat{Z} = \hat{Z}(X^\rho, \mu^i)$ where $X^\rho$ denote a set of off-equilibrium variables. The equillibrium configurations occur at the values $X^\rho_{eq} = X^\rho_{eq}(\mu^i)$ which are solutions of $$\frac{\partial \hat{Z}}{\partial X^\rho} = 0.$$ A stable equilibrium takes place at the point $X^\rho_{eq}$ if and only if it is a local maximum of $\hat{Z}$ at fixed $\mu^i$. Therefore stable solutions have a matrix $$\left. \frac{\partial^2 \hat{Z}(X^\tau, \mu^i)}{\partial X^\rho \partial X^\sigma}\right|_{X^\tau = X^\tau_{eq}(\mu^i)}
 \label{eq: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 (\[eq:matrix\]). What we have to do is to plot a conjugate variable $\beta_a (\mu^a) = \frac{\partial Z_{eq}}{\partial \mu^a}$ against a control parameter $\mu^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 $\beta_a (\mu^a)$ has a vertical tangent. Figure \[fig:turning\] 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 (\[eq:matrix\]) 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.

![Example of stability curves. A conjugate variable $\beta_a$ against a control parameter $\mu_a$ is plotted. The point $A$ is a turning point, where a stability change can occur. At point $B$ a stability change doesn’t occur even if the slope changes its sign there.[]{data-label="fig:turning"}](figure5.eps)

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 (\[eq:Renyi\_entropy\]) and the control parameters are $M$ and $J$. The conjugate variables are $\beta$ and $-\alpha$, where $\beta$ is the inverse of the Rényi temperature (\[eq:Renyi\_temperature\]) and $\alpha = \frac{\Omega}{T_R}$. To study the stability of the black hole we need to plot the stability curves $\beta(M)$ at constant $J$ and $-\alpha(J)$ at constant $M$.

When $J=J_0$ is a constant, the variables $\beta/\sqrt{J_0}$ and $M/\sqrt{J_0}$ can be written as functions of the normalized parameter $h$ with constant $\lambda J_0$. On Fig. \[fig:micro\_bM\] we plotted the stability curves $\beta(M)$ for different values of $\lambda J_0$. For comparison, the stability curve of the standard Kerr black hole ($\lambda = 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 [@Kaburaki:1993ah]), and therefore there is no stability change for any $M$. Like in the $\lambda = 0$ case, the $\lambda>0$ curves also diverge asymptotically as we consider the extreme black hole limit $h\rightarrow 1$. In the large $M$ limit the black holes approach the static Schwarzschild solution ($h=0$). It can also be seen that the $\lambda > 0$ stability curves are similar to the Schwarzschild-Rényi case in the large $M$ region (see Fig. 3 of [@Czinner:2015eyk]).

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.

![Curves of the conjugate variable $\beta(M)$ at fixed $J$ in the microcanonical treatment. The $\lambda = 0$ (black) curve represents the stability curve of the black hole in the standard thermodynamic approach, while the $\lambda J_0 = 0.01$ (green), $\lambda J_0 = 0.02$ (red) and $\lambda J_0 =0.1$ (blue) curves are the stability curves within the Rényi approach. No vertical tangent occurs in either case. By rotating the figure clockwise with $\frac{\pi}{2}$, the stability curves of the canonical treatment can be obtained, i.e. $-M(\beta)$ at fixed $J$. In this case, the $\lambda J_0 = 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.[]{data-label="fig:micro_bM"}](figure6.eps)

By looking at Fig. \[fig:micro\_bM\], we can also see that there are two points where the tangent of the stability curves with smaller $\lambda$ (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 [[@Davies:1978mf]]{} of the standard Kerr black hole case. On Fig. \[fig:hc\_J0\] we plotted the lines of constant $J$ in the normalized parameter space of $(h,k)$. Here the $\lambda 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.

![Plots of $J = \mbox{const.}$ curves for $\lambda J_0 = 0.01$ (green), $\lambda J_0 = 0.02$ (red) and $\lambda J_0 =0.1$ (blue) parameter values on the $(h,k)$ space.[]{data-label="fig:hc_J0"}](figure7.eps)

When $M=M_0$ is a constant, $\alpha$ and $J/M_0^2$ become functions of $h$ with a constant $\lambda M_0^2$. On Fig. \[fig:micro\_maJ\], we plotted the $-\alpha(J)$ stability curves for different values of $\lambda M_0^2$. There is no vertical tangent and hence no stability change occurs at any point for any $\lambda$. The $M=const.$ lines in the parameter space of $(h,k)$ are plotted on Fig. \[fig:hc\_M0\].

![Curves of the conjugate variable $-\alpha(J)$ at fixed $M$. The $\lambda = 0$ (black) curve represents the stability curve of a Kerr black hole in the standard approach within the microcanonical treatment. The $\lambda M_0^2 = 0.01$ (green), $\lambda M_0^2 = 0.08$ (red) and $\lambda M_0^2 =0.3$ (blue) curves are the stability curves of the Rényi approach. No vertical tangents occurs.[]{data-label="fig:micro_maJ"}](figure8.eps)

![Plots of $M = \mbox{const.}$ curves for $\lambda M_0^2 = 0.01$ (green), $\lambda M_0^2 = 0.08$ (red) and $\lambda M_0^2 =0.3$ (blue) parameter values on the $(h,k)$ space.[]{data-label="fig:hc_M0"}](figure9.eps)

Black holes in a heat bath {#sec: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 temperature. The Massieu function in this case is $$\Psi(\beta,J) = S_R - \beta M = -\beta F,$$ where $F = M - T_R S_R$ is the Helmholtz free energy. The conjugate variables of the control parameters are $-M$ and $-\alpha$. To study the stability of the black hole we need to plot the stability curves $-M(\beta)$ at constant $J$ and $-\alpha(J)$ at constant $\beta$.

The stability curves of $-M(\beta)$ at constant $J$ are simply the $\frac{\pi}{2}$ clockwise rotated versions of Fig. \[fig:micro\_bM\]. We can see that there is a vertical tangent along the stability curve in the standard Boltzmann treatment ($\lambda = 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 [@Kaburaki:1993ah], 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 $\lambda J_0$. We can see that there are two turning points on the $\lambda J_0 = 0.01$ stability curve when we rotate the plots of Fig. \[fig:micro\_bM\] with $\pi/2$ clockwise. These turning points disappear when the value of $\lambda J_0$ is larger than a critical value. As a consequence, stability change occurs only when the parameter $\lambda$ 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 [[@Czinner:2015eyk]]{}.

The stability curves of $-\alpha(J)$ at constant $\beta$ are plotted on Fig. \[fig:cano\_maJ\], while the lines of constant $\beta$ in the parameter space of $(h,k)$ are plotted on Fig. \[fig:hc\_beta0\]. One can see again that the stability curves have no vertical tangent when the $\lambda$ parameter is sufficiently large, similar to the case of $\beta(M)$ with constant $J$. Vertical tangents to the stability curves appear however when $\lambda$ and/or $\beta_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 $\lambda \beta_0^2$ is less than some critical value, there are two branches of the stability curves $-\alpha(J)$ at constant $\beta_0$. This result is consistent with the fact that in the Schwarzschild-Rényi case there are two black holes with the same $\beta$, as it can be seen on Fig. 1 in [@Czinner:2015eyk]. The lower branch of the $\lambda \beta_0^2=9$ curve terminates at the $h = 0$ small, static black hole limit. The behavior of this branch is similar to the curve of the $\lambda = 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 [@Kaburaki:1993ah]. The upper branch of the $\lambda \beta_0^2=9$ stability curve belongs to larger mass black holes and terminates at the large, static black hole limit when $h \rightarrow 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.]{}

![Curves of the conjugate variable $-\alpha(J)$ at fixed $\beta$ in the canonical approach. The $\lambda = 0$ (black) curve describes the standard thermodynamic approach. The $\lambda \beta_0^2 = 9$ (brown), $\lambda \beta_0^2 = 12.7$ (green), $\lambda \beta_0^2 = 13.4$ (red) and $\lambda \beta_0^2 = 20$ (blue) curves are the stability curves of the Rényi model. The $\lambda \beta_0^2 = 12.7$ (green) curve has two turning points, while the stability curve of $\lambda \beta_0^2 = 9$ (brown) exhibits two branches with a single turning point in the lower branch.[]{data-label="fig:cano_maJ"}](figure10.eps)

![Plots of $\beta = \mbox{const.}$ curves on the $(h,k)$ space for $\lambda \beta_0^2 = 9$ (brown), $\lambda \beta_0^2 = 12.7$ (green), $\lambda \beta_0^2 = 13.4$ (red) and $\lambda \beta_0^2 = 20$ (blue).[]{data-label="fig:hc_beta0"}](figure11.eps)

Phase Transitions {#sec:phase}
-----------------

We have shown in our earlier paper [@Czinner:2015eyk] that a Hawking–Page phase transition can be observed for static black holes in the Rényi approach. Previously, it had also been shown [@Caldarelli:1999xj; @Tsai:2011gv; @Altamirano:2014tva] 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 $\lambda J_0 =$ 0.01, 0.02 and 0.1 at constant $J$ is displayed on Fig. \[fig:FT\_Jconst\]. 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. \[sec:heat\_bath\]. It is generally assumed that $F\approx 0$ for a thermal gas, so the phase transition occurs around the temperature where the free energy of the black hole becomes zero.

![Free energy of a Kerr black hole in the Rényi model against the temperature for various angular momenta $J_0$, $\lambda J_0 = 0.01$ (green), $\lambda J_0 = 0.02$ (red) and $\lambda J_0 = 0.1$ (blue). Characterristic swallowtail behaviour is observed for $\lambda J_0 = 0.01$ (green), which corresponds to a SBH/LBH phase transition.[]{data-label="fig:FT_Jconst"}](figure12.eps)

An SBH/LBH phase transition can also be observed for Kerr black holes in the Rényi model when we enlarge the $\lambda J_0 =$ 0.01 curve of Fig. \[fig:FT\_Jconst\] on Fig \[fig:FT\_swallow\]. 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 $\lambda$, the swallowtail behavior disappears, as it can be seen on the $\lambda J_0 =$ 0.02 curve on Fig. \[fig:FT\_Jconst\]. This suggests the existence of a critical point where the phase transition becomes second order.

![Close up figure of the free energy of a Kerr black hole in the Rényi model for $\lambda J_0 = 0.01$ on Fig \[fig:FT\_Jconst\]. The intermediate, unstable branch is displayed with a dashed line.[]{data-label="fig:FT_swallow"}](figure13.eps)

Kerr-AdS black holes {#sec:Kerr-AdS}
====================

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 $$\begin{aligned}
 ds^2 & = & - \frac{\Delta_r}{\rho^2} \left( dt - \frac{a \sin^2 \theta}{\Xi} d \phi \right)^2 
 + \frac{\rho^2}{\Delta_r} dr^2 + \frac{\rho^2}{\Delta_\theta} d \theta^2 \nonumber \\
 & & + \frac{\Delta_\theta \sin^2 \theta}{\rho^2} \left( a dt - \frac{r^2 + a^2}{\Xi} d\phi \right)^2,\end{aligned}$$ where $$\Delta_r = (r^2 + a^2) \left( 1 + \frac{r^2}{l^2} \right) - 2 m r, ~~~ \Delta_\theta
= 1 - \frac{a^2 \cos^2 \theta}{l^2}, ~~~$$ $$\rho^2 = r^2 + a^2 \cos^2 \theta, ~~~ \Xi = 1 - \frac{a^2}{l^2}.$$ The thermodynamic quantities are written in terms of $a$, $l$ and the horizon radius $r_{+}$, which is obtained by solving $\Delta = 0$. The Hawking temperature of the horizon is given by $$T = \frac{1}{2 \pi r_{+}} \left( \frac{(a^2 +3 r_{+}^2)(r_{+}^2/l^2+1)}{2(a^2 + r_{+}^2)} -1\right),$$ while the Bekenstein-Hawking entropy of the black hole is $$S =\pi \frac{a^2 + r_{+}^2}{1 -a^2/l^2}.$$ The angular momentum of a Kerr-AdS black hole is $$J = \frac{(r_{+}^2 + a^2)(1 + r_{+}^2/l^2)}{2 r_{+}} \frac{a}{(1 - a^2/l^2)^2},$$ the angular velocity of the horizon is $$\Omega = \frac{a}{l^2} \frac{r_{+}^2 + l^2}{r_{+}^2 + a^2},$$ and the mass-energy parameter of the black hole can be re-expressed as $$M = \frac{(r_{+}^2 + a^2)(1 + r_{+}^2/l^2)}{2 r_{+}} \frac{1}{(1 - a^2/l^2)^2}.$$ The heat capacity at constant angular velocity can be computed as $$C_\Omega = \frac{2 \pi l^2 r_{+}^2( 3 r_{+}^4 +( a^2 + l^2) r_{+}^2 
 - a^2 l^2)}{(l^2 - a^2)(3 r_{+}^4 -(a^2 + l^2) r_{+}^2 - a^2 l^2)},$$ and the heat capacity at constant angular momentum takes the form $$C_J = \frac{2 \pi l^4  \left(a^2+r_{+}^2\right)^2 \left(-a^2 l^2+\left(a^2+l^2\right) r_{+}^2
 +3 r_{+}^4\right)}{(l^2 -a^2) X},$$ where $$\begin{aligned}
 X & = -l^4 r_{+}^4+3 l^2 r_{+}^6+a^6 \left(l^2+r_{+}^2\right) +a^4 \left(3 l^4+13 l^2 r_{+}^2+6 r_{+}^4\right) \\
   & \quad +a^2 \left(6 l^4 r_{+}^2+23 l^2 r_{+}^4+9 r_{+}^6\right).\end{aligned}$$ Here we introduced the normalized parameters $$p \equiv \frac{|a|}{r_{+}}, ~\mbox{and}~~ s \equiv \frac{l}{r_{+}}.$$ The heat capacities $C_\Omega$ 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_\Omega$ and $C_J$ as shown on Fig. \[fig:heat\_capacity\_AdS\].

![Phase diagram of Kerr-AdS black holes in the standard model. On the dashed curve $C_J$, while on the dotted curve $C_\Omega$ diverges. In region I, $C_{\Omega}<0$ and $C_{J} <0$, in region II, $C_{\Omega}<0$ and $C_{J} >0$, and in region III, $C_{\Omega}>0$ and $C_{J} >0$. In region IV, there is no physical solution because $|a| > l$.[]{data-label="fig:heat_capacity_AdS"}](figure14.eps)

Similarly to the Kerr-Rényi case in Sec. \[sec:Stability\], 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, $$\tilde{\beta} = \frac{\beta}{l},~~ 
 \tilde{J} = \frac{J}{l^2}, ~~
 \tilde{M} = \frac{M}{l}.$$

First we consider the microcanonical ensemble. The stability curves $\beta(M)$ at constant $J$ are plotted on Fig. \[fig:micro\_bM\_ads\], while the curves of constant $J$ in the parameter space of $(p,s)$ are plotted on Fig. \[fig:hc\_J\_ads\]. Just like in the Kerr-Rényi case, we can see that there is no turning point of stability. The stability curves $-\alpha(J)$ at constant $M$ and the curves of constant $M$ in the $(p,s)$ space are depicted on Fig. \[fig:micro\_maJ\_ads\] and Fig. \[fig:hc\_M0\_ads\] respectively. The behavior of the stability curves is almost identical to the one of the Kerr-Rényi case.

![Stability curves of $\beta(M)$ at fixed $J$ for Kerr-AdS black holes in the microcanonical treatment. The curves of $\tilde{J} = 0.01$ (green), $\tilde{J} = 0.025$ (red), and $\tilde{J} =0.05$ (blue) are plotted. No vertical tangent occurs in either case. The figure rotated by $\frac{\pi}{2}$ clockwise represents the stability curves of $-M(\beta)$ at fixed $J$ for the canonical ensemble, in which case the $\tilde{J} = 0.01$ (green) curve has two vertical tangents.[]{data-label="fig:micro_bM_ads"}](figure15.eps)

![Plots of $J = \mbox{const.}$ curves for $\tilde{J} = 0.01$ (green), $\tilde{J} = 0.025$ (red) and $\tilde{J} =0.05$ (blue) on the $(p,s)$ space.[]{data-label="fig:hc_J_ads"}](figure16.eps)

![Curves of the conjugate variable $-\alpha(J)$ at fixed $\beta$ for Kerr-AdS black holes in the canonical treatment. The stability curves of $\tilde{M} = 0.3$ (green), $\tilde{M} = 0.4$ (red) and $\tilde{M} = 0.6$ (blue) are plotted. There are no turning points on the diagram.[]{data-label="fig:micro_maJ_ads"}](figure17.eps)

![Plots of $M = \mbox{const.}$ curves for $\tilde{M} = 0.3$ (green), $\tilde{M} = 0.4$ (red) and $\tilde{M} = 0.6$ (blue) on the $(p,s)$ plane.[]{data-label="fig:hc_M0_ads"}](figure18.eps)

For the canonical system, the stability curves of $-M(\beta)$ at constant $J$ can be seen on Fig. \[fig:micro\_bM\_ads\] if we rotate it by $\frac{\pi}{2}$ clockwise. The figure shows the existence of a critical temperature, above which the Kerr-AdS black holes allow a first order SBH/LBH phase transition in the canonical ensemble. On Fig. \[fig:cano\_maJ\_ads\] we plotted the stability curves of $-\alpha(J)$ at constant $\beta$. The curves of constant $\beta$ on the $(p,s)$ space are plotted on Fig. \[fig:hc\_beta0\_ads\]. There are no turning points on the lower temperature (larger $\beta$) curves, but higher temperature curves exhibit turning points. Therefore, a stability change of Kerr-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 $(\tilde{J},-\alpha) = (0,0)$, where the Kerr-AdS black hole reduces to the Schwarzschild-AdS case. A vertical tangent occurs in the small black hole branch only, and no vertical tangent exists in the large black hole branch. From this result we can conclude that small and slowly rotating Kerr-AdS black holes are unstable in the canonical ensemble.

![Curves of the conjugate variable $-\alpha(J)$ at fixed $\beta$ for Kerr-AdS black holes in the canonical ensemble. The curves of $\tilde{\beta} = 3$ (brown), $\tilde{\beta} = 3.63$ (green), $\tilde{\beta} = 3.7$ (red) and $\tilde{\beta} = 4$ (blue) are plotted. The $\tilde{\beta} = 3.63$ (green) curve has two turning points, while the $\tilde{\beta} = 3$ (brown) curve has two branches and the lower branch has a turning point.[]{data-label="fig:cano_maJ_ads"}](figure19.eps)

![Plots of $\beta = \mbox{const.}$ curves on the $(p,s)$ space for $\tilde{\beta} = 3$ (brown), $\tilde{\beta} = 3.63$ (green), $\tilde{\beta} = 3.7$ (red) and $\tilde{\beta} = 4$ (blue).[]{data-label="fig:hc_beta0_ads"}](figure20.eps)

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 $\lambda$ and the AdS curvature parameter $l$ for black holes with identical horizon temperatures [@Biro:2013cra]. 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 equations $$\hat{T}_R (h,k) = \hat{T}_{AdS} (p,s), ~~~ \hat{\Omega}_R (h,k) = \hat{\Omega}_{AdS} (p,s),$$ where we normalized the quantities by the horizon radius $r_{+}$ as $$\hat{T} = T r_{+}, ~~~ \hat{\Omega} = \Omega r_{+}.$$ 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 $$\begin{aligned}
 p & = & \frac{3+h^2+k \pi -h^4 k \pi -Y}{2 h}, \\
 s & = &  \sqrt{\frac{3(1+h^2)}{3-h^2+2 k \pi -2 h^4 k \pi 
 -Y}},\end{aligned}$$ where $$Y = \sqrt{\left(3-h^2\right)^2+2 \left(3+h^2-3 h^4-h^6\right) k \pi 
 +\left(1-h^4\right)^2 k^2 \pi ^2}.$$ These equations provide a very interesting correspondence between the two approaches.

Summary and Conclusions {#sec:summary}
=======================

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 [@Biro:2013cra; @Czinner:2015eyk], 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 temperature-energy 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 $\lambda$, 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 $\lambda$ parameter value for the black hole problem is yet to be developed, but a similar approach has been considered in [@CzM] 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 $\lambda$ 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_{\lambda}(S)=\sqrt{S}$ and for Kerr black holes $H_{\lambda}(S)=\frac{S}{\sqrt{S-a^2\pi}}$ with $\lambda = 0$. The corresponding thermodynamic and stability problems (by also applying the formal logarithm method) has been studied in [@czinner2015black] and [@CzIg], 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 [@Biro:2012bka], 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 $\lambda=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 [@Biro:2013qea].

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. [@Carlip:2014pma] and references therein), and we expect that other parametric situations are also possible which might be connected to the parametric Rényi description.

The research leading to this result was supported by JSPS via an Invitation Fellowship for Research in Japan (Long-term) (No. L14710) and a Grant-in-Aid for Scientific Research (C) (No. 23540319). V.G.Cz thanks to Fundação para a Ciência e Tecnologia (FCT) Portugal, for financial support through Grant No. UID/FIS/00099/2013.

[99]{}

S.W. Hawking and D. Page, *Thermodynamics of Black Holes in anti-De Sitter Space*, *Commun. Math. Phys.* [**87**]{} (1983) 577. J.M. Maldacena, *The Large N limit of superconformal field theories and supergravity*, *Int. J. Theor. Phys.* [**38**]{} (1999) 1113. E. Witten, *Anti-de Sitter space and holography*, *Adv. Theor. Math. Phys.* [**2**]{} (1998) 253. D. Mateos, R.C. Myers and R.M. Thomson, *Holographic Phase Transitions with Fundamental Matter*, *Phys. Rev. Lett.* [**97**]{} (2006) 091601. D. Mateos, R.C. Myers and R.M. Thomson, *Thermodynamics of the brane*, *JHEP* [**05**]{} (2007) 067. J.D. Bekenstein, *Black holes and entropy*, *[Phys. Rev. D]{}* [**7**]{} (1973) 2333. J.M. Bardeen et al., *The Four laws of black hole mechanics*, *Commun. Math. Phys.* [**31**]{} (1973) 161. S.W. Hawking, *Particle Creation by Black Holes*, *Commun. Math. Phys.* [**43**]{} (1975) 199; S.W. Hawking, *Black Holes and Thermodynamics*, *[Phys. Rev. D]{}* [**13**]{} (1976) 191. T.S. Biró and P. Ván, *Zeroth law compatibility of nonadditive thermodynamics*, *[Phys. Rev. E]{}* [**83**]{} (2011) 061147. M.C. Mackey, *Time’s Arrow: The Origins of Thermodynamic Behaviour*, Springer, New York, (1992). C. Tsallis, *Introduction to Non-Extensive Statistical Mechanics: Approaching a Complex World*, Springer (2009). T.S. Biró, *Abstract composition rule for relativistic kinetic energy in the thermodynamical limit*, *Europhys. Lett.* [**84**]{} (2008) 56003. P.T. Landsberg, *Is equilibrium always an entropy maximum?*, *J. Stat. Phys.* [**35**]{} (1984) 159. A.M. Salzberg, *Exact Statistical Thermodynamics of Gravitational Interactions in One and Two Dimensions*, *J. Math. Phys.* [**6**]{} (1965) 158. M.E. Fisher and D. Ruelle, *The Stability of Many‐Particle Systems*, *J. Math. Phys.* [**7**]{} (1966) 260; L.G. Taff, *Celestia Mechanics*, Wiley, New York, (1985) p. 437. W.C. Saslaw, *Gravitational Physics of Stellar and Galactic Systems*, Cambridge University Press, Cambridge, (1985) p. 217. D. Pavon,*Thermodynamics of superstrings*, *Gen. Rel. Grav.* [**19**]{} (1987) 375. J. Binney and S. Tremaine, *Galactic Dynamics*, Princeton University Press, Princeton, (1987) p. 267. H.E. Kandrup, *Mixing and “violent relaxation” for the one-dimensional gravitational Coulomb gas*, *Phys. Rev. A* [**40**]{} (1989) 7265. H.S. Robertson, *Statistical Thermophysics*, Prentice-Hall, Englewood Cliffs, NJ, (1993) p. 96. H. Bacry, *The existence of dark matter in question*, *Phys. Lett. B* [**317**]{} (1993) 523. J.W. Gibbs, *Elementary Principles in Statistical Mechanics–Developed with Especial Reference to the Rational Foundation of Thermodynamics*, C. Scribner’s Sons, New York, (1902) (Yale University Press, New Haven, 1948; OX Bow Press, Woodbridge, Connecticut, 1981), p. 35. C. Tsallis and L.J.L. Cirto, *Black hole thermodynamical entropy*, *Eur. Phys. J. C* [**73**]{} (2013) 2487. P.C.W. Davies, *Proc. R. Soc. Lond. A*, *Thermodynamics of Black Holes* [**353**]{} (1977) 499. P.T. Landsberg and D. Tranah, *Entropies need not to be concave*, *Phys. Lett. A* [**78**]{} (1980) 219. N.T. Bishop and P.T. Landsberg, *The thermodynamics of a system containing two black holes and black-body radiation*, *Gen. Rel. Grav.* [**19**]{} (1987) 1083. D. Pavón and J.M. Rubí, *On some properties of the entropy of a system containing a black hole*, *Gen. Rel. Grav.* [**18**]{} (1986) 1245. J.  Maddox, *When entropy does not seem extensive*, *Nature*, [**365**]{} (1993) 103. G. Gour, *Entropy bounds for charged and rotating systems*, *Class. Quant. Grav.* [**20**]{} (2003) 3403. J. Oppenheim, *Thermodynamics with long-range interactions: From Ising models to black holes*, *[Phys. Rev. E]{}* [**68**]{} (2003) 016108. A. Pesci, *Entropy of gravitating systems: Scaling laws versus radial profiles*, *Class. Quant. Grav.* [**24**]{} (2007) 2283. R.F. Aranha et al., *The Efficiency of Gravitational Bremsstrahlung Production in the Collision of Two Schwarzschild Black Holes*, *Int. J. Mod. Phys. D* [**17**]{} (2008) 2049. P.T. Landsberg, *Thermodynamics and Statistical Mechanics*, Dover, New York (1990). C.E. Shannon, *A mathematical theory of communication*, *Bell. Syst. Tech. J.* [**27**]{} (1948) 379; ibid. 623. A.I. Khinchin, *Mathematical Foundations of Information Theory*, Dover, New York (1957). P. Tempesta, *Beyond the Shannon-Khinchin Formulation: The Composability Axiom and the Universal Group Entropy*, *Ann. Phys.* [**365**]{} (2016) 180. C. Tsallis, *Possible Generalization of Boltzmann-Gibbs Statistics*, *J. Stat. Phys.* [**52**]{} (1988) 479. T.S. Biró and V.G. Czinner, *A $q$-parameter bound for particle spectra based on black hole thermodynamics with Rényi entropy*, *Phys. Lett. B* [**726**]{} (2013) 861. A. Rényi, *On the dimension and entropy of probability distributions*, *Acta. Math. Acad. Sci. Hung.* [**10**]{} (1959) 193. A. Renyi, *Probability Theory*, North Holland, Amsterdam (1970). J.D. Bekenstein, *A Universal Upper Bound on the Entropy to Energy Ratio for Bounded Systems*, *[Phys. Rev. D]{}* [**23**]{} (1981) 287. C. Beck, *Generalized statistical mechanics of cosmic rays*, *Physica A* [**331**]{} (2004) 173. C. Beck, *Superstatistics in high-energy physics. Application to cosmic ray energy spectra and $e^+e^-$ annihilation*, *Eur. Phys. J. A*[**40**]{} (2009) 267. T.S. Biró and K. Ürmössy, *Non-extensive equilibration in relativistic matter*, *J. Phys. G* [**36**]{} (2009) 064044. O. Kaburaki, I. Okamoto and J. Katz, *Thermodynamic stability of Kerr black holes*, *[Phys. Rev. D]{}* [**47**]{} (1993) 2234. J. Katz, I. Okamoto, and O. Kaburaki, *Thermodynamic stability of pure black holes*, *Class. Quant. Grav.* [**10**]{} (1993) 1323. O. Kaburaki, *Critical behavior of extremal Kerr-Newman black holes*, *Gen. Rel. Grav.* [**28**]{} (1996) 843. H. Poincaré, *Sur l’[é]{}quilibre d’une masse fluide anim[é]{}e d’un mouvement de rotation*, *Acta. Math.* [**7**]{} (1885) 259. G. Arcioni and E. Lozano-Tellechea, *Stability and critical phenomena of black holes and black rings*, *[Phys. Rev. D]{}* [**72**]{} (2005) 104021. M. Azreg-Aïnou and M. E. Rodrigues, *Thermodynamical, geometrical and Poincaré methods for charged black holes in presence of quintessence*, *JHEP* [**09**]{} (2013) 146. V.G. Czinner and H. Iguchi, *Rényi entropy and the thermodynamic stability of black holes*, *Phys. Lett. B* [**752**]{} (2016) 306. V.G. Czinner and F.C. Mena, *Relative information entropy in cosmology: The problem of information entanglement*, *Phys. Lett. B* [**758**]{} (2016) 9. A. Chamblin et al., *Charged AdS black holes and catastrophic holography*, *[Phys. Rev. D]{}* [**60**]{} (1999) 064018. A. Chamblin et al., *Holography, thermodynamics and fluctuations of charged AdS black holes*, *[Phys. Rev. D]{}* [**60**]{} (1999) 104026. M.M. Caldarelli, G. Cognola and D. Klemm, *Thermodynamics of Kerr-Newman-AdS black holes and conformal field theories*, *Class. Quant. Grav.* [**17**]{} (2000) 399. Y-D. Tsai, X.N. Wu and Y. Yang, *Phase Structure of Kerr-AdS Black Hole*, *[Phys. Rev. D]{}* [**85**]{} (2012) 044005. N. Altamirano et al., *Thermodynamics of rotating black holes and black rings: phase transitions and thermodynamic volume*, *Galaxies* [**2(1)**]{} (2014) 89. S. Abe, *General pseudoadditivity of composable entropy prescribed by the existence of equilibrium*, *[Phys. Rev. E]{}* [**63**]{} (2001) 061105. http://tsallis.cat.cbpf.br/TEMUCO.pdf. C.A.M. Valdés et al., *Nonextensivity and Tsallis entropy in DNA fragmentation patterns by ionizing radiation*, *Journal of Modern Physics* [**3**]{} (2012) 431. I. Okamoto and O. Kaburaki, *The third law of thermodynamics for Kerr black holes*, *Mon. Not. Roy. Astron. Soc.* [**250**]{} (1991) 300. J. Katz, *On the number of unstable modes of an equilibrium*, *Mon. Not. Roy. Astron. Soc.* [**183**]{} (1978) 765. J. Katz, *On the Number of Unstable Modes of an Equilibrium - Part Two*, *Mon. Not. Roy. Astron. Soc.* [**189**]{} (1979) 817. R. Sorkin, *A Criterion for the onset of instability at a turning point*, *Astrophys. J.* [**249**]{} (1981) 254. O. Kaburaki, *Should entropy be concave?*, *Physics Letters A* [**185**]{} (1994) 21. V.G. Czinner, *Black hole entropy and the zeroth law of thermodynamics*, *Int. J. Mod. Phys. D* [**24**]{} (2015) 1542015. V.G. Czinner and Hideo Iguchi, *A zeroth law compatible model to Kerr black hole thermodynamics*, *Universe* [**3**]{} (2017) 14. T.S. Biró, *Ideal gas provides q-entropy*, *Physica A* [**392**]{} (2013) 3132. T.S. Biró et al., *Quark-gluon plasma connected to finite heat bath*, *Eur. Phys. J. A* [**49**]{} (2013) 110. S. Carlip, *Black Hole Thermodynamics*, *Int. J. Mod. Phys. D* [**23**]{} (2014) 1430023.
