Thermodynamics and Phase Transitions of Dyonic AdS Black Holes in Gauss-Bonnet-Scalar Gravity

: In this paper, by treating the cosmological constant as a thermodynamic pressure, we study the thermodynamics and phase transitions of the dyonic AdS black holes in Gauss-Bonnet-Scalar gravity, where the conformal scalar field is considered. In a more general extended phase space, we first verified the first law of black hole thermodynamics, and find that it is always true. Meanwhile, the corresponding Smarr relation is also obtained. Then, we found that this black hole exhibits interesting critical behaviors in six dimensions, i.e., two swallowtails can be observed simultaneously. Interestingly, in a specific parameter space, we observed the small/intermediate/large black hole phase transitions, with the triple point naturally appearing. Additionally, the small/large black hole phase transition, similar to the liquid/gas phase transition of the van der Waals fluids, can also be found in other parameter regions. Moreover, we note that the novel phase structure composed of two separate coexistence curves discovered in the dyonic AdS black holes in Einstein-Born-Infeld gravity disappears in Gauss-Bonnet-Scalar gravity. This suggests that this novel phase structure may be related to gravity theory, and importantly, it is generally observed that the triple point is a universal property of dyonic AdS black holes. On the other hand, we calculated the critical exponents near the critical points and found that they share the same values as in mean field theory. Finally, it is true that these results will provide some deep insights into the interesting thermodynamic properties of the dyonic AdS black holes in the background of conformal scalar fields.


