sign in sign up
<<
Home , Black hole

OCU-PHYS-536 AP-GR-167

Hawking radiation from squashed Kaluza-Klein holes with quantum effects

Ken Matsuno∗ Department of Mathematics and , Graduate School of , Osaka City University, Sumiyoshi, Osaka 558-8585, Japan

Abstract We study the from the five-dimensional charged static squashed Kaluza- Klein black hole by the tunneling of charged fermions and charged scalar , including the phenomenological effects predicted by the generalized uncertainty principle with the minimal measurable length. We derive corrections of the Hawking to , which are related to the energy of the emitted , the size of the compact extra dimension, the charge of the black hole and the existence of the minimal length in the squashed Kaluza-Klein geometry. We show that the quantum gravity effect may slow down the increase of the Hawking temperature due to the radiation, which may lead to the thermodynamic stable remnant of the order of the Planck after the evaporation of the squashed Kaluza-Klein black hole. We also find that the sparsity of the Hawking radiation may become infinite when the mass of the squashed Kaluza-Klein black hole approaches its remnant mass. arXiv:2104.00891v2 [hep-th] 17 Sep 2021

∗Electronic address: [email protected]

1 I. INTRODUCTION

Hawking radiation is one of the interesting phenomena where both of general relativity and quantum theory play a role. While a black hole does not emit radiation from a classical point of view, it is shown that, using the semiclassical method of quantum theory in curved , a black hole can radiate from an horizon like a blackbody at the Hawking temperature which is proportional to the of the black hole [1]. Hawking radiation and its black hole temperature are derived by several approaches. Parikh and Wilczek proposed an elegant derivation of Hawking radiation on the basis of the tunneling method [2]. The essential idea of the tunneling mechanism is that particle-antiparticle pairs are formed close to the horizon inside a black hole. While the ingoing particles moving toward the center of the black hole are trapped inside the horizon, a part of the outgoing particles escapes outside the horizon by the quantum tunneling effect. If the particle which comes out to our Universe has positive energy, we can regard such a particle outside the horizon as the radiation from the black hole. The Hawking temperature is derived by comparing the tunneling probability of an outgoing particle with the Boltzmann factor in . The tunneling method has also been used to obtain the Hawking radiation of scalar particles and fermions across an in some black hole geometries [3–6]. Then the derivation of Hawking radiation on the basis of the tunneling mechanism has been actively discussed in the literature (see the reference [7] as a review). A generalized uncertainty principle with a minimal measurable length is attracted a lot of interest to consider phenomenological quantum gravity effects. Though the debates on the quantum of black holes are crucial and long standing, there is no systematic study by quantum gravity so far to lay a solid theoretical foundation for the arguments. Then some different models of quantum gravity are effective ways to understand gravity behaviors at a sufficiently small scale. The combination of relativistic and quantum effects implies that the conventional notion of distance would break down the latest at the Planck scale. The basic argument is that the resolution of small distances requires test particles of short wavelengths and thus of high energies. At such small scales, the gravitational effects by the high energies of test particles would significantly disturb the structure which was tried. Then some quantum gravity theories suggest that there would exist a finite limit to the possible resolution of distances, which would be of order the . From

2 several studies in theory, such a minimal measurable length would be obtained by a generalized uncertainty principle, which is a quantum gravity inspired modification to the conventional Heisenberg uncertainty principle. Then various types of generalized uncertainty principles have been heuristically derived from thought experiments and are often taken as phenomenological models that would accommodate a minimal length [8–16], though there exist attempts to make its formulations more mathematically rigorous [17]. Generalized uncertainty principles have been applied to some different systems and played an important role to consider its corrections by supposed quantum gravity theories [18–35]. Motivated by these discussions of generalized uncertainty principles, we are interested in performing in- depth studies of the quantum features of higher-dimensional black hole solutions. Quantum gravity effects predicted by generalized uncertainty principles on the Hawking radiation have been studied by the tunneling of various particles in a variety of spacetimes including , electrovacuum and with a vast array of scalar fields or effective fluids [36–52]. In this paper, we focus on the quantum tunneling radiation of charged particles coming from higher-dimensional black holes based on a generalized uncertainty principle in the spacetime with compact extra dimensions. Higher-dimensional black hole solutions are actively discussed in the context of string theories and braneworld models. Since our observable world is effectively four dimensional, we can regard higher-dimensional black hole solutions with compactified extra dimensions as candidates of realistic models. We call these Kaluza-Klein black holes. The four-dimensional uniquely describes the general relativistic gravitational field in vacuum with spherical symmetry. However, even if we impose asymptotic flatness to the four- dimensional part of the higher-dimensional spacetime model with Kaluza-Klein structure, the metric is not uniquely determined. A family of five-dimensional squashed Kaluza-Klein black hole solutions [53–57] represent fully five-dimensional black holes near the squashed S3 horizons and asymptote to effective four-dimensional spacetimes with a twisted S1 as an extra dimension at infinity. Then we can regard a series of squashed Kaluza-Klein black hole solutions with a twisted compactified extra dimension as one of realistic higher- dimensional black hole models. Several aspects of squashed Kaluza-Klein black holes have been discussed, for example, thermodynamics [58–60], Hawking radiation [61–72], stabilities [73, 74], gyroscope precession [75, 76], thin disk [77], X-ray reflection spectroscopy [78], deflection [79], strong gravitational lensing [80–84] and black hole shadow [85, 86].

3 In this paper, we investigate the Hawking radiation by the tunneling of charged parti- cles and its quantum gravity effects based on the generalized uncertainty principle in the five-dimensional charged static squashed Kaluza-Klein black hole spacetime. To the best our knowledge, Hawking radiation as a tunneling process in the presence of a generalized uncertainty principle has not been discussed in asymptotically Kaluza-Klein spacetimes. In the present work, we extend the derivation of Hawking radiation on the basis of the tunnel- ing mechanism in four-dimensional black hole spacetimes to the case of the five-dimensional squashed Kaluza-Klein black hole in the context of the generalized uncertainty principle with the minimal measurable length. According to the discussion of the Hawking radiation in general relativity, primordial mini black holes in the very early Universe would have been completely evaporated [87]. However, during the final stages of the Hawking evaporation, the semiclassical approach would be expected to break down due to the dominance of quantum gravity effects. If the black holes would not evaporate completely but turn into the remnants at the end of the evaporation process, such remnants might be a constituent of dark , similar to other candidates including weakly interacting massive particles, the axion, the axino, the neutralino and the gravitino [88–94]. Then it is interesting to investigate how higher-dimensional black holes with compact extra dimensions evolve under a Hawking evaporation process modified by a generalized uncertainty principle. In this paper, we study the evaporation of the squashed Kaluza-Klein black hole by the tunneling of particles in the framework of the generalized uncertainty principle as one of the quantum gravity effects in the Hawking radiation. This paper is organized as follows. In the section II, we review the properties of five- dimensional charged static Kaluza-Klein black hole solutions with squashed horizons. In the section III, we study the tunneling radiation of charged fermions from the squashed Kaluza-Klein black hole with the quantum gravity effects in the context of the generalized uncertainty principle. Then we consider the evaporation process of the black hole and the sparsity of the radiation by the tunneling of fermions in the squashed Kaluza-Klein geometry. In the section IV, we study the Hawking radiation by the tunneling of charged scalar particles from the squashed Kaluza-Klein black hole based on the generalized uncertainty principle. The section V is devoted to summary and discussion.

