This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

\UseRawInputEncoding
thanks: e-mail:03090137@sxdtdx.edu.cn,duyzh22@sxdtdx.edu.cn

Topology of the charged AdS black hole in restricted phase space

Han Wang1, Yun-Zhi Du2,3111the corresponding author 1College of Mathematics and Statistics, Shanxi Datong University, Datong 037009, China
2Department of Physics, Shanxi Datong University, Datong 037009, China
3Institute of Theoretical Physics, Shanxi Datong University, Datong, 037009, China
Abstract

The local topological properties of black hole systems can be expressed by the winding numbers as the defects. As so far, AdS black hole thermodynamics is often depicted by the dual parameters of (T,S),(P,V),(Φ,Q)(T,S),~{}(P,V),(\Phi,Q) in the extended phase space, while there is several study on the black hole thermodynamics in the restricted phase space. In this paper, we analyze the topological properties of the charged AdS black holes in the restricted phase space under the higher dimensions and higher order curvature gravities frame. The results show that the topological number of the charged black hole in the same canonical ensembles is a constant and is independent of the concrete dual thermodynamical parameters. However, the topological number in the grand canonical ensemble is different from that in the canonical ensemble for the same black hole system. Furthermore, these results are independent of the dimension dd, the highest order kk of the Lanczos-Lovelock densities.

I Introduction

As we know that after the last century black holes have been proposed to be not only the strong gravitational systems, but also the thermodynamic systems that satisfying the four laws of thermodynamics Bardeen1973 ; Bekenstein1973 . In 1983, Hawking and Page Hawking1983 proved a phase transition between a pure AdS spacetime and a stable large black hole state, i.e., Hawking-Page phase transition. Subsequently, it was explained to a confinement/deconfinement phase transition in the gauge theory Witten1998 . That leads to more attention being paid on black hole. Because of the lack of pressure in the traditional black hole thermodynamics, the authors of Kastor2009 expanded black hole thermodynamics into the expanded phase space by regarding the negative cosmological constant as pressure. Especially the thermodynamical properties of AdS black holes in the expanded phase space had been widely investigated Hendi2017a ; Hennigar2017a ; Frassin ; Kubiznak2012 ; Cai2013 ; Ma2017 ; Banerjee2017 ; Mann1207 ; Wei2015 ; Bhattacharya2017 ; Zeng2017 ; Zhang1502 ; Du2021 ; Zhang2020 . Recently, the holographic thermodynamics Visser2022 ; Ahmed2023 and the restricted phase space thermodynamics Zhao2022 ; Gao2022 ; Gao2022a ; Du2023 of AdS black holes had been proposed to give a holographic interpretation of black hole thermodynamics. On the other hand, based on the ϕ\phi-map topological flow theory Duan1979 the authors of Wei2022 presented that black hole solutions are just the defects that described by winding numbers. The corresponding winding number of the local stable black hole is one, while it is negative one for the local unstable black hole. The topological number is just the sum of all winding numbers to reveal the global topological nature of black hole. The black hole solutions may be classified by topological numbers. Based on these results, we will investigate the topological property of the charged AdS black holes under different gravity frames in the restricted phase space.

Recently, the authors of Kong2022 had checked out the restricted phase space thermodynamics in the higher dimensions and higher order curvature gravities Cai2002 ; Lanczos1932 ; Lovelock1971 . The subclass of Lanczos-Lovelock models with some particular choices of the coupling coefficients that known as the class of black hole scan models was adopted to be a simple example for the application of the extended phase space thermodynamics to the higher dimensions and higher order curvature gravities models with the character parameters of (d,kd,k). The parameter dd stands for the spacetime dimension and kk is the integer with the condition 1k[(d1)/2]1\leq k\leq[(d-1)/2] that represents the highest order of the Lanczos-Lovelock densities appearing in the action. In the restricted phase space, the thermodynamics of three typical models with (d,k)=(5,1),(5,2),(6,2)(d,k)=(5,1),~{}(5,2),~{}(6,2) which are representative of the Einstein-Hilert (EH), Chern-Simons (CS) and Born-Infield (BI) like gravity models, were investigated. From the thermodynamical results, it is proofed that the EH and BI like models seemly belong to the same universality class while the CS like models do not. On the other hand, RN-AdS black hole is the typical solution in the EH gravity model. Thence, in this work we will probe the topology of the RN-AdS black holes and check out whether the EH and BI like models belong to the same class from the perspective of the topology in the different ensembles.

Inspired by these, firstly we investigate the topology of the RN-AdS black hole under the restricted phase space frame in two different ensembles in Sec. II. Then from the perspective of the topology we check out whether two typical models with (d,k)=(5,1),(6,2)(d,k)=(5,1),~{}(6,2) belong to the same class in two different ensembles. A brief summary is given in Sec. IV.

II Topology of RN-AdS black hole in the restricted phase space

In this part we will discuss the topology of the RN-AdS black hole in different ensembles under the restricted phase space.

II.1 In the canonical ensemble

For the four-dimensional RN-AdS black hole, its metric has the following form

ds2=f(r)c2dt2+f1(r)dr2+r2(dθ2+sin2θdϕ2),\displaystyle ds^{2}=-f(r)c^{2}dt^{2}+f^{-1}(r)dr^{2}+r^{2}\left(d\theta^{2}+\sin^{2}\theta d\phi^{2}\right), (1)
h(r)=12GMr+Gq¯2r2+r2l2,Aμ=(ϕ¯(r)/c,0,0,0),Φ¯(r)=q¯r,\displaystyle h(r)=1-\frac{2GM}{r}+\frac{G\bar{q}^{2}}{r^{2}}+\frac{r^{2}}{l^{2}},~{}~{}A^{\mu}=(\bar{\phi}(r)/c,0,0,0),~{}~{}\bar{\Phi}(r)=\frac{\bar{q}}{r}, (2)

where ll is related to the cosmological constant via l2=3/Λl^{2}=-3/\Lambda. As M>q¯/GM>\bar{q}/\sqrt{G}, the function h(r)h(r) has two distinct real zeros at r=r±r=r_{\pm}, r+r_{+} corresponds to the black hole event horizon radius. The black hole mass parameter can be described as

M=r+2G(1+Gq¯2r+2+r+2l2).\displaystyle M=\frac{r_{+}}{2G}\left(1+\frac{G\bar{q}^{2}}{r_{+}^{2}}+\frac{r_{+}^{2}}{l^{2}}\right). (3)

In the restricted phase space for the RN-AdS black hole, the macro states are characterized by the following three pairs of dual thermodynamical parameters (S,T),(Q,Φ),(S,T),~{}(Q,\Phi), and (μ,N)(\mu,N). The corresponding first law of thermodynamics and smarr formula read

dM=TdS+ΦdQ+μdN,M=TS+ΦQ+μN\displaystyle dM=TdS+\Phi dQ+\mu dN,~{}~{}M=TS+\Phi Q+\mu N (4)

with

M(S,Q,N)\displaystyle M(S,Q,N) =\displaystyle= S2+πSN+π2Q22π3/2l(SN)1/2,T(S,Q,N)=3S2+πSNπ2Q24π3/2lS(SN)1/2,\displaystyle\frac{S^{2}+\pi SN+\pi^{2}Q^{2}}{2\pi^{3/2}l(SN)^{1/2}},~{}~{}T(S,Q,N)=\frac{3S^{2}+\pi SN-\pi^{2}Q^{2}}{4\pi^{3/2}lS(SN)^{1/2}}, (5)
Φ\displaystyle\Phi =\displaystyle= (πSN)1/2,μ=S2πSN+π2Q24π3/2lN(SN)1/2,\displaystyle\left(\frac{\pi}{SN}\right)^{1/2},~{}~{}\mu=-\frac{S^{2}-\pi SN+\pi^{2}Q^{2}}{4\pi^{3/2}lN(SN)^{1/2}}, (6)
G\displaystyle G =\displaystyle= l2N,S=πr+2G,Q=qlG.\displaystyle\frac{l^{2}}{N},~{}~{}S=\frac{\pi r_{+}^{2}}{G},~{}~{}Q=\frac{ql}{\sqrt{G}}. (7)

If the parameters S,Q,NS,~{}Q,~{}N are scaled as SλS,QλQ,NλNS\rightarrow\lambda S,~{}Q\rightarrow\lambda Q,~{}N\rightarrow\lambda N, from eq. (7) we can see that the mass parameter will also be scaled as MλMM\rightarrow\lambda M, and T,Φ,μT,~{}\Phi,~{}\mu will not be scaled. Hence the first order homogeneity of the mass parameter and the zeroth order homogeneity of T,Φ,μT,~{}\Phi,~{}\mu are clear as crystal. Note that the zeroth order homogeneous functions are intensive the critical point

Sc=πN6,Tc=63πl,Qc=N6,Fc=6N18l.\displaystyle S_{c}=\frac{\pi N}{6},~{}~{}T_{c}=\frac{\sqrt{6}}{3\pi l},~{}~{}Q_{c}=\frac{N}{6},~{}~{}F_{c}=\frac{\sqrt{6}N}{18l}. (8)

Introducing the relative parameters s=S/Sc,t=T/Tc,q=Q/Qcs=S/S_{c},~{}t=T/T_{c},~{}q=Q/Q_{c} and a Legendre transform of eq. (7), the relative parameter of free energy is given as Gao2022a

fFFc=q2+s2+6s4ts3/24s1/2,F(T,Q,N)=M(T,Q,N)TS.\displaystyle f\equiv\frac{F}{F_{c}}=\frac{q^{2}+s^{2}+6s-4ts^{3/2}}{4s^{1/2}},~{}~{}F(T,Q,N)=M(T,Q,N)-TS. (9)

In order to uncover the thermodynamical topology, the vector field mapping ϕ:X={(s,θ)|0<s<,\phi:~{}X=\{(s,\theta)|0<s<\infty, 0<θ<π}20<\theta<\pi\}\rightarrow\mathbb{R}^{2} is defined as Wei2020 ; Wei2022 ; Wei2022a ; Wu2023 ; Wu2023a

ϕ(s,θ)=(fs,cotθcscθ).\displaystyle\phi(s,\theta)=\left(\frac{\partial f}{\partial s},-\cot\theta\csc\theta\right). (10)

When the parameter ss is used to characterize AdS black hole, it becomes the first parameter of the domain of the mapping ϕ\phi. Another parameter θ\theta serves as an auxiliary function and it is utilized to construct the second component of the mapping. The component ϕθ\phi^{\theta} is divergent at θ=0,π\theta=0,~{}\pi, thus the direction of the vector points outward there. It is obvious that the zero point of ϕ\phi corresponds to the black hole with the temperature of T=τ1T=\tau^{-1} as θ=π/2\theta=\pi/2. Therefore, the zero point of the mapping can be used to characterize the black hole solution with a given parameter τ\tau. Based on the Duan’s ϕ\phi-mapping topological current theory Duan1979 ; Duan1984 , the zero points of the mapping ϕ\phi is linked to the topological number. The topological number can be obtained by the weighted sum of the zero points of the mapping ϕ\phi. The weight of each zero point is determined by its nature. A saddle point has a weight of negative one, and conversely, the weight of an extremum point is one. A topological current can be described as the following form

jμ=12eμνρϵabνnaρnb,μ,ν,ρ=0,1,2,a,b=1,2,\displaystyle j^{\mu}=\frac{1}{2}e^{\mu\nu\rho}\epsilon_{ab}\partial_{\nu}n^{a}\partial_{\rho}n^{b},~{}~{}~{}~{}\mu,\nu,\rho=0,1,2,~{}~{}a,b=1,2, (11)

where nn is the unit vector (ns,nθ)(n^{s},n^{\theta}) with ns=ϕs/ϕn^{s}=\phi^{s}/\parallel\phi\parallel and nθ=ϕθ/ϕn^{\theta}=\phi^{\theta}/\parallel\phi\parallel. And it satisfies the conserve law: μjμ=0\partial_{\mu}j^{\mu}=0. As shown in Refs. Duan1998 ; Fu2000 , the topological current is a δ\delta-function of the field configuration

jμ=δ2(ϕ)Jμ(ϕx),\displaystyle j^{\mu}=\delta^{2}(\phi)J^{\mu}\left(\frac{\phi}{x}\right), (12)

where the three-dimensional Jacobian reads ϵabJμ(ϕx)=ϵμνρνϕaρϕb\epsilon^{ab}J^{\mu}\left(\frac{\phi}{x}\right)=\epsilon^{\mu\nu\rho}\partial_{\nu}\phi^{a}\partial_{\rho}\phi^{b}. As ϕa(xi)=0\phi^{a}(x_{i})=0, the topological current equals to zero. The topological number in a parameter region Σ\Sigma can be calculated by the following expression

W=Σj0d2x=i=1N¯βiηi=i=1N¯wi,\displaystyle W=\int_{\Sigma}j^{0}d^{2}x=\sum\limits^{\bar{N}}_{i=1}\beta_{i}\eta_{i}=\sum\limits^{\bar{N}}_{i=1}w_{i}, (13)

where j0=i=1N¯βiηiδ2(xsi)j^{0}=\sum\limits^{\bar{N}}_{i=1}\beta_{i}\eta_{i}\delta^{2}\left(\overrightarrow{x}-\overrightarrow{s}_{i}\right) is the density of the topological current, βi\beta_{i} is the Hopf index which is always positive. sis_{i} is the ii-th zero point of the mapping ϕ\phi, ηi=signJ0(ϕ/x)si=±1\eta_{i}=signJ^{0}(\phi/x)_{s_{i}}=\pm 1 is the Brouwer degree. For the ii-th zero point of the vector, the winding number wiw_{i} is determined by the stability of black hole state.

