Thermodynamic consequences of well-known regular black holes under modified first law

We modify the first law of thermodynamics of regular black hole of the non-minimal Einstein–Yang–Mill theory with gauge field of magnetic Wu–Yang type and a regular black hole which is associated with the cosmological constant by the surface tensions. The corresponding Smarr relations are also satisfied for them. We calculate the Gibbs free energy to discuss the global stability of the black holes. It is found that Gibbs free energy exhibits the positive behavior in most of the cases which leads to globally thermodynamically stability of both black holes.


Introduction
In modern theoretical physics, black hole (BH) is an interesting subject. At all times, BH is associated with a singularity inside it. Einstein's theory of space, time and gravity empowers the existence of singularities. According to the General Relativity (GR), a singularity is a definite point in the spacetime where explicit quantities even the curvature of spacetime is infinite. The problem of singularities is one of the exceptional problems and GR made the prediction about the singularities in the interior of BHs. Firstly, it was anticipated that the both of these singularities does not naturally exist but occurs as a result of the preparative procedure of finding the harmony in the solutions. Penrose and Hawking's theorems of singularities lead us to believe that they are inevitable [1][2][3][4][5][6]. As the initial exact BH solutions in GR predicts the singularity in the interior of an event horizon. Sakharov [7,8] and Gliner [9] proposed that for the sources like matter which have a de-Sitter core in the heart of the space-time its singularities could have been avoided and then Bardeen availed this theory to put forward the initial solution of the static spherically symmetric regular black hole (RBH), called "Bardeen BH". Bronnikov and Melnikov particularize a e-mails: jawadab181@yahoo.com; abduljawad@cuilahore.edu.pk b e-mail: amna.khawer22@gmail.com the detail of the many types of the RBH [10], explained the geometries of RBH in a nutshell and revealed the actuality of the black-universe (BU) solutions in GR in agreement with the minimally coupled phantom fields. By coupling the GR to non-linear electrodynamics (NED) many different kinds of RBH solutions have been derived [11][12][13][14][15]. It is mention here that RBHs asymptotically behave like an ordinary charged Reissner-Nordstrom BH solutions and all fields and curvature invariants are regular everywhere instead of usual singularities at their center. These solutions and singularity theorems have no contradiction at all between them.
Ayon-Beato and Garcia, in GR, introduced a proposal of exact nonsingular BH solution fulfilling the weak energy condition (WEC), which is produced by the new NED combined to gravity [11]. Bambi and Modesto proposed the solutions for the rotating RBH [16]. The spherically symmetric radiating and rotating solutions have been proposed by using the Newman-Janis algorithm [17]. Fan and Wang, in the gravity model composed those solutions of BHs which bear magnetic charge on them. Strategy of formulating the solutions of the electrically charged BHs has been illustrated and the thermodynamical possessions has been studied along with the formulation of the first law of thermodynamics. Also generalization of the construction of AdS BH solutions has been done [18].
Hawking and Penrose, during the late 1960s, put into practical a new complicated mathematical model which led to both proposing and proving the numerous theorems of singularity [3]. There is a parallelism between the ordinary 2nd law (OSL) of thermodynamics, which is stated as "the entropy of the closed system never decreases" and Hawking's theorem, whose statement is that "the surface area of BH never decreases" [19,20], and Bekenstein was the first to analyzed the resemblance between these two [21]. Then in the 1972, Hawking-Bekenstein debate started until 1974, when Hawking discovered the radiation of BH near event horizon produced by it's quantum effects and named it, Hawking's radiation [22].
The discovery about the area of the event horizon is irreducible lead towards a successful conclusions in GR and made up an essential demonstration that the BHs have thermodynamical properties. Thermodynamical properties of the space-times proved to be very important and grasped the worthwhile attentions [23]. They investigated the properties for the magnetically charged regular black hole (MCRBH). Phase transitions is another well-known phenomenon which is based on the fundamental principles of thermodynamics. In the succeeding work, the phase transitions and the thermodynamical stabilities has been studied for many different BHs [24].
Kastor, Ray and Traschen, for the anti de Sitter (AdS) BH, treated the cosmological constant associated with the pressure such as P = − 8π , while treated its conjugate quantity as thermodynamic volume. The expanded first law of thermodynamics presented which has the new term in the format of effective volume times change in pressure, which indicates that the mass of AdS BH should be seen as the enthalpy of BH [25]. Kubiznak and Mann acquired this concept and investigated the thermodynamics of the charged AdS BH in the extended phase space [26]. In this way, the investigation on thermodynamics of the many BHs has been made and many significative critical phenomena were found in [15,[27][28][29][30][31][32][33][34][35][36][37][38][39][40][41][43][44][45][46][47][48][49][50][51]. Hansen et al. modified the first law of thermodynamics by using the radial Einstein's equation at the Kerr BH by the surface tension [52]. Chen and Zeng, modify the first law of thermodynamics of the Schwarzschild dS-BH at the BH horizon and the cosmological horizon. They have also derived the apparent horizon of the Friedmann-Robertson-Walker cosmology by the surface tension [53]. In recent work, for the dS and AdS space-times, Chen, Qingyu and Tao, modified the first law of thermodynamics by the surface tension at the BH horizon and cosmological horizon [54].
In this paper, our aim is to investigate the thermodynamics at the BH horizon of the two BHs (the non-minimal MCRBH and the RBH with cosmological constant in NED). In Sect. 2, we modify the first law of thermodynamics for the nonminimal MCRBH. We treat the cosmological constant to be fixed and drive the modified first law in the presence of the surface tension at the horizon. The correlated Smarr relation will be evaluated. Cosmological constant will be associated with the pressure as P = − 8π and get modified laws for nonminimal MCRBH in the stretched out phase space. In Sect. 3, we modify the first law of thermodynamics at the BH horizon treating the cosmological constant as fixed by the surface tension and as associated with the pressure (P = −r 2 + 8π(q 2 +r 2 + ) ) of the RBH with cosmological constant in NED. Sect. 4 is devoted for conclusions.

Modification of First Law of Thermodynamics for non-minimal RBH
The establishment of the BH thermodynamics came into being in between the 1960's to 1970's. The discoveries of Bekenstein and Hawking break the grounds for BH thermodynamics [22] and [55]. Recently, Astorino examined the conserved charges and the thermodynamics of the accelerating (RN-BH), where he used the phase space technique for evaluating the mass which satisfies the standard first law of thermodynamics and the corresponding Smarr relation [56].
In the mean while, Chen, Qingyu and Tao make great efforts in modifying the first law of thermodynamics with fixed and varied cosmological constant for RN-AdS and dS BHs [54].

Fixed Cosmological Constant
Consider the line-element of non-minimal MCRBH whose new exact regular spherically symmetric solution was presented by Balakin, Lemos and Zayats [57], which is given by with the metric function (r ), where λ, Q m , M, and r are non-minimal parameter of the theory, magnetic charge of the Wu-Yang gauge field, asymptotic mass of the object which appears to be the constant of integration, cosmological constant and radial variable of MCRBH, respectively. For our study, we consider λ > 0, > 0, ≤ 0, Q 2 m > 0 and M ≥ 0. The limiting case λ = 0 gives the minimal coupled magnetized RN solution along with the cosmological constant as Since, there exist a curvature singularity at r + = 0 for this MCRBH solution. Also, for λ < 0 at finite positive r + , there exists space-time curvature singularity, therefore, we are only considering the case when λ > 0. Balakin, Lemos and Zayats, analyzed the horizons as functions of the parameters [57]. Equation (2) manifest explicitly a four parameter family of exact solutions, which contains solutions along with horizons, based on the corresponding values of these four parameters. Furthermore, for the horizons, all possible solutions depend on the value and sign of the . Generally, Cauchy horizon, event horizon and cosmological horizon become apparent, in few cases these horizons can occur simultaneously with each other or three altogether.
Using Eq. (2), we can obtain the mass of the non-minimal MCRBH in horizon radius as follows For > 0, which is the asymptotically dS, the MCRBH appears to have three horizons, Cauchy, event and cosmological horizons, based on the values above mentioned four parameters. In order to find the number of horizons, an auxiliary function f (r + , λ, , Q m ) is introduced and Eq. (8) takes the form as, and In [57], Fig For = 0, that is, when all space-times are asymptotically flat. In this case, cosmological horizon doesn't exist, only Cauchy and event horizons become apparent. To sum it up, MCRBH solution can have three horizons, Cauchy, event and cosmological horizons for the positive values of cosmological constant and for = 0, only Cauchy and event horizons. Thus, we want to calculate the thermal quantities on the outer horizon and we refer it by r + . The event horizon of the non-minimal MCRBH is obtained by setting (r ) = 0 in Eq. (2), which is located at the outer horizon r + . In order to carry out the analysis of the thermodynamical properties of the non-minimal MCRBH, we use the definition of Misner-Sharp mass [58]. For the non-minimal MCRBH, the Misner-Sharp mass on the horizon can be obtained as The non-minimal MCRBH entropy and temperature are respectively, where (r + ) = ∂ ∂r | r =r + . The value of T is determined by r + , Q m and . Figure 1 illustrates the behavior of the temperature T with respect to r + for specific values of Q m and . The temperature shows the fluctuations between negative and positive values with the passage of horizon radius. For example, when = 0.01, the temperature exhibits the negative behavior for 0.55 ≤ r + ≤ 0.9, however, it is positive and attains maximum value at r + = 1.9 while approaches to zero for higher values of horizon radius (r + ≥ 1). When = 0, the temperature is negative for the small values of horizon radius, positive for the higher values and approaches to T = 0.01. For = − 0.01, the temperature approaches to T = 0.15. We observe that the temperature is higher for positive cosmological constant as compare to the negative cosmological constant and = 0. From Eq. (9), we have the expression for physical mass as By using Eqs. (9) and (10), the temperature leads to For obtaining real positive temperature, the restrictions are made to the other parameters and the conditions are as follows: The concept of horizon thermodynamics emerged from the discovery that Einstein's equations on the BH horizon can be interpreted as a thermodynamical identity. The radial Einstein equation at the horizon is needed to be evaluated essentially in order to investigate the thermodynamics and the surface tension By inserting Eq. (12) into the Eq. (9), the non-minimal MCRBH temperature can be rewritten as follows Since the entropy of a non-minimal MCRBH which is related to the area of the BH horizon is S = πr 2 + . Furthermore, the partial derivative of S is δS = 2πr + δr + . Multiplying the both sides of Eq. (16) by δS, we get The first term on the R.H.S in above equation is depend on the matter while the second term can be identified as the differential form of the Misner-Sharp mass. T r r | r + can be obtained from the Eq. (13) as, Also, Eq. (14) can also be written as, The above equation is the modified first law of thermodynamics at the non-minimal MCRBH horizon. Where σ = surface tension at the horizon and δ A can be identified as the differential form of the horizon area, whereas horizon area is A = 4S. The surface tension σ depends upon three factors such as magnetic charge Q m , non-minimal parameter λ and cosmological constant . The surface tension approaches to zero for r 6 it yields positive surface tension and if r 6 then the surface tension will be negative. The corresponding Smarr relation takes the following form It is found that this relation is satisfied for the values of E, T, S, A. It is also mentioned here that the modified first law of thermodynamics in Eq. (16) is different from the Eq.(3.8) in [26]. Conventional thermodynamics analysis requires an application of the first law of thermodynamics, known as energy analysis. In Conventional thermodynamics, in accordance with the specifications of stable equilibrium, we calculate the thermodynamic variables, such as, the heat capacity at the constant pressure and volume and the Gibbs free energy in order to analyze the local and global thermodynamic stability of the BHs in GR. Moreover, the free energy in the grand canonical ensemble, also called Gibbs free energy. The expression for the Gibbs free energy is We get G = r + 4 , using the expressions for E, T, S, σ and A. Carrying out the differential on the Gibbs free energy yields, As the effective temperature represents the temperature of a black body that would emit the same total amount of radiation, consider the effective temperature be T e f f = T − 4σ = negative. Figure 4 shows that the effective Gibbs free energy decreases with the increase of the effective temperature.

The Pressure
Here, we discuss the pressure at the horizon of MCRBH. Horizon thermodynamics takes a different approach to the problem of pressure in the BH thermodynamics. Since temperature with surface gravity yields T = κ 2π = (r + ) 4π , horizon thermodynamics is based upon the approach that the energy-momentum tensor on the horizon is interpreted as The area and the volume are inter-linked as V = Ar + 3 which generates σ δ A = PδV . By using Eq. (21), the modified law The pressure and surface tensions are directly proportional to each other, thus if pressure is positive, negative or zero then it will be in accordance with positive, negative or zero surface tension, respectively. Using the pressure, P = T r r | r + in the Eq. (13), we get, which leads to The Van der Waals equation correlate with state variables the pressure P, the volume V , the number of particles N and the temperature T , which can be written as where k is the Boltzmann constant, ν is the specific volume, b is the size of molecules of the system and a is a constant which represents the attraction between the molecules of the system. By comparing the Eqs. (23) and (25), we get a = 1 2π , b = 0 and ν = 2r + . Hence the above expression of pressure takes the following form The Gibbs free energy is which behave in accordance with varying the horizon equation of state by treating the pressure P and temperature T as independent thermodynamics quantities then horizon equation of state established a new horizon first law in [59], which is same as Eq. (22) where temperature T , volume V and pressure P were specified and the horizon entropy S and the Gibbs free energy G were the derived concepts. By applying a degenerate Legendre transformation on Eq. (25), first law can be easily regained. This condition prevail over the mystery between the heat and the work terms. In Fig. 5, we observe the behavior of Gibbs free energy for different values of temperature and fixed values of pressure. For the first case, when the cosmological constant is positive, diagram illustrates that the Gibbs free energy values are lower for the high temperatures. As temperature decreases, the Gibbs free energy increases for all the fixed values of pressure. In Fig. 6, when = 0, the Gibbs free energy remains the same with the increase in temperature. When T = 0.016, curves bounce back and Gibbs free energy increases for the lower values of temperature. For the negative value of pressure, the Gibbs free energy approaches to zero while for positive value of pressure and P = 0, the Gibbs free energy increases.
By taking the differential on the entropy S, horizon area A and the volume V , we get, which leads to here, the coefficient of δ A depends upon the magnetic charge By using the transformation P(δV ) = δ(PV ) − V δ P and differential of the Eq. (31), we obtain Thus, Eq. (33) is also the modified first law of thermodynamics of the non-minimal MCRBH. The corresponding Smarr relation leads to

Modifying Laws for RBH with Cosmological Constant
Mo Wen-Juan, Rong-Gen Cai and Su Ru-Keng [60] presented a new kind of RBH solution associated with the cosmological constant in NED. The (r ) of this metric has the following form By using (r ) = 0, we attain the BH horizon. The expression for mass M is , The q represent the electric charge of the given RBH. As we are considering the definition of the Misner-Sharp mass for investigating the thermodynamic properties of our RBH, therefore, the expression for mass is E = . By using the expression of mass M from Eq. (36), we have get a very brief expression as E = r + 2 . The temperature of this RBH with cosmological constant is, Figure 7 illustrates the behavior of the temperature and horizon radius at the fixed values of q and . We have considered the three cases of in Fig. 7. For the positive cosmological constant, the temperature is increasing for  (22), which is the modification of the first law of thermodynamics at horizon of the RBH with cosmological constant, the surface tension at the horizon becomes σ = − 3q 6 + 6q 4 r 2 + + 4q 2 r 4 + + r 4 + (q 2 + r 2 The surface tension σ depends on the electric charge q and the cosmological constant . When r 4 + = − (3q 6 +6q 4 r 2 + +4q 2 r 4 + )l 2 (q 2 +r 2 + ) 2 , the surface tension is zero, r 4 + > − (3q 6 +6q 4 r 2 + +4q 2 r 4 + )l 2 (q 2 +r 2 + ) 2 , the surface tension is positive while negative surface tension is obtained for r 4 . The corresponding Smarr relation is same as Eq. (17). The Gibbs free energy turns out to be Figure 8 represents the behavior of Gibbs free energy for the different values of surface tension, three cases will be discussed, i.e. the surface tension is positive, negative and zero. When the temperature is low, Gibbs free energy is increasing at slow rate. When T = 0.16, the Gibbs free energy increasing rapidly and temperature bounces back for the − 0.02 ≤ σ ≤ 0.02. The effective temperature for the RBH with cosmological constant is also found to be T e f f = T − 4σ = 1 4πr + . Now we will discuss the pressure for the RBH with cosmological constant at the horizon. Since the thermodynamic pressure is P = T r r |r + , therefore, the expression for the pressure is evaluated as P = − 3q 6 + 6q 4 r 2 + + 4q 2 r 4 + + r 4 + (q 2 + r 2 + ) 2 8πr 2 In the presence of pressure, the Gibbs free energy is, Figure 9 represents the behavior of the Gibbs free energy varying along with the temperature at the fixed values of pressures. We observe that when P < 0, the Gibbs free energy increases for higher values of temperature, after certain value it bounces back and temperature decreases. Once the Gibbs free energy reaches its maximum value, it starts decreasing. Both red and green curves follows this behavior. But when P ≥ 0, the Gibbs free energy increase for the different values of temperature. Now, we will consider the cosmological constant, which is correlated with pressure and investigate the thermodynamics of the RBH with the cosmological constant. The first term on the R.H.S of Eq. (14) can be rewritten as 2r + T r r | r + δS = − 3q 6 + 6q 4 r 2 + + 4q 2 r 4 + 4πr + (q 2 + r 2 2r + T r r | r + δS = − 3q 6 + 6q 4 r 2 + + 4q 2 r 4 By using the relations between the entropy, horizon area and volume, we obtain the equal sign in the above equation. It is clear that the first term depends on the electric charge and second term depends on the cosmological constant. Let σ eff = − 3q 6 +6q 4 r 2 + +4q 2 r 4 + 16πr + (q 2 +r 2 + ) 3 and P = −r 2 + 8π(q 2 +r 2 + ) , then also for RBH with cosmological constant, Eq. (14) takes the form of the Eq. (31), which further can be rewritten as Eq. (33). In these equations E o = E + PV is identified as the enthalpy and its expression is evaluated as After making analysis on the E o , we have seen that it is the Misner-Sharp mass for the RBH with cosmological constant and by taking into account the expressions of T, S, σ, A, P and V the corresponding Smarr relation is same as the Eq. (34), that we have obtained for the non-minimal RBH. We obtain the expression for the corresponding Gibbs free energy as In Fig. 10 The corresponding Gibbs free energy is evaluated as, Figure 11 illustrate the behavior of the Gibbs free energy at different values of T o and fixed pressure. When T o = 0.112, the curve bounces back and the Gibbs free energy increases along with the decrease in temperature at fixed values of pressure. When P = 0, the Gibbs free energy increases. We conclude that when 0 ≤ P ≤ 0.5, the Gibbs free energy increase for the positive values of temperature.

Conclusion
The radial Einstein's equation at the BH event horizon and cosmological horizon has been used in this paper in which we have obtained the modified first laws of thermodynamics for two BHs: (1) The non-minimal MCRBH and (2) RBH with cosmological constant, for which all the correlated Smarr relations are also satisfied. We treated the cosmological con- stant in two ways, firstly it has been adopted as fixed and then secondly, as the variable which is correlated to the pressure. The surface tension σ can be negative, zero and positive which value is ascertained by using the horizon radius r + , magnetic charge Q m , non-minimal parameter λ, electric charge q and the cosmological horizon . The first law of thermodynamics has been modified in the presence of σ for MCRBH. For > 0, plot of Gibbs free energy versus temperature in the presence of surface tension (for any value) shows that for the increasing values of temperature, the Gibbs free energy approaches to the same value, but for = 0 and negative surface tension, the Gibbs free energy and temperature, both are negative except for σ ≥ 0, the Gibbs free energy increase with increasing temperature. We have obtained the expression for effective Gibbs free energy by modifying the temperature into effective temperature, also analyzed that the effective Gibbs free energy decreases along with increasing effective temperature. We have also discussed pressure at horizon of MCRBH. For fixed values of pressure, negative, positive and zero, behavior of Gibbs free energy has been studied. For > 0, the Gibbs free energy decreases and approaches to the same value with increase in temperature for all fixed values of pressure. For = 0, we have analyzed that, the Gibbs free energy increases when P = 0 and p > 0, after bouncing back at T = 0.016 but for P < 0, it decrease to negative values.
• For RBH with cosmological constant, the first law of thermodynamics has also been modified as Eq. (22) in the presence of σ , which expression is given by Eq. (38). For = 0.01, we have analyzed that the Gibbs free energy increases for all, negative, zero and positive values of σ along with the increase in temperature, after bouncing back at T = 0.16. We have also obtained the expression for Gibbs free energy in the presence of pressure for RBH with cosmological constant. We have analyzed that for = 0.01, the Gibbs free energy decreases and approaches to negative values for negative pressure whereas, for positive pressure, the Gibbs free energy increases.
• We have obtained the modified temperature T o by using the relation between the Bekenstein entropy area relation. The first law of thermodynamics could have taken the possession of the differential form based on the above replacement for which the corresponding Smarr relation is not obeyed. This type of situation was also explained in [61]. When the surface effects of the curved surface are taken into the consideration then the pressures are generally unsteady. The cosmological constant has been treated as the thermodynamic pressure such as P = − 8π which has been accounted for exploring the thermodynamics. To make further investigation for the RBH with cosmological constant , we treated the cosmological constant as the pressure P = −r 2 + 8π(q 2 +r 2 + ) . For RBH with cosmological constant, we have analyzed that, when the is correlated with the pressure, the Gibbs free energy in the presence of effective surface tension increases for lower negative values of σ eff . For positive fixed values of pressure, we have also analyzed that the Gibbs free energy increases for the increasing values of modified temperature, for the RBH with cosmological constant.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .