Impact of Thermal Fluctuations on Logarithmic Corrected Massive Gravity Charged Black Hole

We investigate the influence of the first-order correction of entropy caused by thermal quantum fluctuations on the thermodynamics of a logarithmic corrected charged black hole in massive gravity. For this black hole, we explore the thermodynamic quantities, such as entropy, Helmholtz free energy, internal energy, enthalpy, Gibbs free energy and specific heat. We discuss the influence of the topology of the event horizon, dimensions and nonlinearity parameter on the local and global stability of the black hole. As a result, it is found that the holographic dual parameter vanishes. This means that the thermal corrections have no significant role to disturb the holographic duality of the logarithmic charged black hole in massive gravity, although the thermal corrections have a substantial impact on the thermodynamic quantities in the high-energy limit and the stability conditions of black holes.


Introduction
The thermodynamics of black holes (BHs) has a major role in BH physics and has gained much attention after the detection of Hawking radiations [1]. Bardeen [2] developed four laws of BH thermodynamics and explained the relationship between the thermodynamics of BH and gravity. This relationship was established by making proportionality between the temperature and surface gravity as well as entropy and surface area of the event horizon. In BH thermodynamics, Bekenstein's area law yields entropy, while the first law of thermodynamics provides the temperature of BH [3,4]. After the discovery of these thermodynamical quantities of BHs, the thermodynamics of BH has become an important part for characterizing the physicality and viability of many BH solutions. The corrected thermodynamics of BH has important features, which bring local and global stability or instability, criticality, holographic duality and many other important aspects of BHs. The thermodynamical corrections due to fluctuations have gained a prominent place in BH physics. The corrected BH thermodynamics in the background of matter field were discussed in [5,6]. The influence of thermal fluctuations on charged anti-de Sitter (AdS) BH has shown excellent results of corrections in AdS BHs [7]. The deep thermodynamical analysis of BHs showed that the quantum approach at small scales to the thermodynamics of BH is inevitable, and it leads to corrections in various thermodynamical quantities. The GUP-corrected thermodynamics for all black objects is one of the important approaches [8].
The logarithmic corrections to Godel BH and quantum corrections to the thermodynamics of BH with the Cardy formula were also studied in [9,10], respectively. Nozari et al. [11] provided a detailed study of the influence of corrections on the BH thermodynamics. The quantum corrections to thermodynamics of quasi-topological BHs were studied in [12]. These thermal corrections for Schwarzschild-Beltrami-de Sitter BH [13] and the BH in MG were discussed in literature works [14]. The pioneering work of Frolov et al. [15] studied The paper is structured as follows. In the next section, we present thermal corrections of logarithmic charged BH in MG. In Section 3, we discuss the holographic duality and thermal stability. In Section 4, we summarize our results in concluding remarks.

Thermal Corrections of Logarithmic Charged Black Hole in Massive Gravity
The (n + 1) dimensional action of Einstein MG with negative cosmological constant in context of LNED is given by the following [50]: where R and Λ = −n(n−1) 2l 2 are the Ricci scalar and the cosmological constant, with l being the AdS space time radius, and F = F µν F µν represents the Maxwell invariant. β is a nonlinear constant; for β → ∞, the LNED changes to linear Maxwell. m 2 is a positive parameter of MG, and (m → 0) represents the translation invariance. In Lagrangian (1), c i , U and Γ are the constants, symmetric polynomials and reference metric, respectively. We consider the line element of (n + 1) dimensional space time to obtain the static charged BH solution.
The metric of BH in context of LNED turns out to be as follows [47]: where Υ = 1 + q 2 β 2 r 2−2n , F is hypergeometric function and m 0 is related to the total mass of the BH. q is a constant, which represents the total charge of BH.
In order to find the mass of logarithmic charged BH, setting ψ(r + ) = 0 in Equation (2), one can obtain the following relation: where r + is the horizon radius. The thermodynamics of BH consists of many important thermodynamical quantities; entropy is one of the major thermodynamical quantities, which is defined as follows: The Hawking temperature of BH can be obtained by using T = ≤ 2π = 1 4π ( ∂ ψ(r) ∂r ). The temperature of BH can also be obtained by using following relation [47]: By using the above relation, Hawking temperature of logarithmic charged BH turns out to be as follows: where we use = (c 1 r 3 ) . The usual relation for BH entropy where A is the area of the BH event horizon in the absence of the correction terms. The quantum corrections of gravity near the Planck scale changes the manifold structure of spacetime. This modifies the holographic principle, which leads to change in the entropy-area law of BH. The corrected entropy-area relation is given by the following [51][52][53]: By using Hawking temperature T and the entropy of zeroth order S o in the above relation, the corrected entropy of logarithmic charged BH turns out to be as follows: Figures 1-4 demonstrate the behavior of corrected entropy of logarithmic charged BH in MG vs r + . These plots show the influence of corrected parameter ξ, topology of event horizon k, dimensions n and nonlinearity parameter β on the corrected entropy of BH. There exists a critical point in the corrected entropy of BH ( Figure 1) at r + = 2.6. The behavior of the curves changes before and after the critical point, i.e., the corrected entropy increases for increasing values of the corrected parameter before the critical point, and it decreases after the critical point. For positive values of ξ, entropy increases, which represents stability. For negative ξ, entropy tends to decrease asymptotically, which shows the instability of BH; this behavior changes after the critical point. Figure 2 shows the behavior of corrected entropy for k = 0 flat, k = 1 spherical and k = −1 hyperbolic surfaces. We can observe that the corrected entropy increases for a hyperbolic surface, while it decreases for the case of a flat surface. Figures 3 and 4 show the influence of different dimensions and the nonlinearity parameter on the corrected entropy of logarithmic charged BH in MG. For higher dimensions, the corrected entropy increases, which represents stability; for smaller dimension, the corrected entropy decreases. The impact of dimensions is significant for larger values of r + . The corrected entropy increases with increasing values of the nonlinearity parameter for logarithmic charged BH in MG. For all the plots, as we increase r + , the BH grows, and this causes the area of the event horizon to grow bigger, and hence, the entropy increases. For logarithmic charged BH in MG, we observe that the entropy increases exponentially.    The electrical charge per unit volume of BH is defined as follows [47]: Another important thermodynamical quantity is volume V = ( ∂M ∂P ) S,Q of the BH which it is related to the mass and pressure of the BH. The relation for thermodynamic volume of the logarithmic charged BH becomes the following: Helmholtz free energy (HFE) is the measure of useful work obtained from a closed thermodynamic system. The free energy F can be obtained from F = M − TS, and for the considered BH, it takes the following form: where in the above relation, ∆ = Figures 5-8 demonstrate the corrected HFE behavior of logarithmic charged BH in massive gravity vs. r + . These plots show the influence of corrected parameter ξ, topology of event horizon k, dimensions n and nonlinearity parameter β on the free energy of BH. From Figure 5, one can clearly see that the positive values of ξ increase the HFE, which indicates the instability. In contrast, negative ξ decreases HFE. This plot helps us to identify the region of instability. A negative value of Helmholtz free energy will not be able to extract any useful work from the concerned BH. For ξ = −1, F changes the sign from positive to negative, which represents a phase transition. Figure 6 shows the behavior of HFE for k = 0 flat, k = 1 spherical and k = −1 hyperbolic surfaces. We can observe that there exists discontinuity in the plot for the case of the flat surface. The strange plot of HFE is due to the overcome of quantum effects by tidal forces at very small values of r + . Figures 7 and 8 show the influence of different dimensions and nonlinearity parameters on the HFE of logarithmic charged BH in MG. HFE increases with increasing dimensions, while it decreases with increasing values of the nonlinearity parameter. The HFE decreases w.r.t r + , and there exists discontinuity, which is due to the overcome of quantum effects by tidal forces.  By using ϕ = ∂M the relation for chemical potential of the considered BH turns out to be the following: Moreover, the internal energy for the logarithmic charged BH can be obtained by using E = F + TS. It takes the following form Figures 9-12 demonstrate the internal energy behavior of logarithmic charged BH in MG vs. r + . For the negative value of ξ, the internal energy decreases, and it increases for the positive value. This is in agreement with the first law of BH thermodynamics, as we have already observed that entropy increases for negative value of ξ. Thus, even after incorporating the quantum correction, the first laws of BH thermodynamics hold. The decrease in internal energy of logarithmic charged BH in massive gravity is due to the quantum corrections, which lead to stability of the BHs. For k = 1 spherical and k = −1 hyperbolic surface, the behavior of internal energy of the considered BH is quite different; this different behavior is due to the quantum corrections in the internal energy. The internal energy increases with increasing the dimension and nonlinearity parameter of BH. The modified pressure can be calculated by using P = −( dF dV ); see Appendix A. Figures 13-15 demonstrate the corrected pressure behavior of logarithmic charged BH in MG vs. r + . One can observe that for positive values of k, the pressure becomes negative and it shifts toward a positive value for negative k. This justifies that negative k increases the stability of the BH, which follows the previous cases of entropy, internal energy and free energy. The impact of different dimensions n is quite interesting on the corrected pressure. The pressure remains positive throughout the range of r + for the lower dimension of n, and it becomes negative; there exists discontinuity for higher dimensions. The influence of the nonlinearity parameter is significant only for smaller BHs, and pressure increases with increasing values of the nonlinearity parameter. The enthalpy is another significant thermodynamical quantity, which measures the energy changes of the system. It also helps us to determine the equilibrium conditions of the system. The enthalpy in BH thermodynamics has become more significant after its inclusion in first law of BH thermodynamics. The mass M of an AdS BH works as enthalpy in classical thermodynamics. The enthalpy can be calculated by using H = E + PV; see Appendix B for mathematical relation. Figures 16-19 demonstrate the enthalpy behavior of logarithmic charged BH in MG vs. r + . From the plots, it is clear that enthalpy increases with increasing the value of r + . The negative correction parameter decreases the enthalpy and hence, induces stability; there exist discontinuities in the plots of enthalpy for spherical and hyperbolic surfaces, which is because of the overcome of quantum effects by tidal forces. The enthalpy of logarithmic charged BH in MG increases with increasing dimensions and the nonlinearity parameter.  Figure 19. The enthalpy of logarithmic charged BH in MG. We set k = 1, c 0 = 1, c 1 = 1, n = 3, c 3 = 1, c 2 = 1, and c 4 = 1.
In order to have the static boundary of BH at a fixed temperature, one should have fixed pressure and temperature. In this scenario, the thermodynamical potential to be utilized is the Gibbs free energy. Once we have relations for HFE, pressure and volume, it is simple to calculate the Gibbs free energy by using G = F + PV; see Appendix C.
To study the impact of parameters on the global stability of BH, we plot the relation of Gibbs free energy in Figures 20-23. G > 0 represents the global stability of the BH, while G < 0 shows the region of global instability. The plots show that for negative correction parameter ξ = −1, the Gibbs free remains negative throughout the horizon radius, which represents the global instability of BH. For ξ = 0, 1, smaller BHs remain stable while larger BHs become instable. The Gibbs free energy remains positive for k = 0 flat, k = 1 spherical and k = −1 hyperbolic surfaces, which is a sign of global stability. There exists discontinuity in the plot for spherical surface, which represents the global instability; this unusual behavior is due to the overcome of quantum effects by tidal forces.

Holographic Duality and Thermal Stability
The van der Waals system is one of the most relevant models for discussing the liquidgas system and its critical characteristics. The modification in ideal gas equation provides the equation of state for this model. The van der Waals model is given as follows [54]: where P, T and V are pressure, temperature and specific volume, respectively, B and a represent the size and strength of attraction of the molecules. Here, k = 1 is the Boltzmann constant. One can obtain the ideal gas law by setting a = B = 0. The relationship between the van der Waals and BH helps us to consider the analogy between the temperature of the fluid and the temperature of the BH. The above model can also be expressed as follows: In order to have holographic duality of charged BH in MG, the condition P = P v must hold. By using the temperature and thermodynamic volume of BH, one can obtain the pressure as follows: Figures 24-26 demonstrate the behavior ∆P = P − P v in terms of V. For large BHs, ∆P → 0, which shows that the thermal correction has no significant role to disturb the holographic duality of the logarithmic charged BH in MG. The influence of the topology of the event horizon k and nonlinearity parameter is opposite.  The stability of BH can be considered in both canonical and grand canonical ensembles. In canonical ensemble, C > 0 along with T > 0 represents the thermal stability, while C < 0 represents the unstable region, and the phase transition occurs in an unstable region to attain stability. The specific heat is defined as follows [47]: Now, by plugging the relations of entropy and temperature of the logarithmic charged BH in MG in the above equation, we can obtain the specific heat (see Appendix D).
For unstable BH, phase transition occurs to gain stability. The roots and points of divergences of C provides the phase transition points. The roots of the heat capacity help us to determine the thermal transitions between the physical and un-physical states of BH. For the considered BH, as the roots and divergences of the heat capacity are not possible analytically, we have shown the behavior specific heat of logarithmic charged BH in MG in Figures 27-30. One can see that in the plot of Figure 27, the divergence in the heat capacity occurs at r + = 1.4. The heat capacity remains positive before and after the point of divergence for ξ = 0, −1, which represent the stable phase. The phase transition occurs at the point of divergence, and the negative range of C shows the region of instability of BH. There exist discontinuities or divergences for k = 0 flat, k = 1 spherical and k = −1 hyperbolic surfaces for smaller BHs-which represents instability, while for larger values of r + , BHs becomes stable for flat and hyperbolic surfaces. The nonlinearity parameter β shows the significant role on the specific heat of the considered BH. For the small BH, the positive heat capacity shows the local stability, while for large BH, C < 0 and discontinuity in heat capacity represent the local instability. One can see that for logarithmic charged BH in higher dimensions, the instability of BH increases.    The determinant of the Hessian matrix H helps us to discuss the local stability. As the charge remains as a fixed parameter in the canonical ensemble, the positive region of heat capacity C is sufficient for the stability of BH. However, we need to examine the sign of H to examine the thermal stability in the grand canonical ensemble. The Hessian matrix is defined as follows [47]: From the above equation, by setting H F X i ,X j = 0, one can easily obtain the following relation: where τ 1 = ∂ 2 F ∂T 2 and τ 2 = ∂ 2 F ∂≡ 2 . The important condition for the stability of BH is τ ≥ 0. Figure 31 demonstrates the behavior of trace of Hessian matrix τ verses r + for different values of ξ. The necessary condition for the stability is τ ≥ 0; one can observe this for all the values of the corrected parameters ξ, τ ≤ 0, which show the local instability of logarithmic charged BH in the MG.

Concluding Remarks
In this work, we considered the logarithmic charged BH in MG with the negative cosmological constant and studied the corrected thermodynamics. We computed the relations of logarithmic corrected S, F, E, H, G and C. In order to investigate the local and global stability and phase transition points, we plotted these relations against horizon radius r + . We provided the deep analysis of the impact of corrected parameter ξ, topology of event horizon k, dimensions n and nonlinearity parameter β on the thermodynamical quantities of the BH.
We found that there exists a critical point in the corrected entropy of BH ( Figure 6) at r + = 2.6, and the behavior of the curve changes before and after the critical point. The corrected entropy increases for increasing values of corrected parameter ξ before the critical point, and it decreases after the critical point. For positive values of ξ, the entropy increases, which represents the stability, while for negative ξ, the entropy tends to decrease asymptotically, which shows the instability of BH. This behavior changes after the critical point. The HFE helped us to identify the region of instability for ξ = −1, as F changes the sign from positive to negative, which represents a phase transition. For negative values of ξ, the internal energy decreases, while it increases for positive values, which is in agreement with the first law of BH thermodynamics. The decrease in internal energy of logarithmic charged BH in MG is due to the quantum corrections, which lead to stability of the BHs. We also observed that for the negative correction parameter ξ = −1, the Gibbs free remains negative throughout the horizon radius, which represents the global instability of BH. For ξ = 0, 1 smaller BHs remain stable, while larger BHs become unstable. The Gibbs free energy remains positive for k = 0 flat, k = 1 spherical and k = −1 hyperbolic surfaces, which is sign of global stability. We analyzed that in the plot of Figure 27 , the divergence in the heat capacity occurred at r + = 1.4. The heat capacity remained positive before and after the point of divergence for ξ = 0, −1, which showed the stable phase. Phase transition occurred at the point of divergence, and the negative range of C showed the region of instability of BH. For large BHs, ∆P → 0, which showed that thermal correction has no significant role in disturbing the holographic duality of the logarithmic charged BH in MG.  Acknowledgments: A.J. and S.C. are thankful to HEC Pakistan for the Indigenous Fellowship phase II.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
For the logarithmic charged BH in MG, the modified pressure becomes: