arxiv: v1 [hep-th] 29 Aug 2018

Transcription

1 Thermodynamics of charged accelerating AdS black holes and holographic heat engines Jialin Zhang, Yanjun Li, and Hongwei Yu arxiv: v1 [hep-th] 29 Aug 2018 Department of Physics and Synergetic Innovation Center for Quantum Effects and Applications, Hunan Normal University, Changsha, Hunan , China Abstract We explicitly formulate the consistent first law of thermodynamics for a charged accelerating Anti-de Sitter (AdS) black hole with a reasonable normalization of the timelike vector. The expression of renormalized thermodynamic mass is exactly derived, which displays an intertwining behavior with the charge in a complex way due to the special asymptotic structure of the gauge field. The thermodynamic length as a potential conjugate to the varying cosmic string tension is also introduced and analyzed in detail. Furthermore, we investigate the properties of holographic heat engines with charged accelerating AdS black holes as the working substance in a benchmarking scheme. We find that the efficiencies of the black hole heat engines can be influenced by both the size of the benchmark circular cycle and the cosmic string tension. More interestingly, the existence of charge may significantly increase the efficiencies of the black hole heat engines and make them be more sensitive to a varying cosmic string tension. A cross-comparison of the efficiencies of different black hole heat engines suggests that the acceleration also increases the efficiency and renders it more sensitive as the charge varies. PACS numbers: Dy,05.70.Ce,05.70.Fh Corresponding author at hwyu@hunnu.edu.cn

2 I. INTRODUCTION Ever since Hawking s pioneering work on the radiation of black holes [1, 2], the black hole thermodynamics has become an important topic of intensive research and been considered as a key to gaining insight into the quantum nature of gravity. In the last decades, the thermodynamics of several classes of black holes has been studied, and remarkable progresses have been made in, for example, connecting the laws of gravitation with the laws of thermodynamics [3, 4], deploying the black holes as a holographically dual system in quantum chromodynamics [5] and condensed matter physics [6, 7], and approaching the thermodynamics of black holes geometrically [8 17]. Recent studies show that the thermodynamical properties of black holes in anti-de Sitter (AdS) space are quite different from those in flat or de Sitter(dS) space. In particular, a subject dubbed black hole chemistry [18, 19], in which the negative cosmologic constant Λ is regarded as a thermodynamic pressure in the extended phase space [20 26], has started to attract a growing deal of interest. To date, the investigation in black hole chemistry has led to the discovery of a lot of interesting thermodynamic phenomena, such as the analogous phase of van der Waals fluid [27], re-entrant phase transitions [28] and the analogous superfluid phase transitions [29]. In the context of black hole chemistry, the concept of holographic heat engines which can extract work with the AdS black holes as the working substance in cycles in the pressurevolume phase plane was proposed by Johnson [30]. Subsequently, in order to better cross compare the efficiency of the heat engines with different black holes as the working material, a circular cycle of the heat engine in the pressure-volume plane has been introduced as a benchmarking prescription by Chakraborty and Johnson [31]. After Johnson s pioneering work, the concept of heat engines has been generalized to other different black holes (for an incomplete list of references, see[32 46]). Particularly, the properties of the holographic heat engine with accelerating AdS black holes as the working substance has drawn considerable interest. The accelerating AdS black holes are described by the so-called C metric [47 50], which can be used to study either one or two accelerating black holes, the pair creation of black holes [51], the construction of the black ring solution in 5-dimensional gravity [52], and the test of gravitational radiation [53]. However, the C metric, despite of the aforementioned various applications, remains some-

