Quintessence background for $5D$ Einstein-Gauss-Bonnet black holes

As we know that the Lovelock theory is an extension of the general relativity to the higher-dimensions, in this theory the first and the second order terms correspond to the general relativity and the Einstein-Gauss-Bonnet gravity, respectively. We obtain a $5D$ black hole solution in Einstein-Gauss-Bonnet gravity surrounded by the quintessence matter, and also analyze their thermodynamical properties. Owing to the quintessence corrected black hole, the thermodynamic quantities have also been corrected except for the black hole entropy, and phase transition is achievable. The phase transition for the thermodynamic stability, is characterized by a discontinuity in the specific heat at $r=r_C$, with the stable (unstable) branch for $r<(>) r_C$.


I. INTRODUCTION
The gravity theory with higher-curvature term, the socalled Lovelock is one of the natural generalization of Einstein's general relativity, introduced originally by Lanczos [1], and rediscovered by David Lovelock [2], the action of it contains higher-order curvature terms and that reduces to the Einstein-Hilbert action in four-dimensions, and its second-order term is the Gauss-Bonnet invariant. The Lovelock theories have some special characteristics, among the larger class of general higher-curvature theories, in having field equations involving not more than second derivatives of the metric. As higher-dimensional members of Einstein's general relativity family, the Lovelock theories allow us to explore several conceptual issues of gravity in a broader setup. Hence, these theories receive significant attention, especially when finding black hole solutions. Besides, the theory is known to be free of ghosts about other exact backgrounds [3]. The theory represents a very interesting scenario to study how higher-order curvature corrections to the black hole physics substantially changes the qualitative features, we know from our experience with black holes in general relativity. Since, its inception, steady attention has been devoted to black hole solutions, including their formation, stability, and thermodynamics. The spherically symmetric static black hole solution for the Einstein-Gauss-Bonnet theory was first obtained by Boulware and Deser [3], and later several authors explored exact black hole solutions and their thermodynamical properties [4][5][6][7][8][9]. The generalization of the Boulware-Desser solution has been obtained with a source as a cloud of strings, in Einstein-Gauss-Bonnet gravity [10,11], and also in Lovelock grav-ity [12][13][14].
An intense activity of studying black hole solutions in Einstein-Gauss-Bonnet theory of gravity is due to the fact that, besides theoretical results, cosmological evidence, e.g., dark matter and dark energy. The quintessence is a hypothetical form of dark energy postulated as an explanation of the observation for an acceleration of the Universe, rather than due to a true cosmological constant. If quintessence exists all over in the Universe; it can also be around a black hole. In this paper, we are interested in a solution to the Einstein equations with the assumption of spherical symmetry, with the quintessence matter obtained by Kiselev [15] and was also rigorously analyzed by himself and others [15][16][17][18][19][20]. In particular, spherically symmetric quintessence black hole solutions [15], which is extended to the higher dimensions [21], to include Narai solutions [22,23], and also to the charged black holes [24]. The black hole thermodynamics for the quintessence corrected solutions were obtained in [25][26][27][28][29] and, quasinormal modes of such solutions are also discussed [30][31][32][33]. The generalization of the spherical quintessential solution [15] to the axially symmetric case, Kerr-like black hole, was also addressed, recently [34,35]. However, the solution of the Einstein-Gauss-Bonnet theory surrounded by the quintessence matter is still not explored, i.e., the black holes surrounded by the quintessence matter in Einstein-Gauss-Bonnet theory is still unknown. It is the purpose of this paper to obtain an exact new five-dimensional (5D) spherically symmetric black holes solution for the Einstein-Gauss-Bonnet surrounded by the quintessence matter. In particular, we explicitly bring out how the effect of a background quintessence matter can alter black hole solutions and their thermodynamics. In turn, we analyze their thermodynamical properties and perform a thermodynamic stability analysis.
The paper is organized as follows. In Sec. II, we derive a Einstein-Gauss-Bonnet solution to the 5D spherically symmetric static Einstein's equations surrounded by the quintessence matter. In Sec. III, we have discussed the thermodynamics of the 5D Einstein-Gauss-Bonnet black holes surrounded by the quintessence matter. The paper ends with concluding remarks in Sec. IV.
We use units which fix the speed of light and the gravitational constant via G = c = 1, and use the metric signature (−, +, +, +, +).

II. QUINTESSENCE MATTER SURROUNDING BLACK HOLE
The Lovelock theory is an extension of the general relativity to higher-dimensions. In this theory the first and second order terms correspond to the general relativity and the Einstein-Gauss-Bonnet, respectively. The action for 5D Einstein-Gauss-Bonnet gravity with matter field reads: with κ 5 = 1. I S denotes the action associated with matter and α is coupling constant that we assume to be non-negative. The Einstein term is L 1 = R, and the second-order Gauss-Bonnet term L GB is Here, R µν , R µνγδ , and R are the Ricci tensor, Riemann tensor, and Ricci scalar, respectively. The variation of the action with respect to the metric g µν gives the Einstein-Gauss-Bonnet equations where G E µν is the Einstein tensor while G GB µν is given explicitly by [36] and T S µν is the energy-momentum tensor of the matter that we consider as a quintessence matter. We note that the divergence of Einstein-Gauss-Bonnet tensor G GB µν vanishes. Here, we want to obtain 5D static spherically symmetric solutions of Eq. (3) surrounded by the quintessence matter and investigate its properties. We assume that the metric has the form [5,14] whereγ ij is the metric of a 3D constant curvature space k = −1, 0, or 1. In this paper, we shall restrict to k = 1. Using this metric ansatz, the Einstein-Gauss-Bonnet Eq. (3) reduces to The energy-momentum tensor of the quintessence matter (see Ref. [15] for further details) gets modified to where B(r) is a quintessential parameter and r a r b = 1 4 r n r n .
Thus, we have the equation of state of the form where for the quintessential matter −1 < ω < 0, which implies −4 < β < 0 within this set up, the parameter B of the energy-momentum tensor reads [21] Hence, the energy-momentum tensor for a quintessence matter takes the form where ρ is the proper density of the quintessence matter. The range of the parameters ω and β gets modified in 5 dimensions.
The Gauss-Bonnet term is the only possibility for the leading correction to Einstein general relativity. The static spherically symmetric black hole solution of Einstein action, modified by the Einstein-Gauss-Bonnet term, was first obtained by Boulware and Deser [3] and demonstrated that the only spherically symmetric solution to the Einstein-Gauss-Bonnet theory is Schwarzschild type solution. Later, Wiltshire [37] included the Maxwell field to the Einstein-Gauss-Bonnet action, found the charged black hole in this theory which was a generalization of the Reissner-Nordström black hole. In general, it is difficult to find solution of Einstein-Gauss-Bonnet field equations (3)  state. Here, we shall find a black hole solution surrounded by the quintessence matter with the equation of state (9). Making use of Eqs. (6) and (11), we deduce the master equation for the Einstein-Gauss-Bonnet gravity as In general, Eq. (12) has one real and two complex solutions. It may have three real solutions as well under some conditions. Here, we consider only the real solution. Interestingly, Eq. (12), for the Einstein-Gauss-Bonnet case, admits a general solution by appropriately relating M and q with integrating constants c 1 and c 2 [11]. Equation (13) is an exact solution of the field equation (12) for an equation of state (9) which in case of no quintessence ω = 0 reduces to Boulware and Deser [3] Gauss-Bonnet black hole solution, and for ω = 1/2 and q = −4Q 2 /3 to a solution mathematically similar to the the charged Gauss-Bonnet black hole due to Wiltshire [37]. When ω = −1, q = Λ/3, the Eq. (13) corresponds to Gauss-Bonnet de Sitter solution. In the limit α → 0, the negative branch of the solution (13) reduces to the general relativity solution. To study the general structure of the solution (13), we take the limit r → ∞ or M = q = 0 in solution (13) to obtain this means the plus (+) branch of the solution (13) is asymptotically de Sitter (anti-de Sitter) depending on the sign of α (±), whereas the minus branch of the solution (13) is asymptotically flat. In the large r limit (or α → 0), Eq. (13) reduces to the 5D Schwarzschild solution surrounded by the quintessence matter. Thus, the negative branch solution (13) is well behaved and it represents short distance correction to 5D black hole solution of general relativity. In similar way, when M = 0, the solution (13) takes the form Obviously, by proper choice of the functions M and q, and parameter ω, one can generate many other known solutions. The above solutions include most of the known spherically symmetric solutions of the Einstein-Gauss-Bonnet field equations (3).

III. THERMODYNAMICS
In this section, we shall discuss the thermodynamical properties of 5D quintessential black hole within Einstein-Gauss-Bonnet framework. Henceforth, we shall restrict ourselves to the negative branch of the solution (13). By definition of a horizon, the value of r = r + is an event horizon when f (r + ) = 0. This is shown in Fig. 1, by plotting f (r) as a function of r. It is interesting to note that the black holes admit only one horizon and the radius of the horizon increases with the value of the quintessence matter parameter ω. Next, we explore the thermodynamics of the black hole solution (13) surrounded by the quintessence matter in Einstein-Gauss-Bonnet framework. The Einstein-Gauss-Bonnet black holes surrounded by the quintessence matter are characterized by their mass (M ) and a quintessence matter parameter (ω). From Eq. (13), the mass of the black hole is obtained in terms of the horizon radius (r + ): The mass of the black hole is plotted in Fig. 2 for various values of parameters ω and α, which shows an increase in the black hole mass with horizon radius r + . To discuss the thermodynamics of the metric (5) with function (13), we start with the Hawking temperature. The Hawking temperature associated with the black hole is defined by T = κ/2π, where κ is the surface gravity defined by [9,14], Hence, the Hawking temperature for the Einstein-Gauss-Bonnet black hole surrounded by the quintessence matter can be calculated as Note that the factor in the numerator of Eq. (18) modifies the Gauss-Bonnet black hole temperature [9,14], and taking the limit q → 0, we recover the Gauss-Bonnet black hole temperature as and when α → 0, it becomes the temperature given by T + = 1 2πr+ . It is interesting to note that for a particular radius of the horizon, the Hawking temperature vanishes. The Hawking temperature diverges in general relativity as r + → 0. However, in Einstein-Gauss-Bonnet it remains finite as shown in Fig. 3. Also, when q = 0 and α = 0, the Hawking temperature has a peak which decreases and shifts as α increases or q increases (cf., Fig. 3).
A black hole behaves as a thermodynamical system; quantities associated with it must obey the first law of thermodynamics [7]: Hence, the entropy [38] can be obtained from the integration Now, the entropy of the Einstein-Gauss-Bonnet gravity black holes surrounded by the quintessence matter, reads However, it is interesting to note that the entropy of the black hole has no effect of the quintessence matter parameter.
Next, let us calculate the Wald entropy for the 5D black hole (5) with f (r) given by Eq. (13). Wald [39] showed that the black hole entropy can be calculated by where dV 3 is the volume element on Σ and the integral is performed on 3D space-like surface Σ. ε ab is the binormal vector to surface Σ normalized as ε ab ε ab = −2, and the L is Lagrangian density as in (1). We note that the integrand can be calculated as On substituting the Eq. (24) into the Eq. (23), one obtain Wald entropy of the 5D black hole (5) as as f − (r + ) = 0 and hence the Wald entropy Eq. (25) has exactly same expression as obtained in Eq. (22). Furthermore, we verify that the 5D black hole (5) satisfy the first law of thermodynamics. The variation of Wald entropy (25) with respect to the radius r + gives Hence, with the help of Eqs. (18), (26) and (27), one can conclude that dM EGB Hence the first law of the black hole thermodynamics holds for 5D black hole (5) with f (r) given by Eq. (13). Finally, we analyze how the quintessence matter influences the thermodynamic stability of the Einstein-Gauss- Bonnet black holes. The heat capacity of the black hole is defined as [7], The heat capacity of Einstein-Gauss-Bonnet black hole surrounded by the quintessence matter, using Eqs. (16), (18), and (29) is given by It is clear that the heat capacity depends on both the Gauss-Bonnet coefficient α, and a quintessence matter parameter ω. When α → 0, it returns to the general relativity case. If in addition q = 0, it becomes which is exactly same as the Einstein-Gauss-Bonnet case [9,14]. We again recall that for C > 0 (C < 0), the black hole is thermodynamically stable (unstable). It is difficult to analyze the heat capacity analytically hence, we plotted it in Fig. 4, for different values of α and ω. Again, we note that there is a change of sign in the heat capacity around r C , and C is discontinuous at r + = r C . The heat capacity is positive for r + < r C and thereby suggesting the thermodynamical stability of a black hole. On the other hand, the black hole is unstable for r + > r C . Thus, the heat capacity of an Einstein-Gauss-Bonnet black hole, for different values of ω and α, is positive for r + < r C , while for r + > r C , it is negative. The phase transition occurs from a lower mass black hole with the negative heat capacity to a higher mass black hole with positive heat capacity.
It may be noted that the critical radius r C changes drastically in the presence of the quintessence matter, thereby affecting the thermodynamical stability. Indeed, the value of critical r C increases with the increase in the quintessence matter parameter ω for a given value of the Gauss-Bonnet coupling constant α. Further, r C is also sensitive to the Gauss-Bonnet parameter α (cf. Fig. 4), and the critical parameter r C also increases with α.

