Thermodynamics of deformed AdS-Schwarzschild black hole

By implementing the gravitational decoupling method, we find the deformed AdS-Schwarzschild black hole solution when there is also an additional gravitational source, which obeys the weak energy condition. We also deliberately choose its energy density to be a certain monotonic function consistent with the constraints. In the method, there is a positive parameter that can adjust the strength of the effects of the geometric deformations on the background geometry, which we refer to as a deformation parameter. The condition of having an event horizon limits the value of the deformation parameter to an upper bound. After deriving various thermodynamic quantities as a function of the event horizon radius, we mostly focus on the effects of the deformation parameter on the horizon structure, the thermodynamics of the solution and the temperature of the Hawking–Page phase transition. The results show that with the increase of the deformation parameter: the minimum horizon radius required for a black hole to have local thermodynamic equilibrium and the minimum temperature below which there is no black hole decrease, and the horizon radius of the phase transition and the temperature of the first-order Hawking–Page phase transition increase. Furthermore, when the deformation parameter vanishes, the obtained thermodynamic behavior of the black hole is consistent with that stated in the literature.


I. INTRODUCTION
Before the 1970s, there was no compelling reason to study the thermodynamics of black holes, until Hawking's area theorem changed this view [1].Then, Bekenstein linked Hawking's area theorem with the second law of thermodynamics by assigning entropy to black holes [2].Using such similarity and assuming the event horizon of black holes as a Killing horizon, the four laws of black hole mechanics were formulated [3].However, with the classical view that black holes absorb all matter and energy and emit no radiation, attributing temperature and entropy to black holes was questionable.But by using quantum effects, Hawking showed that black holes have radiation with the spectrum of a black body with a certain temperature [4].After the formulation of the standard laws of thermodynamics of black holes, several interpretations and investigations were carried out in the thermodynamics committee about black holes, such as entropy and temperature, the first law of thermodynamics, and relations between extensive and intensive quantities, see, e.g., Refs.[5,6] and references therein.
The thermodynamic laws of ordinary materials are associated with a pressure-volume term, however those for black holes do not contain such a term.To overcome this inconsistency, an idea based on the immersion of black holes in a background with a negative cosmological constant was presented [7].In this view, pressure can be considered equal to a negative cosmological constant.In this regard, in the context of black hole thermodynamics, black holes with asymptotic anti-de Sitter (AdS) are more interesting.The idea of using the cosmological constant as the thermodynamic pressure [8,9] generalized the first law of black hole thermodynamics and formed an extended phase-space, see, e.g., Refs.[10,11].Hence, in the presence of the cosmological constant, the mass of a black hole does not represent the internal energy of the thermodynamic system and is interpreted as the gravitational version of the enthalpy [12].
The new perspective of black hole thermodynamics based on the new interpretation of the cosmological constant and black hole mass has led to new phenomena associated with black holes.Accordingly, topics such as thermodynamic equilibrium and phase transition can be addressed.One of the most important thermodynamic behaviors is the phase transition between sufficiently large (compared to the AdS radius) black holes and an environment full of special radiation.This is known as the Hawking-Page phase transition [13], which is considered as a confinement/deconfinement phase transition in boundary conformal field theory [14].
Furthermore, black hole thermodynamics has wide range of applications compared to the Einstein gravity.Even, the investigation of more complex behaviors in the issue of the thermodynamic phase transition of black holes in modified gravitational theories has received much attention, see, e.g., Refs.[15,16].In general, finding a black hole as a solution in a gravitational theory leads to the investigation of thermodynamic behavior from the researchers' point of view.Subsequently, more studies have been conducted and other novel phenomena and phase structures have been observed, see, e.g., Refs.[17][18][19][20][21][22].
The non-linear nature of equations in gravitational the-ories leads to approximate methods and, of course, innovative proposals for finding analytical solutions, which is always an important task.In this regard, one of the attractive approach for searching and analyzing solutions of the gravitational equations is the gravitational decoupling (GD) method, which serves as a useful tool, see, e.g., Refs.[23][24][25] and references therein.Through this method, the known solutions of the standard gravitational action can be extended (with some minimal set of requirements) to additional sources and the domain of modified gravitational theories.The implementation of the GD method makes it possible to decouple the gravitational equations of sources into two parts, one for the standard field equations and one for an additional gravitational source.This method has been applied in the context of the Randall-Sundrum braneworld [26,27] and has also been extended to investigate in other gravitational issues, including new black hole solutions, see, e.g., Refs.[25,[28][29][30][31][32][33][34][35][36][37][38][39][40].In the present work, while employing the GD method, we intend to extend the AdS-Schwarzschild vacuum solution in the presence of a generic gravitational source, which satisfies the weak energy condition.Thereafter, we try to find possible deformed AdS-Schwarzschild black hole solution and then investigate its thermodynamic properties and Hawking-Page phase transition while varying the relevant parameters.
The work is organized as follows.In Sec.II, while briefly reviewing and implementing the GD method, we introduce the desired gravitational action.In Sec.III, we derive the deformed AdS-Schwarzschild black hole.In Sec.IV, we scrutinize the structure of the horizon and calculate the thermodynamic quantities in the extended phase-space to investigate the horizon structure, the thermodynamics of the solution and the temperature of the Hawking-Page phase transition.Finally, we summarize the results in Sec.V.

II. MODIFIED GRAVITATIONAL ACTION AND DECOUPLING FIELD EQUATIONS
We consider the following action in four dimensions in which an additional general Lagrangian term is added to the standard Einstein gravity with the cosmological constant, namely where g is the determinant of the metric, R is the Ricci scalar, Λ is the cosmological constant, L m is the usual matter Lagrangian and L X represents Lagrangian for any other matter or new gravitational sector beyond general relativity, e.g.Lovelock gravity and/or new other scalar/vector/tensor field(s).Meanwhile, κ = 8πG N /c 4 , however through the work, we use the natural units ℏ = 1 = c, and also employ the (−, +, +, +) signature.
The variation of the action with respect to the metric gives the field equations where G µν is the Einstein tensor and T (tot) µν represents the total symmetric energy-momentum tensor as wherein and As an ansatz, we demand the solution to be the spherically symmetric and static spacetime, namely ds 2 = −e ν(r) dt 2 + e µ(r) dr 2 + r 2 (dθ 2 + sin 2 θ dϕ 2 ), (6) in the spherical coordinates.Hence, metric (6) must satisfy the field Eqs.(2), namely where prime denotes derivative with respect to r and T (tot)3 3 = T (tot)2 2 due to the spherical symmetry.However, we effectively identify the energy-momentum tensors as where (ε, ϵ and E) are energy densities, (p r , p r and P r ) are radial pressure densities, and (p t , p t and P t ) are tangential pressure densities.In general, when pt ̸ = pr , these definitions clearly indicate an anisotropy case.
In addition, the conservation equation with metric (6), is a linear combination of Eqs. ( 7)-( 9).However, in terms of those two sources in relation (3), while using relations (10) Now, to solve the system Eqs.( 7)-( 9), we implement the GD method, which is explained in detail in Ref. [24], of course, we mention the necessary steps very briefly.We assume that the solution of Eqs. ( 7)- (9) for the source T (m) µν when T (X) µν = 0 is a general static spherically symmetric one as The energy-momentum tensor T (m) µν is conserved with this metric, i.e.
where ∇ µ is calculated according to metric (15).We also assume that the effects of the presence of the source T (X) µν on solution (15) are in the geometric deformation and where f (r) and g(r) are geometric deformations that alter the radial and temporal metric components, respectively, and the constant α is a free positive parameter, which can somehow adjust the strength of the effects on these components simultaneously.Moreover, having α guarantees that solution (15) is recovered in the limit α → 0. We refer to the parameter α as a deformation parameter.
Then, substituting decompositions (17) and ( 18) into the original field Eqs. ( 7)-( 9) causes those to be separated into two decoupled sets of equations.One set of the field equations is for the standard energy-momentum tensor Another set of the field equations is for the source T (X) which clearly shows that when the deformation parameter is zero, this source vanishes.
At this stage, one can investigate the field Eqs. ( 19)- (21) to determine (T (m) µν , ζ, λ), and then solve the field Eqs. ( 22)- (24) to specify (T (X) µν , g, f ).In other words, both sets, although separated, remain gravitationally connected, i.e. to solve the second set, one needs first to solve the first set.However, since the number of unknowns in each set of equations is more than the number of independent equations, it is necessary to apply additional conditions and/or extra relations, e.g. the equation of state.More explicitly, after solving the first set, five unknowns remain in the three equations corresponding to the additional source.Hence, to solve this set, we need to impose two constraints, which we will perform in the next section.

III. DEFORMATION OF ADS-SCHWARZSCHILD BLACK HOLE
In order to deform AdS-Schwarzschild black hole for action (1), we start in the absence of L X and with the solution of the exterior Schwarzschild (i.e., in vacuum with T (m) µν = 0) in the AdS background.In this regard, the field Eqs. ( 2) give for the region r > R, where R is the surface of the selfgravitating system, M is the ADM mass and l = 3/|Λ| is the AdS radius.
In order to have black holes with a well-defined event horizon structure for action (1), the sufficient condition (i.e., the Kerr-Schild condition) can be imposed [25] on metric (6).Actually, by adding this ansatz, the deformed black hole metric is a metric that respects the symmetries.Accordingly, a direct consequence of Eqs. ( 7) and ( 8) plus condition (26) is that the total source T When T (m) µν = 0, this relation yields the constraint i.e., in this case, with positive energy density of the source T (X) µν , only its negative radial pressure density is allowed.Furthermore, by deriving the conservation of the energymomentum tensor of sources with metric (6) while using decompositions ( 17) and ( 18) as well as using relations (11) and ( 16), the first equality in relation (25), and condition (26), gives 1 This relation, while assures the decoupling of the field equations, shows the exchange of energy between the sources.Moreover, from relation (29), it is obvious that there is no energy exchange in the absence of µν , or in the special case where ϵ = −p r , and/or when the geometric deformation g is constant.
From now on, for simplicity, we consider the background geometry to be the AdS-Schwarzschild vacuum, i.e. when T (m) µν = 0 and we have solution (25), but henceforth we do not mention the subscript.In fact, in the continuation, our aim is to explore the back-reaction of a static, spherically-symmetric energy-momentum tensor on the four-dimensional AdS-Schwarzschild black hole.Of course, it should be noted that in the mentioned case, one can get the solution without resorting to the GD method.However, we utilize the GD method as an alternative approach to solve this case to better observe the effect of moving away from the AdS-Schwarzschild solution and also to benefit from varying the deformation parameter when studying the thermodynamic properties.In other words, the merit of the GD method, in general, is that it is easy any interpretation in terms of energy exchange as well as in terms of superposition of configurations.
Accordingly, in the next step, by substituting condition (26) and the AdS-Schwarzschild solution (25) into decompositions (17) and (18), we obtain Hence, the line element (6) becomes Therefore, to completely specify the line element (6) in this case, we need to specify the behavior of function B(r) (or indeed g(r)).
1 It is instructive to emphasize that the use of the first equality in relation (25) and condition ( 26) is necessary to obtain the result of Eq. ( 29).
Now, using Eqs.( 22) and ( 30), the differential equation governing the function B(r) reads ) Eq. (33) shows that the behavior of function B(r) depends on the function of energy density E, which in turn depends on the choice of T (X) µν as the additional source.
To proceed, we resort to the weak energy condition for the T (X) µν , namely In the special case under consideration, due to constraint (28), condition (35) holds.Considering Eq. ( 33), condition (34) puts the constraint on the function B(r).While using the conservation Eq. ( 14) plus constraint (28), condition (36) imposes the constraint on the function E(r).
Once again, to continue the process, we deliberately choose the energy density function of the additional source as the special monotonic function2 which satisfies3 both condition (34) and constraint (38), and has rapid asymptotic decay.In fact, considering Eq. ( 22) and choosing function (39), we have selected a special function for the geometric deformation f (r) and in turn the function B(r).In function (39), r ̸ = −β, where β is a constant parameter with length dimension that controls the behavior of E(r) at r = 0, i.e. it is necessary to avoid the central singularity.We refer to the parameter β as a control parameter.In addition, the role of the deformation parameter in function ( 39) is consistent with Eq. ( 22).Meanwhile, substituting (39) into Eq.(38) gives and in turn, by using relation (12) and constraint (28), the trace of the energy-momentum tensor T (X) µν is obtained to be Thus, under the mentioned considerations, the energymomentum tensor T (X) µν depends only on a single function of the radial coordinate, its energy density function.Now, by substituting function (39) into Eq.( 33) and solving the differential equation for the function B(r), we obtain where c 1 is an integration constant.However, to have the background AdS-Schwarzschild solution (25) in the limit α → 0 (i.e., when B(r)− −− → α→0 1, and consistent with definition ( 32)), it is necessary to set4 Hence, the metric function reads which preserves the asymptotic behavior in the form of the AdS solution.Moreover, the asymptotic behavior of the term resulting from the additional source is 1/r 2 , i.e. like a Maxwell field.
Therefore, the line element for the region outside a self-gravitating system filled with the energy-momentum tensor T (X) µν , while choosing the behavior of its energy density as function (39), is We refer to the black hole solution of this metric as a deformed AdS-Schwarzschild black hole, and in the next section, investigate its horizon structure and thermodynamics.