Introduction
The thermodynamics of black holes (BHs) constitutes an important research subject in the field of BH studies.In the 1970s, it was realized that BH is a thermodynamic system with temperature and entropy [1][2][3].The intriguing thermodynamic properties of BHs are gradually being discovered, which differ from those of ordinary thermodynamic systems.In particular, Hawking and Page discovered a phase transition between stable BHs and thermal radiation, which is referred to as the Hawking-Page phase transition [4].On the other hand, the anti-de Sitter/conformal field theory (AdS/CFT) correspondence indicates that the thermodynamics of BHs in AdS space can correspond to the thermodynamics of the dual strongly coupled conformal field theory in the boundary of AdS space [5][6][7].Therefore, in the framework of the AdS/CFT correspondence, the Hawking-Page phase transition is interpreted as a confinement/deconfinement phase transition of gauge fields [8].Motivated by this inspiration, the thermodynamics and phase transitions of BHs have been widely investigated [9][10][11][12][13][14].
Recently, there has been increasing attention directed towards studying the thermodynamics of AdS BHs in the extended phase space.It is worth noting that, in this framework, the negative cosmological constant Λ is interpreted as the thermodynamic pressure, while its conjugate quantity is regarded as the thermodynamic volume [15][16][17][18][19][20][21].Subsequently, Kubiz ňák and Mann compared the charged AdS BHs with the van der Waals (vdW) liquid-gas system, and demonstrated that they share the same oscillatory behavior of pressure-volume, critical exponents, and scaling relations [22].Thus, an analogy between charged AdS BHs and vdW systems has been established.Specifically, the small/large BH phase transition in AdS BHs is similar to the liquid/gas phase transition in vdW fluids.This analogy was applied to different types of AdS BHs, and suggested that such a small/large BH phase transition widely existed [23][24][25][26][27][28].Later, it was discovered that AdS BHs exhibit more interesting phase transitions in the extended phase space, such as the reentrant phase transition and the triple point [29][30][31][32][33][34][35][36][37][38][39][40][41].Furthermore, in the context of higher-order Lovelock gravity, there are also some intriguing BHs phase transitions, such as multicritical phase transitions [42][43][44].
By adding higher-order curvature terms to the gravitational action, the extended gravity theories have been established, which include Lovelock gravity and GB gravity [45][46][47].These extended gravity theories have provided new insights into the study of BHs in higherdimensional cases.For example, Wei found that charged AdS BHs in GB gravity exhibits small/intermediate/large BH phase transitions in six-dimensional spacetime, with the appearance of triple point [31].Frassino studied the phase transitions of charged AdS BHs in third-order Lovelock gravity and confirmed the existence of triple point and reentrant phase transitions in higher dimensions [32].More studies have shown that AdS BHs in extended gravity exhibits many intriguing properties [42][43][44][48][49][50].On the other hand, quasitopological electromagnetism is a class of interesting and significant objects in the study of BHs [51][52][53][54][55][56].Recently, Liu et al. have introduced a novel concept of quasitopological electromagnetism, which is defined as the square of the norm of the topological wedge product of the Maxwell field strength of order k (k > 2) [51].Subsequently, Li investigated the phase transitions of the dyonic AdS BHs in Einstein-Born-Infeld (EBI) gravity and obtained some interesting results, such as triple point and novel phase structure composed of two separate coexistence curves [57].Moreover, we also revealed the intriguing thermodynamic properties of dyonic AdS BHs in Einstein-Gauss-Bonnet (EGB) gravity and observed the triple point [58].
Recently, a class of high-dimensional dyonic BH solutions, which include the dyonic BHs in Gauss-Bonnet-Scalar (GBS) gravity, has been derived by coupling Lovelock-Scalar gravity with quasitopological electromagnetism [59].The conformal scalar field is very important for BHs and it is also quite interesting.Recently, Oliva et al. have developed a model for the gravity theory coupled to the real scalar field [60,61].Based on this model, higher dimensional BHs with scalar hair have been studied [62][63][64][65][66].It is believed that the conformal scalar field can affect the thermodynamic stability of hairy BHs [62,63].Interestingly, one also finds that the conformal scalar field has an impact on the local stability of BHs [59].Specifically, the range of horizon radius for stable BHs decreases with the increase in the conformal scalar field parameter H. Clearly, these studies indicate that BHs influenced by the conformal scalar field exhibit many intriguing thermodynamic properties which are worthy of further investigation.In addition, although dyonic AdS BHs in EBI gravity and EGB gravity exhibit some rich phase transitions, such as triple points [57,58], it is still unclear whether these intriguing phase transitions also exist in GBS gravity when considering the conformal scalar field.Therefore, in this paper, we study the thermodynamics and phase transitions of the dyonic BHs in GBS gravity by considering the effect of the conformal scalar field in the extended phase space.We are aim to further reveal the interesting thermodynamic properties of the dyonic BHs in GBS gravity and provide insights into the influence of the conformal scalar field on dyonic AdS BHs.
This paper is organized as follows.In Section 2, we give a review of the dyonic BHs in GBS gravity when considering the conformal scalar field.In Section 3, we study the thermodynamics of the dyonic BHs in the extended phase space.In Section 4, we investigate the phase transitions and phase diagrams of the dyonic AdS BHs, where the conformal scalar field is considered.Section 5 involved the computation of critical exponents near the critical points.Finally, Section 6 concludes with a summary and discussion.

