Hawking radiation as instantons

There have been various interpretations of Hawking radiation proposed based on the perturbative approach, and all have confirmed Hawking's original finding. One major conceptual challenge of Hawking evaporation is the associated black hole information loss paradox, which remains unresolved. A key factor to the issue is the end-stage of the black hole evaporation. Unfortunately by then the evaporation process becomes non-perturbative. Aspired to provide a tool for the eventual solution to this problem, here we introduce a new interpretation of Hawking radiation as the tunneling of instantons. We study instantons of a massless scalar field in Einstein gravity. We consider a complex-valued instanton that connects an initial pure black hole state to a black hole with a scalar field that represents the Hawking radiation at future null infinity, where its action depends only on the areal entropy difference. By comparing it with several independent approaches to Hawking radiation in the perturbative limit, we conclude that Hawking radiation may indeed be described by a family of instantons. Since the instanton approach can describe non-perturbative processes, we hope that our new interpretation and holistic method may shed lights on the information loss problem.


I. INTRODUCTION
One of the most challenging issues in black hole physics is the information loss problem. A black hole evaporates via Hawking radiation [1]. The evaporation depends only on the property of the horizon and hence after the evaporation, it seems that the original information cannot be recovered [2]. The problem has been addressed from various points of view, yet it remains unresolved. One aspect of the problem that may be agreed by most experts is the preservation of unitarity, which is also supported by the notion of AdS/CFT correspondence [3]. Many ideas have been proposed to explain the unitarity maybe preserved during the black hole evaporation (e.g., [4,5]), but no consensus has been reached yet (e.g., [6,7]). A conventional interpretation is the black hole complementarity conjecture [4], according to which information is recovered by Hawking radiation through quantum entanglement [8] without inconsistency. However, some counterexamples and criticisms have recently been reported by several authors [5,9,10], but the AMPS firewall conjecture [5] is also challenged [7].
One difficulty is that there exist only limited tools at hand to deal with the physics of black hole evaporation.
If one only relies on perturbative quantum field theory (e.g., [11]), then the information loss is apparent.
This suggests, under the belief of unitarity, that perturbative methods are insufficient and so non-perturbative methods [12] or new symmetry principles [13,14] may have to be invoked to conclusively address the information loss paradox. Regarding the former approach, we recall that one of the most developed techniques is the Euclidean path-integral approach [15,16].
In many situations, the non-perturbative approach through the Euclidean path-integral makes use of the concept of instantons. We note, however, that Hawking radiation as it was originally derived is a perturbative phenomenon, while an instanton is a non-perturbative process. So at first glance the two notions seem incompatible. Nevertheless, we know that once gravity is taken into account, instantons can describe phenomena associated with vacuum fluctuations that are consistent with that deduced from the perturbative quantum field theory [17,18].
In this paper we will demonstrate that Hawking radiation is describable by Euclidean instantons. To accomplish this, we face two challenges.
1. One issue is the complex-valued nature of instantons. As we will show below, this description requires complex-valued instantons [19]. The Wick-rotation of time requires that all functions should be complexified. The meaning of such instantons has not been completely clarified yet. Nevertheless, complex-valued instantons have been used in quantum cosmology and referred to as fuzzy instantons [20,21], which have led to satisfactory results. Thus we take the view that this approach is legitimate.
2. The other is the regularity of the horizon in Euclidean spacetime. If one chooses the Euclidean time period as the inverse of the Hawking temperature, then the matter distribution between the horizon and the infinity would cause a cusp at the horizon. Unless this cusp is treated by a proper method, one cannot assign a well-defined probability for such an instanton. Fortunately, however, there has been some progress in this issue, and we now have a good reason to believe that the presence of a cusp may be appropriately regularized [22]. Namely, such cusp instantons may be regarded as legitimate, which should contribute to the transition amplitude.
Adopting the above mentioned resolutions, we show that Hawking radiation can be described in terms of Euclidean instantons. Our result implies that there may be additional situations where instantons can play a more useful role than what we previously thought. If Hawking radiation can be described by the tunneling of instantons, it can cover not only the large mass limit, but also the small mass limit when the quantum effects are no more perturbative. This instanton approach may therefore provide a holistic methodology to shed lights on the eventual solution to the information loss paradox. We hope that this work may further promote this research direction [23]. This paper is organized as follows. In SEC. II, we construct scalar field instantons that include Hawking radiation. In SEC. III, we argue that these instantons cover the particle tunneling pictures, including Hartle-Hawking picture and Parikh-Wilczek picture. Finally, in SEC. IV, we summarize our discussions.

