HTML conversions sometimes display errors due to content that did not convert correctly from the source. This paper uses the following packages that are not yet supported by the HTML conversion tool. Feedback on these issues are not necessary; they are known and are being worked on.

  • failed: cmap

Authors: achieve the best HTML results from your LaTeX submissions by following these best practices.

License: CC Zero
arXiv:2403.14730v1 [gr-qc] 21 Mar 2024

Topology of dyonic AdS black holes with quasitopological electromagnetism in Einstein-Gauss-Bonnet gravity

Hao Chen haochen1249@yeah.net School of Physics and Electronic Science, Zunyi Normal University, Zunyi 563006,PR China    Meng-Yao Zhang gs.myzhang21@gzu.edu.cn School of Mathematics and Statistics, Guizhou University, Guiyang, 550025, China    Hassan Hassanabadi h.hasanabadi@shahroodut.ac.ir Faculty of Physics, Shahrood University of Technology, Shahrood, Iran    Bekir Can Lütfüoğlu bekir.lutfuoglu@uhk.cz Department of Physics, University of Hradec Králové,
Rokitanského 62, 500 03 Hradec Králové, Czechia.
   Zheng-Wen Long zwlong@gzu.edu.cn College of Physics, Guizhou University, Guiyang, 550025, China
(March 21, 2024)
Abstract

Abstract: In this study, we employ the thermodynamic topological method to classify critical points for the dyonic AdS black holes with QTE in the EGB background. To this end, we find that there is a small/large BH phase transition in any space-time dimension, a conventional critical point exists with the total topological charge of Qt=1subscript𝑄𝑡1Q_{t}=-1italic_Q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = - 1. The existence of the coupling constant α𝛼\alphaitalic_α gives rise to a more intricate phase structure of the black hole, with the emergence of a triple points requires α0.5𝛼0.5\alpha\geq 0.5italic_α ≥ 0.5 and d=6𝑑6d=6italic_d = 6. Interestingly, the condition for the simultaneous occurrence of small/intermediate and intermediate/large phase transition is that the coupling constant a takes a special value (α=0.5𝛼0.5\alpha=0.5italic_α = 0.5), the two conventional critical points (CP1,CP2)𝐶subscript𝑃1𝐶subscript𝑃2(CP_{1},CP_{2})( italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) of the black hole are (physical) critical point, and the novel critical point that lacks the capability to minimize the Gibbs free energy. The critical point (QCP1=QCP2=1subscript𝑄𝐶subscript𝑃1subscript𝑄𝐶subscript𝑃21Q_{CP_{1}}=Q_{CP_{2}}=-1italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - 1) is observed to occur at the maximum extreme points of temperature in the isobaric curve, while the critical point (QCP3=1)subscript𝑄𝐶subscript𝑃31(Q_{CP_{3}}=1)( italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 1 ), emerges at the minimum extreme points of temperature. Furthermore, the number of phases at the novel critical point exhibits an upward trend, followed by a subsequent decline at the conventional critical points. With the increase of the coupling constant (α=1)𝛼1(\alpha=1)( italic_α = 1 ), although the system has three critical points, only CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is a (physical) critical point, and the CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT serves as the phase annihilation point. This means that the coupling constant α𝛼\alphaitalic_α has a non-negligible effect on the phase structure of the black hole.

Black hole; Dark matter; Dark energy;Phase transition; Joule-Thomson expansion
preprint: APS/-QNM

I Introduction

The official announcement on February 11, 2016 of the detection of gravitational waves by the LIGO represents a groundbreaking achievement that provides significant validation for Einstein’s general relativity. This momentous discovery not only confirms the existence of black holes but also verifies the presence of gravitational waves ch1 ; ch2 ; ch3 ; ch4 . The emergence of black hole thermodynamics stems from the application of general laws in thermodynamics to the general relativity, with the objective of investigating the characteristics inherent to black holes. The fundamental concept involves regarding black holes as a system governed by thermodynamic principles, wherein mass corresponds to energy, entropy is directly proportional to the surface area, and temperature is represented by surface gravity. In 1974, Hawking made a groundbreaking discovery that black holes exhibit thermal radiation ch5 , this revelation not only resolves the inherent contradictions within the thermodynamics of black holes but also successfully establishes a profound connection between gravitational theory, quantum mechanics, and thermodynamics. The phenomenon of phase transition is commonly observed in thermodynamic systems, wherein the Schwarzschild black hole, when considered as a thermodynamic system within the framework of AdS, exhibits an unstable state from a thermodynamic perspective at low temperatures. However, with increasing temperature, the free energy is reduced, and finally the black hole reaches stability. Essentially, the stability of a black hole hinges on the critical temperature known as the Hawking-Page (HP) phase transition ch6 . Based on the AdS/CFT correspondence theory, there has been considerable focus on studying the HP phase transition, which is recognized as a phase transition between confinement and deconfinement in the gauge/gravity duality ch7 .

Kastor and Ray proposed the concept of regarding the cosmological constant ΛΛ\Lambdaroman_Λ as a variable thermodynamic pressure within the black hole, and suggested that the volume of the black hole in AdS space can be considered as its conjugate quantity, while the black hole’s mass is interpreted as the spacetime’s enthalpy, rather than its internal energy ch8 . This notion was further expanded to Lovelock theory and applied to black holes, which demonstrates that Lovelock coupling serves as an extension of the first law of thermodynamics ch9 . Additionally, both the free energy and of ADM mass AdS black holes were derived ch10 . Furthermore, investigations on dilaton-AdS black holes have revealed that their mass represents the system’s enthalpy ch11 , the phenomenon has led to a surge in research on thermodynamic phase transitions in AdS black holes. For instance, the charged RN-AdS black hole in four-dimensional spacetime was examined by Kubiznak and Mann, who identified a significant transition between small and large black holes, which bears resemblance to the gas-liquid phase transition observed in van der Waals fluids ch17 . The phase transition in nonlinear electromagnetic field was investigated by Hendi and Vahidinia through considering the Einstein gravitational solution under the background of generalized Maxwell theory ch18 . Critical phenomena and critical exponents near the critical point of Gauss-Bonnet black hole were studied by Wei et al, who also examined the phenomenon of phase transition in spin Kerr black holes with a high number of dimensions ch19 ; ch20 . Zou et al. explored P-V critical behavior of Born-Infeld Anti-deSitter black hole and Gauss-Bonnet AdS black hole ch21 ; ch22 . It is well known that the typical first-order phase transition has only one critical point. However, the existence of multiple critical points is observed in certain categories of black holes, such as, the occurrence of reentrant phase transition ch24 ; ch25 , the presence of triple points in dynamic AdS black holes under quasi-topological electromagnetism (QTE) ch26 ; ch27 , and the λ𝜆\lambdaitalic_λ phase transition observed in Horava gravity ch28 , among others. A complementary aspect to the investigation of black hole criticality is the concept of thermodynamic topology. The Duan ’s φ𝜑\varphiitalic_φ-mapping theory proposed is utilized to explore the criticality of black holes, resulting in distinct characteristics of topological charge exhibited by critical points within the thermodynamic space, these critical points can be classified into conventional and novel types cq4 . Furthermore, Wei et al. employed the thermodynamic topological method to investigate thermodynamic properties of black holes, this approach enables analysis of stability by considering the winding number associated with defect we1 . The thermodynamic topology classification of different black holes are found in detail in we2 ; we3 ; we4 ; we5 ; we7 ; we9 ; we8 ; we10 ; we12 ; we13 ; we14 ; we15 ; we16 ; we17 ; we18 .

The theory of extended gravity, which encompasses EGB gravity, provides novel perspectives on the characteristics of black holes ch29 ; ch30 ; ch31 ; ch32 ; ch33 ; ch34 ; ch35 . Through an examination of the thermodynamic characteristics of the Born-Infeld AdS black hole within the framework of 4D Einstein-Gauss-Bonnet (EGB) gravity, it has been observed that there exist three critical points within the black hole, these critical points exhibit numerous intriguing properties when subjected to EGB gravity ch25 . Recently, the QTE is characterized by the square norm of the topological 4-form FF𝐹𝐹F\wedge Fitalic_F ∧ italic_F, indicating that the system’s energy-momentum tensor corresponds to a perfect fluid with isotropic properties ch36 . Based on this, the thermodynamic phase transition of dyonic AdS black holes with QTE is discussed in detail. By changing the coupling parameters, a triple point phase structure is observed ch27 . Additionally, through an examination of the interplay between QTE and EGB gravity, Sekhmani et al. derived a solution for dyonic AdS black holes with a high number of dimensions. Furthermore, they conducted a comprehensive analysis on how the coupling constant α𝛼\alphaitalic_α influences both the energy emission rates and black hole shadow ch37 . On this basis, a comprehensive analysis is conducted on the thermodynamic phase transition of the dyonic AdS black holes with QTE in the EGB background, and it is found that the system parameters will appear triple points in special regions cq3 . In this work, we will use the thermodynamic topological method to investigate the critical points of black holes, and further understand the influence of GB coupling constant on the phase structure of dyonic AdS black holes. The structure of the manuscript is as follows: In section 2, we will give a concise overview of the thermodynamic properties associated with dyonic AdS black holes with QTE in the EGB background and present a comprehensive explanation of Duan’s φ𝜑\varphiitalic_φ-mapping theory. In section 3, we will classify the critical points of black holes in different space-time dimensions. In section 4, we uncover the nature of the critical point by the behavior of the isobaric curve. Section 5 contains the remarks and conclusions.

II Thermodynamics of dyonic AdS black holes and topology

We first consider that the d-dimensional action involving EGB gravity and QTE can be formulated in a minimal coupling manner ch37

SD=116πddxg(2Λ+α𝒢+QE),subscript𝑆𝐷116𝜋superscript𝑑𝑑𝑥𝑔2Λ𝛼𝒢subscriptQES_{D}=\frac{1}{16\pi}\int d^{d}x\sqrt{-g}\left(\mathcal{R}-2\Lambda+\alpha% \mathcal{G}+\mathcal{L}_{\mathrm{QE}}\right),italic_S start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 16 italic_π end_ARG ∫ italic_d start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT italic_x square-root start_ARG - italic_g end_ARG ( caligraphic_R - 2 roman_Λ + italic_α caligraphic_G + caligraphic_L start_POSTSUBSCRIPT roman_QE end_POSTSUBSCRIPT ) , (1)

wherein, Λ=(d1)(d2)2l2Λ𝑑1𝑑22superscript𝑙2\Lambda=-\frac{(d-1)(d-2)}{2l^{2}}roman_Λ = - divide start_ARG ( italic_d - 1 ) ( italic_d - 2 ) end_ARG start_ARG 2 italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG represents the cosmological constant within the AdS space, α𝛼\alphaitalic_α is the symbol used to denote the GB coupling constant, which possesses dimensions of [length]22{}^{2}start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT, the GB term can be expressed as

𝒢=24μνμν+μνρσμνρσ.𝒢superscript24superscript𝜇𝜈subscript𝜇𝜈superscript𝜇𝜈𝜌𝜎subscript𝜇𝜈𝜌𝜎\mathcal{G}=\mathcal{R}^{2}-4\mathcal{R}^{\mu\nu}\mathcal{R}_{\mu\nu}+\mathcal% {R}^{\mu\nu\rho\sigma}\mathcal{R}_{\mu\nu\rho\sigma}.caligraphic_G = caligraphic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 caligraphic_R start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT caligraphic_R start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT + caligraphic_R start_POSTSUPERSCRIPT italic_μ italic_ν italic_ρ italic_σ end_POSTSUPERSCRIPT caligraphic_R start_POSTSUBSCRIPT italic_μ italic_ν italic_ρ italic_σ end_POSTSUBSCRIPT . (2)

The matter field, a fundamental concept in the realm of physics, can be characterized by an electromagnetic Lagrangian with quasi-topological properties and represented as

QE=(142+12p!2+βint),subscriptQE14superscript212𝑝superscript2𝛽subscriptint\displaystyle\mathcal{L}_{\rm{QE}}=-(\frac{1}{4}\mathcal{F}^{2}+\frac{1}{2p!}% \mathcal{H}^{2}+\beta\mathcal{L}_{\rm{int}}),caligraphic_L start_POSTSUBSCRIPT roman_QE end_POSTSUBSCRIPT = - ( divide start_ARG 1 end_ARG start_ARG 4 end_ARG caligraphic_F start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 italic_p ! end_ARG caligraphic_H start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_β caligraphic_L start_POSTSUBSCRIPT roman_int end_POSTSUBSCRIPT ) , (3)

where, p=d2𝑝𝑑2p=d-2italic_p = italic_d - 2, 2=μνμνsuperscript2subscript𝜇𝜈superscript𝜇𝜈\mathcal{F}^{2}=\mathcal{F}_{\mu\nu}\mathcal{F}^{\mu\nu}caligraphic_F start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = caligraphic_F start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT caligraphic_F start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT and 2=ρ1ρpρ1ρpsuperscript2subscriptsubscript𝜌1subscript𝜌𝑝superscriptsubscript𝜌1subscript𝜌𝑝\mathcal{H}^{2}=\mathcal{H}_{\rho_{1}\ldots\rho_{p}}\mathcal{H}^{\rho_{1}% \ldots\rho_{p}}caligraphic_H start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = caligraphic_H start_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT … italic_ρ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT end_POSTSUBSCRIPT caligraphic_H start_POSTSUPERSCRIPT italic_ρ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT … italic_ρ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT end_POSTSUPERSCRIPT, and β𝛽\betaitalic_β means the coupling constant associated with the dimension [length]22{}^{2}start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT. Additionally, the interaction term reads

int=δγ1γDλ1λDλ1λ2λ3λDγ1γ2γ3γD.subscriptintsuperscriptsubscript𝛿subscript𝛾1subscript𝛾𝐷subscript𝜆1subscript𝜆𝐷subscriptsubscript𝜆1subscript𝜆2subscriptsubscript𝜆3subscript𝜆𝐷superscriptsubscript𝛾1subscript𝛾2superscriptsubscript𝛾3subscript𝛾𝐷\displaystyle\mathcal{L}_{\rm{int}}=\delta_{\gamma_{1}\ldots\gamma_{D}}^{% \lambda_{1}\ldots\lambda_{D}}\mathcal{F}_{\lambda_{1}\lambda_{2}}\mathcal{H}_{% \lambda_{3}\ldots\lambda_{D}}\mathcal{F}^{\gamma_{1}\gamma_{2}}\mathcal{H}^{% \gamma_{3}\ldots\gamma_{D}}.caligraphic_L start_POSTSUBSCRIPT roman_int end_POSTSUBSCRIPT = italic_δ start_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT … italic_γ start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT … italic_λ start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT end_POSTSUPERSCRIPT caligraphic_F start_POSTSUBSCRIPT italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT caligraphic_H start_POSTSUBSCRIPT italic_λ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT … italic_λ start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT end_POSTSUBSCRIPT caligraphic_F start_POSTSUPERSCRIPT italic_γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT caligraphic_H start_POSTSUPERSCRIPT italic_γ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT … italic_γ start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT end_POSTSUPERSCRIPT . (4)

The above action, referring to a specific mathematical calculation or physical process, leads to the emergence of a BH solution, which can be described by a metric

ds2=f(r)dt2+f(r)1dr2+r2dΩd22,𝑑superscript𝑠2𝑓𝑟𝑑superscript𝑡2𝑓superscript𝑟1𝑑superscript𝑟2superscript𝑟2𝑑superscriptsubscriptΩ𝑑22\displaystyle ds^{2}=-f(r)dt^{2}+f(r)^{-1}dr^{2}+r^{2}d\Omega_{d-2}^{2},\ italic_d italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = - italic_f ( italic_r ) italic_d italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_f ( italic_r ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_d italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d roman_Ω start_POSTSUBSCRIPT italic_d - 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (5)

where dΩd22𝑑superscriptsubscriptΩ𝑑22d\Omega_{d-2}^{2}italic_d roman_Ω start_POSTSUBSCRIPT italic_d - 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT represents the unit sphere in d2𝑑2d-2italic_d - 2) dimensions, The metric function is provided by

f(r)=1+r22α~(11+4α~mrd1+8α~Λ(d2)(d1)ϱ)𝑓𝑟1superscript𝑟22~𝛼114~𝛼𝑚superscript𝑟𝑑18~𝛼Λ𝑑2𝑑1italic-ϱf(r)=1+\frac{r^{2}}{2\tilde{\alpha}}\left(1-\sqrt{1+\frac{4\tilde{\alpha}m}{r^% {d-1}}+\frac{8\tilde{\alpha}\Lambda}{(d-2)(d-1)}-\varrho}\right)italic_f ( italic_r ) = 1 + divide start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 over~ start_ARG italic_α end_ARG end_ARG ( 1 - square-root start_ARG 1 + divide start_ARG 4 over~ start_ARG italic_α end_ARG italic_m end_ARG start_ARG italic_r start_POSTSUPERSCRIPT italic_d - 1 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 8 over~ start_ARG italic_α end_ARG roman_Λ end_ARG start_ARG ( italic_d - 2 ) ( italic_d - 1 ) end_ARG - italic_ϱ end_ARG ) (6)

with

ϱ=2α~d3(Qm2+Qe2ξ),italic-ϱ2~𝛼𝑑3superscriptsubscript𝑄𝑚2superscriptsubscript𝑄𝑒2𝜉\varrho=\frac{2\tilde{\alpha}}{d-3}\left(Q_{m}^{2}+Q_{e}^{2}\xi\right),italic_ϱ = divide start_ARG 2 over~ start_ARG italic_α end_ARG end_ARG start_ARG italic_d - 3 end_ARG ( italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ξ ) , (7)

the hypergeometric function, denoted by ξ𝜉\xiitalic_ξ, can be expressed as

ξ=2F1[1,d32(d2);73d42d;8βQm2Γ(d1)2r2d4]subscript2𝜉subscript𝐹11𝑑32𝑑273𝑑42𝑑8𝛽superscriptsubscript𝑄𝑚2Γsuperscript𝑑12superscript𝑟2𝑑4\xi=\ _{2}F_{1}\left[1,\frac{d-3}{2(d-2)};\frac{7-3d}{4-2d};\frac{-8\beta Q_{m% }^{2}\Gamma(d-1)^{2}}{r^{2d-4}}\right]italic_ξ = start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ 1 , divide start_ARG italic_d - 3 end_ARG start_ARG 2 ( italic_d - 2 ) end_ARG ; divide start_ARG 7 - 3 italic_d end_ARG start_ARG 4 - 2 italic_d end_ARG ; divide start_ARG - 8 italic_β italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_Γ ( italic_d - 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 italic_d - 4 end_POSTSUPERSCRIPT end_ARG ] (8)
α~=(d4)(d3)α~𝛼𝑑4𝑑3𝛼\widetilde{\alpha}=(d-4)(d-3)\alphaover~ start_ARG italic_α end_ARG = ( italic_d - 4 ) ( italic_d - 3 ) italic_α (9)

here,Qesubscript𝑄𝑒Q_{e}italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT and Qmsubscript𝑄𝑚Q_{m}italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT are essential parameters that indicate the electrical and magnetic charges of black holes, respectively. These parameters are critical in comprehending the characteristics of black holes. m𝑚mitalic_m denotes the mass of the solution within a designated parameter domain, serving as the constant of integration derived from the boundary conditions ch37 . In this case, Arnowitt-Deser-Misner mass is given by

M=(d2)𝒱d216πm𝑀𝑑2subscript𝒱𝑑216𝜋𝑚\displaystyle M=\frac{(d-2)\mathcal{V}_{d-2}}{16\pi}mitalic_M = divide start_ARG ( italic_d - 2 ) caligraphic_V start_POSTSUBSCRIPT italic_d - 2 end_POSTSUBSCRIPT end_ARG start_ARG 16 italic_π end_ARG italic_m (10)

with

𝒱d2=2πd12Γ(d12),subscript𝒱𝑑22superscript𝜋𝑑12Γ𝑑12\mathcal{V}_{d-2}=\frac{2\pi^{\frac{d-1}{2}}}{\Gamma\left(\frac{d-1}{2}\right)},caligraphic_V start_POSTSUBSCRIPT italic_d - 2 end_POSTSUBSCRIPT = divide start_ARG 2 italic_π start_POSTSUPERSCRIPT divide start_ARG italic_d - 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG start_ARG roman_Γ ( divide start_ARG italic_d - 1 end_ARG start_ARG 2 end_ARG ) end_ARG , (11)

the symbol 𝒱d2subscript𝒱𝑑2\mathcal{V}_{d-2}caligraphic_V start_POSTSUBSCRIPT italic_d - 2 end_POSTSUBSCRIPT represents the volume of a unit sphere in (d2)𝑑2(d-2)( italic_d - 2 ) dimensions. The concept of thermodynamic pressure is closely related to the cosmological constant ΛΛ\Lambdaroman_Λ, specifically, P=Λ8π𝑃Λ8𝜋P=-\frac{\Lambda}{8\pi}italic_P = - divide start_ARG roman_Λ end_ARG start_ARG 8 italic_π end_ARGch8 . Based on this, the event horizon’s radius can be utilized to ascertain thermodynamic properties of black holes, including Hawking temperature, entropy, volume, and Gibbs free energy

T𝑇\displaystyle Titalic_T =1[8(d2)πr++16πr+1𝒩1]𝒩2[(12+12d(d5)\displaystyle=\frac{1}{\left[8(d-2)\pi r_{+}+16\pi r_{+}^{-1}\mathcal{N}_{1}% \right]\mathcal{N}_{2}}[(12+12d(d-5)= divide start_ARG 1 end_ARG start_ARG [ 8 ( italic_d - 2 ) italic_π italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT + 16 italic_π italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT caligraphic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ] caligraphic_N start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG [ ( 12 + 12 italic_d ( italic_d - 5 ) (12)
+32πPr+2+2d(d5)r+2𝒩1r+62dQm2)𝒩2r+6Qe2]\displaystyle\left.\left.+32\pi Pr_{+}^{2}+2d(d-5)r_{+}^{-2}\mathcal{N}_{1}-r_% {+}^{6-2d}Q_{m}^{2}\right)\mathcal{N}_{2}-r_{+}^{6}Q_{e}^{2}\right]+ 32 italic_π italic_P italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_d ( italic_d - 5 ) italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT caligraphic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 6 - 2 italic_d end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) caligraphic_N start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ]