4 II. SQUASHED KALUZA-KLEIN BLACK HOLES

We consider the charged static Kaluza-Klein black holes with squashed S3 horizons in the five-dimensional Einstein-Maxwell theory [55]. The metric and the Maxwell field are given by

2 2 2 2 K 2 2 2 2 2 2 r∞ 2 ds = −F dt + dρ + ρ K dθ + sin θdφ + 2 (dψ + cos θdφ) , (1) √ F 4K 3ρ ρ A dxµ = ± + − dt, (2) µ 2ρ with (ρ − ρ )(ρ − ρ ) ρ + ρ F = + − ,K2 = 0 , (3) ρ2 ρ 2 where the parameters r∞, ρ± and ρ0 are related as r∞ = 4 (ρ+ + ρ0)(ρ− + ρ0). The coordi- nates run the ranges of −∞ < t < ∞, 0 < ρ < ∞, 0 ≤ θ ≤ π, 0 ≤ φ ≤ 2π and 0 ≤ ψ ≤ 4π. The squashed Kaluza-Klein black hole solution is asymptotically locally flat, i.e., the metric asymptotes to a twisted constant S1 fiber bundle over the four-dimensional -

. The parameter r∞ denotes the size of the compactified extra dimension at infinity. In this paper, to avoid the existence of naked singularities on and outside the horizon, we

restrict ourselves to the ranges of parameters such that ρ+ ≥ ρ− ≥ 0, ρ− + ρ0 > 0. In these parameters, the outer and the inner horizons are located at ρ = ρ+ and ρ = ρ−, respectively. The Komar mass and the charge of the black hole are given by √ πr∞ ρ+ + ρ− 2πr∞ √ ρ+ρ− M = (ρ+ + ρ−) = , |Q| = ρ+ρ− = , (4) G5 2G4 G5 G4

respectively, where the five-dimensional G5 and the four-dimensional

one G4 are related as G5 = 2πr∞G4 [75]. In terms of the mass and the charge, the outer and   p 2 2 the inner horizons are expressed as ρ± = G4 M ± M − Q . The Hawking temperature and the of the black hole associated with the surface gravity of the outer horizon κ and the area of the outer horizon A are obtained as

κ ρ+ − ρ− TKK = = p , (5) 2π 4πρ+ ρ+ (ρ+ + ρ0) p A πρ+ ρ+ (ρ+ + ρ0) SKK = = , (6) 4G5 G4 √ respectively, which satisfy the Smarr-type formula M − 2QA+ 3 = 2TKKSKK, where p . |A+| = 3ρ−/ρ+ 2 is the gauge potential on the outer horizon [58–60]. When ρ− = ρ+,

5 the Hawking temperature TKK vanishes and the metric (1) represents the extremal charged

Kaluza-Klein black hole with the mass M = |Q| [55]. The parameter ρ0 gives the typical

scale of transition from five dimensions to effective four dimensions [75]. In the limit ρ0 → 0, we obtain the metric (1) with K = 1 which represents the four-dimensional Reissner- Nordstr¨omblack hole with a twisted constant S1 fiber. We expect the appearance of the

higher-dimensional corrections, which are related to the parameter ρ0, to the tunneling radiation of particles across the black hole horizon and its Hawking temperature of four- dimensional relativity.

III. TUNNELING OF FERMIONS WITH QUANTUM GRAVITY EFFECTS

A. Derivation of modified Hawking temperature

We consider the Hawking radiation by the tunneling of charged particles from the five- dimensional charged static squashed Kaluza-Klein black hole (1), including the quantum gravity effects predicted by a generalized uncertainty principle. A number of generalized uncertainty principles are proposed as phenomenological models that would describe quan- tum gravity inspired uncertainty in the position of a particle in the vicinity of the black hole horizon [8–17, 26]. In this paper, we focus on the quadratic generalized uncertainty principle with the minimal measurable length in the form [9–15, 23, 24] 1 ∆x∆p ≥ 1 + β (∆p)2 , (7) 2 where ∆x and ∆p respectively represent the uncertainties in the position and the momen- 2  2 tum of a particle, β = β0lp = β0 mp , β0 is a dimensionless positive parameter encoding √ quantum gravity effects on the particle dynamics, lp = ~G4 is the Planck length and p mp = ~/G4 is the Planck mass. When ∆x  lp, one recovers the Heisenberg uncertainty principle as ∆x∆p ≥ 1/2. In standard , ∆x can be made arbitrarily small by letting ∆p grow correspondingly. However, this is no longer the case if the gen- eralized uncertainty relation (7) holds. When ∆p increases with decreasing ∆x, the term β (∆p)2 on the right hand side of the uncertainty relation (7) will eventually grow faster than the left hand side. Then ∆x can no longer be made arbitrarily small and there exists a √ minimal measurable length ∆x0 = lp β0. The generalized uncertainty relation (7) implies a small correction term to the commutation relation in the associative Heisenberg algebra

6 kl  for mirror-symmetric states, i.e., [xi, pj] = i 1 + βg pkpl δij, where xi = x0i is the position jk  operator, pi = p0i 1 + βg p0jp0k is the momentum operator, and x0i and p0i satisfy the canonical commutation relation [x0i, p0j] = iδij [12, 13]. As the consequences of the general- 2 ij 2 ized uncertainty principle (7), using the energy mass shell condition E = g pipj + m , the ˜ 2 ij modified energy operator takes the form E = E (1 − βm − βg pipj), where E = i (∂/∂t) is the energy operator in standard quantum mechanics and m is the mass of a particle [19, 20]. The Hawking radiation emitted by black holes may contain several kinds of particles. The choice of spinors as emitted particles is based on the prediction that the Hawking radiation would favor lower spin, lighter particles and it would be expected to be dominated by fermions [95]. Then, in this section, we consider the tunneling of charged fermions across the horizon of the squashed Kaluza-Klein black hole with the quantum gravity effect. Using the ˜ momentum operator pi and the energy operator E, which leads effectively to a replacement of third order time derivatives by spatial derivatives, we obtain the modified Dirac equation with the Maxwell field up to the first order in the parameter β [19, 20, 22, 25, 37] which describes the behaviors of spinning particles in curved spacetimes as

 t i µ  µ  2 2 jk  i~γ ∂t + i~ γ ∂i + γ Ωµ − eγ Aµ + m 1 − βm + β~ g ∂j∂k Ψ = 0, (8)

µ where ∂µ = ∂/∂x , e is the charge of the particle, Ψ is the modified Dirac field, Ωµ = ab  a a λ ν ν a i ωµ Σab 2, ωµ b = eνeb Γµλ −eb ∂µeν is the spin connection defined by the Christoffel symbol λ a ab  a b Γµν and the vielbein eµ, and Σ = i γ , γ 4 is the Lorentz spinor generator. The a µ µ a µν a b ab vielbeins eµ and the gamma matrices γ = ea γ satisfy the conditions g eµeν = η = µ ν µν diag (−1, 1, 1, 1, 1) and {γ , γ } = 2g , respectively. When Aµ = 0, the equation (8) describes the modified Dirac equation without the Maxwell field [36]. For a spin-1/2 fermion, there are a spin-up and a spin-down states. In this paper, we only consider the tunneling radiation of the fermion with the spin-up state, since the discussion of the spin-down state is fully analogous. For the spin-up state, we employ the following ansatz for the spinor field describing the fermion,   U(xµ)    0   i    µ Ψ =   exp I(x ) , (9) V (xµ) ~   0

7 with the action I and the functions U and V . The corresponding gamma matrices of the squashed Kaluza-Klein spacetime (1) are       0 0 1 0 0 0 1 0 0 0 0 1   √     1  0 0 0 1 F 0 0 0 −1 1 0 0 1 0 t   ρ   θ   γ = √   , γ =   , γ =   , F −1 0 0 0 K 1 0 0 0  ρK 0 1 0 0       0 −1 0 0 0 −1 0 0 1 0 0 0    ˜  0 0 0 −i −2K/r∞ 0 0 iK    ˜  1 0 0 i 0   0 −2K/r∞ −iK 0  φ   ψ   γ =   , γ =   , (10) ρK sin θ 0 −i 0 0   0 iK˜ 2K/r 0     ∞  ˜ i 0 0 0 −iK 0 0 2K/r∞ where K˜ = cos θ /(ρK sin θ) [67, 70]. Substituting the metric (1), the Maxwell field (2), the ansatz (9) and the matrices (10) into the Dirac equation (8), the Wentzel-Kramers-Brillouin approximation to the leading order in ~ yields the equations of motion √ ! 1 F 2K  U √ (∂tI + eAtσ) − σ ∂ρI − σV ∂ψI − m = 0, (11) F K r∞ √ ! 1 F 2K  V √ (∂tI + eAtσ) + σ ∂ρI − σU ∂ψI + m = 0, (12) F K r∞  i  σ ∂ I + (∂ I − cos θ∂ I) = 0, (13) θ sin θ φ ψ where βF β β (∂ I − cos θ∂ I)2 4βK2 σ = 1 − βm2 − (∂ I)2 − (∂ I)2 − φ ψ − (∂ I)2 . (14) 2 ρ 2 2 θ 2 2 2 2 ψ K ρ K ρ K sin θ r∞ It would be difficult to solve the equations (11)-(13) directly. Then, according to the Killing vector fields in the squashed Kaluza-Klein spacetime, ∂/∂t, ∂/∂φ and ∂/∂ψ, we consider the action I in the form

I = −ωt + R(ρ) + Θ(θ) + Jφ + Lψ, (15)

where ω, J and L are the emitted fermion’s energy and the angular momenta in the φ and the ψ direction, respectively. Substituting the action (15) into the equation (13), we have dΘ i  + (J − L cos θ) dθ sin θ " # βF dR2 β dΘ2 β (J − L cos θ)2 4βK2L2 × 1 − βm2 − − − − = 0. (16) 2 2 2 2 2 2 2 K dρ ρ K dθ ρ K sin θ r∞

8 Since the parameter β is a small quantity which represents the quantum gravity effects by the generalized uncertainty principle, the expression inside the square brackets in the equation (16) cannot vanish. Then we obtain

dΘ i + (J − L cos θ) = 0. (17) dθ sin θ

We see that Θ must be a complex function. Substituting the action (15) and the equation (17) into the equations (11) and (12), and cancelling the functions U and V , we have

 2 2 2 2   2 2 2  2 2  2 4K L e At 2 4βK L 2 2ωeAt 4βK L 2 ω 2 − − m 2 + βm − 1 − 2 + βm − 1 − r∞ F r∞ F r∞ F   2 2   2 2 2 2    2 F 4βK L 2 12βK L 2βe At 2 2ωβeAt dR + 2 2 + βm − 1 2 − − βm − 1 − 2 K r∞ r∞ F K dρ 2  2 2 2 2   4 2 3  6 βF 12βK L βe At 2 dR β F dR + 4 2 − + βm − 2 + 6 = 0. (18) K r∞ F dρ K dρ Neglecting the higher-order terms of β in the equation (18) and solving the obtained equation

on the black hole horizon ρ = ρ+ yield the imaginary part of the action, p πρ+ ρ+ (ρ+ + ρ0)(ω − eA+) 2 ImRout = −ImRin = (1 + βΞf ) + O β , (19) ρ+ − ρ− with

3ωm2 ωL2 Ξf = + 2 (ω − eA+) 2ρ+ (ρ− + ρ0)(ω − eA+) 2 ω ρ+ (4ρ+ (ρ+ − 2ρ−) + ρ0 (3ρ+ − 7ρ−)) + eωρ+A+ (4ρ+ρ− + ρ0 (ρ+ + 3ρ−)) + 2 , (20) 2 (ρ+ + ρ0)(ρ+ − ρ−) where the ρ-integral is performed by deforming the contour around the pole at the horizon which lies along the line of integration and gives πi the residue, and Rout and Rin correspond to the outgoing and the ingoing solutions, respectively. Then the tunneling probability amplitude of the charged fermions reads p ! exp (−2ImRout) 4πρ+ ρ+ (ρ+ + ρ0) (1 + βΞf ) Γ ' ' exp − (ω − eA+) , (21) exp (−2ImRin) ρ+ − ρ−

where the contribution from the function Θ cancels out upon dividing the outgoing proba- bility by the ingoing one, since the same solution for Θ is obtained for both the outgoing and the ingoing cases. Thus, by comparing the probability amplitude (21) to the first order

in the energy with the Boltzmann factor Γ = exp (− (ω − eA+)/ T ) in a thermal equilibrium

9 state at the temperature T , we obtain the modified Hawking temperature of the squashed Kaluza-Klein black hole (1) as

2 T = TKK (1 − βΞf ) + O β , (22)

where the temperature TKK and the correction Ξf are given by the equations (5) and (20), respectively. We see that the modified Hawking temperature (22) depends on the energy ω, the mass m, the charge e and the L of the emitted fermion, and is modified by the squashed Kaluza-Klein geometry, the Maxwell field and the generalized uncertainty principle through the parameters ρ0, ρ− and β, respectively. By taking some limits in the equation (22), we obtain the Hawking of some five and four-dimensional black holes. First, when β = 0, the equation (22) coincides with the Hawking temperature of the five-dimensional squashed Kaluza-Klein black hole [58–60]. .  2 p 2 2 Second, when ρ− = 0,L = 0, introducing the new parameters ρ+ = r+ 2 r∞ − r+ . p 2 2 and ρ0 = r∞ − r+ 2, then taking the limit r∞ → ∞, the equation (22) represents the modified Hawking temperature of the five-dimensional Schwarzschild-Tangherlini black hole

obtained by the uncharged fermion tunneling [41, 42]. Lastly, when ρ0 = 0, e = 0,L = 0, the equation (22) represents the modified Hawking temperature of the four-dimensional Reissner-Nordstr¨omblack hole by the uncharged fermion tunneling [38].

B. Evaporation of black holes

We consider the evaporation process of the squashed Kaluza-Klein black hole by the

tunneling of fermions, which is affected by the correction Ξf in the modified Hawking tem- perature T , as one of the quantum gravity effects in the Hawking radiation. The mass of the squashed Kaluza-Klein black hole would decrease due to the radiation. When the black hole mass approaches the order of the Planck mass, the quantum gravity effect could be considered and we may discuss the value of the correction Ξf . First, when Ξf = 0, the effects related to the energy of the fermion, the asymptotically Kaluza-Klein structure and the Maxwell field are canceled. Then the Hawking temperature without the quantum grav-

ity correction TKK appears and results in the complete evaporation. Second, when Ξf < 0, the temperature T is higher than the temperature TKK. Then the black hole accelerates the

evaporation and there is no remnant left. Lastly, when Ξf > 0, the temperature T is lower

10 than the temperature TKK. This implies that the combination of the effects related to the energy of the fermion, the compact extra dimension, the Maxwell field and the generalized uncertainty principle slows down the increase of the temperature due to the radiation. Then the evaporation may cease at the particular mass of the black hole and the black hole may be in a stable balanced state, leading to the remnant mass. Here, we restrict ourselves to the uncharged fermion radiation from the uncharged Kaluza-

Klein black hole (1), i.e., ρ− = 0. In this case, the correction Ξf is positive. Since all particles emitted by the Hawking radiation near the horizon region are effectively massless [2–6], the mass of the emitted fermion does not be taken into account in the following discussion, i.e., m = 0. Further we assume that the fermion has no momentum in the direction of the extra dimension [75], i.e., L = 0. Then, using the saturated form of the uncertainty principle ω ≥ 1/∆x and the uncertainty in the position x for the events near the black hole horizon

∆x ' 2ρ+ [21], we have the Hawking temperature (22) in terms of M, r∞, β0, mp and lp as   p 2 2 mp 10µ + 3 4µ + ν 2 T = 1 − β0   + O β0 , (23) r   2 p 2 2 4π µ 2µ + p4µ2 + ν2 32µ 2µ + 4µ + ν

where µ = M/mp and ν = r∞/lp. We find that, for the large black hole , since the quantum gravity effect is negligible at that scale, the modified temperature (23) asymptotes

to the temperature TKK with ρ− = 0, which monotonically increases with decreasing M/mp

for fixed r∞/lp and is higher than the temperature (23) with β0 6= 0. Using the temperature (23) and the Smarr-type formula, we have the entropy with the quantum gravity effect as   r   p 2 2 p 2 2 10µ + 3 4µ + ν 2 S = 2πµ µ 2µ + 4µ + ν 1 + β0   + O β0 . (24) 32µ2 2µ + p4µ2 + ν2

When the quantum gravity parameter vanishes, i.e., β0 = 0, the equation (24) coincides

with the entropy SKK with ρ− = 0. Using the temperature (23) and the entropy (24), we obtain the modified C = T (∂S/∂T ) as   5 3 2 2 4 2 2 2 p 2 2 C = π 2048µ + 448µ ν − 2β0µν + 1024µ + 96µ ν − 3β0ν 4µ + ν   r   3 2  p 2 2 p 2 2 × 64µ − 10β0µ + 32µ − 3β0 4µ + ν µ 2µ + 4µ + ν h   4 2 2 4 3 2 p 2 2 × 48β0µ 512µ + 136µ ν + 5ν + 256µ + 36µν 4µ + ν  p i−1 −512µ3 128µ4 + 32µ2ν2 + ν4 + 64µ3 + 8µν2 4µ2 + ν2 , (25)

11 1000 0.015 T 0 0.010 C mp -1000 0.005 -2000 0.000 1 2 3 4 5 6 7 8 1 2 3 4 5 6

M/mp M/mp

FIG. 1: Hawking temperatures (23) and heat capacities (25) in various β0 for r∞/lp = 1. β0 = 0

(solid curves), β0 = 15 (Mrm/mp ' 1.36,Mcr/mp ' 2.37, dashed curves) and β0 = 30 (Mrm/mp '

1.93,Mcr/mp ' 3.36, dotted curves).

0.007 6000 0.006 4000 0.005 2000 T 0.004 C 0 mp 0.003 0.0008 -2000 0.002 0.0004 -4000 0.001 0 -6000 2 2.1 2.2 2.3 0.000 -8000 2 4 6 8 10 2 3 4 5 6 7

M/mp M/mp

FIG. 2: Hawking temperatures (23) and heat capacities (25) in various r∞/lp for β0 = 40. r∞ = 0

(Mrm/mp ' 2.24,Mcr/mp ' 3.87, solid curves), r∞/lp = 8 (Mrm/mp ' 2.13,Mcr/mp ' 3.96, dashed curves) and r∞/lp = 25 (Mrm/mp ' 2.02,Mcr/mp ' 4.12, dotted curves). which gives the information on the thermodynamic stability of the present system. We find that the heat capacity (25) with β0 = 0 monotonically increases with decreasing M/mp for

fixed r∞/lp and is negative for M > 0. We show the behaviors of the Hawking temperature

T/mp and the heat capacity C versus M/mp in the figures 1 and 2. From the left panels of the figures 1 and 2, we see that, as the black hole mass decreases, the modified temperature

(23) with β0 6= 0 reaches the local maximum value at the critical mass s 2 p 2 mp 6β0 (3β0 − 5ν ) + 3β0 6β0 (6β0 + 5ν ) Mcr = 2 , (26) 4 2 (3β0 − 2ν )

12 and then decreases to zero at the minimum value of the mass s 2 p 2 mp 2β0 (β0 + 3ν ) + β0 2β0 (2β0 + 3ν ) Mrm = 2 . (27) 4 2 (β0 + 2ν ) The existence of local maximum and vanishing values of the Hawking temperature are related to a transition of the black hole and an evaporation remnant. From the right panels of the figures 1 and 2, we find that there exist three regimes in the heat capacity (25) with

β0 6= 0, i.e., C < 0 for M > Mcr, C ≥ 0 for Mrm ≤ M < Mcr, and C → ∞ for M → Mcr. The discontinuity between negative and positive values of the heat capacity indicates a phase transition from a thermodynamic unstable phase to a stable one. Then we see that, at the

local maximum temperature specified by the mass Mcr, the system undergoes a transition from an unstable negative heat capacity phase to a stable positive heat capacity cooling

down towards a cold extremal configuration with the mass Mrm due to the radiation. At the

Mrm, the black hole no longer radiates since the Hawking temperature (23) vanishes. At the same time, the heat capacity (25) also vanishes, which implies that the black hole may not exchange its energy with the surrounding environment. Then the generalized uncertainty principle prevents the squashed Kaluza-Klein black hole to completely evaporate and results in the thermodynamic stable remnant, similar to the evaporation of the black holes in the noncommutative model and the asymptotically safe gravity [96, 97]. From the left panels of the figures 1 and 2, we see that the mass of the black hole remnant Mrm

increases with increasing β0 for fixed r∞/lp, while it decreases with increasing r∞/lp for

fixed β0. When r∞ = 0, the equation (27) represents the mass of the four-dimensional

Schwarzschild black hole remnant [41]. If the quantum gravity parameter is β0 ' 1 [41] and the size of the extra dimension is r∞ ' 0.1 mm [75], we can estimate that the mass of the −8 squashed Kaluza-Klein black hole remnant is Mrm ' 10 kg, which is of order the Planck mass like the black hole remnants in the minimally geometric deformation model and the quadratic gravity [48, 52]. We consider the impact of the quantum gravity effect onto the sparsity of the Hawking radiation during the evaporation process of the squashed Kaluza-Klein black hole. Sparsity is defined by the average time gap between two successive emissions of the particle over the characteristic timescale of the emission of the individual particle [95, 98]. It has been shown that the Hawking radiation is sparse throughout the whole Hawking evaporation process of the four-dimensional black hole [98]. However, it has also been shown that the

13 sparsity of the Hawking radiation would be modified by some quantum gravity effects in some higher-dimensional black hole spacetimes [30, 31, 33, 97, 99–103]. Further, once the wave effect of the radiation is taken into account, the sparsity becomes a crucial feature that extends the black hole lifetime. Then we consider the sparsity of the Hawking radiation in 2 the form η = λ / Aeff, where λ = 2π /T is the thermal wavelength of the Hawking particle, ˜. Aeff = 27A 4 is the effective area that corresponds to the universal cross section at high frequencies and A˜ is the area of the black hole horizon [30, 31, 33, 102]. If the sparsity is much less than 1, the Hawking radiation is a typical blackbody radiation where its thermal wavelength is much shorter than the size of the emitting body. On the other hand, if the sparsity is much greater than 1, the Hawking radiation is not a continuous emission of particles but a sparse radiation, i.e., most particles are randomly emitted in a discrete manner with pauses in between. Using the temperature (23) and the entropy (24) with the ˜. relation S = A (4G4), we obtain the modified sparsity of the Hawking radiation as  7/2 h   3 11/2 p 2 2 3 2  p 2 2 η = 1048576π µ 2µ + 4µ + ν 27 64µ + 10β0µ + 32µ + 3β0 4µ + ν  2−1 3 2  p 2 2 × 64µ − 10β0µ + 32µ − 3β0 4µ + ν . (28)

We see that the sparsity (28) depends on the black hole mass M, the extra dimension size r∞ and the quantum gravity parameter β0. When r∞ = 0, β0 = 0, the equation (28) coincides with the sparsity of the Hawking radiation from the four-dimensional Schwarzschild black hole, which keeps a constant value 64π3/ 27 during the whole evaporation process [30]. We

show the behaviors of the sparsity η versus M/mp in the figure 3. From the figure 3, we

see that the sparsity increases with decreasing M/mp for fixed r∞/lp and β0, and then the sparsity (28) with β0 6= 0 diverges at M = Mrm. This implies that the average time between the emission of the successive Hawking quanta becomes much larger than the timescales set by the energies of the emitted quanta, so the radiation is very sparse. Then the sparsity of the Hawking radiation is enhanced due to the quantum gravity effect. Thus, similar to the evaporation of the black holes in the noncommutative model and the asymptotically safe gravity [97], the squashed Kaluza-Klein black hole with the generalized uncertainty principle would take an infinite amount of time to radiate a particle at the final stage of the evaporation, and then turn into the remnant when the black hole mass M approaches the mass Mrm. Since the squashed Kaluza-Klein black hole remnant would not radiate and its gravitational interaction would be very weak, it would be difficult to observe the remnants in

14 1000 8000

800 6000 600 η η 4000 400

200 2000

0 0 1 2 3 4 5 6 2.0 2.5 3.0 3.5 4.0

M/mp M/mp

FIG. 3: Sparsities of Hawking radiation (28) in various β0 for r∞/lp = 1 (left panel). β0 = 0 (solid curve), β0 = 15 (Mrm/mp ' 1.36, dashed curve) and β0 = 30 (Mrm/mp ' 1.93, dotted curve).

The same ones in various r∞/lp for β0 = 40 (right panel). r∞ = 0 (Mrm/mp ' 2.24, solid curve), r∞/lp = 8 (Mrm/mp ' 2.13, dashed curve) and r∞/lp = 25 (Mrm/mp ' 2.02, dotted curve). our Universe directly. However, it would be expected that one possible indirect signature of the black hole remnant might be associated with the cosmic background [89, 90].

IV. TUNNELING OF SCALAR PARTICLES WITH QUANTUM GRAVITY EF- FECTS

We consider the Hawking radiation by the tunneling of charged scalar particles in the five-dimensional charged static squashed Kaluza-Klein black hole spacetime (1) with the quantum gravity effects. Taking into account the generalized uncertainty principle (7), we obtain the modified Klein-Gordon equation with the Maxwell field up to the first order in the parameter β [41] as

 tµ g (~∂t + ieAt)(~∂µ + ieAµ) ij 2 2 2 kl  ˜ + g (~∂i + ieAi)(~∂j + ieAj) − m 1 − 2βm + 2β~ g ∂k∂l Ψ = 0, (29) ˜ where Ψ is the modified scalar field. When Aµ = 0, the equation (29) describes the modified Klein-Gordon equation without the Maxwell field [39]. We assume that the wave function of the Klein-Gordon equation (29) takes the form  i  Ψ˜ = exp I˜(xµ) , (30) ~

15 where I˜ is the action of the emitted scalar particle. Substituting the metric (1), the Maxwell field (2) and the ansatz (30) into the Klein-Gordon equation (29), the Wentzel-Kramers- Brillouin approximation to the leading order in ~ yields the equation of motion 1  2 ∂ I˜+ eA +σ ˜ (2βσ˜ − 1) = 0, (31) F t t where 2  ˜ ˜ F  2 1  2 ∂φI − cos θ∂ψI 4K2  2 σ˜ = m2 + ∂ I˜ + ∂ I˜ + + ∂ I˜ . (32) 2 ρ 2 2 θ 2 2 2 2 ψ K ρ K ρ K sin θ r∞ According to the Killing vector fields in the squashed Kaluza-Klein spacetime, we consider the action I˜ in the form

I˜ = −ωt˜ + R˜(ρ) + Θ(˜ θ) + Jφ˜ + Lψ,˜ (33)

whereω, ˜ J˜ and L˜ are the emitted particle’s energy and the angular momenta in the φ and the ψ direction, respectively. Substituting the action (33) into the equation (31), we obtain !2 !4 1 (1 − 4βm2) F dR˜ 2βF 2 dR˜ m2 1 − 2βm2 − (˜ω − eA )2 + − = 0, (34) F t K2 dρ K4 dρ

where we restrict ourselves to the s-wave particles with θ = const, J˜ = 0, L˜ = 0, since the tunneling effect is a quantum one arising within the Planck length near the horizon

region [2–6]. Solving the equation (34) on the black hole horizon ρ = ρ+ for the classically forbidden trajectory yields the imaginary part of the action, p ˜ ˜ πρ+ ρ+ (ρ+ + ρ0) (˜ω − eA+) 2 ImRout = −ImRin = (1 + βΞs) + O β , (35) ρ+ − ρ− with

2 m ρ+ (˜ω − eA+) Ξs = + 2 [˜ω (4ρ+ (ρ+ − 2ρ−) + ρ0 (3ρ+ − 7ρ−)) 2 2 (ρ+ + ρ0)(ρ+ − ρ−)

+ eA+ (2ρ+ (ρ+ + ρ−) + ρ0 (3ρ+ + ρ−))] , (36)

˜ ˜ where Rout and Rin correspond to the outgoing and the ingoing solutions, respectively. Then the tunneling probability amplitude of the charged scalar particles takes the form

 ˜  ! exp −2ImRout p ˜ 4πρ+ ρ+ (ρ+ + ρ0) (1 + βΞs) Γ '   ' exp − (˜ω − eA+) . (37) ˜ ρ+ − ρ− exp −2ImRin

16 Thus, by comparing the probability amplitude (37) to the first order in the energy with ˜  . ˜ the Boltzmann factor Γ = exp − (˜ω − eA+) T in a thermal equilibrium state at the temperature T˜, we obtain the modified Hawking temperature of the squashed Kaluza-Klein black hole (1) as

˜ 2 T = TKK (1 − βΞs) + O β , (38)

where the temperature TKK and the correction Ξs are given by the equations (5) and (36), respectively. We see that the modified Hawking temperature (38) is determined by the mass M and the charge Q of the squashed Kaluza-Klein black hole, the size of the extra dimension

r∞, the quantum gravity parameter β0, and the energyω ˜, the mass m and the charge e of the emitted scalar particle. By taking some limits in the equation (38), we obtain the Hawking temperatures of some five and four-dimensional black holes. First, when β = 0, the equation (38) coincides with the Hawking temperature of the five-dimensional squashed Kaluza-Klein black hole .  2 p 2 2 [58–60]. Second, when ρ− = 0, introducing the new parameters ρ+ = r+ 2 r∞ − r+ . p 2 2 and ρ0 = r∞ − r+ 2, then taking the limit r∞ → ∞, the equation (38) represents the modified Hawking temperature of the five-dimensional Schwarzschild-Tangherlini black hole

obtained by the uncharged scalar particle tunneling [41, 42]. Lastly, when ρ0 = 0, ρ− = 0, the equation (38) represents the modified Hawking temperature of the four-dimensional Schwarzschild black hole by the uncharged scalar particle tunneling [39]. The evaporation process of the squashed Kaluza-Klein black hole by the tunneling ra- ˜ diation of scalar particles depends on the correction Ξs in the modified temperature T . If we consider the quantum tunneling radiation of uncharged massless scalar particles near the horizon region on the basis of the generalized uncertainty principle in the uncharged

Kaluza-Klein black hole spacetime (1), i.e., ρ− = 0, m = 0, the correction Ξs is positive. Then, similar to the Hawking radiation by the tunneling of fermions discussed in the sec- tion III, the equation (38) shows that the generalized uncertainty principle slows down the increase of the Hawking temperature due to the radiation and the evaporation may cease at the particular mass of the black hole, leading to the thermodynamic stable Planck mass remnant.

17 V. SUMMARY AND DISCUSSION

We study the Hawking radiation of charged fermions and charged scalar particles from the five-dimensional charged static squashed Kaluza-Klein black hole on the basis of the tunneling mechanism, including the quantum gravity effects predicted by the quadratic gen- eralized uncertainty principle with the minimal measurable length. We derive the modified Hawking temperature with the corrections related to the energy of the emitted particle, the compact extra dimension, the Maxwell field and the generalized uncertainty principle in the squashed Kaluza-Klein black hole background. Some known Hawking temperatures in five and four-dimensional black hole spacetimes are obtained by taking limits in the modified temperature. We consider the evaporation process of the five-dimensional squashed Kaluza-Klein black hole with the generalized uncertainty principle by the tunneling of particles as one of the quantum gravity effects in the Hawking radiation. We see that the generalized uncertainty principle may prevent the complete evaporation of the squashed Kaluza-Klein black hole. As the black hole mass decreases due to the radiation, the black hole with the negative heat capacity may undergo a phase transition to the one with the positive heat capacity. At the minimum mass of the black hole, both the Hawking temperature and the heat capacity may vanish. This implies that the black hole may not exchange its energy with the surrounding environment. Then the evaporation of the squashed Kaluza-Klein black hole may cease, leading to the thermodynamic stable remnant. Further we consider the sparsity of the Hawking radiation from the squashed Kaluza-Klein black hole in the presence of the quantum gravity effect. We find that the sparsity may increase with decreasing the black hole mass and become infinite at the final stage of the evaporation. This also indicates that the generalized uncertainty principle may stop the Hawking radiation and lead to the black hole remnant. We see that the mass of the squashed Kaluza-Klein black hole remnant increases with increasing the quantum gravity parameter, while decreases with increasing the extra dimension size. If the quantum gravity parameter is of order 1 and the extra dimension size is of order 0.1 mm, the black hole remnant mass is of order the Planck mass. Since such a Planck mass remnant would be a ground state mass of the black hole, the squashed Kaluza-Klein black hole remnant would have no hair. This would be a desirable property as a [48, 52, 88–90, 97]. If the squashed Kaluza-Klein metric would describe the

18 geometry around a and the generalized uncertainty principle considered in this paper would play an important role in the quantum nature of the black hole, the squashed Kaluza-Klein black hole remnants would be a dark matter candidate. Moreover, the variations of the quantum gravity parameter and the extra dimension size provide specific signatures on the quantum features of the squashed Kaluza-Klein black hole solutions with the generalized uncertainty principle which would open the possibility of testing such higher- dimensional models by using future astronomical and astrophysical observations. At present, gravity has been experimentally verified at 9 × 10−5 kg or more [104] and the quantum control has been realized at 5 × 10−11 kg or less [105]. Then it is expected that theoretical and experimental studies of macroscopic quantum phenomena on the Planck mass scale in gravitational fields would bring significant progress in quantum gravity researches. Since the energy of the higher-dimensional black hole remnant would decrease to a detectable range in the , such remnants would be observed in future accelerator experiments [106–110]. If higher-dimensional black holes would be created in an accelerator and we assume that the five-dimensional squashed Kaluza-Klein black hole solutions would describe geometries around such black holes, we expect that our present work would make a contribution to the verification of the Hawking radiation and the extra dimension in asymptotically Kaluza-Klein spacetimes.

Acknowledgments

The author would like to thank Dr. Hideki Ishihara, Dr. Ken-ichi Nakao, Dr. Hirotaka Yoshino and colleagues at the theoretical and gravity group in Osaka City Uni- versity for valuable suggestions and discussions. This work was partly supported by Osaka City University Advanced Mathematical Institute (MEXT Joint Usage/Research Center on Mathematics and Theoretical Physics JPMXP0619217849).

[1] S. W. Hawking, Commun. Math. Phys. 43, 199-220 (1975); 46, 206(E) (1976). [2] M. K. Parikh and F. Wilczek, Phys. Rev. Lett. 85, 5042-5045 (2000). [3] K. Srinivasan and T. Padmanabhan, Phys. Rev. D 60, 024007 (1999). [4] M. Angheben, M. Nadalini, L. Vanzo and S. Zerbini, JHEP 05, 014 (2005).

19 [5] R. Kerner and R. B. Mann, Class. Quant. Grav. 25, 095014 (2008). [6] R. Kerner and R. B. Mann, Phys. Lett. B 665, 277-283 (2008). [7] L. Vanzo, G. Acquaviva and R. Di Criscienzo, Class. Quant. Grav. 28, 183001 (2011). [8] C. A. Mead, Phys. Rev. 135, B849-B862 (1964). [9] D. Amati, M. Ciafaloni and G. Veneziano, Phys. Lett. B 216, 41-47 (1989). [10] K. Konishi, G. Paffuti and P. Provero, Phys. Lett. B 234, 276-284 (1990). [11] M. Maggiore, Phys. Lett. B 304, 65-69 (1993). [12] A. Kempf, G. Mangano and R. B. Mann, Phys. Rev. D 52, 1108-1118 (1995). [13] A. Kempf and G. Mangano, Phys. Rev. D 55, 7909-7920 (1997). [14] F. Scardigli, Phys. Lett. B 452, 39-44 (1999). [15] R. J. Adler and D. I. Santiago, Mod. Phys. Lett. A 14, 1371 (1999). [16] W. Greiner et al., Relativistic Quantum Mechanics, Vol. 2 (Springer, New York, 2000). [17] M. Isi, J. Mureika and P. Nicolini, JHEP 11, 139 (2013). [18] R. J. Adler, P. Chen and D. I. Santiago, Gen. Rel. Grav. 33, 2101-2108 (2001). [19] S. Hossenfelder, M. Bleicher, S. Hofmann, J. Ruppert, S. Scherer and H. Stoecker, Phys. Lett. B 575, 85-99 (2003). [20] U. Harbach, S. Hossenfelder, M. Bleicher and H. Stoecker, Phys. Lett. B 584, 109-113 (2004). [21] A. J. M. Medved and E. C. Vagenas, Phys. Rev. D 70, 124021 (2004). [22] K. Nozari and M. Karami, Mod. Phys. Lett. A 20, 3095-3104 (2005). [23] S. Das and E. C. Vagenas, Phys. Rev. Lett. 101, 221301 (2008). [24] S. Liberati, L. Maccione and T. P. Sotiriou, Phys. Rev. Lett. 109, 151602 (2012). [25] T. L. Antonacci Oakes, R. O. Francisco, J. C. Fabris and J. A. Nogueira, Eur. Phys. J. C 73, 2495 (2013). [26] S. Hossenfelder, Living Rev. Rel. 16, 2 (2013). [27] A. N. Tawfik and A. M. Diab, Int. J. Mod. Phys. D 23, no.12, 1430025 (2014). [28] F. Scardigli and R. Casadio, Eur. Phys. J. C 75, no.9, 425 (2015). [29] Z. W. Feng, S. Z. Yang, H. L. Li and X. T. Zu, Phys. Lett. B 768, 81-85 (2017). [30] A. Alonso-Serrano, M. P. Dabrowski and H. Gohar, Phys. Rev. D 97, no.4, 044029 (2018). [31] Y. C. Ong, JHEP 10, 195 (2018). [32] J. C. S. Neves, Eur. Phys. J. C 80, no.4, 343 (2020). [33] A. Alonso-Serrano, M. P. Dabrowski and H. Gohar, Phys. Rev. D 103, no.2, 026021 (2021).

20 [34] L. Petruzziello and F. Illuminati, Nature Commun. 12, no.1, 4449 (2021). [35] F. Tamburini, F. Feleppa and B. Thid´e,arXiv:2103.13750 [gr-qc]. [36] D. Chen, H. Wu and H. Yang, Adv. High Energy Phys. 2013, 432412 (2013). [37] D. Y. Chen, Q. Q. Jiang, P. Wang and H. Yang, JHEP 11, 176 (2013). [38] D. Chen, H. Wu and H. Yang, JCAP 03, 036 (2014). [39] D. Chen, H. Wu, H. Yang and S. Yang, Int. J. Mod. Phys. A 29, no.26, 1430054 (2014). [40] M. A. Anacleto, F. A. Brito and E. Passos, Phys. Lett. B 749, 181-186 (2015). [41] Z. W. Feng, H. L. Li, X. T. Zu and S. Z. Yang, Eur. Phys. J. C 76, no.4, 212 (2016). [42] Z. Feng, X. Zhu, G. Li, W. Fang and X. Zu, Mod. Phys. Lett. A 31, no.07, 1650022 (2016). [43] H. L. Li and R. Lin, Can. J. Phys. 94, no.11, 1153-1157 (2016). [44]A. Ovg¨un,Int.¨ J. Theor. Phys. 55, no.6, 2919-2927 (2016). [45] I. Sakalli, A. Ovg¨unand¨ K. Jusufi, Astrophys. Space Sci. 361, no.10, 330 (2016). [46] H. L. Li and S. R. Chen, Gen. Rel. Grav. 49, no.10, 128 (2017). [47] R. V. Maluf and J. C. S. Neves, Phys. Rev. D 97, no.10, 104015 (2018). [48] R. Casadio, P. Nicolini and R. da Rocha, Class. Quant. Grav. 35, no.18, 185001 (2018). [49] H. L. Li, Z. W. Feng, S. Z. Yang and X. T. Zu, Eur. Phys. J. C 78, no.9, 768 (2018). [50] D. J. Qi, L. Liu and S. X. Liu, Can. J. Phys. 97, no.9, 1012-1018 (2019). [51] M. Rizwan, M. Z. Ali and A. Ovg¨un,Mod. Phys. Lett. A 34, no.23, 1950184 (2019). [52] I. Kuntz and R. Da Rocha, Eur. Phys. J. C 80, no.5, 478 (2020). [53] P. Dobiasch and D. Maison, Gen. Relativ. Gravit. 14, 231 (1982). [54] G. W. Gibbons and D. L. Wiltshire, Ann. Phys. (N.Y.) 167, 201 (1986); 176, 393(E) (1987). [55] H. Ishihara and K. Matsuno, Prog. Theor. Phys. 116, 417 (2006). [56] C. Stelea, K. Schleich, and D. Witt, Phys. Rev. D 78, 124006 (2008). [57] S. Tomizawa and S. Mizoguchi, Phys. Rev. D 87, 024027 (2013). [58] R. G. Cai, L. M. Cao and N. Ohta, Phys. Lett. B 639, 354-361 (2006). [59] Y. Kurita and H. Ishihara, Class. Quant. Grav. 24, 4525-4532 (2007). [60] Y. Kurita and H. Ishihara, Class. Quant. Grav. 25, 085006 (2008). [61] H. Ishihara and J. Soda, Phys. Rev. D 76, 064022 (2007). [62] S. Chen, B. Wang and R. K. Su, Phys. Rev. D 77, 024039 (2008). [63] S. W. Wei, R. Li, Y. X. Liu and J. R. Ren, Eur. Phys. J. C 65, 281-287 (2010). [64] H. L. Li, EPL 92, no.2, 20003 (2010).

21 [65] D. J. Qi, Int. J. Theor. Phys. 49, 1405-1410 (2010). [66] K. Matsuno and K. Umetsu, Phys. Rev. D 83, 064016 (2011). [67] H. L. Li, Chin. Phys. B 20, 030402 (2011). [68] D. J. Qi, Commun. Theor. Phys. 56, 1171-1174 (2011). [69] P. G. Nedkova and S. S. Yazadjiev, Eur. Phys. J. C 73, no.4, 2377 (2013). [70] M. M. Stetsko, Eur. Phys. J. C 76, 48 (2016). [71] B. Chen, Int. J. Theor. Phys. 55, no.7, 3181-3190 (2016). [72] J. W. Hu, J. H. Wu and X. M. Liu, Int. J. Theor. Phys. 56, no.2, 480-493 (2017). [73] M. Kimura, K. Murata, H. Ishihara, and J. Soda, Phys. Rev. D 77, 064015 (2008); 96, 089902(E) (2017). [74] M. Kimura and T. Tanaka, Classical Quantum Gravity 36, 055005 (2019). [75] K. Matsuno and H. Ishihara, Phys. Rev. D 80, 104037 (2009). [76] M. Azreg-A¨ınou,M. Jamil, and K. Lin, Chin. Phys. C 44, 065101 (2020). [77] S. Chen and J. Jing, Phys. Lett. B 704, 641 (2011). [78] J. Zhu, A. B. Abdikamalov, D. Ayzenberg, M. Azreg-Ainou, C. Bambi, M. Jamil, S. Nam- palliwar, A. Tripathi, and M. Zhou, Eur. Phys. J. C 80, 622 (2020). [79] K. Matsuno, Phys. Rev. D 103, no.4, 044008 (2021). [80] Y. Liu, S. Chen, and J. Jing, Phys. Rev. D 81, 124017 (2010). [81] S. Chen, Y. Liu, and J. Jing, Phys. Rev. D 83, 124019 (2011). [82] J. Sadeghi, A. Banijamali, and H. Vaez, Astrophys. Space Sci. 343, 559 (2013). [83] J. Sadeghi, J. Naji, and H. Vaez, Phys. Lett. B 728, 170 (2014). [84] L. Ji, S. Chen, and J. Jing, J. High Energy Phys. 03 (2014) 089. [85] F. Long, J. Wang, S. Chen, and J. Jing, J. High Energy Phys. 10 (2019) 269. [86] M. Ghasemi-Nodehi, M. Azreg-A¨ınou,K. Jusufi, and M. Jamil, Phys. Rev. D 102, 104032 (2020). [87] M. Sasaki, T. Suyama, T. Tanaka and S. Yokoyama, Class. Quant. Grav. 35, no.6, 063001 (2018). [88] J. H. MacGibbon, Nature 329, 308-309 (1987). [89] P. Chen and R. J. Adler, Nucl. Phys. B Proc. Suppl. 124, 103-106 (2003). [90] P. Chen, New Astron. Rev. 49, 233-239 (2005). [91] F. Scardigli, C. Gruber and P. Chen, Phys. Rev. D 83, 063507 (2011).

22 [92] I. Dalianis and G. Tringas, Phys. Rev. D 100, no.8, 083512 (2019). [93] B. V. Lehmann, C. Johnson, S. Profumo and T. Schwemberger, JCAP 10, 046 (2019). [94] Y. Bai and N. Orlofsky, Phys. Rev. D 101, no.5, 055006 (2020). [95] D. N. Page, Phys. Rev. D 13, 198-206 (1976). [96] Y. S. Myung, Y. W. Kim and Y. J. Park, JHEP 02, 012 (2007). [97] S. Di Gennaro and Y. C. Ong, JCAP 07, 041 (2021). [98] F. Gray, S. Schuster, A. Van-Brunt and M. Visser, Class. Quant. Grav. 33, no.11, 115003 (2016). [99] S. Hod, Phys. Lett. B 756, 133-136 (2016). [100] A. Paul and B. R. Majhi, Int. J. Mod. Phys. A 32, no.16, 1750088 (2017). [101] Y. G. Miao and Z. M. Xu, Phys. Lett. B 772, 542-546 (2017). [102] Z. W. Feng, X. Zhou, S. Q. Zhou and D. D. Feng, Annals Phys. 416, 168144 (2020). [103] S. Schuster, Class. Quant. Grav. 38, no.4, 047002 (2021). [104] T. Westphal, H. Hepach, J. Pfaff and M. Aspelmeyer, Nature 591, no.7849, 225-228 (2021). [105] M. Underwood, D. Mason, D. Lee, H. Xu, L. Jiang, A. B. Shkarin, K. Børkje, S. M. Girvin, and J. G. E. Harris, Phys. Rev. A 92, 061801(R) (2015). [106] B. Koch, M. Bleicher and S. Hossenfelder, JHEP 10, 053 (2005). [107] G. C. Nayak, Phys. Part. Nucl. Lett. 8, 337-341 (2011). [108] J. Mureika, P. Nicolini and E. Spallucci, Phys. Rev. D 85, 106007 (2012). [109] L. Bellagamba, R. Casadio, R. Di Sipio and V. Viventi, Eur. Phys. J. C 72, 1957 (2012). [110] D. M. Gingrich and B. Undseth, Phys. Rev. D 102, no.9, 095020 (2020).

23