Mass-induced instability of SAdS black hole in Einstein-Ricci cubic gravity

We perform the stability analysis of Schwarzschild-AdS (SAdS) black hole in the Einstein-Ricci cubic gravity. It shows that the Ricci tensor perturbations exhibit unstable modes for small black holes. We call this the mass-induced instability of SAdS black hole because the instability of small black holes arises from the massiveness in the linearized Einstein-Ricci cubic gravity, but not a feature of higher-order derivative theory giving ghost states. Also, we point out that the correlated stability conjecture holds for the SAdS black hole by computing the Wald entropy of SAdS black hole in Einstein-Ricci cubic gravity.


Introduction
The study of higher-derivative gravity theories has attracted critical attention in quantum gravity. Stelle's seminal work [1] has shown that the fourth-order gravity is renormalizable and has a finite Newtonian potential at origin. However, this gravity belongs to a nonunitary theory because it has a massive spin-2 pole with negative residue which could be interpreted as a state of negative norm (ghost). It turns out that the infinite derivative gravity (nonlocal gravity) is ghost-free and renormalizable around the Minkowski spacetime background when one chooses the exponential form of an entire function [2,3]. The most general class of theories that are ghost-free on any background is Lovelock gravity [4] whose terms of order k in the curvature are topological in d = 2k dimensions and vanish identically for d < 2k. Quasi-topological gravities provide additional example in higher dimensions than four [5,7,6], but all quasi-topological theories are trivial in four dimensions.
Recently, the Einsteinian cubic gravity of L EC = R − 2Λ 0 − λP/6 with P Riemann polynomials was introduced to indicate that it is neither topological nor trivial in four dimensions [8]. It was shown that black hole solutions of this gravity have a number of interesting properties [9,10,11], but these belong to either numerical or approximate solutions. That is, the Einstein equation cannot be solved analytically. An obstacle to studying these black holes is the lack of an analytic solution. At this stage, we remind the reader that at the critical points, the Einsteinian cubic gravity admits AdS black boles in four and five dimensions [12].
On the other hand, it is important to note that Ricci polynomials are much more manageable, compared to Riemann polynomials. The Ricci cubic gravity [13] can be composed of three terms from six cubic invariants in four dimensions, which has a similar property to the Ricci quadratic gravity (fourth-order gravity). It is known that the linearized theory of any higher-order gravity around a maximally symmetric background can be mapped into the linearized theory of fourth-order gravity [14]. This may imply that if one performs the linear stability analysis for a black hole obtained from a general quadratic gravity, these results could apply to analyzing the stability of the same black hole obtained from any higher-order gravity. A crucial benefit of Ricci cubic gravity is that the Schwarzschild black hole to Einstein gravity is a solution to this theory. In this case, the solution represented by mass r 0 = 2M describes the gravitational field outside of a static matter distribution, because its linearized theory reduces to that of Einstein gravity. It was proved that the only theories susceptible of admitting solutions with g tt g rr = −1 and representing the exterior field of a spherically symmetric distribution of mass are those that only propagate a massless spin-2 mode with 2 DOF (degrees of freedom) on the vacuum [15]. We note that the Schwarzschild-AdS (SAdS) black hole to Einstein gravity with a cosmological constant is a solution to the theory, but the SAdS is not a solution to the Einsteinian cubic gravity. Importantly, one can construct a covariant linearized gravity on the SAdS black hole in Ricci cubic gravity, but the covariant linearized theory of Riemann polynomials is allowed only on a maximally symmetric vacuum of AdS 4 spacetimes. However, the SAdS solution does not describe the exterior field of a spherically symmetric distribution of mass in Einstein-Ricci cubic gravity because its linearized equation (29) describing 5 DOF could not reduce to that of Einstein gravity.
Also, the (in)stability of SAdS black hole is known for fourth-order gravity [16] and thus, this result will be compared to that of Ricci cubic gravity. These are a few of reasons why we wish to introduce the Ricci cubic gravity as a higher-order gravity in the study of a black hole.
If a SAdS black hole is obtained from the Ricci cubic gravity, one has to ask what this all mean for "physical black hole"? A physical black hole can be selected by the stability analysis. If it is stable against the metric perturbation, one accepts it as a physical black hole. If not, one has to reject it. This is a long-standing issue since 1957. First of all, the linearized equation around the Schwarzschild black hole is given by δR µν (h) = 0 in Einstein gravity. Then, the metric perturbation h µν is classified depending on the transformation property under parity, namely odd and even. Using the Regge-Wheeler [17] and Zerilli gauge [18], one obtains two distinct perturbations: odd and even parities. It turns out that the Schwarzschild black hole is stable against the metric perturbation [19,20]. Investigating the stability analysis of the SAdS black hole in Einstein gravity with a cosmological constant, one might use the linearized Einstein equation δG µν (h) = 0. It turns out to be stable by following the Regge-Wheeler prescription [21,22].
We would like to mention that the Regge-Wheeler prescription is limited to the secondorder gravity. Thus, one could not implement the Regge-Wheeler prescription to perform the stability analysis of black hole found in higher-order gravity. For a higher-order gravity, Whitt [23] has argued that provided both massive spin-0 and spin-2 gravitons are nontachyonic, the Schwarzschild black hole is classically stable in fourth-order gravity when using the linearized Ricci tensor equation. In this case, the linearized Ricci tensor could represent a massive spin-2 field. Considering an auxiliary field formulation for decreasing a fourth-order gravity to a second-order theory of gravity [24], one found that the linearized equation for Ricci tensor δR µν is transformed exactly into that for auxiliary field ψ µν . Hence, one does not worry about the ghost (an unhealthy massive spin-2 field) problem arising from the fourth-order gravity because the linearized Ricci tensor δR µν as a healthy massive spin-2 field satisfies a second-order equation [25]. Visiting this stability issue again, it has shown that the small black hole in Einstein-Weyl gravity is unstable against s(l = 0)-mode Ricci tensor perturbation, while the large black hole is stable against s-mode perturbation [26].
Actually, this was performed by comparing the linearized Ricci tensor equation with the linearized metric equation around the five-dimensional black string where the Gregory-Laflamme instability appeared [27].
In addition, one observes that there was a close connection between thermodynamic instability and classical [Gregory-Laflamme] instability for the black strings/branes. This Gubser-Mitra proposal [28] was referred to as the correlated stability conjecture (CSC) [29].
The CSC states that the classical instability of a black string/brane with translational symmetry and infinite extent sets in precisely, when the corresponding thermodynamic system becomes thermodynamically unstable (that is, either Hessian matrix1 < 0 or heat capacity < 0). Here the additional assumption of translational symmetry and infinite extent has been added to ensure that finite size effects do not spoil the thermodynamic nature of the argument and to exclude a well-known case of the Schwarzschild black hole. A famous example of holding in the CSC is the five-dimensional black string. It is known that the Schwarzschild black hole is classically stable, but thermodynamically unstable because of its negative heat capacity. Also, the SAdS black hole is stable against the metric perturbation, whereas the small (large) black hole with r + < r * = ℓ/ √ 3(r + > r * ) is thermodynamically unstable (stable) because of negative (positive) heat capacity. Therefore, the last two examples show violation of the CSC and inapplicability of s-mode perturbation because a massless spin-2 mode perturbation starts from l = 2.
However, the situation is changed when one analyzes a black hole found in a fourthorder gravity (a massive gravity) where the infinite extent with translational symmetry is absent clearly. An important thing being different from the Einstein gravity is the appearance of a massive spin-2 mode. A massive spin-2 mode allows us to define a mass squared M 2 and to analyze the black hole with s(l = 0)-mode perturbation. Considering a setting e Ω r 0 t e i k r 0 z for a black string perturbation [27], there exists a critical wave number k c where for k < k c (k > k c ), the black sting is unstable (stable) against metric perturbations.
There is an unstable (stable) mode for any wavelength large (smaller) than the critical wavelength λ GL = 2πr 0 k c . Here, the mass M of a massive spin-2 mode plays a role of k/r 0 because the mass operator is given by ∂ z . This implies that the massiveness (M 2 = 0) takes over a black string located in z direction effectively. The Gregory-Laflamme instability is an s-wave unstable mode from the four-dimensional perspective [25]. In this respect, the dRGT massive gravity [30] having a Schwarzschild solution is subject directly to an s-wave instability [31,32]. Here, we pay our attention to an equivalence between 4D black hole in linearized massive gravity and 5D black string in linearized Einstein gravity from the four-dimensional perspective.
Furthermore, the Gregory-Laflamme instability condition (massiveness) picks up the small AdS black hole with r + < r * which is thermodynamically unstable in fourth-order gravity. It is known that the CSC holds for the SAdS black hole in Einstein-Weyl gravity by establishing a connection between the thermodynamic instability and the Gregory-Laflamme instability [33]. Also, the CSC holds for the BTZ black hole regardless of the horizon radius r + in three-dimensional new massive gravity. Hence it is quite interesting to check whether the CSC holds for a black hole found in higher-order gravity theory.
In this work, we will investigate classical instability and thermodynamics of SAdS black holes in Einstein-Ricci cubic gravity. Hereafter, we rename Ricci cubic gravity as Einstein-Ricci cubic gravity for a precise definition. We perform the stability analysis of SAdS black hole in Einstein-Ricci cubic gravity by introducing the Gregory-Laflamme scheme.
Computing the Wald entropy, we derive other thermodynamic quantities by making use of the first-law of thermodynamics in Einstein-Ricci cubic gravity. Finally, we wish to establish a connection between the Gregory-Laflamme instability and the thermodynamic instability of SAdS black holes in α = −3β Einstein-Ricci cubic gravity. This will provide another example for which the CSC holds.