with

𝒩1=(d4)(d3)(d2)αsubscript𝒩1𝑑4𝑑3𝑑2𝛼\mathcal{N}_{1}=(d-4)(d-3)(d-2)\alphacaligraphic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = ( italic_d - 4 ) ( italic_d - 3 ) ( italic_d - 2 ) italic_α (13)
𝒩2=r+2d+8r+4βQm2Γ(d1)2,subscript𝒩2superscriptsubscript𝑟2𝑑8superscriptsubscript𝑟4𝛽superscriptsubscript𝑄𝑚2Γsuperscript𝑑12\mathcal{N}_{2}=r_{+}^{2d}+8r_{+}^{4}\beta Q_{m}^{2}\Gamma(d-1)^{2},caligraphic_N start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_d end_POSTSUPERSCRIPT + 8 italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_β italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_Γ ( italic_d - 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (14)
S𝑆\displaystyle Sitalic_S =0rhT1(Hr)P,Qe,Qm𝑑rabsentsuperscriptsubscript0subscript𝑟superscript𝑇1subscript𝐻𝑟𝑃subscript𝑄𝑒subscript𝑄𝑚differential-d𝑟\displaystyle=\int_{0}^{r_{h}}T^{-1}\left(\frac{\partial H}{\partial r}\right)% _{P,Q_{e},Q_{m}}dr= ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( divide start_ARG ∂ italic_H end_ARG start_ARG ∂ italic_r end_ARG ) start_POSTSUBSCRIPT italic_P , italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT , italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_d italic_r (15)
=πd12r+d2[1+2(d3)(d2)αr+2]2Γ(d12),absentsuperscript𝜋𝑑12superscriptsubscript𝑟𝑑2delimited-[]12𝑑3𝑑2𝛼superscriptsubscript𝑟22Γ𝑑12\displaystyle=\frac{\pi^{\frac{d-1}{2}}r_{+}^{d-2}\left[1+2(d-3)(d-2)\alpha r_% {+}^{-2}\right]}{2\Gamma\left(\frac{d-1}{2}\right)},= divide start_ARG italic_π start_POSTSUPERSCRIPT divide start_ARG italic_d - 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d - 2 end_POSTSUPERSCRIPT [ 1 + 2 ( italic_d - 3 ) ( italic_d - 2 ) italic_α italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT ] end_ARG start_ARG 2 roman_Γ ( divide start_ARG italic_d - 1 end_ARG start_ARG 2 end_ARG ) end_ARG ,
V=(HP)S,Qe,Qm=2πd12r+d1(d1)Γ(d12),𝑉subscript𝐻𝑃𝑆subscript𝑄𝑒subscript𝑄𝑚2superscript𝜋𝑑12superscriptsubscript𝑟𝑑1𝑑1Γ𝑑12\displaystyle V=\left(\frac{\partial H}{\partial P}\right)_{S,Q_{e},Q_{m}}=% \frac{2\pi^{\frac{d-1}{2}}r_{+}^{d-1}}{(d-1)\Gamma\left(\frac{d-1}{2}\right)},italic_V = ( divide start_ARG ∂ italic_H end_ARG start_ARG ∂ italic_P end_ARG ) start_POSTSUBSCRIPT italic_S , italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT , italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT end_POSTSUBSCRIPT = divide start_ARG 2 italic_π start_POSTSUPERSCRIPT divide start_ARG italic_d - 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d - 1 end_POSTSUPERSCRIPT end_ARG start_ARG ( italic_d - 1 ) roman_Γ ( divide start_ARG italic_d - 1 end_ARG start_ARG 2 end_ARG ) end_ARG , (16)
G𝐺\displaystyle Gitalic_G =2(d2)r+2d+2𝒩3+2r+2d𝒩1𝒩3+32πP𝒩3r+2d+4d1absent2𝑑2superscriptsubscript𝑟2𝑑2subscript𝒩32superscriptsubscript𝑟2𝑑subscript𝒩1subscript𝒩332𝜋𝑃subscript𝒩3superscriptsubscript𝑟2𝑑4𝑑1\displaystyle=2(d-2)r_{+}^{2d+2}\mathcal{N}_{3}+2r_{+}^{2d}\mathcal{N}_{1}% \mathcal{N}_{3}+\frac{32\pi P\mathcal{N}_{3}r_{+}^{2d+4}}{d-1}= 2 ( italic_d - 2 ) italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_d + 2 end_POSTSUPERSCRIPT caligraphic_N start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT + 2 italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_d end_POSTSUPERSCRIPT caligraphic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT caligraphic_N start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT + divide start_ARG 32 italic_π italic_P caligraphic_N start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_d + 4 end_POSTSUPERSCRIPT end_ARG start_ARG italic_d - 1 end_ARG (17)
𝒩3((d2)r+2+2𝒩1)𝒩2[2(d3)(d2)α+r+2]subscript𝒩3𝑑2superscriptsubscript𝑟22subscript𝒩1subscript𝒩2delimited-[]2𝑑3𝑑2𝛼superscriptsubscript𝑟2\displaystyle-\frac{\mathcal{N}_{3}}{\left((d-2)r_{+}^{2}+2\mathcal{N}_{1}% \right)\mathcal{N}_{2}}\left[2(d-3)(d-2)\alpha+r_{+}^{2}\right]- divide start_ARG caligraphic_N start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_ARG start_ARG ( ( italic_d - 2 ) italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 caligraphic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) caligraphic_N start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG [ 2 ( italic_d - 3 ) ( italic_d - 2 ) italic_α + italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ]
×[(2𝒩4r+2dr+8Qm2)𝒩2r+2d+8Qe2]absentdelimited-[]2subscript𝒩4superscriptsubscript𝑟2𝑑superscriptsubscript𝑟8superscriptsubscript𝑄𝑚2subscript𝒩2superscriptsubscript𝑟2𝑑8superscriptsubscript𝑄𝑒2\displaystyle\times\left[\left(2\mathcal{N}_{4}r_{+}^{2d}-r_{+}^{8}Q_{m}^{2}% \right)\mathcal{N}_{2}-r_{+}^{2d+8}Q_{e}^{2}\right]× [ ( 2 caligraphic_N start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_d end_POSTSUPERSCRIPT - italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) caligraphic_N start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_d + 8 end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ]
+r+8(Qm2+Qe2ξr+)𝒩3d3superscriptsubscript𝑟8superscriptsubscript𝑄𝑚2superscriptsubscript𝑄𝑒2subscript𝜉subscript𝑟subscript𝒩3𝑑3\displaystyle+\frac{r_{+}^{8}\left(Q_{m}^{2}+Q_{e}^{2}\xi_{r_{+}}\right)% \mathcal{N}_{3}}{d-3}+ divide start_ARG italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT ( italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ξ start_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) caligraphic_N start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_ARG start_ARG italic_d - 3 end_ARG

with

𝒩3=πd32r+5d16Γ(d12),subscript𝒩3superscript𝜋𝑑32superscriptsubscript𝑟5𝑑16Γ𝑑12\mathcal{N}_{3}=\frac{\pi^{\frac{d-3}{2}}r_{+}^{-5-d}}{16\Gamma\left(\frac{d-1% }{2}\right)},caligraphic_N start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT = divide start_ARG italic_π start_POSTSUPERSCRIPT divide start_ARG italic_d - 3 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 5 - italic_d end_POSTSUPERSCRIPT end_ARG start_ARG 16 roman_Γ ( divide start_ARG italic_d - 1 end_ARG start_ARG 2 end_ARG ) end_ARG , (18)
𝒩4=6r+2+d(d5)r+2+16Pπr+4+(d5)𝒩1.subscript𝒩46superscriptsubscript𝑟2𝑑𝑑5superscriptsubscript𝑟216𝑃𝜋superscriptsubscript𝑟4𝑑5subscript𝒩1\mathcal{N}_{4}=6r_{+}^{2}+d(d-5)r_{+}^{2}+16P\pi r_{+}^{4}+(d-5)\mathcal{N}_{% 1}.caligraphic_N start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT = 6 italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_d ( italic_d - 5 ) italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 16 italic_P italic_π italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT + ( italic_d - 5 ) caligraphic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT . (19)

To gain a deeper comprehension of the phase structure, we will employ the thermodynamic topology approach to classify the critical points. Drawing inspiration from cq4 , which implies the allocation of topological charge to these critical points,we denote the temperature of a black hole as T=T(S,P,xi)𝑇𝑇𝑆𝑃superscript𝑥𝑖T=T\left(S,P,x^{i}\right)italic_T = italic_T ( italic_S , italic_P , italic_x start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ). Based on this premise, the determination of the critical point of the black hole can be achieved as follows:

(ST)P,xi=0,(S,ST)P,xi=0.formulae-sequencesubscriptsubscript𝑆𝑇𝑃superscript𝑥𝑖0subscriptsubscript𝑆𝑆𝑇𝑃superscript𝑥𝑖0\left(\partial_{S}T\right)_{P,x^{i}}=0,\quad\left(\partial_{S,S}T\right)_{P,x^% {i}}=0.( ∂ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT italic_T ) start_POSTSUBSCRIPT italic_P , italic_x start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT end_POSTSUBSCRIPT = 0 , ( ∂ start_POSTSUBSCRIPT italic_S , italic_S end_POSTSUBSCRIPT italic_T ) start_POSTSUBSCRIPT italic_P , italic_x start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT end_POSTSUBSCRIPT = 0 . (20)

The elimination of thermodynamic pressure leads to the derivation of a novel temperature function from Eq.(20), the next step entails the construction of a Duan’s potential cq4

Ψ=1sinθT(S,xi),Ψ1𝜃𝑇𝑆superscript𝑥𝑖\Psi=\frac{1}{\sin\theta}T\left(S,x^{i}\right),roman_Ψ = divide start_ARG 1 end_ARG start_ARG roman_sin italic_θ end_ARG italic_T ( italic_S , italic_x start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ) , (21)

note that the inclusion of 1/sinθ1𝜃1/\sin\theta1 / roman_sin italic_θ serves as an additional factor that aids in the topological examination. The application of Duan’s ϕitalic-ϕ\phiitalic_ϕ-mapping theory is employed to introduce a new vector field denoted as ϕ=(ϕS,ϕθ)italic-ϕsuperscriptitalic-ϕ𝑆superscriptitalic-ϕ𝜃\phi=(\phi^{S},\phi^{\theta})italic_ϕ = ( italic_ϕ start_POSTSUPERSCRIPT italic_S end_POSTSUPERSCRIPT , italic_ϕ start_POSTSUPERSCRIPT italic_θ end_POSTSUPERSCRIPT ), where

ϕS=(SΨ)θ,xi,ϕθ=(θΨ)S,xi.formulae-sequencesuperscriptitalic-ϕ𝑆subscriptsubscript𝑆Ψ𝜃subscript𝑥𝑖superscriptitalic-ϕ𝜃subscriptsubscript𝜃Ψ𝑆subscript𝑥𝑖\phi^{S}=(\partial_{S}\Psi)_{\theta,x_{i}},\phi^{\theta}=(\partial_{\theta}% \Psi)_{S,x_{i}}.italic_ϕ start_POSTSUPERSCRIPT italic_S end_POSTSUPERSCRIPT = ( ∂ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT roman_Ψ ) start_POSTSUBSCRIPT italic_θ , italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT , italic_ϕ start_POSTSUPERSCRIPT italic_θ end_POSTSUPERSCRIPT = ( ∂ start_POSTSUBSCRIPT italic_θ end_POSTSUBSCRIPT roman_Ψ ) start_POSTSUBSCRIPT italic_S , italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT . (22)

Inspired by cq5 , in the case where θ𝜃\thetaitalic_θ equals π2𝜋2\frac{\pi}{2}divide start_ARG italic_π end_ARG start_ARG 2 end_ARG, the ϕitalic-ϕ\phiitalic_ϕ component of the vector field remains consistently at zero, allowing for a description of the topological flow as follows:

Jμ=12πϵμνλϵabνnaλnb,superscript𝐽𝜇12𝜋superscriptitalic-ϵ𝜇𝜈𝜆subscriptitalic-ϵ𝑎𝑏subscript𝜈superscript𝑛𝑎subscript𝜆superscript𝑛𝑏J^{\mu}=\frac{1}{2\pi}\epsilon^{\mu\nu\lambda}\epsilon_{ab}\partial_{\nu}n^{a}% \partial_{\lambda}n^{b},italic_J start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 italic_π end_ARG italic_ϵ start_POSTSUPERSCRIPT italic_μ italic_ν italic_λ end_POSTSUPERSCRIPT italic_ϵ start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT ∂ start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT , (23)

here ν=xνsubscript𝜈superscript𝑥𝜈\partial_{\nu}=\frac{\partial}{\partial x^{\nu}}∂ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT = divide start_ARG ∂ end_ARG start_ARG ∂ italic_x start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT end_ARG and xν=(t,r,θ)superscript𝑥𝜈𝑡𝑟𝜃x^{\nu}=(t,r,\theta)italic_x start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT = ( italic_t , italic_r , italic_θ ), the unit vector n𝑛nitalic_n is represented by n=(n1,n2)𝑛superscript𝑛1superscript𝑛2n=(n^{1},n^{2})italic_n = ( italic_n start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT , italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ), where n1=ϕSϕsuperscript𝑛1superscriptitalic-ϕ𝑆normitalic-ϕn^{1}=\frac{\phi^{S}}{\left\|\phi\right\|}italic_n start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT = divide start_ARG italic_ϕ start_POSTSUPERSCRIPT italic_S end_POSTSUPERSCRIPT end_ARG start_ARG ∥ italic_ϕ ∥ end_ARG, n2=ϕθϕsuperscript𝑛2superscriptitalic-ϕ𝜃normitalic-ϕn^{2}=\frac{\phi^{\theta}}{\left\|\phi\right\|}italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = divide start_ARG italic_ϕ start_POSTSUPERSCRIPT italic_θ end_POSTSUPERSCRIPT end_ARG start_ARG ∥ italic_ϕ ∥ end_ARG. The fulfillment of μJμ=0subscript𝜇superscript𝐽𝜇0\partial_{\mu}J^{\mu}=0∂ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_J start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT = 0 is a prerequisite for the topological current. In addition, the topological charge linked to a parameter region \sum can be expressed as follows:

Qt=j0d2x=i=1Nwi,subscript𝑄𝑡subscriptsuperscript𝑗0superscript𝑑2𝑥superscriptsubscript𝑖1𝑁subscript𝑤𝑖Q_{t}=\int_{\sum}j^{0}d^{2}x=\sum_{i=1}^{N}w_{i},italic_Q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = ∫ start_POSTSUBSCRIPT ∑ end_POSTSUBSCRIPT italic_j start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x = ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT italic_w start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , (24)

the winding number wisubscript𝑤𝑖w_{i}italic_w start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT determines the topological charge Q𝑄Qitalic_Q at critical points (ϕa(xi)=0)superscriptitalic-ϕ𝑎superscript𝑥𝑖0(\phi^{a}(x^{i})=0)( italic_ϕ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT ( italic_x start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ) = 0 ), with negative or positive values corresponding to a winding number of 11-1- 1 or +11+1+ 1 respectively. The summation of charges at all critical points within a thermodynamic system defines the total topological charge Qtsubscript𝑄𝑡Q_{t}italic_Q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT, this fundamental concept holds great significance in the analysis of thermodynamic systems as it aids in comprehending their global characteristics and facilitating classification. By aggregating the charges associated with each critical point, we can derive a quantitative indicator to depict the system’s topology structure. The classification of distinct thermodynamic systems is determined based on the total topological charge. The aforementioned approach will be employed to investigate dyonic AdS black holes with QTE in the context of EGB gravity.

Refer to caption
Refer to caption
Figure 1: d=4𝑑4d=4italic_d = 4: Up panel: The red arrows represent the vector field n𝑛nitalic_n for dyonic AdS black holes with QTE in the EGB background. The green dot are (r,θ)=(2.449490,π/2)𝑟𝜃2.449490𝜋2(r,\theta)=(2.449490,\pi/2)( italic_r , italic_θ ) = ( 2.449490 , italic_π / 2 ), which represent the critical points CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT. Down panel: The deflection angle Ω(ϑ)Ωitalic-ϑ\Omega(\vartheta)roman_Ω ( italic_ϑ ) as a function of ϑitalic-ϑ\varthetaitalic_ϑ for contours C1subscript𝐶1C_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (green curve), C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (brown curve).
Refer to caption
Refer to caption
Figure 2: d=5𝑑5d=5italic_d = 5: Up panel: The red arrows represent the vector field n𝑛nitalic_n for dyonic AdS black holes with QTE in the EGB background. The blue dot are (r,θ)=(2.658528,π/2)𝑟𝜃2.658528𝜋2(r,\theta)=(2.658528,\pi/2)( italic_r , italic_θ ) = ( 2.658528 , italic_π / 2 ), which represent the critical points CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT. Down panel: The deflection angle Ω(ϑ)Ωitalic-ϑ\Omega(\vartheta)roman_Ω ( italic_ϑ ) as a function of ϑitalic-ϑ\varthetaitalic_ϑ for contours C1subscript𝐶1C_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (blue curve), C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (brown curve).
Refer to caption
Refer to caption
Figure 3: d=6,α=0.1formulae-sequence𝑑6𝛼0.1d=6,\alpha=0.1italic_d = 6 , italic_α = 0.1: Up panel: The red arrows represent the vector field n𝑛nitalic_n for dyonic AdS black holes with QTE in the EGB background. The green dot are (r,θ)=(1.711584,π/2)𝑟𝜃1.711584𝜋2(r,\theta)=(1.711584,\pi/2)( italic_r , italic_θ ) = ( 1.711584 , italic_π / 2 ), which represent the critical points CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT. Down panel: The deflection angle Ω(ϑ)Ωitalic-ϑ\Omega(\vartheta)roman_Ω ( italic_ϑ ) as a function of ϑitalic-ϑ\varthetaitalic_ϑ for contours C1subscript𝐶1C_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (green curve), C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (brown curve).

III Topology of critical points

In this section, we will extensively examine the topological classification of critical points in black holes in different dimensions. By analyzing the metric of black holes in equation (5), we can deduce the thermodynamic function representing temperature without pressure

Refer to caption
Refer to caption
Figure 4: d=6(α=0.5,β=0.1,Qe=1,Qm=2)𝑑6formulae-sequence𝛼0.5formulae-sequence𝛽0.1formulae-sequencesubscript𝑄𝑒1subscript𝑄𝑚2d=6(\alpha=0.5,\beta=0.1,Q_{e}=1,Q_{m}=2)italic_d = 6 ( italic_α = 0.5 , italic_β = 0.1 , italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT = 1 , italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 2 ): Up panel: The red arrows represent the vector field n𝑛nitalic_n for dyonic AdS black holes with QTE in the EGB background. The colored dots are (r,θ)=(1.489048,π/2)𝑟𝜃1.489048𝜋2(r,\theta)=(1.489048,\pi/2)( italic_r , italic_θ ) = ( 1.489048 , italic_π / 2 ), (r,θ)=(2.706456,π/2)𝑟𝜃2.706456𝜋2(r,\theta)=(2.706456,\pi/2)( italic_r , italic_θ ) = ( 2.706456 , italic_π / 2 ) and (r,θ)=(1.873569,π/2)𝑟𝜃1.873569𝜋2(r,\theta)=(1.873569,\pi/2)( italic_r , italic_θ ) = ( 1.873569 , italic_π / 2 ) which represent the critical points CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT, respectively. Down panel:The deflection angle Ω(ϑ)Ωitalic-ϑ\Omega(\vartheta)roman_Ω ( italic_ϑ ) as a function of ϑitalic-ϑ\varthetaitalic_ϑ for contours C1subscript𝐶1C_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (green curve), C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (brown curve) and C3subscript𝐶3C_{3}italic_C start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT (blue curve) and C4subscript𝐶4C_{4}italic_C start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT (black curve)
Refer to caption
Refer to caption
Figure 5: d=6(α=1,β=0.1,Qe=1,Qm=2)𝑑6formulae-sequence𝛼1formulae-sequence𝛽0.1formulae-sequencesubscript𝑄𝑒1subscript𝑄𝑚2d=6(\alpha=1,\beta=0.1,Q_{e}=1,Q_{m}=2)italic_d = 6 ( italic_α = 1 , italic_β = 0.1 , italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT = 1 , italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 2 ): Up panel: The red arrows represent the vector field n𝑛nitalic_n for dyonic AdS black holes with QTE in the EGB background. The colored dots are (r,θ)=(1.008427,π/2)𝑟𝜃1.008427𝜋2(r,\theta)=(1.008427,\pi/2)( italic_r , italic_θ ) = ( 1.008427 , italic_π / 2 ), (r,θ)=(3.609964,π/2)𝑟𝜃3.609964𝜋2(r,\theta)=(3.609964,\pi/2)( italic_r , italic_θ ) = ( 3.609964 , italic_π / 2 ) and (r,θ)=(3.284875,π/2)𝑟𝜃3.284875𝜋2(r,\theta)=(3.284875,\pi/2)( italic_r , italic_θ ) = ( 3.284875 , italic_π / 2 ) which represent the critical points CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT, respectively. Down panel: The deflection angle Ω(ϑ)Ωitalic-ϑ\Omega(\vartheta)roman_Ω ( italic_ϑ ) as a function of ϑitalic-ϑ\varthetaitalic_ϑ for contours C1subscript𝐶1C_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (green curve), C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (brown curve) and C3subscript𝐶3C_{3}italic_C start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT (blue curve).
ΨΨ\displaystyle\Psiroman_Ψ =14πsinθr+(6(d4)(d3)α+r+2)absent14𝜋𝜃subscript𝑟6𝑑4𝑑3𝛼superscriptsubscript𝑟2\displaystyle=\frac{1}{4\pi\sin\theta r_{+}\left(6(d-4)(d-3)\alpha+r_{+}^{2}% \right)}= divide start_ARG 1 end_ARG start_ARG 4 italic_π roman_sin italic_θ italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ( 6 ( italic_d - 4 ) ( italic_d - 3 ) italic_α + italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG (25)
×[2(3d)(2(d5)(d4)α+r+2)\displaystyle\times\left[2(3-d)\left(2(d-5)(d-4)\alpha+r_{+}^{2}\right)\right.× [ 2 ( 3 - italic_d ) ( 2 ( italic_d - 5 ) ( italic_d - 4 ) italic_α + italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT )
+r+82d(Qm2+r+4dQe2(r+2d+8βΓ(d1)2r+4Qm2)2)].\displaystyle\left.+r_{+}^{8-2d}\left(Q_{m}^{2}+\frac{r_{+}^{4d}Q_{e}^{2}}{% \left(r_{+}^{2d}+8\beta\Gamma(d-1)^{2}r_{+}^{4}Q_{m}^{2}\right)^{2}}\right)% \right].+ italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 8 - 2 italic_d end_POSTSUPERSCRIPT ( italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 italic_d end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_d end_POSTSUPERSCRIPT + 8 italic_β roman_Γ ( italic_d - 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ] .

The vector field ϕitalic-ϕ\phiitalic_ϕ is a spatial vector function that can be separated into its radial component, denoted as ϕr+superscriptitalic-ϕsubscript𝑟\phi^{r_{+}}italic_ϕ start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT end_POSTSUPERSCRIPT, and its component related to the polar angle, denoted as ϕrθsuperscriptitalic-ϕsubscript𝑟𝜃\phi^{r_{\theta}}italic_ϕ start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_θ end_POSTSUBSCRIPT end_POSTSUPERSCRIPT. This decomposition aids in gaining a deeper understanding of the properties exhibited by vector fields. Through appropriate manipulation, we can derive the outcome represented by

ϕr+=cscθ4π(6(d4)(d3)α+r+2)2[4(d3)α\displaystyle\phi^{r_{+}}=\frac{\csc\theta}{4\pi\left(6(d-4)(d-3)\alpha+r_{+}^% {2}\right)^{2}}[4(d-3)\alphaitalic_ϕ start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT end_POSTSUPERSCRIPT = divide start_ARG roman_csc italic_θ end_ARG start_ARG 4 italic_π ( 6 ( italic_d - 4 ) ( italic_d - 3 ) italic_α + italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ 4 ( italic_d - 3 ) italic_α (26)
+r+2)+2r+82d(r+4dQe2(r+2d+8βΓ(d1)2r+4Qm2)2\displaystyle\left.+r_{+}^{2}\right)+2r_{+}^{8-2d}\left(\frac{r_{+}^{4d}Q_{e}^% {2}}{\left(r_{+}^{2d}+8\beta\Gamma(d-1)^{2}r_{+}^{4}Q_{m}^{2}\right)^{2}}\right.+ italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + 2 italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 8 - 2 italic_d end_POSTSUPERSCRIPT ( divide start_ARG italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 italic_d end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_d end_POSTSUPERSCRIPT + 8 italic_β roman_Γ ( italic_d - 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG
+Qm2)2(6(d4)(d3)α+r+2)\displaystyle\left.+Q_{m}^{2}\right)-2\left(6(d-4)(d-3)\alpha+r_{+}^{2}\right)+ italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) - 2 ( 6 ( italic_d - 4 ) ( italic_d - 3 ) italic_α + italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT )
×(62d(d4)r+6+4dQe2(r+2d+8βΓ(d1)2r+4Qm2)3\displaystyle\times\left(6-2d-\frac{(d-4)r_{+}^{6+4d}Q_{e}^{2}}{\left(r_{+}^{2% d}+8\beta\Gamma(d-1)^{2}r_{+}^{4}Q_{m}^{2}\right)^{3}}\right.× ( 6 - 2 italic_d - divide start_ARG ( italic_d - 4 ) italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 6 + 4 italic_d end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_d end_POSTSUPERSCRIPT + 8 italic_β roman_Γ ( italic_d - 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG
+r+62dQm2(+8dβΓ(d1)2r+4+4dQe2(r+2d+8βΓ(d1)2r+4Qm2)3\displaystyle+r_{+}^{6-2d}Q_{m}^{2}\left(+\frac{8d\beta\Gamma(d-1)^{2}r_{+}^{4% +4d}Q_{e}^{2}}{\left(r_{+}^{2d}+8\beta\Gamma(d-1)^{2}r_{+}^{4}Q_{m}^{2}\right)% ^{3}}\right.+ italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 6 - 2 italic_d end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( + divide start_ARG 8 italic_d italic_β roman_Γ ( italic_d - 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 + 4 italic_d end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_d end_POSTSUPERSCRIPT + 8 italic_β roman_Γ ( italic_d - 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG
+4d))+1r+2(6(d4)(d3)α+r+2)\displaystyle+4-d))+\frac{1}{r_{+}^{2}}\left(6(d-4)(d-3)\alpha+r_{+}^{2}\right)+ 4 - italic_d ) ) + divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ( 6 ( italic_d - 4 ) ( italic_d - 3 ) italic_α + italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT )
×(2(d3)(2(d5)(d4)α+r+2)\displaystyle\times\left(-2(d-3)\left(2(d-5)(d-4)\alpha+r_{+}^{2}\right)\right.× ( - 2 ( italic_d - 3 ) ( 2 ( italic_d - 5 ) ( italic_d - 4 ) italic_α + italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT )
+r+82d(Qm2+r+4dQe2(r+2d+8βΓ(d1)2r+4Qm2)2))],\displaystyle+\left.\left.r_{+}^{8-2d}\left(Q_{m}^{2}+\frac{r_{+}^{4d}Q_{e}^{2% }}{\left(r_{+}^{2d}+8\beta\Gamma(d-1)^{2}r_{+}^{4}Q_{m}^{2}\right)^{2}}\right)% \right)\right],+ italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 8 - 2 italic_d end_POSTSUPERSCRIPT ( italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 italic_d end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_d end_POSTSUPERSCRIPT + 8 italic_β roman_Γ ( italic_d - 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ) ] ,
ϕθsuperscriptitalic-ϕ𝜃\displaystyle\phi^{\theta}italic_ϕ start_POSTSUPERSCRIPT italic_θ end_POSTSUPERSCRIPT =cotθcscθ4πr+(6(d4)(d3)α+r+2)[2(3d)\displaystyle=\frac{\cot\theta\csc\theta}{4\pi r_{+}\left(6(d-4)(d-3)\alpha+r_% {+}^{2}\right)}\left[2(3-d)\right.= divide start_ARG roman_cot italic_θ roman_csc italic_θ end_ARG start_ARG 4 italic_π italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ( 6 ( italic_d - 4 ) ( italic_d - 3 ) italic_α + italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG [ 2 ( 3 - italic_d ) (27)
×2(d5)(d4)α+r+2)+r+82d\displaystyle\left.\times 2(d-5)(d-4)\alpha+r_{+}^{2}\right)+r_{+}^{8-2d}× 2 ( italic_d - 5 ) ( italic_d - 4 ) italic_α + italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 8 - 2 italic_d end_POSTSUPERSCRIPT
×(Qm2+r+4dQe2(r+2d+8βΓ(d1)2r+4Qm2)2))].\displaystyle\left.\left.\times\left(Q_{m}^{2}+\frac{r_{+}^{4d}Q_{e}^{2}}{% \left(r_{+}^{2d}+8\beta\Gamma(d-1)^{2}r_{+}^{4}Q_{m}^{2}\right)^{2}}\right)% \right)\right].× ( italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 italic_d end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_d end_POSTSUPERSCRIPT + 8 italic_β roman_Γ ( italic_d - 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ) ] .

The normalized vector field is defined by

n=(ϕrϕ,ϕθϕ).𝑛superscriptitalic-ϕ𝑟normitalic-ϕsuperscriptitalic-ϕ𝜃normitalic-ϕn=(\frac{\phi^{r}}{\left\|\phi\right\|},\frac{\phi^{\theta}}{\left\|\phi\right% \|}).italic_n = ( divide start_ARG italic_ϕ start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT end_ARG start_ARG ∥ italic_ϕ ∥ end_ARG , divide start_ARG italic_ϕ start_POSTSUPERSCRIPT italic_θ end_POSTSUPERSCRIPT end_ARG start_ARG ∥ italic_ϕ ∥ end_ARG ) . (28)

Next, we investigate the topological charge linked to the critical point across various dimensions. Through analysis and examination of the topological structure, the topological charge Q𝑄Qitalic_Q is assigned a specific value for each critical point, when encircling a single critical point CP𝐶𝑃CPitalic_C italic_P with a contour C𝐶Citalic_C, it is observed that its non-zero topological charge Q𝑄Qitalic_Q signifies significant information associated with the enclosed area. Conversely, if contour C𝐶Citalic_C lies outside of critical point CP𝐶𝑃CPitalic_C italic_P, its topological charge Q𝑄Qitalic_Q becomes zero cq6 ; cq7 . If there are many critical points within the contour, the total topological charge Q𝑄Qitalic_Q is obtained by summing up the individual topological charges associated with each critical point. To determine the topological charge at a critical point, in the orthogonal plane, we employ the method of Duan’s φ𝜑\varphiitalic_φ-mapping theory to ensure that the constructed contour exhibits a positive value. We can transform the critical point on the (S,θ)𝑆𝜃(S,\theta)( italic_S , italic_θ ) plane to the (r,θ)𝑟𝜃(r,\theta)( italic_r , italic_θ ) plane using reparameterization. In this case, the critical point is situated at (r0,π2)subscript𝑟0𝜋2(r_{0},\frac{\pi}{2})( italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , divide start_ARG italic_π end_ARG start_ARG 2 end_ARG )

{r=acosϑ+r0θ=bsinϑ+π2.\left\{\begin{matrix}r=a\cos\vartheta+r_{0}\\ \theta=b\sin\vartheta+\frac{\pi}{2}.\end{matrix}\right.{ start_ARG start_ROW start_CELL italic_r = italic_a roman_cos italic_ϑ + italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_θ = italic_b roman_sin italic_ϑ + divide start_ARG italic_π end_ARG start_ARG 2 end_ARG . end_CELL end_ROW end_ARG (29)

where ϑ(0,2π)italic-ϑ02𝜋\vartheta\subset(0,2\pi)italic_ϑ ⊂ ( 0 , 2 italic_π ). After this transformation, we can express the deflection angle Ω(ϑ)Ωitalic-ϑ\Omega(\vartheta)roman_Ω ( italic_ϑ ) is given by

Ω(ϑ)=0ϑϵabnaϑnbdϑ.Ωitalic-ϑsuperscriptsubscript0italic-ϑsubscriptitalic-ϵ𝑎𝑏superscript𝑛𝑎subscriptitalic-ϑsuperscript𝑛𝑏𝑑italic-ϑ\Omega(\vartheta)=\int_{0}^{\vartheta}\epsilon_{ab}n^{a}\partial_{\vartheta}n^% {b}d\vartheta.roman_Ω ( italic_ϑ ) = ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ϑ end_POSTSUPERSCRIPT italic_ϵ start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT ∂ start_POSTSUBSCRIPT italic_ϑ end_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT italic_d italic_ϑ . (30)

Naturally, the topological charge can be determined by the deflection angle as follows:

Q=Ω(2π)2π.𝑄Ω2𝜋2𝜋Q=\frac{\Omega(2\pi)}{2\pi}.italic_Q = divide start_ARG roman_Ω ( 2 italic_π ) end_ARG start_ARG 2 italic_π end_ARG . (31)

Now, let us investigate the classification of critical points pertaining to AdS black holes. Research has indicated that black holes demonstrate characteristic swallowtail behavior across various dimensions of space-time and undergo small/large BH phase transitions resembling those observed in van der Waals liquid/gas systems cq3 . Interestingly, in the case of d=6𝑑6d=6italic_d = 6, when the coupling constant α𝛼\alphaitalic_α falls within a particular parameter range (α0.5)𝛼0.5(\alpha\geq 0.5)( italic_α ≥ 0.5 ), a noteworthy phenomenon occurs where the black hole experiences a phase transition characterized by small/intermediate/large BH phases. Additionally, this leads to the emergence of a triple point naturally, signifying that the black hole possesses three critical points. In the case where the coupling constant α=1𝛼1\alpha=1italic_α = 1, it can be observed that despite the presence of three critical points, there are indications of two first-order phase transitions within the system as suggested by the existence of two swallowtails.

Table 1: Critical values at Qe=1,Qm=2,β=0.1formulae-sequencesubscript𝑄𝑒1formulae-sequencesubscript𝑄𝑚2𝛽0.1Q_{e}=1,Q_{m}=2,\beta=0.1italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT = 1 , italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 2 , italic_β = 0.1 for dyonic AdS black holes with QTE in the EGB background.
CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT
rcsubscript𝑟𝑐r_{c}italic_r start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 2.449490 - -
(α=0.5,d=4formulae-sequence𝛼0.5𝑑4\alpha=0.5,d=4italic_α = 0.5 , italic_d = 4) Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 0.043317 - -
Pcsubscript𝑃𝑐P_{c}italic_P start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 0.003316 - -
rcsubscript𝑟𝑐r_{c}italic_r start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 2.658528 - -
(α=0.5,d=5formulae-sequence𝛼0.5𝑑5\alpha=0.5,d=5italic_α = 0.5 , italic_d = 5) Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 0.063461 - -
Pcsubscript𝑃𝑐P_{c}italic_P start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 0.006193 - -
rcsubscript𝑟𝑐r_{c}italic_r start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 1.711584 - -
(α=0.1,d=6formulae-sequence𝛼0.1𝑑6\alpha=0.1,d=6italic_α = 0.1 , italic_d = 6) Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 0.138928 - -
Pcsubscript𝑃𝑐P_{c}italic_P start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 0.027903 - -
rcsubscript𝑟𝑐r_{c}italic_r start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 1.489048 2.706456 1.873569
(α=0.5,d=6formulae-sequence𝛼0.5𝑑6\alpha=0.5,d=6italic_α = 0.5 , italic_d = 6) Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 0.064736 0.064872 0.064642
Pcsubscript𝑃𝑐P_{c}italic_P start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 0.006536 0.006576 0.006353
rcsubscript𝑟𝑐r_{c}italic_r start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 1.008427 3.609964 3.284875
(α=1,d=6formulae-sequence𝛼1𝑑6\alpha=1,d=6italic_α = 1 , italic_d = 6) Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 0.055925 0.045933 0.045928
Pcsubscript𝑃𝑐P_{c}italic_P start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 0.050611 0.003311 0.003308
rcsubscript𝑟𝑐r_{c}italic_r start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 0.929668 - -
(α=0.5,d=7formulae-sequence𝛼0.5𝑑7\alpha=0.5,d=7italic_α = 0.5 , italic_d = 7) Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 0.113123 - -
Pcsubscript𝑃𝑐P_{c}italic_P start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 0.288091 - -

However, the intermediate BH branch is unable to participate in the phase transition.

We aim to topologically classify the critical points for the dyonic AdS black holes with QTE in the EGB background. Based on the data presented in TABLE I, it is clear that when d𝑑ditalic_d takes values of 4, 5, and 7, there exists only one critical point for the black hole. The system undergoes a small/large phase transition of the black hole. For the case of d=4𝑑4d=4italic_d = 4, the characteristics of θ𝜃\thetaitalic_θ and radius r+subscript𝑟r_{+}italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT are presented in relation to the vector field n𝑛nitalic_n in FIG. 1. Identify a critical point located at (θ,r+)𝜃subscript𝑟(\theta,r_{+})( italic_θ , italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT )=(π2,2.449490)𝜋22.449490(\frac{\pi}{2},2.449490)( divide start_ARG italic_π end_ARG start_ARG 2 end_ARG , 2.449490 ). Create two contours, namely C1subscript𝐶1C_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, where C1subscript𝐶1C_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT should enclose critical point CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, while no critical points should be enclosed by C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (refer to TABLE II for contour parameters). In addition, the values of Ω(ϑ)Ωitalic-ϑ\Omega(\vartheta)roman_Ω ( italic_ϑ ) for contours C1subscript𝐶1C_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are observed to be 2π2𝜋-2\pi- 2 italic_π and 00, respectively. The topological charge of critical point CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT was found to be negative. In contrast, contour C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT did not enclose any critical point, resulting in a topological charge of zero. Inspired by cq4 , CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is classified as a conventional critical point possessing a topological charge of Q=QCP1=1𝑄subscript𝑄𝐶subscript𝑃11Q=Q_{CP_{1}}=-1italic_Q = italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - 1. The presence of only one critical point in the system, which is called as a physical critical point. Similarly, in the case of space-time dimensions d=5𝑑5d=5italic_d = 5 and 7777, we observe analogous outcomes as for d=4𝑑4d=4italic_d = 4. The distinction lies in the fact that the coordinates of the physical critical points differ (d=5,(θ,r+)=(π2,2.658528),d=7,(θ,r+)=(π2,0.929668))formulae-sequence𝑑5formulae-sequence𝜃subscript𝑟𝜋22.658528formulae-sequence𝑑7𝜃subscript𝑟𝜋20.929668(d=5,(\theta,r_{+})=(\frac{\pi}{2},2.658528),d=7,(\theta,r_{+})=(\frac{\pi}{2}% ,0.929668))( italic_d = 5 , ( italic_θ , italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) = ( divide start_ARG italic_π end_ARG start_ARG 2 end_ARG , 2.658528 ) , italic_d = 7 , ( italic_θ , italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) = ( divide start_ARG italic_π end_ARG start_ARG 2 end_ARG , 0.929668 ) ). Further details can be found in FIG.2 and FIG.6.

Table 2: Contours at Qe=1,Qm=2,β=0.1formulae-sequencesubscript𝑄𝑒1formulae-sequencesubscript𝑄𝑚2𝛽0.1Q_{e}=1,Q_{m}=2,\beta=0.1italic_Q start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT = 1 , italic_Q start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 2 , italic_β = 0.1 for dyonic AdS black holes with QTE in the EGB background.
C1subscript𝐶1C_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT C3subscript𝐶3C_{3}italic_C start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT C4subscript𝐶4C_{4}italic_C start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT
a𝑎aitalic_a 0.07 0.07 - -
(α=0.5,d=4formulae-sequence𝛼0.5𝑑4\alpha=0.5,d=4italic_α = 0.5 , italic_d = 4) b𝑏bitalic_b 0.4 0.4 - -
r0subscript𝑟0r_{0}italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT 2.449490 3.4 - -
a𝑎aitalic_a 0.07 0.07 - -
(α=0.5,d=5formulae-sequence𝛼0.5𝑑5\alpha=0.5,d=5italic_α = 0.5 , italic_d = 5) b𝑏bitalic_b 0.4 0.4 - -
r0subscript𝑟0r_{0}italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT 2.658528 3.8 - -
a𝑎aitalic_a 0.15 0.15 - -
(α=0.1,d=6formulae-sequence𝛼0.1𝑑6\alpha=0.1,d=6italic_α = 0.1 , italic_d = 6) b𝑏bitalic_b 0.4 0.4 - -
r0subscript𝑟0r_{0}italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT 2.57753 2.4 - -
a𝑎aitalic_a 0.15 0.15 0.15 1.1
(α=0.5,d=6formulae-sequence𝛼0.5𝑑6\alpha=0.5,d=6italic_α = 0.5 , italic_d = 6) b𝑏bitalic_b 0.4 0.4 0.4 0.5
r0subscript𝑟0r_{0}italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT 1.489048 2.706456 1.873569 2.09775
a𝑎aitalic_a 0.15 0.15 0.15 -
6(α=1,d=6formulae-sequence𝛼1𝑑6\alpha=1,d=6italic_α = 1 , italic_d = 6) b𝑏bitalic_b 0.4 0.4 0.4 -
r0subscript𝑟0r_{0}italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT 1.008427 3.609964 3.284875 -
a𝑎aitalic_a 0.15 0.15 - -
6(α=0.5,d=7formulae-sequence𝛼0.5𝑑7\alpha=0.5,d=7italic_α = 0.5 , italic_d = 7) b𝑏bitalic_b 0.4 0.4 - -
r0subscript𝑟0r_{0}italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT 0.929668 2 - -

Next, let’s analyze the situation where d=6.𝑑6d=6.italic_d = 6 .
(i)𝑖(i)( italic_i ): α=0.1,β=0.1formulae-sequence𝛼0.1𝛽0.1\alpha=0.1,\beta=0.1italic_α = 0.1 , italic_β = 0.1
The presence of a conventional critical point for a black hole can be observed at coordinates (1.711584,π/2)1.711584𝜋2(1.711584,\pi/2)( 1.711584 , italic_π / 2 ), with a topological charge of (Q=QCP1=1)𝑄subscript𝑄𝐶subscript𝑃11(Q=Q_{CP_{1}}=-1)( italic_Q = italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - 1 ). This result exhibits remarkable similarities to findings in other dimensions (d=4,5,7)𝑑457(d=4,5,7)( italic_d = 4 , 5 , 7 ).
(ii)𝑖𝑖(ii)( italic_i italic_i ): α=0.5,β=0.1formulae-sequence𝛼0.5𝛽0.1\alpha=0.5,\beta=0.1italic_α = 0.5 , italic_β = 0.1
According to the findings cq3 , it is explicitly stated that the LBH/ IBH/LBH phase transition occurs occurs at α=0.5𝛼0.5\alpha=0.5italic_α = 0.5 (refer to TABLE.I for other parameter values). The contours C1subscript𝐶1C_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, and C3subscript𝐶3C_{3}italic_C start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT surround the critical points CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (1.489048,π/2)1.489048𝜋2(1.489048,\pi/2)( 1.489048 , italic_π / 2 ), CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (2.706456,π/2)2.706456𝜋2(2.706456,\pi/2)( 2.706456 , italic_π / 2 ), CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT (1.873569,π/2)1.873569𝜋2(1.873569,\pi/2)( 1.873569 , italic_π / 2 ), respectively. And contour C4subscript𝐶4C_{4}italic_C start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT encircle three critical points. We can find from the deflection angle diagram that QCP1=QCP2=1subscript𝑄𝐶subscript𝑃1subscript𝑄𝐶subscript𝑃21Q_{CP_{1}}=Q_{CP_{2}}=-1italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - 1, QCP3=1subscript𝑄𝐶subscript𝑃31Q_{CP_{3}}=1italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 1. Thus, this means that CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are the conventional critical points, while CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT is the novel critical point. In addition, we see that C4subscript𝐶4C_{4}italic_C start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT is equal to -2π𝜋\piitalic_π, meaning that the topological charge of the three critical points surrounded by C4subscript𝐶4C_{4}italic_C start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT is 11-1- 1. In addition, the total topological charge can also be calculated as Q=QCP1+QCP2+QCP3=1𝑄subscript𝑄𝐶subscript𝑃1subscript𝑄𝐶subscript𝑃2subscript𝑄𝐶subscript𝑃31Q=Q_{CP_{1}}+Q_{CP_{2}}+Q_{CP_{3}}=-1italic_Q = italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT + italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT + italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - 1. The triple point phenomenon is intriguing due to its manifestation of the distinctive characteristics associated with black hole phase transition. For this occurrence to take place, it is essential for the CP1𝐶subscript𝑃1{CP_{1}}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and CP2𝐶subscript𝑃2{CP_{2}}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT to be physical. This implies that under specific temperature and pressure ranges, a phase transition occurs in black holes, resulting in the coexistence of three distinct phases simultaneously. However, noted that the novel critical point CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT is considered non-physical as it cannot be attained by manipulating external conditions such as temperature and pressure. Unlike the conventional critical point, achieving minimal Gibbs free energy at the novel critical point CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT is unattainable.
(iii)𝑖𝑖𝑖(iii)( italic_i italic_i italic_i ): α=1,β=0.1formulae-sequence𝛼1𝛽0.1\alpha=1,\beta=0.1italic_α = 1 , italic_β = 0.1
We find that there are two swallowtails curves in the G-T diagram cq3 in the range of pressure CP3<P<CP2𝐶subscript𝑃3𝑃𝐶subscript𝑃2CP_{3}<P<CP_{2}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT < italic_P < italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, we will investigate from a topological point of view why the system has three critical points and only small/large BH phase transition occurs. The first observation is that the contours C1subscript𝐶1C_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, and C3subscript𝐶3C_{3}italic_C start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT enclose the critical points CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (1.008427,π/2)1.008427𝜋2(1.008427,\pi/2)( 1.008427 , italic_π / 2 ), CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (3.609964,π/2)3.609964𝜋2(3.609964,\pi/2)( 3.609964 , italic_π / 2 ), CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT (3.284875,π/2)3.284875𝜋2(3.284875,\pi/2)( 3.284875 , italic_π / 2 ) respectively within the vector field in FIG.5. And the contours C1subscript𝐶1C_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, and C3subscript𝐶3C_{3}italic_C start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT corresponding deflection angle Ω(ϑ)=2π,2πΩitalic-ϑ2𝜋2𝜋\Omega(\vartheta)=-2\pi,-2\piroman_Ω ( italic_ϑ ) = - 2 italic_π , - 2 italic_π and 2π2𝜋2\pi2 italic_π for ϑ=2πitalic-ϑ2𝜋\vartheta=2\piitalic_ϑ = 2 italic_π, respectively, this means that the black holes associated topological charge of QCP1=subscript𝑄𝐶subscript𝑃1Q_{CP_{1}}=-italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = -1, QCP2=1subscript𝑄𝐶subscript𝑃21Q_{CP_{2}}=1italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 1, and QCP3=1subscript𝑄𝐶subscript𝑃31Q_{CP_{3}}=1italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 1. As a result, the total topological charge can be determined by Q=QCP1+QCP2+QCP3=1𝑄subscript𝑄𝐶subscript𝑃1subscript𝑄𝐶subscript𝑃2subscript𝑄𝐶subscript𝑃31Q=Q_{CP_{1}}+Q_{CP_{2}}+Q_{CP_{3}}=-1italic_Q = italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT + italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT + italic_Q start_POSTSUBSCRIPT italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - 1. Naturally, the conventional critical point is denoted as CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, whereas the novel critical point is represented by CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT. As the system undergoes first-order transition, it can be inferred that only CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT represents the conventional (physical) critical point. Furthermore, CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT are associated with non-physical critical points. In addition, the findings indicate that the conventional critical point is considered as an indicator of a first-order phase transition , whereas the novel critical point does not possess the ability to minimize Gibbs free energy and thus is defined as the critical point associated with non-physical characteristics ni1 . However, the conventional critical point CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT lacks the ability to minimize the Gibbs free energy and is defined as a non-physical critical point. This suggests that the aforementioned conclusion cannot be universally applied to all systems.

Refer to caption
Refer to caption
Figure 6: d=7𝑑7d=7italic_d = 7: Up panel: The red arrows represent the vector field n𝑛nitalic_n for dyonic AdS black holes with QTE in the EGB background. The blue dot are (r,θ)=(0.929668,π/2)𝑟𝜃0.929668𝜋2(r,\theta)=(0.929668,\pi/2)( italic_r , italic_θ ) = ( 0.929668 , italic_π / 2 ), which represent the critical points CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT. Down panel: The deflection angle Ω(ϑ)Ωitalic-ϑ\Omega(\vartheta)roman_Ω ( italic_ϑ ) as a function of ϑitalic-ϑ\varthetaitalic_ϑ for contours C1subscript𝐶1C_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (blue curve), C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (green curve).
Refer to caption
Refer to caption
Refer to caption
Figure 7: a isobaric curves (colored solid curves) for dyonic AdS black holes with QTE in the EGB background. shown in the Tr+𝑇subscript𝑟T-r_{+}italic_T - italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT plane. The black curve is for the extremal points of the temperature. Up panel:(α=0.1)𝛼0.1(\alpha=0.1)( italic_α = 0.1 ),Middle panel:(α=0.5)𝛼0.5(\alpha=0.5)( italic_α = 0.5 ), Down panel:(α=1)𝛼1(\alpha=1)( italic_α = 1 )

IV Nature of Critical points

In the scenario where (d=6,α=0.1)formulae-sequence𝑑6𝛼0.1(d=6,\alpha=0.1)( italic_d = 6 , italic_α = 0.1 ) is considered, it can be observed that the black hole experiences small/large phase transition with only one critical point. Similar conclusions can be drawn for d=4𝑑4d=4italic_d = 4, 5555, and 7777. Hence, this section will primarily focus on examining the critical point of the black hole at d=6𝑑6d=6italic_d = 6. To this end, the equation of state ( 12 ) is employed to generate the isobaric curves(represented by colored curve) and the critical temperature (depicted as a black curve) in FIG.7. We can find some interesting results as follows:
(i)𝑖(i)( italic_i ) α=0.1,β=0.1formulae-sequence𝛼0.1𝛽0.1\alpha=0.1,\beta=0.1italic_α = 0.1 , italic_β = 0.1
Up panel: Each isopiestic curve corresponds to distinct maximum extreme points of temperature, with the red curve representing P=0.02𝑃0.02P=0.02italic_P = 0.02, the green curve representing P=0.025𝑃0.025P=0.025italic_P = 0.025, the blue curve corresponding to PCP1=0.027903subscript𝑃𝐶𝑃10.027903P_{CP1}=0.027903italic_P start_POSTSUBSCRIPT italic_C italic_P 1 end_POSTSUBSCRIPT = 0.027903, and the orange curve corresponding to P=0.03𝑃0.03P=0.03italic_P = 0.03. Additionally, it is noteworthy that the conventional critical point (QCP1=1)subscript𝑄𝐶𝑃11(Q_{CP1}=-1)( italic_Q start_POSTSUBSCRIPT italic_C italic_P 1 end_POSTSUBSCRIPT = - 1 ) coincides with the maximum temperature extreme point. By employing Maxwell’s equal area law, we can eliminate the unstable region near CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (referred to as the intermediate black hole branch). Moreover, our findings reveal that the isobaric curves are divided into distinct regions with varying phase numbers by the critical temperature, and the arrangement of phases in proximity to the critical point is determined based on these regions. Interestingly, an increase in pressure (redorange)𝑟𝑒𝑑𝑜𝑟𝑎𝑛𝑔𝑒(red\rightarrow orange)( italic_r italic_e italic_d → italic_o italic_r italic_a italic_n italic_g italic_e ) results in a reduction of isopiestic lines near CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT. Consequently, CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT can be regarded as a phase annihilation point.
(ii)𝑖𝑖(ii)( italic_i italic_i ) α=0.5,β=0.1formulae-sequence𝛼0.5𝛽0.1\alpha=0.5,\beta=0.1italic_α = 0.5 , italic_β = 0.1
Middle panel: We can observe that the isobaric curves PCP3=0.06353(purple)subscript𝑃𝐶𝑃30.06353𝑝𝑢𝑟𝑝𝑙𝑒P_{CP3}=0.06353(purple)italic_P start_POSTSUBSCRIPT italic_C italic_P 3 end_POSTSUBSCRIPT = 0.06353 ( italic_p italic_u italic_r italic_p italic_l italic_e );P=0.006469(orange)𝑃0.006469𝑜𝑟𝑎𝑛𝑔𝑒P=0.006469(orange)italic_P = 0.006469 ( italic_o italic_r italic_a italic_n italic_g italic_e ); PCP1=0.006536(blue)subscript𝑃𝐶𝑃10.006536𝑏𝑙𝑢𝑒P_{CP1}=0.006536(blue)italic_P start_POSTSUBSCRIPT italic_C italic_P 1 end_POSTSUBSCRIPT = 0.006536 ( italic_b italic_l italic_u italic_e );PCP2=0.0065757(red)subscript𝑃𝐶𝑃20.0065757𝑟𝑒𝑑P_{CP2}=0.0065757(red)italic_P start_POSTSUBSCRIPT italic_C italic_P 2 end_POSTSUBSCRIPT = 0.0065757 ( italic_r italic_e italic_d ) exhibits temperature extreme points. We have observed that the maximum extreme temperature is associated with the green dots, which represent conventional critical points CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT. On the other hand, the blue point (novel critical point CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT) accurately represents the minimum extreme temperature. By applying Maxwell’s equal area law, we can effectively eliminate unstable regions near CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT; however, this approach may not be successful for unstable regions near CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT. These findings align perfectly with the results cq4 . Furthermore, as pressure increases (indicated by orange-red), we notice a decrease in the number of regions surrounding CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, as well as a decrease in the number of regions near CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT. We designate CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT as the phase generation point while referring to critical points CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT as phase annihilation points.
(iii)𝑖𝑖𝑖(iii)( italic_i italic_i italic_i ) α=1,β=0.1formulae-sequence𝛼1𝛽0.1\alpha=1,\beta=0.1italic_α = 1 , italic_β = 0.1
Down panel: We can find than the isobaric curves (PCP3(P_{CP3}( italic_P start_POSTSUBSCRIPT italic_C italic_P 3 end_POSTSUBSCRIPT=0.03308(purple)0.03308𝑝𝑢𝑟𝑝𝑙𝑒0.03308(purple)0.03308 ( italic_p italic_u italic_r italic_p italic_l italic_e ); P𝑃Pitalic_P=0.003309(orange)0.003309𝑜𝑟𝑎𝑛𝑔𝑒0.003309(orange)0.003309 ( italic_o italic_r italic_a italic_n italic_g italic_e ); PCP2subscript𝑃𝐶𝑃2P_{CP2}italic_P start_POSTSUBSCRIPT italic_C italic_P 2 end_POSTSUBSCRIPT=0.003311(blue)0.003311𝑏𝑙𝑢𝑒0.003311(blue)0.003311 ( italic_b italic_l italic_u italic_e ); P=0.025(green);(PCP1=0.050611(red))𝑃0.025𝑔𝑟𝑒𝑒𝑛subscript𝑃𝐶𝑃10.050611𝑟𝑒𝑑P=0.025(green);(P_{CP1}=0.050611(red))italic_P = 0.025 ( italic_g italic_r italic_e italic_e italic_n ) ; ( italic_P start_POSTSUBSCRIPT italic_C italic_P 1 end_POSTSUBSCRIPT = 0.050611 ( italic_r italic_e italic_d ) ) exhibit temperature extreme points. Similarly, the green dots (representing conventional critical points CP1𝐶subscript𝑃1CP_{1}italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT) correspond to the point of maximum extreme temperature point, while the blue dot (novel critical point CP3𝐶subscript𝑃3CP_{3}italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT) accurately corresponds to the point of minimum extreme temperature point. Based on research findingscq3 , it has been observed that there is no occurrence of a first-order phase transition in close proximity to the conventional critical point CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, which can also be found in we2 ; we7 . Notably, the application of Maxwell’s law of equal area can effectively remove the unstable region near CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT. However, it is not suitable for indicating a first-order phase transition. Therefore, an alternative classification should be employed, it can be inferred that CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT may serve as a phase annihilation point.

V Conclusions

In summary, we employ the Duan ’s φ𝜑\varphiitalic_φ-mapping theory to examine the topological structure of thermodynamics in black holes. The categorization of critical points for black holes is carried out across diverse dimensions of space-time. In the scenario where d=6𝑑6d=6italic_d = 6, the black hole will have triple point when α𝛼\alphaitalic_α specific value is assumed by the GB coupling constant. Our focus lies on assessing how this coupling constant impacts the critical point. This investigation contributes to a more profound comprehension of EGB modified black holes. The primary highlights can be summarized as follows:
(i)𝑖(i)( italic_i ) Black holes have a critical point in any space-time dimension, there is a small/large BH phase transition, which is a conventional (physical) critical point due to the total topological charge 11-1- 1. Different d𝑑ditalic_d will lead to different positions of the critical point in vector field n𝑛nitalic_n. Initially, as d increases, the critical point also increases before eventually decreasing again in TABLE I.
(ii)𝑖𝑖(ii)( italic_i italic_i ) For d=6𝑑6d=6italic_d = 6, we will find that when the coupling constant satisfies α=0.5𝛼0.5\alpha=0.5italic_α = 0.5 and 1111, the system will appear triple points, but the phase structure is completely different. More specifically, the three critical points can be divided into two conventional ( physical ) critical points (CP1,CP2)𝐶subscript𝑃1𝐶subscript𝑃2(CP_{1},CP_{2})( italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) and a novel critical point (CP3)𝐶subscript𝑃3(CP_{3})( italic_C italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ), the novel critical point is not a physical point because it does not minimize the Gibbs free energy. When a = 1, the system has only one traditional ( physical ) critical point (CP1)𝐶subscript𝑃1(CP_{1})( italic_C italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ). In both cases, the total topological charge is 11-1- 1.
(iii)𝑖𝑖𝑖(iii)( italic_i italic_i italic_i ) From the isobaric curve and critical temperature in FIG.7, we can observe that the conventional critical point is located at the maximum extreme point of temperature, and the unstable region near it can be eliminated by Maxwell ’s equal area law. For the case of coupling constant α=1𝛼1\alpha=1italic_α = 1 , it is found that the critical point CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is no longer applicable. It should be noted that the novel critical point appears at the minimum extreme point of temperature. The results also reveal a decrease in the number of phases along the isobaric curve at the conventional critical point as pressure increases (α=1)𝛼1(\alpha=1)( italic_α = 1 ), while the number of phases increases at the novel critical point. Therefore, we define the CP2𝐶subscript𝑃2CP_{2}italic_C italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT as the phase annihilation point (α=1)𝛼1(\alpha=1)( italic_α = 1 ).

Acknowledgements.
This work is supported by the Doctoral Foundation of Zunyi Normal University of China (BS [2022] 07, QJJ-[2022]-314), and the National Natural Science Foundation of China (Grant Nos. 12265007 and 12264061)

References

  • (1) B. P. Abbott et al. (LIGO Scientific, Virgo), Phys. Rev. Lett., 116 (2016) 241103.
  • (2) B. P. Abbott et al. (LIGO Scientific, Virgo), Phys. Rev. Lett., 116 (2016) 061102.
  • (3) B. P. Abbott et al. (LIGO Scientific, Virgo), Phys. Rev. Lett., 119 (2017) 141101.
  • (4) B. P. Abbott et al. (LIGO Scientific, Virgo), Phys. Rev. Lett., 119 (2017) 161101.
  • (5) S.W. Hawking. Commun, Math. Phys., 43 (1975) 199.
  • (6) S.W. Hawking and D.N. Page, Commun. Math. Phys., 87 (1983) 577.
  • (7) E.Witten, Adv. Theor. Math. Phys., 2 (1998) 505.
  • (8) D. Kastor, S. Ray and J. Traschen, Class. Quant. Grav., 26 (2009)195011.
  • (9) D. Kastor, S. Ray and J. Traschen, Class. Quant. Grav., 27 (2010) 235014.
  • (10) D. Kastor, S. Ray and J. Traschen, Class. Quant. Grav., 28 (2011) 195022.
  • (11) D. Kastor, S. Ray and J. Traschen, Class. Quant. Grav., 36 (2018) 024002.
  • (12) D. Kubiznak and R. B.Mann, J. High Energy Phys., 1207 (2012) 033.
  • (13) S. H. Hendi and M. H. Vahidinia, Phys. Rev. D 88 (2012) 084045.
  • (14) S. W. Wei and Y. X. Liu, Phys. Rev. D 87 (2013) 044014.
  • (15) S. W. Wei, P. Cheng and Y. X. Liu, Phy. Rev. D 93 (2016) 084015.
  • (16) D. C. Zou, S. J. Zhang and B. Wang, Phy. Rev. D 89 (2014) 044002.
  • (17) D. C. Zou, Y. Liu and B. Wang, Phy. Rev. D 90 (2014) 044063.
  • (18) S. Gunasekaran, R. B. Mann and D. Kubiznak, J. High Energy Phys., 11 (2012) 110.
  • (19) N. Altamirano, D. Kubiznak and R. B. Mann, Phys. Rev. D 88 (2013) 101502.
  • (20) C.M. Zhang, D. C. Zou and M. Zhang, Phys. Lett. B 811 (2020)135955.
  • (21) M. D. Li, H. M. Wang and S. W. Wei, Phys. Rev. D 105 (2022) 104048.
  • (22) R. A. Hennigar and R. B. Mann, Phys. Rev. Lett., 118 (2017) 021301.
  • (23) S.W. Wei, Y .X. Liu, Phys. Rev. D 105 (2022) 104003.
  • (24) S.W. Wei, Y.X. Liu and R. B. Mann, Phys. Rev. Lett., 129 (2022) 191101.
  • (25) P.Ku. Yerra and C. Bhamidipati, Phys.Rev.D 105 (2022) 104053.
  • (26) N.-C. Bai, L. Li and J. Tao, Phys. Rev. D 107 (2023) 064015.
  • (27) P.Ku. Yerra and C. Bhamidipati, Phys. Lett B 835 (2022) 137591.
  • (28) P.Ku. Yerra and C. Bhamidipati, Phys. Rev. D 106 (2022) 064059.
  • (29) C. Fang, J. Jiang and M. Zhang J. High Energy Phys., 01 (2023) 102.
  • (30) N.J. Gogoiand P. Phukon, Phys. Rev. D 107 (2023) 106009.
  • (31) M.R. Alipour, Mohammad Ali S. Afshar, S. N.Gashti and J. Sadeghi, Phys.Dark Univ. 42 (2023) 101361.
  • (32) D.Wu, Eur.Phys.J.C 83 (2023) 365.
  • (33) D.Wu, Eur.Phys.J.C 83 (2023) 589.
  • (34) D.Wu, Phys.Rev.D 107 (2023) 084002.
  • (35) D.Wu, Phys.Rev.D 107 (2023) 024024.
  • (36) D.Wu, Phys.Rev.D 108 (2023) 084041.
  • (37) D. Chen, Y. He, J. Tao and W. Yang Eur. Phys. J. C 84 (2024) 96.
  • (38) D. Chen, Y. He and J. Tao Eur. Phys. J. C 83 (2023) 872.
  • (39) M.-Y. Zhang, H. Chen, H. Hassanabadi,Z.-W. Long and H. Yang, Eur.Phys. J. C 83 (2023) 773.
  • (40) R. G. Cai, L. M. Cao, L. Li and R. Q. Yang, J. High Energy Phys. 09 (2013) 005.
  • (41) J. X. Mo and W. B. Liu, Phys. Rev. D 89 (2014) 084057.
  • (42) S. G. Ghosh, M. Amir and S. D. Maharaj, Eur. Phys. J. C 77 (2017) 530.
  • (43) D. Glavan and C. Lin, Phys. Rev. Lett. 124 (2020) 081301.
  • (44) K. Yang, B. M. Gu, S. W. Wei and Y. X. Liu, Eur. Phys. J. C 80 (2020) 662.
  • (45) Z. C. Lin, K. Yang, S. W. Wei, Y. Q. Wang and Y. X. Liu, Eur. Phys. J. C 80 (2020) 1033.
  • (46) D. V. Singh, B. K. Singh and S. Upadhyay, Annals Phys. 434 (2021) 168642.
  • (47) H. S. Liu, Z. F. Mai, Y. Z. Li and H. Lu, Sci. China Phys. Mech. Astron. 63 (2020) 240411.
  • (48) Y. Sekhmani, H. Lekbich, A. El Boukili and M. B. Sedra, Eur. Phys. J. C 82 (2022) 1087.
  • (49) P.-H. Mou, Q.-Q. Jiang, K.-J. He and G.-P. Li, Triple points and phase transitions of D-dimensional dyonic AdS black holes with quasitopological electromagnetism in Einstein-Gauss-Bonnet gravity, arXiv: 2310.08010
  • (50) R.G. Cai, Y .P . Hu, Q.Y . Pan, Y .L. Zhang, Phys. Rev. D 91 (2015) 024032.
  • (51) E. Babichev and A. Fabbri, Phys. Rev. D 90 (2014) 084019.
  • (52) M. Chabab, H. E.Moumni, S. Iraoui and K. Masmar, Eur. Phys. J. C 79 (2019) 342.
  • (53) A. Dehghani and S.H. Hendi, Phys. Rev. D 104 (2021) 024025.