A complementary third law for black hole thermodynamics

There are some examples in the literature, in which despite the fact that the underlying theory or model does not impose a lower bound on the size of black holes, the final temperature under Hawking evaporation is nevertheless finite and nonzero. We show that under some loose conditions, the black hole is necessarily an effective remnant, in the sense that its evaporation time is infinite. That is, the final state that there is nonzero finite temperature despite having no black hole remaining cannot be realized. We discuss the limitations, subtleties, and the implications of this result, which is reminiscent of the third law of black hole thermodynamics, but with the roles of temperature and size interchanged. We therefore refer to our result as the “complementary third law” for black hole thermodynamics.


Introduction: the issue with temperature of black hole remnants
In the usual picture of Hawking evaporation, an asymptotically flat Schwarzschild black hole evaporates completely in finite time, although the time scale is extremely long for a stellar mass black hole. 1 Since the Hawking temperature is inversely proportional to the mass, the black hole becomes hotter as it shrinks. Eventually the energy scale becomes so high that new physics could potentially enter and affect the subsequent evolution. In particular, novel quantum gravity effect may put a stop on Hawking evaporation, thus resulting in a black hole "remnant". The idea of a black hole remnant can be traced back to the work of Aharonov, Casher and Nussinov [1]. It has been suggested that black hole remnants could help to ameliorate the black hole information paradox, though there are arguments 1 A solar mass non-rotating neutral black hole takes 10 67 years to evaporate, which far exceeds the current age of the Universe ∼ 10 10 years. a e-mail: ycong@yzu.edu.cn against the very existence of remnants. See [2] for a comprehensive discussion of the debate and a review of various remnant scenarios.
A popular way of obtaining a black hole remnant is via the generalized uncertainty principle (GUP), which incorporates the effect of gravity into the Heisenberg's uncertainty principle. Since GUP arises from various general considerations involving gravity and quantum mechanics, as well as string theory [3][4][5][6][7][8][9][10][11], it is usually treated as a phenomenological approach to study various properties of quantum gravity. The simplest form of GUP is given by From here onwards, we seth = c = G = k B = 1, unless when explicitly restored for clarity. Note that if α ∼ O(1), as is usually considered in theoretical calculation, then the correction term becomes important at Planck scale. It has been argued that this leads to a correction in the Hawking temperature, resulting in a black hole remnant [12]. In this scenario however, as the evaporation stops at some finite mass, the temperature also stops at a finite, nonzero value, see Fig.  1 below. This is somewhat peculiar: a positive temperature seems to suggest that the black hole continues to emit particles, how then does the evaporation completely stop? A possible interpretation is as follows: since the remnant heat capacity vanishes, there is no thermodynamical interaction with its environment. Therefore, it is thermodynamically inert and behaves like an elementary particle [13]. The finite remnant "temperature" should therefore be interpreted as energy of the remnant (via E = k B T ). If one takes α < 0 in Eq. (1), we would find that there is no lower bound for black hole mass, so in principle the black hole can evaporate completely. However the final temperature is finite and nonzero [14,15], also see Fig. 1.  Such a choice of sign of α may seem unusual, but it is consistent with some quantum gravity models in which physics at the Planck scale "classicalized" and becomes deterministic (as the RHS of the GUP equation goes to zero when p · c is equal to the Planck energy) [13][14][15][16]. Specifically, we have seen that a lattice "spacetime crystal" gives rise to such a GUP [13]. Negative GUP parameter is also required if one accepts that Wick-rotation can be applied to obtain GUPcorrected black hole temperature from a Schwarzschild-like black hole (with higher order terms) [17]. More recently, it has also been shown that non-commutative geometry [18], as well as corpuscular gravity, give rise to negative GUP parameter [19].
In addition, even without GUP, the situation that x p ∼ 0 is compatible with other scenarios that have been proposed in the literature, such as the possibility that h is a dynamical field that flows to zero in high energy limit [20,21]; and with Planck mass fixed in 4-dimensions, Asymptotically Safe Gravity with G → 0 is equivalent tō h → 0 [22] (see also a similar proposal in f (R) gravity [23]). Incidentally, although some of the very first GUP scenarios came from string theory, in which α is naturally positive, it is not clear that negative α is incompatible with string theory. For example, low energy limit of string theory gives rise to charged dilaton black hole, with string coupling being weak near the singularity. Though we have no right to trust this solution near the singularity, as Horowitz put it [24], it is The Hawking temperature of an asymptotically flat Schwarzschild black hole with α = − 1, here in black dash-dotted curve, as a function of time, shows that the temperature tends to a constant value. The mass of the black hole, in red solid curve, tends to zero asymptotically. In order to display both curves in the same diagram, we have multiplied the Hawking temperature by a factor of 120, so that the temperature curve tends to 120T * = 120/4π ≈ 9.549 tempting to speculate whether "contrary to the usual picture of large quantum fluctuations and spacetime foam near the singularity, quantum effects might actually be suppressed." This would then be, at least naively, compatible with the x p ∼ 0 scenario at high energy. This can be seen as a virtue of GUP as a phenomenological tool: by taking different signs of α, it can accommodate different kinds of quantum gravity models. The question is: How does one make sense of a nonzero black hole temperature if the black hole has completely evaporated? A possible interpretation is that this is the temperature of the Hawking radiation at the final moment just before the black hole disappears [14]. This parallels the explanation for the α > 0 case, but instead of interpreting the temperature as energy of the remnant (since there is no radiation), one now takes it to be the temperature of the final emission of radiation (since there is no black hole). An alternative interpretation as "vacuum fluctuation" is discussed further in Sect. 2.
However, in [15], a more detailed study of the evaporation process reveals that the α < 0 GUP-corrected black hole actually takes an infinite amount of time to evaporate completely, so there is no need to resort to the aforementioned interpretation; the black hole simply continues to evaporate indefinitely (indeed its heat capacity is always negative and so it interacts thermodynamically with the environment), with temperature asymptotes to a finite nonzero value. This value is T * = 1/(4π √ |α|) [14,15] (note the typo in [15]).
Thus even though black hole mass is not bounded below, the black hole behaves effectively as a meta-stable remnant at late times. See Fig. 2. This behavior is due to the fact that dM/ dt, though always negative, is not monotonic, and tends to zero as M → 0.
Indeed, the fact that the black hole lifetime is infinite can be shown analytically [15]. The evolution equation is (a minus sign in missing in [15]): dM dt = − M 6 (4|α|π) 4 (2) so as M becomes sufficiently small, we have which leads to where M 0 is the "initial" (small) mass. This leads to a natural question: how generic is this behavior of having an infinite evaporation time when the final temperature is finite and nonzero? Are there any condition required to ensure this? Black holes are known to behave like thermodynamical systems, in particular the third law states that a black hole (of course of nonzero mass) cannot reach zero temperature state in finite number of steps. Here we are claiming a complementary result: a black hole cannot reach a state with nonzero temperature but zero mass. This is the complementary third law of black hole thermodynamics.
We found that this behavior is in fact rather general, and can be stated as Theorem Consider an n-dimensional neutral static black hole spacetime, with areal radius r , and horizon at r = r h . Assume that the Hawking temperature T and the black hole mass M are analytic functions of r h . Suppose dM/ dt = −C AT n , where C > 0 is a constant, and T → T * ∈ (0, ∞) as r h → 0, then r h → 0 only if t → ∞, provided that the k-th derivative M (k) , for k < n − 1, do not all vanish when r h = 0.
We have assumed that there is no problem with convergence of the series expansion. In particular, since T and M are defined only in the domain [0, ∞), all the associated limits and differentiability at 0 are to be understood as being onesided (r h → 0 + ). Note that Hawking temperature of the This result reminded us of the ("Nernst version" of) third law of black hole thermodynamics: zero temperature (extremal) black hole, which is of nonzero size, is unattainable in finite number of steps. Here we have the opposite scenario, zero mass/size 2 black hole is unattainable in 2 Zero mass and zero size are not always interchangeable, see Sect. 4. finite time 3 under Hawking evaporation if the temperature is nonzero. We therefore refer to this theorem as a "complementary third law". See Sect. 4 for more discussions.
Here we assume that the Stefan-Boltzmann law for arbitrary spacetime dimension n 4 holds during the entire evolution [25,26], which of course need not be the case; see, e.g., [27][28][29][30], for a different viewpoint. Note that we assumed that the Hawking evaporation is governed by the simple Stefan-Boltzmann law only. This means, for example, we do not study asymptotically de-Sitter spacetimes, in which Gibbons-Hawking temperature from the cosmological horizon would contribute. Likewise, in an asymptotically locally anti-de Sitter spacetime, the usual reflective boundary condition would complicate the situation (see, however, Sect. 4 for more discussions). A denotes the horizon area, with the constant C incorporating the Boltzmann constant and the greybody factor. The effect of greybody factor, as well as the discreteness of the Hawking radiation (the sparsity [31][32][33]) can be ignored since their effects would result in an even longer evaporation time [15]. In the geometric optic approximation, it is the geometric optic cross section that goes into the Stefan-Boltzmann law, but we also absorb this correction into C, as it will not affect the qualitative behavior of the solution.
In the following, we will first illustrate the theorem with another concrete example, before proving the theorem for the general case. Unlike the GUP-corrected black hole studied in [14] discussed above, the following example is obtained from classical modified gravity (the Hawking radiation itself is of course semi-classical).

