Maxwell's equal-area law for Gauss-Bonnet Anti-de Sitter black holes

Interpreting the cosmological constant \Lambda as a thermodynamic pressure and its conjugate quantity as a thermodynamic volume, we study the Maxwell's equal area law of higher dimensional Gauss-Bonnet-AdS black holes in extended phase space. These black hole solutions critically behave like Van der Waals systems. It has been realized that below the critical temperature T_c the stable equilibrium is violated. We show through calculations that the critical behaviors for the uncharged black holes only appear in d=5. For the charged case, we analyse solutions in d = 5 and d = 6 separately and find that, up to some constrains, the critical behaviors only appear in the spherical topology. Using the Maxwell's construction, we also find the isobar line for which the liquid-gas-like phases coexist.


Introduction
Recently, the study of thermodynamical properties of the black holes using techniques explored in statistical physics and fluids has received a special interest [1,2,3]. These researches have brought new understanding of the fundamental physics associated with the critical behaviors of several black holes in various dimensions using either numerical or analytic methods [1,4,5,6,7,8,9]. In particular, Ads black holes in arbitrary dimensions have been extensively investigated in many works [9,10,11,12,13,14,15,16,17,18,19,20,21,22].
More precisely, the state equations P = P(T, v) have been established by considering the cosmological constant as the thermodynamic pressure and its conjugate as the thermodynamic volume. In fact, this issue has opened a new way to study the behavior of the RN-AdS black hole systems using the physics of Van der Waals fluids. Indeed, it has been shown that the corresponding P-V criticality can be linked to the liquid-gas systems of statistical physics.
More recently, a special interest has been devoted to study of Maxwell's equal for some black hole solutions. More precisely, the state equations of various AdS black hole in P − v diagrams have been worked out showing the existence of undesirable negative pressure.
They appear also thermodynamic unstable regions associated with the condition ∂P ∂v > 0. In this way, the corresponding system can be contracted and expanded automatically [18,19,20,23,24,25]. In van der Waals systems, such problems have been overcame by using the Maxwell equal area approach.
The aim of this work is to contribute to these researches by studying higher dimensional uncharged and charged Gauss-Bonnet-Anti-de Sitter black holes in extended phase space.
Considering the cosmological constant Λ as a thermodynamic pressure and its conjugate quantity as a thermodynamic volume, we present the Maxwell's equal area law of higher Gauss-Bonnet-AdS black holes in extended space. These black hole solutions involve critical behaviors like Van der Waals gas. It has been shown that below the critical temperature T c the stable equilibrium is violated. Based on the Maxwell's construction, we obtain the isobar line for which the liquid-gas phases coexist.
The paper is organized as follow. In section 2, we give an overview on the thermodynamics of higher dimensional Gauss-Bonnet black holes in AdS geometry. In section 3, we first study the equal area law of Gauss-Bonnet-AdS black hole in extended phase space for uncharged solutions. Next, we extend the analysis to the charged case in five dimensions.
Then, we present the results for higher dimensions. The last section is devoted to conclusion.

Thermodynamics of Gauss-Bonnet black holes in AdS space
In this section, we give an overview on thermodynamics of Gauss-Bonnet black holes in AdS space. This matter is based on many works given by other authors [16,26,27,28]. Indeed, we start by considering a d-dimensional Einstein-Maxwell theory in the presence of the Gauss-Bonnet terms and a cosmological constant . The corresponding action reads where α GB is the Gauss-Bonnet coefficient with dimension [length] 2 . In this action, F µν is the Maxwell field strength given by Roughly speaking, the discussion will be given here corresponds to the case with a positive Gauss-Bonnet coefficient, namely, α GB ≥ 0. As shown in [26], dynamical solutions appear only in higher dimensional theories (d ≥ 5). For this reason, it should be interesting to concentrate on such models. In this way, the above action produces a static black hole solution with the following metric where h ij dx i dx j is the line element of a (d − 2)-dimensional maximal symmetric Einstein space with the constant curvature (d − 2)(d − 3)k and volume Σ k . It is noted that k takes three values 1, 0 and −1, corresponding to the spherical, Ricci flat and hyperbolic topology of the black hole horizon, respectively. According to [22,27,29,30], the metric function f takes the following form where α = (d − 3)(d − 4)α GB . In this solution, M represents the black hole mass, Q is linked to the charge of the black hole and P = − Λ 8π . For well-defined vacuum solution associated with M = Q = 0, the effective Gauss-Bonnet coefficient α and pressure P must satisfy the following constraint It is recalled that the mass M can be given in terms of the horizon radius r h of the black hole determined by the largest real root of the equation f (r h ) = 0. It takes the following form Similarly, the Hawking temperature of the black hole reads as It is worth noting that in the discussion of the thermodynamics of the black hole in the extended phase space by considering the pressure P = − Λ 8π , the black hole mass M should be identified with the enthalpy H ≡ M rather than the internal energy of the gravitational system [31]. It turns out that many other thermodynamic quantities can be obtained using thermodynamical equations. For instance, the entropy S, thermodynamic volume V and electric potential (chemical potential) Φ take the following forms These thermodynamic quantities satisfy the following differential form where is the conjugate quantity to the Gauss-Bonnet coefficient α, being considered as a variable.
By the scaling argument, the generalized Smarr relation for the black holes can be written as Moreover, many solutions treating other black holes are given in [13,31]. This class of black hole solutions has volume V ∼ Σ k r d−1 h and an area given by A ∼ Σ k r d−2 h . In this way, the black hole horizon has a scalar curvature R h ∼ k/r 2 h . The Gauss-Bonnet term on the horizon is R GB ∼ k 2 /r 4 h . In fact, the first term in (11) can be put in the form VR GB , while the second term takes the form TAR h . It is observed that both terms vanish in the flat solution corresponding to k = 0. It is noted that the first term in (11) is nothing but the second term in the black hole mass (5), while the second term in (11) is identified with the second term of the black hole entropy (7) multiplied by Hawking temperature T. Using Legendre transformations, the Gibbs free energy and the Helmholtz free energy read as However, it is recalled that the Helmholtz free energy F can be obtained by removing the contribution of the background of the AdS vacuum solution. In this way, only the situation associated with a negative F can be regarded as the black hole solution considered as a thermodynamically AdS vacuum solution. In the case associated with the hyperbolic horizon (k = −1) [32], the black hole entropy (7) could be negative. In fact, a negative entropy has no meaning in statistical physics. For these reasons, the constaints should be added. It is noted that in the case k = −1, the metric function f (3) shows the existence of a minimal horizon radius r 2 h ≥ 2α. However, the non-negative values of the black hole entropy (7) produce a more strong constraint on the horizon radius given by

The equal area law of Gauss-Bonnet-AdS black hole in extended phase space.
Having given the essentials on the thermodynamics of the Gauss-Bonnet-Anti-de Sitter black holes, we move to investigate the corresponding Maxwell's equal area law.

The construction of equal area law in P − v diagram
In the study of der Waals system, for an isotherm below the critical temperature T < T c the two points, of the plan (P, v), solving the equation indicate the the stability limit of the system. It is observed that the critical point corresponding to the largest volume is interpreted as the stability limit of the gaseous phase. However, the critical point associated with the smallest volume corresponds to the stability limit of the liquid phase. The chemical potential of such thermodynamic systems should satisfy In the isotherm transformation, the difference of chemical potential between two states with the pressure P and P 0 should have the following form It is realized that the Gauss-Bonnet-Ads black holes present similar behaviors as illustrated in figure 1. It observed from this figure that, at the point E, the black hole lies at "gas" phase. However, at point "A", the black hole lies at "liquid" phase completely. Moreover, the region between A and E can be considered as a coexistence phase. However, the oscillating part of the curve between A and E cannot be the coexistence line due to the fact that the part BD violates the equilibrium conditions. At A and E, the chemical potentials µ A (T, P) and This relation produces the thermodynamic condition for the phase equilibrium. Indeed, the equation (17) gives showing that the area EDC and ABC are equal.
In what follows, the main aim is to investigate the position of the points A and E for Gauss-Bonnet-AdS black holes in higher dimensions. First, we discuss the uncharged case then we will present numerical solutions for the charged black holes in the next section.