II. CONSTRUCTION OF INSTANTONS
The scattering amplitude from an in-state (defined at the past null infinity, say (h in ab , φ in )) to an out-state (defined at the future null infinity, say (h out ab , φ out )) ( Fig. 1) is formally defined by the path-integral, where we sum over all g µν and φ that connects (h in ab , φ in ) and (h out ab , φ out ). Although this integral is not easy to evaluate in practice, extrapolating the flat space quantum field theory, we may assume that it can be evaluated by the analytic continuation to the Euclidean time t = −iτ [24]: which can be further approximated by invoking two assumptions. First, we assume a certain spacetime symmetry and restrict it to the mini-superspace. If the system includes a black hole, then the spherical symmetry is the simplest ansatz for the metric. Second, we approximate the path-integral by using the steepest-descent approximation, where we sum over all on-shell solutions, i.e., the instantons, that connects (h in ab , φ in ) and (h out ab , φ out ). Now let us consider the meaning of on-shell solutions. First, we note that the Euclidean geometry connects the past null infinity, (h in ab , φ in ), and the future null infinity, (h out ab , φ out ). For our purpose the initial condition may include the formation of the black hole, e.g., via a thin-shell collapse or a star collapse. While this adds complications to the problem, it is a classical process and hence as long as we study quantum effects of a black hole, it may be reasonable to start with the maximally extended static black hole (i.e., Schwarzschild) solution ( Fig. 2); and consider instantons on the solution. no scalar field and the future (upper) part has a non-trivial scalar field solution. As we magnify around the symmetric point (Einstein-Rosen bridge), we can approximate that any solution is a superposition of in-going and out-going modes.
We can choose a condition that there are only out-going modes. As we impose the reality condition at (III), the solution is complex-valued for (I) and (II). Since the out-going pulse has energy δM , the green colored region should have mass To cancel out the boundary term at infinity, we define the Euclidean time period as τT = 8πM . Therefore, the horizon becomes a cusp singularity.
The most trivial manifold for this setup is just the Lorentzian solution itself, but in this case the on-shell solutions have no quantum effects. The next simplest solution is the Euclidean Schwarzschild solution that is analytically continued at a certain time slice Σ[t] (say, at t = 0, where t is the conventional Schwarzschild time coordinate). Then we can paste the Euclidean manifold to the past part as well as to the future part (Fig. 3).
Although two Euclidean manifolds are apparently disconnected [25], such linkage between the initial state and the final state is considered consistent [26]. Now we consider solutions with non-trivial field configurations on the manifold in Fig. 3. If we assume that the initial quantum state (in-state) is a vacuum, then at the on-shell level there should be no particles present in the past region. However, there exists a nonzero probability to find a non-vanishing field configuration in the future region (Upper of Fig. 3).
At this juncture we have a couple of questions: 1. Since the solution is in general complex-valued, can one be sure that the solution satisfies the classicality, i.e., reality, condition at future null infinity? 2. Is the probability for each instanton solution consistent with Hawking's result?
Below we answer these questions,. We first define our model in Subsection II A, where we consider Einstein gravity with a massless scalar field. Then in Subsection II B, we answer the first question and clarify that there exists a continuous family of non-trivial instanton solutions that are real at future null infinity. Finally, in Subsection II C, we answer the second question by showing that such classicalized instantons indeed correspond to the Hawking radiation.

A. Model
We consider Einstein gravity with a massless scalar field, where R is the Ricci scalar, K is the Gibbons-Hawking boundary term, and K o is the corresponding boundary term for the periodically identified flat space [27]. Note that there is no potential term and so φ is a free scalar field. Since the Einstein equation gives R = 8π(∇φ) 2 , the volume integration of the on-shell action vanishes.
Therefore for any field configuration, the on-shell Euclidean action becomes Note that the probability of a process mediated by this instanton is P ∼ e −2B , where

B. Two Wick rotations and the classicality of outgoing modes
Let us construct an instanton solution as shown in Fig. 3. Namely, the initial state is a pure black hole, and the final state is a black hole plus a scalar field.
Initially the spacetime is Schwarzschild, where −∞ < t ≤ 0 on the physical side (i.e., the right-hand side of the spacetime in Fig. 3). This is analytically continued to the Euclidean Schwarzschild solution at t = 0, where τ = it + const., and it is periodic with the period τ T = 8πM . As shown in Fig. 3, half of this solution is matched to the initial Lorentzian Schwarzschild. For the sake of convenience, we put τ = it−τ T /2 (equivalently, t = −i(τ + τ T /2)). So its range is −τ T /2 ≤ τ ≤ 0, and it matches to the unphysical side (i.e., the left-hand side in Fig. 3) of the same initial Lorentzian Schwarzschild solution at τ = 0. Now we consider another Euclidean Schwarzschild solution but with a non-trivial scalar field configuration, with the range of the Euclidean time 0 ≤ τ ≤ τ T /2, which is eventually analytically continued to the Lorenztian These two Euclidean solutions are disconnected, but the same formula Eq. (6) can be consistently interpreted as having connected at infinity (r → ∞), where the action S E (background) corresponds to the first pure Schwarzschild instanton and S E (solution) to the second one with a non-trivial field configuration.
To obtain this instanton, we first consider its analytic continued solution in the Lorentzian regime. The solution of the Klein-Gordon equation can be expressed as where R is the conventional Schwarzschild coordinate and r is the tortoise coordinate, and hence −∞ < r < ∞. The corresponding Klein-Gordon equation for each mode becomes There is a potential barrier, but it disappears for R → 2M (r → −∞) as well as for R → ∞ (r → ∞).
Since we are interested in describing Hawking radiation, we focus on the field configurations in the near horizon limit (R → 2M ) in the Euclidean solution. Thus in the regions of the Lorentzian solution close to the matching surface, the waves are also confined to the near horizon limit. In this limit the equation is well approximated by Since the equation is insensitive up to ℓ, to avoid unnecessary complications, we further focus on the ℓ = 0 modes, i.e., spherically symmetric configurations. In addition, since we are interested in the radiation that is emitted to the future null infinity, we consider waves propagating near the black hole (or future) horizon. Therefore, Now we match this to the Euclidean solution through the time slice at t = 0. With t = −i(τ − τ T /2), the above gives This is matched to the Lorentzian solution at τ = 0 to the unphysical, t < 0, region of the solution.
Now since t runs backward in time on the unphysical side of the solution, this describes waves propagating out of the black hole in the past direction. Thus seeing from the physical side, they are the negative energy waves that compensate the energy carried out by waves propagating along and outside the black hole horizon.
Keeping the above discussion in mind, now we impose the physical requirements. First, we demand the asymptotic classicality, namely the out-going mode should contain only real particles. Second, we impose the energy conservation law. Thus if there is an out-going energy flux, then the black hole mass (Misner-Sharp mass) should be smaller than the ADM mass at infinity.
The reality condition implies Note that this implies that the solution in the Euclidean time, Eq. (14), is complex: As discussed in the Introduction, we take the view that complex-valued instantons are legitimate solutions that contribute to the transition amplitude.
To track the energy carried by the radiation, we express the field as where A ω = √ 2ωa out ω . Since the Hamiltonian, H, around the horizon is the number of out-going particles is N ω ∝ |A ω | 2 . Hence, the Misner-Sharp mass at the horizon should be M ′ = M − δM due to the energy conservation.

C. Probabilities
Now we proceed to compute the probability of the process mediated by the instanton. Based on Eq. (5), since there is no contribution from the volume integration, the probability interpretation is the same as that of the thermal thin-shell instanton [31]. The only contributions are from the boundary terms; one from the boundary at infinity and the other at the horizon.
Regarding the boundary term at infinity, since we consider matching the two disconnected instantons at infinity, naturally they cancel each other. Thus the only possible contribution is from the horizon. If there were no cusp, there would be no contribution from the horizon either. But in our case we do have a cusp because of the difference between the period associated with the black hole in the out state and the period fixed at infinity.
Thus we expect to have some contribution from the cusp singularity after some proper regularization [28]. To regularize the cusp, we can apply the same procedure used for thin-shell instantons [22]. We obtain where A and M are the areal radius and mass of the initial black hole, respectively, while the primed ( ′ ) quantities denote those for the final black hole. One can calculate the probability for the same process by using the Hamiltonian approach [32] and can obtain the same result [25]. Hence, we obtain a transition probability that depends only on the area change.
To evaluate the probability, we first consider an emission of a single quantum with energy ω of the order of Thus the exponent is perfectly consistent with Hawking radiation, with the Hawking temperature identified as T H = (8πM ) −1 . Following the usual procedure [33], the total partition function may be evaluated as where Ω is the volume, T is the time, K is determined from the perturbations around the instanton, and As usual this gives the decay rate, Summing over all possible numbers of quanta (ω, 2ω, 3ω, ...), we then readily recover the Planck distribution with T = 1/8πM as we advertised [34].

D. Analogy to the Bogoliubov transformation
Our approach is not the same as the canonical approach based on the Bogoliubov transformation [1,35].
However, we can appreciate their analogy through the following argument.
The scalar field corresponds to particles for ω > 0 and antiparticles for ω < 0. Then before the Wick-rotation, the solution looks like Note that the coordinate time of the left side of the Penrose diagram flows oppositely. In order to identify this with the causal time, one can substitute t → −t, and hence, Now we identify this with the in-state of the Bogoliubov transformation, where a −ω corresponds the particle modes and a ω corresponds the antiparticle modes. In order to avoid the confusion, we redefine the mode functions such that by identifying a −ω =b ω and a ω =b −ω .
After the Wick-rotation, the field becomes From the investigations of the previous subsection, we obtain a linear transformation between (b ω ,b −ω ) and The reality condition after the Wick-rotation requires that b ω = b * ω , i.e., a −ω = e −ωτT a * ω . By using this relation, we can present the equivalent transformation such that frequencies, then the Bogoliubov transformation from the in-going modes to the out-going modes will be with the normalization condition |α ω | 2 − |β ω | 2 = 1. This transformation matrix should include all contributions from different paths and will be approximated by (Fig. 4) where I is the identity matrix (which represents the trivial classical mode propagation) and Q is the contribution from the Euclidean path-integral that is approximated by instantons. Then, the Bogoliubov matrix becomes where N = 1 − e −ωτT is the normalization factor. Finally, we obtain the Hawking-like relation such that where t ω is the transmission coefficient for the modes that propagate in large r region. This is of course not a bona fide quantum field theoretical derivation, since we did not regard all coefficients as real operators. Rather, we heuristically identify the relation between mode coefficients. However, it is worthwhile to see that our results interestingly reproduce that of Hawking radiation.

III. HAWKING RADIATION REVISITED
Our result reveals that instanton approach can indeed recapitulate the Hawking radiation with a thermal spectrum. Traditionally, there have been at least three approaches to Hawking radiation.
First, one can derive Hawking radiation by using the Bogoliubov transformation, which we have checked and demonstrated the consistency between the two approaches in the previous subsection, at least heuristically.
Second, one can derive the evaporation of a black hole by using the renormalized energy-momentum tensor [35]. In order to obtain the renormalized energy-momentum tensor, one needs to invoke a regularization method. The simplest example of exact solution was obtained in 2-dimensions [36], which may be extended to higher dimensional cases using the S-wave approximation [37]. It has been shown that the Euclidean approach may reproduce the approximate forms of the regularized energy-momentum tensors [38].
Third, one can derive it by using the concept of tunneling [39,40]. One can interpret Hawking emission as the tunneling of a particle with energy ω, where the decay rate Γ ∼ e −ω/T and the temperature T can be estimated. In this tunneling approach, one obtains the Boltzmann distribution. After summing over all quanta, one recovers the Planck distribution [34].
In this section we further discuss this tunneling picture [39,40], since the tunneling approach has a direct correspondence to the instanton picture. A. Hartle-Hawking method Although the original derivation by Hawking used the Bogoliubov transformation method [1], soon an explanation based on the pair-creation of particles was given by Hartle and Hawking [39]. They considered the path-integral of a particle that starts from inside the black hole, goes backward in time, reaches outside the event horizon, and then goes forward in time to the asymptotic future infinity (Left of Fig. 5). Although such a trajectory is not allowed in classical physics, it is possible in quantum mechanics.
Hartle and Hawking calculated the decay rate by taking the following logical steps (Right of Fig. 5). - Step 1: For a black hole formed by gravitational collapse, the tunneling process can be described by the analytic continuation of the static black hole solution.
-Step 2: The process to be computed is the tunneling in the left diagram in Fig. 5. The emission probability can be calculated from the propagator, i.e., the path-integral, of the particle. One then invokes the Euclidean analytic continuation to find out an equivalent integration.
-Step 3: In order to calculate the path-integral or the propagator, one can rely on the following three properties.  Fig. 5).
2. The path-integral has the time-translation symmetry for a given r. Therefore, Path (1) in Fig. 5 does not contribute to the propagator.
3. The path-integral has the time-reversal symmetry. Hence, Path (3) in Fig. 5 is equivalent to the absorption probability for a given particle.
-Step 4: Hartle and Hawking chose the equivalent path as that indicated in the Right diagram in Fig. 5.
Therefore, the probability is contributed by two parts: the Boltzmann factor from the Euclidean time evolution (a half period of the Euclidean manifold) and the absorption probability of the particle. However, this approach has several problems. First, the Euclidean analytic continuation cannot form a Euclidean manifold with the signatures (+, +, +, +); rather, it forms (+, +, −, −). Since there is no welldefined manifold for such a signature, it is not easy to impose back-reactions; Hartle and Hawking's description therefore only holds either for emissions with negligible energy ω ≪ M , or for a perfect thermal equilibrium system with negligible black hole mass decrease [41]. One therefore needs find an equivalent steepest-descent method that is applicable to the Euclidean manifold that include processes with non-negligible emission energy ω. Note that it is conceptually easy to see that the Hartle-Hawking approach is analytically equivalent to the instanton approach (Fig. 6), while the latter approach only invokes the Euclidean manifold that satisfies the energy conservation consistently.