IV. STRUCTURE OF HORIZON AND THERMODYNAMICS OF DEFORMED ADS-SCHWARZSCHILD BLACK HOLE
First of all, the metric function (44) in the limit β → 0 becomes the metric function corresponding to the wellknown charged AdS black hole, i.e.
Therefore, √ α can be interpreted as the electric charge of black hole.In this particular case, the thermodynamic behavior of the black hole is reminiscent of that investigated in Refs.[7,41].However, in these references, with the electric charge as a conserved charge of the theory, the possibility of variable electrostatic energy with an electric potential has been considered, while we consider the deformation parameter only as a free parameter. 5Moreover, we emphasize that in the case of β = 0, relation (41) indicates that the trace of the additional energy-momentum tensor is zero, and utilizing constraint (28) and relation (40) gives P eff ≡ 1 3 (P r + P t + P t ) = 1 3 E.These results plus the second footnote confirm that the additional energy-momentum tensor is the one of a Maxwell field and its solution is simply the electrically charged static AdS black hole.
On the other hand, we deliberately put the control parameter, i.e. β, by hand in function (39) so that β ̸ = 0 prevents the divergence of the energy density E. Therefore, it is better to check the effect of the control parameter on the behavior of the metric function.For this purpose, we expand the metric function (44) around the zero of its argument to obtain (47) Hence, with the value of the metric function F (r) has no singularity at r = 0. Further investigation of this particular case may lead to a regular black hole, however we leave this investigation to another work.That is, although relation (48) is necessary to avoid the divergence of the metric function, we do not limit our investigation to it in this work.However, by the simple transformation r → β + r, the control parameter is turned on.Hence, the energy-momentum tensor T (X) µν no longer describes a pure Maxwell field and the corresponding metric is modified.Actually, we are interested to have the black hole solution of metric (45) in the presence of non-vanishing β.For this purpose, the equation governing the radius of the event horizon r h is determined by the larger root of equation which is a sixth order equation.Its solution is complicated and not useful for our purpose.Hence, instead of solving it analytically, we restrict ourselves to specifying its important consequences by resorting to its various corresponding diagrams in FIG. 1 with fixed values of its constants.For each of the curves in this figure, the intersection of the metric function with the horizontal axis determines its corresponding event horizon radius.
In this regard, this figure shows that the condition of having an event horizon limits the value of the deformation parameter to an upper bound, α max .In general, the value of α max depends on the values of other parameters, namely M , l and β.For instance, FIG. 1 indicates that when M = 1, l = 1 and β = 1, the radius of the event horizon decreases with increasing the value of α, until we obtain 6 α max ≈ 6.In the continuation of this work, we pay attention to this limitation.In other words, we comply the value 0 < α < α max as a condition for having an event horizon for the deformed AdS-Schwarzschild black hole.
Alternatively, we have plotted the metric function F (r), relation (44), with respect to r, while varying the control parameter β and fixing the other parameters including the deformation parameter α in FIG. 2. This figure shows that when M = 1, l = 1 and α = 4, the radius of the event horizon increases with increasing the value of β.This effect is expected because the β parameter shifts the radial coordinate as mentioned earlier.In turn, increasing the radius of the event horizon affects the thermodynamic properties, and this indicates the influence of the control parameter on them.
To investigate the thermodynamic behavior of the obtained solution, while using Eq.(49), we first express the 6 However, using Eq. ( 49) with these mentioned values gives α mass of the black hole in terms of the radius of the event horizon, i.e.
Furthermore, the entropy of the deformed AdS-Schwarzschild black hole can be obtained from the Bekenstein-Hawking formula [2,42] as a quarter of the area of the event horizon, i.e.
where the Planck length L Pl is considered in the natural units and the second equality is due to the static and spherical symmetry.Also, the definition of black hole pressure (density) P in AdS space (in the natural units with a negative cosmological constant, see, e.g., Ref. [7]) is Then by obtaining r h and l from relations (51) and (52), and substituting into relation (50), the mass of the deformed AdS-Schwarzschild black hole can also be expressed in terms of thermodynamic quantities S and P as Hence, for the deformed AdS-Schwarzschild black hole, we can write the (generalized) first law of black hole thermodynamics in the extended phase-space as where T and V are the Hawking temperature and the thermodynamic volume, respectively.Of course, in relation (54), the presence of V dP instead of P dV indicates that the mass M is the enthalpy of black hole instead of the internal energy.Now, using relation (53), let us derive these thermodynamic quantities.In this regard, the temperature reads Of course, this relation can also be derived using the definition of the Hawking temperature in terms of the radius of the event horizon, i.e.
where we have substituted M from relation (53).Then, taking r h and l from relations (51) and ( 52), and substituting into relation (56), it reads the same as relation (55) as expected.Also, the thermodynamic volume is [43,44] that, in terms of the radius of the event horizon, obviously reads To continue investigation of the thermodynamic behavior of the deformed AdS-Schwarzschild black black hole, we apply the thermodynamic machinery suggested in Ref. [7].We assume that the black hole occurs in a canonical ensemble.We also assume that each related extended phase-space contains a fixed value of the deformation parameter.In this case, the Gibbs free energy is8 which in terms of the radius of the event horizon becomes In this regard, it is known that the condition of better thermal equilibrium globally corresponds with more negative values of the Gibbs free energy.Also, the criterion of the phase transition is related to the vanishing Gibbs free energy of the black hole.
Another useful quantity in the thermodynamic study of a black hole is the specific heat capacity at constant pressure, C P , which determines the local thermodynamic stability of the black hole.In the case at hand, it is , (61) where we have used the chain rule and relation (56).However, the local thermodynamic stability corresponds to positive C P values.
On the other hand, among the various black hole phase transitions, one of the most significant is the Hawking-Page phase transition [13].This is a study of the thermodynamics between an AdS-Schwarzschild black hole and the thermal AdS space.In this respect, we consider the constant values of β = 1 and9 l = 1 and choose several constant values of the deformation parameter within its allowed range to investigate its effects on the stability of the black hole and particularly on the Hawking-Page phase transition.
First, let us investigate the thermodynamic stability of the black hole locally through the behavior of the specific heat capacity at constant pressure.In this regard, FIG. 3 indicates the behavior of the C P function with respect to the radius of the event horizon for several values of the deformation parameter in its allowed range.The resulting curves contain a discontinuity at a certain horizon radius, such that for the chosen values it is approximately in the range of 0.55 < r h < 0.6.Moreover, the black holes with horizon radii located in the region of negative C P values are thermodynamically unstable, and a minimum value of horizon radius is required for a black hole to have local thermodynamic equilibrium.However, as the deformation parameter increases, such minimum required radius decreases.FIG. 3 also shows that black holes in the C P < 0 region have small horizon radii (which we refer to as small black holes (SBHs)) and transform into other thermodynamically allowed states.However, black holes with larger horizon radii in the C P > 0 region (which we refer to as large black holes (LBHs)) have a clear local thermodynamic stability.Hence, LBHs are thermodynamically preferred over SBHs.In addition, according to this figure, as α increases, the horizon radius of locally stable black holes decreases.Furthermore, in the region of LBHs, for all curves, as the horizon radius increases, the corresponding value of C P first decreases to a minimum and then increases.To scrutinize the thermodynamic behavior of black holes globally, we investigate the behavior of the Gibbs free energy function.
In this regard, first we have depicted the G function versus r h with fixed values of β = 1 and l = 1, and different chosen values of α in FIG. 4. It is clear that the global thermodynamic equilibrium is better achieved with less negative values of the G function.Here, as an important aspect, this figure illustrates that LBHs have less negative values of the G function.Hence, FIG. 4 confirms that in addition to local thermodynamic stability, this group of black holes also has global thermodynamic stability, while SBHs do not.Thus once again, LBHs are thermodynamically preferred over SBHs.This figure also shows that the function G is first ascending and then descending for all allowed values of the deformation parameter.However, as the deformation parameter increases, these maximum values of the G function and also the horizon radius of the phase transition (wherein G = 0) increase.Nevertheless, increasing the deformation parameter causes the thermodynamic stability of a black hole to be disturbed compared to its previous position.Also, each maximum point of the G function represents a horizon radius and a temperature.However, to better realize the phase transition and the global thermodynamic stability of black holes, we probe the behavior of the G function versus the horizon temperature in the following figure.Although it is better and instructive to first plot the temperature versus the radius of the event horizon.
In this respect, FIG. 5 shows the behavior of the temperature of horizon with respect to r h for fixed values of β = 1 and l = 1, and different chosen values of α.This figure indicates that as r h increases, the temperature of horizon first decreases and then increases for all values of α.The minimum temperature in this figure exactly corresponds to the maximum Gibbs function in FIG. 4. In general, the changing behavior of the temperature function is in accord with the behavior of the G function.Now, employing relations (56) and (60), FIG.6 shows the behavior of the Gibbs function with respect to the horizon temperature for fixed values of β = 1 and l = 1, and different chosen values of the deformation parameter.The cusp of each curve in this figure represents the maximum of the G function (as shown in FIG. 4) and the minimum temperature (as shown in FIG. 5).At temperatures below the minimum temperature, there are no black holes.FIG.6 indicates that as the temperature increases from its minimum, there are two branches of black holes.In both branches the G function decreases, which is in accordance with FIG. 4. The right/upper branches contain SBHs with C P < 0, while the left/lower branches have C P > 0 and include black holes with intermediate horizon radii (which we refer to as intermediate black holes (IBHs)) with positive values of the G function, the intersection of curves with the horizontal axis (the phase transition with G = 0), and LBHs with negative values of the G function.
Since thermodynamically, smaller (and even negative) values of the G function are preferred, there are actually two global thermodynamically stable phases.The first phase in which G = 0 represents the immersion environment of black holes and includes the thermal radiation (the thermal bath filled with the cosmological constant).This region is preferred with respect to IBHs.The sec-

V. CONCLUSIONS
Initially, we considered the Einstein-Hilbert action with the presence of the cosmological constant and a standard matter source plus any additional gravitational source.The aim of this work is to find the black hole solution(s) for such an action and to search for the corresponding thermodynamic behaviors.To proceed, we have employed the GD method, which serves as a useful tool for searching solutions to the gravitational equations.Then, we have taken the background as the AdS-Schwarzschild vacuum solution, and have looked for the static spherically symmetric solution(s) when the additional source is also present.By implementing the GD method, we have used two geometric deformation functions to alter the radial and temporal components of the background metric.Meanwhile, in process of this method, there is a common positive parameter that can adjust the strength of the effects on these components simultaneously, and we refer to it as a deformation parameter.When this parameter vanishes, the background solution is recovered.Through the GD method, the field equations are separated into two decoupled sets of equations for each source, although they remain gravitationally connected.
In continuation, after solving the set of equations for the background, five unknowns remain in the three equations corresponding to the additional source.Hence, to solve this set, we imposed two constraints.Actually, we have assumed that the additional source obeys the weak energy condition, and we have also deliberately chosen its energy density function to be a special monotonic function proportional to the inverse of the distance to the fourth power.Also for consistency, this function is proportional to the deformation parameter and includes a control parameter to prevent it from diverging at the center.Actually, in the absence of the control parameter, we showed that the additional energy-momentum tensor is the one of a Maxwell field.Moreover, to have black holes with a well-defined event horizon structure, we have imposed the Kerr-Schild condition that the radial and temporal components of the solution to be inverses of each other regardless of their signs in the signature.We refer to the obtained solution as a deformed AdS-Schwarzschild black hole.The solution found, although complicated, turns out to be analytical and with quite interesting features.Then, the focus of the work has been to investigate the horizon structure, the thermodynamics of the solution and the Hawking-Page phase transition mainly by varying the deformation parameter.However, within the work (including figures), we considered several values of the deformation parameter with fixed values of other constants (hence, actually at constant pressure).
Since we intended to consider the black hole solution, we first plotted the metric function versus the distance.For each of the curves in this figure, the intersection of the metric function with the horizontal axis specifies that there is an event horizon radius.However, the figure shows that the condition of having an event horizon limits the value of the deformation parameter to an upper bound, as we have also shown through the corresponding derivation.Hence, we confined the investigation to vary the deformation parameter up to its upper bound value.
Next, we assumed that the black hole occurs in a canonical ensemble and wrote the first law of black hole thermodynamics in the extended phase-space and then derived various thermodynamic quantities as a function of the event horizon radius.We also assumed that each related extended phase-space contains a fixed value of the deformation parameter as a free parameter and not a thermodynamic quantity.Thereafter, to determine the thermodynamic stability of the black hole locally, we plotted the obtained heat capacity at constant pressure versus the radius of the event horizon.The figure shows that a discontinuity occurs between the negative and positive values of the heat capacity at constant pressure for all its curves.The negative region contains SBHs that are locally thermodynamically unstable.Whereas, the positive region contains LBHs that are locally thermodynamically stable.In other words, for a black hole to have local thermodynamic equilibrium, a minimum value of horizon radius is required, and as the deformation parameter increases, this minimum required radius decreases.
In addition, in the region of LBHs, for all curves, as the horizon radius increases, the corresponding value of C P first decreases to a minimum and then increases, while as the deformation parameter increases, the horizon radius of locally stable black holes decreases.Furthermore, to determine the thermodynamic stability of the black hole globally, we plotted the obtained Gibbs free energy versus the radius of the event horizon.This figure shows that the Gibbs function starts from the region of positive values, and with the increase of the horizon radius, it first increases to a maximum and then decreases to more negative values after crossing the horizontal axis (i.e., the phase transition point).The figure also illustrates that LBHs have less negative values of the Gibbs free energy.Hence, in addition to local thermodynamic stability, this group of black holes also has global thermodynamic stability and is thermodynamically preferred over SBHs.Moreover, with the increase of the deformation parameter, the maximum values of the Gibbs function and the horizon radius of the phase transition increase.Nevertheless, increasing the deformation pa-rameter causes the thermodynamic stability of a black hole to be disturbed compared to its previous position.
Then, to better realize the phase transition and global thermodynamic stability of black holes, we plotted the temperature versus the event horizon radius and the Gibbs free energy versus the horizon temperature.The figures indicate that there is a minimum temperature that exactly corresponds to the maximum of the Gibbs function.As the deformation parameter increases, the minimum temperature decreases and no black hole exists below this minimum temperature.The second figure illustrates that SBHs are in one branch of the curves, and in the other branch are IBHs, the phase transition point (i.e., thermal radiation) and LBHs.The proximity of the phase transition between the thermal radiation medium and the deformed AdS-Schwarzschild LBHs is known as the Hawking-Page phase transition.As the deformation parameter increases, the temperature of the first-order Hawking-Page phase transition increases.
In the special case of vanishing the deformation parameter, the obtained thermodynamic behavior of the black hole is consistent with that stated in Ref. [7].Also, we showed that in the special case of vanishing the control parameter, the obtained metric function corresponds to the charged AdS black hole, which was investigated in Refs.[7,41] with the square of the deformation parameter as the role of electric charge.In addition, we showed that increasing the control parameter increases the radius of the event horizon, which in turn affects the thermodynamic properties.Furthermore, we indicated that for a certain non-zero value of the control parameter, the obtained metric function has no singularity at the center.Further investigation of this particular case may lead to a regular black hole, however we leave these investigations to another work.
µν must satisfy the equation of state pr = −ε.

FIG. 1 .
FIG. 1. [color online] Using relation (44), the figure schematically (i.e., scale-free) illustrates the behavior of the function F (r) with respect to r for fixed values of M = 1, l = 1 and β = 1, and different chosen values of α.As α increases, the radius of the event horizon decreases, until αmax ≈ 6.

FIG. 2 .
FIG. 2. [color online] Using relation (44), the figure schematically (i.e., scale-free) illustrates the behavior of the function F (r) with respect to r for fixed values of M = 1, l = 1 and α = 4, and different chosen values of β.As β increases, the radius of the event horizon increases.

FIG. 3 .
FIG. 3. [color online]Using relation (61), the figure schematically (i.e., scale-free) illustrates the behavior of the CP function with respect to r h for fixed values of β = 1 and l = 1, and different chosen values of α.In the CP < 0 region, the black holes do not have local thermodynamic equilibrium.In the CP > 0 region, as α increases, the black holes with smaller horizon radii have local thermodynamic equilibrium.Also, in this region for all curves, as r h increases, the corresponding value of CP first decreases to a minimum and then increases.

FIG. 4 .
FIG. 4. [color online] Using relation (60), the figure schematically (i.e., scale-free) illustrates the behavior of the function G with respect to r h for fixed values of β = 1 and l = 1, and different chosen values of α.The intersection points of the curves with the horizontal axis (i.e., G = 0) indicate the horizon radius of the phase transition.As α increases, the r h of the phase transition and the maximum value of the G function also increase.More interestingly, in the G < 0 region, black holes with larger r h have less negative G function.

FIG. 5 .
FIG. 5. [color online] Using relation (56), the figure schematically (i.e., scale-free) illustrates the behavior of the temperature of horizon with respect to r h for fixed values of β = 1 and l = 1, and different chosen values of α.As r h increases, the horizon temperature first decreases to a minimum and then increases.

FIG. 6 .
FIG. 6. [color online] The figure schematically (i.e., scale-free) illustrates the behavior of the G function with respect to the horizon temperature for fixed values of β = 1 and l = 1, and different chosen values of α.A deformed AdS-Schwarzschild black hole, like a Schwarzschild black hole, exhibits a phase transition with a thermal radiation medium.The region of the thermal radiation is indicated by an arrow for each curve, and the end of each arrow indicates the temperature of the first-order Hawking-Page phase transition.