Uncharged solutions
In this subsection, we discus the neutral case corresponding to Q = 0. Identifying the pressure with the cosmological constant and the corresponding conjugate thermodynamic volume, the Hawking temperature (6) can be explored to write down the state equation. The latter reads as To make contact with the van der Waals equation, we should use a series expansion with the inverse of specific volume v. The calculation leads to Identify the specific volume v with the horizon radius of the black holes as proposed in [13,9,21] the specific volume v can be related to the horizon radius r h . Using (22), one can recover all the results in terms of v. In this way, the Eq. (20) becomes To obtain the equal area isobar, P = P 0 , one use the following relations: and Pdv.
The equal area law requires the following equality leading to In the isothermal curves, the points v 1 and v 2 should satisfy From these equations, one can derive the relations and Putting x = v 1 v 2 ,(0 ≤ x ≤ 1), we get the following relations and Using Eqs (32),(33) and(34), we get a polynomial equation admitting v 2 2 as a real positive root where the coefficients a, b and c read as Solving the above polynomial equation, we get In the limit In what follows, we discuss the convergence of this limit in terms of the topology and the dimension of the space-time. In fact, three situations can appear. They are classified as follows: • In the case of k = 0 corresponding to flat topology, this limit diverges for any dimensions d. This shows the absence of the critical points. Indeed, this can be seen from the equation of state eq. (23). The letter reduces to P = T v indicating that no phase transition can happen.
• In the case of k = −1 associated with the hyperbolic topology, a negative value of v 2 2 appears in d = 5. For d ≥ 6, the limit also diverges. This implies that there does not exist any phase transition.
• In the case of k = 1 corresponding to the spherical topology, this limit diverges only when d ≥ 6. In the case of d = 5, we have It is noted that a similar result has been found in [26]. Substituting the (38) in (33) to eliminate v 2 and setting T 0 = χT c , the critical temperature reads as We can also obtain the following relation In the table (1), we present the numerical values of the v 1,2 and P 0 for different values of the constant α.   For different values of α, we observe that the isobar in the isotherm will be short by increasing the temperature. When the temperature reaches the critical one, the boundaries of the isobar coincide v 1 = v 2 = v c .

Charged black hole solutions
This section concerns the case of charged Gauss-Bonnet-AdS black holes in the AdS geometry.
In this case, the state equation can be written as A close inspection on this equation shows that one may give a separate study. First, we deal with the five dimensional case. Then, we study models associated with d ≥ 6.

Five dimensional case
In d = 5, the equation of the state (40) reduces to Using a similar analysis used in the case of uncharged solutions, the equations (32), (33) and (34) become respectively and Using the above equations, we obtain also the polynomial equation of v 2 . It is given by where the coefficients a, b and c take now the following form The limit This shares similarities with the one obtained in [26]. The critical specific volume, which is the positive real root of (50), takes the following form The corresponding critical temperature reads as As in the uncharged case, we discuss the existence of the critical points for different topologies. In fact, we have three situations: • for the flat topology, the apparition of k, in fourth term in (51), reveals the absence of the critical behavior.
• for the hyperbolic one, the constraint on the positivity of the temperature and the specific volume does not show critical behavior.
• for the spherical topology, the existence of the critical point is controlled by the following constraint The Maxwell's construction can be obtained using the same analysis of the previous section.
This gives rise to an equation depend only on x. To derive the values of v 1,2 and the P 0 , we should find x. For the charged case, we list all these results on the table (2).
showing that the charge and α can be interrelated [26]. In what follows, we plot the isotherms in the (P, v) diagram and show the equal Maxwell's area.

Higher dimensional cases
The complete form of the equation of states (23) gives also a polynomial equation of v 2 .
Similar calculations can be done for higher dimensional black holes. In the limit x → 1, the specific volume is the real positive solution with the following polynomial form The critical temperature reads as It is observed that the Ricci flat topology and the hyperbolic one do not allow the existence of critical behaviors. Up to some conditions, the critical behaviors appear only in the case of the spherical topology corresponding to k = 1. In fact, α|Q| should not be too large [26].
For simplicity reasons, we restrict our study to d = 6. In particular, we present the numerical results in d = 6 for x, v 1,2 and P 0 in the following table, (1, 1)  Table 3: Numerical values of x, v 1 , v 2 and P 0 at different values of α in six dimensional spherical topology with non vanishing charge.
It is observed from the table (3) that P 0 increases with Q and χ. However, the specific volume v 2 decreases with Q and χ. The comparison between table (2) and (3) shows the effect of the black hole dimensions. In fact, the P 0 and v 2 increases and decreases resepectively with the dimension. To illustrate this effect, we plot these results in the fig.(4) showing the Maxwell equal area. Figure 4: The P − v diagram of the charged Gauss-Bonnet-AdS black holes in six dimensions, the dashed blue curve corresponds to the critical temperature T c , the red one corresponds to isotherm with 0.8T c and the green one is asscociated with 0.7t c .

Conclusion
In this paper, we have studied the the Maxwell's equal area law of higher dimensional Gauss-BonnetAnti-de-Sitter black holes. The corresponding critical behaviors share similarity with van der Waals one. We have shown that this construction can be used to eliminate the region of the violated stable equilibrium ∂P ∂v > 0. In particular, we have found the isobar line in which the two real phases coexist. It has been realized that this construction can be viewed as a simple way to drive the coordinates of the critical points. We have presented numerical calculations showing that the critical behaviors for the uncharged black holes appear only when d = 5. For the charged case, we have studied solutions in d = 5 and d = 6 separately and showed that, up to some constrains, the critical behaviors appear only in the spherical topology.