Another example: a black hole remnant in massive gravity
Note that the theorem is quite generic: it does not need the underlying theory to be general relativity. Here for explicitness we show an example in the context of dRGT (de Rham-Gabadadze-Tolley) massive gravity [34][35][36][37], in which graviton has nonzero mass. A dyonic black hole solution in this theory was found in [38], with metric coefficient given by where m 0 is related to the physical mass of the black hole. Likewise, q E and q M are related to the electric and magnetic charges, respectively. In addition, m is essentially the graviton mass, c is a positive constant, while c 1 , c 2 and k (the sectional curvature of the horizon) can be positive or negative. The Hawking temperature for this black hole is [38] T = 1 4π We emphasize that these black holes are "physical", in the sense that they have well-behaved thermodynamical properties, as shown in [38]. Of course, dRGT massive gravity has two metric tensors: a dynamical one, g μν , and a fixed background fiducial metric f μν . The first and foremost requirement (and the entire purpose of dRGT theory) is that the theory should be ghost-free. This is not guaranteed for an arbitrary choice of the fiducial metric. In most applications in the context of holography, a degenerate fiducial metric is chosen (as opposed to a Minkowski metric). Such metric was also applied in [38] to obtain the black hole solution above. Such a choice of f μν has been demonstrated to also lead to a ghost-free massive gravity theory [39]. Thus, the theory (and hence the solution) is also physical in the sense that there is no ghost.
In the absence of magnetic charge, the authors showed that by tuning the various parameters such that where E is the electric potential, one could have a solution with Hawking temperature of the form where P = − /(8π) is the pressure term in the extended black hole thermodynamics in an asymptotically locally antide Sitter spacetime. In the limit of vanishing horizon r h → 0, we see that 0 < T = m 2 cc 1 /(4π) < ∞. The authors interpreted this as a fluctuation in the temperature of the background spacetime after the black hole has evaporated. That is to say, there is a trace of the black hole left behind if we look at the energy fluctuation of the vacuum. This is of course tiny because graviton mass m is miniscule, although m of order unity was studied in [38] in the context of holography.
(Of course, large AdS black holes with the usual reflective boundary condition will not evaporate in the first place; but see below.) The simplest way to achieve Eq. (7) is by considering the neutral case and set E = 0. Take also = 0 so that the pressure term vanishes. We choose k = 1 (this is enforced by a topological theorem in general relativity once = 0 [40], but this may not be the case for massive gravity). Since c 2 can be negative, one can set c 2 = −1/ m 2 c 2 . This yields, surprisingly, a constant Hawking temperature regardless of the black hole size, namely T ≡ m 2 cc 1 /(4π), independent of r h .
To further simplify the calculation, we set the numerical values m = c = 1, so c 2 = −1. We remind the readers that our purpose is only to illustrate the aforementioned theorem. It is possible that with this choice of the parameter values the black hole becomes unstable or other problems might arise. 4 The readers are referred to [38] for detailed study of the black hole solutions. (Massive gravity also admits a more conventional black hole remnant that tends to zero temperature with finite size [45].) The physical mass (the mass that appears in the first law of thermodynamics) is [38] which, with our choice of the parameter values, reduces to M = r 2 h /4. In Fig. 3, we set the initial condition r h = 1000, and obtain a plot which shows that the horizon tends to zero size only asymptotically.
The analytic proof is straightforward: with M = r 2 h /4 and T = 1/(4π), we have which yields, withC = 2C/(4π) 4 , where r h (0) is the initial horizon size, and t * is the time at which the horizon has shrunk to ε. Integrating yields Therefore, in order to shrink to zero size, ε → 0, one must have an infinite evaporation time t * → ∞.
The caveat here is that in this particular example, we have assumed that the area that appears in the Stefan-Boltzmann law is the horizon area. However due to the metric function being −g tt = r/2 − 2M/r ∼ r/2 at large r , the asymptotic structure is not flat. We know that in AdS, the effective emitting area of black holes with genus g 1 (k = − 1, 0) is a constant (essentially the square of AdS length scale) and  The black hole also asymptotes to zero size as time goes to infinity independent of the black hole mass [46]. Thus, for our example with unusual asymptotic structure, a similar study along the line of [46] should be carried out to determine its effective emitting area, which might not simply be the horizon area up to some factors. In other words, our example above may not be correct, but we use it to illustrate how the theorem would work if the emission surface is indeed the horizon area.
We can also choose other values for the various parameters so that the black hole is asymptotically AdS-like, and then modify the boundary conditions to allow the black hole to evaporate (see Sect. 4 for details). An example with m = 1 = c, c 2 = − 1, k = l = 1 are shown in Fig. 4. Its evolution time is also infinite, but the rate of Hawking evaporation is different from the previous example without the cosmological constant term (the "pressure" term). The difference will become clear in the next section after the theorem is proved.

Proof of the theorem
We now proceed to prove the theorem. As in the previous example, we will work with r h in place of M, since M is an increasing function of r h . (Though this does not imply that r h = 0 ⇐⇒ M = 0). The dimensionality of spacetime is n 4. We assume that r is the areal radius (if not, change to an appropriate coordinate system under which this is true), thus A ∝ r n−2 h . Absorb the proportional constant into C. It suffices to consider the late stages of the evolution. Assuming analyticity of the Hawking temperature, we can Taylor expand around r h = 0 to obtain where T (0) = T * ∈ (0, ∞) in the statement of the theorem, and prime denotes derivative with respect to r h . Similarly, Taking the derivative yields To lowest order in r h , the differential equation is Since M increases with r , M > 0. In particular, M (0) > 0.
where t 0 1 is the initial condition at a sufficiently late time where the series approximation is valid. Integrating yields It is now clear that if ε → 0, t * must tend to infinity (since n 4).
In particular, in 4-dimensions, we get, This is the case for the GUP-corrected black hole in [15] with negative GUP parameter α, since r h = 2M, the same as the usual Schwarzschild black hole; c.f. Eq. (4). If M (0) = 0, Eq. (16) would lead to, to lowest order of r h , the differential equation This leads to ε r h (0) The evaporation time is clearly infinite for n 5. In 4dimensions, ε is exponential in −t * . Thus, we again obtain an infinite evaporation time. This is the case for the massive gravity black hole example illustrated in Fig. 3.
However, if both M (0) and M (0) vanish, and M (0) = 0, then the integration gives Note that in 4-dimensions, the corresponding result is for some constantK . Therefore t * is now finite as ε → 0: In n 5 the evaporation time is still infinite. Indeed, in general, one can see that as long as the lowest order of nonzero M (k) (0) is k = n − 1, the evaporation time will be finite. This completes the proof. In particular, this implies that in 4-dimensions, the evaporation is infinite if M (0) and M (0) do not both vanish.

