Prospecting Black Hole Thermodynamics with Fractional Quantum Mechanics

This paper investigates whether the framework of fractional quantum mechanics can broaden our perspective of black hole thermodynamics. Concretely, we employ a {\it space-fractional} derivative \cite{Rie} as our main tool. Moreover, we restrict our analysis to the case of a Schwarzschild configuration. From a subsequently modified Wheeler-DeWitt equation, we retrieve the corresponding expressions for specific observables. Namely, the black hole mass spectrum, $M$, its temperature $T$, and entropy, $S$. We find that these bear consequential alterations conveyed through a fractional parameter, $\alpha$. In particular, the standard results are recovered in the specific limit $\alpha=2$. Furthermore, we elaborate how generalizations of the entropy-area relation suggested by Tsallis and Cirto \cite{Tsallis} and Barrow \cite{Barrow} acquire a complementary interpretation in terms of a fractional point of view. A thorough discussion of our results is presented.


Introduction
Black hole (BH) physics constitutes an enthusiastic research domain for this century. On the one hand, gravitational waves resonating from colliding BHs have been measured [4] and are now a regular astronomical probing tool. Outstandingly, astrophysicists have also captured the first-ever image of a BH at the center of galaxy M87 [5]. Furthermore, there is a collection of additional observational data [6] hinting at the presence of an event horizon for BHs. In general, BH candidates may be stellar-mass objects in a X-ray binary system or supermassive at the center of typical galaxies [7,8, a e-mail: shahram.jalalzadeh@ufpe.br b e-mail: filipe.rodrigues@ufpe.br c e-mail: pmoniz@ubi.pt 9]. All these achievements have been recently conveyed and recognized within the 2020 Nobel prize in Physics [10].
On the other hand, innovative work on BH physics in the past forty years or so have also brought us to the forefront of unanticipated research directions. In the early 70's of the last century, Bekenstein proposed that a BH entropy is proportional to its horizon area [11]. Then, the evaporation of BHs (through the seminal use of quantum fields in curved spaces) was formulated by Hawking [12]. Testing this property observationally has been recently appraised, bearing promising results [13]. Thereafter, physicists realized an intimate connection between geometrical horizons, thermodynamic temperature and quantum mechanics [14]. Moreover, it was advanced [15] that the BH horizon area may be quantized, and the corresponding eigenvalues would be given by A n = γL 2 P n, n = 1, 2, 3, ... , where γ is a dimensionless constant of order one and L P = √ G is the Planck length. The literature has since then been enriched with contributions strengthening in favor of the area spectrum (1), including informationtheoretic considerations [16], ranging from string theory arguments [17] to the periodicity of time [18], plus e.g., a Hamiltonian quantization of a dust collapse [19].
Within the viewpoint outlined in the previous paragraph, we suggest a complementary direction of exploration. Let us be more specific. Recently, a set of pertinent arguments and results have been put forward to apply fractional calculus [20,21,22,23] in quantum physics. Such framework is known as Fractional Quantum Mechanics (FQM), see e.g. [24]. Fractional calculus has been embraced mainly within the last century. In essence, it follows from extending the meaning of derivatives to the case where the order is any number, i.e., irrational, fractional or complex. Its peculiar interest notwithstanding, the obstacles have been fruitfully addressed. In particular, it is currently known that non-integer order systems can describe the dynamical behavior of specific classes of materials and processes over different time and frequency scales. Fractional calculus has assisted in scattering theory, diffusion, probability, potential theory and elasticity. Therefore, it was only sensible to embrace and explore it within quantum physics.
FQM was built upon Feynman and Hibbs [25] assertion about quantum mechanics and path integrals. Moreover, Nelson [26] clarified about classical Brownian motion and quantum mechanical features, whereas Abbott and Wise [27] determined that in 1-dimensional settings, the fractal (or Hausdorff) dimension of those paths is 2. The essential motivation for FQM emerges in that if we restrict the path integral (Feynman) description of quantum mechanics to Brownian paths only, it will be challenging to explain a few other pertinent quantum phenomena [28]. Such difficulties have led to consider a generalization of the Feynman path integral, specifically by replacing the Gaussian probability distribution by Lévy's [29]; The Hausdorff dimension of the Lévy path is then equal to the fractional parameter α.
Extended versions of the Schrödinger equation (SE) are then retrieved, namely from such broader path integral. More concretely, by including non-Brownian trajectories in the path integral formalism of quantum physics, either the space-fractional [30], time-fractional [31], and space-time-fractional [32] versions of the standard SE can be elaborated. Let us be precise and clarify that it is in space-fractional quantum mechanics, where the Feynman path integral method is modified such that the Gaussian probability distribution is replaced by Lévy's [29], that we obtain the following and broadly used modified SE. More particularly, starting from the Hamiltonian a generalization [24,28,30] is produced by means of in which D α is a coefficient carrying dimension [D α ] = erg 1−α cm α sec −α . The parameter α is known as Lévy's fractional parameter and is associated to the concept of Lévy path and its induced fractal dimension. In the Feynman path integral, the measure is generated by the process of the Brownian motion, and the path's corresponding (fractal) dimension is d For more details, please browse through, e.g., [1,20,21,22,24,28,30] and other references therein.
Therefore, choosing a space representation withp → −i ∇ andr → r, we obtain the space-fractional SE In the above fractional SE, the fractional Riesz derivative [1,24,33], (− 2 ∆) α/2 , is defined in terms of the Fourier transformation F The infinite-well example was one of the first solutions of space-fractional SE, which Laskin solved [28]. Despite its simplicity, this problem is critical since it is the prototype of a quantum detector with internal degrees of freedom.
In what concerns the time-fractional SE, the time evolution is given by the Caputo fractional derivative [34]. In this case, the associated Hamiltonian is non-Hermitian and not local in time. The space-time fractional SE has been introduced by Wang and Xu [32], where they employed a combination of space and timefractional models to establish their equation. The spacetime fractional SE is a generalized version of the equation (4): where 1 < α ≤ 2, 0 < β ≤ 1, β and D α,β are two scale coefficients with physical dimensions [ β ] = erg.sec β and [D α,β ] = erg 1−α .cm α .sec −αβ . In addition, ∂ β t denotes the left Caputo fractional derivative [35] of order β: whereḟ (τ ) = df (τ ) dτ . Proceeding towards a broader context, we recall that Wheeler [36] seminally suggested a foam structure for spacetime on the Planck scale. Hence, allowing that a fractal quality could therein be natural. Thus, contemplating a fractional WDW equation could unveil interesting aspects regarding the gravitational domain [37,38]. In fact, in the last few years, FQM has indeed been proposed and employed as a tool to explore features within quantum cosmology and quantum gravity [39]. It 3 has pointed to exciting opportunities and sidelong connections between unexpected mathematics and physics domains. See [37,38], for a recent survey about engaging FQM within a canonical route towards quantum gravity and cosmology. Additionally, a few more other references can be found in [40].
Thus, in our paper we propose to apply space-FQM to investigate thermodynamic properties of the Schwarzschild BH. The Bekenstein-Hawking entropy expression for a Schwarzchild BH with mass M is given by where A = 16πG 2 M 2 is the BH horizon area. A suitably extended fractional version of the Wheeler-DeWitt (WDW) equation for the Schwarzschild BH will enable us to retrieve altered expressions for the entropy and other specific BH observables. Moreover, Tsallis and Cirto recently proposed a new formulation for the Schwarzschild BH's horizon entropy and the entropy-area relation [2]. They asserted that since the Bekenstein-Hawking entropy is proportional to the horizon area (instead of its volume), it would hint that Boltzmann-Gibbs entropy would be unsuitable to describe BHs. These authors have proposed a modified form upon the BH entropy given by Eq. (8), such that an extended expression, using a generalized non-additive entropy, would be where γ is an unspecified constant and δ denotes the non-additivity parameter. The above equation suggests that the entropy is a power-law function of its area. In addition, by introducing a fractal structure for the horizon surface of the BH, Barrow [3] conceived a 3-dimensional spherical analog of a Koch Snowflakea sphereflake-using a hierarchy of touching spheres around the event horizon. This fractal structure of the BH surface changes its actual area, which in turn leads us to a new entropy relation, namely, where A P is the Planck area and △ would suggestively denote an induced deformation of the horizon, with △ = 0 reproducing the conventional Bekenstein-Hawking entropy (simplest horizon structure) and with △ = 1 corresponding to the most intricate structure. Note that the Barrow modified entropy resembles Tsallis and Cirto's non-additive entropy (9); nevertheless, the involved foundations and physical principles are entirely different. As an additional application of BH prospecting by means of FQM, we will retrieve these two extended entropy expressions but without using generalized non-additive features.
In essence, we propose and employ herewith a spacefractional formulation of the WDW equation for the Schwarzschild BH. Hence, let us then mention that the remaining of this paper is organized as follows. In section 3, we summarize the canonical quantization of the Schwarzschild BH, having briefly reviewed the corresponding Hamiltonian setting in section 2. In Section 4, we hypothesize a fractional calculus setting for the corresponding WDW equation as applied to the Schwarzschild BH. Specifically, we consider the observables reviewed in section 3 and present them (e.g., the Schwarzschild BH entropy) from within the framework of FQM. Section 5 contains our conclusions and a thorough discussion. Throughout this paper we shall work in natural units, = c = k B = 1.

Canonical quantization of a Schwarzchild black hole
We recall in this section results concerning the canonical quantization of the Schwarzschild BH that will be of relevance for our analysis in subsequent sections. The spherical symmetric ADM line element is where dΩ 2 is the line element for the unit two sphere S 2 . We follow Kuchař's fall-off conditions [41]. Firstly, it confirms that the coordinates r and t are extended to the Kruskal manifold, −∞ < r, t < ∞, and secondly, that the spacetime is asymptotically flat as well. Likewise, the fall-off conditions guarantee that the 4momentum at the infinities, r → ±∞, has no spatial component, which indicates that the BH is at rest concerning the left and right asymptotic Minkowski spacetimes. By fixing the asymptotic values of the lapse function, N , at infinities, r → ±∞ to be t-dependent quantities (denoted by N ± (t), respectively), the Hamiltonian form of the Einstein-Hilbert action functional, with appropriate boundary terms, reads then (12) in which the conjugate momenta of Λ and R are M P = 1/ √ G is the Planck mass,ḟ = ∂ t f , f ′ = ∂ r f and the quantities M ± (t) are defined by the asymptotic fall-off of the configuration variables. Note that on a classical solution, M ± are equal to the Schwarzschild mass of the BH. Furthermore, H and H r are the super-Hamiltonian and the radial super-momentum constraints, respectively, given by Following Kuchař [41] introducing the following two (M, Π M ) and (R, Π R ) pairs of canonical transformations the action (12) becomes where the new super-Hamiltonian, H, and the supermomentum, H r , are The variation of ADM action (17) with respect to the lapse function N and N r imposes the Hamiltonian and momentum constraints or equivalently The constraint M ′ = 0 means that M is homogeneous M = M (t). Now, by substituting Π R ≈ 0 and M = M (t) back into the action (17) and in addition the new conjugate momenta of M defined by we obtain Note that the new conjugate variables (M, P ) obey the Poisson bracket {M, P } = 1. Following Louko and Mäkelä [42], if we picked the right-hand side asymptotic Minkowski time as the observer time parameter, we should restrict N + = 1 and N − = 0. Next, the reduced action (21) reduces to the following simple form where H(M ) = M is the reduced Hamiltonian of the BH. One can easily show that the field equations' solution is M = const. and P = −t, as expected. The BH's mass constancy follows Birkhoff's theorem, which declares that the mass is the only time-independent and coordinate invariant solution. Moreover, the conjugate momenta, P , describes the asymptotic time coordinate at the spacelike slice. Once again, by using the canonical transformations (M, P ) → (x, p), introduced by Louko and Mäkelä [42] the reduced action (22) takes the form where H = M is the Hamiltonian and is given (utilizing the second transformation of (23)), as It is of interest to acknowledge the following at this point. Transformations (23) map the BH solution into a wormhole solution, in which x represents the wormhole throat [42]. The time evolution of this corresponding wormhole throat is given by Hamilton's equationṡ with a solution for the wormhole throat as As we can see from the above solutions, it is necessary to define transformations (23) so that the time parameter t (or equivalently the conjugate momentum P ) is restricted into the finite interval where T H is, remarkably, the Hawking temperature of the BH.

Thermodynamic implications from the quantum Schwarzschild black hole
In the coordinate representationp → −id/dx,x → x, the canonical quantization procedure upon the previous section gives us a suitable time-independent WDW equation, for the simple one-dimensional minisuperspace of the Schwarzschild BH As usual, there is an operator-ordering problem in addressing the above WDW equation, but as we are interested in BH states whose mass M ≫ M P , its particular resolution will not surpass our semiclassical considerations. Hence, we adopt the following wide enough factor-ordering [43] where i+j +k = −1. Note that because in the left hand side of (30) we have p 2 /x, so, in its right side the power of x is also −1. We thus write the WDW equation (29) where q = (ij + ik + jk − 2)/3. If we redefine the wave function ψ(x) → √ xψ(x), and choose ordering in which q = −3/4, then Eq. (31) will reduce to which is a SE for the harmonic oscillator of the Planck's mass M P , the Planck's angular frequency ω P = 1/t P defined in terms of Planck's time t P = 1/M P . The domain of definition for x is x ≥ 0, and consequently, the Hamiltonian operator of the harmonic oscillator in (32) is defined on a dense domain C ∞ (0, +∞). Hence, H is not an essentially self-adjoint operator. It may constitute an Hermitian operator if A necessary and sufficient condition for the validity of this condition is where a prime symbol, ′, denotes the derivative of ψ(x) with respect to x. As pointed out by Tipler [44], the constant γ (with dimension of length) would be a new fundamental constant of theory. To avoid such, we set it to be zero. Hence, we assume which constitutes the DeWitt boundary condition [45].
In addition, we are interested in square-integrable wave functions in the interval 0 ≤ x < +∞, which implies the second boundary condition ψ(x → +∞) = 0. The general solution of the WDW equation (32) is where N is a normalization constant. Applying the boundary condition (34) on the obtained wave function (38) gives us The solution of the above equation for M , gives us the mass spectrum of the BH.
We can see that the spectrum of M 2 is not equally spaced for small values of ν. On the other hand, for large values of ν (which is defined by if we use the asymptotic relation (37) in (39) we find which provides us the mass spectrum where n is an integer and n ≫ 1. Note that the above result assumes that M ≫ M P , i.e., ν ≫ 1 and is therefore essentially semiclassical. To see this, let us apply the Bohr-Sommerfeld quantization rule to the Hamiltonian (25). The classical turning points of (25) are x = 0 and x = 2M/M 2 P . Hence, the conventional treatment of the Bohr-Sommerfeld quantization rule yields which allows to extract the mass spectrum (41). Bekenstein [15] firstly found a similar mass spectrum. Generally, the proportionality constant for the square root of n in (41) is model dependent. Since then, several authors [15,46] have used different arguments for a quantum BH spectrum of the type (41). Before proceeding, let us keep in mind equation (41), as this and others in this section will bear alterations that will be brought from the intrinsic features of FQM, as applied to the Schwarzshild BH (see next section). S. Hawking [12,47] showed in 1974 that due to quantum fluctuations, BHs emit black-body radiation, consistently with the corresponding entropy being one fourth of the event horizon area, namely A = 16πG 2 M 2 . Following Refs. [48] and [49], let us assume that the Hawking radiation of a massive BH where M ≫ M P and n ≫ 1, is emitted when the BH system spontaneously jumps from the state n + 1 towards the closest lower state level, i.e., n, as described by (41). Let us now denote the frequency of the emitted thermal radiation as ω 0 . Then which agrees with the classical BH oscillation frequencies which scales as 1/M . We thus find a BH to radiate with a characteristic temperature T ∝ M 2 P /M , matching the Hawking temperature.
The characteristic BH time (the lifetime of the BH at the state M (n + 1) before decaying into the lower state M (n)) can be defined [49] as whereṀ = dM/dt is the mass loss of the BH because of its evaporation; in the second equality we used the definition of ω 0 expressed in (43). As discussed in [48,50], because of the interaction of the BH with the vacuum of the quantum fields, the width of the states, W n , is not zero. The width of state n can be estimated [48,50] as where β ≪ 1 is a numerical dimensionless factor. Then, by inserting (43) into the uncertainty relation W n τ n ≃ 1 and eliminating τ n in resulting equation, with the assistance of (44) we can finḋ If we further assume, on the one hand, that the origin of the Hawking radiation emerges from the highly blueshifted modes just outside the horizon, and, on the other hand, take the BH as a black-body, then the radiated power is given by the Stefan-Boltzmann law [51,52]  where σ S = π 2 /60 is the Stefan-Boltzmann constant and A = 16πM 2 /M 4 P is the horizon area. Eliminating the mass loss of the BH, using Eqs. (46) and (47), gives us the effective temperature The BH entropy can then be expressed as By choosing β = 1/15360π [53], Eqs. (48) and (49) yields where S B-H = 4πGM 2 is the Bekenstein-Hawking entropy. The logarithmic correction to the Bekenstein-Hawking entropy is obtained using other methods [53] except that the overall factor 2π is model dependent.

Fractional quantum mechanics and Schwarzschild black hole thermodynamics
In this section, we will establish and then investigate the alterations brought from employing FQM features towards a Schwarzschild BH. This will provide broader expressions for some thermodynamic observables, specifically in terms explicitly dependent on the fractional parameter, 1 < α ≤ 2, since the limit α = 2 will conveys us to the standard case (summarized in sections 2 and 3).
Concretely, a fractional WDW equation for a BH will be employed, in a corresponding minisuperspace. Our approach stems from Laskin's seminal papers [28,54], leading to an Hamiltonian, which includes a fractional kinetic term in terms of the quantum Riesz fractional operator; this method has been extended to the WDW equation in Refs. [37,38].
The fractional extension of the WDW equation (32) is given by where z = x − M/M 2 P is the new coordinate in the 1-dimensional minisuperspace, ∆ = d 2 /dz 2 , (−∆) α 2 is the Riesz fractional derivative [1,20,21,22,55] and 1 < α ≤ 2. Unfortunately, there is no known general solution, explicitly bearing a dependence on α for the above fractional WDW equation. Therefore, we may resort to employ the Bohr-Sommerfeld quantization rule. By replacing |p| = (−∇) 1 2 , the fractional WDW equation leads to the following relation which is, in fact, the fractional version of (25). Hence, the classical turning points (where |p| = 0) are z = where y = (M α+2 P /M 2 ) 1 α z and B(a, b) is the beta function. Thus, the application of the standard Bohr-Sommerfeld quantization rule in this setting gives us the following semi-classical mass spectrum For α = 2 we recover spectrum (41) as written in the previous section. Fig. (1) shows the BH mass spectrum for three values of α. As we show in this figure, the mass of the BH increases with n but at a faster rate for α = 2, namely the standard case. Similarly to the previous section 3, if we write the frequency of the radiation emitted by the BH, ω 0 , we 8 now obtain . Note that the energy of the emitted radiation is a function of the fractional-order α. We then find that for a massive BH, where M ≫ M P , the energy of emitted radiation is minimal for α → 1, and it will increases up to a maximum at α = 2. Additional elements are conveyed within Fig. (2). Interestingly, the mass spectrum (54) (which includes the standard case (41)) may induce observable signatures that gravitational waves [56] may inform about.
To be more clear, let us contemplate the mass spectrum as in (54) and further discuss it within a wider analysis as follows. Take the interaction of a gravitational wave with a BH. Suppose the BH is initially at state n 1 . Regarding the mass spectrum in (54), we can also consider process whereby the BH could alternatively absorb a gravitational wave whose frequency, f GW satisfies where ∆n = n 1 − n 2 and n 2 denotes the final state of the BH after absorption of the gravitational wave. As known in classical GR, the effective potential that describes the motion of a test particle has a maximum at the 3/2 Schwarzschild radius, called photon sphere. Essentially, this (potential barrier) screens the nearhorizon region from external observers [57]. The gravitational radiation proceeding towards the horizon will be scattered at the horizon if its frequency does not match, e.g., (56). Otherwise, radiation with the frequencies such as (56) in our example, would be absorbed.
The previous assertions notwithstanding, if then we consider back-reaction quantum effects, BHs are not perfect absorbers of gravitational radiation and part of the radiation will be reflected from the horizon area.
The reflected part of the radiation then interacts with the potential barrier at the photon sphere. Then, it will partially be transmitted, and the other portion will be reflected again back towards the horizon. This outcome generates a series of so-called gravitational-wave echoes [58]. This feature could eventually be detected either in the inspiral stage of BH binaries or during the last stages of BH relaxation following a merger. For further details and detection methods please see [59]. Therefore, BHs may act as "magnifying lenses" in the sense that they could bring new features of the BH horizonarea within the realm of gravitational waves observations. Furthermore, using the formula of the characteristic time of an evaporating BH, τ −1 =Ṁ /ω 0 , the uncertainty relation for width of states, W n τ n ≃ 1 and the relation W n = βω 0 , we computė which will increase as mass is lost. Inserting the above relation into the Stefan-Boltzmann law (47), it gives the corresponding temperature of the BH 1 T = Note that for α = 2 we will recover equations of the previous section. Moreover, using (49), we find the fractional entropy of the BH to be To further contrast a BH with FQM features with a standard Schwarzschild BH, let us examine Eqs. (55), (57), (58) and (60) and for a stellar BH with mass approximately 10M ⊙ . For the limiting α = 2 and α → 1 values of the Lévy index, the values of S(α), T (α), ω 0 (α) andṀ (α) are given in Table (1). Its first row,  Table 1 The fractional entropy, S(α), the temperature, T (α), the mass loss,Ṁ(α) and the frequency of the emitted thermal radiation, ω 0 (α) of the Schwarzschild black hole given by Eqs. (55), (57), (58) and (60)  together with Fig. (2), shows that the entropy of fractional BH can be much larger than in the standard Schwarzschild case. The second row shows that a stellar BH temperature (with mass 10M ⊙ ) in fractional formalism is always less than 10 −9 K. On the other hand, in the extreme case α → 1, the temperature of a fractional BH is approximately zero: T (α → 1) ≃ 10 −48 K. Moreover, the manner the mass decayṀ (α) proceeds is indicated by by the third row, in agreement with the temperature range: it is quite negligible for the extreme case, α → 1. Since the average temperature of the Universe at the present epoch is about 2.7 K, all stellar BHs with M = 10M ⊙ and 1 < α ≤ 2 are absorbing more matter and radiation than they emit Hawking radiation and will not begin to evaporate until the Universe has expanded and cooled below their corresponding temperature. Therefore, in our setting, the BHs with fractional features could be almost eternal within α → 1. Let us now remark the following. Tsallis and Cirto [2] investigated the entropy of a Schwarzschild BH, using appropriate non-additive generalizations for d−dimensional systems and suggested a generalized entropy. In their study [2], a non-additive entropy is defined by for a set of N discrete states, where δ > 0 denotes the non-additivity parameter and p i is a probability distribution [61]. For δ = 1 we recover the standard Boltzmann-Gibbs entropy. Then, as Tsallis and Cirto demonstrated [2], the generalized BH entropy can be written as If we identify the non-additivity parameter of Tsallis and Cirto as δ = 2+α 2α , the entropy (60) retrieved from fractional quantum mechanical methods can be rewritten as where This shows that the leading term of the FQM computed BH entropy can be expressed in terms of the Tsallis and Cirto result for a BH. Noteworthy, for 1 < α ≤ 2 the fractional entropy of a BH varies in the corresponding interval as A natural question that arises is the meaning of the entropy relation, ranging from Eqs. (60) to (65). To assist in addressing this question, let us first remind the fractal nature of the Lévy path [29] in FQM. As we know, in the Feynman path integral representation of quantum mechanics, the measure is generated by the process of the Brownian motion, and the fractal dimension of the Feynman's path is d (Feynman) fractal = 2. In addition,, FQM is based on the fractional path integral [28,29]. The Lévy path integral corresponding to the SE (51) is where V (z(t)) is the potential as a functional of the Lévy path and z i (z f ) denotes the initial (final) point.
In this case, the fractional path integral measure is .
The measure (67) indicates that a length increment, ∆z = z j − z j−1 , and a time increment, ∆t, satisfy the fractional scaling relation The above scaling relation implies that the fractal dimension of the Lévy path is d (Lévy) fractal = α. Thus, fractional features within FQM could be interpreted as being generated by a Lévy stochastic process [29] and suggesting a fractal structure 2 . Hence, it may be possible to theorize within the FQM framework herewith proposed, that the horizon of the BH may unveil a fractal structure whose relevance will be α-dependent. In this context, let us rewrite the first term on the right-hand side of (60) as follows by means of which we hence define A fractal . Being more concrete, the entropy is proportional to a surface area, but one explicitly given by Furthermore, the above expression additionally suggests a fractal dimension [62] of the BHs surface, namely 2 According to Mandelbrot "A fractal is by definition a set for which the Hausdorff-Besicovitch dimension strictly exceeds the topological dimension" [62]. For example, the Hausdorff dimension of a regular Brownian surface is 2.79 and a triangular von Koch fractal Surface is 1 + log 2 (3) = 2.5849 .
Moreover, equation (60) will take the following form Let us add further to the question set a few paragraphs above, namely on the meaning of entropy within Eqs. (60) to (65). By means of a very interesting paper, J. Barrow recently conjectured that BHs might exist with such extremely wrinkled surfaces (referred to as a rough horizon), so that the event horizon would be a fractal surface [3]. He proposed that the fractal area of the BH, A fractal is related to the ordinary area, where △ = 0 is corresponding to the ordinary nonfractal horizon, and △ = 1 to the most complex structure. Inserting the above area definition into entropy (60) we find where we defined △ = 2/α − 1. Note that the extreme case of Barrow △ = 1, (α = 1), does not exist in our model. Therefore, the above equation clarifies that the Barrow entropy can be obtained from within FQM. We can elaborate more on the physical consequences of the above entropy (or equivalently from (60)) formula. If we assume that the number of degrees of freedom, N , in the horizon is [63] then, combining Eqs. (58), (75) and (76) plus some algebra, we retrieve a modified equipartition theorem where and T P is the Planck temperature. The heat capacity, C of a fractional BH can be computed from the expression where a prime means a derivative relative to the BH mass, M . Substituting (75) into the above definition which shows that for any value of 0 < △ < 1, the heat capacity in our semiclassical model, M ≫ M P , is negative and consequently the BH is unstable. Thus, from the geometrical BH surface, we may admit that there is more to unveil. In particular, we may speculate that a behavior, suggesting a broader dimensional dynamics at the horizon, can be explored in the limit when Lévy's parameter approaches α → 1, supporting the conjecture proposed in [3]. We may, therefore, be allowed to speculate about the possibility that at some scale, a broader surface area of a BH can be considered, encompassing the standard value 16πG 2 M 2 , because of the fractal structure of the horizon. In this context, let us recall Eq. (73) regarding how the entropy varies and how the Bekenstein-Hawking expressions fits in. It brings an interesting perspective such that the entropy of a BH (according to FQM) could be proportional to a power of its fractal surface, depending on the choice of α.

Discussion and Outlook
The purpose in this paper was to employ FQM to discuss specific BH features. It is now well demonstrated that BHs can form [10] and a set of properties has been widely described. Moreover, BHs encompass scenarios where strong gravity and quantum mechanics can meet. New insights need therefore to be proposed. Thus, this paper was set up to explore whether thermodynamic observables might be altered if an additional ingredient is present (say, imported from a quantum mechanical framework) and if it could reveal itself in special BH situations.
Hence, we started by briefly summarizing how concrete thermodynamic properties can be retrieved from a simple quantum mechanical BH (Schwarschild). This was done in sections 2 and 3. We then proceeded to an analysis employing a space-fractional derivative in a corresponding WDW equation. A broader perspective about BH thermodynamics within FQM was therefore presented in section 4, specifically in terms explicitly dependent on the fractional parameter, 1 < α ≤ 2. We retrieved altered expressions for physical quantities that were then contrasted with those of the standard case (α = 2) and summarized in sections 2 and 3.
Then in section 4 we started by elaborating on the mass spectrum of a BH and how it varies according to α. Whereas the frequency of gravitational waves generated by a standard BH is proportional to the inverse mass, in a fractional scenario the frequency would instead allow a dependence as 1/M 4 α −1 . As pointed out in [59], the horizon behaves as a filter of gravitational waves with a set of absorption lines. Admitting the possibility of the fractal structure of the horizon, the horizon radius is different from the Schwarzschild radius [3], and the effective radius is wider, r = 2M G(1 + ǫ), in which ǫ ≪ 1 depends on how the surface is correspondingly altered. Any subsequent gravitational waves may bear a delay due to the time it takes the signal to transit from the horizon (with such possible fractal structure) to the photon sphere. This delay time scales logarithmically with ǫ [59] or, as we suggested herein, the fractality of the horizon.
In addition, we have shown how the entropy, the temperature, and the mass loss of the BH are highly sensitive to the value of α. The Bekenstein-Hawking entropy of supermassive black holes [64] is of the order of 10 104 and the corresponding fractional entropy of the supermassive BHs with α ∼ 1 will be 10 156 . If we extend a FQM perspective to the cosmic event horizon (CEH), then we can compute that S CEH ∝ (GH 2 ) −(α+1) 2α ≃ 10 183 . This assumption has been interestingly employed to discuss the late-time acceleration of the Universe [68]. Furthermore, we retrieved a formula for the Schwarzschild BH in which the entropy is a polynomial function of its area. We then explained that the surface area of such BH can exceed the standard value 16πG 2 M 2 of a BH, as suggested by a possible fractal structure. In this context, Tsallis and Cirto's formalism [2] must be mentioned, constituting an extension of Boltzmann-Gibbs statistical theory, representing a non-extensive, i.e., a non-additive entropy 3 We have indicated herewith that from a FQM perspective we can alternatively, obtain Tsallis and Cirto's results without committing to the extended entropy features. In addition, we proposed that the inherent fractal nature of Lévy paths leads to the extension of Barrow entropy formula [3]. Barrow's modified entropy formally coincides with Tsallis and Cirto's definition; notwithstanding, the involved foundations and interpretation of physical results are different. This is not a surprising result because, over the last several years, it has been shown that a fractal structure in properties of physical systems could be a possible origin for non-additive statistics [69].
Having summarized and discussed our results, it is also pertinent to suggest an outlook of subsequent directions to prospect. A particular line in our sight is to explore BH entropy and information loss within FQM. We plan to extend the context of paper [70] with the assistance of space-fractional derivative operators in the corresponding equations therein. The results in [70] indicated that for the final stage of a BH evaporation, which included back reaction from Hawking radiation, we get a strong entangled state (between a mimicked BH and the radiation). So, given how the entropy may vary with α, what would FQM add within a suitable set up about information and BHs?
It would also be of interest to investigate BHs within FQM and in a Brans-Dicke or even a scalar-tensor theory, where more parameters are allowed to juggle with. Or employ other dynamical configurations instead, beyond the Schwarzschild case. It is possible using other fractional derivatives as well, besides space-fractional. Most likely, the use of numerical methods would be mandatory.