This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Non-trivial quantum fluctuations in asymptotically non-flat black-hole space-times

Manu Srivastava manu$˙[email protected] Department of Physics, Indian Institute of Technology Bombay, Mumbai 400076, India    S. Shankaranarayanan [email protected] Department of Physics, Indian Institute of Technology Bombay, Mumbai 400076, India
Abstract

Hawking radiation remains a crucial theoretical prediction of semi-classical gravity and is considered one of the critical tests for a model of quantum gravity. However, Hawking’s original derivation used quantum field theory on a fixed background. Efforts have been made to include the space-time fluctuations arising from the quantization of the dynamical degrees of freedom of gravity itself and study the effects on the Hawking particles. Using semi-classical analysis, we study the effects of quantum fluctuations of scalar field stress-tensors in asymptotic non-flat spherically symmetric black-hole space-times. Using two different approaches, we obtain a critical length-scale from the horizon at which gravitational interactions become large, i.e., when the back reaction to the metric due to the scalar field becomes significant. For 4-D Schwarzschild AdS (SAdS) and Schwarzschild de Sitter (SdS), the number of relevant modes for the back-reaction is finite only for a specific range of values of M/LM/L (where MM is the mass of the black-hole, and LL is related to the modulus of the cosmological constant). For SAdS (SdS), the number of relevant modes is infinite for M/L1M/L\sim 1 (0.2<M/L<1/(33)0.2<M/L<1/(3\sqrt{3}). We discuss the implications of these results for the late stages of black-hole evaporation.

I Introduction

Hawking radiation is one of the most striking effects widely agreed to arise from quantum theory and general relativity [1, 2]. Hawking radiation arises from quantizing matter fields in a fixed background space-time with an event horizon. In other words, black-holes constantly radiate particles with energies of the order of Hawking temperature (which is proportional to the surface gravity of the horizon) [3, 4, 5]. Hawking’s original derivation makes two key assumptions: First, the radiation which emerges at late times is strongly redshifted as it climbs out of the black-hole’s gravitational well. In other words, the derivation requires incoming vacuum modes with frequencies far above the Planck scale [6, 7, 8, 9]. Second, the derivation ignores the interactions of the radiation with matter [10].

In the modern viewpoint, general relativity is an effective description valid below a particular cutoff momentum scale [11, 12, 13]. In the effective field theory picture, an effective Lagrangian (valid below some scale E<ΛE<\Lambda)

effg=MPl22R+aRμνRμν+bR2+dRμνλρRμνλρ+eR+cΛ2R3+-\frac{\mathcal{L}_{\mathrm{eff}}}{\sqrt{-g}}=\frac{M_{\rm Pl}^{2}}{2}R+aR_{\mu\nu}R^{\mu\nu}+bR^{2}+dR_{\mu\nu\lambda\rho}R^{\mu\nu\lambda\rho}+e\square R+\frac{c}{\Lambda^{2}}R^{3}+\cdots (1)

accounts for all tree and loop level corrections from particles of masses greater or equal to Λ\Lambda [14]. a,b,c,d,a,b,c,d, and ee are dimensionless constants and MPlM_{\rm Pl} is the Planck mass. However, these do not include the gravitational effects generated from quantum loops of known particles. In principle, these can be computed and are known to be finite and free from renormalization ambiguities [15]. While these effects are, of course, highly suppressed in space-times without horizons, it has been shown that these corrections can be large close to the horizons [16, 17, 18, 19, 20, 21].

Recently, Almheiri et al. [22] used a similar argument and proposed a firewall near the horizon to resolve the inconsistency in black-hole complementarity [10]. At infinity, Hawking radiation is dominated by low angular momentum modes, since high angular momentum modes are trapped behind a potential barrier. More specifically, they showed that an infalling particle would encounter a divergent stress tensor at the stretched horizon [10] and, hence, the need for a firewall. However, in these calculations, it is assumed that the matter fluctuations do not significantly alter the outgoing modes that carry the thermal radiation to distant observers.

In an attempt to address this issue, in this work, we evaluate the quantum fluctuations of the massless scalar field in asymptotically non-flat black-hole space-times. In comparison, the effect of quantum fluctuations of the massless scalar field in an asymptotically flat black-hole space-times has been performed earlier [16, 21], such an analysis has not been performed for asymptotic AdS or dS space-times. In this work, we obtain a length scale from the horizon at which gravitational interactions become large, i.e., when the back reaction to the metric due to the scalar field becomes significant. Within this length scale, Hawking’s assumption [1, 2] of a free field propagating on a given classical background breaks down. In the 4-D Schwarzschild AdS case, we show that the semiclassical analysis to account for the back-reaction breaks down in the limit M/L1M/L\sim 1 where MM is the mass of the black-hole and LL is related to the modulus of the cosmological constant. [This is the limit considered in AdS/CFT correspondence where the large black holes in AdS are dual to the high-temperature phase of the dual field theory [23].] In the 4-D Schwarzschild-de Sitter case, we show the analysis can only be employed for the event horizon and breaks down in the Nariai limit [24].

We calculate the back-reaction of the fluctuations on the horizon using two different approaches. In the first approach, we use the statistical mechanical properties of the perturbations and calculate the change in the mass of the black-hole due to these fluctuations. In the second approach, we calculate the effect of a shockwave (due to an infalling massless particle) on the outgoing Hawking particle [25, 26]. While the two approaches use two different features of the fluctuations, we show that they provide similar results when the semiclassical analysis is valid. We discuss the limitation of our approach near the cosmological horizon in the case of Schwarzschild de Sitter (SdS) and when MLM\sim L for Schwarzschild Anti-de Sitter (SAdS).

In Sec. (II), we discuss the kinematical properties of SAdS and SdS. In Sec. (III), we introduce the setup that we use to calculate the back-reaction of the quantum fluctuations on a 4-D spherically symmetric black-hole space-time. In Secs. (IV) and (V), we obtain the back-reaction using two different approaches. In Sec. (VI), we discuss the limitations of the analysis in M/L1M/L\sim 1 limit. The two Appendices contain the details of the calculations of Approach 2. We use (,+,+,+)(-,+,+,+) signature for the 4-D space-time metric [27]. We use the geometric units G=c=1G=c=1, and to match the dimensions, we retain the Boltzmann constant (kBk_{B}) and Planck’s constant (\hbar).

II Kinematical properties of the black-hole space-times

We consider the following general 4-dimensional spherically symmetric metric:

ds2=f(r)dt2+dr2f(r)+r2dΩ2ds^{2}=-f(r)dt^{2}+\frac{dr^{2}}{f(r)}+r^{2}d\Omega^{2} (2)

with dΩ2d\Omega^{2} being the usual 22-sphere element. In this work, we consider asymptotically non-flat space-times. Specifically, we consider two cases: Schwarzschild Anti-de Sitter (SAdS) and Schwarzschild de Sitter (SdS). In the former we have [23]:

f(r)=12Mr+r2L2=1+r2L2rhr(rh2L2+1)f(r)=1-\frac{2M}{r}+\frac{r^{2}}{L^{2}}=1+\frac{r^{2}}{L^{2}}-\frac{r_{h}}{r}\left(\frac{r_{h}^{2}}{L^{2}}+1\right) (3)

For SdS, we have [24]:

f(r)=12Mrr2L2f(r)=1-\frac{2M}{r}-\frac{r^{2}}{L^{2}} (4)

where MM is the mass of the black-hole and ±L2\pm L^{-2} corresponds to the measure of curvature in the asymptotic limit of dS4dS_{4} and AdS4AdS_{4}. In the limit of LL\to\infty, the above two line-elements reduce to Schwarzschild black-hole (which is a spherically symmetric black-hole in asymptotically flat space-time).

Using the property that the event horizon is a null hypersurface, the horizon (rh)r_{h}) is determined by the condition gμνμNνN=0g^{\mu\nu}\,\partial_{\mu}N\,\partial_{\nu}N=0. For the line-element (2), NN is a function of rr and gμνμNνN=0g^{\mu\nu}\,\partial_{\mu}N\,\partial_{\nu}N=0 leads to f(r)=0f(r)=0. Thus, the location of the event horizon is given by the roots of the equation f(r)=0f(r)=0. SAdS line-element (3) supports only one real positive root and is given by

χ=rhL=23sinh[13sinh1(33ML)]\chi=\frac{r_{h}}{L}=\frac{2}{\sqrt{3}}\sinh\left[\frac{1}{3}\sinh^{-1}\left(3\sqrt{3}\frac{M}{L}\right)\right] (5)

In the limit of LL\to\infty, the above event horizon reduces to rh=2Mr_{h}=2M. The horizon radius is denoted as rhr_{h}.