3 what esoteric in its thermodynamic prescription partly due to the fact that it has an uncommon asymptotic behavior and at least one non-removable conical singularity on the azimuthal axis of symmetry. As a result, discrepant and conflicting results have been obtained in the investigation of the thermodynamics of the accelerating AdS black holes [54 58]. For example, a set of thermodynamic variables was obtained in Refs. [54 56], where corresponding thermodynamic mass and temperature are associated with the coordinate time rather than with the normalized timelike Killing vector. However, later, the corresponding first law of thermodynamics was found to be incorrect [59]. Alternate expressions for the thermodynamic mass and temperature have also been posited by assuming integrability of a scaling of the mass and temperature in Ref. [57], but the cosmic string tension, which should be considered as a thermodynamic parameter relevant to the varying conical deficit, was not included in the first law. Furthermore, no physical interpretation was given either for the scaling or for why the energy content of the tension is thermodynamically irrelevant. Recently, A. Anabalon et al re-examined the thermodynamics of uncharged accelerating AdS black holes by means of the mass and temperature generated by a normalized time that just corresponds to the time of an asymptotic observer [58]. Such generated thermodynamic mass has been verified to be consistent with the mass obtained by holographic computing, and hence it reconciles discrepancies and conflicts in the previous investigations and confirms at this point that such choice of normalization of timelike Killing vector is applicable and well grounded. More recently, based on the correct thermodynamics of uncharged accelerating AdS black holes stated in Ref. [58], the properties of the corresponding benchmarking holographic heat engines have been studied in detail [46]. However, a formulation of the first law of thermodynamics and an examination of the holographic heat engines when the accelerating black hole is charged, where the additional special asymptotic structure of the gauge field renders the thermodynamics quite challenging, are quite challenging [58] and remain open. The purpose of the present paper is to fill this gap, i.e., we will first try to formulate the first law of thermodynamics for the charged accelerating AdS black holes and then study the influence of charge on the corresponding holographic heat engines. This paper is organized as follows. We begin in Sec.II by seeking the explicit expressions for the thermodynamic parameters of charged accelerating AdS black holes with an appropriate normalized time. In Sec.III, we investigate the properties of the benchmarking

4 holographic heat engines with charged accelerated AdS black holes as working substances. We end with a summary and a conclusion in Sec.IV. II. THE CHARGED ACCELERATING ADS BLACK HOLE AND ITS THERMO- DYNAMICS A charged accelerating AdS black hole is described by the following C metric [50, 54, 55] ds 2 = 1 [ ( f(r)dt 2 dr2 dθ 2 )] Ω 2 f(r) r2 g(θ) +g(θ)sin2 θ dφ2, (1) K 2 where and the gauge potential Ω = 1+Arcosθ, (2) f(r) = (1 A 2 r 2 ) (1 2m ) r + e2 + r2 r 2 l, 2 (3) g(θ) = 1+2mAcosθ+e 2 A 2 cos 2 θ, (4) B = e dt. (5) r Here, the parameters m, e and A are related to the mass, the electric charge and the magnitude of acceleration of the black hole respectively, l represents the AdS radius, and K characterizes the conical deficit of the spacetime. Specifically, the conical deficits on the spatial two poles, θ + = 0 (the north pole ) and θ = π (the south pole ) are [ δ ± = 2π 1 g(θ ] ±). (6) K Such conical deficits correspond to the tension of the cosmic string [54, 55] µ ± = δ ± 8π = 1 4 1±2mA+e2 A 2 4K. (7) In order to avoid the occurrence of negative tension defects, clearly we should require 0 µ + µ 1/4. To obtain the correct thermodynamics, we should first select an appropriate normalized time [60]. As far as metric (1) is concerned, such a normalization of the time coordinate can be chosen such that it corresponds to the time of an asymptotic observer, τ = αt, with α = 1 A 2 l 2 [58]. In the following discussions, we will try to formulate the exact

5 expressions of the thermodynamic parameters for charged accelerating black holes with the normalized time τ. The temperature T can be found by using the conventional Euclidean method associated with the normalized time, T = f (r + ) 4απ = m e2 + A2 m 2παr+ 2 2παr+ 3 2απ A2 r + 2απ + r + 2παl, (8) 2 where r + represents the black hole event horizon. Here, note that in order to have a welldefined temperature associated with the horizon of the black hole, we assume that the acceleration is small such that the acceleration horizon can be negated by a negative cosmological constant and only the black hole horizon exists [55]. As usual, the entropy S of the accelerating black hole is one quarter of the horizon area S = πr 2 + K(1 A 2 r 2 +). (9) Meanwhile, the thermodynamic electric charge Q of the black hole can be obtained by Q = 1 db = e 4π K, (10) and the thermodynamic pressure in extended black hole thermodynamics reads Ω=0 P = Λ 8π = 3 8πl 2. (11) With regard to the thermodynamic mass and the electric potential, the influence of the asymptotic structure of the gauge field should be taken into account. In fact, there is an electric flux from pole to pole at the AdS boundary (i.e., Arcosθ = 1) db = easinθdt dθ. (12) Such fluxes make the black hole mass and charge inextricably intertwined, therefore, a renormalization of the thermodynamic mass and the electrostatic potential is needed [55]. In order to obtain the correct thermodynamic variables for charged accelerating AdS black holes, we begin by taking the following variation δf = f δr + + f f f f δm+ δe+ δa+ δl = 0. (13) r + m e A l In addition, according to Eqs. (7),(9) and (10), the variation of some thermodynamic parameters yields δs = 2πr + 2πr K(1 A 2 r+) 2 2δr A πr 2 K(1 A 2 r+) 2 2δA + δk (1 A 2 r+) 2 K, (14) 2

