\UseRawInputEncoding
thanks: e-mail:[email protected], [email protected]

Barrow’s non-linear charged Anti-de Sitter black hole and stability

Yun-Zhi Du1,2111the corresponding author and Qiang Gu1 Shanxi Datong University, Datong 037009, China
State Key Laboratory of Quantum Optics and Quantum Optics Devices, Shanxi University, Taiyuan, Shanxi 030006, China
Abstract

As we know that the horizon area of a black hole will increase when it absorbs matters. While based on Barrow’s idea of fractal black hole horizon, ones [Phys. Lett. B 831 (137181) 2022] had proposed that for a spherically fractal structure the minimal increase of the horizon area is the area of the smallest bubble sphere. And the corresponding black hole entropy is of a logarithmic form, which is similar to that of Boltzmann entropy under a certain condition. Based on these, we re-derive the entropy of the Barrow’s Einstein-power-Yang-Mills (EPYM) AdS black hole, and calculate the temperature and heat capacity of the Barrow’s EPYM AdS black hole. There exists an interesting phenomena that the ratio between the Barrow’s temperature and the Hawking temperature of the EPYM AdS black hole is fully consistent with that of other Schwarzschild-like black holes. The Barrow’s temperature and Hawking temperature with the certain range of ΛΛ\Lambdaroman_Λ are monotonically increasing and the corresponding heat capacities are all positive, which means these black holes are thermodynamically stable. Besides, for the Barrow’s EPYM AdS black hole its heat capacity has a Schottky anomaly-like behavior, which may reflect the existence of the discrete energy level and the microscopical degree of freedom.

I Introduction

As we known that there arises the non-renormalizability which is a well-known puzzle when gravity is combined with quantum theory. In 1916191619161916 Einstein published the general relativity theory in which a concept about acceleration differences between non-inertial frames was proposed. It sheds light on the idea that the space and the time are both distinct features, and the spacetime is curved. However, the spacetime curvature depends on the presence of gravity, which is also capable of wrapping light, radio waves and many other things. In 1939193919391939, Oppenheimer and Snyder published a paper titled as “On Continuing Gravitational Contraction” on black holes Oppenheimer1939 . They predicted the continued contraction of a star, forming a body with an extreme attraction force that could not escape even light from it. However, Einstein established the general relativity and predicted the presence of gravitational waves, although he did not have any experimental evidence to prove the existence of gravitational waves. Over nearly half of a century, when LIGO physically detected the undulations in spacetime induced by gravitational waves, the so-called “GW150914” phenomenon, is a proof of Einstein’s theory Abbatial2016 . This is the first clear observation of a pair of black holes that merge to form a single black hole (BH). Black holes are possibly the most ideal thermal structures that have played an increasingly important role in the general relativity interdisciplinary fields, the information theory, the quantum mechanics, and the statistical physics.

The four laws of BH mechanics formalized its behavior and the microscopic properties, closely implies the thermal properties of BHs Bardeen1973 ; Wald2001 . After the work of Bekenstein, Stephen Hawking was the first to demonstrate the presence of BH radiations, which is defined as a tunneling process caused by vacuum fluctuations near the horizon. A small amount of heat called Hawking Radiation recorded by Hawking Hawking1975 is radiated due to irregularities of quantum mechanics. BHs were found to possess thermodynamic properties such as entropy and temperature in the Pioneering work Bekenstein1972 ; Bekenstein1973 . These findings show that a BH is a thermodynamic system with a Hawking temperature proportional to its surface gravity on the horizon area and obeys the four BH thermodynamic laws. AdS BH thermodynamics is the most impressive instrument for the study of quantum gravity in spacetime containing horizon Birrell1982 . The quantum scenario for Hawking radiation provides a valuable insight that a black hole temperature is proportional to its horizon area and that its entropy is related to its area gravity Abreu2020 . This analysis revealed that there is a strong relationship between thermodynamics and gravity. The Barrow entropy Barrow2020 arises from the notion that the surface of a BH can be altered as a result of quantum gravitational effects. Recently Barrow has shown that the area of a BH can be represented as a discontinuous and possibly fractal structure of the geometry of the horizon Barrow2020 . Starting from the idea of “Koch snowflake” Barrow proposed a model in which the event horizon of a Schwarzschild black hole radius is surrounded by N𝑁Nitalic_N smaller spheres whose radii have a ratio λ𝜆\lambdaitalic_λ with respect to the original sphere. After n𝑛nitalic_n iterations, the radius is rn+1=λrnsubscript𝑟𝑛1𝜆subscript𝑟𝑛r_{n+1}=\lambda r_{n}italic_r start_POSTSUBSCRIPT italic_n + 1 end_POSTSUBSCRIPT = italic_λ italic_r start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT. By adding smaller spheres to the surface, the total area and volume after an infinite number of steps should be the sum of all the intricate structures, which can be finite and infinite. The corresponding entropy can be very large. Subsequently, the Barrow’s idea has been further extended to the investigation of the dark energy Saridakis2020 ; Moradpour2020 , the cosmology Salehi2023 ; Komatsu2024 ; Okcu2024 and black hole thermodynamics Abreu2020a ; Ma2022 ; Ladghami2024 ; Rani2023 .

In Ref. Ma2022 the authors proposed a new way to investigate the effect of a minimal length on the entropy and the stability of black holes with fractal structures by considering the generalized uncertainty principle. Based on this idea, we apply the corresponding Barrow’s entropy formula to the Einstein-power-Yang-Mills (EPYM) AdS black hole, and analyze the stability of Barrow’s black holes. In Sec. II, we briefly review the thermodynamic quantities and stability of the standard EPYM AdS black hole. In Sec. III, we review Barrow’s fractal black hole and the Barrow’s entropy formula with the consideration of the generalized uncertainty principle. Then based on the fractal structure we calculate the temperature and heat capacity of the AdS black hole with the non-linear source. A brief summary is given in Sec. IV.

II Non-linear charged AdS black hole Solution

The action for four-dimensional Einstein-power-Yang-Mills (EPYM) gravity with a cosmological constant ΛΛ\Lambdaroman_Λ was given by Zhang2015 ; Corda2011 ; Mazharimousavi2009 ; Lorenci2002 (8πG=18𝜋𝐺18\pi G=18 italic_π italic_G = 1)

I=12d4xg(R2Λγ)𝐼12superscript𝑑4𝑥𝑔𝑅2Λsuperscript𝛾\displaystyle I=\frac{1}{2}\int d^{4}x\sqrt{g}\left(R-2\Lambda-\mathcal{F}^{% \gamma}\right)italic_I = divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ italic_d start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x square-root start_ARG italic_g end_ARG ( italic_R - 2 roman_Λ - caligraphic_F start_POSTSUPERSCRIPT italic_γ end_POSTSUPERSCRIPT ) (1)

with the Yang-Mills (YM) invariant \mathcal{F}caligraphic_F and the YM field Fμν(a)superscriptsubscript𝐹𝜇𝜈𝑎F_{\mu\nu}^{(a)}italic_F start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT

\displaystyle\mathcal{F}caligraphic_F =\displaystyle== Tr(Fμν(a)F(a)μν),Trsubscriptsuperscript𝐹𝑎𝜇𝜈superscript𝐹𝑎𝜇𝜈\displaystyle\operatorname{Tr}(F^{(a)}_{{\mu\nu}}F^{{(a)\mu\nu}}),roman_Tr ( italic_F start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT italic_F start_POSTSUPERSCRIPT ( italic_a ) italic_μ italic_ν end_POSTSUPERSCRIPT ) , (2)
Fμν(a)superscriptsubscript𝐹𝜇𝜈𝑎\displaystyle F_{\mu\nu}^{(a)}italic_F start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT =\displaystyle== μAν(a)νAμ(a)+12ξC(b)(c)(a)Aμ(b)Aν(c).subscript𝜇superscriptsubscript𝐴𝜈𝑎subscript𝜈superscriptsubscript𝐴𝜇𝑎12𝜉superscriptsubscript𝐶𝑏𝑐𝑎superscriptsubscript𝐴𝜇𝑏superscriptsubscript𝐴𝜈𝑐\displaystyle\partial_{\mu}A_{\nu}^{(a)}-\partial_{\nu}A_{\mu}^{(a)}+\frac{1}{% 2\xi}C_{(b)(c)}^{(a)}A_{\mu}^{(b)}A_{\nu}^{(c)}.∂ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT - ∂ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 italic_ξ end_ARG italic_C start_POSTSUBSCRIPT ( italic_b ) ( italic_c ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT italic_A start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_b ) end_POSTSUPERSCRIPT italic_A start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_c ) end_POSTSUPERSCRIPT . (3)

Here, Tr(Fμν(a)F(a)μν)=a=13Fμν(a)F(a)μνTrsubscriptsuperscript𝐹𝑎𝜇𝜈superscript𝐹𝑎𝜇𝜈subscriptsuperscript3𝑎1subscriptsuperscript𝐹𝑎𝜇𝜈superscript𝐹𝑎𝜇𝜈\operatorname{Tr}(F^{(a)}_{\mu\nu}F^{(a)\mu\nu})=\sum^{3}_{a=1}F^{(a)}_{\mu\nu% }F^{(a)\mu\nu}roman_Tr ( italic_F start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT italic_F start_POSTSUPERSCRIPT ( italic_a ) italic_μ italic_ν end_POSTSUPERSCRIPT ) = ∑ start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a = 1 end_POSTSUBSCRIPT italic_F start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT italic_F start_POSTSUPERSCRIPT ( italic_a ) italic_μ italic_ν end_POSTSUPERSCRIPT, R𝑅Ritalic_R and γ𝛾\gammaitalic_γ are the scalar curvature and a positive real parameter, respectively; C(b)(c)(a)superscriptsubscript𝐶𝑏𝑐𝑎C_{(b)(c)}^{(a)}italic_C start_POSTSUBSCRIPT ( italic_b ) ( italic_c ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT represents the structure constants of three-parameter Lie group G𝐺Gitalic_G; ξ𝜉\xiitalic_ξ is the coupling constant; and Aμ(a)superscriptsubscript𝐴𝜇𝑎A_{\mu}^{(a)}italic_A start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_a ) end_POSTSUPERSCRIPT represents the SO(3)𝑆𝑂3SO(3)italic_S italic_O ( 3 ) gauge group Yang-Mills (YM) potentials defining by the Wu-Yang (WY) ansatz Balakin2016 ; Balakin2007 ; Balakin2007a . For this system, the EPYM black hole solution with the negative cosmological constant ΛΛ\Lambdaroman_Λ in the four-dimensional spacetime is obtained by adopting the following metric Yerra2019 :

ds2=f(r)dt2+f1dr2+r2dΩ22,𝑑superscript𝑠2𝑓𝑟𝑑superscript𝑡2superscript𝑓1𝑑superscript𝑟2superscript𝑟2𝑑superscriptsubscriptΩ22\displaystyle ds^{2}=-f(r)dt^{2}+f^{-1}dr^{2}+r^{2}d\Omega_{2}^{2},italic_d italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = - italic_f ( italic_r ) italic_d italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_f start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_d italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d roman_Ω start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (4)
f(r)=12MrΛ3r2+(2q2)γ2(4γ3)r4γ2,𝑓𝑟12𝑀𝑟Λ3superscript𝑟2superscript2superscript𝑞2𝛾24𝛾3superscript𝑟4𝛾2\displaystyle f(r)=1-\frac{2M}{r}-\frac{\Lambda}{3}r^{2}+\frac{\left(2q^{2}% \right)^{\gamma}}{2(4\gamma-3)r^{4\gamma-2}},italic_f ( italic_r ) = 1 - divide start_ARG 2 italic_M end_ARG start_ARG italic_r end_ARG - divide start_ARG roman_Λ end_ARG start_ARG 3 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG ( 2 italic_q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_γ end_POSTSUPERSCRIPT end_ARG start_ARG 2 ( 4 italic_γ - 3 ) italic_r start_POSTSUPERSCRIPT 4 italic_γ - 2 end_POSTSUPERSCRIPT end_ARG , (5)

where dΩ22𝑑superscriptsubscriptΩ22d\Omega_{2}^{2}italic_d roman_Ω start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT is the metric on unit 2222-sphere with volume 4π4𝜋4\pi4 italic_π and q𝑞qitalic_q is the YM charge. Note that this solution is valid for the condition of the non-linear YM charge parameter γ0.75𝛾0.75\gamma\neq 0.75italic_γ ≠ 0.75, and the power YM term holds the weak energy condition (WEC) for γ>0𝛾0\gamma>0italic_γ > 0 Corda2011 . In the extended phase space, ΛΛ\Lambdaroman_Λ was interpreted as the thermodynamic pressure P=Λ8π𝑃Λ8𝜋P=-\frac{\Lambda}{8\pi}italic_P = - divide start_ARG roman_Λ end_ARG start_ARG 8 italic_π end_ARG. For γ=1𝛾1\gamma=1italic_γ = 1, the YM charge term in the metric has the same form as the Maxwell charge term for the Einstein-Maxwell-Yang-Mills (EMYM) theory. Thus the contribution of the YM charge term on the thermodynamic property should be same as that of the Maxwell charge term in EMYM theory. The only difference is that they have different gauge groups: the YM field is of SO(3)𝑆𝑂3SO(3)italic_S italic_O ( 3 ), while the EM field is of U(1)𝑈1U(1)italic_U ( 1 ). These comments are consistent with that in Ref. Moumni2018 . The black hole event horizon locates at f(r+)=0𝑓subscript𝑟0f(r_{+})=0italic_f ( italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) = 0. The parameter M𝑀Mitalic_M represents the ADM mass of the black hole and it reads

M(r+,q,γ)=r+2[1Λr+23+(2q2)γ2(4γ3)r+4γ2].𝑀subscript𝑟𝑞𝛾subscript𝑟2delimited-[]1Λsuperscriptsubscript𝑟23superscript2superscript𝑞2𝛾24𝛾3superscriptsubscript𝑟4𝛾2\displaystyle M(r_{+},q,\gamma)=\frac{r_{+}}{2}\left[1-\frac{\Lambda r_{+}^{2}% }{3}+\frac{\left(2q^{2}\right)^{\gamma}}{2(4\gamma-3)r_{+}^{4\gamma-2}}\right].italic_M ( italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT , italic_q , italic_γ ) = divide start_ARG italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG [ 1 - divide start_ARG roman_Λ italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 3 end_ARG + divide start_ARG ( 2 italic_q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_γ end_POSTSUPERSCRIPT end_ARG start_ARG 2 ( 4 italic_γ - 3 ) italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 italic_γ - 2 end_POSTSUPERSCRIPT end_ARG ] . (6)

And in our set up it is associated with the enthalpy of the system. The black hole temperature and Hawking entropy were given by Zhang2015

T=14πr+(1Λr+2(2q2)γ2r+4γ2),S=πr+2.formulae-sequence𝑇14𝜋subscript𝑟1Λsuperscriptsubscript𝑟2superscript2superscript𝑞2𝛾2superscriptsubscript𝑟4𝛾2𝑆𝜋superscriptsubscript𝑟2\displaystyle T=\frac{1}{4\pi r_{+}}\left(1-\Lambda r_{+}^{2}-\frac{\left(2q^{% 2}\right)^{\gamma}}{2r_{+}^{4\gamma-2}}\right),~{}~{}~{}~{}~{}S=\pi r_{+}^{2}.italic_T = divide start_ARG 1 end_ARG start_ARG 4 italic_π italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT end_ARG ( 1 - roman_Λ italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG ( 2 italic_q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_γ end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 italic_γ - 2 end_POSTSUPERSCRIPT end_ARG ) , italic_S = italic_π italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (7)

In order to probe the stability of the EPYM AdS black hole, we present the heat capacity at the constant charge q𝑞qitalic_q, which is a counterpart of the heat capacity at the constant volume, as the following form

C=TST=2πr+2(1Λr+2(2q2)γ2r+4γ2)1+Λr+2(2γ1)(2q2)γ2r+4γ2.𝐶𝑇𝑆𝑇2𝜋superscriptsubscript𝑟21Λsuperscriptsubscript𝑟2superscript2superscript𝑞2𝛾2superscriptsubscript𝑟4𝛾21Λsuperscriptsubscript𝑟22𝛾1superscript2superscript𝑞2𝛾2superscriptsubscript𝑟4𝛾2\displaystyle C=T\frac{\partial S}{\partial T}=-\frac{2\pi r_{+}^{2}\left(1-% \Lambda r_{+}^{2}-\frac{\left(2q^{2}\right)^{\gamma}}{2r_{+}^{4\gamma-2}}% \right)}{1+\Lambda r_{+}^{2}-\frac{(2\gamma-1)\left(2q^{2}\right)^{\gamma}}{2r% _{+}^{4\gamma-2}}}.italic_C = italic_T divide start_ARG ∂ italic_S end_ARG start_ARG ∂ italic_T end_ARG = - divide start_ARG 2 italic_π italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - roman_Λ italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG ( 2 italic_q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_γ end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 italic_γ - 2 end_POSTSUPERSCRIPT end_ARG ) end_ARG start_ARG 1 + roman_Λ italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG ( 2 italic_γ - 1 ) ( 2 italic_q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_γ end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 italic_γ - 2 end_POSTSUPERSCRIPT end_ARG end_ARG . (8)

The behaviours of the temperature and the heat capacity for the EPYM AdS black hole are shown in Fig. 1. The results indicate that there exist a critical value of the cosmological constant Λc=γ1/2(2q2γ2)γ2γ1(4γ1)12γ1subscriptΛ𝑐𝛾12superscript2superscript𝑞2superscript𝛾2𝛾2𝛾1superscript4𝛾112𝛾1\Lambda_{c}=-\frac{\gamma-1/2}{(2q^{2}\gamma^{2})^{\frac{\gamma}{2\gamma-1}}(4% \gamma-1)^{\frac{1}{2\gamma-1}}}roman_Λ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = - divide start_ARG italic_γ - 1 / 2 end_ARG start_ARG ( 2 italic_q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG italic_γ end_ARG start_ARG 2 italic_γ - 1 end_ARG end_POSTSUPERSCRIPT ( 4 italic_γ - 1 ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 italic_γ - 1 end_ARG end_POSTSUPERSCRIPT end_ARG, when ΛΛcΛsubscriptΛ𝑐\Lambda\geq\Lambda_{c}roman_Λ ≥ roman_Λ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT the temperature is monotonically increasing with r+subscript𝑟r_{+}italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT, the heat capacity that has a Schottky Anomaly-like behaviour and it is always positive. Below the critical value, the temperature have some extremes, and the corresponding heat capacity is positive only in the interval where the slope of the temperature is positive. This configuration corresponds to the first-order phase transition between the high-potential black hole and the low-potential one as shown in Ref. Du2021 . In the next part, we will investigate the temperature and heat capacity of the Barrow’s EPYM AdS black hole.

Refer to caption
Figure 1: Plots of T,Cr+𝑇𝐶subscript𝑟T,C-r_{+}italic_T , italic_C - italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT with the parameters q=0.85𝑞0.85q=0.85italic_q = 0.85, γ=0.8𝛾0.8\gamma=0.8italic_γ = 0.8, and Λ=0.04474Λ0.04474\Lambda=-0.04474roman_Λ = - 0.04474 (the thicked solid lines), Λ=0.08947Λ0.08947\Lambda=-0.08947roman_Λ = - 0.08947 (the thin solid lines), Λ=0.1342Λ0.1342\Lambda=-0.1342roman_Λ = - 0.1342 (the dashed lines).

III Barrow’s EPYM AdS entropy and stability

From Barrow’s idea, the fractal structure of Schwarzschild black hole horizon is that there are many smaller spheres that can attach on the black hole horizon, then more smaller spheres attaching to these spheres and so on. At each step there are N𝑁Nitalic_N spheres and the horizon radius is λ𝜆\lambdaitalic_λ times smaller than that of the sphere in the previous step. For no cut off at some small finite scale, the corresponding area of the horizon and volume of the black hole were given in Ref.

A=4πrh2n=0(Nλ2)n,V=4πrh33n=0(Nλ3)n,formulae-sequencesubscript𝐴4𝜋subscriptsuperscript𝑟2subscriptsuperscript𝑛0superscript𝑁superscript𝜆2𝑛subscript𝑉4𝜋subscriptsuperscript𝑟33subscriptsuperscript𝑛0superscript𝑁superscript𝜆3𝑛\displaystyle A_{\infty}=4\pi r^{2}_{h}\sum^{\infty}_{n=0}\left(N\lambda^{2}% \right)^{n},~{}~{}V_{\infty}=\frac{4\pi r^{3}_{h}}{3}\sum^{\infty}_{n=0}\left(% N\lambda^{3}\right)^{n},italic_A start_POSTSUBSCRIPT ∞ end_POSTSUBSCRIPT = 4 italic_π italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT ∑ start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n = 0 end_POSTSUBSCRIPT ( italic_N italic_λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT , italic_V start_POSTSUBSCRIPT ∞ end_POSTSUBSCRIPT = divide start_ARG 4 italic_π italic_r start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT end_ARG start_ARG 3 end_ARG ∑ start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n = 0 end_POSTSUBSCRIPT ( italic_N italic_λ start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT , (9)

where rhsubscript𝑟r_{h}italic_r start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT is the black hole horizon radius and λ𝜆\lambdaitalic_λ is the ratio of radii with respect to the original sphere. As λ2<N<λ3superscript𝜆2𝑁superscript𝜆3\lambda^{-2}<N<\lambda^{-3}italic_λ start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT < italic_N < italic_λ start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT the black hole horizon area will be infinite, while the volume is finite. On the other hand, Barrow proposed that the entropy can take the form as S(Ah/Apl)1+/2𝑆superscriptsubscript𝐴subscript𝐴𝑝𝑙12S\approx(A_{h}/A_{pl})^{1+\triangle/2}italic_S ≈ ( italic_A start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT / italic_A start_POSTSUBSCRIPT italic_p italic_l end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 1 + △ / 2 end_POSTSUPERSCRIPT with 0<<1010<\triangle<10 < △ < 1. Here Ahsubscript𝐴A_{h}italic_A start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT stands for the standard horizon area and Aplsubscript𝐴𝑝𝑙A_{pl}italic_A start_POSTSUBSCRIPT italic_p italic_l end_POSTSUBSCRIPT is the Planck area. As =00\triangle=0△ = 0 it reduces to the standard black hole entropy, and it corresponds to one of the most intricate horizon when =11\triangle=1△ = 1.

In the frame of the generalized uncertainty principle, the lowest limit of the length scale is the Planck length lpsubscript𝑙𝑝l_{p}italic_l start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT. The radius of the smallest sphere in the fractal structure can not exactly be lpsubscript𝑙𝑝l_{p}italic_l start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT. From the recurrence relation rn+1=λrhsubscript𝑟𝑛1𝜆subscript𝑟r_{n+1}=\lambda r_{h}italic_r start_POSTSUBSCRIPT italic_n + 1 end_POSTSUBSCRIPT = italic_λ italic_r start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT, at some cut-off step n1subscript𝑛1n_{1}italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT it only needs to satisfy rn1lp,λrn1<lpformulae-sequencesubscript𝑟𝑛1subscript𝑙𝑝𝜆subscript𝑟𝑛1subscript𝑙𝑝r_{n1}\geq l_{p},~{}\lambda r_{n1}<l_{p}italic_r start_POSTSUBSCRIPT italic_n 1 end_POSTSUBSCRIPT ≥ italic_l start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT , italic_λ italic_r start_POSTSUBSCRIPT italic_n 1 end_POSTSUBSCRIPT < italic_l start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT. In this case, the Barrow’s black hole entropy was given in Ref. Ma2022

SB=ln2c1ln(Ah/Apl)+c0,subscript𝑆𝐵𝑙𝑛2subscript𝑐1𝑙𝑛subscript𝐴subscript𝐴𝑝𝑙subscript𝑐0\displaystyle S_{B}=\frac{ln2}{c_{1}}ln(A_{h}/A_{pl})+c_{0},italic_S start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT = divide start_ARG italic_l italic_n 2 end_ARG start_ARG italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG italic_l italic_n ( italic_A start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT / italic_A start_POSTSUBSCRIPT italic_p italic_l end_POSTSUBSCRIPT ) + italic_c start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , (10)

where c1=λ2n1(1Nλ2)1(Nλ2)n1+1subscript𝑐1superscript𝜆2subscript𝑛11𝑁superscript𝜆21superscript𝑁superscript𝜆2subscript𝑛11c_{1}=\frac{\lambda^{2n_{1}}(1-N\lambda^{2})}{1-(N\lambda^{2})^{n_{1}+1}}italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = divide start_ARG italic_λ start_POSTSUPERSCRIPT 2 italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ( 1 - italic_N italic_λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 1 - ( italic_N italic_λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 1 end_POSTSUPERSCRIPT end_ARG is a positive constant and c0subscript𝑐0c_{0}italic_c start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is a integrate constant. Although the Barrow entropy does not satisfy the standard area law, or the logarithmic correction to area law, it seams to be the well-known Boltzmann entropy S=kBlnΩ𝑆subscript𝑘𝐵𝑙𝑛ΩS=k_{B}ln\Omegaitalic_S = italic_k start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT italic_l italic_n roman_Ω when Ah/AplΩsubscript𝐴subscript𝐴𝑝𝑙ΩA_{h}/A_{pl}\leftrightarrow\Omegaitalic_A start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT / italic_A start_POSTSUBSCRIPT italic_p italic_l end_POSTSUBSCRIPT ↔ roman_Ω and ln2c1=1𝑙𝑛2subscript𝑐11\frac{ln2}{c_{1}}=1divide start_ARG italic_l italic_n 2 end_ARG start_ARG italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG = 1, c0=0subscript𝑐00c_{0}=0italic_c start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0. And the corresponding parameter Aplsubscript𝐴𝑝𝑙A_{pl}italic_A start_POSTSUBSCRIPT italic_p italic_l end_POSTSUBSCRIPT stands for the area occupied by one microscopic state, Ah/Aplsubscript𝐴subscript𝐴𝑝𝑙A_{h}/A_{pl}italic_A start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT / italic_A start_POSTSUBSCRIPT italic_p italic_l end_POSTSUBSCRIPT is just the number of microscopic states of the black hole.

Since as a thermodynamic system, the thermodynamic quantities of black holes should satisfy the thermodynamical law, dM=TBdSB+𝑑𝑀subscript𝑇𝐵𝑑subscript𝑆𝐵dM=T_{B}dS_{B}+...italic_d italic_M = italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT italic_d italic_S start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT + …. Hence the Barrow temperature of the EPYM AdS black hole reads

1TB=SBM=4ln2c1r+(1Λr+2(2q2)γ/2rh4γ2).1subscript𝑇𝐵subscript𝑆𝐵𝑀4𝑙𝑛2subscript𝑐1subscript𝑟1Λsuperscriptsubscript𝑟2superscript2superscript𝑞2𝛾2superscriptsubscript𝑟4𝛾2\displaystyle\frac{1}{T_{B}}=\frac{\partial S_{B}}{\partial M}=\frac{4ln2}{c_{% 1}r_{+}\left(1-\Lambda r_{+}^{2}-(2q^{2})^{\gamma}/2r_{h}^{4\gamma-2}\right)}.divide start_ARG 1 end_ARG start_ARG italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT end_ARG = divide start_ARG ∂ italic_S start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_M end_ARG = divide start_ARG 4 italic_l italic_n 2 end_ARG start_ARG italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ( 1 - roman_Λ italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( 2 italic_q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_γ end_POSTSUPERSCRIPT / 2 italic_r start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 italic_γ - 2 end_POSTSUPERSCRIPT ) end_ARG . (11)

When Λ=0Λ0\Lambda=0roman_Λ = 0 and q=0𝑞0q=0italic_q = 0, this temperature reduces to the one of the Barrow’s Schwarzschild black hole; it also can reduce to the temperature of the Barrow’s RN black hole with the parameters Λ=0Λ0\Lambda=0roman_Λ = 0, γ=1𝛾1\gamma=1italic_γ = 1. Furthermore the temperature of the Barrow’s RN-AdS black hole can also be recovered with the parameter γ=1𝛾1\gamma=1italic_γ = 1. There exists an interesting phenomena, the ratio between two entropies τTB/T=c1πr+2/ln2𝜏subscript𝑇𝐵𝑇subscript𝑐1𝜋superscriptsubscript𝑟2𝑙𝑛2\tau\equiv T_{B}/T=c_{1}\pi r_{+}^{2}/ln2italic_τ ≡ italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT / italic_T = italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_π italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_l italic_n 2 is fully consistent with that of other Barrow’s black holes. It’s obviously that the Barrow’s temperature and Hawking temperature both become zero in the extremal limit, Λ=0Λ0\Lambda=0roman_Λ = 0 and r+=2(γ1)/(4γ2)qγ/(2γ1)subscript𝑟superscript2𝛾14𝛾2superscript𝑞𝛾2𝛾1r_{+}=2^{(\gamma-1)/(4\gamma-2)}q^{\gamma/(2\gamma-1)}italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT = 2 start_POSTSUPERSCRIPT ( italic_γ - 1 ) / ( 4 italic_γ - 2 ) end_POSTSUPERSCRIPT italic_q start_POSTSUPERSCRIPT italic_γ / ( 2 italic_γ - 1 ) end_POSTSUPERSCRIPT.

In order to investigate the stability of the Barrow’s EPYM AdS black hole, the corresponding heat capacity at the constant charge can be calculated as the following

CB=MTB=2ln2(1Λr+2(2q2)γ2r+(4γ2))c1(13Λr+2(4γ3)(2q2)γ2r+4γ2).subscript𝐶𝐵𝑀subscript𝑇𝐵2𝑙𝑛21Λsuperscriptsubscript𝑟2superscript2superscript𝑞2𝛾2superscriptsubscript𝑟4𝛾2subscript𝑐113Λsuperscriptsubscript𝑟24𝛾3superscript2superscript𝑞2𝛾2superscriptsubscript𝑟4𝛾2\displaystyle C_{B}=\frac{\partial M}{\partial T_{B}}=\frac{2ln2\left(1-% \Lambda r_{+}^{2}-\frac{\left(2q^{2}\right)^{\gamma}}{2r_{+}^{(4\gamma-2)}}% \right)}{c_{1}\left(1-3\Lambda r_{+}^{2}-\frac{(4\gamma-3)\left(2q^{2}\right)^% {\gamma}}{2r_{+}^{4\gamma-2}}\right)}.italic_C start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT = divide start_ARG ∂ italic_M end_ARG start_ARG ∂ italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT end_ARG = divide start_ARG 2 italic_l italic_n 2 ( 1 - roman_Λ italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG ( 2 italic_q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_γ end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 4 italic_γ - 2 ) end_POSTSUPERSCRIPT end_ARG ) end_ARG start_ARG italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( 1 - 3 roman_Λ italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG ( 4 italic_γ - 3 ) ( 2 italic_q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_γ end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 italic_γ - 2 end_POSTSUPERSCRIPT end_ARG ) end_ARG . (12)

The corresponding behaviours of two temperatures and the heat capacity as the functions of horizon radius are shown in Fig. 2(a). For the EPYM AdS black hole with the any given value of ΛΛ\Lambdaroman_Λ, the Barrow’s temperature TBsubscript𝑇𝐵T_{B}italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT is monotonically increasing with the increasing of r+subscript𝑟r_{+}italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT, as well as the Hawking temperature T𝑇Titalic_T with ΛΛcΛsubscriptΛ𝑐\Lambda\geq\Lambda_{c}roman_Λ ≥ roman_Λ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT, while as Λ<ΛcΛsubscriptΛ𝑐\Lambda<\Lambda_{c}roman_Λ < roman_Λ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT the Hawking temperature T𝑇Titalic_T first increase to a maximum and then decreases monotonically. With the increasing of r+subscript𝑟r_{+}italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT, the heat capacity CBsubscript𝐶𝐵C_{B}italic_C start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT is positive and it also has the behaviour of the Schottky anomaly as like the standard black hole system with ΛΛcΛsubscriptΛ𝑐\Lambda\geq\Lambda_{c}roman_Λ ≥ roman_Λ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT. And CBsubscript𝐶𝐵C_{B}italic_C start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT will approach to zero when TB0subscript𝑇𝐵0T_{B}\rightarrow 0italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT → 0, in the limitation of TBsubscript𝑇𝐵T_{B}\rightarrow\inftyitalic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT → ∞ the Barrow’s heat capacity CBsubscript𝐶𝐵C_{B}italic_C start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT tends to the constant 2ln2/(3c1)2𝑙𝑛23subscript𝑐12ln2/(3c_{1})2 italic_l italic_n 2 / ( 3 italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ). From the viewpoint of the heat capacity, the existence of the peak indicates the Barrow’s EPYM AdS black hole should have discrete energy levels microscopically. If regarding the Barrow’s EPYM AdS black hole mass as the internal energy, the mass parameter should be the statistical average of these energy levels. At a enough high temperature of the Barrow’s EPYM AdS black hole, the mass parameter is nearly proportional to the Barrow’s temperature (see Fig. 3). The heat capacity keeps a constant value thereafter.

Refer to caption
(a) T,TBr+𝑇subscript𝑇𝐵subscript𝑟~{}~{}T,T_{B}-r_{+}italic_T , italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT - italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT
Refer to caption
(b) CB,TBr+subscript𝐶𝐵subscript𝑇𝐵subscript𝑟~{}~{}C_{B},T_{B}-r_{+}italic_C start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT , italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT - italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT
Figure 2: The parameters are set to c1=0.02,q=0.85formulae-sequencesubscript𝑐10.02𝑞0.85c_{1}=0.02,~{}q=0.85italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 0.02 , italic_q = 0.85 and γ=0.8𝛾0.8\gamma=0.8italic_γ = 0.8. In the left, the dashed line stands for the behaviour of TBr+subscript𝑇𝐵subscript𝑟T_{B}-r_{+}italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT - italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT, and the solid line is the curve of Tr+𝑇subscript𝑟T-r_{+}italic_T - italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT. In the right, the solid line stands for the behaviour of CBr+subscript𝐶𝐵subscript𝑟C_{B}-r_{+}italic_C start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT - italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT, and the dashed line is the curve of TBr+subscript𝑇𝐵subscript𝑟T_{B}-r_{+}italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT - italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT.
Refer to caption
Figure 3: The plots of the black hole mass vs the Barrow’s temperature with the parameters c1=0.02,q=0.85,γ=0.8formulae-sequencesubscript𝑐10.02formulae-sequence𝑞0.85𝛾0.8c_{1}=0.02,~{}q=0.85,~{}\gamma=0.8italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 0.02 , italic_q = 0.85 , italic_γ = 0.8.

IV Discussions and Conclusions

Based on the Barrow’s idea of the fractal structure for black holes, the authors Ma2022 proposed the area of the smallest bubble sphere is the minimal increase of horizon area. The entropy is of a logarithmic form for the static spherically symmetric black holes, which can be related to the Boltzmann entropy that can reflect the microscopic structure of black holes. Based on these we found that the Barrow’s EPYM AdS black hole with a any given value of ΛΛ\Lambdaroman_Λ is always stable and there is no phase transition as well as the standard EPYM AdS black hole with ΛΛcΛsubscriptΛ𝑐\Lambda\geq\Lambda_{c}roman_Λ ≥ roman_Λ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT, and as Λ<ΛcΛsubscriptΛ𝑐\Lambda<\Lambda_{c}roman_Λ < roman_Λ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT the standard EPYM AdS black hole will undergo a first-order phase transition in which a unstable intermediate-potential black hole can be survived. On the other hand, for the Barrow’s EPYM AdS black hole its heat capacity exhibits a Schottky anomaly behavior, which has also been found and discussed in other gravitational system. This maybe due to the existence of discrete energy levels and restricted microscopic degrees of freedom. For low temperatures, the Barrow’s heat capacity increases rapidly with the Barrow’s temperature, while at high enough temperature its behaviour is like that of any usual thermodynamical systems.

Acknowledgements

We would like to thank Prof. Ren Zhao and Meng-Sen Ma for their indispensable discussions and comments. This work was supported by the National Natural Science Foundation of China (12375050, 12075143), the Natural Science Foundation of Shanxi Province (202203021221209, 202303021211180), the Teaching Reform Project of Shanxi Datong Universtiy (XJG2022234), and the Program of State Key Laboratory of Quantum Optics and Quantum Optics Devices (KF202403).

References

  • (1) J. R. Oppenheimer and H. Snyder, On Continued Gravitational Contraction, Rev. Lett. 56 (1939) 455.
  • (2) B. P. Abbatial, et al, Properties of the Binary Black Hole Merger GW150914, Phys. Rev. Lett. 116 (2016) 241102.
  • (3) J. M. Bardeen, B. Carter, and S. W. Hawking, The Four Laws of Black Hole Mechanics, Commun. Math. Phys. 31 (1973) 161.
  • (4) R. M. Wald, The Thermodynamics of Black Holes, Living Rev.4 (2001) 6.
  • (5) S. W. Hawking, Particle creation by black holes, Commun. Math. Phys. 43 (1975) 199.
  • (6) J. D. Bekenstein, Black holes and the second law, Lett. Nuovo Cimento 4 (1972) 737.
  • (7) J. D. Bekenstein, Black Holes and Entropy, Phys. Rev. D 7 (1973) 2333.
  • (8) N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Spacetime, Cambridge University Press, Cambridge, 1982.
  • (9) E. M. C. Abreu and J. A. Neto, Barrow black hole corrected-entropy model and Tsallis nonextensivity,Phys. Lett. B 810 (2020) 135805.
  • (10) J. D. Barrow, The area of a rough black hole, Phys. Lett. B 808 (2020) 135643.
  • (11) E. N. Saridakis, Barrow holographic dark energy, Phys. Rev. D 102 (2020) 123525.
  • (12) H. Moradpour, A. H. Ziaie, and M. Kord Zangeneh, Generalized entropies and corresponding holographic dark energy models, Eur. Phys. J. C 80 (2020) 732.
  • (13) A. Salehi, Test of Barrow entropy using a model independent approach, Eur.Phys.J.C 83 (2023) 11.
  • (14) N. Komatsu, Cosmological model based on both a holographiclike connection and Padmanabhan s holographic equipartition law, Phys. Rev. D 109 (2024) 2.
  • (15) O. Okcu, Investigation of generalised uncertainty principle effects on FRW cosmology, Nucl. Phys. B 1004 (2024) 116551.
  • (16) E. M. C. Abreu, J. A. Neto, and E. M. Barboza, Barrow’s black hole entropy and the equipartition theorem, EPL (Europhysics Letters) 130 (2020) 40005.
  • (17) L.-H Wang and M.-S Ma, Barrow black holes and the minimal length, Phys. Lett. B 831 (137181) 2022.
  • (18) Y. Ladghami, B. Asfour, A. Bouali, A. Errahmani, and T. Ouali, Barrow entropy and AdS black holes in RPS thermodynamics, Phys.Dark Univ. 44 (2024) 101470.
  • (19) S. Rani, A. Jawad, and M. Hussain, Impact of barrow entropy on geometrothermodynamics of specific black holes, Eur. Phys. J. C 83 (2023) 8.
  • (20) M. Zhang, Z.-Y. Yang, D.-C. Zou, W. Xu, and R.-H. Yue, PV𝑃𝑉P-Vitalic_P - italic_V criticality of AdS black hole in the Einstein-Maxwell-power-Yang-Mills gravity, Gen. Rel. Grav. 47, 14 (2015), arXiv:1412.1197.
  • (21) C. Corda and H.J. Mosquera Cuesta, Inflation from R2superscript𝑅2R^{2}italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT gravity: a new approach using nonlinear electrodynamics, Astropart. Phys. 34 (2011) 587, arXiv:1011.4801.
  • (22) S. H. Mazharimousavi and M. Halilsoy, Lovelock black holes with a power-Yang-Mills source, Phys. Lett. B 681 (2009) 190, arXiv:0908.0308.
  • (23) V. A. De Lorenci, R. Klippert, M. Novello, and J. M. Salim, Nonlinear electrodynamics and FRW cosmology, Phys. Rev. D 65 (2002) 063501.
  • (24) A. B. Balakin, J. P. S. Lemos, and A.E. Zayats, Regular nonminimal magnetic black holes in spacetimes with a cosmological constant, Phys. Rev. D 93 (2016) 024008, arXiv:1512.02653.
  • (25) A. B. Balakin, H. Dehnen, and A. E. Zayats, Nonminimal Einstein-Yang-Mills-Higgs theory: associated, color, and color-acoustic metrics for the Wu-Yang monopole model, Phys. Rev. D 76 (2007) 124011, arXiv:0710.5070.
  • (26) A. B. Balakin and A. E. Zayats, Non-minimal Wu-Yang monopole, Phys. Lett. B 644 (2007) 294, gr-qc/0612019.
  • (27) P. K. Yerra and B. Chandrasekhar, Heat engines at criticality for nonlinearly charged black holes, Mod. Phys. Lett. A 34 (2019) 1950216, arXiv:1806.08226.
  • (28) H. El Moumni, Revisiting the phase transition of AdS-Maxwell-power-Yang-Mills black holes via AdS/CFT tools, Phys. Lett. B 776 (2018) 124.
  • (29) Y.-Z Du, H.-F Li, F. Liu, R. Zhao, and L.-C Zhang, Phase transition of non-linear charged Anti-de Sitter blackh holes, Chinese Phys. C 45 (2021) 11.