For 0<ML<1330<\frac{M}{L}<\frac{1}{3\sqrt{3}}, SdS line-element (4) supports two real positive zeroes [24], and are given by:

χb=rbL=23cos(π+ψ3)χc=rcL=23cos(πψ3),\chi_{b}=\frac{r_{b}}{L}=\frac{2}{\sqrt{3}}\cos\left(\frac{\pi+\psi}{3}\right)\quad\quad\quad\chi_{c}=\frac{r_{c}}{L}=\frac{2}{\sqrt{3}}\cos\left(\frac{\pi-\psi}{3}\right)\,, (6)

where

ψ=cos1(33ML)\psi=\cos^{-1}\left(3\sqrt{3}\frac{M}{L}\right) (7)

rbr_{b} corresponds to the event horizon of the black-hole, while rcr_{c} refers to the cosmological horizon. In the range 0<ML<1330<\frac{M}{L}<\frac{1}{3\sqrt{3}}, χb\chi_{b} takes the values in (0,13\frac{1}{\sqrt{3}}).

Since we want to obtain the length-scale in which the back-reaction of the fluctuations on the horizon can not be ignored, evaluating the quantities in terms of χ,χc,χb\chi,\chi_{c},\chi_{b} (instead of other parameters like M/LM/L) enables exact calculations. In the next section, we will discuss this aspect.

To compare and contrast our results with asymptotically flat space-time black-holes, we define the ratios of quantities for the SAdS (SdS) black-holes to the quantities for the Schwarzschild black-hole. Physically, this corresponds to the two black-holes with identical mass. We denote these quantities with a subscript RR (for the ratio) and a superscript SAdS (SdS). For instance, the dimensionless surface gravity for SAdS is given by:

κRSAdS=κSAdSκS=(1+χ2)(3χ2+1)\kappa^{SAdS}_{R}=\frac{\kappa^{SAdS}}{\kappa^{S}}=(1+\chi^{2})(3\chi^{2}+1) (8)

where the surface gravity for the black-hole line-element (2) is given by:

κ=12(df(r)dr)r=rh\kappa=\frac{1}{2}\left(\frac{df(r)}{dr}\right)_{r=r_{h}} (9)

and the corresponding Hawking temperature TT is

T=κ2πkBT=\frac{\hbar\kappa}{2\pi k_{B}} (10)
Refer to caption
(a) SAdS
Refer to caption
(b) SdS
Figure 1: Plot of ration of surface gravity with χ\chi for SAdS and SdS.

For the SdS case, there exists 2 horizons and the effective surface gravity is given by [24]:

κeff1=|κb1|+|κc1|\kappa_{eff}^{-1}=|\kappa^{-1}_{b}|+|\kappa^{-1}_{c}| (11)

where κb\kappa_{b} and κc\kappa_{c} are calculated using (9) at the black-hole horizon and the cosmological horizon, respectively. In Fig 1, we have plotted κRSAdS\kappa^{SAdS}_{R} and κeffSdS\kappa^{SdS}_{\rm eff} as a functions of χ\chi. Note that as χ\chi increases, the black-hole horizon size increases.

Our aim is to obtain a length scale at which the back-reaction due to the quantum scalar field becomes significant. In order to do that we need to obtain a line-element near the horizon for any spherically symmetric space-time. To do that, we make the coordinate transformation (t,r)(t,γ)(t,r)~{}\to~{}(t,\gamma), given by:

γ1κf,dγ=12κdrffdr,\gamma\equiv\frac{1}{\kappa}\sqrt{f},\qquad d\gamma=\frac{1}{2\kappa}\frac{d_{r}f}{\sqrt{f}}\,dr\,, (12)

where κ\kappa is given by (9). In the above local transformation, the horizon (rhr_{h}) is at γ=0\gamma=0. The line-element (2) becomes [28]:

ds2=κ2γ2dt2+4κ2(drf)2dγ2+r2dΩ2,ds^{2}=-\kappa^{2}\gamma^{2}dt^{2}+4\,\frac{\kappa^{2}}{(d_{r}f)^{2}}d\gamma^{2}+r^{2}d\Omega^{2}\,, (13)

and hence, near the horizon, we have

ds2κ2γ2dt2+dγ2+r02dΩ2.ds^{2}\to-\kappa^{2}\gamma^{2}dt^{2}+d\gamma^{2}+r_{0}^{2}\,d\Omega^{2}\,. (14)

For space-times with non-degenerate horizons (like SAdS and SdS for which κ0\kappa\neq 0), f(r)f(r) can be expanded around rhr_{h} as f(r)=f(rh)(rrh)f(r)=f^{\prime}(r_{h})~{}(r-r_{h}). Near the horizon, the invariant distance (12) becomes:

γ2f(rh)rrh=(2(rrh)κ)1/2\gamma\sim\frac{2}{\sqrt{f^{\prime}(r_{h})}}\sqrt{r-r_{h}}=\left(\frac{2(r-r_{h})}{\kappa}\right)^{1/2} (15)

The line-element (14) describes the near horizon geometry for a generic spherically symmetric black-hole space-time with any asymptotic structure. This relation is crucial for one of the approaches to obtain the critical invariant distance [16]. We will discuss more on this in Sec. (V).

III Setup

In this section, we provide the setup to evaluate the invariant distance near the horizon where the back-reaction due to the quantum field becomes significant. The quantum field we consider is a massless, minimally coupled scalar field propagating in the spherically symmetric background (2). The evolution equation of the scalar field is given by:

Φ=1gμ(ggμννΦ)=0\square\Phi=\frac{1}{\sqrt{-g}}\partial_{\mu}\left(\sqrt{-g}g^{\mu\nu}\partial_{\nu}\Phi\right)=0 (16)

The rotational symmetry of the line-element (2), allows us to decompose into the normal modes ul(r,ω)u_{l}(r,\omega) (where, ω=E/\omega=E/\hbar) of the field Φ\Phi as follows:

ϕlm=ul(r,ω)rYlm(θ,ϕ)eiωt\phi_{lm}=\frac{u_{l}(r,\omega)}{r}Y_{lm}(\theta,\phi)e^{-i\omega t} (17)

Putting this back in the wave equation we get:

(d2dr2+ω2V(r))ul(r,ω)=0\left(\frac{d^{2}}{dr_{*}^{2}}+\omega^{2}-V(r)\right)u_{l}(r,\omega)=0 (18)

where

r=drf(r)r_{*}=\int\frac{dr}{f(r)} (19)

denotes the tortoise coordinate, and the effective 1-D potential (often referred to as Regge-Wheeler potential) V(r)V(r) is:

V(r)=f(r)[1rdf(r)dr+l(l+1)r2]V(r)=f(r)\left[\frac{1}{r}\,\frac{df(r)}{dr}+\frac{l(l+1)}{r^{2}}\right] (20)
Refer to caption
(a) SAdS
Refer to caption
(b) SdS
Figure 2: Plot of effective 1-D potential (4M2V)(4M^{2}V) for SAdS and SdS black-holes (M/L=0.1{M}/{L}=0.1).

The potential V(r)V(r) has maxima near event horizon for SAdS and SdS cases and can be treated as an effective 1-D quantum mechanical problem of a particle in a potential barrier. Fig. 2 is the plot of dimensionless potential (4M2V)(4M^{2}V) with dimensionless radial coordinate (r/(2M))(r/(2M)). We have set M/L=0.1{M}/{L}=0.1. For both the cases, near the event horizon, the barrier height increases with increasing ll. Thus, large ll modes get trapped near the horizon.

Our analysis rests upon the fact that there is a peak (positive barrier) in the potential profile that traps the high angular momentum modes close to the horizon. For SAdS black-hole, there is no maximum for χ1\chi\sim 1 (also M/L1M/L\sim 1). In SdS, there is no maximum near the cosmological horizon (rc)(r_{c}) for any value of M/LM/L. We relegate the implication of this condition to Sec. VI. We assume that the effective 1-D potential V(r)V(r) has a barrier for the rest of the analysis.

The modes that contribute to the energy or mass content of the black-hole (as seen from infinity) are the scattering modes of (18). The condition for such modes is:

ω2>V(r)\omega^{2}>V(r) (21)

As mentioned earlier, the black-holes constantly radiate particles with energies of the order of Hawking temperature. Hence, we take ωκ\omega\sim\kappa. Using the relation (9), we have:

(f(rh)2)2>f(r)[1rdf(r)dr+l(l+1)r2].\left(\frac{f^{\prime}(r_{h})}{2}\right)^{2}>f(r)\left[\frac{1}{r}\,\frac{df(r)}{dr}+\frac{l(l+1)}{r^{2}}\right]\,. (22)

Using the near horizon approximation, we have

