BTZ black holes in massive gravity

BTZ black holes in massive gravity BTZ black holes in massive gravity Sumeet Chougule• , Sanjib Dey∗ , Behnam Pourhassan◦ and Mir Faizal†§ • School of Physics, University of Hyderabad, Hyderabad, India Department of Physics, Indian Institute of Science Education and Research Mohali, Sector 81, SAS Nagar, Manauli 140306, India ◦ School of Physics, Damghan University, Damghan, 3671641167, Iran † Irving K. Barber School of Arts and Sciences, University of British Columbia -Okanagan, Kelowna, British Columbia V1V 1V7, Canada § Department of Physics and Astronomy, University of Lethbridge, Lethbridge, Alberta, T1K 3M4, Canada E-mail: sumeetchougule@live.com, dey@iisermohali.ac.in, b.pourhassan@umz.ac.ir, mirfaizalmir@googlemail.com arXiv:1809.00868v1 [gr-qc] 4 Sep 2018 ∗ Abstract We analyze certain aspects of BTZ black holes in massive theory of gravity. The black hole solution is obtained by using the Vainshtein and dRGT mechanism, which is asymptotically AdS with an electric charge. We study the Hawking radiation using the tunneling formalism as well as analyze the black hole chemistry for such system. Subsequently, we use the thermodynamic pressure-volume diagram to explore the efficiency of the Carnot heat engine for this system. Some of the important features arising from our solution include the non-existence of quantum effects, critical Van der Walls behaviour, thermal fluctuations and instabilities. Moreover, our solution violates the Reverse Isoperimetric Inequality and, thus, the black hole is super-entropic, perhaps which turns out to be the most interesting characteristics of the BTZ black hole in massive gravity. 1 Introduction Astronomical observations suggest that our universe is expanding, in fact, the observational data from the type-I supernovae proposes that our universe is expanding in an accelerating rate [1–3]. Theoretically, this accelerated expansion is due to the creation of a negative pressure implying a positive vacuum energy density which is caused due to a cosmological constant term in the Einstein’s field equation. Such a cosmological constant can be associated with the zero point energy of quantum field theories, however, the value of the cosmological constant obtained from the quantum field theoretical calculations is several orders of magnitude greater than that arising from the observational astronomy. Various models have been proposed to explain the origin of the cosmological constant [4–6] including some large distance modifications of general relativity [7]. All of these modifications are mainly constrained in such a way that they are compatible with the theory of general relativity at a scale at which general relativity has been observed [8,9]. Alternative schemes also exist, for instance, it is possible to obtain a long distance infrared modification of the general relativity along with massive gravitons, where the mass of the graviton can be obtained from the observational data [10]. It may be noted that by adding a small Fierz-Pauli like mass term to the action of general relativity, we do not obtain a stable theory consistent with the zero mass limit [11]. In fact, such type of straightforward modified theories are not physical as they violate the known experimental bounds obtained from solar system tests [8,9]. The inclusion of non-linearities in some theories, where the Vainshtein mechanism helps to recover General Relativity at the solar system scales [12, 13], gives rise to the Boulware-Deser ghost [14], which can be removed by introducing higher-order terms in the massive action as it was done in the case of the dRGT theory of Massive 1 BTZ black holes in massive gravity gravity [15], where a subclass of massive potentials is considered [16, 17]. Within this framework, the higher order-term when grouped with the ghost like term becomes a total derivative and, thus, the action is resumed having two free parameters. Mass terms can be produced by using many other mechanisms, for instance, by breaking the Lorentz symmetry of the system [18]. Nevertheless, mass terms are very useful and they have been utilized to study various interesting models; such as, in Gauss-Bonnet massive gravity [19]. It has been observed that the massive gravitons can lead to interesting modification of black hole thermodynamics. The modification to the behavior of black hole by the inclusion of graviton mass has also been analyzed in the extended phase space [10] in order to study the phase transition of black holes [20]. Besides, the cosmological solutions [21] and the initial value constraint [22], etc., have been explored in the context of massive gravity. BTZ black hole is another interesting object which was introduced [23]. It is possible to construct a BTZ black hole in massive gravity. In fact, an asymptotically AdS charged BTZ black hole has been constructed in a massive theory of gravity and various different aspects of such a solution have been studied [24]. In this work, we study both the dynamical and thermodynamical aspects of such a solution. In fact, using the geometrical thermodynamic approaches Weinhold, Ruppeiner and Quevedo metrics have been studied for this system. We analyze the Hawking radiation [25] for such a system using the tunneling formalism [26–29]. We also analyze the black hole chemistry for this BTZ black hole [30, 31]. This is done by relating each black hole parameter with a chemical equivalent using the first law of thermodynamics. The cosmological constant is considered a thermodynamic parameter related to the pressure of the system [30, 31]. Subsequently, using the pressure-volume terms, we study the heat engines and their efficiency [32, 33]. 2 Massive gravity Let us start by recollecting the notions of three dimensional massive gravity [15]. First, we define the constants √ for massive gravity ci and the symmetric polynomials of the eigenvalues Ui of the d × d µ matrix Kν = g µα fαν U1 = [K], U2 = [K]2 − [K2 ], U3 = [K]3 − 3[K][K2 ] + 2[K3 ], U4 = [K]4 − 6[K2 ][K]2 + 8[K3 ][K] + 3[K2 ]2 − 6[K4 ], (2.1) √ √ with [K] = Kµµ and ( A)µν ( A)νλ = Aµλ . The three dimensional action of massive gravity with an abelian U (1) gauge field and negative cosmological constant is known to have the form [24] " # Z 4 X √ 1 S=− ci Ui (g, f ) , (2.2) d3 x −g R − 2Λ + L(F) + M̃ 2 16π i=1 where L(F) is the Lagrangian for the vector gauge field, Λ stands for the cosmological constant, R represents the scalar curvature and M̃ , f are the mass term and fixed symmetric tensor, respectively. Also, F = Fµν F µν is the Maxwell invariant, with Fµν = ∂µ Aν − ∂ν Aµ being the Faraday tensor and Aµ being the gauge potential. By using the action (2.2), one can utilize the variational principle to obtain the field equations for gravitation, as computed in [24]. Nevertheless, in order to obtain a static solution of the charged AdS black hole in 3D, we can start with the following ansatz for the metric [24, 34] dr2 + r2 dφ2 , (2.3) ds2 = −f (r)dt2 + f (r) where f (r) is an arbitrary function of the radial coordinate. An exact solution of the metric (2.3) can be obtained by choosing a reference metric as given by fµν = diag(0, 0, c2 hij ), 2 (2.4) BTZ black holes in massive gravity with c being a positive constant. With the given ansatz (2.4), Ui ’s can easily be computed as U1 = c/r, U2 = U3 = U4 = 0, which indicates that the contribution of massive gravity is arising only from the U1 . Furthermore, keeping in mind that we are going to study a linearly charged BTZ black hole, we can choose the Lagrangian of Maxwell field as L(F) = −F. In addition, by considering a gauge potential related to the radial electric field to be of the form Aµ = h(r)δµt , and by following the procedure explained in [24], one obtains an exact form of the radial function f (r) as given by [24] r f (r) = −Λr2 − m − 2q 2 ln ( ) + M̃ 2 cc1 r, (2.5) l from which one obtains the exact expression of the metric for the massive gravity in the given scenario. Here m = 8M and q = 2Q, with M, Q being the mass and electric charge of the black hole, respectively. Here, l is an arbitrary constant having the dimension of the length, which is arising from the fact that the logarithmic arguments should be dimensionless. In what follows, we shall consider the cosmological constant Λ = −1/l2 , since Λ has a dimension of inverse squared length. However, it should be noted that the metric corresponding to (2.5) can also be constructed by using other methods available in the literature. For instance, in [35–38], the authors have explored a procedure by using the Stückelberg method, where the Stückelberg fields can be considered to be in a unitary gauge, so that the corresponding fiducial metric becomes the Minkowskian. This is the simplest case that one can consider. However, the fiducial metric coming out of such theories may not be unique, as it depends on the choice of the gauge field. If the Stückelberg fields are not in unitary gauge then one obtains an associated fiducial metric also but it is not a Minkowskian anymore. This method is quite simple and, surely, it has its own beauty, however, in this paper, we have constructed the metric from a slightly different procedure as explained earlier in this section, where we have used a reference metric for the purpose. 3 Hawking Radiation as Tunneling Among many approaches of analyzing the Hawking radiation of a black hole, in this manuscript we consider the method of quantum tunneling [26–29]. There are many reasons for this. For instance, other approaches for deriving Hawking radiation deal with the principle of detailed balance, the background geometry is considered fixed and the energy conservation is not enforced during the emission process in general cases. More precisely, in a general cases like massive gravity, the energy conservation is not valid in the usual sense, however, this is because the time-like Killing vector in massive gravity is not defined in the same direction of the ordinary time-coordinate. Actually, the fact is that the black-hole radiation is related to the way how one defines the time (vacuum), therefore, in order to restore the notion of energy conservation for those cases, one can redefine the time, as indicated in [35]. There are many ways to redefine such time coordinate, such as, the path integral method [36, 37], Bogolibov method [38], etc. However, in tunneling formalism the energy conservation is utilized to obtain non-thermal corrections to the spectrum of particles and, thus, this method shows the conservation of energy in a more explicit way. Moreover, as the tunneling process takes place at the horizon, the coordinate system is required to be non-singular at the horizon. Thus, Painlevé coordinates are useful for such analysis [39]. In this formalism it is argued that when a classical stable system becomes quantum mechanically unstable, it is natural to consider tunneling. The Hawking radiation occurs due to the tunneling of virtual particles. The idea is to consider the vacuum fluctuations near the horizon which creates a pair of particle and anti-particle. When a pair is created just inside the horizon, the positive energy particle tunnels across the horizon and escapes to infinity as a Hawking radiation. While the black hole absorbs the negative energy particle and its mass is decreased. Similarly, for the pair which is created outside the horizon, the anti-particle tunnels inside the black hole before it is annihilated. Thus, in both of the cases the black hole absorbs the negative energy particle by decreasing its mass, while the positive energy particle escapes to infinity to be observed as Hawking radiation. Under this formalism, the 3 BTZ black holes in massive gravity probability of tunneling is given as [26] Γ ∼ e−2ImS , (3.1) where S is the action of the trajectory. The barrier for the tunneling is provided by the outgoing particle itself. Black holes lose energy due to the radiation and, thus, it shrinks in order to conserve the energy. Consequently, the horizon is contracted with respect to its original size and the amount of contraction depends on the energy of the outgoing particle. In this way, the outgoing particle itself provides the barrier. Now, for the case of a massive BTZ black hole the form of the metric is given by (2.3), which reduces to a form with [24] f (r) = r2 − m − 2q 2 ln(r) − M̃ 2 cc1 r ≃ r2 − m − 2q 2 (r − 1) − M̃ 2 cc1 r, for a constant value of l = 1. The horizon for this metric is at q 2q 2 + M̃ 2 cc1 ± (2q 2 + M̃ 2 cc1 )2 − 4(2q 2 − m) r± = , 2 (3.2) (3.3) so that we can write the metric (2.3) in the following modified form ds2 = −(r2 − m − 2q 2 (r − 1) − M̃ 2 cc1 r)dt2 + (r2 − m − 2q 2 (r − 1) − M̃ 2 cc1 r)−1 dr2 + r2 dφ2 . (3.4) Notice that, at r+ there is a coordinate singularity, therefore, in order to study the physics across the horizon we need to change the coordinate system again such that the metric is well behaved at the horizon. Therefore, we use the Painlevé time t [39], which defines a new time coordinate with respect to the Schwarzschild time ts with an arbitrary function R̃(r) t = ts − R̃(r), dt = dts − R̃′ (r)dr, dt2 = dt2s + R̃′2 (r)dr2 − 2R̃′ (r)drdts , with which we can rewrite (3.4) as follows h i ds2 = −f (r)dt2s + f −1 (r) − f (r)R̃′2 (r) dr2 + 2f (r)R̃′ (r)drdts + r2 dφ2 . (3.5) (3.6) Since, R̃(r) has been considered as an arbitrary function, we have the freedom p to specify it in such a way that the coefficient of dr2 in (3.6) becomes unity and, thus, R̃′ (r) = 1 − f (r)/f (r), so that (3.6) further reduces to p (3.7) ds2 = −f (r)dt2s + 2 1 − f (r)drdts + dr2 + r2 dφ2 . Correspondingly, the radial null geodesic is given by p 0 = −f (r)dt2s + 2 1 − f (r)drdts + dr2   p dr dr 2 0 = −f (r) + 2 1 − f (r) + , dts dts (3.8) so that ṙ = ±1 − p 1 − f (r), (3.9) where the upper (lower) sign corresponds to the outgoing (ingoing) geodesics with the assumption that the time increases towards the future. Let us now consider the pair production inside the horizon at rin ≃ r+ . If, ω be the energy of the particle created, the mass of the black hole after the emission of the particle becomes m − ω and, hence, the horizon contracts from rin = 2q 2 + M̃ 2 cc1 + q q    (2q 2 + M̃ 2 cc1 )2 − 4(2q 2 − m) /2 to rout = 2q 2 + M̃ 2 cc1 + (2q 2 + M̃ 2 cc1 )2 − 4(2q 2 − m + ω) /2. 4 BTZ black holes in massive gravity The difference between rout and rin acts as a barrier of potential V for the particle tunneling. In this region, ω < V and, therefore, the action is imaginary which can be written as follows Z rout Z rout Z pr Z rout Z m−ω dH ′ ImS = Im dr, (3.10) pr dr = Im dpr dr = Im ṙ rin 0 rin rin m where we use the Hamilton’s equation to replace dp′r by dH/ṙ, followed by a change of variable from momentum to energy. Subsequently, by considering the case of outgoing geodesic in (3.9) we obtain Z rout Z ω −dω ′ p ImS = Im dr (3.11) 1 − f (r) 0 1− rin Z rout Z ω −dω ′ q = Im dr, (3.12) 0 1− rin 1 − r2 + m + 2q 2 (r − 1) + M̃ 2 cc1 r − ω ′ where H = m − ω ′ . Note that while the self-gravitation of the system is taken into account, the mass of the black hole decreases from m to m − ω ′ and, thus, we replace m by m − ω ′ in (3.12). Now, considering u = 1 − r2 + m + 2q 2 (r − 1) + M̃ 2 cc1 r − ω ′ , we have du = −dω ′ , therefore, we can write Z rout Z u(ω) du √ dr, ImS = Im (3.13) u(0) (1 − u) rin which has a simple pole at u = 1 and, thus, the residue at u = 1 is −2. Therefore, (3.13) becomes Z rout ImS = −Im 4πdr = −4π(rout − rin ). (3.14) rin Correspondingly, the transmission probability is given by Γ(ω) ≃ e−2ImS = e8π(rout −rin ) = e8πσ , (3.15) with 8ω 2 −4ω +h σ = rout − rin = q i3/2 , 2 2 2 2 2 2 2 2 (2q + M̃ cc1 ) − 4(2q − m) (2q + M̃ cc1 ) − 4(2q − m) (3.16) where we have considered the binomial series upto second order in ω. The second term in the exponential in (3.15) is the non-thermal correction to Hawking radiation, whereas the first order term corresponds to the Boltzmann factor exp[− Tω ], such that the Hawking temperature TH turns out to be q (2q 2 + M̃ 2 cc1 )2 − 4(2q 2 − m) 2r+ − 2q 2 − M̃ 2 cc1 = , (3.17) TH = 32π 32π where the effect of the massive parameter is clearly visible. 4 Black Hole Chemistry It is customary that every black hole parameter is associated with a chemical equivalent compatible with the first law of thermodynamics [30, 31]. Therefore, we can write dE = T dS + V dP + work terms, κ dM = dA + ΩdJ + Φdq, 8π 5 (4.1) (4.2) BTZ black holes in massive gravity and compare the mass M with the internal energy E, surface gravity κ with the temperature T and the horizon area A with entropy S, however, we do not have any gravitational analogue for pressure P and volume V in space-time with Λ = 0. But, for space-time with non-zero cosmological constant it is possible to find an analogue to pressure-volume terms. The basic idea of black hole chemistry is to regard Λ as a thermodynamical variable in analogy to the pressure in the first law. The mass M is then considered to be the gravitational analogue of chemical enthalpy, which we denote by M ′ . Under this framework, the pressure P is related to cosmological constant Λ as P =− (D − 2)(D − 1) Λ = , 8π 16πl2 (4.3) where D is the dimension of the system. The most general Smarr formula [40] for D < 4 for a charged singly-rotating black hole is given by [30, 31] (D − 3)GD M = (D − 2)T S + (D − 2)ΩJ − 2V P + (D − 3)Φq, (4.4) where J is the angular momentum, Ω represents angular velocity and GD stands for the Ddimensional Newton’s constant. For the case of charged non-massive black holes, (4.4) and the first law of thermodynamics dM ′ = T dS + V dP + Φdq, (4.5) hold [30], however, our motivation is to test whether both of them are satisfied for the charged massive BTZ black hole with mass term M̃ . Since, in this case, the metric is given by (2.3), we obtain the temperature as follows [30] T = f ′ (r+ ) M̃ 2 cc1 r+ q2 = 2− − . 2π πl πr+ 2π (4.6) We also compute the entropy S, pressure P and enthalpy M ′ as given in the following 1 1 S = πr+ , P = , 2 8πl2   r2 M̃ 2 cc1 r+ q2 r+ − ln . M ′ = +2 − 4l 2 l 4 By using the above equations (4.7) and (4.8), we can rewrite the enthalpy as   q2 32P S 2 8S 2 P ′ − ln − 4l2 M̃ 2 cc1 SP, M (S, P, q) = π 4 π (4.7) (4.8) (4.9) so that the volume V and the electric potential Φ turn out to be V = Φ = dM ′ dP dM ′ dq (S,q) (S,P ) 2 − 2q 2 πl2 − 2πr+ l2 M̃ 2 cc1 , = 2πr+ = −q ln   r+ . l (4.10) (4.11) Note that, the volume defined by (4.10) is a thermodynamical volume and is not the usual geometrical volume. Nevertheless, it is easy to cross check that all of our results satisfy the first law of thermodynamics (4.5), which ensure the fact that all of our calculations are indeed correct. However, the Smarr relation (4.4) is not satisfied in the given case as expected and it is indicated already in some articles, for instance in [24, 30], that for massive black holes the usual Smarr relation may be violated and one may need to modify it accordingly. In order to preserve the Smarr relation let us introduce an extra parameter G in the Enthalpy M ′ (4.8) as follows     2 r+ q2 M̃ 2 cc1 r+ 8S 2 P q2 4l2 M̃ 2 cc1 SP r+ 32P S 2 ′ ln = − ln , (4.12) M = 2− − − 4l 2 l 4G π 4 π G 6 BTZ black holes in massive gravity so that the electric potential remains the same as given by (4.11), however, the volume is modified as dM ′ 2πr+ l2 M̃ 2 cc1 M̃ 2 cc1 S 8S 2 q2 2 V = = 2πr+ − 2q 2 πl2 − = − − . (4.13) dP (S,q,G) G π 4P 2πP G Because of the insertion of G, we are forced to introduce another parameter K K= dM ′ dG = (S,P,q) M̃ 2 cc1 r+ , 4G 2 (4.14) which is the thermodynamic conjugate to G. In order to conserve the first law of thermodynamics, (4.12) enforces us to introduce the parameter G into the metric (2.3) also, so that (2.5) is modified to r M̃ 2 cc1 r f (r) = −Λr2 − m − 2q 2 ln ( ) + , (4.15) l G and, thus, the temperature obtains a new form T = r+ q2 M̃ 2 cc1 f ′ (r+ ) = 2− − . 2π πl πr+ 2πG (4.16) With these new formalism, we propose a generalization of the Smarr formula (4.4) (D − 3)GD M = (D − 2)T S + (D − 2)ΩJ − 2V P − 2KG + (D − 3)Φq, (4.17) which can be applied to the BTZ black holes including the massive case in D < 4. It is straightforward to verify that the generalized Smarr relation (4.17) is satisfied with the expressions given in (4.7), (4.12), (4.13), (4.14) and (4.16). Surely, the new results satisfy the first law of thermodynamics (4.5), which ensures that there is no violation of the basic principle even after introducing the new parameter G. Let us now set out to calculate the Gibbs free energy and its derivative with respect to r+   2 r+ q2 q2 r+ ′ ln (4.18) + , G = M − TS = − 2 − 4l 2 l 2 r 2 + q 2 l2 dG =− + < 0. (4.19) dr+ 2r+ l2 Note that the Gibbs free energy G does not have any additional effect arising from the new parameter G, viz., (4.18) remains invariant with any of the formalisms considered, and so the second equation (4.19). In rest of the article we shall only explore the case where the parameter G is considered, since the other case does not satisfy the Smarr relation and we consider that case to be non-physical. Nevertheless, the second equation (4.19) implies that the massive charged BTZ black hole does not admit any critical Van der Waals behavior. While verifying the validity of the Reverse Isoperimetric Inequality [30], we find s r q 2 l2 l2 M̃ 2 cc1 2π V = 1− 2 − < 1, with A = 4S, (4.20) R= A 2π Gr+ r+ which means that in our case for all non-zero q and M̃ the Reverse Isoperimetric Inequality (4.20) is violated and, hence, the entropy exceeds the expected thermodynamic maximum giving rise to a super-entropic black hole. People sometimes try to restore the Reverse Isoperimetric Inequality by introducing a new thermodynamical variable, for instance, in [30] the authors include a renormalization length scale to make the isoperimetric ratio R ≥ 1. However, in our case we notice that because of the presence of the massive parameter M̃ it is never possible to make R ≥ 1 and, thus, the BTZ black hole in massive gravity remains super-entropic always. It should be noted that the super-entropic solution is sometimes even more interesting than the ordinary AdS black holes satisfying Reverse Isoperimetric Inequality. They give rise to non-compact event horizons with finite area, for further details on the interesting facts about the super-entropic black holes, one may see, for instance [41]. 7 BTZ black holes in massive gravity 5 Heat Engines In this section, we explore the behavior of a massive charged BTZ black hole (4.4) as a heat engine [32, 33]. Our main interest is to calculate the efficiency η of the cycle η= QC W =1− , QH QH (5.1) where W is the net output work and QH , QC are the net input and output heat flow, respectively. We will now consider a rectangular cycle in the P − V plane, which is familiar as the Carnot cycle. This rectangle will be in (P, V ) coordinates with the corner points described by (Pi , Vj ), where i = T, B stand for top and bottom and j = L, R denote the left and right. The efficiency of the Carnot cycle is [32] η= W (PT − PB )(VR − VL ) = ′ . QH M (PT , VR ) − M ′ (PB , VL ) − VL (PT − PB ) (5.2) The total amount of work done is given by the area enclosed by the rectangle. As we see that S in (4.7) and V in (4.13) are not independent, thus, the adiabats and isochores are the same [33]. Consequently, the Carnot and the Stirling cycles coincide. Nevertheless, in order to compute 2 the q efficiency (5.2), we first solve S from (4.13) as S = (M̃ cc1 ± A1 )/(32P G), with A1 (P, V ) = M̃ 4 c2 c21 + 32πP G(q 2 + 4P V ) and then substitute the solution in (4.12), so that we obtain the enthalpy as  
 M̃ 2 cc A ± M̃ 2 cc   8πP l2 G − 1 2 1 1 1 q G ∓ . M ′ (P, V ) = P V + 1 + ln(32πP ) + 2ln 4 64πP G 2 A1 ± M̃ 2 cc1 (5.3) ′ ′ Now, M (PT , VR ) and M (PB , VL ) can be obtained by replacing P by PT , PB and V by VR , VL in (5.3), respectively, so that we can compute the efficiency (5.2) as η= 64πPT PB G 2 (PT − PB )(VL − VR ) , N (5.4) with N = 16πPT PB G 2 " )# PB (A2 ± M̃ 2 cc1 )2 + M̃ 4 c2 c21 (PT − PB ) 4PT (VL − VR ) + q ln PT (A3 ± M̃ 2 cc1 )2    + M̃ 2 cc1 PB ±8πl2 PT G(A2 − A3 ) ∓ A2 ± A3 PT , 2 ( (5.5) where A2 = A1 (PT , VR ) and A3 = A1 (PB , VL ). Therefore, we notice that the mass term M̃ has a significant effect on the efficiency of the heat engine. The form of the efficiency in (5.4) along with (5.5) stands for a complete general expression for the massive BTZ black hole, which one can analyze further. 6 Concluding remarks We have studied a black hole solution of a charged massive BTZ black hole by using the Vainshtein and dRGT mechanism and analyzed the effects of such black hole in different scenarios. In particular, we have studied the Hawking radiation form this solution by utilizing the tunneling formalism. The solution that we obtain violates the Reverse Isoperimetric Inequality and, thus, they are superentropic. They also do not admit any critical Van der Walls behavior. Moreover, we have explored the black hole chemistry followed by the efficiency of a heat engine from the P − V diagram for this system. Our results show that unlike the ordinary BTZ black holes, the thermal fluctuation is 8 BTZ black holes in massive gravity absent in the massive gravity case. Thus, the BTZ black holes in massive gravity scenario are free from instabilities which is perhaps one of the greatest shortcomings for ordinary BTZ black holes. Furthermore, since the thermal fluctuations may be interpreted as quantum effects, the massive gravity does not feel important quantum effects. There are many interesting aspects which can be followed from our results. Firstly, it would be interesting to analyze the holographic entanglement entropy and holographic complexity dual to the BTZ AdS black hole. Secondly, the holographic conductivity has also been analyzed in certain aspects of massive gravity [42]. Therefore, it would be interesting to analyze such effects in our case also. Finally, our solution can be used to construct a CFT dual by using the standard formulation of AdS/CFT correspondence [43–47]. Acknowledgements: S.D. acknowledges the support of research grant (DST/INSPIRE/04/2016/ 001391) from the Department of Science and Technology, Govt. of India. References [1] A. G. Riess et al., Observational evidence from supernovae for an accelerating universe and a cosmological constant, Astron. J. 116, 1009 (1998). [2] S. Perlmutter et al., Discovery of a supernova explosion at half the age of the universe, Nature 391 , 51–54 (1998). [3] J. L. Tonry et al., Cosmological results from high-z supernovae, Astrophys. J. 594, 1 (2003). [4] P. J. E. Peebles and B. Ratra, The cosmological constant and dark energy, Rev. Mod. Phys. 75, 559 (2003). [5] E. J. Copeland, M. Sami and S. Tsujikawa, Dynamics of dark energy, Int. J. Mod. Phys. D 15, 1753 (2006). [6] J. A. Frieman, M. S. Turner and D. Huterer, Annu. Rev. Astron. Astrophys. 46, 385–432 (2008). [7] M. Khurshudyan, B. Pourhassan and A. Pasqua, Higher derivative corrections of f (R) gravity with varying equation of state in the case of variable G and Λ, Can. J. Phys. 93, 449–455 (2015). [8] H. van Dam and M. Veltman, Massive and mass-less Yang-Mills and gravitational fields, Nucl. Phys. B 22, 397–411 (1970). [9] Y. Iwasaki, Consistency condition for propagators, Phys. Rev. D 2, 2255 (1970). [10] S. Upadhyay, B. Pourhassan and H. Farahani, P − V criticality of first-order entropy corrected AdS black holes in massive gravity, Phys. Rev. D 95, 106014 (2017). [11] W. Pauli and M. Fierz, On relativistic field equations of particles with arbitrary spin in an electromagnetic field, Helv. Phys. Acta 12, 297 (1939). [12] A. I. Vainshtein, To the problem of nonvanishing gravitation mass, Phys. Lett. B 39, 393 (1972). [13] E. Babichev and C. Deffayet, An introduction to the Vainshtein mechanism, Class. Quantum Grav. 30, 184001 (2013). [14] D. G. Boulware and S. Deser, Can gravitation have a finite range? Phys. Rev. D 6, 3368 (1972). [15] K. Hinterbichler, Theoretical aspects of massive gravity, Rev. Mod. Phys. 84, 671 (2012). 9 BTZ black holes in massive gravity [16] C. de Rham, G. Gabadadze and A. J. Tolley, Resummation of massive gravity, Phys. Rev. Lett. 106, 231101 (2011). [17] S. F. Hassan, R. A. Rosen and A. Schmidt-May, Ghost-free massive gravity with a general reference metric, J. High Energy Phys. 2012, 026 (2012). [18] A. H. Chamseddine and V. Mukhanov, Massive hermitian gravity, J. High Energy Phys. 2012, 036 (2012). [19] S. H. Hendi, S. Panahiyan and B. E. Panah, Charged black hole solutions in Gauss-Bonnetmassive gravity, J. High Energy Phys. 2016, 129 (2016). [20] S. H. Hendi et al., Phase transition of charged black holes in massive gravity through new methods, Ann. Phys. (Berlin) 528, 819 (2016). [21] M. Wyman, W. Hu and P. Gratia, Self-accelerating massive gravity: time for field fluctuations, Phys. Rev. D 87, 084046 (2013). [22] M. S. Volkov, Stability of Minkowski space in ghost-free massive gravity theory, Phys. Rev. D 90, 024028 (2014). [23] M. Banados, C. Teitelboim and J. Zanelli, Black hole in three-dimensional spacetime, Phys. Rev. Lett. 69, 1849 (1992). [24] S. H. Hendi, B. E. Panah and S. Panahiyan, Massive charged BTZ black holes in asymptotically (a)dS spacetimes, J. High Energy Phys. 2016, 029 (2016). [25] S. W. Hawking, Particle creation by black holes, Commun. Math. Phys. 43, 220 (1975). [26] M. K. Parikh and F. Wilczek, Hawking radiation as tunneling, Phys. Rev. Lett. 85, 5042 (2000). [27] H. Kim, Hawking radiation as tunneling from charged black holes in 0A string theory, Phys. Lett. B 703, 94 (2011). [28] A. Yale, Exact Hawking radiation of scalars, fermions, and bosons using the tunneling method without back-reaction, Phys. Lett. B 697, 398 (2011). [29] K. Nozari and S. H. Mehdipour, Hawking radiation as quantum tunneling from a noncommutative Schwarzschild black hole, Class. Quantum Grav. 25, 175015 (2008). [30] A. M. Frassino, R. B. Mann and J. R. Mureika, Lower-dimensional black hole chemistry, Phys. Rev. D 92, 124069 (2015). [31] D. Kubiznak, R. B. Mann and M. Teo, Black hole chemistry: thermodynamics with Lambda, Class. Quantum Grav. 34, 063001 (2017). [32] R.A. Hennigar et al., Holographic heat engines: general considerations and rotating black holes, Class. Quantum Grav. 34, 175005 (2017). [33] C. V. Johnson, Holographic heat engines, Class. Quantum Grav. 31, 205002 (2014). [34] R.-G. Cai, Y.-P. Hu, Q.-Y. Pan and Y.-L. Zhang, Thermodynamics of black holes in massive gravity, Phys. Rev. D 91, 024032 (2015). [35] I. Arraut, Komar mass function in the de Rham–Gabadadze–Tolley nonlinear theory of massive gravity, Phys. Rev. D 90, 124082 (2014). 10 BTZ black holes in massive gravity [36] I. Arraut, On the apparent loss of predictability inside the de Rham-Gabadadze-Tolley nonlinear formulation of massive gravity: The Hawking radiation effect, Europhys. Lett. 109, 10002 (2015). [37] I. Arraut, Path-integral derivation of black-hole radiance inside the de-Rham-Gabadadze-Tolley formulation of massive gravity, Euro. Phys. J. C 77, 501 (2017). [38] I. Arraut, The black hole radiation in massive gravity, Universe 4, 27 (2018). [39] P. Painlevé, The classical mechanics and the theory of the relativity, Compt. Rend. Aca. Sci. (Paris), 173, 677–680 (1921). [40] L. Smarr, Mass formula for Kerr black holes, Phys. Rev. Lett. 30, 71 (1973). [41] R. A. Hennigar, D. Kubizk and R. B. Mann, Entropy inequality violations from ultraspinning black holes, Phys. Rev. Lett. 115, 031101 (2015). [42] L. Alberte and A. Khmelnitsky, Stability of massive gravity solutions for holographic conductivity, Phys. Rev. D 91, 046006 (2015). [43] A. Sinha, On the new massive gravity and AdS/CFT, J. High Energy Phys. 2010, 061 (2010). [44] J. Sadeghi and B. Pourhassan, Drag force of moving quark at N = 2 supergravity, J. High Energy Phys. 2008, 026 (2008). [45] B. Pourhassan and J. Sadeghi, STU–QCD correspondence, Can. J. Phys. 91, 995 (2013). [46] J. Sadeghi, B. Pourhassan and S. Heshmatian, Application of AdS/CFT in quark-gluon plasma, Adv. High Energy Phys. 2013, 759804 (2013). [47] V. Niarchos, Multi-string theories, massive gravity and the AdS/CFT correspondence, Fortsch. Phys. 57, 646 (2009) . 11