IV. CONCLUSION
The Einstein-Gauss-Bonnet theory has a number of additional nice properties than Einstein's general relativity that are not enjoyed by other higher-curvature theories. Hence, Einstein-Gauss-Bonnet theory received significant attention, especially when finding black hole solutions. We have obtained an exact 5D static spherically symmetric black hole solutions to Einstein-Gauss-Bonnet gravity surrounded by the quintessence matter. We then proceeded to find exact expressions, in Einstein-Gauss-Bonnet gravity, for the thermodynamical quantities like the black hole mass, Hawking temperature, entropy, specific heat and in turn also analyzed the thermodynamical stability of black holes. It turns out that due to quintessence correction to the black hole solution, the thermodynamical quantities are also getting corrected except for the entropy which does not depend on the background quintessence. The entropy of a black hole in Einstein-Gauss-Bonnet gravity doesn't obeys the area law.
The phase transition is characterized by the divergence of specific heat at a critical radius r C which is changing with Gauss-Bonnet parameter α as well as with w. In particular, the black hole is thermodynamically stable, with a positive heat capacity for the range 0 < r < r C and unstable for r > r C . It would be important to understand how these black holes with positive specific heat (C > 0) would emerge from thermal radiation through a phase transition. We also discuss the phase transition of the black holes. The results presented here are the generalization of the previous discussions, on the Einstein-Gauss-Bonnet black hole, in a more general setting, and the possibility of a further generalization of these results to Lovelock gravity is an interesting problem for future research. On the other hand, the quintessence matter alone (M = 0) has horizon placed at r + = (q) 1/4ω+2 . Obviously, in the limit q → 0, the above solution will reduce to 5D general relativity black hole in which case R = R = 0. Next, we shall discuss the thermodynamics of the 5D black hole surrounded by the quintessence matter. We note that the gravitational mass of a black hole is determined by g rr (r + ) = 0, which, from Eq. (A2), reads (A7) Eq. (A7) takes the form of the 5D Schwarzschild black hole mass M = r 2 + , when q → 0 [5]. Accordingly, the Hawking temperature of the black hole at outer horizon r + , reads Then, we can easily see that the temperature is positive. Taking the limit q → 0, we recover the temperature for 5D general relativity T + = 1 2πr+ [9,14], which shows that Hawking temperature diverges as r + → 0. Next, we turn to calculate the Wald entropy associated with the black hole horizon from Eq. (25), we arrive at Thus, we note that the quantity 4πr 3 + /3 of Eq. (A9) is just an area of the black hole horizon. In proper units the Eq. (A9) may be written as S + = A/4G (see Appendix B). Thus, we conclude that the entropy of the 5D black hole obeys the area law, and the entropy of a black hole has no influence of the quintessence matter. Next, we turn our attention to the stability of the black holes by calculating the specific heat of the black hole solution (A2) and to discuss the effect of the quintessence matter. The Schwarzschild black hole and higher-dimensional Schwarzschild-Tangherlini always have negative heat capacity, indicating thermodynamic instability of these black holes [6]. Inserting Eqs. (A7) and (A8) in (29), we obtain It is well known that the thermodynamic stability of the system is related to the sign of the heat capacity (C). If the heat capacity is positive (C > 0), then the black hole is stable; when it is negative (C < 0), the black hole is said to be unstable. It is clear from Eq. (A10) that the heat capacity depends on the quintessence matter. When q → 0, one gets C = −4πr 3 + , which means 5D general relativity black holes are thermodynamically unstable [7,14]. Next, we analyze the effect of the quintessence matter on thermodynamical stability of a black hole. We plot specific heat (C) with radius r in Fig. 5 for different values of the parameter ω and q. It is seen that the heat capacity discontinuous at r = r C , for each ω, and for a given q. We observe that the heat capacity is C > 0 (C < 0) for r + < r C (r + > r C ). Thus, the 5D black hole is thermodynamically stable for r < r C , as the black hole has the positive heat capacity and unstable for r > r C (cf. Fig. 5). The black hole mass increases with increase in r + . Hence, the phase transition occurs from a lower mass black hole with negative heat capacity to a higher mass black hole with positive heat capacity.

Appendix B: Thermodynamics of d-dimensional black holes
The metric of d-dimensional spherically symmetric black hole surrounded by quintessence reads [21] where G (d) is the d-dimensional Newton constant, dΩ (d−2) is the line element on a unit (d − 2) sphere, and Γ (d−2) is the volume of the unit (d − 2) sphere. The black hole has a horizon at r = r + , which is solution of f (r + ) = 0 The black hole mass in terms of horizon r + reads The Hawking temperature is By applying the Wald entropy formula (25), one can deduce the general expression for the entropy as and it satisfy the area law.