In the canonical ensemble, substituting eq. (9) into eq. (10), we have

ϕs=fs=q28s3/2+3s1/28+34s1/21τ,ϕθ=cotθcscθ.\displaystyle\phi^{s}=\frac{\partial f}{\partial s}=-\frac{q^{2}}{8s^{3/2}}+\frac{3s^{1/2}}{8}+\frac{3}{4s^{1/2}}-\frac{1}{\tau},~{}~{}\phi^{\theta}=-\cot\theta\csc\theta. (14)

The zero points are determined by the expression of ϕs=0\phi^{s}=0. So we have

1τ=q28s3/2+3s1/28+34s1/2.\displaystyle\frac{1}{\tau}=-\frac{q^{2}}{8s^{3/2}}+\frac{3s^{1/2}}{8}+\frac{3}{4s^{1/2}}. (15)

It is obvious that there are two extremes: τmin,max=2(11q2)3/23(11q2)q2\tau_{min,max}=\frac{2\left(1\mp\sqrt{1-q^{2}}\right)^{3/2}}{3(1\mp\sqrt{1-q^{2}})-q^{2}}. As q=1q=1, τmin=τmax=τc=1\tau_{min}=\tau_{max}=\tau_{c}=1. Note that the generation point satisfies the constraint conditions: τs=0,2τs2>0\frac{\partial\tau}{\partial s}=0,~{}\frac{\partial^{2}\tau}{\partial s^{2}}>0, and the annihilation point obeys the conditions: τs=0,2τs2<0\frac{\partial\tau}{\partial s}=0,~{}\frac{\partial^{2}\tau}{\partial s^{2}}<0. The zero points of ϕs\phi^{s} in the diagram of τs\tau-s is displayed in Fig. 1. τmin\tau_{min} is related with the generation point, and τmax\tau_{max} is of the annihilation point. As τmin<τ<τmax\tau_{min}<\tau<\tau_{max}, there are three intersection points for the RN-AdS black hole in the canonical ensemble. The black, red dashed, blue lines are for the low-potential black hole (LPSB), intermediate-potential black hole (IPBH), and high-potential black hole (HPBH), respectively. They are just the stable, unstable, and stable black hole states. The corresponding winding numbers are 1, -1, 1, so the topological number is W=11+1=1W=1-1+1=1 as τmin<τ<τmax\tau_{min}<\tau<\tau_{max}. The intersection points exactly are satisfying the equation τ=1/T\tau=1/T. For τ=τmin,max\tau=\tau_{min,max}, the three intersection points for the RN black hole coincide. And they will disappear as τ<τmin\tau<\tau_{min} or τ>τmax\tau>\tau_{max}, the three branches black holes reduce to the one stable black hole and the topological number is still one. These results are consistent with that in Ref. Fan2023 .

Refer to caption
Figure 1: Zero points of ϕs\phi^{s} in the diagram of τs\tau-s with q=0.4q=0.4 for the RN-AdS black hole in the canonical ensemble.

II.2 In the grand canonical ensemble

Now we proceed to investigate the topology of the RN-AdS black hole in the grand canonical ensemble. For the given values of Φ\Phi and NN, from eq. (7) the minimum temperature reads Gao2022a

Tmin=32(1l2Φ2)1/2πl,\displaystyle T_{min}=\frac{\sqrt{3}}{2}\frac{\left(1-l^{2}\Phi^{2}\right)^{1/2}}{\pi l}, (16)

the corresponding minimum entropy, mass parameter, and free energy are

Smin=πN3(1l2Φ2),Mmin=N(2+l2Φ2)33l(1l2Φ2)1/2,Fmin=N(1+5l2Φ2)63l(1l2Φ2)1/2.\displaystyle S_{min}=\frac{\pi N}{3}\left(1-l^{2}\Phi^{2}\right),~{}~{}M_{min}=\frac{N\left(2+l^{2}\Phi^{2}\right)}{3\sqrt{3}l}\left(1-l^{2}\Phi^{2}\right)^{1/2},~{}~{}F_{min}=\frac{N\left(1+5l^{2}\Phi^{2}\right)}{6\sqrt{3}l}\left(1-l^{2}\Phi^{2}\right)^{1/2}. (17)

By introducing the relative parameters t=T/Tmin,s=S/Smin,m=M/Mmin,f=F/Fmint=T/T_{min},~{}s=S/S_{min},~{}m=M/M_{min},~{}f=F/F_{min}, we can obtain the express

f=s1/2[3+s+(3s)l2Φ2]1+5l2Φ23ts(1l2Φ2)1/21+5l2Φ2,\displaystyle f=\frac{s^{1/2}\left[3+s+(3-s)l^{2}\Phi^{2}\right]}{1+5l^{2}\Phi^{2}}-\frac{3ts\left(1-l^{2}\Phi^{2}\right)^{1/2}}{1+5l^{2}\Phi^{2}}, (18)

so the component of the mapping ϕs\phi^{s} reads

ϕs=3[1+s(1l2Φ2)+l2Φ2]2s1/2(1+5l2Φ2)3t(1l2Φ2)1/21+5l2Φ2.\displaystyle\phi^{s}=\frac{3\left[1+s(1-l^{2}\Phi^{2})+l^{2}\Phi^{2}\right]}{2s^{1/2}(1+5l^{2}\Phi^{2})}-\frac{3t\left(1-l^{2}\Phi^{2}\right)^{1/2}}{1+5l^{2}\Phi^{2}}. (19)

Hence, we have

τ=2s1/2(1l2Φ2)1/2[1+s(1l2Φ2)+l2Φ2].\displaystyle\tau=\frac{2s^{1/2}\left(1-l^{2}\Phi^{2}\right)^{1/2}}{\left[1+s(1-l^{2}\Phi^{2})+l^{2}\Phi^{2}\right]}. (20)

There exist one extreme: τextr=(1l4Φ4)(1+l4Φ4)\tau_{extr}=\sqrt{\frac{\left(1-l^{4}\Phi^{4}\right)}{\left(1+l^{4}\Phi^{4}\right)}}. The zero points of ϕs\phi^{s} in the diagram of τs\tau-s is displayed in Fig. 2. τextr\tau_{extr} is related with the annihilation point. As τ<τextr\tau<\tau_{extr}, there are two intersection points for the RN-AdS black hole in the grand canonical ensemble. The red and blue lines are for the low-potential black hole (LPSB) and high-potential black hole (HPBH), respectively. They are just the stable and unstable black hole states. The corresponding winding numbers are 1, -1, so the topological number is W=11=0W=1-1=0 as τ<τextr\tau<\tau_{extr}. For τ=τextr\tau=\tau_{extr}, the two intersection points for the RN black hole coincide. Furthermore, there is no black hole as τ>τextr\tau>\tau_{extr}.