Einstein-Ricci cubic gravity
We start with the Einstein-Ricci cubic gravity (ER) in four dimensional spacetimes [13] with κ = 1/G the inverse of Newtonian constant, Λ 0 the bare cosmological constant, and (e 1 , e 2 , e 3 ) three cubic parameters. Here we observe from the second term of (1) that the fourth-order gravity is embedded into the Einstein-Ricci cubic gravity. From the action (1), the Einstein equation is derived to be where the P -tensor is defined by Explicitly, it takes the form At this stage, we propose that a SAdS black hole solution to the Einstein gravity, is also a solution to Eq.(2). Here the metric function is given by with ℓ the curvature radius of AdS 4 spacetimes. The effective cosmological constant Λ is related to the bare cosmological constant as We note that a black hole mass parameter is determined as which is not surely the horizon radius r + . Hereafter we denote all background quantities with the "overbar". The background Ricci tensor and Ricci scalar are given bȳ In this case, one notes that the background P -tensor takes a maximally symmetric form It is easy to show that the SAdS black hole (5) to the Einstein equation of G µν − Λḡ µν = 0 is also the solution to the Einstein-Ricci cubic gravity when one substitutes (10) together with (9) into (2). However, it is important to note that the background Riemann tensor for the SAdS black hole is not given by the AdS 4 -curvature tensor which means that the SAdS spacetimes (5) is not a maximally symmetric vacuum.

Linearized Einstein-Ricci cubic gravity
To perform the stability analysis, we introduce the metric perturbation around the SAdS black hole as Then, we may define the linearized Ricci tensor and scalar as where with h = h ρ ρ . In this case, the linearized Einstein tensor can be written by Introducing two new parameters α = 4e 2 + 3e 3 and β = 2(6e 1 + e 2 ), the linearized Einstein equation can be rewritten compactly as where the background Lichnerowicz operators are defined by acting on scalar and tensor, respectively,∆ The linearized equation (17) is a second-order equation for δG µν and δR, but it becomes a fourth-order equation for h µν . This is surely our expectation that the linearized theory of any higher-order gravity is a fourth-order theory of gravity. Here, we have a fourth-order equation (17) because of the inclusion of the cosmological constant Λ 0 in the action (1).
Putting Λ = 0(Λ 0 = 0) yields Ricci-flat spacetimes on which cubic curvature terms give no contributions to the linearized equations, leading to δR µν = 0. This is one reason why we included the cosmological constant in the beginning action (1), compared to the fourthorder gravity. In addition, it is worth noting that Eq.(17) leads to (2.15) in Ref. [14] when replacingR µνρσ byR AdS 4 µνρσ in (11). Taking the trace of Eq.(17) leads to the linearized Ricci scalar equation One notes that Eq. (17) is a coupled second-order equation for δR µν and δR, which seems to be difficult to be solved. One way to avoid this difficulty is to split Eq. (17) into the traceless and trace parts by choosing α and β appropriately. For this purpose, we introduce a traceless Ricci tensor as Then, Eqs. (17) and (19) lead to where the mass squared µ 2 2 for spin-2 mode, and the mass squared µ 2 0 for spin-0 mode are given by First of all, decoupling of all massive modes requires either Λ = 0 or α = β = 0. The former case corresponds to the Ricci-flat spacetimes on which cubic curvature tensor gives no contribution to the linearized Ricci tensor equation. On the other hand, the latter case yields the quasi-topological Ricci cubic gravity whose linearized equation is exactly the same form as in the Einstein gravity with a cosmological constant [13]. Hence, these cases lead to a linearized second-order gravity.
A nontrivial decoupling between traceless and trace parts may occur when choosing a condition of α = −2β. This case corresponds to the linearized fourth-order (Einstein-Ricci quadratic) gravity [34] because one parameter e 3 can be represented by e 1 and e 2 . In this case, Eqs. (21) and (22) lead to the massive spin-2 and massive spin-0 equations, separately, Here the mass squared µ 2 2 and µ 2 0 are given by However, the number of degrees of freedom (DOF) for δR µν cannot be reduced to five because the contracted Bianchi identity of∇ µ δG µν = 0 does not imply the transverse due to δR = 0 in α = −2β Einstein-Ricci cubic gravity.
A promising choice may be done by requiring the non-propagation of the Ricci scalar.
From Eq. (19), imposing the condition of α = −3β, we have a constraint of non-propagating Ricci scalar Also, one finds from (23) that the mass squared µ 2 0 of massive spin-0 blows up, which means that the massive spin-0 is decoupled from the theory. This case corresponds to the linearized Einstein-Weyl gravity [16,33] because the linearized Ricci scalar is decoupled from the theory. Considering δG µν in (16) together with δR = 0, Eq.(17) leads to the massive spin-2 equation for δR µν as with the mass squared Here, one requires M 2 2 > 0 to avoid the tachyonic instability of δR µν propagating on the SAdS black hole [35]. Taking into account Eq.(28), the contracted Bianchi identity provides a desired transverse condition∇ µ δR µν = 0.
Hence, the DOF of δR µν becomes five from the counting of 10 − 4 − 1 = 5 in α = −3β Einstein-Ricci cubic gravity. Performing the stability analysis of the SAdS black hole in α = β = 0 Einstein-Ricci cubic gravity (quasi-topological Ricci cubic gravity) [13], one might use the linearized equation δG µν (h) = 0 with δG µν (h) in (16). It turns out to be stable by following the Regge-Wheeler prescription [21,22,36]. In these cases, the s(l = 0)-mode analysis is not necessary to show the stability of the Schwarzschild and SAdS black holes because the massless spin-2 mode starts from l = 2.

