The effect of Rényi entropy on Hawking radiation

It is widely believed that Hawking radiation originates from excitations near the horizons of black holes (Hawking in Commun Math Phys 43:199, 1975, Giddings in Phys Lett B 754:39, 2016, Mathur in Class Quantum Gravity 26:224001, 2009). However, Giddings (Phys Lett B 754:39, 2016) proposed that the Hawking radiation spectrum that characterizes evaporating semi-classical black holes originates from a quantum “atmosphere”, which extends beyond the horizon of a black hole. Although several research projects have been conducted in this field, they have not yet taken into account the effect of Rényi entropy. In the present article, we will therefore consider the effect of Rényi entropy on Hawking radiation power. We assume that if the effect of Rényi entropy is very small, we suggest that the Hawking radiation should originate from the quantum “atmosphere” which extends beyond the black hole’s horizon for finite dimensions. That is, that Giddings’ suggestion is the more likely of the above possibilities. However, for infinite dimensions, both suggestions are equally credible. We briefly consider the very large effect of Rényi entropy on Hawking radiation power as well. We find that if the effect of Rényi entropy is very large and ω/TBH\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega /T_{BH}$$\end{document} is very small, then the power spectral density SR\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_R$$\end{document} is proportional to the power spectral density SBH\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_{BH}$$\end{document}.


Introduction
Based on Hawking's original paper, the evaporation process of black holes is characterized by a non-unitary evolution of quantum fields in curved spacetimes [1]. More specifically, matter fields in a pure quantum state will collapse to form a black hole that eventually evaporates into a mixed thermal state [1]. It is widely believed that the semi-classical Hawking radiation spectra should be modified to ensure quantum unitarity since unitary temporal evolution is one of the fundamental principles of quantum mechanics [2,3].
It is generally believed that Hawking radiation spectra originate from quantum excitations near the horizon region, namely, r = r −r H r H [2,3]. Therefore, it was widely expected that modified Hawking radiation spectra should also be characterized by a relatively short length scale r r H , where r H is the radius of the black hole's horizon [3].
However, Giddings suggested that the radiation spectra of a black hole originates from an effective quantum "atmosphere" which extends well beyond the horizon of the black hole [2]. By comparing the numerical results for the Hawking radiation power P B H of an evaporating (3 + 1)-dimensional Schwarzschild black hole of horizon radius r H with the familiar Stefan-Boltzmann radiation power P B B = σ AT 4 of a (3 + 1)-dimensional flat space perfect blackbody emitter of radius r A [4,5], Giddings concluded that the source of Hawking radiation was a quantum region which was located outside the black hole and that the effective radius r A was described by the following relation (1.1) According to Ref. [2], Eq. (1.1) is consistent with the existence of an effective emitting "atmosphere" which extends well beyond the horizon of a black hole, i.e., "atmosphere" is outside and far from the horizon of a black hole.
In considering the correctness of Giddings' conclusion, one interesting question that should be addressed is if the radiation spectra of a black hole in any (D + 1)dimension originates from an effective quantum "atmosphere" which extends well beyond the horizon of the black hole? In order to answer this question, in Ref. [6], the author studied the Hawking radiation spectra of Schwarzschild black holes in a (D + 1)-dimension. Following Giddings's method, the author defined the effective radii r A (D) of the black-hole quantum atmospheres by equating the Hawking radiation powers of the (D+1)-dimensional black holes to the corresponding Stefan-Boltzmann radiation powers of flat space perfect blackbody emitters [6]. The dimensionless radii r A /r H parameters, which characterize the effective black-hole quantum atmospheres, are functions which monotonically decrease as the number (D + 1) increases [6]. In particular, in the large D 1 regime, Ref. [6] suggests that radiating (D + 1)-dimensional Schwarzschild black holes are characterized by the relation (r A − r H )/r H 1. Following Giddings' argument, other studies on black-hole quantum atmospheres have been conducted [7][8][9]. For example, in Ref. [7], Yan-Gang Miao and Zhen-Ming Xu investigated the Hawking radiation cascade from a five-dimensional charged black hole with a scalar field coupled to higher-order Euler densities in a conformally invariant manner. Myungseok Eune and Wontae Kim [8] showed that the outer temperature vanishes at the horizon and has a peak at a scale whose radial extent is set by the horizon radius which then decreases to the Hawking temperature at infinity. Ramit Dey, Stefano Liberati and Daniele Pranzetti [9] have shown that the Hawking quanta originate from a quantum atmosphere around the black hole with energy density and fluxes of particles peaking at about 4MG, which is contrary to the popular belief that these quanta originate from ultra high energy excitations very close to the horizon. At present, the theory of a "black-hole quantum atmosphere" is a frontier area of research which is still far from reaching any final widely accepted conclusion. Results on this topic which are not directly related to the present article are not listed here for reasons of brevity.
One important element that appears to have received little attention in previous articles is the following. In 1902 Gibbs pointed out that, in a system where the partition function diverges, the standard Boltzmann-Gibbs theory cannot be applied, and large-scale gravitational systems are known to fall within this class [10]. Then Tsallis generalized standard statistical mechanics (which arises from the hypothesis of weak probabilistic correlations and their connection to ergodicity) to nonextensive one, which can be applied in all cases, and still possessing standard Boltzmann-Gibbs theory as a limit [10]. Therefore, the usual Boltzmann-Gibbs additive entropy must be generalized to the nonextensive, i.e., non-additive entropy (the entropy of the whole system is not necessarily the sum of the entropies of its sub-systems), which is named Tsallis entropy [10]. Non-extensive statistical mechanics does not only include one example, i.e., Tsallis entropy. Rényi entropy is another main example of such a nonextensive entropy. In this paper, we will therefore principally discuss the effect of Rényi entropy on Hawking radiation power, in particular on the effective radius of a Schwarzschild black hole quantum atmosphere. In Sect. 2, we introduce the basic concepts of Hawking radiation spectra of Schwarzschild black holes and Rényi entropy. In Sect. 3, we derive the effective radius of a black-hole quantum atmosphere taking Rényi entropy into consideration. In Sect. 4, numerical results of the effective radius of the black-hole quantum atmosphere are obtained. In Sect. 5, we briefly discuss when the effect of Rényi entropy is very large, how it impacts on Hawking radiation power. In Sect. 6, we discuss the results we have obtained.

The Hawking radiation spectra of (D+1)-dimensional Schwarzschild black holes
For (D +1)-dimensional Schwarzschild black holes, the semi-classical Hawking radiation power for one bosonic degree of freedom is described as [4,5,11,12] Here j is the angular harmonic index of the emitted field modes and ω is the emitted frequency of field. The Bekenstein-Hawking temperature of the black hole is given by We have set G = c = k B = = 1. r B H is the horizon radius of the black hole [6,13]. The greybody factors of the black-hole-field system composition are denoted as = (ω; j, D), which are dimensionless coefficients which quantify the interaction of the emitted fields with curved black-hole spacetime [4,5].

An area law prescription for Rényi entropy
One of the most remarkable discoveries in fundamental physics is that the entropy of a black hole has a value equal to a quarter of the horizon area in Planck units [3,14]: Here G N is the Newton's constant. Ryu and Takayanagi [14,15] generalized the relationship eq.(2.3) in the context of gauge/gravity duality. They proposed that the von Neumann entropy is determined by the area of a codimensional-2 minimal surface in the dual spacetime [14,15]: The above discussions of area laws eq.(2.3) and (2.4) are limited to the von Neumann entropy. However, according to Gibbs' argument [10], the Boltzmann-Gibbs (BG) theory cannot be applied in systems with divergence in the partition function, such as in a gravitational system [10]. Therefore, thermodynamic entropy of such nonstandard systems cannot be described by an extensive entropy but must instead be generalized to a non-extensive entropy [10]. In Ref. [14], the author has generalized the area laws to the case of Rényi entropy, in which a label index q is used and the entropy is defined in terms of the density matrix ρ of the entangling region as [14]: When q → 1, the von Neumann entropy S = −T r(ρ ln ρ) is recovered. Rényi entropy contains richer physical information about the entanglement structure of a quantum state. The properties of Rényi entropy have been extensively studied using by numerical methods [16], in spin chains [17], in tensor networks [18], in free field theories [19], in two-dimensional CFTs [20] or higher [21], and in the context of gauge/gravity duality [22]. These results have also been generalized to charged [23] and supersymmetric cases [24]. The generalized area-law prescription for holographic Rényi entropy can be derived by applying the replica trick in the context of gauge/gravity duality [14,25]. Rényi entropy (2.5) of integer q > 1 is determined by the partition function of the QFT on a branched cover M q , defined by taking q copies of the original Euclidean spacetime M 1 on which the QFT lives with a cut along the entangling region and gluing them along the cuts in a cyclic order. This can be written as [14,25]:  [14,25]. In the large N limit where the bulk physics is classical, we have [14]: where I bulk [B q ] denotes the on-shell action of the bulk solution B q . If we assume the Z q replica symmetry, then one can obtain that [14]: whereB q = B q /Z q and the bulk action I bulk is the Einstein action Here X μ , G μν and R are the coordinates, metric and Ricci scalar in the bulk [14]. As a result, one can obtain a generalized area law: where the cosmic brane is analogous to the Ryu-Takayanagi minimal surface [14]. The above arguments can be generalized to non-integer q [14]. More details can be found in Ref. [14].

Basic formulae for Rényi entropy
A long-standing problem of non-extensive thermodynamics is deriving a formulation which is compatible with the zeroth law of thermodynamics [26]. Based on the concept of composability, Abe showed [26,27] that the most general non-additive entropy composition rule which is compatible with a homogeneous equilibrium has the form where H λ is a differentiable function of S and λ is a real constant parameter. S 1 , S 2 and S 12 are the entropies of the subsystems and the total system, respectively. Biró and Ván [28] showed that, the most general entropy function compatible with the zeroth law for homogeneous systems has the following form: which is additive for composition, namely, The corresponding temperature is then where we have assumed the principle of additivity holds for the energy composition [26]. For classical black holes, the Bekenstein-Hawking formula satisfies the Eq. (2.3), i.e., where H λ (S) = √ S and λ → 0 [26,27]. The Rényi entropy can be defined as which is equivalent to the choices of H λ (S) = S and λ = 1 − q in eq.(2.4) [26], if the original entropy functions satisfy the following condition where q is a real parameter called the non-extensivity parameter. It can be shown that the formal logarithm of Tsallis entropy (2.18) provides the Rényi entropy, namely, When q approaches to 1 (λ → 0), both Rényi and Tsallis entropy give the Shannon entropy S = − i p i ln p i [26].

The effective radius of a black-hole quantum atmosphere taking Rényi entropy into consideration
In this section, in order to determine the location of effective radius of a Schwarzchild black hole, we will derive the formulae for black hole thermodynamics considering the perturbations of Rényi entropy.

The Rényi temperature of a (D+1)-dimensional Schwarzschild black hole
Considering Eq. Based on Ref. [6], the horizon radius of a (D + 1)-dimensional Schwarzschild black hole of mass M is given by Then the thermodynamic entropy 1 of (D + 1)-dimensional Schwarzschild black hole is We have pointed out that the entropy of large-scale gravitational system cannot be described by Gibbs-Boltzmann statistical mechanics and should be replaced by Tsallis entropy [10]. Then if we regard the entropy of a black hole as being Tsallis entropy, namely, S(r B H ) = S T [26,29], then the Rényi entropy of a (D + 1)-dimensional Schwarzschild black hole is then where λ = 1 − q. We can therefore conclude that irrespective of the value of λ, the value of S R will be positive. The derivative of Rényi entropy can then be obtained as So that the Rényi temperature 2 of a (D + 1)-dimensional Schwarzschild black hole is When D = 3 and λ → 0, the Rényi temperature becomes which is the well-known result for a (3 + 1)-dimensional Schwarzschild black hole [6]. (1 + x), where x = λ π D/2 r D−1 B H 2 (D/2) . In Sect. 3, we consider only the case of |x| 1 and a modification to order x.

The effective radius of a black-hole quantum atmosphere taking Rényi entropy into account
On the one hand, for (D + 1)-dimensional Schwarzschild black holes, the semiclassical Hawking radiation power for one bosonic degree of freedom should be modified as (3.10) So that Eq. (3.9) can be written as On the other hand, the scalar radiation power of a spherically-symmetric blackbody (BB) at temperature T and having radius R in (D + 1) spacetime dimensions is given by the generalized Stefan-Boltzmann radiation law [6] where σ is the Stefan-Boltzmann constant in (D + 1)-dimensions: and A D−1 (R) is the surface area of an emitting body in (D + 1)-dimensions: In this article, we suggest that the Stefan-Boltzmann law (3.14) still holds since it is a non-gravitational system which can be described by the standard statistical mechanics. However, it should be modified as, where T R and r are the Rényi temperature and radius in (D + 1)-dimensions, respectively. The reason can be explained as the following. If we repeat the derivation of Stefan-Boltzmann's law [30], we can know that the only thermodynamic quantity is temperature T , while other physical quantities are mechanical quantities. Now that we have considered the effect of Rényi entropy and Rényi temperature, the only change in Stefan-Boltzmann's law is to change the ordinary temperature T into Rényi temperature T R , and the other parts remain unchanged. Giddings pointed out that the effective radius of a black-hole atmosphere can be determined by equating the radiation power P B B of a flat space perfect blackbody emitter with the Hawking radiation power P B H [2,6]. Therefore, following Ref. [2,6], we can define the effective radius r A of the black-hole quantum atmosphere by the following relation: , T B H ). (3.18) In this article, the eq.(3.18) is modified to, (3.19) where T R (T B H ) means that T R is a function of T B H . Since T R = T B H (1 + x) and |x| 1, therefore  (r B H , T B H ), (3.22) where C D+1 is a constant in (D + 1)-dimensions. Considering Eqs. (3.14) and (3.18), we can therefore derive that, (3.23) which means that our assumption is self-consistent. Combining Eqs. (3.13), (3.17), (3.19) and (3.22), we can obtain the effective radius r AR of a black-hole quantum atmosphere as being, where we have defined the scaled Hawking radiation power of a black hole in (D + 1)dimensions [6] as,P If we define the effective radius r A of the black-hole quantum atmosphere without taking into account the effect of Rényi entropy [6] as (3.26) and the entropy difference as, (3.27) When x = 0, then F = 1, and the effective radius r AR of the black-hole quantum atmosphere, taking into account the effect of Rényi entropy, can be written as There are two reasons why we only expand 1 exp(ω/T R )−1 to order x. Firstly, if we expand the term to order x then the effective radius r AR can be written as the product of the standard effective radius r A and the difference F, which makes for an easier comparison of the results with that given by standard statistical mechanics. Secondly, standard statistical mechanics can describe the universe accurately, therefore, the effect of Rényi entropy can be neglected completely under normal circumstances. In fact, as we have pointed out in Sect. 2.3, when x is zero, the standard results are recovered.
which characterizes the effective quantum atmospheres of the radiating black holes in (D + 1) dimensions. Considering Eq. (3.28), we can rewrite Eq. (4.1) as, wherer A = r A /r B H . Based on Hawking's and Giddings' arguments [2,6], there are two possible values ofr AR in (D + 1) dimensions, namely,r AR ∼ O(0) orr AR ∼ O(1). We will calculate the numerical results for the two possibilities respectively. Forr AR ∼ O(0), if we setr AR = 0, we have, Forr AR ∼ O(1), if we setr AR = 1, we have, where the subscript D denotes the dimensionality of space.

The (3 +1)-dimensional case
The Hawking radiation power of scalar quanta from a (3 + 1)-dimensional Schwarzschild black hole is given by [6,31]  Then we haver A = 2.679 for the effective radius of the black-hole quantum atmosphere. Forr AR ∼ O(0), from Eq. (4.3), we have x 3 = 1.0661, which is O(1) and inconsistent with the assumption in Sect. 3, namely, |x| 1. We cannot adjust the value ofr AR to make |x| 1. Therefore, for the (3 + 1)-dimensional case, Hawking's suggestion is less feasible.
Forr AR ∼ O(1), from Eq. (4.4), we have x 3 = 0.1786. In fact, we can adjust the value ofr AR so that x can be very closed to zero. Therefore, for the (3+1)-dimensional case, Giddings' suggestion is more feasible.

The (D +1)-dimensional cases
In the previous subsection, we have seen that the effective quantum atmosphere for a (3 + 1)-dimensional Schwarzschild black hole is described by the relation  Table 3 The values of x D for the radiating black holes in (D + 1) dimensions using Giddings' suggestion D + 1 5 6 7 8 9 1 0 1 1 Based on the numerical results in (3 + 1) dimensions, we can deduce, that Hawking's suggestion is less feasible than that of Giddings'. We shall now calculate the values of the dimensionless radiir AR , which describe the effective quantum atmospheres of the radiating Schwarzschild black holes, and the value of x, which quantifies the effect of Rényi entropy in (D + 1)-dimensions. In Ref. [6,32], the Hawking radiation powers of Schwarzschild black holes in (D + 1)dimensions have been obtained numerically. We list the results of the dimensionless radii r A (D) in Table 1 [6].

The Giddings' cases for intermediate D-values
Forr AR ∼ O(1), from eq.(4.4) and  Table 3. Based on Table 3, we find that ifr AR ∼ O(1), then the values of x D can be very closed to zero and we can therefore conclude that for any dimensions, Giddings' suggestion appears to be feasible.
Combining this with the conclusions in Sect. 4.1 and 4.2.1, we can deduce that if we consider the effect of Rényi entropy, for any dimensions, Hawking's suggestions appear to be much less feasible that those of Gidding. Therefore, Hawking radiation most likely originates from a quantum atmosphere which extends beyond the horizons of black holes. We can also infer that the above conclusion is a general result and true independently of the value of λ in non-extensive statistical mechanics.

The large D regime
In this section, we will consider the case for in the large D regime to see whether the conclusions which we have obtained in Sects. 4.1, 4.2.1 and 4.2.2 are still established. The Hawking radiation spectrum for a Schwarzschild black holes in (D + 1) dimensions is characterized by ω D /(exp(ω/T R ) − 1) (see Eq. (3.9)). This frequency dependent function has a peak at [6] ω peak where W (x) is the Lambert function and T R is the Rényi temperature of the black hole as given in Eq. (3.7). Inserting Eq. (3.7) into Eq. (4.7), we then have λ peak where λ peak is the peak of wave length. Since |x| 1, the modification term does not affect the above conclusion, namely, that in the large D regime the characteristic wavelength is much shorter than the horizon radius r B H of a radiating Schwarzschild black hole, which is consistent with the result obtained without taking Rényi entropy into consideration.
References [6] and [32] have pointed out that the inequality (4.8) implies that the corresponding Hawking emission spectra of these higher-dimensional black holes are described extremely well by a short wavelength approximation. In particular, the effective radiating r A of a black hole in the large D regime is determined by the highenergy (short wavelengths) absorptive radius of the black hole [6,12,[32][33][34][35][36][37][38], namely, In the large D 1 regime, we find that ,  1)x ≈ −1, which is the necessary condition for Giddings' arguments. We also find that λ < 0.
Since D 1, we know that |x| is extremely small for the two cases, so that the effect of Rényi entropy can be neglected completely. In addition, these two suggestions are both feasible.

Very large effect of Rényi entropy on Hawking radiation power
In this section, we will briefly mention the properties of Hawking radiation power in the situation where the effect of Rényi entropy is extremely large. Although this is a situation which does not exist in the world we have explored, it is nevertheless of interest. on Hawking radiation power. If the effect of Rényi entropy is very large and ω/T B H is very small, we found that the power spectral density S R is proportional to S B H , that is, S R = x S B H , where x is defined as in Sect. 3.1 and includes the effect of Rényi entropy.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.