Refer to caption
Figure 2: Zero points of ϕs\phi^{s} in the diagram of τs\tau-s with l2Φ2=0.1l^{2}\Phi^{2}=0.1 for the RN-AdS black hole in the grand canonical ensemble.

III Topology of the AdS black hole in the higher dimensional and higher curvature gravity

Recently, the restricted phase space thermodynamics is shown to be applicable to a large class of the higher dimensional and higher curvature gravity models with the coupling to the Maxwell field Kong2022 ; Gogoi2023 , which are known as black hole scan models Crisostomo2000 and are labeled by the spacetime dimension dd and the highest order kk of the Lanczos-Lovelock densities appearing in the action. Three typical example cases with (d,k)=(5,1)(d,k)=(5,1), (5,2)(5,2), and (6,2)(6,2) are chosen as example cases and studied in some detail. These cases are representatives of Einstein-Hilbert, Chern-Simons and Born-Infield like gravity models. However, the Chern-Simons like (5,2)(5,2)- model behaves quite differently. This seems to indicate that the Einstein-Hilbert and Born-Infield like models belong to the same universality class while the Chern-Simons like models do not. In this part we will check out whether the EH and BI like models do belong the same universality class from topology.

To be more concrete, the metric function of the black holes in the higher dimensions and higher curvature gravities with the coupling to Maxwell’s field was given in Refs. Kong2022 ; Gogoi2023 ; Crisostomo2000

ds(d,k)2\displaystyle ds^{2}_{(d,k)} =\displaystyle= f(d,k)(r)dt2+f(d,k)1(r)dr2+r2dΩd22,\displaystyle-f_{(d,k)}(r)dt^{2}+f_{(d,k)}^{-1}(r)dr^{2}+r^{2}d\Omega^{2}_{d-2},
f(d,k)(r)\displaystyle f_{(d,k)}(r) =\displaystyle= 1+r2l2(2G(d,k)M+δd2,k,1rd2k1G(d,k)Q2(d3)r2(dk2))1/k,Φ=Q(d3)rd3,\displaystyle 1+\frac{r^{2}}{l^{2}}-\left(\frac{2G_{(d,k)}M+\delta_{d-2,k,1}}{r^{d-2k-1}}-\frac{G_{(d,k)}Q^{2}}{(d-3)r^{2(d-k-2)}}\right)^{1/k},~{}~{}\Phi=\frac{Q}{(d-3)r^{d-3}}, (21)

where G(d,k)G_{(d,k)} is the Newton constant, and δd2,k,1\delta_{d-2,k,1} is the Kronecker delta function. QQ and MM are the charge and mass parameters of black hole. Note that for k=1k=1, the three dimensional black hole Banados1992 ; Banados1993 and Schwarzschild-AdS solutions of the higher dimensional Einstein-Hilbert action with the negative cosmological constant both can be recovered. And when setting k=[d12]k=\left[\frac{d-1}{2}\right], we also can obtain the black hole solutions corresponding to Born-Infeld (BI) and Chern-Simons (CS) theories Banados1994 .

The black hole horizon rhr_{h} is located at one of the zero points of f(d,k)(r)f_{(d,k)}(r). By solving the equation of f(d,k)(rh)=0f_{(d,k)}(r_{h})=0, the black hole mass parameter can be obtained as the following form

M(rh,G(d,k),Q)=δd2,k,12G(d,k)+rhd2k1(1+rh2/l2)k2G(d,k)+Q22(d3)rhd3.\displaystyle M(r_{h},G_{(d,k)},Q)=-\frac{\delta_{d-2,k,1}}{2G_{(d,k)}}+\frac{r_{h}^{d-2k-1}\left(1+r_{h}^{2}/l^{2}\right)^{k}}{2G_{(d,k)}}+\frac{Q^{2}}{2(d-3)r_{h}^{d-3}}. (22)

The temperature of the black hole can be evaluated by the Euclidean period method, which gives

T=14πkrh2d1(1+rh2/l2)1k[(d2k1)rh2d(1+rh2/l2)kG(d,k)Q2rh2k+4]+rh2πl2.\displaystyle T=\frac{1}{4\pi kr_{h}^{2d-1}}\left(1+r_{h}^{2}/l^{2}\right)^{1-k}\left[(d-2k-1)r_{h}^{2d}\left(1+r_{h}^{2}/l^{2}\right)^{k}-G_{(d,k)}Q^{2}r_{h}^{2k+4}\right]+\frac{r_{h}}{2\pi l^{2}}. (23)

The entropy of the black hole was given in Ref. Gogoi2023 as the following form

S=2πkG(d,k)0rhrd2k+1(1+r2/l2)k1𝑑r.\displaystyle S=\frac{2\pi k}{G_{(d,k)}}\int_{0}^{r_{h}}r^{d-2k+1}\left(1+r^{2}/l^{2}\right)^{k-1}dr. (24)

For establishing the restricted phase space thermodynamics (RPST) of the (d,k)(d,k)-models, we introduce two thermodynamical parameters: the effective number of the microscopic degrees of freedom NN and the chemical potential μ\mu, which are defined as

N=Ld2kG(d,k),μ=(MTSΦQ)/N.\displaystyle N=\frac{L^{d-2k}}{G_{(d,k)}},~{}~{}\mu=(M-TS-\Phi Q)/N. (25)

Here, LL is a constant that is introduced to make NN dimensionless. The corresponding re-scaled charge and potential are

Q~=QL(d2k)/2G(d,k),Φ~=QG(d,k)(d3)L(d2k)/2rhd3.\displaystyle\widetilde{Q}=\frac{QL^{(d-2k)/2}}{\sqrt{G_{(d,k)}}},~{}~{}\widetilde{\Phi}=\frac{Q\sqrt{G_{(d,k)}}}{(d-3)L^{(d-2k)/2}r_{h}^{d-3}}{\color[rgb]{1,0,0}.} (26)

The thermodynamical parameters are satisfied the following expressions

M=TS+Φ~Q~+μN,dM=TdS+Φ~dQ~+μdN.\displaystyle M=TS+\widetilde{\Phi}\widetilde{Q}+\mu N,~{}~{}dM=TdS+\widetilde{\Phi}d\widetilde{Q}+\mu dN. (27)

Note that the above equations indicate that the RPST formalism holds on for the (d,k)(d,k)-models with d>3d>3.

III.1 Topology of RPST for (5,1)(5,1)-model: EH gravity coupled to Maxwell’s field in the canonical ensemble

Setting (d,k)=(5,1)(d,k)=(5,1), the eq. (21) becomes

f(5,1)=1+r2l2+G(5,1)Q22r42G(5,1)Mr2.\displaystyle f_{(5,1)}=1+\frac{r^{2}}{l^{2}}+\frac{G_{(5,1)}Q^{2}}{2r^{4}}-\frac{2G_{(5,1)}M}{r^{2}}. (28)

This is the black hole solution in the Einstein-Hilbert gravity theory, and the entropy is S=A/4GS=A/4G with the Newton constant G=3πG(5,1)/4G=3\pi G_{(5,1)}/4. Through the transformations

rh=[3SL3/(2πN)]1/3,G(5,1)=L3/N,Q=Q~/N,S~=(3SL3/N)1/3,\displaystyle r_{h}=[3SL^{3}/(2\pi N)]^{1/3},~{}G_{(5,1)}=L^{3}/N,~{}Q=\widetilde{Q}/\sqrt{N},~{}\widetilde{S}=(3SL^{3}/N)^{1/3}, (29)

the temperature and mass parameter of the black hole can be rewritten as

T\displaystyle T =\displaystyle= S~2(2π)2/3(2+4S~2(2π)2/3/l2(2π)4/3L3Q~2N2S~4),\displaystyle\frac{\widetilde{S}}{2(2\pi)^{2/3}}\left(2+4\widetilde{S}^{2}(2\pi)^{-2/3}/l^{2}-\frac{(2\pi)^{4/3}L^{3}\widetilde{Q}^{2}}{N^{2}\widetilde{S}^{4}}\right), (30)
M\displaystyle M =\displaystyle= S~2N2(2π)2/3L3(1+(2π)2/3S~2l2+(2π)4/3L3Q~22N2S~4).\displaystyle\frac{\widetilde{S}^{2}N}{2(2\pi)^{2/3}L^{3}}\left(1+\frac{(2\pi)^{-2/3}\widetilde{S}^{2}}{l^{2}}+\frac{(2\pi)^{4/3}L^{3}\widetilde{Q}^{2}}{2N^{2}\widetilde{S}^{4}}\right). (31)

It is transparent that when S~λS~,Q~λQ~,NλN\widetilde{S}\rightarrow\lambda\widetilde{S},~{}\widetilde{Q}\rightarrow\lambda\widetilde{Q},~{}N\rightarrow\lambda N, the mass parameters behaves as MλMM\rightarrow\lambda M, while the temperature is kept unchanged under these re-scalings. From these thermodynamical parameters, we can obtain the critical point

S~c=(2π)1/3l3,Q~c=2135l2NL3/2,Tc=435πl,Mc=7l2N30L3,Fc=l2N18L3.\displaystyle\widetilde{S}_{c}=\frac{(2\pi)^{1/3}l}{\sqrt{3}},~{}~{}\widetilde{Q}_{c}=\sqrt{\frac{2}{135}}\frac{l^{2}N}{L^{3/2}},~{}~{}T_{c}=\frac{4\sqrt{3}}{5\pi l},~{}~{}M_{c}=\frac{7l^{2}N}{30L^{3}},~{}~{}F_{c}=\frac{l^{2}N}{18L^{3}}. (32)
Refer to caption
Figure 3: Zero points of ϕs\phi^{s} in the diagram of τs\tau-s with q=0.8q=0.8 for the higher dimensions and higher order curvature gravity model with (d,k)=(5,1)(d,k)=(5,1) in the canonical ensemble.

With the relative parameters t=T/Tc,s=S/Sc,q=Q/Qc,ϕ~=Φ~/Φ~c,m=M/Mc,f=F/Fct=T/T_{c},~{}s=S/S_{c},~{}q=Q/Q_{c},~{}\widetilde{\phi}=\widetilde{\Phi}/\widetilde{\Phi}_{c},~{}m=M/M_{c},~{}f=F/F_{c}, the relative free energy reads

f=mMc/FctsTcSc/Fc=s4/3+3s2/3+q25s2/316s5τ.\displaystyle f=mM_{c}/F_{c}-tsT_{c}S_{c}/F_{c}=s^{4/3}+3s^{2/3}+\frac{q^{2}}{5s^{2/3}}-\frac{16s}{5\tau}. (33)

From the above equation, the components of the mapping ϕ\phi can be written as

ϕs=4s1/33+2s1/32q215s5/3165τ,ϕθ=cotθcscθ.\displaystyle\phi^{s}=\frac{4s^{1/3}}{3}+\frac{2}{s^{1/3}}-\frac{2q^{2}}{15s^{5/3}}-\frac{16}{5\tau},~{}~{}\phi^{\theta}=-\cot\theta\csc\theta. (34)

The zero points of the mapping ϕ\phi can be calculated by ϕs=0\phi^{s}=0. We solve it and get

τ=5s1/312+58s1/3q224s5/3.\displaystyle\tau=\frac{5s^{1/3}}{12}+\frac{5}{8s^{1/3}}-\frac{q^{2}}{24s^{5/3}}. (35)

In the Fig. 3, we can see that there are three branches as τmin<τ<τmax\tau_{min}<\tau<\tau_{max}, the large black hole branch for τ<τmin\tau<\tau_{min}, and the small black hole branch for τ>τmax\tau>\tau_{max}. That means the existence of a phase transition in the region of τmin<τ<τmax\tau_{min}<\tau<\tau_{max}. The winding numbers of the large black hole branch for τ<τmin\tau<\tau_{min} and small black hole branch for τ>τmax\tau>\tau_{max}, which are stable, both read w=1w=1. And the corresponding topological number equals to one. While for the system in the region of τmin<τ<τmax\tau_{min}<\tau<\tau_{max}, there exist the stable large and small black hole, and the unstable intermediate black hole, thus the corresponding winding number are w=1,1,1w=1,1,-1. The topological number is W=1W=1. These results are consistent with that for the RN-AdS black hole in the canonical ensemble.

III.2 Topology of RPST for (5,1)(5,1)-model: EH gravity coupled to Maxwell’s field in the grand canonical ensemble

In the grand canonical ensemble, using the eqs. (25), (26), and (29), we can obtain

Q2=4NΦ~2(2π)4/3(3L3SN)4/3,S~=L(3SN)1/3,\displaystyle Q^{2}=\frac{4N\widetilde{\Phi}^{2}}{(2\pi)^{4/3}}\left(\frac{3L^{3}S}{N}\right)^{4/3},~{}~{}\widetilde{S}=L\left(\frac{3S}{N}\right)^{1/3}, (36)

then substituting them into eqs. (30) and (31), the temperature and mass parameter become

T\displaystyle T =\displaystyle= (3π2NL32S)1/3(2+4L2l2[3S2πN]2/34Φ~2L3),\displaystyle\left(\frac{3\pi^{2}NL^{3}}{2S}\right)^{1/3}\left(2+\frac{4L^{2}}{l^{2}}\left[\frac{3S}{2\pi N}\right]^{2/3}-4\widetilde{\Phi}^{2}L^{3}\right), (37)
M\displaystyle M =\displaystyle= N2L(3S2πN)2/3(1+L2l2[3S2πN]2/3+2Φ~2L3).\displaystyle\frac{N}{2L}\left(\frac{3S}{2\pi N}\right)^{2/3}\left(1+\frac{L^{2}}{l^{2}}\left[\frac{3S}{2\pi N}\right]^{2/3}+2\widetilde{\Phi}^{2}L^{3}\right). (38)

From the above equation, the minimum temperature is determined by the expression TS=0\frac{\partial T}{\partial S}=0, yielding

Tmin=2(12Φ~2L3)πl,\displaystyle T_{min}=\frac{\sqrt{2(1-2\widetilde{\Phi}^{2}L^{3})}}{\pi l}, (39)

the corresponding minimum entropy, mass and free energy read

Smin=πNl3(12Φ~2L3)3/232L3,Mmin=Nl2(12L3Φ~2)(3+2L3Φ~2)8L3,Fmin=N(1+5l2Φ~2)(1Φ~2l2)1/263l.\displaystyle S_{min}=\frac{\pi Nl^{3}\left(1-2\widetilde{\Phi}^{2}L^{3}\right)^{3/2}}{3\sqrt{2}L^{3}},~{}~{}M_{min}=\frac{Nl^{2}(1-2L^{3}\widetilde{\Phi}^{2})(3+2L^{3}\widetilde{\Phi}^{2})}{8L^{3}},~{}~{}F_{min}=\frac{N(1+5l^{2}\widetilde{\Phi}^{2})\left(1-\widetilde{\Phi}^{2}l^{2}\right)^{1/2}}{6\sqrt{3}l}. (40)
Refer to caption
Figure 4: Zero points of ϕs\phi^{s} in the diagram of τs\tau-s with L3Φ~2=0.1L^{3}\widetilde{\Phi}^{2}=0.1 for the AdS black hole under the (5,1)-model frame in the grand canonical ensemble.

With the definitions f=F/Fmin=(MTS)/Fmin,1/τ=t=T/Tmin,s=S/Sminf=F/F_{min}=(M-TS)/F_{min},~{}1/\tau=t=T/T_{min},~{}s=S/S_{min}, we can obtain the components of the mapping ϕ\phi

ϕs=fs=1Fmin(MsTminSminτ),ϕθ=cotθcscθ.\displaystyle\phi^{s}=\frac{\partial f}{\partial s}=\frac{1}{F_{min}}\left(\frac{\partial M}{\partial s}-\frac{T_{min}S_{min}}{\tau}\right),~{}~{}\phi^{\theta}=-\cot\theta\csc\theta. (41)

The zero points of the mapping ϕ\phi can be calculated by ϕs=0\phi^{s}=0. We solve it and get

τ=2s1/3(12L3Φ~2)1+2L3Φ~2+s2/3(12L3Φ~2).\displaystyle\tau=\frac{\sqrt{2}s^{1/3}(1-2L^{3}\widetilde{\Phi}^{2})}{1+2L^{3}\widetilde{\Phi}^{2}+s^{2/3}(1-2L^{3}\widetilde{\Phi}^{2})}. (42)

As s=(1+2L3Φ~2)3/2/(12L3Φ~2)3/2s=(1+2L^{3}\widetilde{\Phi}^{2})^{3/2}/(1-2L^{3}\widetilde{\Phi}^{2})^{3/2}, τ=τmax=2([1+2L3Φ~212L3Φ~2]1/3+[1+2L3Φ~212L3Φ~2]2/3)\tau=\tau_{max}=\frac{\sqrt{2}}{\left(\left[\frac{1+2L^{3}\widetilde{\Phi}^{2}}{1-2L^{3}\widetilde{\Phi}^{2}}\right]^{1/3}+\left[\frac{1+2L^{3}\widetilde{\Phi}^{2}}{1-2L^{3}\widetilde{\Phi}^{2}}\right]^{2/3}\right)}. The zero points of ϕs\phi^{s} in the diagram of τs\tau-s is displayed in Fig. 4. τmax\tau_{max} is related with the annihilation point. As τ<τmax\tau<\tau_{max}, there are two intersection points for the AdS black hole under the (5,1)(5,1)-model frame in the grand canonical ensemble. The red and blue lines are for the low-potential black hole (LPSB) and high-potential black hole (HPBH), respectively. They are just the stable and unstable black hole states. The corresponding winding numbers are 1, -1, so the topological number is W=11=0W=1-1=0 as τ<τmax\tau<\tau_{max}. For τ=τmax\tau=\tau_{max}, the two intersection points for this type AdS black hole coincide. Furthermore, there is no black hole as τ>τmax\tau>\tau_{max}.

III.3 Topology of RPST for (6,2)(6,2)-model: BI gravity coupled to Maxwell’s field in the canonical ensemble

When adopting (d,k)=(6,2)(d,k)=(6,2), the Wald entropy for the charged spherically symmetric AdS black hole solution was given in Ref. Kong2022 as

S=π8Mrh32Q2/πG(6,2)l4rh4G(6,2)l2\displaystyle S=\pi\sqrt{\frac{8Mr_{h}^{3}-2Q^{2}/\pi}{G_{(6,2)}l^{4}}}-\frac{r_{h}^{4}}{G_{(6,2)}l^{2}} (43)

With the definitions S¯=L2S/N+πl2,Q¯=Q/N\bar{S}=L^{2}S/N+\pi l^{2},~{}\bar{Q}=Q/N, the corresponding radiation temperature and mass parameter are

T\displaystyle T =\displaystyle= 5S¯214π1/2lS¯3/2+13πl2S¯4π3/2l3S¯1/2π2L2Q¯28π5/4l3/2S¯1/2(S¯1/2π1/2l)5/2,\displaystyle\frac{5\bar{S}^{2}-14\pi^{1/2}l\bar{S}^{3/2}+13\pi l^{2}\bar{S}-4\pi^{3/2}l^{3}\bar{S}^{1/2}-\pi^{2}L^{2}\bar{Q}^{2}}{8\pi^{5/4}l^{3/2}\bar{S}^{1/2}\left(\bar{S}^{1/2}-\pi^{1/2}l\right)^{5/2}}, (44)
M\displaystyle M =\displaystyle= N3S¯26π1/2lS¯3/2+3πl2S¯+π2L2Q¯26π5/4l3/2L2(S¯1/2π1/2l)3/2.\displaystyle N\frac{3\bar{S}^{2}-6\pi^{1/2}l\bar{S}^{3/2}+3\pi l^{2}\bar{S}+\pi^{2}L^{2}\bar{Q}^{2}}{6\pi^{5/4}l^{3/2}L^{2}\left(\bar{S}^{1/2}-\pi^{1/2}l\right)^{3/2}}. (45)

From these thermodynamical parameters, the critical point reads