SAdS black hole in Einstein-Weyl gravity
The Regge-Wheeler prescription is no longer suitable for performing the stability analysis of the SAdS black hole in fourth-order gravity, since it focused on a linearized second-order gravity. One may explore the black hole stability by means of the second-order equation for the linearized Ricci tensor δR µν , instead of the fourth-order equation for the metric perturbation h µν [23,25]. The s-mode analysis is an essential tool to detect the instability of small SAdS black holes [26] obtained from Einstein-Weyl gravity given by with C µνρσ the Weyl tensor. Considering thermodynamics of the SAdS black hole in (5), we usually denote the small (large) SAdS black holes by the condition of r + < r * (r + > r * ) with r * = ℓ/ √ 3 and ℓ 2 = −3/Λ. Here r * is a position where the heat capacity blows up [see Fig. 4]. Its linearized equation is given by with the mass squared of massive spin-2 mode Here, one requires M 2 > 0 to avoid the tachyonic instability of δG µν propagating on the SAdS black hole background [35]. We depict M as a function of b with κ = 1 in Fig. 1. It is worth noting that M 2 is zero at b = b * = ℓ 2 /6 = 16.6, which corresponds to the critical gravity. The tachyon-free condition implies the allowed range for b as However, M 2 > 0 is not a necessary and sufficient condition to guarantee a stable SAdS black hole. If M 2 < 0(b > b * ), one does not need to consider a further analysis because it implies the tachyonic instability.
Even though one starts with 4 DOF, they are related to each other when one uses the Bianchi identity of∇ µ δG µν = 0 together with δG = −δR = 0. Hence, we derive one decoupled second-order equation for δG tr , where A, B and C were given by (20) in [37,33]. See O(1)/r 0 | r + =1,2,4 = 0.84, 0.40, 0.18 whose horizontal lines appear in Fig. 1. From the observation of Fig. 2 with O(1) ≃ 0.85, we find unstable modes for given r + = 1, 2, 4 as As the horizon size r + increases, the instability region becomes narrow and narrow. We call this instability as the Gregory-Laflamme instability [27,37] because the four-dimensional linearized equation for h µν around the five-dimensional black string background leads to (33) when replacing h µν by δG µν . We check that for r + = 6 > r * = 5.7, the maximum value of Ω is less than 10 −4 , which implies that there is no unstable modes for large black hole with r + > r * . The case of M = 0 yields the critical gravity avoiding massive spin-2 and spin-0 modes when choosing the transverse-traceless gauge for h µν [35].
It show clearly that the Gregory-Laflamme instability of small black holes in the Einstein-Weyl gravity is due to the massiveness of M ∈ (0, M c ), but not a feature of fourth-order gravity giving ghost states. Taking into account the number of degrees of freedom (DOF), it is helpful to show that the SAdS black hole is physically stable in the Einstein gravity [21,22], whereas the small SAdS black hole is unstable in the Einstein-Weyl gravity. The number of DOF of the metric perturbation is two in the Einstein gravity, while the number of DOF of massive spin-2 δG µν is five in the Einstein-Weyl gravity. The s(l = 0)-mode analysis of the massive spin-2 with five DOF shows the Gregory-Laflamme instability. The s-mode analysis is useful for handling the massive spin-2 mode in the Einstein-Weyl gravity, but is not suitable for the massless spin-2 mode in the Einstein gravity.

