Total spectral distributions from Hawking radiation

Taking into account the time dependence of the Hawking temperature and finite evaporation time of the black hole, the total spectral distributions of the radiant energy and of the number of particles have been explicitly calculated and compared to their temporary (initial) blackbody counterparts (spectral exitances).


Introduction
One of the most famous features of the Hawking black hole (BH) radiation is its blackbody thermal (Planckian) spectrum [1][2][3][4][5]. However, it is also well known that the BH blackbody thermality is modified in a number of ways [6,7]. Actually, the shape of the BH spectrum is only approximately Planckian, with certain key modifications, especially for small BHs. A direct consequence of the Hawking radiant flux from the BH is its evaporation [1,2], i.e., the mass M of the BH monotonically decreases in time. Moreover, since the BH blackbody temperature T H (the Hawking temperature) depends on the reciprocal of M (see (11)), T H is a growing function of time [1,2]. The author of [6] has presented an elegant and exhaustive discussion of various modifications and limitations of the thermality of the Hawking radiation. According to [6] the Planckian shape of the Hawking spectrum will be modified by at least three distinct physical effects: graybody factors, adiabaticity constraints, and available phase space. From our perspective, somewhat arbitrarily, we could add to the list the following (not quite independent, nor new) modifications of the thermality: the BH radiation lasts for a finite period of time t e (the evaporation time) [1,2], and the total number N of the particles emitted by the BH is finite. Interestingly, it appears that the particles are being emitted from the BH (in a sense) rather rarely [8]. The mentioned (finite) total number N (strictly speaking, its average) of particles (photons) emitted has been calculated quite recently in [9], a e-mail: bobroda@uni.lodz.pl and especially in [10]. The finiteness of N is not quite unexpected, but it is not obvious for gapless particles, e.g. for photons. Moreover, this fact (i.e., the finiteness of N ) could possibly also influence discussion around the BH information paradox [5]. Since we deal with a finite process (t e < +∞), and the Hawking temperature T H is a (non-constant) function of time, what could be physically more representative for the BH evaporation process than the initial or temporary blackbody spectrum expressed by the spectral radiant exitance e ω (see (2)) is the total spectral distribution of energy E ω coming from the total particle production from the BH. Therefore, the aim of the present work is to calculate the total spectral distribution of energy, E ω , and its actinometric counterpart, N ω (the total spectral distribution of the number of particles), both emitted during the whole BH evaporation time t e . In other words, we assume complete evaporation (the final mass M f = 0). This assumption implicitly ignores possible modifications of the evaporation process for the temporary mass of the order of the Planck mass, when unknown quantum gravity effects should probably be taken into account. For simplicity, we also assume that at any instant of time the radiation is given by the Hawking blackbody formula for the Schwarzschild BH. Moreover, we limit our discussion to one species of massless particles, i.e., to photons. Thanks to these simplifications, we are able to express both formulas in closed analytical forms as polylogarithm functions. The both calculated distributions are not only interesting in itself, but also they give rise to the approximate notion of the (time) average(d) temperature T (see (30) and (31)), which we aim to estimate, as a quantity alternative to the (temporary) Hawking temperature T H .
According to standard terminology, the notion radiometric refers to the quantities corresponding to radiant energy, whereas the notion actinometric (photometric or photonic) refers to the quantities corresponding to the number of radiant particles (photons) [11]. In the case of spectral (i.e., ωdependent) quantities, both possibilities are multiplicatively related with the Planck multiplierhω. The simplest relation of this type is where e ω is the spectral radiant exitance, in the case of the blackbody, given by the well-known Planck law [12] (cf. [13,14]) and N ω is the corresponding spectral particle (photon) exitance assuming the standard form ( Here,h is the reduced Planck constant (h/2π ), c the speed of light, k B the Boltzmann constant, ω the angular frequency, T the blackbody temperature, and Li 0 is the polylogarithm function Li of the zeroth order (see the appendix). The notation assumed in the paper is not quite standard-see Chapts. 34 and 36 in [15] for the standard notation. Except where necessary, for clarity, we confine ourselves only to radiometric quantities (actinometric quantities can easily be reproduced according to (10)). Nevertheless, we aim to treat both types of quantities equally and complementing each other.
The total radiant exitance, denoted in our paper by e, is given (for the blackbody) by the famous Stefan-Boltzmann law, where σ is the Stefan-Boltzmann constant, whereas the corresponding total particle (photon) exitance, with ζ(3) ≈ 1.20206. For processes being considered in a finite time interval, we can introduce the notion of the total energy, expressed by integration within this time interval, and the notion of the total number of particles, where A denotes the area of the surface of the source. Analogously, the total spectral distribution of energy, we are most interested in, is defined by whereas the total spectral distribution of the number of particles defined by is directly related to (8) by the relation (cf. (1))

Total number of particles
For the Schwarzschild BH, the spectral radiant exitance e ω , and the spectral particle (photon) exitance N ω is given by the standard thermodynamic formulas (2), and (3), respectively, with T equal the Hawking BH temperature T H , i.e., where G is the Newton gravitational constant. Analogously, the total radiant exitance e, and the total particle (photon) exitance N is given by another couple of standard thermodynamic formulas, (4) and (5), respectively, with T = T H .
The total radiant energy (see (6)) is directly given by the Einstein mass-energy equivalence formula whereas the total number N of particles (photons) emitted by the BH has been determined quite recently in [9,10]. Since, in the next section, for technical reasons, we will need in our calculations the formula expressing the mass flow rate dM/dt in terms of the mass M, involved in the derivation of the evaporation time t e , we rederive it in the present section. As a byproduct of our derivation, we compute the total number N of particles (photons) emitted during the whole BH evaporation process. By virtue of the equivalence formula (12), and of the Stefan-Boltzmann law (4), the mass flow rate where dE/dt can be interpreted as the velocity of evaporation. Inserting to (13) the standard BH horizon surface area formula, and the Hawking BH temperature, Eq. (11), we obtain the differential equation expressing the mass flow rate in terms of M This equation yields the evaporation time where M 0 denotes the initial mass of the BH. As a byproduct of our derivation (15), we can now calculate the total number N of particles (photons) emitted during the whole BH evaporation process. To this end, we will apply a technical trick which consists in using the chain rule to the number flow rate (given by differential form of (7)), and next performing some further formal manipulations (insertions). Thus, first, we get Inserting the mass flow rate (15) into the central part of (17), and next the area formula (14), the total particle exitance (5), the Hawking temperature (11) into the RHS of (17), respectively, after simple rearrangements, we obtain, in agreement with earlier derivations given in [9,10], the differential equation where the Riemann zeta function ζ(3) ≈ 1.202. The solution of (18) yields the total number of particles emitted Thus, the total radiant energy E, and the total number N of particles (photons) emitted, both related to complete evaporation, are determined.

Total spectral distributions
The total spectral distribution E ω of the total energy E emitted during the whole BH evaporation process is determined by the spectral power (given by the differential form of (8)) or, equivalently, using the chain rule as a technical trick, Making use of the mass flow rate equation (15), the area equation (14) and the Planck law (2), we obtain from (21) the (differential) element of the total spectral distribution of energy or in terms of an (auxiliary) dimensionless BH mass Now, we should integrate out the RHS of (24) with respect to the mass x over the finite mass interval [x 0 , 0], corresponding to the total (initial) mass M 0 of the BH, where is a dimensionless frequency. Utilizing the integral identity (36), we obtain our final formula for the total spectral distribution of the total energy, Here, the Riemann zeta function ζ(5) (≈ 1.03693) is the only contribution from 0 (the upper limit of the integral, corresponding to the final mass M f = 0). More precisely, for k = 4 in the sum in the identity (36) (see the Riemann zeta function identity (35)), and there are no other contributions to (36) from 0. In fact, for k = 1, 2, 3, we have (see (35)) lim x→0

Discussion and conclusions
It seems that the best way to discuss the shape of the total spectral distributions E ω and N ω for the BH, given by (26) and (10), is by comparison to analogous temporary (initial) spectral exitances for the corresponding blackbody, i.e., to e ω and N ω , respectively. Differences and similarities between respective spectral functions can be directly observed in Fig. 1. Their maxima, x max (F ω ), have been numerically determined, and they are given in the table in Fig. 1. In particular, we can observe that the values of the arguments for the maxima of the total distributions are (a bit) greater than the values of the arguments for the maxima of their temporary (initial) blackbody counterparts, i.e., x max (E ω ) > x max (e ω ) and x max (N ω ) > x max (N ω ). Since the shapes of the total spectral distributions in their central parts (in particular, we can ignore IR and UV tails, which are a bit unphysical; see [6]) are similar to their blackbody counterparts, making use of the above mentioned shifts of the maxima, we can introduce the approximate notion of the (time) average(d) temperature T . More precisely, actually we deal with two a bit differing average(d) temperatures. One temperature, T E , corresponds to energy, and the second one, T N , corresponds to the number of particles, respectively, and they both are multiplicatively related to T H ; namely, by and respectively. Here, the Hawking temperature T H is the temporary blackbody temperature of the initial BH. Thus, on average, the BH temperature is some 1.3-1.4 times greater than its initial Hawking temperature T H . Since the Hawking radiation of a BH is a blackbody radiation with the temperature T H , which changes in time, and the BH evaporation time t e is finite (see (16)), in this paper, we have focused on the total spectral distributions corresponding to the complete evaporation process, rather than on temporary (initial) blackbody quantities (spectral exitances). The total spectral distribution of energy, E ω , and the total spectral distribution of the number of particles, N ω , have been explicitly calculated and compared to their initial blackbody counterparts.
Integrating by parts, we can verify the following integral identity [11], which is very useful in our calculations: x n Li 0 (e −x )dx = const. − n! n k=0 x n−k (n − k)! Li k+1 (e −x ), valid for any non-negative integer n.