Review of the Dyonic BHs
The action of high curvature gravity coupled to the conformal scalar field and matter sources can be expressed as [59] where L qt represents the Lagrangian density of matter, a p is a Lovelock coupling constant, and b p is a conformal coupling constant.Here, δ , where ξ describes the scalar field [60,61].As a certain class of Lovelock coupling gravity, the dimensionally continued gravity can be obtained when assuming [67][68][69] (2) It should be noted that l is related to the cosmological constant, i.e., Λ = − (d−1)(d−2) . Meanwhile, n is associated with the dimensionality parameter d, where d = 2n + 1 is for odd dimensions and d = 2n + 2 is for even dimensions.The gravitational field equation in Equation ( 1) can be determined by the action principle, i.e., where ς is the energy-momentum tensor corresponding to the material source, and which is defined by One should note that I qt represents the action associated with matter.And, the energymomentum tensor of the conformal scalar field is Based on the action principle, the scalar field equation is constructed as Moreover, the Lagrangian density of the material source is where η is a coupling parameter, and L int takes the form of In addition, the energy-momentum tensor corresponding to the Lagrangian density (7) can be expressed as To construct the hairy dyonic BH solutions, one applies condition (2) and the static spherical symmetry line element where dΥ 2 d−2 is the metric of the (d − 2)-dimensional hyper-surface of curvature (d − 2)(d − 3)γ, which is defined as The magnetic field of this (d − 2)-dimensional hyper-surface can be expressed as where q is related to the magnetic charge and Σ represents the volume of the (d − 2)dimensional hyper-surface.For the purely electric case, the corresponding Maxwell tensor is of the form where prime denotes differentiation with respect to r.Based on Equations ( 12) and ( 13), the following equation can be obtained, Integrating the above equation, it gives It should be noted that the constant of integration Q in the above equation is related to the electric charge.The configuration of a scalar field can be defined as When the conditions and are satisfied, it shows that ξ(r) is a solution of Equation (6).It can be shown that there is an unknown X in Equations ( 17) and ( 18), so that one of these two equations must constitute a constraint on the constant b p 's.Therefore, by considering the gravitational field equations given in Equation (3), choosing the Lovelock parameter a p arbitrarily, and based on the scalar field ξ(r) subject to the constraints of Equations ( 17) and ( 18), the energy-momentum tensor (9), and (15), an independent equation of motion can be obtained, which takes the form of where , with the conditions p ≥ 2. By solving Equation (19), one obtains where is a parameter related to the conformal scalar field.Therefore, by setting α 0,1,2 ̸ = 0 and α p = 0 for p ≥ 3 in Equation ( 21), a spherical (γ = 1) 1 dyonic BHs with the effect of conformal scalar field can be obtained in GBS gravity, which is [59] where is a hypergeometric function.Here, M is the BH mass, H is the parameter related to the conformal scalar field, Σ is the volume of the (d − 2)-dimensional hyper-surface, and l corresponds to the negative cosmological constant as . Moreover, Q and q represent the electric and magnetic charges of the BH, respectively.

Thermodynamics of the Dyonic BHs
In this subsection, we study the thermodynamics of the dyonic BHs in the extended phase space, where the conformal scalar field is considered.In this case, the negative cosmological constant Λ is regarded as the thermodynamic pressure P = − Λ 8π [15].Furthermore, by solving the equation f (r h ) = 0, the outer horizon radius r h of BH can be obtained, which is determined by the largest root of this equation.Thus, we can express the BH mass in terms of the horizon radius r h as Based on the definition of Hawking temperature T = f ′ (r h ) 4π , the BH temperature can be determined, Furthermore, in the extended phase space, the BH mass should be regarded as the enthalpy rather than the internal energy, i.e., H ≡ M. Therefore, we can calculate the other thermodynamic quantities of BH, such as the thermodynamic volume V, entropy S, electric potential Φ Q , and magnetic potential Φ q as follows, Considering the characteristics of the parameters H, α 2 , and η, we regard them as the new thermodynamic variables.Therefore, it can be verified that these thermodynamic quantities satisfy the follow differential form where is the conjugate quantity to H, is the conjugate quantity to α 2 , is the conjugate quantity to η.In addition, the corresponding Smarr relation can be obtained, which is of the form Through the above discussion, we find that both the first law of BH thermodynamics and the Smarr relation hold in a more general extended phase space which H is considered as the new thermodynamic variable.
The Gibbs free energy, as a quantity describing the global stability of a BH system, is given by G = H − TS.Therefore, we obtain The appearance of a swallowtail in a Gibbs free energy temperature (G − T) diagram indicates the occurrence of a BH phase transition.Therefore, we study BH phase transitions by analyzing the swallowtail observed in the G − T diagram.

Phase Transitions and Phase Diagrams of the Dyonic BHs
In this section, we would like to study the phase transitions and phase diagrams of the dyonic AdS BHs in GBS gravity.In particular, this paper focuses on the BH phase transition in six-dimensional spacetime 2 .Naturally, based on the temperature Equation ( 26), the equation of state for the BH can be obtained as Moreover, considering the thermodynamic volume V ∝ r d−1 h , the critical point can be determined by the following conditions, As is well known, the local thermodynamic stability of BHs is measured by the heat capacity with a constant pressure, denoted as C P .When C P is positive, it indicates that the system is locally stable.On the contrary, when C P is negative, it suggests that the system is locally unstable.For a constant pressure P, the heat capacity C P is defined as By substituting thermodynamic quantities into the calculation, C P can be further represented as In this paper, we follow the conditions of T > 0 and S > 0. Thus, in the T − r h diagram, the BH branches with positive and negative slopes correspond to stable and unstable phases, respectively, and their corresponding C P values are positive and negative.
Next, we proceed to study BH phase transitions and phase diagrams, with a particular focus on the triple point.In addition, we also focus on the effect of the conformal scalar field on BH phase transitions.Therefore, in this paper, we set the electric charge Q = 10, magnetic charge q = 5, and parameter Σ = 1, and vary the conformal scalar field parameter H, the coupling parameters α 2 , and η to investigate the phase transitions and phase diagrams of the dyonic BHs.

Phase Transitions by Fixing H and η While Varying α 2
In this subsection, we study the BH phase transitions and phase diagrams by fixing the conformal scalar field parameter H = 0.01 and the coupling parameter η = 0.01, while varying α 2 as 1, 5, and 10.

α 2 = 1
In this case, according to Equations ( 37) and ( 38), a critical point can be determined, which is The behaviors of temperature T with respect to horizon radius r h and Gibbs free energy G with respect to temperature T are plotted in Figure 1a,b, respectively.It should be noted that isobaric curves of the same color in (a) and (b) correspond to same pressure values.In Figure 1a, when P < P c , two extremal points appear on the isobaric curves (represented by red and green), which divide these two curves into three branches: the stable small BH branch, the unstable intermediate BH branch, and the stable large BH branch.Among them, stable branches are represented by solid curves with a positive value of C P , while unstable branches are represented by dashed curves with a negative value of C P .When P > P c , there is no extremal point on the isobaric curve, and T increases monotonically with r h .Now, let us focus on the behavior of Gibbs free energy in Figure 1b.For P < P c , a swallowtail appears on each isobaric curve, which suggests the occurrence of a first-order small/large BH phase transition.It should be noted that the nonsmooth points on the isobaric curves in the G − T diagram correspond to the extremal points on the isobaric curves in the T − r h diagram.Regarding the red and green isobaric curves, it can be found that the system initially exhibits a small BH phase and turns into a large BH phase near the intersection as temperature increases.Furthermore, by comparing the red with the green isobaric curves, it can be discovered that the size of the swallowtail decreases as the pressure increases.When the pressure reaches P c , the swallowtail disappears.For P > P c , the Gibbs free energy decreased monotonically with temperature, which indicates that no phase transition occurs in the system.The phase diagram for the dyonic BHs, as shown in Figure 1c.It can be observed that the pressure increases monotonically with temperature and terminates at the critical point (P c , T c ).The region of small BHs is located above the coexistence curve, while the region of large BHs lies below it.This is a typical small/large BH phase transition, which is similar to the vdW liquid/gas phase transition.
Firstly, we focus on the behavior of temperature, as shown in Figure 2a.When P < P c1 , two extremal points appear on the blue isobaric curve, which divide it into three branches: the stable small BH branch, the unstable intermediate BH branch, and the stable large BH branch.Interestingly, when P c1 < P < P c2 , four extremal points emerged on the red isobaric curve, which divide it into five branches: the stable small BH branch, the unstable small BH branch, the stable intermediate BH branch, the unstable large BH branch, and the stable large BH branch.For P = P t = 0.00337487, it is easy to utilize Maxwell equal area laws to construct two pairs of equal area regions in the T − S diagram, as illustrated in Figure 2b.These two pairs of regions have the same temperature, i.e., T = T t = 0.0485154.This implies that the BH undergoes two phase transitions simultaneously at such pressure and temperature.In fact, this result indicates the existence of a triple point, where small, intermediate, and large BH phases can coexist.When P increases to P c2 , the BH system undergoes a second-order BH phase transition.For P c2 < P < P c3 , two extremal points on the orange isobaric curve divide it into three branches: the stable small BH branch, the unstable intermediate BH branch, and the stable intermediate BH branch.Moreover, the first-order phase transition turns to a second-order phase transition as P approaches to P c3 .When P > P c3 , the temperature T increases monotonically with r h , which implies that there is only one BH branch.We have plotted the behavior of Gibbs free energy G with respect to temperature T, as shown in Figure 3.When P < P c1 , a swallowtail appears in the G − T diagram, which indicates a small/large BH phase transition.As the pressure increased to P c1 < P < P t , two swallowtails appear in Figure 3b, which suggests the potential existence of two BH phase transitions.However, the intermediate BH branch is suppressed by the lower free energy branch and does not participate in the phase transition.Therefore, in this case, only the small/large BH phase transition occurs.Increasing the pressure to P = P t , it can be discovered that the intersection points of two swallowtails appear at the same point, which suggests the occurrence of small/intermediate/large BH phase transitions.For P t < P < P c2 , two swallowtails can be observed, and all BH branches can participate in the phase transitions.This indicates that in this case, the system can undergo small/intermediate and intermediate/large BH phase transitions simultaneously.When the pressure is increased to P c2 < P < P c3 , one swallowtail can be noticed, which implies that there exists only one small/intermediate BH phase transition.When P > P c3 , the Gibbs free energy decreases monotonically with temperature, which indicates that no phase transition occurs in the system.
The phase diagram is illustrated in Figure 4, where Figure 4b is a local magnification near the triple point.As can be seen from the phase diagram, the system undergoes the small/large BH phase transition at P < P t .Interestingly, at the triple point (P t , T t ), the small, intermediate, and large BHs phases coexist, and the system undergoes the small/intermediate/large BH phase transitions.When P t < P < P c2 , the small/intermediate and intermediate/large BH phase transitions appear simultaneously.When the pressure increases to P c2 < P < P c3 , the system undergoes only the small/intermediate BH phase transition.In summary, these results suggest the presence of a rich variety of phase transition types in this parameter region.
Firstly, we would like to analyze the behavior of temperature T with respect to r h , as shown in Figure 5.When P < P c1 , two extremal points appear on the blue isobaric curve, which divide it into three branches: the stable small BH branch, the unstable intermediate BH branch, and the stable large BH branch.As the pressure increases to P c1 < P < P c2 , the red isobaric curve is divided into five branches by four extremal points: the stable small BH branch, the unstable small BH branch, the stable intermediate BH branch, the unstable large BH branch, and the stable large BH branch.For P c2 < P < P c3 , there are two extremal points on the orange isobaric curve, which suggests the existence of three BH branches.When the pressure increases to P > P c3 , the temperature increases monotonically with r h , which implies that there is only one BH branch.Next, let us analyze the behavior of Gibbs free energy G with respect to T. When P < P c1 , the appearance of a swallowtail in Figure 6a actually indicates the small/large BH phase transition.When the pressure increased to P c1 < P < P c2 , two swallowtails can be observed in Figure 6b.However, the intermediate BH branch is suppressed by the BH branches with lower Gibbs free energy and does not participate in the phase transition.Therefore, in this pressure range, only the small/large BH phase transition occurs.For P c2 < P < P c3 , the presence of a swallowtail indicates the occurrence of the small/large BH phase transition.When P > P c3 , the monotonic decrease in Gibbs free energy suggests that no phase transition occurs in the system.Finally, the phase diagram is illustrated as shown in Figure 7.The coexistence curve originates from the origin and terminates at the critical point (P c3 , T c3 ).From the phase diagram, it is evident that while there are three critical points in this parameter region, only the small/large BH phase transitions occur.

Phase Transitions by Fixing H and α 2 While Varying η
In this subsection, by setting the conformal scalar field parameter H = 0.01 and the coupling parameter α 2 = 6, while varying η as 0.01, 0.1, and 5, we study the BH phase transitions and phase diagrams.
In this case, three critical points can be obtained, which are T c2 = 0.0455153, P c2 = 0.00317372, ( 49) We have plotted the behavior of temperature T with respect to r h , as shown in Figure 8a.For P c1 < P < P c2 , there are four extremal points on the isobaric curve, which suggests the existence of five BH branches.When P = P t = 0.00292831, we have utilized Maxwell equal area law to construct two pairs of equal area regions with the same temperature T = T t = 0.0447573, as shown in Figure 8b.In fact, this predicts the existence of a triple point.Therefore, we plot the behavior of Gibbs free energy with respect to temperature, as shown in Figure 9.As anticipated, the three BH branches intersect at a point when P = P t , which indicates the occurrence of the small/intermediate/large BH phase transitions.Finally, the phase diagram is plotted in Figure 10.In the phase diagram, the triple point where small, intermediate, and large BHs coexist can be observed.4.2.2.η = 0.1 Like the previous case, three critical points can be obtained, which are T c2 = 0.0455649, P c2 = 0.00319072, T c3 = 0.0463009, P c3 = 0.00429228.
The T − r h diagram is illustrated in Figure 11, where four extremal points appear on the red isobaric curve at P = 0.00316.However, the intermediate BH is suppressed and does not participate in phase transitions, as shown in Figure 12b.This indicates that the system only undergoes the small/large BH phase transitions, similar to the case discussed in Section 4.1.3.Finally, the phase diagram, as illustrated in Figure 13, indicates that only the small/large BH phase transitions occur in this parameter region.
In Figure 14, there exists an isobaric curve with four extremal points, which seems to indicate complex phase transition behavior.However, by analyzing the behavior of Gibbs free energy in Figure 15, it can be found that only the small/large BH phase transitions occur, which is similar to the previous case.The phase diagram in Figure 16 also supports our analysis, which indicates the small/large BH phase transitions occur in this parameter region.
The behavior of temperature T with respect to r h is shown in Figure 17a.For P c1 < P < P c2 , four extremal points appear on each isobaric curve, which indicates a rich variety of phase transitions.Obviously, when P = P t = 0.00342036, we have utilized the Maxwell equal area law to construct two pairs of equal area regions, as shown in Figure 17b.Then, we plotted the behavior of Gibbs free energy in Figure 18.In particular, when P = P t = 0.00342036, it can be discovered that three BH branches intersect at one point, which suggests the occurrence of small/intermediate/large BH phase transitions.Finally, in the phase diagram shown in Figure 19, a triple point can be observed where small, intermediate, and large BH phases coexist.
The behaviors of temperature and Gibbs free energy are plotted in Figures 20 and 21.After careful analysis, it can be determined that only small/large BH phase transitions occur in this parameter region.Finally, the phase diagram in Figure 22 further shows this small/large BH phase transition.T vs. r h , where the blue, purple, red, black, orange, green, and gray curves represent P = 0.003, 0.00326926, 0.0035, 0.00383222, 0.008, 0.0161884, and 0.02, respectively.The solid and dashed curves represent stable and unstable branches, respectively.

Critical Exponents
It is widely believed that the critical exponents offer a valuable method for describing the behavior of physical quantities near the critical point, and they do not depend on the details of the physical system.Therefore, in this section we would like to calculate the critical exponents in the vicinity of the critical points.
For convenience, we start by defining some reduced parameters, whereby we set where τ = T T c is defined as the reduced thermodynamic temperature, we set where ν = V V c is defined as the reduced thermodynamic volume, and we define as the reduced thermodynamic pressure.
Next, let us review the definitions of the critical exponents α, β, γ, and δ near the critical point [22]: (1) Exponent α determines the behavior of the specific heat at constant volume, (2) Exponent β describes the behavior of the order parameter η 1 = V l − V s (the difference between the volumes of the coexisting large and small BHs) on a given isotherm (3) Exponent γ governs the behavior of the isothermal compressibility κ T (4) Exponent δ reflected the following behavior on the critical isotherm T = T c Subsequently, we will proceed to compute these critical exponents.From Equation ( 28), it can be found that the entropy S is independent of T, which leads Thus, the critical exponent is α = 0.In addition, by substituting the reduced parameters introduced in Equations ( 63)-( 65) into Equation (37), we can calculate the corresponding state equations.Therefore, it is easy to express the reduced pressure near the critical point as The values of the expanded coefficients in Equation ( 71) for different parameters were calculated, and the corresponding results are listed in Table 1.As shown in Table 1, the coefficient A 0 = 1, the coefficient B 0 is positive, and the coefficients A 3 and B 1 are negative.Moreover, the coefficients A 1 and A 2 are absent in the BH system, and thus they have not been listed in this table.Therefore, the reduced pressure can be re-expressed as first law of BH thermodynamics and derive the corresponding Smarr relation in a more general extended phase space.Then, we study the BH phase transitions by analyzing the characteristic behaviors of temperature and Gibbs free energy in six-dimensional spacetime.
On the other hand, to analyze the effect of the conformal scalar field on the BH phase transition, we fix and vary the values of the conformal scalar field parameter H and the coupling parameters α 2 and η to study the BH phase transition.We first consider the case where H = 0.01 and η = 0.01, while α 2 takes the values of 1, 5, and 10.For α 2 = 1, we observed a typical small/large BH phase transition, which is similar to the vdW liquid/gas phase transition.When α 2 = 5, we discovered four extremal points on the isobaric curve in the T-r diagram, as well as two swallowtails in the G-T diagram.This actually indicates a rich variety of phase transitions beyond the small/large BH phase transition.As a result, the small/intermediate/large BH phase transitions can be found in this case.Additionally, the triple point where the small, intermediate, and large BHs can coexist is obtained, i.e., (P t = 0.00337487, T t = 0.0485154).For α 2 = 10, although there are three critical points, as the intermediate BH branch has a higher free energy and is suppressed by the lower free energy BH branch, thus only the small/large BH phase transition can be discovered.Then, we fixed H = 0.01 and α 2 = 6 while we varied the parameter η to study the phase transition of BHs.When η = 0.01, we observed the small/intermediate/large BH phase transitions, as well as the triple point, which is located at (P t = 0.00292831, T t = 0.0447573).For the case of η = 0.01 and 5, which is similar to the case of H = 0.01, η = 0.01, and α 2 = 10, only the small/large BH phase transition occurs.Finally, we set α 2 = 5 and η = 0.01 and varied the parameter H to study the BH phase transition.For H = 0.1, the small/intermediate/large BH phase transitions can be found.As expected, a triple point appears, which is located at (P t = 0.00342036, T t = 0.0486811).While for H = 1, only the small/large BH phase transitions can be observed.Based on the above discussions, it can be found that the conformal scalar field has some significant impact on the BH thermodynamics and phase transitions.On the other hand, we also note that the novel phase structure composed of two separate coexistence curves, discovered in EBI gravity [57], is absent in GBS gravity.This is consistent with the results obtained in EGB gravity [58], further suggesting that this novel phase structure is related to the gravity theory.In fact, these results also demonstrate that the triple point, where small, intermediate, and large BHs can coexist, is a universal feature of dyonic AdS BHs.
Moreover, we calculated the critical exponents near the critical points and obtained results for α = 0, β = 1 2 , γ = 1, and δ = 3.This implies that these critical exponents share the same values as in mean field theory and are consistent with those obtained in other BH systems.Finally, it can be summarized that our conclusions will provide important insights for a deep understanding of the intriguing thermodynamic properties of the dyonic AdS BHs in GBS gravity.
In addition, holographic duality provides a valuable method for the study of the thermodynamic properties of AdS BHs, i.e., it reveals these characteristics from an alternative perspective [70,71].Therefore, it is worth utilizing holographic duality to further investigate the thermodynamics of this dyonic AdS BHs, which may reveal more intriguing thermodynamic properties.This will also be a part of our future work.