S¯c\displaystyle\bar{S}_{c} =\displaystyle= (0.9058+π)l2,Q¯=0.0358l2/L,Tc=0.1747/l,\displaystyle(0.9058+\pi)l^{2},~{}\bar{Q}=0.0358l^{2}/L,~{}T_{c}=0.1747/l, (46)
Φ~c\displaystyle\widetilde{\Phi}_{c} =\displaystyle= 0.2406/(lL),Mc=0.241lN/L2,Fc=0.0828lN/L2.\displaystyle 0.2406/(lL),~{}M_{c}=0.241lN/L^{2},~{}F_{c}=0.0828lN/L^{2}. (47)
Refer to caption
Figure 5: Zero points of ϕs\phi^{s} in the diagram of τs\tau-s with q=0.4q=0.4 for the higher dimensions and higher order curvature gravity model with (d,k)=(6,2)(d,k)=(6,2) in the canonical ensemble.

By introducing the relative parameters f=F/Fc=(MTS)/Fc,1/τ=t=T/Tc,s=S¯/S¯c,q=Q¯/Q¯cf=F/F_{c}=(M-TS)/F_{c},~{}1/\tau=t=T/T_{c},~{}s=\bar{S}/\bar{S}_{c},~{}q=\bar{Q}/\bar{Q}_{c}, the component of the mapping ϕs\phi^{s} reads

ϕs=s3(0.9058+π)2s26π(0.9058+π)3/2s3/2+3π(0.9058+π)s+(0.0358qπ)20.4968π5/4[(0.9058+π)sπ]3/21.91115τ.\displaystyle\phi^{s}=\frac{\partial}{\partial s}\frac{3(0.9058+\pi)^{2}s^{2}-6\sqrt{\pi}(0.9058+\pi)^{3/2}s^{3/2}+3\pi(0.9058+\pi)s+(0.0358q\pi)^{2}}{0.4968\pi^{5/4}\left[\sqrt{(0.9058+\pi)s}-\sqrt{\pi}\right]^{3/2}}-\frac{1.91115}{\tau}. (48)

The zero points of the mapping ϕ\phi can be calculated by ϕs=0\phi^{s}=0. In the Fig. 5, we can see that there are three branches as τmin<τ<τmax\tau_{min}<\tau<\tau_{max}, the large black hole branch for τ<τmin\tau<\tau_{min}, and the small black hole branch for τ>τmax\tau>\tau_{max}. That means the existence of a phase transition in the region of τmin<τ<τmax\tau_{min}<\tau<\tau_{max}. The winding numbers of the large and small black hole branches, which are stable, both read w=1w=1. And the corresponding topological number equals to one. While for the system in the region of τmin<τ<τmax\tau_{min}<\tau<\tau_{max}, there exist the stable large and small black hole, and the unstable intermediate black hole, thus the corresponding winding number are w=1,1,1w=1,1,-1. The topological number is W=1W=1. These results are consistent with that for the RN-AdS black hole and the AdS black hole under the (5,1)(5,1)-model in the canonical ensemble.

III.4 Topology of RPST for (6,2)(6,2)-model: BI gravity coupled to Maxwell’s field in the grand canonical ensemble

In this part we will consider the topology of the charged AdS black hole for the case of (6,2)(6,2)-model in the grand canonical ensemble. From eq. (24), we have

rh2=l2(1+1+SG(6,2)πl2).\displaystyle r_{h}^{2}=l^{2}\left(-1+\sqrt{1+\frac{SG_{(6,2)}}{\pi l^{2}}}\right). (49)

With the above equation, eqs. (25), (26), and the definition S~=S¯/(πl2)=L2Sπl2N+1,Q¯=Q/N\widetilde{S}=\bar{S}/(\pi l^{2})=\frac{L^{2}S}{\pi l^{2}N}+1,~{}\bar{Q}=Q/N, the eqs. (44) and (45) can be rewritten as

M\displaystyle M =\displaystyle= Nl[S~22S~3/2+S~+3Φ~2N2l2L2(S~1)3]2L2(S~1)3/2,\displaystyle\frac{Nl\left[\widetilde{S}^{2}-2\widetilde{S}^{3/2}+\widetilde{S}+3\widetilde{\Phi}^{2}N^{2}l^{2}L^{2}(\sqrt{\widetilde{S}}-1)^{3}\right]}{2L^{2}(\sqrt{\widetilde{S}}-1)^{3/2}}, (50)
T\displaystyle T =\displaystyle= 5S~214S~3/2+13S~4S~1/29Φ~2N2l2L2(S~1)38πlS~(S~1)5/2.\displaystyle\frac{5\widetilde{S}^{2}-14\widetilde{S}^{3/2}+13\widetilde{S}-4\widetilde{S}^{1/2}-9\widetilde{\Phi}^{2}N^{2}l^{2}L^{2}(\sqrt{\widetilde{S}}-1)^{3}}{8\pi l\sqrt{\widetilde{S}}(\sqrt{\widetilde{S}}-1)^{5/2}}. (51)
Refer to caption
Figure 6: Zero points of ϕS~\phi^{\widetilde{S}} in the diagram of τS~\tau-\widetilde{S} with Φ~2N2l2L2=0.8\widetilde{\Phi}^{2}N^{2}l^{2}L^{2}=0.8 for the AdS black hole under the (6,2)-model frame in the grand canonical ensemble.

Through the redefinitions f(S~)=FŁ2/(lN),m(S~)=MŁ2/(lN),1/τ=t(S~)=TlNf(\widetilde{S})=F\L^{2}/(lN),~{}m(\widetilde{S})=M\L^{2}/(lN),~{}1/\tau=t(\widetilde{S})=TlN, the components of the mapping ϕ\phi read

ϕS~=f(S~)S~=m(S~)S~1τ,ϕθ=cotθcscθ,\displaystyle\phi^{\widetilde{S}}=\frac{\partial f(\widetilde{S})}{\partial\widetilde{S}}=\frac{\partial m(\widetilde{S})}{\partial\widetilde{S}}-\frac{1}{\tau},~{}~{}\phi^{\theta}=-\cot\theta\csc\theta, (52)

From the zero-point equation of the mapping ϕS~=0\phi^{\widetilde{S}}=0, we can obtain

1τ=5S~+S~(9Φ~2N2l2L24)9Φ~2N2l2L28S~S~1.\displaystyle\frac{1}{\tau}=\frac{5\widetilde{S}+\sqrt{\widetilde{S}}(9\widetilde{\Phi}^{2}N^{2}l^{2}L^{2}-4)-9\widetilde{\Phi}^{2}N^{2}l^{2}L^{2}}{8\sqrt{\widetilde{S}}\sqrt{\sqrt{\widetilde{S}}-1}}. (53)

The zero points of ϕS~\phi^{\widetilde{S}} in the diagram of τS~\tau-\widetilde{S} is displayed in Fig. 6. τmax\tau_{max} is related with the annihilation point. As τ<τmax\tau<\tau_{max}, there are two intersection points for the AdS black hole under the (6,2)(6,2)-model frame in the grand canonical ensemble. The red and blue lines are for the low-potential black hole (LPSB) and high-potential black hole (HPBH), respectively. They are just the stable and unstable black hole states. The corresponding winding numbers are 1, -1, so the topological number is W=11=0W=1-1=0 as τ<τmax\tau<\tau_{max}. For τ=τmax\tau=\tau_{max}, the two intersection points for this type AdS black hole coincide. Furthermore, there is no black hole as τ>τmax\tau>\tau_{max}.