6 m m ) K δa = δµ + +δµ Aδ(, K (15) (1 e 2 A 2 ) δk 2K = 2 A2 eδq e2 A 2 ( m )+ m δ K m δµ + + m+ae2 m δµ. (16) Eq. (13) can be rewritten, after a substitution of Eqs. (14)-(16) into it and some algebraic manipulations, as ( αm δ K where and ) = TδS+V H δp+φ H δq ( r + 1+Ar + Al 2) δµ + α ( r + 1 Ar + +Al 2) δµ α +mδk 2K 2 α, (17) V H := 4π [ r 3 ] + 3Kα (1 A 2 r+) 2 2 +ma2 l 4 (18) Φ H := e αr +. (19) Eq. (17) shows that as long as the accelerating black hole is charged, αm/k can not be identified as the thermodynamic mass due to the presence of variation of K. Although αm/k is consistent with the expression of mass derived in Ref. [58] by using the method of conformal completion [61, 62], and even satisfies the relation αm/k = 2TS +QΦ H 2PV H, (20) the polar electric field at the AdS boundary prevents us from carrying out the method of conformal completion, and a renormalization of the thermodynamic mass is still compulsory. Without loss of generality, it is convenient to assume that the renormalized thermodynamic mass, electric potential and volume respectively satisfy M = γ(e,a,l) m K, Φ = Φ H Φ 0, V = V H V e, (21) where γ(e,a,l), Φ 0 and V e are the corresponding auxiliary functions. In the extended phase space, the thermodynamic mass M usually is identified with the enthalpy rather than the internal energy. Thus, we have ( m δm = γδ + K) m ( e γδe+ A γδa+ l γδl ) K ( m = γδ +m e γδq+ K) me eγ δk + m Aγ K 2 K δa+ m lγ δl. (22) K

7 Moreover, according to Eq. (17), we can obtain ( αm TδS +VδP +ΦδQ λ + δµ + λ δµ = δ K [ [ λ + r + Al 2 A 2 l 2 r + α(1+ar + ) ] δµ + ) V e δp Φ 0 δq mδk 2αK 2 λ r + +Al 2 A 2 l 2 r + α(1 Ar + ) ] δµ, (23) where λ ± are the thermodynamic lengths introduced as conjugates to the string tensions in Refs. [55, 58]. Combining Eq. (22) and Eq. (23), we have after some algebra δm (TδS +VδP +ΦδQ λ + δµ + λ δµ ) ( = [γ α A A γ A2 l 2 + e2 +2e 3 α e γ )] ( m δ α 1 A 2 e K) [ 2 + 3V e 4πl + mla2 3 Kα + m ] lγ δl K [ ] A 2 em + Φ 0 + (1 A 2 e 2 )α + m(1+a2 e 2 ) 1 A 2 e 2 e γ δq [ r + + λ + A γ (1+Ar + )α + (1+2eα ] eγ) (1 A 2 e 2 )α (m Ae2 ) δµ + [ r + + λ + A γ (1 Ar + )α + (1+2eα ] eγ) (1 A 2 e 2 )α (m+ae2 ) δµ. (24) If the first law holds, the right-hand side of Eq. (24) must vanish. As a result, the analytic forms of the relevant thermodynamic parameters can be obtained in principle if γ(e, A, l) is given. Now, in order to get the exact expression of γ(e,a,l), we first check the Smarr formula M = γ m K = 2TS 2PV +QΦ = αm K +2PV e QΦ 0, (25) where Eq. (20) was used in the last equality. So Eq. (25) can be taken as another additional constraint equation in determining V e and Φ 0. From the condition that the right-hand of Eq. (24) vanishes, the expression of γ(e,a,l) can be explicitly obtained. This has been done in appendix A, where we have ( ) 1 A2 l γ(e,a,l) = 2 Ae 2Ae + 1+A 2 e 2 1+A 2 e 2E i arcsin (l+a 2 e 2 l) 2 1+A 2 e 2 4e ( ) 2 Ae 2Ae 1+A 2 e 2F i arcsin (l+a 2 e 2 l) 2 1+A 2 e 2, (26) 4e 2 where E i (ϕ m) and F i (ϕ m) are the corresponding incomplete elliptic integrals E i (ϕ m) = ϕ 0 (1 msin 2 θ) 1/2 dθ,f i (ϕ m) = ϕ 0 (1 msin 2 θ) 1/2 dθ. (27)

8 Then, the following thermodynamic parameters follow from the fact that each term in the right-hand side of Eq. (24) must vanish individually 3 4πr + V = V H V e = 3Kα(1 A 2 r +2 ) 4πl3 m γ 2 3K l, (28) Φ = Φ H Φ 0 = e αr + + λ ± = and they certainly obey the first law A 2 em α(1 A 2 e 2 ) + m(1+a2 e 2 ) 1 A 2 e 2 γ e, (29) r + α(1±ar + ) ± γ A (1+2eα eγ) (1 A 2 e 2 )α (m Ae2 ), (30) δm = TδS +VδP +ΦδQ λ + δµ + λ δµ. (31) For the case of vanishingly-small acceleration (i.e., Al 1 and ea 1), the proposed thermodynamic parameters can be approximated as αm M [1 2 ] (1+A 2 e 2 )K 3 A4 e 2 l 2, V 4π [ r+ 3 3Kα (1 A 2 r+) + A2 l 4 m ], A 2 e 2 Φ 1 [ e ] A 2 em(2α 1), α r + λ ± 1 [ r ] + m A(l 2 +e 2 )+A 2 e 2 m(4α 1). (32) α 1±Ar + In order to gain a better understanding of the thermodynamic behavior, let us recall the expression of the conical deficits in Eq. (6). As we see that the conical deficits in this spacetime are unequal on different axes, so it is impossible to choose a fixed K to simultaneously regularize the metric at both two poles. Such kind of irregularity along an axis is just precisely a definition of a conical singularity which encodes the information about the presence of cosmic string. Actually, the regularity of the metric at one pole demands K ± = g(θ ± ) = 1±2mA+e 2 A 2. (33) Here, we choose K = K + = 1+2mA+e 2 A 2 so that the metric is just regular on the north pole (i.e., µ + = 0), leaving a conical deficit on the south pole, then the string tension on the south pole yields µ = ma/k. We can see, from Fig.(1), that when the renormalized thermodynamic mass M is comparable to AdS radius l, the thermodynamic length λ becomes negative. And the thermodynamic length is not always a monotonic decreasing function of M in the region of M/l < 1,

9 (a) (b) FIG. 1: The thermodynamic length of the south pole is plotted as a function of the renormalized thermodynamic mass M for (a) the varying tension with a fixed Q/l = 0.10, and (b) the varying Q/l with a fixed tension µ = All the solid circle points represent the correspondingextremal black holes. andthe functional behavior depends onthe concrete values of µ andq. Furthermore, when M becomes very large for a fixed tension µ, the influence of the charge variation on the thermodynamic length almost disappears. Here, notice that if M/l is smaller than that of the extremal black hole, the temperature becomes negative, which is usually considered as unphysical. Moreover, it should be pointed out that as Al is drawing near one, the presence of the acceleration horizon makes the thermodynamic temperature ambiguous [54, 55, 58]. So, the yellow curve (µ = 0.15) in Fig. (1(a)) does not continue to arbitrarily small M/l (i.e., M/l or equivalently Al ) in order for the thermodynamic temperature of accelerating black holes to make sense. Neither does the red curve (Q/l = 0) in Fig. (1(b)) with M/l (or equivalently Al ). In Fig. (2), we have displayed the variation of the free energy F = M TS with the temperature for varying tension and charge in the canonical ensemble. For the case of a fixed Q/l = 0.10, as the tension µ is increased, the swallowtail style will be suppressed, leaving only two wings, which means that the increasing tension renders the charged accelerating black holes far away from extremality, and a spontaneous phase transition from a small black hole to a large black hole will occur. For a fixed tension, the prominent swallowtail exists in

10 (a) (b) FIG. 2: The free energy is plotted as a function of temperature for varying string tension with fixed Q/l = 0.10 in (a), and varying Q/l with fixed µ = 0.10 in (b). theregionwhereq/lbecomescomparabletoµ. AsQ/lisgraduallyincreasedordecreased, the swallowtail will disappear. Especially for large enough Q/l, only a single branch is left at last, which means that there is a critical point for Q/l at which the intermediate regime with negative specific heat for the charged accelerating black hole just disappears. III. BENCHMARKING BLACK HOLE HEAT ENGINES AND EFFICIENCY With the first law of charged accelerating AdS black holes derived, now we turn our attention to the features of the corresponding holographic heat engine in the benchmarking scheme. It is well-known that the efficiency for a valid heat cycle is defined by η = 1 Q C Q H, (34) where Q C denotes a net output heat flow in one cycle, and Q H represents a net input heat flow. In the benchmarking scheme [35, 45], a circular or elliptical cycle has been suggested to allow for cross-comparison of the efficiencies of holographic heat engines with different black holes as the working substances. In the P V plane, a general circle cycle can be written as P(θ) = P 0 +Rsinθ,V(θ) = V 0 +Rcosθ, (35)

11 where the circle is centered at (V 0,P 0 ) with a radius of R. Note that the relation R < P 0 and R < V 0 should hold for a physical heat cycle. As we have shown in the proceeding section, all the thermodynamic variables in the first law are certain combinations of the black hole parameters and the specific heat capacity at constant volume is nonzero (C V 0). Therefore, it is difficult to obtain the analytic expression for the efficiency of benchmarking black hole heat engines and a numerical integration of TdS will be needed. In order to better understand how the holographic heat engine of the accelerating charged black holes behaves, let us first review the influence of a charge on the temperature of the black holes. According to Eq. (8), it is easy to see that the contribution of a charge to the temperature is negative. Therefore, there exists an upper-bound of Q for the black hole thermodynamics. If Q approaches the upper-bound, the temperature becomes zero and the black hole becomes extremal. Since a valid cycle contour is forbidden to cross the area of non-positive temperature, a zero isothermal curve cannot be tangent to or intersecting with the circular cycle in the P V plane. Forabenchmarking circle in the P V plane, there must be a point of minim temperature on the circle. If we suppose an extremal black hole at the point of minim temperature on the cycle, then the corresponding value of Q of the extremal black hole can be considered as an extremal upper-bound of charge for the benchmarking circle. This upper-bound is denoted by Q max in our following discussions. It is not difficult to deduce that if the thermodynamic charge is larger than Q max, negative temperature occurs and the cycle process stops. Therefore, a benchmarking heat engine of black holes with charge can make physical sense only when Q < Q max. In Fig. (3), we have displayed how Q max changes as a function of µ and R. We can see that the upper-bound Q max is a decreasing function of µ and R. Let us return to the efficiency in the benchmarking circular cycle in Eq. (35), which will be influenced by the thermodynamic variables Q, µ and the circle radius R. We first present thebehavior of theefficiency withtheincreasing circle radiusrforfixed Qandµ in Fig. (4). Here, the Carnot efficiency is classically defined by η C = 1 T C /T H with maximum temperaturet H andminimtemperaturet C intheentirecycleprocess, andη D = 2π/(π+4)is proposed as the universal upper-bound for the efficiencies of benchmarking holographic heat engines on circular cycles in Ref. [45]. As we can see, the efficiency of the benchmarking

12 (a) (b) FIG. 3: Assuming the benchmarking circle is at (V 0,P 0 ) = (100,20) in the P V plane, the upper-bound of charge Q max is plotted as a function of (a)µ for fixed R = 10 and (b) R for fixed µ = 0.2. black hole heat engine η increases with the increase of R, and the corresponding Carnot efficiency η C is always larger than the efficiency η. However, due to the presence of charge, η can approximate to the upper-bound η D so long as R draws near rather than 20 in the uncharged case. A comment is in order here. According to Eq. (8) and Eq. (35), it is not difficult to derive that Q max, the upper-bound of charge, will approximate to 100 for µ = 0.2 and R = If R goes beyond , the Q max will become smaller than 100 (see Fig. (3(b))), then the circular cycle for Q = 100 will cross the region of negative temperature, whichmakesnosenseinphysics. Therefore, inthesettingofq = 100,µ = 0.2 and (V 0,P 0 ) = (100,20), R cannot go beyond , and the efficiency will approximate to the upper-bound 2π/(π+4) in the limit of R approaching In Fig. (5), the efficiencies are plotted as a function of the circle radius R for fixed cosmic string tension and fixed charge respectively. It is easy to find that the efficiency is still an increasing function of R. For a fixed R, the larger the Q or the larger the µ, the higher the efficiency. In Fig (6), we give a representative plot of the efficiency as a function of the cosmic string tensionµ forsomedifferent valuesofchargeq, which illustrates howthevariationofcharge influences the curved contour describing the behavior of the efficiency with increasing µ for a fixed benchmarking circle. Onecan see that when the valueof chargeqissmall a variation of µ does not cause significant changes in the efficiency as one would expect. However,

13 when the value of Q is large, the efficiency becomes an obvious increasing function of the cosmic string tension µ in contrast to the uncharged case [46]. Therefore, the presence of a charge may significantly increase the efficiency of the benchmarking heat engine with slowly accelerating AdS black holes as working substances. Now, we cross-compare the efficiency of different black hole heat engines in the benchmarking scheme. The black hole heat engines we choose are those of the slowly accelerating charged AdS black holes we just studied, the Einstein-Maxwell-AdS black holes (Reissner- Nordstom-like black hole) and the Born-Infeld black holes. We plot the efficiency of the benchmarking heat engine as a function of charge Q with different black holes as working substances in Fig. (7). It is worth noting that both Einstein-Maxwell-AdS black hole and Born-Infeld black hole have a vanishing specific heat at constant volume(i.e., C V = 0, the detailed thermodynamic quantities can be found in Ref. [35, 63 65]), and the holographic heat engine efficiency has been studied in Ref. [31, 45] FIG. 4: The efficiency of a heat engine of charged accelerating black holes is plotted as a function of R with fixed Q = 100 and µ = 0.2, where η C denotes Carnot efficiency and η D = 2π/(π +4) means the upper-bound value of efficiency. We have assumed the centre of the cycle satisfies (V 0,P 0 ) = (100,20). Here, for a valid heat engine cycle, R should be smaller than since the upper-bond Q max = 100 for R and µ = 0.2.

14 (a) (b) FIG. 5: Theefficiency of thebenchmarkingheat engine as a functionof theradiusrfor thecircular cycle is centered at (V 0,P 0 ) = (100,20). (a)thermodynamic charge Q = {1,50,100} with fixed cosmic string tension µ = 0.2. (b) The tension µ = {0.01,0.1,0.2} with fixed Q = FIG. 6: The efficiency of the holographic heat engine as a function of cosmic string tension µ with Q = {1,50,100}. The center of the circle is at (V 0,P 0 ) = (100,20) with a fixed radius R = 10. Fig. (7) reveals that the existence of charge in general increases the work efficiency to some extent. More interestingly, one finds that with acceleration, the holographic heat engine efficiency is usually larger than that without, and moreover it increases more rapidly

15 as the charge grows FIG. 7: The efficiencies vs electric charge Q with the center of benchmarking cycle localized at (V 0,P 0 ) = (100,20) and R = 10. The three solid curves respectively represent the efficiencies of three different families of black holes as the working material, including accelerating charged AdS (µ = 0.2), Born-Infeld (β = 1) and Einstein-Maxwell-AdS black holes. Note that the extreme charge for accelerating charged AdS black holes in this benchmarking cycle reads Q max , of which the efficiency will draw near η D = 2π/(4+π). Here, η C denotes the Carnot efficiency of accelerating charged black hole in this circular cycle. IV. CONCLUSION We have explicitly formulated the first law of thermodynamics for a charged accelerating AdS black hole within the framework where the time of an asymptotic observer is chosen as the normalized time. Due to the presence of a polar electric field at the AdS boundary, an additional renormalization of the thermodynamic mass and the electric potential is carried out. Through some algebraic manipulations in a tactic variation in the position of the horizon, we obtain the exact expression of the renormalized thermodynamic mass in the first

16 law, as well as other thermodynamic quantities, such as, the electric potential, the volume and the thermodynamic length. With the first law, the properties of the thermodynamic length and the possible phases of the charged accelerating black holes are also discussed. Moreover, we have examined the features of the heat engines with charged accelerating black holes as the working substances in the benchmarking scheme of holographic heat engines. Due to the fact that the corresponding thermodynamic variables in the first law are non-linear combinations of the black hole parameters r +, m, e, l and A, the numerical estimation is used to demonstrate the behavior of the efficiency of the holographic heat engines. We found that the efficiency of the holographic heat engine can be influenced by the cosmic string tension, the size of circular cycle and the thermodynamic charge. When the cosmic string tension is fixed, the efficiency of the holographic heat engine is increased by the existence of charge. For a fixed size of the circular cycle, the existence of charge may significantly increase the efficiency of the holographic heat engine as the string tension grows. A cross-comparison of the holographic heat engines with slowly accelerating charged AdS black holes and the Einstein-Maxwell-AdS black holes in a fixed benchmarking cycle shows that the presence of acceleration also increases the efficiency, and moreover, it makes the efficiency more sensitive to a varying charge Q. Acknowledgments This work was supported by the National Natural Science Foundation of China under Grants No and No Appendix A: Derivation of γ(e,a,l) In this appendix we will derive the explicit expression of γ(e,a,l) defined in Eq. (21). For this purpose, let us follow the Smarr formula Eq. (25) to get (γ α) m K = 2PV e QΦ 0. (A1) From the first three terms of the right-hand of Eq. (24), we find ( γ α A A γ A2 l 2 + e2 +2e 3 α e γ ) = 0, (A2) α 1 A 2 e 2

17 3V e 4πl + mla2 3 Kα + m lγ K = 0, A 2 em Φ 0 + (1 A 2 e 2 )α + m(1+a2 e 2 ) 1 A 2 e 2 e γ = 0. Substituting Eq. (A4) into Eq. (A1) leads to [ γ α V e = 4πl2 m 3K (A3) (A4) A 2 e 2 α(1 A 2 e 2 ) e(1+a2 e 2 ) e γ ]. (A5) 1 A 2 e 2 Combining Eq. (A3) and Eq. (A5) to eliminate parameter V e, we have α[(1 A 2 e 2 )(γ l l γ) e(1+a 2 e 2 ) e γ] 1 = 0. (A6) Thus, according to Eq. (A2) and Eq. (A6), we can straightforwardly derive the following results and A γ A = e γ e +l γ l, γ = l γ l A2 l 2 (1 A 2 e 2 ) + e(1+a2 e 2 ) γ 1 A 2 e 2 e (A7) (A8) where we have used α = 1 A 2 l 2. Here, based on the fact that the γ(e,a,l) is dimensionless, it is convenient to introduce two dimensionless parameters x = Al,y = Ae. So, we have γ A = l γ x +e γ y, γ l = A γ x, γ e = A γ y, (A9) and Eq. (A7) will automatically hold. Then Eq. (A8) will lead to one partial differential equation (xy 2 x) γ x = 1 1 x 2 +(y2 1)γ +y(1+y 2 ) γ y. (A10) The solutions for the above equation can be directly obtained as 1 x 2 γ = + y ( 1+y 2 1+y 2E i arcsin 2y ) x 2 (1+y 2 ) 2 1+y 2 4y 2 y ( 1+y 2F i arcsin 2y ) x 2 (1+y 2 ) 2 1+y 2, (A11) 4y 2 where E i (ϕ m) and F i (ϕ m) denote the incomplete elliptic integrals. Therefore, Eq. (26) can be straightforwardly verified by substituting x = Al and y = Ae into Eq. (A11). [1] S. Hawking, Commun. Math. Phys. 43, (1975).

18 [2] S. Hawking, Phys.Rev. D 13, (1976). [3] T. Jacobson, Phys. Rev.Lett. 75, 1260 (1995). [4] T. Padmanabhan, Rept. Prog. Phys. 73, (2010). [5] P.K. Kovtun, D. T. Son and A. O. Starinets, Phys. Rev. Lett. 94, (2005). [6] S. A. Hartnoll, P. K. Kovtun, M. Muller and S. Sachdev, Phys. Rev. B 76, (2007). [7] S. A. Hartnoll, C. P. Herzog and G. T. Horowitz, Phys. Rev. Lett. 101, (2008). [8] J.E. Aman, I. Bengtsson and N. Pidokrajt, Gen.Rel. Grav. 35, 1733 (2003). [9] H. Quevedo, Gen. Rel. Grav. 40, 971 (2008). [10] G. Ruppeiner, Springer Proc. Phys. 153, 179 (2014). [11] S.A.H. Mansoori and B. Mirza, Eur. Phys. J. C 74, 2681 (2014). [12] J. Suresh, R. Tharanath, N. Varghese and V.C. Kuriakose, Eur. Phys. J. C 74, 2819 (2014). [13] J.-L.Zhang, R.-G. Cai, and H. Yu, JHEP 02,143 (2015). [14] J.-L. Zhang, R.-G. Cai, and H. Yu, Phys. Rev. D 91, (2015). [15] S. H. Hendi, A. Sheykhi, S. Panahiyan, B. Eslam Panah, Phys. Rev. D 92, (2015). [16] C. Gruber, O. Luongo and H. Quevedo, arxiv: [gr-qc]. [17] R. Banerjee, B. R. Majhi and S. Samanta, Phys. Lett. B 767,25 (2017). [18] D. Kubiznak and R. B. Mann, Can. J. Phys. 93, 999 (2015). [19] R. B. Mann, Springer Proc. Phys. 170, 197 (2016). [20] C. Teitelboim, Phys.Lett. B 158,293 (1985). [21] M. M. Caldarelli, G. Cognola and D. Klemm, Class. Quant. Grav.17, 399 (2000). [22] D. Kastor, S. Ray and J. Traschen, Class.Quant.Grav. 26, (2009). [23] B. P. Dolan, Class. Quant. Grav. 28, (2011). [24] B. P. Dolan, Class. Quant. Grav. 28, (2011). [25] M. Cvetic, G. W. Gibbons, D. Kubiznak and C. N. Pope, Phys. Rev. D 84, (2011). [26] D. Kubiznak, R. B. Mann, and M. Teo, Class. Quant. Grav. 34, (2017). [27] D. Kubiznak and R. B. Mann, JHEP 07, 033 (2012). [28] N. Altamirano, D. Kubiznak, and R. B. Mann, Phys.Rev. D 88, (2013). [29] R. A. Hennigar, R. B. Mann, and E. Tjoa, Phys. Rev.Lett. 118, (2017). [30] C. V. Johnson, Class. Quant. Grav. 31, (2014). [31] A. Chakraborty and C. V. Johnson, arxiv: [hep-th]. [32] A. Belhaj, M. Chabab, H. EL Moumni, K. Masmar, M. B. Sedra, and A. Segui, JHEP 05,

19 149 (2015). [33] E. Caceres, P. H. Nguyen, and J. F. Pedraza, JHEP 09, 184 (2015). [34] M. R. Setare and H. Adami, Gen. Rel. Grav. 47, 133 (2015). [35] C. V. Johnson, Class. Quant.Grav. 33, (2016). [36] C. V. Johnson, Class. Quant. Grav. 33, (2016). [37] M. Zhang and W.-B. Liu, Int. J. Theor. Phys. 55, 5136 (2016). [38] J. Sadeghi and K. Jafarzade, Int. J. Theor. Phys. 56, 3387 (2017). [39] S.-W. Wei and Y.-X. Liu, arxiv: [gr-qc]. [40] S. H. Hendi, B. Eslam Panah, S. Panahiyan, H. Liu, and X. H. Meng, Phys. Lett. B 781, 40 (2018). [41] C. Bhamidipati and P. K. Yerra, Phys. Lett. B 772, 800 (2017). [42] H. Xu, Y. Sun, and L. Zhao, Int.J.Mod.Phys. D 26, (2017). [43] H. Liu and X.-H. Meng, Eur. Phys. J. C 77, 556 (2017). [44] J.-X. Mo, F. Liang, and G.-Q. Li, JHEP 03, 010 (2017). [45] R. A. Hennigar, F. McCarthy, A. Ballon, and R. B. Mann, Class. Quant. Grav. 34, (2017). [46] J. Zhang, Y. Li and H. Yu, Eur. Phys. J. C 78, 645 (2018). [47] W. Kinnersley and M. Walker, Phys. Rev. D 2, 1359 (1970). [48] J. F. Plebanski and M. Demianski, Ann. Phys.(N.Y.) 98, 98 (1976). [49] O. J. C. Dias and J. P. S. Lemos, Phys. Rev. D 67, (2003). [50] J. B. Griffiths and J. Podolsky, Int. J. Mod. Phys. D 15, 335(2006). [51] F. Dowker, J. P. Gauntlett, D. A. Kastor, and J. H. Traschen, Phys. Rev. D 49, 2909 (1994). [52] R. Emparan and H. S. Reall, Phys. Rev. Lett. 88, (2002). [53] J. Podolsky, M. Ortaggio, and P. Krtous, Phys.Rev. D 68, (2003). [54] M. Appels, R. Gregory and D. Kubiznak, Phys. Rev. Lett. 117, (2016). [55] M. Appels, R. Gregory and D. Kubiznak, JHEP 05, 116 (2017). [56] R. Gregory, J. Phys. Conf. Ser. 942, (2017). [57] M. Astorino, Phys. Rev. D 95, (2017). [58] A. Anabalon, M. Appels, R. Gregory, D. Kubiznak, R. B. Mann, A. Ovgun, arxiv: [hep-th]. [59] M. Appels, Thermodynamics of Accelerating Black Holes, Ph.D. thesis, Durham University

20 (2018). [60] G. W. Gibbons, M. J. Perry and C. N. Pope, Class. Quant. Grav. 22, 1503 (2005). [61] A. Ashtekar and S. Das, Class. Quant. Grav. 17, L17 (2000). [62] S. Das and R. Mann, JHEP 08, 033 (2000). [63] A. Chamblin, R. Emparan, C. V. Johnson, and R. C. Myers, Phys. Rev. D 60, (1999). [64] R.-G. Cai, D.-W. Pang, and A. Wang, Phys. Rev.D 70, (2004). [65] T. K. Dey, Phys. Lett. B 595, 484(2004).