The complementary third law: applicability and subtleties
To summarize our findings so far: in this work, we investigated the conditions for a black hole to have "left-over" nonzero and finite temperature at the end of Hawking evaporation, at which point the black hole shrinks to zero size. To our knowledge (and that of the authors of [38]), the massive gravity black hole discussed in Sect. 2 is the only known example in classical modified gravity with such a property. The GUP corrected black hole studied in [14][15][16] provided another example from a quantum gravitational correction. This result parallels the third law of black hole thermodynamics in which (nonzero mass) extremal black hole with zero temperature cannot be attained in finite time, but with the role of mass and temperature reversed, we therefore dubbed it the "complementary third law of black hole thermodynamics".
To be more specific, we found that if the first and the second derivative of M(r h ) do not both vanish at r h = 0, then the evaporation time in 4-dimensions is actually infinite (analogously in higher dimensions), and so the black hole behaves as an effective, meta-stable remnant. This result does not assume the underlying theory to be general relativity. It is simply a consequence of the mathematical properties of the Stefan-Boltzman differential equation. Our proof relies on the assumption that M(r h ) and T (r h ) are both real analytic functions, which are infinitely differentiable at r h = 0. Although most calculations in physics literature make use of Taylor expansion almost ubiquitously, this assumption might be too strong. Perhaps one could relax it and the theorem would still remain true. An alternative proof without the use of series expansion would be welcomed so the issues of convergence can be avoided altogether. Now we shall discuss a few more aspects of the complementary third law.

