###### Abstract

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.

Mass-induced instability of SAdS black hole

in Einstein-Ricci cubic gravity

[10mm]

Yun Soo Myung^{*}^{*}*e-mail address: ysmyung@inje.ac.kr

[8mm]

Institute of Basic Sciences and Department of Computer Simulation, Inje University Gimhae 50834, Korea

[0pt]

## 1 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 (non-local 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 in the curvature are topological in dimensions and vanish identically for . 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 with 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 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 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 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 in Einstein gravity. Then, the metric perturbation 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 . 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 second-order 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 non-tachyonic, 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 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 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 -mode Ricci tensor perturbation, while the large black hole is stable against -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 or heat capacity ). 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 is thermodynamically unstable (stable) because of negative (positive) heat capacity. Therefore, the last two examples show violation of the CSC and inapplicability of -mode perturbation because a massless spin-2 mode perturbation starts from .

However, the situation is changed when one analyzes a black hole found in a fourth-order 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 and to analyze the black hole with -mode perturbation. Considering a setting for a black string perturbation [27], there exists a critical wave number where for , the black sting is unstable (stable) against metric perturbations. There is an unstable (stable) mode for any wavelength large (smaller) than the critical wavelength . Here, the mass of a massive spin-2 mode plays a role of because the mass operator is given by . This implies that the massiveness () takes over a black string located in direction effectively. The Gregory-Laflamme instability is an -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 -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 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 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 Einstein-Ricci cubic gravity. This will provide another example for which the CSC holds.

## 2 Einstein-Ricci cubic gravity

We start with the Einstein-Ricci cubic gravity (ER) in four dimensional spacetimes [13]

(1) | |||||

with the inverse of Newtonian constant, the bare cosmological constant, and 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

(2) |

where the -tensor is defined by

(3) |

Explicitly, it takes the form

(4) | |||||

At this stage, we propose that a SAdS black hole solution to the Einstein gravity,

(5) |

is also a solution to Eq.(2). Here the metric function is given by

(6) |

with the curvature radius of AdS spacetimes. The effective cosmological constant is related to the bare cosmological constant as

(7) |

We note that a black hole mass parameter is determined as

(8) |

which is not surely the horizon radius . Hereafter we denote all background quantities with the “overbar”. The background Ricci tensor and Ricci scalar are given by

(9) |

In this case, one notes that the background -tensor takes a maximally symmetric form

(10) |

It is easy to show that the SAdS black hole (5) to the Einstein equation of 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-curvature tensor

(11) |

which means that the SAdS spacetimes (5) is not a maximally symmetric vacuum.

## 3 Linearized Einstein-Ricci cubic gravity

To perform the stability analysis, we introduce the metric perturbation around the SAdS black hole as

(12) |

Then, we may define the linearized Ricci tensor and scalar as

(13) |

where

(14) | |||||

(15) |

with . In this case, the linearized Einstein tensor can be written by

(16) |

Introducing two new parameters and , the linearized Einstein equation can be rewritten compactly as

(17) |

where the background Lichnerowicz operators are defined by acting on scalar and tensor, respectively,

(18) |

The linearized equation (17) is a second-order equation for and , but it becomes a fourth-order equation for . 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 in the action (1). Putting yields Ricci-flat spacetimes on which cubic curvature terms give no contributions to the linearized equations, leading to . This is one reason why we included the cosmological constant in the beginning action (1), compared to the fourth-order gravity. In addition, it is worth noting that Eq.(17) leads to (2.15) in Ref.[14] when replacing by in (11).

Taking the trace of Eq.(17) leads to the linearized Ricci scalar equation

(19) |

One notes that Eq.(17) is a coupled second-order equation for and , 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

(20) |

Then, Eqs.(17) and (19) lead to

(21) | |||

(22) |

where the mass squared for spin-2 mode, and the mass squared for spin-0 mode are given by

(23) |

First of all, decoupling of all massive modes requires either or . 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 . This case corresponds to the linearized fourth-order (Einstein-Ricci quadratic) gravity [34] because one parameter can be represented by and . In this case, Eqs.(21) and (22) lead to the massive spin-2 and massive spin-0 equations, separately,

(24) | |||

(25) |

Here the mass squared and are given by

(26) |

However, the number of degrees of freedom (DOF) for cannot be reduced to five because the contracted Bianchi identity of does not imply the transverse condition,

(27) |

due to in 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 , we have a constraint of non-propagating Ricci scalar

(28) |

Also, one finds from (23) that the mass squared 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 in (16) together with , Eq.(17) leads to the massive spin-2 equation for as

(29) |

with the mass squared

(30) |

Here, one requires to avoid the tachyonic instability of propagating on the SAdS black hole [35]. Taking into account Eq.(28), the contracted Bianchi identity provides a desired transverse condition

(31) |

Hence, the DOF of becomes five from the counting of in Einstein-Ricci cubic gravity.

## 4 Black hole stability

### 4.1 SAdS black hole in Einstein gravity

First of all, in Einstein-Ricci cubic gravity (Einstein gravity), the linearized equation around the Schwarzschild black hole is given by with in (14). Then, the metric perturbation is classified depending on the transformation properties under parity, namely odd and even. Using the Regge-Wheeler [17] and Zerilli gauge [18], one obtains two distinct perturbations: odd with 2 DOF and even with 4 DOF. This implies that even though one starts with 6 DOF under the Regge-Wheleer gauge, the propagating DOF is two for a massless spin-2 metric tensor . It turned out that the Schwarzschild black hole is stable against the metric perturbation [19, 20].

Performing the stability analysis of the SAdS black hole in Einstein-Ricci cubic gravity (quasi-topological Ricci cubic gravity) [13], one might use the linearized equation with in (16). It turns out to be stable by following the Regge-Wheeler prescription [21, 22, 36]. In these cases, the -mode analysis is not necessary to show the stability of the Schwarzschild and SAdS black holes because the massless spin-2 mode starts from .

### 4.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 , instead of the fourth-order equation for the metric perturbation [23, 25]. The -mode analysis is an essential tool to detect the instability of small SAdS black holes [26] obtained from Einstein-Weyl gravity given by

(32) |

with 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 with and . Here is a position where the heat capacity blows up [see Fig. 4]. Its linearized equation is given by

(33) |

with the mass squared of massive spin-2 mode

(34) |

Here, one requires to avoid the tachyonic instability of propagating on the SAdS black hole background [35]. We depict as a function of with in Fig. 1. It is worth noting that is zero at , which corresponds to the critical gravity. The tachyon-free condition implies the allowed range for as

(35) |

However, is not a necessary and sufficient condition to guarantee a stable SAdS black hole. If , one does not need to consider a further analysis because it implies the tachyonic instability.

Here, one has which means that there is no massive spin-0 (Ricci scalar) propagation in Einstein-Weyl gravity. Therefore, the linearized Einstein tensor satisfies the transverse-traceless condition of and . Hence, its DOF is five from counting of . The even-parity metric perturbation in Einstein gravity is used for a single -mode analysis in the Einstein-Weyl gravity. Its form is given by and as displayed in the matrix form

(36) |

Even though one starts with 4 DOF, they are related to each other when one uses the Bianchi identity of together with .

Hence, we derive one decoupled second-order equation for ,

(37) |

where and were given by (20) in [37, 33]. See Fig. 2 that is obtained by solving Eq.(37) numerically. We note that the small black holes of horizon radii correspond to with , respectively. ends up at whose horizontal lines appear in Fig. 1. From the observation of Fig. 2 with , we find unstable modes for given as

(38) |

As the horizon size 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 around the five-dimensional black string background leads to (33) when replacing by . We check that for , the maximum value of is less than , which implies that there is no unstable modes for large black hole with . The case of yields the critical gravity avoiding massive spin-2 and spin-0 modes when choosing the transverse-traceless gauge for [35]. At the critical point, the massless spin-2 modes have zero energy whereas the massive spin-2 modes are replaced by the log modes. The presence of log modes implies another instability of the SAdS black hole in critical gravity.

For given , from Figs. 1 and 2, the stable condition of the SAdS black hole in Einstein-Weyl gravity is given by

(39) |

It show clearly that the Gregory-Laflamme instability of small black holes in the Einstein-Weyl gravity is due to the massiveness of , 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 is five in the Einstein-Weyl gravity. The -mode analysis of the massive spin-2 with five DOF shows the Gregory-Laflamme instability. The -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.

### 4.3 SAdS black hole in Einstein-Ricci cubic gravity

First of all, we consider the Einstein-Ricci cubic gravity. We wish to solve Eq.(24) for the traceless Ricci tensor . Actually, it is observed that Eq.(24) becomes Eq.(33) when substituting and by and . However, we note that is a transverse-traceless tensor in Einstein-Weyl gravity, while is a traceless tensor in the Einstein-Ricci cubic gravity. Hence, the DOF of are five, whereas the DOF of is nine. The non-transversality for does not make a further progress on the -mode analysis of stability.

Now, let us consider the 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 and by and . Importantly, the traceless Ricci tensor satisfies the transverse condition (31). Hence, its DOF is determined to be five as for in Einstein-Weyl gravity. At this stage, the even-parity metric perturbation could be chosen for a single -mode analysis in the Einstein-Ricci cubic gravity and whose form is given by

(40) |

Once one starts with 4 DOF, they are related to each other when using the Bianchi identity of together with . Thus, we derive one decoupled second-order equation for ,

(41) |

where and were given by (20) in [33]. A physical mode of grows exponentially in time as with , spatially vanishes at the AdS infinity, and regular at the future horizon [37]. Here is the mass of massive spin-2 mode given by

(42) |

with . We depict as a function of in Fig. 3. The massive spin-2 mass is zero at , where the critical gravity appears. We require to avoid the tachyonic instability of propagating on the SAdS black hole background [35]. The tachyon-free condition implies the allowed range for as

(43) |

If , we do not need to perform a further analysis for the stability because it indicates the tachyonic instability.

It is emphasized that 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 by , we note that ends up at with whose horizontal lines appear in Fig. 3. Here, we find unstable modes for small black holes with as

(44) |

As the horizon size increases, the instability region starting from the origin becomes narrow and narrow. For given , three horizontal lines which are ending points split unstable () and stable () black holes.

It indicates that the Gregory-Laflamme instability of small black holes in the Einstein-Ricci cubic gravity is due to the massiveness () 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 , giving seven DOF. Here the massive ghost (an unhealthy massive spin-2 mode) does not appear because we used the linearized Ricci tensor 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 fourth-order linearized equation when expressing in terms of . 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 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 has nothing to do with the ghost instability and it reflects a feature of massive gravity described by the linearized equation (29).

## 5 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.

### 5.1 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

(45) |

with 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]

(46) |

where the thermodynamic quantities including the Hawking temperature for a SAdS black hole in Einstein gravity take the forms

(47) |

We check that the first-law of thermodynamics is satisfied in Einstein-Weyl gravity as

(48) |

as the first-law is indeed satisfied in Einstein gravity

(49) |

where denotes the differentiation with respect to the horizon radius only.

We briefly sketch the thermodynamic stability of the SAdS black hole found in Einstein-Weyl gravity [33]. First we consider the case of which is dominantly described by the Einstein-Hilbert term. Since the heat capacity blows up at [see Fig. 4], we divide the black hole into the small black hole with and the large black hole with . We know that the small black hole is thermodynamically unstable because the heat capacity , while the large black hole is thermodynamically stable because . For the other case of which is dominantly described by the Weyl term, the situation reverses. The small black hole is thermodynamically stable because , while the large black hole is thermodynamically unstable because . 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 and , 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 is unstable against the -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.

### 5.2 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]:

(50) |

where is the volume element on and is the binormal vector to normalized as . We note that the superscript (0) denotes that the functional derivative with respect to is evaluated on-shell. The background (on-shell) -tensor is given by (10). Now, the Wald entropy takes the form

with and . In the case of Einstein-Ricci cubic gravity, the Wald entropy takes the form

(52) |

Up to now, we know the Wald entropy and the Hawking temperature

(53) |

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 Einstein-Ricci cubic gravity as

(54) |

Using (54) together with (52) and (53), we derive the mass of the SAdS black hole

(55) |

Other thermodynamic quantities of heat capacity and free energy are computed to be

(56) | |||||

(57) |

Now we are in a position to mention the thermodynamic stability of the SAdS black hole in Einstein-Ricci cubic gravity. First, we consider the case of which is dominantly described by the Einstein-Hilbert term. Since the heat capacity blows up at [see Fig. 4], we divide still the black hole into the small black hole with and the large black hole with . Then, it suggests that the small black hole is thermodynamically unstable because , while the large black hole is thermodynamically stable because . For the other case of which is dominantly described by the Ricci cubic terms, the thermodynamic stability reverses. The small black hole is thermodynamically stable because , while the large black hole is thermodynamically unstable because . However, we note that the mass squared 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 , a small (large) black hole with is unstable (stable) against the -mode massive spin-2 perturbation . The Gregory-Laflamme instability condition picks up the small SAdS black hole which is thermodynamically unstable in Einstein-Ricci cubic gravity. This indicates that the CSC proposed by Gubser-Mitra [28] holds for the SAdS black hole found in Einstein-Ricci cubic gravity too. However, the other case of corresponds to the tachyonic instability and its thermodynamic stability contradicts to the conventional one. Hence, the CSC does not hold for .

## 6 Discussions

We would like to mention the Einstein-Ricci cubic gravity according to the relation between and .

i) 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 in Einstein-Ricci cubic gravity because its linearized equation takes the form of without any mass terms. However, the small black hole with is thermodynamically unstable while the large black hole with is thermodynamically unstable. The CSC does not hold for the Einstein-Ricci cubic gravity which is a ghost-free theory.

ii) case.
This case allows a decoupling of massive spin-0 mode from massive spin-0 mode.
This corresponds to the fourth-order gravity. Here is only a traceless tensor whose DOF is nine. We note that the absence of transversality condition makes the stability analysis of SAdS black hole difficult. Its Wald entropy is given by . All thermodynamic quantities are similar to those for Einstein-Ricci cubic gravity except replacing by . For , the small black hole with is thermodynamically unstable while the large black hole with 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 Einstein-Ricci cubic gravity.

iii) case.
This is related to the Einstein-Weyl gravity. We show by adopting the Gregory-Laflamme scheme that the small SAdS black hole with is classically unstable, while the large SAdS black hole with is classically stable. It indicates clearly that the Gregory-Laflamme instability arises from the massiveness () 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 which is dominantly described by the Einstein-Hilbert term, one finds that the small black hole with is thermodynamically unstable, whereas the large black hole with is thermodynamically stable. This shows that the CSC proposed by Gubser-Mitra [28] works for the SAdS black hole obtained from the Einstein-Ricci cubic gravity. The other case of corresponds to the tachyonic instability and its thermodynamic stability is unacceptable. Hence, the CSC does not hold for .

iv) case with .
In this case, considering the transverse and traceless gauge of and , the linearized Einstein tensor takes the from of . Then, its linearized equation is given by 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 Einstein-Ricci cubic gravity passes the stability analysis of SAdS black hole. The 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 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 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].

Acknowledgments

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MOE) (No. NRF-2017R1A2B4002057).

## References

- [1] K. S. Stelle, Phys. Rev. D 16, 953 (1977). doi:10.1103/PhysRevD.16.953
- [2] L. Modesto, Phys. Rev. D 86, 044005 (2012) doi:10.1103/PhysRevD.86.044005 [arXiv:1107.2403 [hep-th]].
- [3] T. Biswas, E. Gerwick, T. Koivisto and A. Mazumdar, Phys. Rev. Lett. 108, 031101 (2012) doi:10.1103/PhysRevLett.108.031101 [arXiv:1110.5249 [gr-qc]].
- [4] D. Lovelock, J. Math. Phys. 12, 498 (1971). doi:10.1063/1.1665613
- [5] R. C. Myers and B. Robinson, JHEP 1008, 067 (2010) doi:10.1007/JHEP08(2010)067 [arXiv:1003.5357 [gr-qc]].
- [6] A. Cisterna, L. Guajardo, M. Hassaine and J. Oliva, JHEP 1704, 066 (2017) doi:10.1007/JHEP04(2017)066 [arXiv:1702.04676 [hep-th]].
- [7] J. Oliva and S. Ray, Class. Quant. Grav. 27, 225002 (2010) doi:10.1088/0264-9381/27/22/225002 [arXiv:1003.4773 [gr-qc]].
- [8] P. Bueno and P. A. Cano, Phys. Rev. D 94, no. 10, 104005 (2016) doi:10.1103/PhysRevD.94.104005 [arXiv:1607.06463 [hep-th]].
- [9] R. A. Hennigar, D. Kubiznak and R. B. Mann, Phys. Rev. D 95, no. 10, 104042 (2017) doi:10.1103/PhysRevD.95.104042 [arXiv:1703.01631 [hep-th]].
- [10] J. Ahmed, R. A. Hennigar, R. B. Mann and M. Mir, JHEP 1705, 134 (2017) doi:10.1007/JHEP05(2017)134 [arXiv:1703.11007 [hep-th]].
- [11] P. Bueno and P. A. Cano, Phys. Rev. D 94, no. 12, 124051 (2016) doi:10.1103/PhysRevD.94.124051 [arXiv:1610.08019 [hep-th]].
- [12] X. H. Feng, H. Huang, Z. F. Mai and H. Lu, Phys. Rev. D 96, no. 10, 104034 (2017) doi:10.1103/PhysRevD.96.104034 [arXiv:1707.06308 [hep-th]].
- [13] Y. Z. Li, H. S. Liu and H. Lu, arXiv:1708.07198 [hep-th].
- [14] P. Bueno, P. A. Cano, V. S. Min and M. R. Visser, Phys. Rev. D 95, no. 4, 044010 (2017) doi:10.1103/PhysRevD.95.044010 [arXiv:1610.08519 [hep-th]].
- [15] P. Bueno and P. A. Cano, Class. Quant. Grav. 34, no. 17, 175008 (2017) doi:10.1088/1361-6382/aa8056 [arXiv:1703.04625 [hep-th]].
- [16] Y. S. Myung, Phys. Lett. B 728, 422 (2014) doi:10.1016/j.physletb.2013.12.019 [arXiv:1308.1455 [gr-qc]].
- [17] T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957). doi:10.1103/PhysRev.108.1063
- [18] F. J. Zerilli, Phys. Rev. Lett. 24, 737 (1970). doi:10.1103/PhysRevLett.24.737
- [19] C. V. Vishveshwara, Phys. Rev. D 1, 2870 (1970). doi:10.1103/PhysRevD.1.2870
- [20] S. Chandrasekhar, in The Mathematical Theory of Black Holes (Oxford University, New York, 1983).
- [21] V. Cardoso and J. P. S. Lemos, Phys. Rev. D 64, 084017 (2001) doi:10.1103/PhysRevD.64.084017 [gr-qc/0105103].
- [22] A. Ishibashi and H. Kodama, Prog. Theor. Phys. 110, 901 (2003) doi:10.1143/PTP.110.901 [hep-th/0305185].
- [23] B. Whitt, Phys. Rev. D 32, 379 (1985). doi:10.1103/PhysRevD.32.379
- [24] S. Mauro, R. Balbinot, A. Fabbri and I. L. Shapiro, Eur. Phys. J. Plus 130, no. 7, 135 (2015) doi:10.1140/epjp/i2015-15135-0 [arXiv:1504.06756 [gr-qc]].
- [25] K. S. Stelle, Int. J. Mod. Phys. A 32, no. 09, 1741012 (2017). doi:10.1142/S0217751X17410123
- [26] Y. S. Myung, Phys. Rev. D 88, no. 2, 024039 (2013) doi:10.1103/PhysRevD.88.024039 [arXiv:1306.3725 [gr-qc]].
- [27] R. Gregory and R. Laflamme, Phys. Rev. Lett. 70, 2837 (1993) doi:10.1103/PhysRevLett.70.2837 [hep-th/9301052].
- [28] S. S. Gubser and I. Mitra, hep-th/0009126.
- [29] T. Harmark, V. Niarchos and N. A. Obers, Class. Quant. Grav. 24, R1 (2007) doi:10.1088/0264-9381/24/8/R01 [hep-th/0701022].
- [30] C. de Rham, G. Gabadadze and A. J. Tolley, Phys. Rev. Lett. 106, 231101 (2011) doi:10.1103/PhysRevLett.106.231101 [arXiv:1011.1232 [hep-th]].
- [31] E. Babichev and A. Fabbri, Class. Quant. Grav. 30, 152001 (2013) doi:10.1088/0264-9381/30/15/152001 [arXiv:1304.5992 [gr-qc]].
- [32] R. Brito, V. Cardoso and P. Pani, Phys. Rev. D 88, no. 2, 023514 (2013) doi:10.1103/PhysRevD.88.023514 [arXiv:1304.6725 [gr-qc]].
- [33] Y. S. Myung and T. Moon, JHEP 1404, 058 (2014) doi:10.1007/JHEP04(2014)058 [arXiv:1311.6985 [hep-th]].
- [34] H. Lu, A. Perkins, C. N. Pope and K. S. Stelle, Phys. Rev. D 96, no. 4, 046006 (2017) doi:10.1103/PhysRevD.96.046006 [arXiv:1704.05493 [hep-th]].
- [35] H. Liu, H. Lu and M. Luo, Int. J. Mod. Phys. D 21, 1250020 (2012) doi:10.1142/S0218271812500204 [arXiv:1104.2623 [hep-th]].
- [36] T. Moon, Y. S. Myung and E. J. Son, Eur. Phys. J. C 71, 1777 (2011) doi:10.1140/epjc/s10052-011-1777-0 [arXiv:1104.1908 [gr-qc]].
- [37] T. Hirayama and G. Kang, Phys. Rev. D 64, 064010 (2001) doi:10.1103/PhysRevD.64.064010 [hep-th/0104213].
- [38] D. M. Eardley, D. L. Lee and A. P. Lightman, Phys. Rev. D 8, 3308 (1973). doi:10.1103/PhysRevD.8.3308
- [39] Y. S. Myung and T. Moon, JCAP 1410, no. 10, 043 (2014) doi:10.1088/1475-7516/2014/10/043 [arXiv:1403.5433 [gr-qc]].
- [40] H. Lu and C. N. Pope, Phys. Rev. Lett. 106, 181302 (2011) doi:10.1103/PhysRevLett.106.181302 [arXiv:1101.1971 [hep-th]].
- [41] R. M. Wald, Phys. Rev. D 48, no. 8, R3427 (1993) doi:10.1103/PhysRevD.48.R3427 [gr-qc/9307038].
- [42] V. Iyer and R. M. Wald, Phys. Rev. D 50, 846 (1994) doi:10.1103/PhysRevD.50.846 [gr-qc/9403028].
- [43] T. Jacobson, G. Kang and R. C. Myers, Phys. Rev. D 49, 6587 (1994) doi:10.1103/PhysRevD.49.6587 [gr-qc/9312023].
- [44] O. J. Kwon, Y. D. Kim, Y. S. Myung, B. H. Cho and Y. J. Park, Phys. Rev. D 34, 333 (1986). doi:10.1103/PhysRevD.34.333
- [45] R. A. Hennigar, M. B. J. Poshteh and R. B. Mann, arXiv:1801.03223 [gr-qc].