Phase structure of charged AdS black holes surrounded by exotic fluid with modified Chaplygin equation of state

By considering the concept of the modified Chaplygin gas (MCG) as a single fluid model unifying dark energy and dark matter, we construct a static, spherically charged black hole (BH) solution in the framework of General Relativity. The $P-V$ criticality of the charged anti-de Sitter (AdS) BH with a surrounding MCG is explored in the context of the extended phase space, where the negative cosmological constant operates as a thermodynamical pressure. This critical behavior shows that the small/large BH phase transition is analogous to the van der Waals liquid/gas phase transition. Accordingly, along the $P-V$ phase spaces, we derive the BH equations of state and then numerically evaluate the corresponding critical quantities. Similarly, critical exponents are identified, along with outcomes demonstrating the scaling behavior of thermodynamic quantities near criticality into a universal class. The use of \emph{geometrothermodynamic} (GT) tools finally offers a new perspective on discovering the critical phase transition point. At this stage, we apply a class of GT tools, such as Weinhold, Ruppeiner, HPEM, and Quevedo classes I and II. The findings are therefore non-trivial, as each GT class metric captures at least either the physical limitation point or the phase transition critical point. Overall, this paper provides a detailed study of the critical behavior of the charged AdS BH with surrounding MCG.


I. INTRODUCTION
General relativity (GR) provides the best available description of gravity so far [1].Among its most enthralling predictions, gravitational waves (GWs) and black holes (BHs) deserve special mention, as they represent the ultimate confirmation of Einstein's theory.On the one hand, the first direct detection of GWs dates back about a decade ago with the emission of signals from a binary BH merger and the subsequent ringdown of the single resulting BH [2].On the other hand, Xrays coming from superheated material swirling around a dark object were recognized as distinctive evidence of a central BH (Cygnus X-1, located within the Milky Way) in the early sixties, although we had to wait another half-century to capture the first ever picture of a similar spacetime oddity -M87* [3].Since then, GWs and BHs have been extensively addressed as a new way to probe the cosmos at a fundamental level.In particular, it is commonly accepted that BH physics could offer valuable insight into the unification of GR, quantum theory, and statistical mechanics [4][5][6][7], opening a novel route to quantum gravity.
A class of solutions of Einstein's equations that have been attracting growing interest in recent years are antide Sitter (AdS) BHs.The observation that asymptotically AdS BHs can be modeled in the language of dual thermal field theory has motivated a fluid-like description of the underlying microphysics.Interestingly enough, such studies have revealed that the equilibrium thermodynamics of BHs can be investigated by considering the geometric properties of their event horizons and other relevant spacetime features [8][9][10][11][12][13][14][15][16][17][18][19].In the ensuing geometrothermodynamic (GT) picture, the scalar curvature of the BH metric represents the thermodynamic (microstructure) interaction, with positive curvature indicating prevailing repulsion and vice versa.Moreover, singularities represent the breakdown of the classical theory because they correspond to an infinite gravitational interaction which cannot be treated properly in GR.The microscopic behavior of BHs has been analyzed in the geometrothermodynamic framework for a wide class of systems [9,13,17,18,20] and in various entropic scenarios [21][22][23][24][25][26].
Further suggestive features of AdS BHs are phase transitions and critical phenomena, which were first demonstrated in the phase space of Schwarzschild-AdS BH non-rotating and uncharged [7].This seminal discovery has opened up a new line of research in the field of BH thermodynamics.In this context, non-trivial results have been achieved for spinning branes [27,28] and charged Reissner-Nordstr öm (RN) BHs [29,30], whose first-order phase transitions display a critical behavior analogous to a van der Waals (vdW)-like (i.e.liquidgas) change of phase.Recently, the correspondence between BHs and condensed matter systems has been explored further by incorporating the variation of the cosmological constant Λ in the first law of BH thermodynamics [31][32][33][34][35][36][37] (see also [38] for more discussion), which allows maintaining consistency with the Smarr relation [32].In this scenario, the BH mass is identified with enthalpy rather than with internal energy.Furthermore, since Λ corresponds to pressure, it is natural to consider the thermodynamic volume of BH as its conjugate variable.
Taking the above arguments seriously, various thermodynamic quantities (such as adiabatic compressibility, specific heat at constant pressure, etc.) have been computed by using standard thermodynamic machinery [33][34][35].Remarkably, the critical behaviour of AdS BHs has been reconsidered in an extended phase space, including pressure and volume as thermodynamic variables [34].As a result of this approach, the P = P(V, T) equation of state has been analyzed for a rotating charged AdS BH, emphasizing analogies with the vdW P − V diagram.The recent study in [38] has gained further progress in the identification of charged BH firstorder transitions with the standard liquid-gas phase transitions by examining the behavior of the Gibbs free energy of a RN-AdS BH in the canonical ensemble (i.e.fixed charge).
Predicting critical points in an BH system is crucial for phase transition analysis.In this way, GT provides an alternative description using a thermodynamic metric, governed by a thermodynamic potential and its derivative with respect to extensive parameters.This metric is implemented in the thermodynamic equilibrium phase space.In this regard, the Weinhold geometry is considered the first attempt to analyze the critical phase transition of a thermodynamic system [39,40].An al-ternative approach, applied to the exploration of critical thermodynamic characteristics, is the Ruppeiner geometry [41,42].As such, these two approaches are known to be closely linked via the inverse of temperature [43].In fact, the Legendre transformation is vital in the context of GT due to its invariance on the thermodynamic potential [44].This invariance enables the use of the Legendre-invariant metric in thermodynamic phase space [45].By the way, the Quevedo metric is a better example of the Legendre-invariant metric [46,47].Furthermore, Hendi et al. [48][49][50][51][52] invented another representation base of a thermodynamic system, a metric that provides a perfect correlation between the curvature scalar and the phase transition point of the heat capacity.Similarly, Mansoori et al. [53][54][55] currently developing another thermodynamical metric.However, there were a few limitations to these thermodynamic metrics.In general, the main motivation behind the use of GT is that it offers an independent view of thermodynamic systems, aiding in the inspection of the bound points, patterns, and stability of phase transitions.It provides microscale behavior and the Ricci scalar sign, indicating repulsive or attractive interactions, i.e., negative or positive, respectively, along the transition curve, while R = 0 indicates the absence of interaction [41].
In the framework of modern cosmology, one of the most conspicuous problems is the explanation of the matter/energy content of the Universe.It is a fact that the amount of invisible dark sectors is about 95% of the total density [56].This dominance has stimulated great efforts to uncover the origin of these mysterious entities.A challenging possibility is that dark components may surround (or even be created inside) BHs.Along this direction, special focus has been devoted to analyzing static spherically symmetric BHs surrounded by quintessence matter (see [57][58][59][60][61][62] and references therein).Other solutions are in [63][64][65][66][67][68][69].Concretely, attempting to describe the acceleration of the universe attributable to a kind of exotic negative-pressure fluid, widely known as dark energy (DE), is a challenge in the fields of astrophysics and theoretical physics.Based on this challenge, various dark energy model candidates provide phenomenological and theoretical predictions to explain the acceleration of the universe [70][71][72][73][74][75].In addition to the familiar types of dark-energy model, there are other types of cosmological dark-energy model deemed interesting, such as the cosmological constant and the universe full of exotic fluids.Keep in mind that the acceleration phase in the late universe is broadly attributed to the cosmological constant Λ, or a constant energy density ρ Λ = 7.0210 −24 g/m 3 , dubbed dark energy.Indeed, the problem puts a frame on looking for a unification of the dark sector of the universe and considers it as a single component that behaves as both dark energy and dark matter.
Recently, new models that mix dark matter and dark energy have been proposed as candidates for the dark components.Among the suggested unified dark fluid models, the Chaplygin gas [74] and related generalizations [76,77] have been largely adopted to explain the accelerated expansion of the Universe [70,[78][79][80].Further applications of the Chaplygin dark fluid (CDF) appear in relation to the Hubble tension [81] and the growth of cosmological perturbations [82], respectively.In [83] the analytical solution and the related thermodynamic quantities have been addressed for a charged static spherically-symmetric BH is surrounded by CDF in the Lovelock gravity theory.Such a model has been later extended to the modified Chaplygin gas (MCG) to study the stability of the MCG-surrounded BHs in Einstein-Gauss-Bonnet [84] and Lovelock [85] gravity.On more thermodynamic shores, the phase transitions and critical behavior of the static spherically-symmetric AdS BHs surrounded by CDF in GR have been the subject of investigation in [86], while a preliminary geometro-thermodynamic analysis has been conducted in [87].All these studies indicate that the (generalized) CDF could not only be a theoretical model, but a naturally existing fluid deserving of further consideration.
Starting from the above premises, in this work, we aim to study the phase structures and transitions of a charged AdS BH surrounded by MCG.Toward this end, the analysis is structured as follows: the next Sec.II is devoted to discussing the charged BH solution surrounded by MCG.In Sec.III we examine the general features of BH chemistry for the obtained solution, while Sec.IV contains the analysis of thermodynamic stability.The study of critical behavior and phase transitions on the basis of the geometrothermodynamic tools are the main objectives of Sec.V and VI, respectively.Conclusions and perspectives are finally summarized in Sec.VII.

II. CHARGED BLACK HOLE SOLUTION SURROUNDED BY MODIFIED CHAPLYGIN GAS
This section presents an examination of the MCG structure within the framework of GR theory.Thus, it is carefully noted that pure Chaplygin gas is said to have an exotic equation of a negative pressure state In addition, the extended version of this gas is the generalized Chaplygin gas (GCG), which is typified by the generalized equation of state with 0 ≤ β ≤ 1 being a universal constant.In practical terms, MCG is a major extension of Chaplygin gas, which is an exotic fluid considered perfect, obeying the calorific EoS p = Aρ − B/ρ β where A ≥ 0, β ≥ 0 and the pair A and B are considered constant through an adiabatic process [88,89].Parametrically speaking, conditions β = 1 and A = 0 lead to a characteristic of the lifting forces on a plane wing in the aerodynamic process.On the other hand, constraining β > 0 and A = 0 define the GCG.On top of that, there is plenty of evidence that the MCG can be rebuilt using k-quintessence (kinetic quintessence) [90,91] and fquintessence (fermionic quintessence) [92].Inspired by the exciting results of the MCG-specific BH chemistry method by assuming the negative cosmological constant as a thermodynamical pressure, revealing the critical behavior from the MCG background is needed [93,94].
In this regard, Ubbanah used a specific BH known as Chaplygin's BH and assumed the negative cosmological constant as the thermodynamic pressure.As a consequence, the work deserves further contributions, and the system is considered to be a heat engine [93].Similarly, modified cosmic Chaplygin gas is examined in the context of the AdS feature by considering BH's chemistry formula.This work highlights several thermodynamic aspects, including the heat engine and Joule-Thomson expansion [94].
In what follows, we consider the system of a charged source with a MCG structure in the context of GR described by the following action: where R is the Ricci scalar, g = det(g µν ) is the determinant of the metric tensor g µν , ℓ is the AdS length, I M is the matter contribution arising from the MCG background, and is the field strength of the electromagnetic field with A µ is the gauge potential.Moreover, κ = 8πG, where G is the Newtonian gravitational constant.Henceforth, we consider G = c = 1.Varying the action (3) leads to the following field equations: where G µν is the Einstein tensor, T CMG µν is the energymomentum tensor for MCG, and T EM µν is the energymomentum tensor for the electromagnetic field which is explicitly given by We consider a static, spherically symmetric, fourdimensional space-time given by in this metric, f (r) is the metric function that depends on r and dΩ 2 = dθ 2 + sin 2 θdϕ 2 .The way to reveal the structure of the electromagnetic field, in particular an electrically charged BH solution, is to assume a radial electric field.Thus, in this choice, the gauge potential is expressed as Exploiting Eqs. ( 6) and ( 7) provides a second-order differential equation given by where the prime and double prime are the first and second derivatives with respect to r.To proceed to solve the differential equation yields a compact solution in the form of where q is an integration constant acting like an electric charge.Taking into account the structure of the gauge field (10) in conjunction with the space-time metric (7), the field strength component appears as follows: We assign the MCG the equation of state p = Aρ − Bρ −β [88,89], where A, B are positive parameters and the parameter β runs in the range 0 ≤ β ≤ 1.In a 4D spherically symmetric spacetime scenario, the related energy-momentum tensor components of the MCG are expressed in the following way: It should be noted that the expression for the energymomentum tensor was first given by Kiselev when studying the static spherically symmetric quintessence surrounding a BH.Due to the fact that spacetime is considered static and spherically symmetric, the r − r component of the energy-momentum tensor should be equal to the t − t component, so On the basis of the isotropic average over the angles, the following where one can obtain It is noteworthy that the given energy density and parameter set require an explicit and constraining statement.To this purpose, taking into account Eqs. ( 13) and ( 15) simplify the process and yield the following expressions: Furthermore, by exploiting equations ( 13) and (15), the parameters α i and β i are constrained to be given by Thus, it is now simple to define the angular components of the energy-momentum tensor devoted to the MCG structure by By considering the spacetime metric (7) together with the field equations (4), the Einstein tensor components can be given by whereby combining the previous set with that of Eqs. ( 13)-( 20) produces a pair of differential equations such as: To properly look for the expression of the energy density, consider the set previously described above by using equations ( 13)- (20) and, in a way, by imposing the conservation condition on the energy-momentum tensor.This consideration step produces an interesting result, as , (25) where an integration constant, γ > 0, is used.It is important to keep in mind that the energy density is restricted to ∼ B 1+A 1 1+β at the asymptotic limit for r.As a result, MCG behaves as a cosmological constant far away from the BH, and as it approaches the BH, it gravitationally grows more densely.
Substituting Eq. ( 25) into Eq.( 23), one can obtain the analytical solution for f (r) in the following compact form where M is the physical parameter standing for the mass of the BH.Whereas, the integration constant q is mainly related to the physical parameter charge Q over a twosphere with an infinite radius according to the conservation law as where ϵ tr = n t σ r = 1 for a static spherically spacetime with n µ is the unit normal and σ ν is the unit normal to a unit two-sphere with an infinite radius.This proves that the physical electric charge Q in our charged BH solution is the same as in Maxwell's equation.Also, the electric potential evaluated at infinity by an observer with relation to the event horizon r + can be given as where χ = C ∂ t is the null generator of the event horizon, while C acts as a fixed constant.Here, the hypergeometric function 2 F 1 [α, ν; λ; ξ] representing the regular solution of the hypergeometric differential equation, is defined for |ξ| < 1 by a power series of the form with (n) k is the (rising) Pochhammer symbol [95].Furthermore, the parameters set (α, ν, λ, ξ) are given, respectively, by .
To explore the asymptotic behavior of the metric function f (r), the limit r → ∞ is taken into account, giving which implies that in the asymptotic limit, the behavior of the solution is controlled by means of the AdS length ℓ (the cosmological constant Λ) and the MCG parameters background.As a result, one can constrain 3 . It may be further pointed out that the metric function can be close to RN-AdS BHs once conditions β = 0, γ = 0, and A > 0 are set after taking the limit B → 0.
To provide an examination regarding the behavior of the metric function for our BH solution, Fig. 1 depicts this behavior.For a certain choice in the parameter space, the metric function is being schematically considered.In particular, for all sets of the parameter space, the root of f (r) is divided into two classes: the smallest root corresponds to two BH horizons, and the largest root is associated with a cosmological horizon.It appears that the two BH horizons, whether event or inner, are closely located at the same horizon radius for the variation of the parameter set (A, B, β).On the other hand, varying the charge parameter at the level of the metric function generates the apparition of multiple event horizons.A closer observation shows that the behavior of the metric function is disproportional with respect to the pairs (A, Q) and (B, β).In what follows, our BH solution will be examined in part by studying its singular behavior and, secondly, by analyzing its violation or satisfaction with respect to the constraints of the Energy Conditions (EC).
To properly verify the singularity and uniqueness of our BH solution, we need to provide an analysis based on scalar invariants.These are the Ricci scalar, the Ricci square, and the Kretschmann scalar.The Ricci scalar for the corresponding metric is given by The Ricci squared is given by Finally, the Kretschmann scalar is found to be where A closer look at the expressions (31), ( 32) and (33) shows that the BH solution represented by this metric is singular for any permissible value of the parameters A and β.Practically speaking, the existence of the singularity results from the mass and charge terms in the BH metric.By imposing the constraints β < 0 and A > 0, the singularity will consequently disappear due to the prod-uct constant γ.Nevertheless, to rule out the singularity made from the mass and the charge terms, it may be useful to describe a procedure with a non-linear charge distribution function similar to Ref. [96].Throughout this work, we will not be thinking about such a situation and will stick to the metric function (26) for the rest of the analysis.A few remarks concerning the mentioned scalars show that the Ricci scalar is not a function of the BH charge Q.In contrast, Ricci squared and Kretschmann scalars are functions of the BH charge Q, so any variation in the BH charge may induce meaningful variations in these same scalars.Therefore, the considered scalars demonstrate that our BH solution is unique, and the AdS background, together with the MCG structure, changes the BH spacetime significantly.
In what follows, we take care to study the energy conditions for our BH solution [97][98][99][100].The elements of the stress-energy tensor T µν governed by Einstein equations (4) for charged AdS BHs surrounding the CMG are as follows: • The weak energy condition (WEC) requires that T µν t µ t ν ⩾ 0 everywhere, for any time vector t µ , which is equivalent to and so The WEC is satisfied since A, B and γ are positive parameters and 0 ⩽ β ⩽ 1. Fig. 2 clearly shows that the energy density ρ along with ρ + P θ is positive for a span of horizon radii.
• The zero energy condition (NEC) stipulates that T µν t µ t ν ⩾ 0 in the overall spacetime for any null vector t µ .The NEC predicts ρ + P r ⩾ 0 which is identically zero, and ρ + P θ ⩾ 0 which is satisfied for equation (41) whenever |A| > 1.
• The strong energy condition (SEC) asserts that T µν t µ t ν ⩾ 1/2 T µν t ν t ν globally, for any time vector t µ which assumes that Consequently, it is obvious that the MCG does not satisfy or violate the SEC, as shown in Fig. 2.This is similar to the situation with the quintessence of dark energy.
Indeed, a violation of the SEC is interpreted as a violation of the attractive behavior of gravity, as evidenced by the dark energy that accelerates the expansion of the universe in cosmological studies, along with the matter content of the background of a regular BH, whose singularity has been superseded by a Sitter core.

III. BLACK HOLE CHEMISTRY WITH MODIFIED CHAPLYGIN GAS
Investigating the thermodynamic quantities of our BH solution is the main task of this section.Therefore, at the radius of the r + horizon of the metric solution (26), the mass of the BH is expressed in such a way as Subsequently, in order to determine the Hawking temperature, it is first necessary to consider the surface gravity [101], which is provided by with ξ µ = ∂/∂t is a killing vector.So, the formula T + = κ/2π is the Hawking temperature, expressed in terms of the BH system parameters as It is worth mentioning that all the parameters set in the BH solution effectively contribute to the behavior of the Hawking temperature.
To describe the behavior of the Hawking temperature, Fig. 3 depicts graphically the Hawking temperature against the horizon radius.It is worth noting that all the parameter variations in the BH system contribute to the effect on the Hawking temperature.Remarkably, the Hawking temperature rises to a maximum at T max + for a specific horizon radius r + .In particular, the maximum related to the Hawking temperature spanned the interval 1.04269 < T max + < 1.06535 for all the considered parameter spaces, and that maximum grows significantly with the charge parameter variation, which is To carry out a nice generation of the corresponding thermodynamic quantities, it is useful to apply the first law of BH thermodynamics.In particular, the first law can be given as [102] where µ i are the chemical potentials corresponding to the conserved charges N i .Holding parameters constant, with the exception of entropy, one can find which is similar compared to multiple background studies [5,6,103,104].
On the other hand, by considering the extended phase space, certain critical processes as a thermodynamic aspect are achieved.At this point, the P − V criticality, which will be a main feature of this work, is treated once the negative cosmological constant behaves like pressure, also known as "black hole chemistry".So, one has [32,105] In light of this objective, the first modified law of thermodynamics is formulated as follows [33,106]: where mass resembles enthalpy.Whereas, to construct the thermodynamic phase space framework, in par- ticular, the parameter function f (r + , M, P, Q) must always vanish under any transformation of the parameter.Further remarks on this subject make similar arguments to consider the constraints f (r + , M, P, Q) = 0 and δ f (r + , M, P, Q) = 0 on the evolution along the space of parameters.Nevertheless, an alternative way is considered: taking the mass parameter, M, as a function of the parameters M(r + , ℓ, Q) too.
The thermodynamic parameters are S, P, and Q.It is then convenient to redefine M = M(S, P, Q) for the possibility of explicitly obtaining dM = ∂M ∂S This is similar to such a differential 1-form in the space of parameters.Accordingly, all components are nothing more than thermodynamic quantities expressed in the context of extended phase space.So, one has In an alternative way, the same results can be expressed in accordance with the variation along the space of the parameters of the condition described by f (r + , M, P, Q), reshape another term for dM giving as follows: which must be in conformity with equation (49).Eq. ( 55) embraces the presence of temperature, which is geometrically defined as and which is a well-known finding, providing It should be pointed out that this expression can also be derived using Wald's formalism; basically, δS = δ ∂L ∂R as long as d f = 0 is satisfied.
Furthermore, the thermodynamic volume and the conjugate potential are defined by the following for-mula: In the case of the BH system, the enthalpy is defined by the total mass of the system.Consequently, in terms of the parameters BH system and in the context of the extended phase space, the thermodynamic volume and the electric potential can thus be formulated as while the Hawking temperature is given by According to Euler's theorem [32,107], with M(S, P, Q), the Smarr formula can be constructed for the charged source in the framework of GR as in which another combination arises between thermodynamic quantities.These amounts in the classical limit (Q → 0) are consistent with the corresponding amounts for the Schwarzchild AdS BH [108].

IV. BLACK HOLE STABILITY
Within the canonical ensemble, heat capacity is a further thermodynamical quantity providing information on the thermal state of the BHs.indeed, the heat capacity involves three specific and fascinating pieces of information.Firstly, the discontinuous behavior of this quantity means the existence of possible thermal phase transitions that the system may undergo.Second, another feature deals with the sign of heat capacity.In essence, the sign shows whether the system is thermally stable or not.In other words, positivity generates thermal stability, whereas the opposite indicates instability [109,110].Third, the corresponding possible set of roots is useful given leads to indicate the sign change, which may show stable/unstable states or bound points.As a consequence of these points, this section and the next are devoted to computing the relevant heat capacity in the extended phase space and inspecting the relation between this quantity and the GT tools.One can compute the heat capacity according to the following form [101] : Likewise, the function is expressed in terms of the parameter space as where Within the realm of BH physics, it is claimed that the associated roots of the heat capacity (C P = T = 0) indicate a one-dimensional line between physical (T > 0) and non-physical (T < 0) BHs, known as the physical limit.In this region, the system exhibits a sign change in heat capacity.On the other hand, the set of divergence points of the heat capacity represents the phase transition critical points of BHs [111].Thus, both the phase transition critical and limitation points of the BHs are explicitly computed with consideration of the following constraints: = 0 phase transition critical points So as to find the physical limitation and phase transition critical points, we consider Eq. ( 43) and solve the following equations for the entropy or in terms of the horizon radius: Attempting to solve these equations appears to be analytically difficult.This is why we apply the numerical approach to find physical limit points and phase transi-tion points for a given parameter space.The presentation of the heat capacity root point and singularity point is illustrated for two given parameter spaces in Tabs.I-I I.  To get an appropriate description of the heat capacity behavior at constant pressure, Fig. 4 presents the variation of C P against the horizon radiusr + .The appearance of divergent points and intersecting points both correspond to a change of sign, which indicates a phase transition, either a second-order phase transition or a firstorder one, respectively.Graphically, the analysis of the behavior of the corresponding heat capacity is present in Fig. 4. Thus, the upper panel concerns a fixed-valued parameter space that generates the behavior of heat capacity as a function of horizon radius.It should be noted that for pressure less than the critical pressure (P < P c ), especially P = 0.0306, there are two roots for the heat capacity (Tab.I).In this case, these two roots are in fact two critical horizon radii, i.e., r ⋆ 1 and r ⋆ 2 (the second root exists only for the upper panel with a certain parameter space), at which, for r + < r ⋆ 1 , the heat capacity of BHs is negative.So, the BHs with a small horizon radius are unstable.In addition, for r ⋆ 2 > r + , the heat capacity is positive which claims that the BHs are stable.It can be clearly observed that there is an important region due to r ⋆ 1 < r + < r div 1 , where r div 1 is a divergence point.The heat capacity of BHs in this region is positive.In other words, BHs have thermal stability when their horizon radius lies within the interval r ⋆ 1 < r + < r div 1 .Similarly, for the region r div 1 < r + < r ⋆ 2 , the heat capacity is negative, which means that our BHs are unstable.On the other hand, for (P < P c ) in the bottom panel with fixed values in parameter space, there are two roots and two critical phase transition points generating four intervals along the horizon radius.Clearly, BHs are thermally stable for regions r ⋆ 1 < r + < r div 1 and r + > r div 2 , while they are thermally unstable for regions r + < r ⋆ 1 and r div 1 < r + < r div 2 .Furthermore, as the pressure closes to critical pressure (P = P c ), the number of physical limitation points and divergent points reduces, and the heat capacity is positive, which means that the BHs are thermally stable and have a one-phase transition.This behavior is similar for the two sets of parameter spaces regarding the heat capacity behavior.For (P > P c ), the critical behavior is no longer present due to the absence of such a discontinuity provided by the divergent point.This situation remains the BHs thermally stable.By contrast, physical limitations are always present.Roughly speaking, our BH solution clearly remains in a thermally stable state.TABLE II: The phase transition critical points with A = 0.1.

V. P − V CRITICALITY
In this section, special care is taken to study the critical P − V for charged AdS BHs surrounded by an MCG structure.As such, by exploiting the expression for the Hawking temperature with (48), we can obtain the corresponding equation of state as follows: Of course, the parameters β, A, B, γ and Q potentially affect this equation.As intended, we can define the specific volume v = 2r + , wherewith the pressure is cast in the standard form P = 1 v T + O(v).Moreover, since the thermodynamic volume V ∝ r 3 + the critical point can be inspected considering the following constraints: or alternatively, The previous constraints set the triply critical point (T c , P c , r c ) with the unknown critical horizon radius r c which needs numerical solving.one has and In Tab.I I I, the numerical sets of critical points are shown as a function of various BH system parameter  values.Thus, the universal constant ratio P c r c T c clearly varies with the variation of the BH system parameter, being proportional to the set of parameters (Q, B, and γ) and disproportionate to the β parameter.
To highlight the key properties of P-V criticality, the T-r + and P-r + diagrams are useful at this stage, Fig. 5 shows the isobaric curve on the T − r + diagram, which is classified according to multiple pressure values.Indeed, the appropriate case P < P c implies two extreme points along the sides of three branches, namely the small BH branch, the intermediate BH branch, and the large BH branch.The small and large BH branches are characterized by a positive slope, meaning that the heat capacity is positive, and the system is thermally stable.While the intermediate BH branch stands on a negative slope, the BH is thermally unstable due to the negativity of the heat capacity.On the other hand, the case P > P c indicates the absence of any extremal point which proves that the system displays one stable BH branch in this situation.Alternatively, Fig. 6 displays, throughout a spectrum of temperature T, the P-r + criticality behavior via the isotherm curve P-r + .It is obvious all too well that the diagram shown is very similar to the VdW liquid-gas system (6).In this way, the situation T < T c predicts the existence of a small-large BH phase transition underlying the liquid-gas phase transi-tion of the VdW.As the temperature reaches the critical temperature (T = T c ), the number of phase transitions is limited to one, for which the first-order phase transition merges with the second-order phase transition like in a real gas system.On the other hand, once (T > T c ), the system displays one-phase behavior, and it is in fact an ideal with an empty set of phase transitions.
The present step applies a different thermodynamic potential to achieve perfect disclosure of the phase transition of the BH system.In particular, the Gibbs free energy is a thermodynamic quantity calculated from the Euclidean action with an appropriate limit term.A fruitful feature resulting from the sign of the Gibbs free en-ergy enables a global stability analysis.In the extended phase space, the thermodynamic potential is now the Gibbs free energy G = M − T S = H − T S. In practice, it should be noted that any discontinuous behavior in the first-or second-order derivatives of the Gibbs energy leads to a first-or second-order phase transition in the system.The Gibbs free energy is therefore given by To provide an appropriate analysis of the behavior of the Gibbs free energy, either as a function of temperature or in terms of the horizon radius, Figs.7 and 8 easily show this behavior.To start with, it is well-known to mention that Figs. 5 and 7 represent the same BH branches, where the same spectrum of values corresponds to the same pressures.Hence, the case (P < P c ) generates a classic swallow-tail phenomenon on the G-T diagram (Fig. 7), implying a first-order small/large BH phase transition.In other words, the non-smooth points on the isobaric curves related to the G − T diagram are those of the extreme points on the isobaric curve in the T − r diagram.In the monitoring phase of the considered scenario, the swallowtail gradually decreases in size, eventually disappearing as the pressure increases.Due to the increasing pressure, the extreme points of the T on the isobaric curve move closer to each other to converge at the critical pressure (P = P c ), and hence, any kind of first-order phase transition in the system disappears completely.For P > P c , G becomes a monotonic func-  tion of T, indicating that no phase transition occurs in the system.Some meaningful information from Gibbs Free Energy is predicting the local stable BH system.By the way, the sign of G claims the system is locally stable when (G < 0) and unstable when (G > 0).For that reason, Fig. 8 shows graphically the behavior of Gibbs free energy as a function of the horizon radius.Thus, it is observed for all the parameter space values that the BH system, in the overall view of the horizon radius, is locally stable and, in contrast, unstable for small values of the horizon radius.Alternatively, one can similarly inspect the same critical study from Fig. 7 as previously mentioned.

A. Critical exponents
Critical exponents perfectly describe the behavior of physical quantities in the vicinity of the critical point.Indeed, critical exponents are independent of such physical systems and can be considered quasi-universal parameters.It is handy to introduce the following notations: in which the critical thermodynamic volume V c is associated with the critical event horizon radius r c by The critical exponents are explicitly defined in the following way: The exponent α describes the behavior of specific heat at a constant volume.It is easy to conclude that the entropy S is independent of the Hawking temperature T, thus hence it can be deduced that α = 0.The exponent β describes the behavior of the order parameter near the critical point.It is then possible to expand the equation of state near the critical point in where F = (A + 1)(3A(β + 1) Accordingly, in terms of numerical results, the dependencies of the C i coefficients on the Q and β parameters are presented in Tab.I I I.
As the pressure remains constant during the phase transition, one can have where ω s and ω l are the reduced volumes of the small and large BHs, respectively.In addition, Maxwell's equal area law is simply given by the following formula and considering the first derivative so that Exploiting Eqs. ( 87) and ( 88) yield the following finding: from which, with Eq. ( 86), it is possible to find an explicit link between ω l and ω s in the following form: where the argument under the square root function remains positive.A quick look at Eq. ( 90) yields the desired results, namely which provides λ = 1/2.The exponent γ describes the critical behavior of the isothermal compressibility κ T given explicitly by giving rise to γ = 1.The exponent δ is in charge of describing the critical behavior of Eq. ( 79) on the critical isotherm T = T c .So, the shape of the critical isotherm is defined at t = 0 providing the following finding: ) which easily proves δ = 3.
From the aforementioned findings, it is clear that the four critical exponents are precisely the ones obtained previously for charged AdS BHs.This in fact demonstrates that the MCG does not alter the critical exponents, similar to the quintessential dark energy.Thus, the universality profile of VdW-like phase transitions and the values of critical exponents for AdS BHs have been ascertained.

VI. THERMAL GEOMETRIES
Based on the arguments set out in the introduction, this section focuses on providing an appropriate analysis of the phase transition for our charged BH solution using GT tools.Attention is paid to examining the shape of the metric in regard to Ruppeiner, Weinhold, HPEM and Quevedo theoretic.So, the Weinhold shape can be given in mass representation as [112] g For a charged AdS BH, the line element appears as follows: Or, in terms of mass matrix representation, the formulation is given explicitly as Similarly, in Ruppeiner formalism one considers entropy as basic thermodynamic potential, It is worth noting that the Ruppeiner metric is linked to the Weinhold metric through a conformal transformation, yielding the following defining expression [113]: Or, in terms of mass matrix representation, one has It should be noted that the treatment will be carried out in terms of pressure, which is related to the derivative of the AdS length by the following expression : In order to properly examine an interesting discussion of the GT tools in the essence of our BH solution in extended phase space, Figs.9-10 provide the necessary analysis in this regard.To begin with, the Weinhold and Ruppeiner structures provide an analysis concerning a certain gauge parameter space, which is illustrated in Fig. 9.It is observed that the zero of the heat capacity yields a divergent scalar curvature, as in the case where the Weinhold scalar curvature has a negative divergent point.This negative divergent point coincides with the physical limitation point of the heat capacity at r ⋆ 1 = 1.63453.Moreover, the application of the Weinhold geometry is compatible with the realm interpretation of the GTs tools.However, the Ruppeiner scalar curvature has three divergent points, one of which is positive.With care, the positive divergence point is shown to coincide with a physical limitation point (the root of the heat capacity at 0) at r div 1 = 1.94108.Furthermore, another coincidence is demonstrated between the Ruppeiener scalar curvature and the phase transition critical point at r ⋆ 1 = 1.63453.Therefore, the Ruppeiener structure offers valuable information.
To get a better approach to the BH phase transition analysis, it might be useful to apply geometrical thermodynamics, such as the HPEM and Quevedo tools.In particular, the Quevedo metric is expressed as follows where Within this structure, Φ, I b , and E a represent the thermodynamic potential and intense and extended variables, respectively, On the other hand, the generalized HPEM metric labeled by n extended variables is given in such a way as [48,50,113] (103) Here, ξ(ξ ̸ = S), M S = ∂M ∂S , and M SS = ∂ 2 M ∂S 2 are extensive parameters, respectively.So, the background of metrics is easily expressed as follows [48,50,113]: By the way, Ricci scalars employed in making these mea-surements have the following denominator [50,113] denom To examine, at a graphical level, the behavior of the three applying GT tools, Fig. 10 presents the relevant results.With regard to the application of the HPEM geometry, the results seem very interesting.This fascinating observation is due to the presence of three divergent points of HPEM scalar curvature, one of which is negative.It is clear that the negative divergent point is perfectly aligned with the phase transition critical point.Furthermore, it is clearly demonstrated that the HEPM scalar curvature reaches zeros of the heat capacity at r.In turn, the application of the Quevedo class provides two divergent points, one of which is negative.Unlike the HPEM scalar curvature, the divergence of the scalar curvature is at the zero point of the heat capacity.Moreover, the negative divergent point coincides with the phase transition critical point.What remains to be discussed is

VII. CONCLUSION
A number of cosmological challenges aim to reveal the unified dark fluid as an unknown energy component of the universal dark sector.What is more, the unified dark fluid is a hybrid of dark matter and dark energy, making it a candidate for the so-called Chaplygin gas.Spurred by this insight, we constructed an exact static, spherically symmetrically charged AdS BH solution with a surrounding MCG.Taking advantage of the MCG equation of state p = Aρ − B/ρ β , we established its energy density with respect to the radial coordinate.To verify the singularity and uniqueness of the BH, we have performed an examination based on scalar invariants.This showed that the BH solution is singular for a given parameter space and absolutely unique.Furthermore, we have investigated the classical energy conditions for the MCG and concluded that it admits the null and weak energy conditions but violates the strong energy condition.Further, the BH solution is graphically represented according to the permitted values on the parameter space.
As part of this work, we have inspected, within the extended phase space, different thermodynamic quantities, and the modified first law of thermodynamics, along with the Smarr relations, have been established.The frame of the extended phase space, together with the statement of BH chemistry, is helpful to consider the P − V criticality.Therefore, we have managed the P − V criticality for the charged AdS BH endowed with MCG structure and discovered a generic small/large BH phase transition, which is analogous to the liquid/gas phase transition of VdW.Furthermore, the critical exponents are shown to be similar to those of the VdW fluid.Furthermore, the heat capacity at constant pressure is said to show the thermal local stability of the BH solution.Indeed, we have shown, according to the sign of the heat capacity, that our BH solution is locally thermally stable.
Among ordinary predictions of the phase transition critical point, GT tools are alternatively used to predict certain correlations.In this way, we have applied in the extended phase space some tools from the GT background, such as Weinhlold, Ruppeiner, HPEM, and Quevedo classes I and II.The examination showed consistent results between the GT scalar curvature and either the physical limitation point (the root of the heat capacity at 0) or the critical phase transition point of the heat capacity (divergent point).In essence, the microstructures are thereby defined precisely for all classes of GT tools and according to the curvature scalar sign.This work raises certain questions and will be a subject for revealing thermodynamic topology or discovering aspects such as the Joule-Thomson expansion.On the other hand, quasi-normal modes and shadow behaviors, as well as the angle of deflection, are topics that could be addressed in the future.

B = 0. 1 BMFIG. 1 :
FIG.1: Variation of the BH metric function(26) with respect to r for various values of the parameter space.

FIG. 3 :
FIG. 3: Variation of the Hawking temperature T + (45) with respect to the horizon radius r + .

FIG. 4 :
FIG. 4: Variation of heat capacity C P (65) as a function of r + for different values of pressure and for fixed values in parameter space.

10 FIG. 5 :
FIG.5: Isobaric curve T-r + diagrams of the BH system for certain values of the parameter space.

20 FIG. 6 :
FIG.6: Isotherme curve P-r + diagrams of the BH system for certain values of the parameter space.

10 FIG. 7 :
FIG. 7:The G − T diagram of the black hole system for certain values of the parameter space.

20 FIG. 8 :
FIG.8:The G − r + diagram of the black hole system for certain values of the parameter space.

TABLE I :
The physical limitation points with A = 0.1.

TABLE III :
Numerical sets for critical physical quantities and coefficients C and C i in P − V critical behavior with A