SAdS black hole in Einstein-Ricci cubic gravity
First of all, we consider the α = −2β Einstein-Ricci cubic gravity. We wish to solve Eq.(24) for the traceless Ricci tensor δR µν . Actually, it is observed that Eq.(24) becomes Eq.(33) when substituting δR µν and µ 2 2 by δG µν and M 2 . However, we note that δG µν is a transverse-traceless tensor in Einstein-Weyl gravity, while δR µν is a traceless tensor in the α = −2β Einstein-Ricci cubic gravity. Hence, the DOF of δG µν are five, whereas the DOF of δR µν is nine. The non-transversality for δR µν does not make a further progress on the s-mode analysis of stability. Now, let us consider the α = −3β Einstein-Ricci cubic gravity. Its linearized equation is given by Eq. (29) with the mass squared (30). Eq.(29) becomes Eq.(33) exactly when substituting δR µν and M 2 2 by δG µν and M 2 . Importantly, the traceless Ricci tensor δR µν satisfies the transverse condition (31). Hence, its DOF is determined to be five as for δG µν in Einstein-Weyl gravity. At this stage, the even-parity metric perturbation could be chosen for a single s-mode analysis in the α = −3β Einstein-Ricci cubic gravity and whose form is given by Once one starts with 4 DOF, they are related to each other when using the Bianchi identity where A, B and C were given by (20) in [33]. A physical mode of δR tr grows exponentially in time as e Ωt with Ω > 0, spatially vanishes at the AdS infinity, and regular at the future horizon [37]. Here M 2 is the mass of massive spin-2 mode given by with κ = 1. We depict M 2 as a function of α in Fig. 3. The massive spin-2 mass M 2 is zero at α = α * = ℓ 4 /3 = 3333, where the critical gravity appears. We require M 2 2 > 0 to avoid the tachyonic instability of δR µν propagating on the SAdS black hole background [35]. The tachyon-free condition implies the allowed range for α as If M 2 2 < 0(α > α * ), we do not need to perform a further analysis for the stability because it indicates the tachyonic instability.
It is emphasized that M 2 2 > 0 is not a necessary and sufficient condition to obtain a stable SAdS black hole. We need to follow Gregory-Laflamme scheme to distinguish between stable and unstable black holes by solving Eq.(41) numerically. Observing Fig. 2 when replacing M by M 2 , we note that M 2 ends up at M c 2 ≡ O(1)/r 0 | r + =1,2,4 = 0.84, 0.40, 0.18 with O(1) = 0.85 whose horizontal lines appear in Fig. 3. Here, we find unstable modes for small black holes with r + = 1, 2, 4 as As the horizon size r + increases, the instability region starting from the origin becomes narrow and narrow. For given r + = 1, 2, 4, three horizontal lines M c 2 which are ending points split unstable (M 2 < M c 2 ) and stable (M 2 > M c 2 ) black holes. It indicates that the Gregory-Laflamme instability of small black holes in the α = −3β Einstein-Ricci cubic gravity is due to the massiveness (0 < M 2 < M c 2 ) of massive spin-2 mode, but not a feature of higher-order gravity giving a ghost. This ghost may appear only when expressing the linearized equation (29) in terms of the metric perturbation h µν , giving seven DOF. Here the massive ghost (an unhealthy massive spin-2 mode) does not appear because we used the linearized Ricci tensor δR µν to represent a healthy massive spin-2 mode.
However, this does not mean that our perturbation analysis misses the ghost instability.
The ghost is present in the spectrum of the theory because the theory provides a fourthorder linearized equation when expressing in terms of h µν . The matter is how to represent a massive spin-2 mode propagating on the SAdS black hole spacetimes. Expressing the linearized equation in term of the linearized Ricci tensor instead of the metric perturbation, it becomes a second-order linearized equation. This is one tip to handle a fourth-order linearized gravity. Considering an auxiliary field formulation for decreasing a fourth-order gravity to a second-order theory of gravity [24], one found that the linearized equation for Ricci tensor δR µν is transformed exactly into that for auxiliary field ψ µν . Using an auxiliary field enables one to distinguish the perturbation related to a massive ghost. Furthermore, we wish to point out that the linearized Ricci tensor could represent a massive spin-2 mode with five DOF [38,39]. We stress here that the Gregory-Laflamme instability based on δR µν has nothing to do with the ghost instability and it reflects a feature of massive gravity described by the linearized equation (29).

Wald entropy and black hole thermodynamics
It is known that the correlated stability conjecture proposed by Gubser-Mitra [28] does not hold for the SAdS black hole found in Einstein gravity, but it holds for the SAdS black hole found in Einstein-Weyl gravity [33]. In order to confirm the classical instability found in the previous section, one has to explore the thermodynamic property of the SAdS black hole obtained from the Einstein-Ricci cubic gravity.

Einstein-Weyl gravity
We start with the Einstein-Weyl gravity because the thermodynamic quantities of the SAdS black hole was completely computed by employing the Abbot-Deser-Tekin method [40]. It was well-known that the Wald entropy of Einstein-Weyl gravity (32) is given by with S BH = πr 2 + the Bekenstein-Hawking entropy of the SAdS black hole in the Einstein gravity. The other thermodynamic quantities of mass, heat capacity, and free energy are given by [40,33] where the thermodynamic quantities including the Hawking temperature for a SAdS black hole in Einstein gravity take the forms We check that the first-law of thermodynamics is satisfied in Einstein-Weyl gravity as as the first-law is indeed satisfied in Einstein gravity where "d" denotes the differentiation with respect to the horizon radius r + only. capacity. This picture persists to C EW with b < ℓ 2 /6 and C α=−3β with α < ℓ 4 /3.
We briefly sketch the thermodynamic stability of the SAdS black hole found in Einstein-Weyl gravity [33]. First we consider the case of b < ℓ 2 /6(M 2 > 0) which is dominantly described by the Einstein-Hilbert term. Since the heat capacity C SAdS blows up at r + = r * = ℓ/ √ 3 = 5.7 [see Fig. 4], we divide the black hole into the small black hole with r + < r * and the large black hole with r + > r * . We know that the small black hole is thermodynamically unstable because the heat capacity C EW < 0, while the large black hole is thermodynamically stable because C EW > 0. For the other case of b > ℓ 2 /6(M 2 < 0) which is dominantly described by the Weyl term, the situation reverses. The small black hole is thermodynamically stable because C EW > 0, while the large black hole is thermodynamically unstable because C EW < 0. In this case, the mass squared of massive spin-2 is negative, which implies the tachyonic instability. So, this case should be excluded from the consideration. We would like to mention that there is no connection between classical stability and thermodynamic instability for small SAdS black hole in Einstein gravity. However, let us see how things are improved in Einstein-Weyl gravity. For a small black hole with r + < r * and b/ℓ 2 /6, the heat capacity is negative which means that it is thermodynamically unstable. On the other hand, we observe from (38) that a small black hole with r + < r * is unstable against the s-mode massive spin-2 perturbation.
The Gregory-Laflamme instability condition picks up the small SAdS black hole which is thermodynamically unstable in Einstein-Weyl gravity. This implies that the CSC [28] holds for the SAdS black hole found in Einstein-Weyl gravity.

Einstein-Ricci cubic gravity
First of all, we wish to compute the Wald entropy of the SAdS black hole in Einstein-Ricci cubic gravity. The Wald entropy is defined by the following integral performed on 2-dimensional spacelike bifurcation surface Σ [41,42,43]: where dV 2 2 = r 2 sin θdθdφ is the volume element on Σ and ǫ µν is the binormal vector to Σ normalized as ǫ µν ǫ µν = −2. We note that the superscript (0) denotes that the functional derivative with respect to R µνρσ is evaluated on-shell. The background (on-shell) P -tensor P µνρσ is given by (10). Now, the Wald entropy takes the form with A = 4πr 2 + and κ = 1/G = 1. In the case of α = −3β Einstein-Ricci cubic gravity, the Wald entropy takes the form Up to now, we know the Wald entropy S α=−3β W and the Hawking temperature However, we do not know the other thermodynamic quantities of the SAdS black hole in Einstein-Ricci cubic gravity. We propose that the first-law should be satisfied in α = −3β Einstein-Ricci cubic gravity as Using (54) together with (52) and (53), we derive the mass of the SAdS black hole Other thermodynamic quantities of heat capacity and free energy are computed to be Now we are in a position to mention the thermodynamic stability of the SAdS black hole in α = −3β Einstein-Ricci cubic gravity. First, we consider the case of α < ℓ 4 /3 < 1(M 2 2 > 0) which is dominantly described by the Einstein-Hilbert term. Since the heat capacity C α=−3β blows up at r + = r * = ℓ/ √ 3 = 0.57ℓ [see Fig. 4], we divide still the black hole into the small black hole with r + < r * and the large black hole with r + > r * . Then, it suggests that the small black hole is thermodynamically unstable because C α=−3β < 0, while the large black hole is thermodynamically stable because C α=−3β > 0. For the other case of α > ℓ 4 /3(M 2 2 < 0) which is dominantly described by the Ricci cubic terms, the thermodynamic stability reverses. The small black hole is thermodynamically stable because C α=−3β > 0, while the large black hole is thermodynamically unstable because C α=−3β < 0. However, we note that the mass squared M 2 2 of massive spin-2 mode is negative. This corresponds to the tachyonic instability and thus, this case is unacceptable.
One finds from (44) that for M 2 2 > 0, a small (large) black hole with r + < r * (r + > r * ) is unstable (stable) against the s-mode massive spin-2 perturbation δR µν . The Gregory-Laflamme instability condition picks up the small SAdS black hole which is thermodynamically unstable in α = −3β Einstein-Ricci cubic gravity. This indicates that the CSC proposed by Gubser-Mitra [28] holds for the SAdS black hole found in α = −3β Einstein-Ricci cubic gravity too. However, the other case of M 2 2 < 0 corresponds to the tachyonic instability and its thermodynamic stability contradicts to the conventional one. Hence, the CSC does not hold for M 2 2 < 0.
i) α = β = 0 case. This case describes a quasi-topological Ricci cubic gravity because its linearized equation around the SAdS black hole reduces to that of Einstein gravity with cosmological constant. There is no counterpart in fourth-order gravity. The absence of massive spin-0 and spin-2 modes implies that the fourth-order terms give no contribution to this linearized theory as if they were purely topological, implying that this linearized theory is ghost-free. The SAdS black hole is stable regardless of the horizon radius r + in α = β = 0 Einstein-Ricci cubic gravity because its linearized equation takes the form of δG µν (h) = 0 without any mass terms. However, the small black hole with r + < r * is thermodynamically unstable while the large black hole with r + > r * is thermodynamically unstable. The CSC does not hold for the α = β = 0 Einstein-Ricci cubic gravity which is a ghost-free theory.
ii) α = −2β case. This case allows a decoupling of massive spin-0 mode from massive spin-0 mode. This corresponds to the fourth-order gravity. Here δR µν is only a traceless tensor whose DOF is nine. We note that the absence of transversality condition∇ µ δR µν = 0 makes the stability analysis of SAdS black hole difficult. Its Wald entropy is given by All thermodynamic quantities are similar to those for α = −3β Einstein-Ricci cubic gravity except replacing [1 − αΛ 2 /3] by [1 − αΛ 2 ]. For α < ℓ 4 /9, the small black hole with r + < r * is thermodynamically unstable while the large black hole with r + > r * is thermodynamically unstable. In this case, we may not discuss the CSC because we could not explore the stability issue of the SAdS black hole explicitly in α = −2β Einstein-Ricci cubic gravity.
iii) α = −3β case. This is related to the Einstein-Weyl gravity. We show by adopting the Gregory-Laflamme scheme that the small SAdS black hole with r + < r * is classically unstable, while the large SAdS black hole with r + > r * is classically stable. It indicates clearly that the Gregory-Laflamme instability arises from the massiveness (0 < M 2 < M c 2 ) of massive spin-2 mode, but not from a feature of Einstein-Ricci cubic (fourth-order) gravity giving ghost states. On the other hand, its thermodynamic quantities are computed by making use of the Wald entropy and the first-law of thermodynamics. For the case of M 2 2 > 0 which is dominantly described by the Einstein-Hilbert term, one finds that the small black hole with r + < r * is thermodynamically unstable, whereas the large black hole with r + > r * is thermodynamically stable. This shows that the CSC proposed by Gubser-Mitra [28] works for the SAdS black hole obtained from the α = −3β Einstein-Ricci cubic gravity. The other case of M 2 2 < 0 corresponds to the tachyonic instability and its thermodynamic stability is unacceptable. Hence, the CSC does not hold for M 2 2 < 0. iv) α = −3β case with M 2 2 = 0(α = α * ). In this case, considering the transverse and traceless gauge of∇ µ h µν = 0 and h = 0, the linearized Einstein tensor takes the from of δG µν = −(∆ L − 2Λ)h µν /2. Then, its linearized equation is given by (∆ L − 2Λ) 2 h µν = 0 that corresponds to the critical gravity. It turned out be an unstable theory, even though massive spin-2 and massive spin-0 modes are decoupled from the theory [35]. On the other hand, all thermodynamic quantities disappear except the Hawking temperature. Therefore, we could not discuss the CSC for this case.
Consequently, we have obtained the SAdS black hole solution (5) from the Einstein-Ricci cubic gravity (1). A physical black hole could be found by performing the analysis of stability which means a perturbation analysis based on physical field around the black hole because the black hole is a solution in the curved spacetimes [20,44]. In this work, we have shown that the α = −3β Einstein-Ricci cubic gravity passes the stability analysis of SAdS black hole. The α = β = 0 case provides a quasi-topological Ricci cubic gravity whose linearized theory is just the linearized Einstein gravity, implying that the black hole solution (5) becomes a physical black hole in α = β = 0 Einstein-Ricci cubic gravity.
However, we could not establish the (in) stability for other relations between α and β. In addition, it is worth noting that the α = −3β Einstein-Ricci cubic gravity has a similar property to the Einstein-Weyl gravity because any higher-order gravity could be mapped into the quadratic gravity at the linearized level.
Finally, if one considers the Riemann cubic gravity (Einsteinian cubic gravity), covariant linearized gravity is possible only in maximally symmetric vacua. Also, black hole solutions are numerical or approximate solutions in the Einsteinian cubic gravity and an obstacle to studying these black hole is the lack of an analytic solution [45].