l2<rh2γ2leading tolmax2rh2γ2.l^{2}<\frac{r_{h}^{2}}{\gamma^{2}}\quad\mbox{leading to}\quad l_{\rm max}^{2}\sim\frac{r_{h}^{2}}{\gamma^{2}}\,. (23)

Thus, the total number of relevant modes are:

N(γ)=l=0lmax(2l+1)lmax2=rh2γ2N(\gamma)=\sum_{l=0}^{l_{\rm max}}(2l+1)\sim l_{\rm max}^{2}=\frac{r_{h}^{2}}{\gamma^{2}} (24)

As mentioned above, our analysis rests upon the fact that the effective 1-D potential has a maximum that traps the high angular momentum modes close to the horizon. In the rest of this work, using two different approaches, we evaluate the length scale at which the back-reaction of the quantum modes on the background metric becomes significant [16]. We refer to the two approaches as Statistical mechanical (approach 1) and Quantum Field Theory (approach 2).

In approach 1, we obtain the critical length by evaluating the back-reaction of the metric fluctuations and several thermodynamical quantities. In approach 2, we calculate the effect of a shockwave (due to an infalling particle) on the outgoing Hawking particles. The two approaches differ in terms of how the critical length is evaluated. Also, the back-reaction of the fluctuations and the critical length at which the Hawking particles change state vary.

IV Approach 1: Statistical Mechanics

This section evaluates the physical quantities from statistical mechanics and shows that the quantum scalar field causes an fluctuation in the horizon radius. We show that the thermal modes of the field back-react on the metric contributing to the energy of the black-hole. The statistical fluctuations lead to the smearing of the horizon. We then calculate the invariant distance at which the back reaction due to the field on the metric becomes significant. We do this by including the effects of the statistical fluctuation in the wave equation (18).

Using N(γ)N(\gamma) obtained in Eq. (24), we can calculate the average energy contributed by these modes, which is given by [29]:

EkBTN(γ)=κ2π(rhγ)2\left\langle E\right\rangle\sim k_{B}\,T\,N(\gamma)=\frac{\hbar\kappa}{2\pi}\left(\frac{r_{h}}{\gamma}\right)^{2} (25)

Similarly, the statistical entropy of these thermal modes is [29]:

SkBN(γ)=kB(rhγ)2S\sim k_{B}\,N(\gamma)\,=\,k_{B}\left(\frac{r_{h}}{\gamma}\right)^{2} (26)

Note that as we move closer to the horizon (γ0\gamma\to 0), the average energy and the scalar field’s statistical entropy diverge. The fluctuations in the energy of these modes are given by:

Δ(E)TkBET(κ2π)N(γ)=(κ2π)rhγ\Delta(E)\sim T\sqrt{k_{B}\frac{\partial\left\langle E\right\rangle}{\partial T}}\sim\left(\frac{\hbar\kappa}{2\pi}\right)\sqrt{N(\gamma)}=\left(\frac{\hbar\kappa}{2\pi}\right)\frac{r_{h}}{\gamma} (27)

The first part of the above equation is a result from statistical mechanics. If the partition function is Zexp(βE)Z\sim\sum{\exp(-\beta E)}, then average energy is given by E=Eexp(βE)Z=lnZβ\langle E\rangle=\frac{\sum{E\exp{(-\beta E)}}}{Z}=-\frac{\partial lnZ}{\partial\beta} and the average of energy squared is given by E2=E2exp(βE)Z=2Zβ21Z\langle E^{2}\rangle=\frac{\sum{E^{2}\exp{(-\beta E)}}}{Z}=\frac{\partial^{2}Z}{\partial\beta^{2}}\frac{1}{Z}. If we differentiate the average energy expression again with respect to β\beta, we get Eβ=1Z2(Zβ)21Z2Zβ2=E2E2=(ΔE)2\frac{\partial\langle E\rangle}{\partial\beta}=\frac{1}{Z^{2}}\left(\frac{\partial Z}{\partial\beta}\right)^{2}-\frac{1}{Z}\frac{\partial^{2}Z}{\partial\beta^{2}}=\langle E\rangle^{2}-\langle E^{2}\rangle=-(\Delta E)^{2}. Converting the β\beta derivative to temperature (T) derivative gives the first part of Eq 27. The subsequent parts of Eq. (27) are substitutions of different quantities from equations (24), (25), (10). We want to make a couple of remarks regarding the above result: First, the above result is generic for any spherically symmetric space-time. Second, like energy and entropy, the energy fluctuations also diverge close to the horizon. The fluctuation in energy leads to a fluctuation in the location of the black-hole horizon.

In the rest of this section, we use the above expression to obtain the fluctuation in the horizon and calculate the length-scale at which the back-reaction becomes relevant. We then apply the results specifically to SAdS and SdS cases.

IV.1 Horizon fluctuations

Let us consider a point r0r_{0} outside the horizon such that,

δ(r0)=r0rh>0andΔδ(r0)=Δrh(r0).\delta\left(r_{0}\right)=r_{0}-r_{h}>0\quad\mbox{and}\quad\Delta\delta\left(r_{0}\right)=\Delta r_{h}(r_{0})\,. (28)

The fluctuation in the black-hole energy leads to the fluctuation in the distance from the horizon (Δδ(r0)\Delta\delta\left(r_{0}\right)). (Note that, throughout this subsection, we deal with the absolute values of Δδ\Delta\delta and Δrh\Delta r_{h} and hence have dropped the ve-ve sign in the right hand equation in (28).) Since δ(r0)\delta(r_{0}) depends on the energy of the field modes, we have

Δδ(r0)=Δrh(r0)=(Δrh)(ΔE)ΔE\Delta\delta\left(r_{0}\right)=\Delta r_{h}(r_{0})=\frac{\partial(\Delta r_{h})}{\partial(\Delta E)}\Delta E (29)

In the case of Schwarzschild black-hole, ΔEΔM\Delta E\sim\Delta M (where MM is the mass of the black-hole), we have (Δrh)/(ΔE)2{\partial(\Delta r_{h})}/{\partial(\Delta E)}\sim 2. [The quantum scalar field modes contribute to the total energy leading to a change in the energy. For an asymptotic observer, the system’s energy content without the scalar field modes is the black-hole mass MM. Including the scalar field modes, the system’s total energy for the observer at infinity is M+Field’s energy contentM+\mbox{Field's energy content}. As a consequence of the Birkhoff theorem, this is the new ADM mass of the black-hole space-time. Hence, we can interpret the fluctuation in the energy content of the field as the fluctuation in the ADM mass of the black-hole. Hence, the relation ΔEΔM\Delta E\sim\Delta M.]

For r0r_{0} to be certainly outside the horizon (including the contribution from the field modes), the following condition must be satisfied [16]:

Δδ(r0)<δ(r0).\Delta\delta\left(r_{0}\right)<\delta\left(r_{0}\right)\,. (30)

To be certain that δ(r0)\delta\left(r_{0}\right) is positive and observer independent, it necessary to relate it to the invariant distance. Using (15), we have:

δ=γ2κ/2\delta=\gamma^{2}\kappa/2 (31)

Substituting Eqs. (31) and (27) in Eq. (29), we get:

γ3>rh(Δrh)(ΔE)\gamma^{3}>\hbar r_{h}\,\frac{\partial(\Delta r_{h})}{\partial(\Delta E)} (32)

Therefore, the minimum length up to which the fluctuations smear the horizon is:

γlim=[(rh)(Δrh)(ΔE)]1/3.\gamma_{{}_{\rm lim}}=\left[(\hbar r_{h})\frac{\partial(\Delta r_{h})}{\partial(\Delta E)}\right]^{1/3}\,. (33)

A closed-form expression can be obtained only for the Schwarzschild black-hole since Δrh\Delta r_{h} and ΔE\Delta E are related by a simple relation. For SAdS or SdS black-hole, Δrh\Delta r_{h} is related to ΔE\Delta E via transcendental functions and can only be obtained (perturbatively) for a narrow range of parameters. To overcome this, we now use the effective metric approach to evaluate the critical length.

IV.2 Back-reaction on the Metric and Critical Length Scale

The above analysis provides a minimum length up to which there is a large uncertainty in defining the location of the horizon. The semi-classical Einstein equation describes the effect of the quantum fields on the background space-time [15]:

Gμν=8πTμνG_{\mu\nu}=8\pi\langle T_{\mu\nu}\rangle (34)

where the source Tμν\langle T_{\mu\nu}\rangle is the quantum expectation value of the matter field stress tensor operator. Note that the expectation value is only defined after suitable regularization and renormalization [30]. The procedure is well-defined away from the Planck energy [30]. The quantum expectation of the matter field will change the original background metric (2) to:

ds2=f~(r)dt2+dr2f~(r)+r2dΩ2,ds^{2}=-\widetilde{f}(r)dt^{2}+\frac{dr^{2}}{\widetilde{f}(r)}+r^{2}d\Omega^{2}\,, (35)

where

f~(r)=f(r)+h(r),\widetilde{f}(r)=f(r)+h(r)\,, (36)

and h(r)h(r) is the change in the background metric due to the quantum modes. We want to make a couple of remarks regarding the above-modified line-element: First, we have made a specific gauge choice of the perturbed metric. This choice of gauge is to keep the calculations tractable. It has also been shown that about 90%90\% of Hawking radiation is in ss-waves [31, 32]. Second, if the back-reaction on the metric is small, h(r)h(r) will be small compared to f(r)f(r). Naively, this implies that if h(r)f(r)h(r)\sim f(r), then the back-reaction on the background is significant. In the rest of this sub-section, we quantify this precisely by obtaining an effective 1-D equation for the quantum field Φ\Phi propagating in the above-modified line-element (35).

Substituting the above modified metric in the scalar field equation (16) and decomposing the modes into spherical harmonics, we obtain the following effective 1-D differential equation:

A(r)r2u~l(r,ω)+[ω2+C+(B2A)ArB24A12Br]u~l(r,ω)=0A(r)\partial_{r_{*}}^{2}\tilde{u}_{l}(r,\omega)+\left[\omega^{2}+C+\left(\frac{B}{2A}\right)\frac{\partial A}{\partial r_{*}}-\frac{B^{2}}{4A}-\frac{1}{2}\frac{\partial B}{\partial r_{*}}\right]\tilde{u}_{l}(r,\omega)=0 (37)

where

A\displaystyle A =\displaystyle= (1+h(r)f(r))2\displaystyle\left(1+\frac{h(r)}{f(r)}\right)^{2} (38a)
B\displaystyle B =\displaystyle= (1+h(r)f(r))(dh(r)drh(r)f(r)df(r)dr)\displaystyle\left(1+\frac{h(r)}{f(r)}\right)\left(\frac{dh(r)}{dr}-\frac{h(r)}{f(r)}\frac{df(r)}{dr}\right) (38b)
C\displaystyle C =\displaystyle= (f(r)+h(r))(1rddr[f(r)+h(r)]+l(l+1)r2)\displaystyle-\left(f(r)+h(r)\right)\left(\frac{1}{r}\frac{d}{dr}\left[f(r)+h(r)\right]+\frac{l(l+1)}{r^{2}}\right) (38c)

Note that the above expression is exact. Let us now compare the (37) with Eq. (18). As mentioned above, in the limit h(r)/f(r)1{h(r)}/{f(r)}\ll 1, the two equations are identical, and the effect of the fluctuation on the wave equation can be treated as negligible. This is also evident if we compare the coefficient of r2u~l(r,ω)\partial^{2}_{r*}\tilde{u}_{l}(r,\omega) in both the equations. This implies that if the following condition:

2h(r)f(r)1\frac{2h(r)}{f(r)}\sim 1 (39)

is satisfied, the back-reaction of the quantum field on the metric will be significant.

We want to make the following two comments about the result: First, we have assumed that the metric perturbation h(r)h(r) is time-independent. However, the back-reaction of the field modes on the background metric is generally time-dependent. In Appendix (C), we have discussed the more general case in which the perturbation of the metric is time-dependent, i.e. h(r,t)h(r,t). The time-average of A(r,t)A(r,t) is identical to A(r)A(r) in Eq. (38a). Hence, the analysis in this section carries through for time-averaged metric perturbations. Second, this expression is valid for any spherically symmetric space-time.

In the case of 4-dimensional Schwarzschild, Reissner-Nordström, SdS and SAdS black-holes, the perturbation of the quantum modes on the metric takes the following simple form:

h(r)=2ΔErh(r)=\frac{-2\Delta E}{r} (40)

where ΔE\Delta E is given by (27). Substituting the above expression in (39) and using (27), we have:

4(κ2πrf)(rhγ)14\left(\frac{\hbar\kappa}{2\pi rf}\right)\left(\frac{r_{h}}{\gamma}\right)\sim 1 (41)

The assumption of a free scalar field (causing no backreaction on the metric) on a classical background will no longer hold when this condition is satisfied. This condition will give a length-scale in which the original Hawking’s derivation can not be trusted. In the rest of this section, we obtain the critical distance, by evaluating the above expression (41), for three different cases :

  1. 1.

    Schwarzschild: In this case, Eq. (41) translates to:

    γcritS[8Mπ]1/3{\gamma^{S}_{crit}}\sim\left[\frac{8M\hbar}{\pi}\right]^{1/3} (42)

    where MM is the mass of the black-hole. The above expression matches with the results of Casher et al. [16]. As mentioned earlier, for Schwarzschild, Eq. (33) can be evaluated, and it matches with the above result (42).

  2. 2.

    SAdS: In this case, (41) translates to:

    γcrit(SAdS)[4π(rh3(rhL)2+1)]1/3=[4Lπ(χ3χ2+1)]1/3{\gamma^{(SAdS)}_{crit}}\sim\left[\frac{4}{\pi}\left(\frac{\hbar r_{h}}{3\left(\frac{r_{h}}{L}\right)^{2}+1}\right)\right]^{1/3}=\left[\frac{4L}{\pi}\left(\frac{\hbar\chi}{3\chi^{2}+1}\right)\right]^{1/3} (43)

    As expected, in the limit LL\rightarrow\infty, the critical length for the SAdS black-hole reduces to the Schwarzschild result. Since, for SAdS, we can not explicitly obtain a closed-form expression (33), we can not compare with the above expression. However, note that both expressions have the same power of γ\gamma and do point the similarity of both the expressions.

  3. 3.

    SdS: As mentioned earlier, SdS admits two (event and cosmological) horizons. We study the backreaction effect near the event horizon only. In Sec. (VI), we discuss the results for the cosmological horizon. In this case, the condition (41) translates to:

    γcrit(4Lπ)1/3(κeff|κb|)1/3(χb13χb2)1/3,\gamma_{crit}\sim\left(\frac{4\hbar L}{\pi}\right)^{1/3}\left(\frac{\kappa_{eff}}{|\kappa_{b}|}\right)^{1/3}\left(\frac{\chi_{b}}{1-3\chi_{b}^{2}}\right)^{1/3}\,, (44)

    where,

    χb(0,13),κeff|κb|=(2+χbχc)(4+χbχc+χcχb).\chi_{b}~{}\in~{}\left(0,\frac{1}{\sqrt{3}}\right)\,,\,\quad\quad\frac{\kappa_{eff}}{|\kappa_{b}|}=\frac{\left(2+\frac{\chi_{b}}{\chi_{c}}\right)}{\left(4+\frac{\chi_{b}}{\chi_{c}}+\frac{\chi_{c}}{\chi_{b}}\right)}\,. (45)

    Using (6), we have:

    χc=134χb2χb2\chi_{c}=\sqrt{1-\frac{3}{4}\chi_{b}^{2}}-\frac{\chi_{b}}{2} (46)

    χb13\chi_{b}\rightarrow\frac{1}{\sqrt{3}} corresponds to the the Nariai limit [24]. The Nariai limit corresponds to the degenerate horizon and the near horizon assumption being Rindler (14) is not valid. From Eq. (44) it is evident that γcrit\gamma_{crit} diverges in this limit. In Sec. (VI), we show that our analysis breaks down much before the Nariai limit.

Note that in obtaining the critical length for SAdS and SdS, we have not made any assumption about the value of χ\chi. To understand the difference between SAdS and Schwarzschild, we now compute the thermodynamical quantities in the two cases for the same MM. [While we provide expressions for SAdS, we plot for both SAdS and SdS.] To do that we first write the ratio of γcrit\gamma_{\rm crit}’s:

γR=ρcritSAdSρcritS=[1(1+χ2)(1+3χ2)]1/3\gamma_{R}=\frac{\rho^{\rm SAdS}_{\rm crit}}{\rho^{\rm S}_{\rm crit}}=\left[\frac{1}{(1+\chi^{2})(1+3\chi^{2})}\right]^{1/3} (47)

Using Eq. (26), the dimensionless entropy at the critical length scale is:

SRSAdS=SSAdSSS=[(3χ2+1)(1+χ2)2]2/3=NSAdSNSS^{\rm SAdS}_{R}=\frac{S^{\rm SAdS}}{S^{S}}=\left[\frac{(3\chi^{2}+1)}{(1+\chi^{2})^{2}}\right]^{2/3}=\frac{N^{\rm SAdS}}{N^{S}} (48)

Fig. 3 is the plot of SRSAdS(SRSdS)S^{\rm SAdS}_{R}(S^{\rm SdS}_{R}) as a function of χ\chi. It is interesting to note that the entropy for the SAdS black-hole is the same as the Schwarzschild black-hole in the obvious case of χ0\chi\rightarrow 0 but also in the case χ1\chi\rightarrow 1 when MM and LL are similar. It is important to note that the entropy has a maximum close to χ0.6\chi\sim 0.6.

Refer to caption
(a) SAdS
Refer to caption
(b) SdS
Figure 3: Plot of ratio of Entropies with χ\chi for SAdS and SdS.

Similarly, using Eqs. (25) and (27), we obtain dimensionless energy and fluctuations in the energy at the critical invariant distance:

ERSAdS\displaystyle\left\langle E\right\rangle^{SAdS}_{R} =\displaystyle= ESAdSES=κSAdSNSAdSκSNS=[(3χ2+1)5(1+χ2)]1/3\displaystyle\frac{\left\langle E\right\rangle^{SAdS}}{\left\langle E\right\rangle^{S}}=\frac{\kappa^{SAdS}N^{SAdS}}{\kappa^{S}N^{S}}=\left[\frac{(3\chi^{2}+1)^{5}}{(1+\chi^{2})}\right]^{1/3} (49)
(ΔE)RSAdS\displaystyle(\Delta E)^{SAdS}_{R} =\displaystyle= (ΔE)SAdS(ΔE)S=κSAdSκS(NSAdSNS)1/2=[(3χ2+1)4(1+χ2)]1/3\displaystyle\frac{(\Delta E)^{SAdS}}{(\Delta E)^{S}}=\frac{\kappa^{SAdS}}{\kappa^{S}}\left(\frac{N^{SAdS}}{N^{S}}\right)^{1/2}=\left[(3\chi^{2}+1)^{4}(1+\chi^{2})\right]^{1/3} (50)
Refer to caption
(a) SAdS
Refer to caption
(b) SdS
Figure 4: Plot of ratio of energies with χ\chi for SAdS and SdS
Refer to caption
(a) SAdS
Refer to caption
(b) SdS
Figure 5: Plot of ratio of fluctuations in Energy with χ\chi for SAdS and SdS

From Fig. (4) and Fig. (5) , we see that the energy and fluctuation in the energy of SAdS increase monotonically with χ\chi. In the case of SdS, the energy and fluctuation in energy each has one maximum. In Sec. (VI), we discuss the implications of these for our results.

While we were able to obtain a critical length scale using the fluctuations of thermal modes’ energy, one of the drawbacks is that we treated h(r)h(r) in Eq.(35) as a classical perturbation rather than a quantum fluctuation. We will overcome this in Approach 2, which we discuss in the next section.

V Approach 2: Quantum Field Theory

This section discusses the second approach to obtain a critical length scale at which the gravitational interactions become significant. We study the interaction of the thermal particles in the black-hole atmosphere, with a point-like massless particle falling radially into the black-hole (ss-wave). (The reason for this choice is because it has been shown numerically that about 90%90\% of Hawking radiation is in ss-waves [31, 32].) A massless point-like particle propagating in the Minkowski space-time produces a shock-wave with which the thermal particles of the atmosphere interact [25, 26].

As discussed in Sec. (II), we will use the shock-wave analysis in the Minkowski space-time to study the near horizon properties in any non-degenerate spherically symmetric black-hole space-time. This is because the near-horizon geometry can be approximated as a Rindler, which has a one-to-one mapping with the Minkowski space-time [33, 34]. We show that the shock-wave interaction with the thermal particles becomes significant as the infalling particle reaches a critical invariant distance from the horizon. We assume that the infalling particle (ss-wave), with momentum pp, to move along the vv-direction (u=u0=u=u_{0}= constant). The line-element for the space-time with the shock-wave is [25, 16]:

ds2=du[dv+2pln(x~2rh2)δ(uu0)du]+dx2+dy2ds^{2}=-du\left[dv+2p\ln\left(\frac{\tilde{x}^{2}}{r_{h}^{2}}\right)\delta\left(u-u_{0}\right)du\right]+dx^{2}+dy^{2} (51)

where x~2=x2+y2\widetilde{x}^{2}=x^{2}+y^{2} and rhr_{h} is the horizon radius. We would like to make a few remarks regarding the above line-element: First, the presence of rhr_{h} in the line-element is to provide a length-scale to compare with x~\widetilde{x}. Thus, in principle, rhr_{h} can be replaced by mass or energy of the black-hole as these also provide a length-scale of the system. Second, the above line element is different from the Minkowski line-element (72) just due to a shift in the vv-direction:

vv+2pln(x~2rh2)Θ(uu0)v\rightarrow v+2p\ln\left(\frac{\tilde{x}^{2}}{r_{h}^{2}}\right)\Theta(u-u_{0}) (52)

Hence, the shock-wave line-element remains the same as Minkowski everywhere except in the u=u0u=u_{0} hyperplane. See Appendix (A) for details on the shock-wave interaction with thermal particles.

Our interest is to obtain the critical length where the interaction between the shock-wave and thermal particles is significant and can not be neglected. One way to calculate the critical length is to calculate the probability of the thermal particles to be in the same state after interaction [16] with the shockwave. Thus, the probability will be close to unity if the interaction is weak; however, it will be negligible if it is strong. To go about calculating this probability, we describe the thermal particles by a wavepacket (g(k)g(k)) with a small spread [that depends on ll, cf. (23)]:

ϕ(x)=x|ϕ=N𝑑kg(k)eik(xT)\phi(x)=\langle x|\phi\rangle=N\int dk\,g(k)\,e^{ik(x-T)} (53)

g(k)g(k) can, in principle, be any function which fixes the momentum width of the wave packet. Using |ϕ|ϕ|2=1|\left\langle\phi|\phi\right\rangle|^{2}=1, we fix NN to be 12π\frac{1}{\sqrt{2\pi}}. The effect of the shockwave on the thermal particles can be described by the shift in vv coordinate given by:

2pln(x~2M2)p,2p\ln\left(\frac{\widetilde{x}^{2}}{M^{2}}\right)\sim p\,,

where pp is the momentum of the incoming particle (along the vv-direction) and is a constant of motion (For details, see Appendix (A).) Note that the log\log term will not contribute much in the final expression for the probability. Hence, the wave packet after crossing the shockwave would be:

ϕ(x)=x|ϕN𝑑kg(k)eik(xT+p),\phi^{\prime}(x)=\langle x|\phi^{\prime}\rangle\simeq N\int dk\,g(k)\,e^{ik\left(x-T+p\right)}\,, (54)

Using the fact that the momentum of a thermal particle near the horizon is 1/γ01/\gamma_{0} (Appendix (B)) before it interacts with the shockwave, we choose g(k)g(k) to be a Gaussian wave-packet:

|g(k)|2=12π(Δk)2e(kk¯)22(Δk)2|g(k)|^{2}=\frac{1}{2\pi(\Delta k)^{2}}e^{-\frac{-(k-\bar{k})^{2}}{2(\Delta k)^{2}}} (55)

where

Δk=k¯=1γ0\Delta k=\bar{k}=\frac{1}{\gamma_{0}} (56)

We have chosen Δk\Delta k to be of the same order as k¯\bar{k}. Thus, the probability of the thermal particle to be in the same state even after interacting with the shockwave is given by:

P=|ϕ|ϕ|2exp(4λ2κ2γ04)P=\left|\left\langle\phi^{\prime}|\phi\right\rangle\right|^{2}\sim\exp\left(-\frac{4\lambda^{2}}{\kappa^{2}\gamma_{0}^{4}}\right) (57)

where λ\lambda is related to the momentum of the infalling massless particle (See Eq. (90)):

λpκγ02andγ0rhl.\lambda\sim\frac{p\,\kappa\,\gamma_{0}}{2}\quad\mbox{and}\quad\gamma_{0}\sim\frac{r_{h}}{l}\,. (58)

Note that γ0\gamma_{0} is the maximum invariant distance that any thermal particle can reach. Probability of the entire thermal atmosphere to not change state is:

Ptotexp(4λ2κ2l(2l+1)γ0(l)4)exp(4λ2κ2γ4N(γ))P_{\rm tot}\sim\exp\left(-\frac{4\lambda^{2}}{\kappa^{2}}\sum_{l}\frac{(2l+1)}{\gamma_{0(l)}^{4}}\right)\sim\exp\left(-\frac{4\lambda^{2}}{\kappa^{2}\gamma^{4}}N(\gamma)\right) (59)

In obtaining the second expression, we approximate the summation by treating γ\gamma to be a free parameter. Eq. (59) is a vital result regarding which we would mention the following points: First, the result we have obtained is different from the one obtained by Casher et al. [16]. We choose a particular g(k)g(k) to get the exact result for the probabilities instead of making approximations. Our results match with that in Ref. [16] in the limit when the exponential is approximated to the first two terms in the series expansion. Second, We have obtained the probability for a general spherically symmetric metric.

From the above expression (59), we can obtain the critical length at which the interactions become strong. When the interaction is weak, most thermal particles will remain in the same state. This corresponds to:

4λ2κ2γ4N(γ)<<1.\frac{4\lambda^{2}}{\kappa^{2}\gamma^{4}}N(\gamma)<<1\,. (60)

This implies that if the following condition:

4λ2κ2γ4N(γ)1\frac{4\lambda^{2}}{\kappa^{2}\gamma^{4}}N(\gamma)\sim 1 (61)

is satisfied, the interaction between the shockwave and the thermal particles become significant. This tells us that up to the following invariant distance from the event horizon:

γ(2λrhκ)1/3(rh2pl)1/3\gamma\sim\left(\frac{2\lambda r_{h}}{\kappa}\right)^{1/3}\sim\left(\frac{r_{h}^{2}\,p}{l}\right)^{1/3} (62)

the interactions cannot be neglected. Note that (62) is related to pp (constant of motion). Physically, for a black-hole with horizon rhr_{h}, pp can not take arbitrarily small values. More specifically, the de Broglie wavelength of the incoming particle can not be larger than the event horizon. If the de Broglie wavelength is larger than the horizon radius, there will not be any interaction. Thus, a black-hole with event horizon rhr_{h} sets the minimum de Broglie wavelength of the incoming particle. This leads to the following condition:

hp<rhp>(2πrh)\frac{h}{p}<r_{h}\,\implies\,p>\left(\frac{2\pi\hbar}{r_{h}}\right) (63)

Substituting the minimum possible value of pp in Eq. (62), ignoring prefactor of π/l\pi/l, we have:

γcrit(2rh)1/3\gamma_{\rm crit}\sim(2\hbar r_{h})^{1/3} (64)

This is an important result regarding which we would like to stress the following points: First, the critical distance we have obtained above is for any spherically symmetric space-time (2) with the event horizon at rhr_{h}. In approach 1, we could obtain a closed-form expression only for a class of spherically symmetric black-hole space-times. However, in approach 2, we do not any make any such assumption. Second, recently, such a feature has emerged in describing the fluid-gravity correspondence leading to the quantization of horizon area [35, 36].

To compare and contrast the two approaches, in the rest of this section, we evaluate the critical distance (64) for the three black-hole space-times:

  1. 1.

    Schwarzschild: In this case, Eq. (64) reduces to:

    γcritS(4M)1/3{\gamma^{S}_{\rm crit}}\sim\left(4\hbar M\right)^{1/3} (65)

    where MM is the mass of the black-hole. This matches with the expression (42) (besides the overall constant). Also, the above expression matches with the corresponding results of Casher et al. [16].

  2. 2.

    SAdS: In this case, Eq. (64) translates to:

    γcritSAdS(2L)1/3χ1/3\gamma_{\rm crit}^{\rm SAdS}\sim(2\hbar L)^{1/3}\chi^{1/3} (66)
  3. 3.

    SdS: Near the event horizon, Eq.  (64) reduces to:

    γcritSdS(2rb)1/3(2L)1/3χb1/3\gamma_{\rm crit}^{\rm SdS}\sim(2\hbar r_{b})^{1/3}\sim(2\hbar L)^{1/3}\chi_{b}^{1/3} (67)

To get a numerical intuition of the critical length scales, let us calculate the numerical value of γcrit\gamma_{crit} using Eq. (64) for a supermassive Schwarzschild blackhole with a billion solar mass. In SI units, Eq. (64) can be rewritten as:

γcrit=(2rh)1/3=(2lpl2rh)1/3.\gamma_{crit}=(2\hbar r_{h})^{1/3}=(2l_{\rm pl}^{2}r_{h})^{1/3}\,.

where lpl=G/c3l_{pl}=\sqrt{\hbar G/c^{3}}. For M=109MM=10^{9}M_{\odot}, we have

rh=2GMc22.95× 1012m.r_{h}=\frac{2GM}{c^{2}}\sim 2.95\,\times\,10^{12}\,m\,.

Substituting lpl1035ml_{\rm pl}\sim 10^{-35}m and rhr_{h} in γcrit\gamma_{crit}, we get γcrit1.541× 1019m\gamma_{crit}\sim 1.541\,\times\,10^{-19}\,m. The value is extremely small (order 1019m10^{-19}\,m) because γcrit=(2rh)1/3=(2lpl2rh)1/3\gamma_{crit}=(2\hbar r_{h})^{1/3}=(2l_{pl}^{2}r_{h})^{1/3} in terms of the Schwarzschild radius rhr_{h}, is suppressed by the Planck length (lpll_{pl}). Thus, we see that for a supermassive black-hole of a billion solar masses, the critical length scale (width of the shell) at which the back-reactions become significant is of the order of 1019m10^{-19}\,m.

Before we proceed with the limitations of the analysis, we want to highlight the following points: First, in obtaining the critical length for SAdS and SdS, we have not assumed the value of χ\chi. Second, these results are consistent with the results obtained by Sorkin using a different procedure [37]. Third, both the approaches can not be trusted beyond γcrit\gamma_{crit}. Around γcrit\gamma_{crit}, the gravitational interactions become strong, and the shock-wave analysis [25] ceases to hold. Lastly, the critical length for both the approaches match for Schwarzschild, they differ for SAdS (43) and SdS (44). In the next section, we identify the cause of the deviation between the two approaches for SAdS and SdS and point to the limitations of the analysis.

VI Limitations of the analysis

Using two different approaches, for general spherically symmetric black-hole space-times, we obtained a critical scale at which the quantum scalar field’s backreaction to the metric becomes significant. While the critical length in the two approaches matches for Schwarzschild black-hole, the functional dependence of the critical length differs in the two approaches for SAdS and SdS black-holes. The plots in Fig. (6) highlight the difference between the critical length (γcrit\gamma_{\rm crit}) as a function of χ\chi in the two approaches for SAdS and SdS.

Refer to caption
(a) SAdS
Refer to caption
(b) SdS
Figure 6: Plot of the critical length γcrit/(4L)1/3\gamma_{\rm crit}/(4\hbar L)^{1/3} as a function of χ\chi (χb)\chi_{b}) in the two approaches for SAdS and SdS.

This is a crucial observation regarding which we want to highlight the following points: First, for SAdS and SdS, the two approaches seem to have different profiles. In the small χ\chi (χb)(\chi_{b}) limit, they agree. Second, in the case of SAdS, while approach 1 predicts finite critical length for all χ\chi, approach 2 predicts a growing critical length for large values of χ\chi. Note that this is the limit (χ1\chi\to 1) considered in AdS/CFT correspondence where the large black holes in AdS are dual to the high-temperature phase of the dual field theory [23]. Third, in the case of SdS, while approach 2 predicts finite critical length for all possible values of χ\chi, approach 1 predicts a diverging critical length for χb1/3\chi_{b}\to 1/\sqrt{3}.

Refer to caption
Figure 7: No Peak in SAdS Potential near horizon for ML=4\frac{M}{L}=4, (l=10l=10)

This leads to the important question: Why is there a discrepancy between the two approaches for SAdS and SdS? As mentioned in Sec. (III), our analysis rests on the fact that there is a peak (positive barrier) in the potential profile V(r)V(r) that localizes the high angular momentum modes close to the horizon. However, this assumption may not be true for all values of χ\chi. In Fig. (7), we have plotted the potential (4M2V)(4M^{2}V) for SAdS for specific values of M/L=4M/L=4 and l=10l=10. As we can see, the peak in the potential is not present, and the potential increases monotonically. Compare this with the plots in Fig. (2), where the potential has a peak at a finite value from the horizon. Specifically, we observe the absence of a peak near the horizon in the potential for M/L1.5{M}/{L}\gtrsim 1.5 and any l<100l<100. Thus, as we increase the value of M/LM/L (corresponding to increasing values of χ\chi), we observe that the peak in the potential slowly vanishes.

Refer to caption
Figure 8: SdS Potential (is neagative) near horizon for ML=0.193\frac{M}{L}=0.193, for different values of l

Fig. (8) is the plot of the potential (4M2V)(4M^{2}V) for the event horizon of SdS for M/L=0.193M/L=0.193. The plot highlights that the potential is negative for any value of ll. Hence, Eq. (21) will not provide any finite value of lmaxl_{\rm max}. Beyond a range of M/L values for SdS, the number of relevant modes is not finite.

This leads to the following: First, there is no upper bound on the angular momentum implying that all the angular momentum modes have to be considered. This means that γ0\gamma\to 0 and the semi-classical analysis is bound to breakdown. Second, for Approach 1, from Eqs. (25, 27), we have:

ΔEE1N0,\frac{\Delta E}{\langle E\rangle}\sim\frac{1}{\sqrt{N}}\sim 0\,, (68)

implying that the backreaction on the metric is negligible. Third, for approach 2, this limit implies that the condition (60) is never true. Physically, this means that gravitational interactions are strong everywhere. When the gravitational interactions are strong everywhere, the shock-wave metric (51) is not a valid metric in this limit. Hence, our analysis will not hold in this limit.

As mentioned earlier, both the approaches require a finite value of lmaxl_{\rm max} (cf. (23)) which is true only if the effective 1-D potential V(r)V(r) (20) has a maximum near the horizon. We have shown that the two approaches give similar critical length if the conditions (21, 23) are satisfied.

VII Conclusions and Discussions

Hawking’s original derivation ignores the radiation’s interaction with matter and requires the incoming vacuum modes with frequencies far above the Planck scale. One way to address both these assumptions is not to make the horizon a fixed surface, but a dynamical one that changes due to the backreaction of the quantum modes [16, 37]. Such an approach is also consistent within the framework of effective field theory [11, 13].

In this work, we aimed to ask until which length scale the quantum modes’ backreaction can be ignored. We identified the length scale where the assumption of a free field on a classical background to explain Hawking radiation is not valid, and the quantum modes smear the horizon. While there have been attempts in the literature, our analysis is different for two specific reasons: First, we have obtained the length scale for an arbitrary 4-dimensional spherically symmetric space-time. We also showed that for a supermassive Schwarzschild black-hole of a billion solar masses, the critical length scale (width of the shell) at which the back-reactions become significant is of the order of 1019m10^{-19}\,m. Second, we have identified the limitations of the approach in the asymptotic non-flat space-times.

We used two different approaches to evaluate the critical length. In the first approach, we used the statistical mechanical properties of the perturbations and calculated the change in the mass of the black-hole due to these fluctuations. In the second approach, we calculated the effect of a shock-wave (due to the thermal particles) on the outgoing Hawking particle. While the two approaches use two different features of the fluctuations, both the approaches assume a finite value of lmaxl_{\rm max}. We showed that as long as this condition is satisfied, both the approaches provide similar results. The results of approach 2 are consistent with Ref. [37].

We showed that this condition is explicitly violated for SAdS and SdS for a finite value of M/LM/L. In SdS, as M/LM/L is increased, the potential is negative for any value of ll. Hence, our results do not apply to the cosmological horizon and in the Nariai limit. In the case of SAdS, we showed that as M/LM/L is increased (corresponding to increasing values of χ\chi), the peak in the potential vanishes. The difference in the Free energy of AdSd+1\mathrm{AdS}_{d+1} with the BH\mathrm{BH} solution (FbhF_{\rm bh}) and free energy of AdSd+1\mathrm{AdS}_{d+1} without black-hole (FthF_{\rm th}) is given by [38]

FbhFth=rhd1VSd1drh2+(d2)L2(L2rh2)F_{\rm bh}-F_{\rm th}=\frac{r_{h}^{d-1}V_{S^{d-1}}}{dr_{h}^{2}+(d-2)L^{2}}\left(L^{2}-r_{h}^{2}\right) (69)

where VSd1V_{S^{d-1}} is the volume of the (d1)(d-1)-sphere. The above expression implies that rh=Lr_{h}=L is the critical size. For rh<Lr_{h}<L, the thermal AdS is less energetic solution while the opposite is true in the limit rh>Lr_{h}>L. Our results imply that critical distance where the back-reaction becomes strong can be computed in the limit of M/L<1M/L<1 (rh/L<1r_{h}/L<1).

Counter-intuitively, our results are more accurate in the late stages of the evaporation process. As black-holes lose mass due to evaporation, M/LM/L decrease, and both the approaches give similar results. However, for the Planck-size black-holes, General Relativity itself needs to be modified to include higher-order Riemann tensor, Ricci tensor and scalar terms (1). In that scenario, the background metric line-element (2) itself may be modified. It is important to note that our calculations do not use any of the physical inputs of the Firewall proposal [22], however, our results may have implications for the Firewall proposal. This is currently under investigation.

One of the reasons for choosing a massless scalar field is that the gravitational perturbations about the background metric have the same form [27, 28]. For the modified gravity theories, we need to obtain the corresponding equation to study the back-reaction. This is currently under investigation.

Our results bring attention to the following interesting questions: Are there modifications to the Hawking temperature, especially in the limit M/L1M/L\sim 1 (for SAdS)? How does the Hawking rate change for these space-times? How is the Page time [39] affected in the limit M/L1M/L\sim 1 (for SAdS)? We hope to return to study some of these questions soon.

Acknowledgements

The authors thank M. Bojowald, T. Dray, Ted Jacobson, Sudipta Mukherji, T. Padmanabhan, P. Ramadevi, A. Virmani, and U. Yajnik for discussions. We are grateful to Mahendra Mali for his help during the initial stages of this work. We thank Saurya Das and S. Mahesh Chandran for helpful comments on an earlier draft. This work is part of the Undergraduate project of MS. The work of SS is supported by the MATRICS SERB grant (MTR/2019/000077) and Homi Bhabha Fellowship.

Appendix A Shockwave Analysis

Here we discuss the effect of a point-like massless particle propagating in the Minkowski background. The Minkowski line-element is given by:

ds2=dT2+dz2+dx2+dy2,ds^{2}=-dT^{2}+dz^{2}+dx^{2}+dy^{2}\,, (70)

where

u=T+zv=Tz.u=T+z\quad\quad\quad\quad v=T-z\,. (71)

This leads to the following form of the Minkowski metric:

ds2=dudv+dx2+dy2ds^{2}=-dudv+dx^{2}+dy^{2} (72)

The above line-element can describe the region near the event horizon of 4-D spherically symmetric black-hole space-time described by the line-element (2). To see this, we define the following:

u=γeκtv=γeκtu=\gamma e^{\kappa t}\quad\quad\quad v=-\gamma e^{-\kappa t} (73)

where κ\kappa is the surface gravity defined in (9) and γ\gamma is given by (12). Under this transformation, the above line-element (72) leads to:

ds2=κ2γ2dt2+dγ2+dx2+dy2ds^{2}=-\kappa^{2}\gamma^{2}dt^{2}+d\gamma^{2}+dx^{2}+dy^{2} (74)

Note that the above line-element corresponds to the near-horizon line-element of a spherically symmetric black-hole obtained in (14) in the xtx-t plane. The 22-Sphere part of the metric can also be mimicked (See [33, 34]). Dray and t’Hooft [25] showed that a massless particle moving in Minkowski space produces a shockwave. Due to the mapping described above, the analysis of Dray and t’Hooft [25] can be extended to any horizon. In the rest of this Appendix, we give the key results which are essential for our calculations. The gravitational field of a massless particle in Minkowski space moving along the vv direction (u=u0u=u_{0}), with momentum pp, is described by the line element [16]:

ds2=du[dv+2pln(x~2rh2)δ(uu0)du]+dx2+dy2ds^{2}=-du\left[dv+2p\ln\left(\frac{\tilde{x}^{2}}{r_{h}^{2}}\right)\delta\left(u-u_{0}\right)du\right]+dx^{2}+dy^{2} (75)

with x~2=x2+y2\widetilde{x}^{2}=x^{2}+y^{2} and rhr_{h} is the horizon radius. rhr_{h} just provides a scale (in the system concerned) to compare x~\widetilde{x} with. rhr_{h} in (75) could, in principle, be replaced by the mass or the energy of the black-hole as these also provide a scale in the system.. Note that the above line element is different from the Minkowski case only by a shift in the vv direction given by:

vv+cΘ(uu0)v\rightarrow v+c\,\Theta(u-u_{0}) (76)

with

c=2pln(x~2rh2)c=2p\ln\left(\frac{\tilde{x}^{2}}{r_{h}^{2}}\right) (77)

Therefore the space remains same as Minkowski everywhere except in the hyperplane u=u0u=u_{0}. The effect of such a shockwave on other particles can be studied using the Hamilton Jacobi formalism. The Hamilton Jacobi equation is given by (SS is action):

4(Su)(Sv)+8pln(x~2rh2)δ(uu0)(Sv)2+(Sx)2+(Sy)2=0-4\left(\frac{\partial S}{\partial u}\right)\left(\frac{\partial S}{\partial v}\right)+8p\ln\left(\frac{\tilde{x}^{2}}{r_{h}^{2}}\right)\delta\left(u-u_{0}\right)\left(\frac{\partial S}{\partial v}\right)^{2}+\left(\frac{\partial S}{\partial x}\right)^{2}+\left(\frac{\partial S}{\partial y}\right)^{2}=0 (78)

At uu0u\neq u_{0} (absence of shockwave), the solution is given:

S0=[kxx+kyy+kvv+(kx2+ky24kv)u]S_{0}=\hbar\left[k_{x}x+k_{y}y+k_{v}v+\left(\frac{k_{x}^{2}+k_{y}^{2}}{4k_{v}}\right)u\right] (79)

Here the kk’s are the wave vectors of the thermal particles before encountering the shockwave. In the presence of shockwave, it is natural to use (76) in (79) to get an action upto the leading correction, which gives:

S=S0+2pvln(x~2rh2)kvθ(uu0)+Subleading termsS=S_{0}+2p^{v}\ln\left(\frac{\tilde{x}^{2}}{r_{h}^{2}}\right)\hbar k_{v}\theta\left(u-u_{0}\right)+\mbox{Subleading terms} (80)

This action can then be used to calculate shifts (refractions) in the xx, yy directions as:

kx(u)=Sx=kx+4pvx~2xkvθ(uu0)+Subleading termsk_{x}(u)=\frac{\partial S}{\partial x}=k_{x}+\frac{4p^{v}}{\tilde{x}^{2}}xk_{v}\theta\left(u-u_{0}\right)+\mbox{Subleading terms} (81)
ky(u)=Sx=ky+4pvx~2ykvθ(uu0)+Subleading termsk_{y}(u)=\frac{\partial S}{\partial x}=k_{y}+\frac{4p^{v}}{\tilde{x}^{2}}yk_{v}\theta\left(u-u_{0}\right)+\mbox{Subleading terms} (82)

Appendix B Calculating the momentum of the Hawking Particles before interaction with the Shockwave

Here, we approximate the momentum of the thermal particles before they interact with the shockwave. In line-element (14), the massless particles follow the geodesics (with gμνV0μV0ν=0g_{\mu\nu}V_{0}^{\mu}V_{0}^{\nu}=0):

Xμ=X0μ+αV0μ.X^{\mu}=X^{\mu}_{0}+\alpha V_{0}^{\mu}\,. (83)

Using a suitable rotation, boost and translation, the above geodesic of the massless particle takes the following form [16]:

Xμ(α)=(T=γ0+α,y=0,x=γ0+α,z=ρ0)0<α<2γ0,X^{\mu}(\alpha)=\left(T=-\gamma_{0}+\alpha,\quad y=0,\quad x=-\gamma_{0}+\alpha,\quad z=\rho_{0}\right)\quad\quad 0<\alpha<2\gamma_{0}\,, (84)

where γ0rhl\gamma_{0}\sim\frac{r_{h}}{l} is the maximum invariant distance that any thermal particle can reach. While the above expression is general and can be applied to any physical situation, our interest is near-horizon. We have used Eq. (23) to relate γ0\gamma_{0} with the orbital angular momentum number \ell and the critical distance.

As in Ref. [16], we assume that the ingoing particle (that produces the shockwave, which is also massless) has momentum pp only along the vv-direction. In other words, it travels at a fixed uu-coordinate.

Now we calculate the γ\gamma at which a thermal particle crosses the shockwave. At that γ\gamma, we will try to estimate the momentum of the thermal particle. Using the thermal particle’s trajectory (84) for a fixed uu equal to that of the ingoing particle, we get (using (71) and (73)):

γcross2=2γ0uu2\gamma^{2}_{cross}=2\gamma_{0}u-u^{2} (85)

We make the assumption that most thermal particles that can reach the maximum γ=γ0\gamma=\gamma_{0}, cross the shockwave at γcross=γ0\gamma_{cross}=\gamma_{0} itself. I.e. most of the thermal particles that can reach γ0\gamma_{0} have already reached γ=γ0\gamma=\gamma_{0} before the shockwave arrives. Using this approximation and the above equation, we get:

γcrossγ0u\gamma_{cross}\sim\gamma_{0}\sim u (86)

Now we calculate the Schwarzschild energies (λ\lambda, negative of the time component of the four momentum covariant vector) of any particle is our spacetime:

λ=pt=κ2γ2pt\lambda=-p_{t}=\kappa^{2}\gamma^{2}p^{t} (87)

which after a short calculation in terms of uu and vv gives:

λ=pt=κ(vpvupu)\lambda=-p_{t}=\kappa(vp_{v}-up_{u}) (88)

For the infalling particle (s-wave):

pv=p;pu=0;pv=0;pu=p2p^{v}=p;\quad p^{u}=0;\quad p_{v}=0;\quad p_{u}=-\frac{p}{2} (89)

Thus we get a relation between the momentum of the particle and its Schwarzschild energy as:

p=2λκu2λκγ0p=\frac{2\lambda}{\kappa u}\sim\frac{2\lambda}{\kappa\gamma_{0}} (90)

For thermal particles we approximate the Schwarzschild energy to be equal to the kBTκk_{B}T\sim\hbar\kappa. Also for the thermal particles, near the horizon the uu and vv direction very nearly coincide (horizon is a null surface), we have ku=kv\hbar k_{u}=\hbar k_{v}. Using this in (88), we get:

ku1γ0k_{u}\sim\frac{1}{\gamma_{0}} (91)

Equation (91) gives an estimate of the momentum of a thermal particle before it crosses the horizon. To fix the width and the peak of the wave packet of the thermal particle before interaction with the shockwave we use (91) as an estimate.

Appendix C Effect of the time-dependent metric perturbation on the critical length scale

In Sec. (IV), we obtained the critical length at which the back-reaction of the quantum field on the metric will be significant assuming that the metric perturbation is time-independent. In this Appendix, we relax the condition and assume that the metric perturbation is time-dependent:

f~(r,t)=f(r)+h(r,t)\tilde{f}(r,t)=f(r)+h(r,t) (92)

We can use the following ansatz for the generic time-dependent field configuration (not necessarily harmonic):

ϕlm=ul(r,t)rYlm(θ,ϕ)\phi_{lm}=\frac{u_{l}(r,t)}{r}Y_{lm}(\theta,\phi) (93)

Substituting the above ansatz in the wave equation (16), we have:

(tf~tul+t2ul)+r~2ulf~(l(l+1)r2+1rrf~)ul=0-(\partial_{t}\tilde{f}\partial_{t}u_{l}+\partial_{t}^{2}u_{l})+\partial_{\tilde{r}^{*}}^{2}u_{l}-\tilde{f}\left(\frac{l(l+1)}{r^{2}}+\frac{1}{r}\partial_{r}\tilde{f}\right)u_{l}=0 (94)

where r~\tilde{r}^{*} is:

r~=drf~\tilde{r}^{*}=\int\frac{dr}{\tilde{f}} (95)

Substituting the above form of the metric (92) in the above wave-equation, we have

Unperturbed Eq. t2ul+r2ulf(l(l+1)r2+1rrf)ul=0\displaystyle~{}~{}~{}-\partial_{t}^{2}u_{l}+\partial_{r^{*}}^{2}u_{l}-f\left(\frac{l(l+1)}{r^{2}}+\frac{1}{r}\partial_{r}f\right)u_{l}=0 (96)
Perturbed Eq. t2u~l(r,t)+A(r,t)r2u~l(r,t)\displaystyle~{}~{}~{}-\partial_{t}^{2}\tilde{u}_{l}(r,t)+A(r,t)\partial_{r^{*}}^{2}\tilde{u}_{l}(r,t)
+[C+(B2A)ArB24A12Br+122ht2+14(ht)2]u~l(r,t)=0\displaystyle+\left[C+\left(\frac{B}{2A}\right)\frac{\partial A}{\partial r^{*}}-\frac{B^{2}}{4A}-\frac{1}{2}\frac{\partial B}{\partial r^{*}}+\frac{1}{2}\frac{\partial^{2}h}{\partial t^{2}}+\frac{1}{4}\left(\frac{\partial h}{\partial t}\right)^{2}\right]\tilde{u}_{l}(r,t)=0

where AA, BB, and CC, defined in Eq (38), are time-dependent. Note that the above equation is similar to (37) except for

122ht2ul+14(ht)2ul\frac{1}{2}\frac{\partial^{2}h}{\partial t^{2}}u_{l}+\frac{1}{4}\left(\frac{\partial h}{\partial t}\right)^{2}u_{l}

terms. Our interest is to study the effect of the fluctuations over a long period. Mathematically, this corresponds to replacing the above expression with time averages (which would be only functions of rr). Thus the time-dependence leads to a change in the potential of the perturbed wave equation; however, it does not change the coefficients of r2u~l\partial_{r^{*}}^{2}\tilde{u}_{l} terms, and the time average of A(r,t)A(r,t) corresponds to A(r)A(r) in (38). Hence, the analysis in Sec. (IV) follows.

References