Critical behavior of AdS Gauss–Bonnet massive black holes in the presence of external string cloud

Following previous study about AdS-Schwarzschild black holes minimally coupled to a cloud of strings in the context of massive gravity (Ghanaatian et al. in Effects of the external string cloud on the Van der Waals like behavior and efficiency of AdS-Schwarzschild black hole in massive gravity, arXiv:1906.00369 [hep-th]) and inspired by strong connection between Gauss–Bonnet Gravity and heterotic string theory, in this paper, we first take into account the Gauss–Bonnet term and we study thermodynamics and critical behavior of these black holes in the extended phase space. The effects of Gauss–Bonnet, massive, and string cloud parameters on the criticality of these black holes has been investigated. It can be seen that the Gauss–Bonnet and massive parameters have opposite effects on the criticality and phase transition of the solutions. We also observe that the increase in the value of the string cloud parameter above a critical value, eliminates the van der Waals like behavior of these solutions. Also, the Joule–Thomson effect is not observed. Then we examine thermal stability of these black holes in canonical ensemble by calculating the heat capacity. In addition, we explore critical behavior in extended phase space by employing heat capacity and consequently, we observe that the results are in agreement with the previous results from the usual method in Sect. 3.


Introduction
After the introduction of Einstein's theory of general relativity, at first a series of mathematical singularities appeared as solutions to Einstein's general equation of relativity. The first settled solution of general relativity as a black hole was found a e-mail: fhranji@gmail.com b e-mail: mehdi.sadeghi@abru.ac.ir c e-mail: m_ghanaatian@pnu.ac.ir d e-mail: forozani@pnu.ac.ir by Schwarzschild in 1916. For a long time, black holes were considered to be a mathematical peculiarity, until in 1973 Bekenstein introduced black holes as interesting thermodynamic systems that follow the laws of usual thermodynamics [2]. This progress made us one step closer to a better understanding of quantum gravity [3].
On the other hand, Einstein's equation generally described an expanding universe, and Einstein first added a cosmological constant to his equation to describe a static universe. But after discovering the expansion of the universe by Hubble in 1931, Einstein ignored the cosmological constant, until 1998, it was discovered that the expansion of the universe is accelerating, imposing a positive value for the cosmological constant. If the cosmological constant is positive, the associated negative pressure will derive an accelerated expansion of the universe, as observed from the Planck Collaboration [4]. So far, cosmological constant is considered as a constant parameter, but on the contrary, in black hole thermodynamics is regarded as thermodynamic pressure which can vary. From this perspective, an extended phase space appears with a new dimension added, and the negative cosmological constant determines a positive varying thermodynamic pressure in this new framework to review a few examples of work in the context of extended phase space, see [5][6][7][8].
Among the modified gravity theories, one of the interesting proposals is Lovelock theory, which is reduced to Einstein's theory in three and four dimensions, but in five dimensions and more, it includes higher curvature terms. The first additional term is Gauss-Bonnet (GB) gravity where is quadratic in curvature tensors and leads to field equations that are second-order unto metric derivatives and do not engage in ghosts [9,10]. It was also shown that the GB theory can be derived from the low-energy limit of heterotic string theory [11,12]. One of the promising theories in modern theoretical cosmology is the scalar-Einstein-Gauss-Bonnet gravity theory [13,14] which is motivated by string theory and shows how the string theory affects the primordial acceleration of the universe. Many other aspects of GB theory have been studied in literature [15][16][17][18][19][20][21][22][23].
One of the most consistent theory of gravity that modifies GR by endowing the graviton with a nonzero mass, is dRGT massive gravity [24]. This theory is ghost free and avoids discontinuities in the limit where the graviton mass goes to zero. Also, the theory of general relativity is modified by massive gravity at large distances, explains the accelerated expansion of the universe, without resorting to the concept of dark energy. In addition , the cosmological solutions of massive gravity and its expanded types, such as bimetric gravity [25], can reflect late-time acceleration compatibility with observations [26][27][28].
However, string theory predicts the existence of a graviton, but we do not have a successful quantum theory of gravity. Superstring theory is a theory that considers the particles and fundamental forces of nature as vibrations of tiny supersymmetric strings. The idea of taking fundamental particles as vibration modes of one-dimensional string objects plays an essential role in these theories. The fast accelerated expansion of the universe in the period of inflation can be attributed to the stretch of such cosmic strings that penetrated everywhere in our observed universe [29]. A cloud of strings, an aggregation of one dimensional objects in a certain geometrical frame, was proposed by Letelier [30]. The gravitational effects of matter in the form of sting cloud are studied [30][31][32].
On the other side, when a control parameter such as temperature is changed in a thermodynamic system, the system may change to a different macroscopic state that is more stable, which this mutation is called phase transition. In a phase transition, a thermodynamic potential such as free energy becomes non-analytic. Although the AdS black holes in the radiation can remain stable in the heat stability, but at a certain critical temperature, there is a phase transition called Hawking-Page phase transition [33]. This phase transition is mainly observed in Einstein's general relativity family such as Gauss-Bonnet gravity [34] and other theories like dilaton gravity [35], magnetic black brane [36], black Dp-branes and R-charged black holes with an IR Cutoff [37] and BTZ black hole [38]. There is also a resemblance between "small black hole/large black hole"(SBH/LBH) phase transition and the liquid-gas phase transition, that this critical behavior is called Van der Waals-like behavior. It can be seen that this behavior is commonly found in charged AdS black holes and in massive gravity theories, as well as the coupled theory with a cloud of strings [39][40][41][42][43][44].
Recently, the effects of a cloud of strings on the extended phase space of Einstein-Gauss-Bonnet AdS black hole is studied and the Van der Waals-like behavior in absence of the GB term is observed [45]. Also, in the previous work [1], we investigated the effects of the external string cloud on the Van der Waals like behavior of AdS-Schwarzschild black holes in massive gravity. With these explanation, in this paper, we study critical behavior of AdS Gauss-Bonnet massive black holes in the presence of external string cloud. One of the motivations of this study is to investigate the simultaneous effects of the massive, Gauss-Bonnet and cloud of strings terms on critical points and other is study the effect of each of these parameters on critical behavior by keeping the other two parameters in constant.
This paper is organized as follows. In Sect. 2, the solutions of AdS-GB black holes in dRGT massive gravity minimally coupled to a cloud of strings are introduced and the metric function and its diagrams in different modes are investigated. In Sect. 3, we examine the first law of thermodynamics of these black holes and employ extended phase space thermodynamics to explore critical points. Also, we study the behavior of system along the coexistence line by plotting isothermal curves in P − T diagrams. We then investigate the possibility of the Joule-Thomson effect in our model by drawing isenthalpic curves in T −P plan. In the following, the critical exponents are calculated in our model. The thermal stability of the solutions in canonical ensemble are studied in Sect. 4, and more, a search for critical behavior in extended phase space has been performed using heat capacity. Finally, in last section we will present our conclusions.