Charged and Dilaton black holes: nonzero mass vs. nonzero size
It is worth emphasizing that this phenomenon is somewhat opposite to that of the third law of black hole thermodynamics, in which the black hole size remains finite (and nonzero) and T = 0 cannot be achieved; whereas here the temperature remains finite (and nonzero) and r h = 0 cannot be achieved. The main difference is that our theorem only concerns neutral black holes, whereas for the third law, it applies to charged and rotating black holes (the only way to get zero temperature in general relativity 5 ). In the presence of electrical charges and other gauge fields, more analysis would be required to study whether the complementary third law holds, since the evolution under Hawking evaporation would be considerably more complicated. Even for asymptotically flat Reissner-Nordström black holes, Hiscock and Weems showed that there are charge loss and mass loss regimes, so the ratio Q/M is not necessarily monotonic in time (the evolution is governed by coupled differential equations in certain range of the parameters) [50]. However, one special charged black hole solution -the GHS black hole (see below) -is worth a separate mention, since it gives us the opportunity to point out that the complementary third law is really stating that the final state with zero mass but nonzero temperature cannot be attained, not zero size. Usually zero mass also coincides with zero size, but this is not always the case, and therefore the distinction is important.
The charged dilatonic "GHS" (Garfinkle-Horowitz-Strominger) black hole [51][52][53], obtained in a low energy limit of string theory, with metric tensor in a Schwarzschild-like coordinate system {t, r, θ, ϕ}: The Hawking temperature is always T = 1/(8π M) independent of the charge, even in the zero size limit. Note that the horizon stays fixed at r = 2M, but r is not an areal radius, so one has to transform via r 2 → R 2 := r (r − Q 2 /M) to a coordinate {t, R, θ, ϕ} in which the areal radius is R = 4M 2 − 2Q 2 , which goes to zero in the extremal limit |Q| = √ 2M, with nonzero temperature (the temperature is not affected under this change of coordinate). This is very different from extremal Reissner-Nordström black hole which has zero temperature.
Despite dM/ dt being independent of Q, charge loss can occur from spontaneous charge particle emission à la Schwinger, and as shown by Hiscock and Weems [50], this may affect the evolution under Hawking evaporation just like in the Reissner-Nordström case. Thus our theorem does not strictly apply.
Nevertheless, if we ignore these subtleties, 6 and consider the mass loss of GHS black hole to be governed only by Stefan-Boltzmann law, then our theorem does apply. Numerically the results are shown in Fig. 5. Note the peculiarity that R → 0 but M tends to a finite value, as the black hole approaches a null singularity of zero size, which somehow supports the mass. It might be illuminating to also discuss here some nonexamples that do not satisfy the premise of the complementary third law (note that this is not the same as "violating" the law). In view of the popularity of holography, asymptotically locally AdS spacetimes are an important class of solutions. As we mentioned in the Introduction, we assumed that the evolution is governed by the simple Stefan-Boltzmann law. In asymptotically locally AdS spacetimes, our analysis can be carried over by replacing the usual reflective boundary condition with an absorbing boundary condition, which allows large black holes to evaporate (see, e.g., [55][56][57]). In holography, this can be done by coupling the boundary field theory with an auxiliary system ("AUX" [56]), such as another field theory. It was previously found that regardless of their horizon topologies, neutral AdS black holes in such spacetimes take about the same amount of time to evaporate down to the same size of order l, the AdS length scale [58,59]. For positively curved (k = 1) case, the black hole takes about the same amount of time to completely evaporate regardless of its initial mass [58,59]. For flat (k = 0) case, the evaporation time is infinite.
Nevertheless these black holes do not satisfy the premise of the theorems. This is because their temperature, in ndimensional spacetime, is given by [60] which either goes to zero (for k = 0) or diverges (for k = 1) as r h → 0. (For k = −1 case, the black hole tends to a min-imum size as the temperature goes to zero.) [59] In addition, for k ∈ {− 1, 0}, the effective emitting surface area of the black hole is independent of the mass, and is completely fixed by the cosmological constant [46]. The complementary third law therefore does not apply to AdS black holes, at least the simplest ones addressed here.

Comparison with conventional thermodynamics
It came at a surprise when Bardeen, Carter and Hawking [61] discovered that black holes satisfy properties that are analogous to thermodynamics. Subsequent discovery that black holes do radiate when quantum mechanics is taken into account, established that black holes are thermodynamical system. In particular, black holes satisfy a form of third law: zero temperature configuration cannot be reached in a finite number of step. This parallels the "Nernst version" of third law in conventional thermodynamics. What about the complementary third law? Is there an analogous phenomenon in other thermodynamical system? By definition, in classical thermodynamics, the temperature is related to the entropy by 1/T = (∂ S/∂U )| V,N , where U is the internal energy of the system, and the volume V , as well as particle number N , are held fixed. It would seem that nonzero temperature (with nonzero internal energy) always corresponds to nonzero entropy. In the context of black holes, nonzero entropy means nonzero area.
The complementary third law essentially says that black holes with zero entropy yet with nonzero temperature cannot be attained. This is therefore consistent with conventional thermodynamics.
Note that if we consider quantizing the system, the entropy can be zero for small temperature if the first excited state of the system has energy higher than k B T . This does not concern us since black hole thermodynamics is analogous to classical ordinary thermodynamics.

Discussions: some remaining puzzles
The question remains for the case M (0) = M (0) in 4dimensions (and analogously in higher dimensions). Black holes that satisfy this property will evaporate in a finite time, and leave behind a nonzero finite temperature. What is the correct interpretation for such a temperature? Is it the temperature of the last bit of radiation, or an energy fluctuation of the ambient spacetime (now without a black hole -so it would be a kind of "vacuum memory")?
Presumably if the final temperature is very low, the second interpretation (to our knowledge, first proposed in [38]) is plausible, but what if the temperature is "high", say 1 • C (note that the theorem does not constrain the value of T * other that it is nonzero and finite), how can this be a "fluctuation" of the vacuum? It would be good to have an explicit black hole solution of this type for a detailed study, if it exists.
Let us speculate on a possibility: since black holes contain an enormous volume [62,63] (which does not decrease even as the black hole evaporates [64][65][66]), the end state of the evolution could be a baby universe that pinches off from the original universe. Perhaps such a pinch-off leaves a finite temperature signature behind in the parent universe. On the other hand, we could turn this around and say that, if the complementary third law is true generically, then maybe there is no 4-dimensional black hole solution that would satisfy M (0) = M (0), and similarly for higher dimensions.
Another issue concerns the singularity of the black hole. As black hole evaporates, does it leave behind a naked singularity? In the usual Schwarzschild case, since temperature becomes extremely high, one may invoke new physics and hope that the even if singularity wasn't already cured by quantum gravity in the first case, will evaporate away together with the horizon. However, if the final temperature remains mild, it leaves open the possibility that naked singularity may form, thus violating the cosmic censorship conjecture.
To summarize, the complementary third law is more restricted than the standard third law in two ways: firstly, we only study Hawking evaporation, not other physical processes. We restricted our study to the static case in n 4 spacetime dimensions, and only to neutral black holes (though the result might hold in more general cases). Secondly, if we express the black hole mass M as a function of its horizon M(r h ), the complementary third law can be violated if the derivatives M (k) (0) = 0 for k < n − 1, but then the standard third law may also be violated under certain circumstances [67][68][69][70]. In a way, the shortcoming of the complementary third law is a virtue since it gives a clear condition for its violation. A further study into these conditions and their physical interpretations could yield a deeper understanding into black hole thermodynamics.