IV Discussions and conclusions

In this work, under the restricted phase space frame we investigated the thermodynamical topology of the RN-AdS black hole as well as two typical models of (d,k)=(5,1),(6,2)(d,k)=(5,1),~{}(6,2) in the higher dimensions and higher order curvature gravity. The results showed that the topology numbers for these black holes in the canonical and grand canonical ensembles are different. For the charged AdS black holes (both the RN-AdS black hole and two others with (d,k)=(5,1),(6,2)(d,k)=(5,1),~{}(6,2)) in the canonical ensemble, the topological number is positive one, while it is zero for the charged AdS black holes in the grand canonical ensemble. This also proves that the EH and BI like models belong to the same gravitational classification from the thermodynamical topology perspective.

Acknowledgements

We would like to thank Prof. Ren Zhao for his indispensable discussions and comments. This work was supported by the National Natural Science Foundation of China (Grant No. 12075143, Grant No. 12375050), and the Natural Science Foundation of Shanxi Province, China (Grant No. 202203021221209, Grant No. 202303021211180).

References

  • (1) J. M. Bardeen, B. Carter, and S. W. Hawking, Commun. Math. Phys. 31, 161 (1973).
  • (2) J. D. Bekenstein, Phys. Rev. D 7, 949 (1973).
  • (3) S. W. Hawking and Don N. Page, Commun. Math. Phys. 87, 577-588 (1983).
  • (4) E. Witten, Adv. Theor. Math. Phys. 2, 505 (1998), arXiv:hep-th/9803131.
  • (5) D. Kastor, S. Ray, and J. Traschen, Class. Quant. Grav. 26 (2009) 195011, arXiv:0904.2765.
  • (6) R. A. Hennigar, E. Tjoa, and R. B. Mann, J. High Energ. Phys (2017) 70, arXiv:1612.06852.
  • (7) R. A. Hennigar and R. B. Mann, Phys. Rev. Lett. 118, 021301 (2017), arXiv:1609.02564.
  • (8) A. M. Frassino, D. Kubiznak, R. B. Mann, et al., J. High Energ. Phys, 09 (2014) 080, arXiv:1406.7015.
  • (9) D. Kubiznak and R. B. Mann, J. High Energ. Phys 1207 (2012) 033, arXiv:1205.0559.
  • (10) R.-G Cai, L.-M Cao, L. Li, et al., J. High Energ. Phys, 2013 (9), arXiv:1306.6233.
  • (11) M.-S Ma, R. Zhao, and Y.-S Liu, Class. Quant. Grav. 34 (2017) 165009, arXiv:1604.06998.
  • (12) R. Banerjee, B. R. Majhi, and S. Samanta, Phys. Lett. B 767 (2017) 25-28, arXiv:1611.06701.
  • (13) D. Kubiznak and R. B. Mann, J. High Energ. Phys 1207 (2012) 033, arXiv:1205.0559.
  • (14) S.-W Wei and Y.-X Liu, Phys. Rev. Lett 115 (2015) 111302.
  • (15) K. Bhattacharya, B. R. Majhi, and S. Samanta, Phys. Rev. D 96, 084037 (2017).
  • (16) X.-X Zeng and L.-F Li, Phys. Lett. B 764 (2017) 100.
  • (17) J.-L Zhang, R.-G Cai, and H.-W Yu,
  • (18) Y.-Z Du, H.-F Li, F. Liu, et al., Chin. Phys. C 45 (2021) 11, arXiv:2112.10403.
  • (19) Y. Zhang, W.-Q Wang, Y.-B Ma, et al., Adv. H. E. Phys. 2020 (2020) 7263059, arXiv:2004.06796.
  • (20) M. R. Visser, Phys. Rev. D 105, 106014 (2022).
  • (21) M. B. Ahmed, W. Cong, D. Kubiznak, et al., Phys. Rev. Lett. 130, 181401 (2023).
  • (22) L. Zhao, Chin. Phys. C 46, 055105 (2022).
  • (23) Z.-Y Gao, X. Kong, and L. Zhao, Eur. Phys. J. C 82, 112 (2022).
  • (24) Z. Y. Gao and L. Zhao, Class. Quant. Grav. 39, 075019 (2022).
  • (25) Y.-Z Du, H.-F Li, Y. Zhang, et al., Entropy 25 (2023) 4, arXiv:2210.02006.
  • (26) Y.-S Duan and M. L. Ge, Sci. Sin. 9, 1072 (1979).
  • (27) S.-W Wei, Y.-X Liu, and R. B. Mann, Phys. Rev. L 129, 191101 (2022).
  • (28) X.-Q Kong, T Wang, Z.-Y Gao, et al., Entropy 24 (2022) 8, arXiv:2208.07748.
  • (29) R.-G Cai, Phys. Rev. D, 2002, 65 (8):084014, arXiv:hepth/0109133.
  • (30) C. Lanczos,Z. Phys., 1932, 73(3):147-168.
  • (31) D. Lovelock, J. Math. Phys., 1971, 12(3):498-501.
  • (32) S.-W. Wei, Phys. Rev. D 102, 064039 (2020).
  • (33) S.-W. Wei and Y.-X. Liu, Phys. Rev. D 105, 104003 (2022).
  • (34) D Wu, Eur. Phys. J. C 83 (2023) 589,arXiv:2306.02324[gr-qc]
  • (35) D Wu, Eur. Phys. J. C 83 (2023) 365, arXiv:2302.01100.
  • (36) Y.-S. Duan, The structure of the topological current, Report No. SLAC-PUB-3301 (1984).
  • (37) Y.-S. Duan, S. Li, and G.-H. Yang, Nucl. Phys. B514, 705 (1998).
  • (38) L.-B. Fu, Y.-S. Duan, and H. Zhang, Phys. Rev. D 61, 045004 (2000).
  • (39) Z.-Y. Fan, Phys. Rev. D 107, 044026 (2023).
  • (40) N. J. Gogoi and P. Phukon, Phys. Rev. D 108 (2023) 6, arXiv:2304.05695.
  • (41) J. Crisostomo, R. Troncoso, and J. Zanelli, Phys. Rev. D, 62 (2000) 8, arXiv:hep-th/0003271.
  • (42) M. Banados, C. Teitelboim, and J. Zanelli, Phys. Rev. Lett. 69 (1992) 1849, arXiv:hep-th/9204099.
  • (43) M. Banados, M. Henneaux, C. Teitelboim, et al., Phys. Rev. D 48 (1993) 1506, arXiv:gr-qc/9302012.
  • (44) M. Banados, C. Teitelboim, and J. Zanelli, Phys. Rev. D 49 (1994) 975.