Figure 2 .
Figure 2. (a) T vs. r h .The blue, red, black, orange, green and gray curves correspond to P = 0.0024, 0.00337487, 0.00374402, 0.004, 0.0044492 and 0.005, respectively.Solid curves indicate stable branches, while dashed curves represent unstable branches.(b) T − S diagram of two pairs of equal area regions at pressure P = P t = 0.00337487.The horizontal line has a temperature T = T t = 0.0485154.

Figure 6 .
Figure 6.G vs. T.The red, green and orange curves represent the small, intermediate and large BHs, while the dashed curves indicate unstable BHs.

Figure 8 .Figure 9 .
Figure 8.(a) T vs. r h .The blue, red, black, orange, green, and gray curves correspond to P = 0.002, 0.00292831, 0.00317372, 0.004, 0.00583728, and 0.007, respectively.The solid and dashed curves represent stable and unstable branches, respectively.(b) T − S diagram of two pairs of equal area regions at a pressure of P = P t = 0.00292831.The horizontal line has a temperature of T = T t = 0.0447573.

Figure 15 .
Figure 15.G vs. T.The red, green and orange curves represent the small, intermediate and large BHs, while the dashed curves indicate unstable BHs.

Figure 17 .
Figure 17.(a) T vs. r h .The blue, red, black, orange, green, and gray curves corresponds to P = 0.0025, 0.00342036, 0.00375132, 0.0043, 0.00494838, and 0.0055, respectively.The solid and dashed curves represent stable and unstable branches, respectively.(b) T − S diagram of two pairs of equal area regions at pressure P = P t = 0.00342036.The horizontal line has a temperature T = T t = 0.0486811.

Author Contributions:
Conceptualization, P.M. and G.L.; methodology, P.M.; software, P.M.; validation, P.M., Z.Y. and G.L.; formal analysis, P.M.; investigation, P.M.; resources, P.M.; data curation, P.M.; writing-original draft preparation, P.M.; writing-review and editing, P.M., Z.Y. and G.L.; visualization, P.M.; supervision, Z.Y. and G.L.; project administration, Z.Y. and G.L.; funding acquisition, Z.Y. and G.L. All authors have read and agreed to the published version of the manuscript.Funding: This work is supported by the National Natural Science Foundation of China (Grant No. 11903025), by the starting fund of the China West Normal University (Grant No. 18Q062), by the Sichuan Science and Technology Program (2023ZYD0023), by the Sichuan Youth Science and Technology Innovation Research Team (21CXTD0038), and by the Natural Science Foundation of Sichuan Province (2022NSFSC1833).
•••σ p λ p is the generalized Kronecker delta and R σλ µν represents the components of the Riemann tensor.The tensor S σλ µν is defined as Figure 18.G vs. T.The red, green, and orange curves represent the small, intermediate, and large BHs, while the dashed curves indicate unstable BHs.

Table 1 .
The values of the expanded coefficients in Equation (71) for different parameters.