Black hole in GB-massive gravity minimally coupled to a cloud of strings
Let us start with an AdS GB-massive gravity in 5-dimensions in the presence of external string cloud. The action is as follows, where R is the scalar curvature, = −6 L 2 is cosmological constant with L as the cosmological constant scale, f is a fixed rank-2 symmetric tensor known as reference metric and m is the massive parameter. The last part called a Nambu-Goto action, in which N P is a positive quantity and is related to the tension of string, (λ 0 , λ 1 ) is a parametrization of the world sheet and χ is the determinant of the induced metric [46][47][48] L gb is the Gauss-Bonnet term of gravity with λ gb its dimensionless coupling. L gb is given by where R μν and R μνρσ are Ricci and Riemann tensors respectively.
In massive term, c i 's are constants and U i 's are symmetric polynomials of the eigenvalues of the 5 λ and the rectangular brackets denote traces.
A generalized version of f μν was proposed in [49,50] with The values of U i 's are calculated as below, The energy-momentum tensor for a cloud of strings is given by where ρ is the proper density of a string cloud and μ ν is the spacetime bivector in which ε ab is Levi-Civita tensor. Conservation of the energy-momentum tensor, ∇ ν T μν = 0 results in, The equation of motion is obtained by variation of the action with respect to the metric tensor g μν where G μν is the Einstein tensor, H μν and X μν are We consider the following metric ansatz for a five-dimensional static spherically symmetric metric Inserting this ansatz into the Eq. (8) yields two solutions for f (r ) in which b 4 is an integration constant and a is a real positive constant known as string cloud parameter. Since black hole should have an event horizon we pick up f − (r ). To calculate the radius of the event horizon, where m 0 is related to the total mass of the black hole with M = 3V 3 m 0 16π where V 3 = 4π 3 , is volume of the three dimensional unit sphere as plane or hyperbola.
To provide the estimated value for the dimensionless coupling coefficient associated with the GB-term, the welldefined constraint of the vacuum solution (b 4 = 0, a = 0 and m = 0) leads to 0 ≤ λ gb ≤ L 2 8 . Besides, the causality and positive requirement of the boundary energy density in holography requires that −7.72 ≤ λ gb ≤ 9.20, [45], then we can estimate L 8.579044 for maximum value of the cosmological constant scale in vacuum solution. However we see that if r → 0 then f (r ) approaches k so there is no r = 0 singularity in vacuum solution. Also in non-vacuum regime we observe that if r → 0 then f (r ) approaches k − m 0 2λ gb , it means adding the Gauss-Bonnet term to the action causes the causal singularity to be removed.
In order to investigate the effects of cloud string parameter and GB parameter on the metric function of gravitational theory we plot f (r ) − r diagrams in the following (Fig. 1). We observe that, due to the addition of the GB parameter, singularity is gone. By comparing the diagrams (a) and (b), we observe that in flat topology, the number of horizons decreases by decreasing the value of the GB parameter. It is also clear from the difference between the two diagrams (c) and (d) that the number of horizons increases with increasing GB parameter. Contrary to this, the comparison of the diagrams (a) and (c) shows that the increase of the string cloud parameter, in contrast to the GB parameter, has a decreasing role in the number of horizons. All in all, it can be said, the existence of a maximum of 3 roots in f (r ) − r diagrams is obvious. It should be noted that in diagrams (a) through (d), the change in the parameters of λ gb and a, except for the effect on the number of horizons, does not have a significant effect on the metric function, especially on asymptotical behavior, but for the negative value for λ gb , the diagrams (e) and (f) show that for radii less than a certain value, the metric function has no value and so is not defined. Also, for the negative values of λ gb , the asymptotical behavior is different from that of the positive values of λ gb .

Thermodynamics and critical behavior in the extended phase space
In this section, we first calculate the thermodynamic quantities of AdS black holes in GB-massive gravity minimally coupled to a cloud of strings in the extended phase space and we present the first law of thermodynamics and corresponding Smarr relation. Then we study the phase transition points and P − V criticality of these black holes. In extended phase space, the mass appears not only as an internal energy, but as an enthalpy of the thermodynamic system, as below The Hawking temperature of the black hole is obtained by applying the definition of surface gravity, is as follows The entropy of the black hole is given by using Wald's formula which clearly shows the correction of the area law with the GB-term for non-flat topology (k = 0). It is also obvious that the entropy of a string cloud and massive term does not affect the black hole entropy.
With these definitions, one can obtain the first law of thermodynamics in extended phase space in the following form where V , conjugating thermodynamical variable corresponding to pressure is thermodynamical volume of the black hole and A, B, C i 's represent for physical quantities conjugated to the parameters λ gb , a and c i 's respectively. By a dimensional argument one can present the Smarr relation as We observe that the massive term c 2 has scaling weight 0 and it is constant in the metric function, so, does not appear in Smarr relation. By putting the thermodynamic quantities presented above in Smarr relation or directly from Eq. (16), we obtain the equation describes the state of the thermodynamic system in the extended phase space, called the equation of state To investigate critical behavior, one can calculate critical values by using the inflection point properties, ( ∂ P ∂r 0 ) T = ( ∂ 2 P ∂r 0 2 ) T = 0, of the equation of state, as follow where r c , T c and P c are called horizon radius, temperature and pressure, respectively. The other thermodynamic quantity that can be calculated in the extended phase space for the considered black hole, is the Gibbs free energy given as In a phase transition, always a thermodynamic potential such as Gibbs free energy becomes non-analytic(discontinuous).
To further investigation of critical behavior and effects of GB, massive and cloud string parameters on the criticality, one can plot P −r 0 and G − T diagrams (Figs. 2, 3, 4, 5, 6, 7) and observe Van der Waals-like behavior clearly. Various diagrams are based on the increase of the GB, massive and cloud string parameters. From P − r 0 diagrams, it can be seen that by decreasing the value of the GB parameter and increasing massive parameter, Van der Waals-like behavior becomes more apparent. Also, the increase of the value of the cloud string parameter in comparison with the two abovementioned parameters leads to a diminution of critical behavior. It can be shown that there is no critical behavior and phase transition for values greater than a critical value of the cloud string parameter, a c , which satisfies the following condition that is, the above inequality holds for a ≤ a c . In addition, the appearance of discontinuity and swallow tail in G − T diagrams corroborates critical behavior and the occurrence of phase transition. It is obvious that the incremental effect of the massive parameter on criticality is in contrast to the GB and cloud string parameters. In Fig. 5, an significant increase in the value of GB parameter has reduced the criticality, which is more evident in flat topology. In Fig. 6, we see that increasing the massive parameter from zero to a common value causes the emergence of criticality and phase transition. Finally, in Fig. 7, it can be seen that the increase in the value of the string cloud parameter exceeds the critical value, eliminates the criticality.
So far we have only investigated the critical behavior of the system by fixing temperature. One can study the critical behavior of black hole temperature in P − r 0 plane and investigate the effect of GB, massive and a cloud of strings parameters on that (Figs. 8, 9, 10). We observe that when the temperature is above the critical value, a behavior similar to that of an ideal gas appears, which is referred to as the ideal gas phase transition. But for temperatures below the critical temperature, three branches are seen, representing small, medium, and large black holes. Except for the latter, which is unstable, the other two are stable and consistent with the Van der Waals liquid/gas phase transition. As we can see in Fig. 9, we find that for temperatures below the critical temperature, the pressure decreases with mass reduction being negative. Also, as seen in Fig. 10a, there exist a particular temperature for which we have ∂ P ∂r 0 = P = 0, similar to what we have seen in Van der Waals fluid before.
Another way to illustrate the phase transition is to use the P − T diagram for two different phases where the black hole phase transition is between the two so that both phases have the same Gibbs free energy. This phase transition is of the first order and occurs where two surfaces of Gibbs free energy intersect, known as coexistence line in P − T diagrams. At any point on this line, the following equations exist between the two phases mentioned, where the indices 1 and 2 correspond to the two different phases of the black hole. The temperature equilibrium of these two phases indicates the isothermal phase transition. We plot equation of pressure with respect to temperature for some values of λ gb , m and a parameters and we observe the effect of changing these parameters on P − T diagrams in Fig. 11. In this figure, p is the value of pressure per unit of critical pressure and τ is the value of temperature per unit of critical temperature. In addition, plotting T − P diagrams and examining the behavior of the thermodynamic system would be helpful. This method reveals a process known as Joule-Thomson expansion, which describes the change in system temperature relative to pressure at a constant enthalpy. That means we will have a isenthalpic process that can display heating and cooling phases. To identify the phase of the system, we need to denote the Joule-Thomson coefficient as μ J T = ( ∂ T ∂ P ) H . If μ J T > 0, it indicates the cooling process in which pressure decreases during the expansion and μ J T < 0 shows heating in which pressure increases. To plot the T − P curves and find the process type in our model, we use Eqs. (15) and (16) for different values of constant mass M (Fig. 12). As can be seen from the diagrams, our process has only one cool-    ing phase and never enters a heating phase as shown in [51] which follows a heating-cooling process in a Joule-Thomson expansion.
Finally, we study the behavior of physical quantities near the critical point. In order to calculate the critical exponents characterizing the behavior of physical quantities in the vicinity of the critical point, it is advisable to use rescaled quantities ν = υ υ c , τ = T T c and p = P P c where υ = 4 3 r 0 is specific volume. This simplification makes the equation of state (26) as follows where is called as law of corresponding state. We can now search for the thermodynamical behavior of the system near the critical points by redefining parameters t, ω That is, the ν, τ and p parameters are expanded around one, so the law of corresponding state would be approximated as To characterize the critical behavior near the critical point, one can introduce the critical exponents as [52] As is clear from the above definitions, the exponents α, β, γ , and δ describe the behavior of specific heat with fixed volume, the order parameter η, the isothermal compressibility coefficient κ T , and the critical isotherm, respectively. The subscripts l and s represent the large black hole and the small black hole, respectively, in the phase transition process. The entropy S does not depend on the Hawking temperature T , so the specific heat at constant volume C υ is equal to zero, consequently the corresponding critical exponent vanishes (α = 0). To estimate the second exponent β, one can evaluate υ l and υ s to calculate the order parameter. During the phase transition the pressure of the black hole keeps unchanged. It results that the large black hole pressure equals the small black hole pressure, p l = p s for which On the other hand, from the Maxwell's equal area law, one can further obtain With two above equations, one can get So the order parameter can be derived as where this leads to the conclusion that β = 1 2 . The isothermal compressibility can be estimated as follows From this one can conclude that γ = 1. The critical isotherm is an isotherm process at critical temperature T = T c or t = 0. Then we can conclude that p − 1 = ω 3 that leads to δ = 3. We see that, the values of critical exponents are independent of GB, massive and cloud string parameters. The critical exponents in our model are the same as those mentioned in other articles [52][53][54], and all the models reviewed have the same scaling laws.

Heat capacity in canonical ensemble vs in the extended phase space
In this section, we first study the stability of the black hole by using heat capacity in canonical ensemble. Then we investigate criticality in extended phase space by employing heat capacity. When the black hole heat capacity is positive, we say that the black hole is thermally stable. On the other side, an unstable black hole may turn into a stable state, which this transition is called phase transition. If we calculate the heat capacity of this black hole, using it, we can determine the type of phase transition. In this way, the presence of roots and divergence points for the heat capacity will represents the type one and the type two phase transition. When the heat capacity is in the form of a fraction, in order to obtain the roots, we set the numerator to zero and for calculation the divergence points, we set the denominator to zero. The heat capacity is calculated as follows It is carefully observed in these equations, since the GB parameter is coupled to the curvature factor, in flat topology, the behavior of heat capacity is independent of GB gravity.
To check the black hole's thermal stability, T −r 0 and C −r 0 diagrams are ploted (see Figs. 13,14,15,16,17,18). For a region of the event horizon, which temperature is negative, the solution is non-physical, and we have removed it from T − r 0 diagrams. It is observed that, for spherical (k = 1) and flat (k = 0) topologies, temperature have one root and for hyperbolic (k = −1) have two roots. If the radius of the black hole event horizon is called the size of the black hole, there will be a black hole with minimum size and non- zero entropy at zero temperature. As one can see, increasing the value of the massive parameter leads to the formation of extrema on the T − r 0 diagrams, indicating type two phase transition (Fig. 13). It is also seen in Fig. 14 that the decrease of the GB parameter has a direct effect on the formation of type two phase transition. Finally, the increase of the string cloud parameter above the critical value results in the loss of the type two phase transition (Fig. 15). It is observed clearly in Fig. 16 that, for spherical topology, with the decrease of the GB parameter, the black hole's instability domain becomes larger. In Fig. 17, the C −r 0 diagrams are ploted in three modes in terms of the massive parameter increase. For small values of massive parameter (Fig. 17a), there is no divergence points, ie, the massive term has a direct effect on the presence of type two phase transition. When the massive parameter increases sufficiently enough (Fig. 17b,  c), the heat capacity has divergencies, indicating type two phase transition, ie, the black hole in these points changes the phase between stable and unstable states. As the massive parameter increases, the black hole's instability domain extends. Finally, in Fig. 18, we observe that with the increase of the cloud string parameter, the black hole stability domain grows, as long as the string cloud parameter exceeds the critical value, the solutions become completely stable.
As mentioned earlier, to calculate the divergence points, we must denote the denominator of the fraction of the heat capacity equal to zero. In extended phase space, the critical values in which phase transition takes place are the same as divergencies in heat capacity. Therefore, as stated in reference [55], by substituting (1) in the denominator of the heat capacity and solving it with respect to pressure, the following relation is obtained  Fig. 19, we can say that when λ gb → 0, the pressure has maximum and the phase transition takes place. Therefore, by using this method, one can calculate critical horizon radius and pressure without referring to the complex formulas given in the usual way.

Conclusion
In this paper, we have studied GB-massive black holes in the presence of external string cloud. It is observed that the solutions of the gravitational theory modified by GB term are protected from the causal singularity. It is worth noting that in all of the computed quantities throughout the paper, apart from the metric function, the GB parameter is coupled to the curvature factor, which means that in the case of flat topology (k = 0), GB gravity has no effect. We observe that the first law of the black hole thermodynamics is modified because the black hole entropy gains a contribution from the GB term. We also investigated the criticality and phase transition of the solution in the extended phase space by employing the equation of state. It is important to note that the massive parameter plays a crucial role in creating critical behavior. Also, when λ gb → 0, the criticality becomes more obvious.
It was shown that, the effects of GB and massive parameters on criticality are opposite of each other. Also, the effect of the cloud string parameter on the criticality, like the effect of the GB parameter, is a decreasing effect. In other words, the cloud of strings coupling with the desired gravitational theory has a positive effect on the critical behavior of the theory's solutions, provided that the value of the string cloud parameter does not exceed its critical value. That is, when the cloud string parameter becomes more than critical value, the solutions are completely stable and the phase transition disappears.
By studying critical behaviour of the black hole temperature, we conclude that when temperature is above the critical value, Van der Waals-like phase transition would be disappeared and an ideal gas-like behaviour recovers. We also plotted the coexistence line in P − T diagrams in which two phases are in equilibrium and Gibbs free energy and Hawking temperature keep unchanged during transition.
We investigated T − P diagrams to find out the sign of the Joule-Thomson coefficient. It is shown that the system has only one cooling phase and never enter to a heating phase, that is, the Joule-Thomson effect does not happen. Investigations on the critical exponents showed that the GB, massive and cloud string parameters do not affect them.
Then, we examined thermal stability of these black holes by calculating heat capacity. As expected, it was observed that, the effects of the above three parameters on stability was the opposite of their effects on the criticality. Finally, using the method mentioned in reference [55], showed that this method is much simpler than the usual method for investigating critical behavior.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The data are available in the article and there is no external data.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .