#### Abstract

In this paper, using Hamilton-Jacobi method, we address the tunnelling of fermions in a 4-dimensional Schwarzschild spacetime. Based on the generalized uncertainty principle, we introduce the influence of quantum gravity. After solving the equation of motion of the spin-1/2 field, we derive the corrected Hawking temperature. It turns out that the correction depends not only on the black hole’s mass but also on the mass (energy) of emitted fermions. It is of interest that, in our calculation, the quantum gravity correction decelerates the temperature increase during the radiation explicitly. This observation then naturally leads to the remnants in black hole evaporation. Our calculation shows that the residue mass is , where is the Planck mass and is a dimensionless parameter accounting for quantum gravity effects. The evaporation singularity is then avoided.

#### 1. Introduction

Hawking radiation is described as a quantum tunnelling effects of particles at horizons of black holes [1–9]. With the consideration of the background variation in black hole evaporation, Parikh and Wilczek studied the tunnelling behaviors of massless scalar particles [3]. They derived the modified emission spectra for spherically symmetric black holes. The leading corrections to the Hawking temperature are found to be also dependent on the energy of emitted particles. This work was extended to massive and charged scalar particles. The Hawking radiation of general black holes was studied in [4, 5]. For an outgoing massive particle, the equation of motion is different from that of a massless particle. The trajectory of massless particles is a null geodesic, while the massive particle's motion satisfies de Broglie wave and is the phase velocity of outgoing particles. Subsequently, the tunnelling behaviors of fermions were carefully investigated with Hamilton-Jacobi method by Kerner and Mann [10, 11]. The Hawking temperatures were recovered by the fermions tunnelling. The extension of this work to complicated spacetimes is referred to in [12–15]. The noncommutative spacetimes bring a similar consequences in the spirit of black holes thermodynamics. In [16], the tunnelling process in the noncommutative Schwarzschild black holes was researched by Parikh-Wilczek tunnelling method. They derived an interesting result that information might be preserved by a stable black hole remnant.

Various theories of quantum gravity predict the existence of a minimum measurable length [17–21]. This length can be approached from the generalized uncertainty principle (GUP). Through the modified fundamental commutation relation [22] the expression of GUP is derived as , where . is the Planck mass. is a dimensionless parameter. From simple electroweak consideration, it is readily to find an upper limit . and are defined by and , respectively. and satisfy the canonical commutation relations . The modification of the fundamental commutation relation is not unique. Other modifications are referred to in [23–25] and references therein.

[1] These modifications are widely applied to gain some information about the quantum properties of gravity. Black holes are effective modes to explore the effects of quantum gravity. Incorporating effects of quantum gravity into black hole physics by GUP, some interesting implications and results were achieved [26–32]. It is shown in [26] that a small black hole is unstable. Moreover, the constraint for a large black hole comparable to the size of the cavity in connection with the critical mass is needed. In [27], the tunnelling radiation in the noncommutative higher spacetime was discussed with GUP. They found that information may be preserved in a stable black hole remnant. The characteristic size in the absorption process, represented by the black hole irreducible mass, was gotten in [28]. The remnant mass and corrections to the area law and heat capacity were obtained in [29]. In [30], following Parikh-Wilczek tunnelling method, based on GUP, the radiation of massless scalar particles in the Schwarzschild black hole was discussed. The commutation relation between the radial coordinate and the conjugate momentum are modified with GUP. The authors treat the natural cutoffs as a minimal length, a minimal momentum, and a maximal momentum. They addressed the tunnelling rate of black holes. The corrected Hawking temperature was obtained and related to the energy of emitted particles.

[2] The purpose of this paper is to investigate fermions' tunnelling behavior cross the event horizon of a 4-dimensional Schwarzschild black hole, where the effects of quantum gravity are taken into account. We first modify the Dirac equation in curved spacetime to reflect the influence of quantum gravity. The model we adopt is the generalized uncertainty principle. We use the Hamilton-Jacob method to solve the equation of motion of the spinor field. Then, the tunnelling rate and Hawking temperature are calculated. Our results show that the quantum correction to the Hawking temperature is dependent not only on the black hole's mass but also on the mass and energy of emitted fermions. Moreover, the correction slows down the temperature increase during the evaporation. This in turn leads to the remnants in black hole evaporation and prevents the existence of the thermodynamics singularity.

[3] The rest of the paper is organized as follows. In Section 2, taking into account the effects of quantum gravity, we modify the Dirac equation in curved spacetime by GUP and get a generalized Dirac equation. In Section 3, the fermion tunnelling behavior in the Schwarzschild black hole is addressed and the corrected Hawking temperature is derived. Section 4 is devoted to our discussion and conclusion.

#### 2. [4] Generalized Dirac Equation in Curved Spacetime

[5] To take into account the effects of quantum gravity, we adopt the generalized commutation relation in [22] to modify the Dirac equation. In (1), the momentum operators are defined by The square of momentum operators is where in the last step, we only keep the leading order term of . To account for the effects from quantum gravity, the frequency is generalized as [33] with the energy operator . Substituting the mass shell condition , we get the generalized expression of energy [34, 35]: The tunnelling of massless scalar particles in the Schwarzschild black hole was studied in detail and the corrected Hawking temperature was derived in [30]. On the other hand, the Dirac equation with the consequence of GUP in flat spacetime has been investigated in [34].

We start with the Dirac equation in curved spacetime, where is the spin connection defined by the tetrad and ordinary connection The Latin indices live in the flat metric while Greek indices are raised and lowered by the curved metric . The tetrad can be constructed from Back in (6), 's are the Lorentz spinor generators defined by Then, one can construct the 's in curved spacetime as To get the generalized Dirac equation in curved spacetime, we rewrite (6) as where denotes the spatial coordinates. The left-hand side of the equation above is related to the energy. Using the generalized expression of energy (5) and the square of momentum operators (3), only keeping the leading order term of , we get Therefore, the generalized Dirac equation in curved spacetime can be written as This is the equation we are going to solve in the next section.

#### 3. Fermion’s Tunnelling with Effects of Quantum Gravity

In this section, we address the tunnelling behavior of spin-1/2 fermions across the event horizon of the Schwarzschild black hole. Effects of quantum gravity are taken into account. The metric is given by with , and is the black hole's mass. We have set . The event horizon is located at . The fermion motion is determined by the generalized Dirac equation (13). For a spin-1/2 particle, there are two states corresponding, respectively, to spin up and spin down. Following the standard ansatz, to describe the motion semiclassically, we assume the wave function of the spin up state as where , , and are functions of coordinates , , , and and is the action of the emitted fermions. The process of spin down is the same as that of spin up. To solve (13), one should choose appropriate gamma matrices by exploiting (8)–(10). It is straightforward to guess a tetrad for the metric (14): Then, our gamma matrices are given by with and . 's are the Pauli matrices with .

Our task is to find the solutions of (13). First substitute the wave function (15) and the matrices (17) in the generalized Dirac equation (13) and cancel the exponential factor. Since we are working with WKB approximation, the contributions from , , and high orders of are neglected. We finally obtain decoupled four Hamilton-Jacobi equations: To find the relevant solution, since the metric has a time-like killing vector, we perform the separation of variables as follows: where turns out to be the energy of the emitted particle. We insert (22) into (18)–(21) and first focus on the last two equations. They are identical after dividing, respectively, by and and can be rewritten as follows: In the equation above, the value in the square bracket cannot vanish since is a small quantity representing the effects from quantum gravity. Therefore, the expression in the round brackets is zero and yields the solution of . In the previous work, though has a complex solution (other than the trivial one ) and gives rise to a contribution to the imaginary part of the action, it has no contribution to the tunnelling rate. Therefore, (23) is simplified as After cancelling and , (18) and (19) are identical and give rise to with Using (24), we find . Neglecting the higher orders of and solving the above equations at the event horizon yields (In [36–38], the authors argued that the invariance under canonical coordinate transformation requires that the integral to calculate the imaginary part of should be a loop rather than an open one-way integral. However, in our calculation, only the difference between the imaginary parts matters.) In the above equation, . The real part is irrelevant to the tunnelling rate. The sign corresponds to outgoing/ingoing wave. Then, the tunnelling rate [39] of the spin-1/2 fermion crossing the horizon is This is the Boltzmann factor with Hawking temperature: where is the original Hawking temperature. It shows that there is a small correction to the Hawking temperature and the correction value is dependent not only on the black hole's mass but also on the mass and energy of emitted fermions. This property has been obtained in the literature. In [3], energy conservation is enforced by dynamical geometry and the tunnelling rate is found to be . Then, the corrected Hawking temperature is , where the leading correction to the Hawking temperature is related to the energy of emitted particles. To address the effects of quantum gravity, the authors of [30, 40] adopted the modified commutation relation between the radial coordinate and the conjugate momentum. They studied the quantum tunnelling of scalar particles in the Schwarzschild black hole. The tunnelling rate was derived as . Thus, the correction to the Hawking temperature is also related to the black hole's mass and the particle's energy.

It is of interest to note that in (29) the quantum correction slows down the increase of the temperature during the radiation. This correction therefore causes the radiation cease at some particular temperature, leaving the remnant mass. To estimate the residue mass, it is enough to consider massless particles. The temperature stops increasing when Then with the observation and , where is the Planck mass and [41, 42] is a dimensionless parameter marking quantum gravity effects, we can get where we have assumed that the maximal energy of the radiated particle is . This result is consistent with those obtained in [29, 43–45]. Compared with previous results, our calculation explicitly shows how the residue mass of black holes arises due to quantum gravity effects. The singularity of black hole evaporation is then prevented by the quantum gravity correction.

#### 4. Discussion and Conclusion

In this work, we modified the Dirac equation in curved spacetime to include the quantum gravity influence. To fulfill this purpose, we employed the generalized uncertainty principle model. This model is derived from the existence of minimal length which arises when combining quantum and gravity. We calculated the radiation of spin-1/2 particles in the 4-dimensional Schwarzschild spacetime with Hamilton-Jacob method. The tunnelling rate and Hawking temperature were presented.

We found that the quantum gravity correction is related not only to the black hole's mass but also to the mass (energy) of emitted fermions. More interestingly, our result shows that the quantum gravity correction explicitly retards the temperature rising in the black hole evaporation. Therefore, at some point during the evaporation, the quantum correction balances the traditional temperature rising tendency. This leads to the existence of the remnants. We showed that the remnants is , where is the Planck mass and from simple electroweak consideration. Therefore, the classical thermodynamics singularity can be avoided and a residue temperature of black holes exists.

We use the 4-dimensional Schwarzschild metric in this work. It is known that the WKB type of approximation is basically the same as working with a 1+1-dimensional spacetime [46–48]. As a consequence, all large nonextremal black holes look basically the same (like Rindler space). It is also true for the charged and rotating black holes. Therefore, the remnants can be also found in other geometries. It is of interest to employ these geometries in the studies. In our calculation, we keep only the leading order of and . It is expected that higher orders of corrections may give more information in the future work.

#### Acknowledgments

Deyou Chen and Houwen Wu are very grateful for Professor S. Q. Wu and Prof. Professor Wang for their useful discussions. This work is supported in part by the NSFC (Grants nos. 11205125, 11175039, and 11375121), Sichuan Province Science Foundation for Youths (Grant no. 2012JQ0039), and China West Normal University Foundation (Grant no. 11B005).