B. Parikh-Wilczek method
Parikh and Wilczek [40] developed a more direct description of this process without using the Euclidean analytic continuation. In this approach the Hawking evaporation is treated as the tunnelling of a particle from inside to outside the horizon (Left of Fig. 7). When we consider tunneling of particles with energy ω, the decay rate is where T is the temperature and S is the Lorentzian action of the particle. Such a tunneling between two null geodesics (one inside and the other outside the event horizon) is not allowed classically. Quantum mechanically, however, the action of the particle can acquire an imaginary part such that where τ T = 1/T is the period of the Euclidean time.
It is interesting to compare this with the instanton picture. The dominant contribution of the tunneling between two null geodesics is equivalent to the tunneling between the left side and the right side of the Einstein-Rosen bridge (Middle of Fig. 7). These two geodesics can be connected by the Euclidean manifold (Right of Fig. 7).

IV. DISCUSSION
In this paper, we derived Hawking radiation using the instanton approach. We interpret the Hawking radiation process as the sum over all possible paths of instantons with certain physical constraints such as the energy conservation and the asymptotic classicality. The probability so obtained is consistent with Hawking's result. It is thus reasonable to think that the black hole evaporation can be regarded as the sum of all possible trajectories of a series of instantons. Even though we discussed 4-dimensional massless scalar field cases, it should be very easy to extend it to other dimensions and other fields.
It is important to reemphasize that the instanton method is applicable to much broader areas than Hawking radiation. Not only can it cover perturbative effects (Hawking radiation) but also non-perturbative effects that may even explain, for example, the transition to a trivial geometry. Although in that case their transition probabilities are suppressed, they have nevertheless non-zero probabilities. These non-perturbative effects may in principle be invoked to resolve the information loss problem, especially in the ω ≃ M limit. The present paper reasonably indicates that the existence of non-perturbative channels are indeed very generic for various gravity models, since they manifest themselves even in the case of a free scalar field. If one takes into account all non-perturbative effects that allow for the transition toward a trivial geometry, then the entire wave function is no more semi-classical and one should be able to see new features of the information loss problem beyond what has been considered. This, however, is beyond the scope of the present paper. In addition, non-perturbative contributions should dominate the late stage of the black hole evolution; then, instantons may provide important clues toward the understanding of the late stage black hole evaporation.
In this work there are several directions that need to be improved. We considered instantons only near the horizon; also, we assumed that the metric back-reactions are negligible. These assumptions are reasonable for ω ≪ M , but if ω M , it becomes more complicated. The probability interpretation should be the same since field and the metric at infinity is not clear. This should be confirmed by more detailed numerical investigations.
Although the Euclidean approach has several limitations, we believe that this is the right way of thinking and we hope that this instanton approach would shed lights on the information loss problem.