Next Article in Journal
Correction to Yao, H.; Qiao, J.-W.; Gao, M.C.; Hawk, J.A.; Ma, S.-G.; Zhou, H. MoNbTaV Medium-Entropy Alloy. Entropy 2016, 18, 189
Next Article in Special Issue
Effects Induced by the Initial Condition in the Quantum Kibble–Zurek Scaling for Changing the Symmetry-Breaking Field
Previous Article in Journal
Microstructures of Al7.5Cr22.5Fe35Mn20Ni15 High-Entropy Alloy and Its Polarization Behaviors in Sulfuric Acid, Nitric Acid and Hydrochloric Acid Solutions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hawking-Like Radiation from the Trapping Horizon of Both Homogeneous and Inhomogeneous Spherically Symmetric Spacetime Model of the Universe

by
Subenoy Chakraborty
1,
Subhajit Saha
2 and
Christian Corda
3,4,5,6,*
1
Department of Mathematics, Jadavpur University, Kolkata 700032, India
2
Department of Physical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur 741246, India
3
Research Institute for Astronomy and Astrophysics of Maragha (RIAAM), Maragha 55134-441, Iran
4
Dipartimento di Fisica, Scuola Superiore di Studi Universitari e Ricerca “Santa Rita”, via Trasaghis 18/E, Roma 00188, Italy
5
Austro-Ukrainian Institute for Science and Technology, Institut fur Theoretishe Physik, Technische Universität, Wiedner Hauptstrasse 8-10/136, Wien A-1040, Austria
6
International Institute for Applicable Mathematics & Information Sciences (IIAMIS), B.M. Birla Science Centre, Adarsh Nagar, Hyderabad 500-463, India
*
Author to whom correspondence should be addressed.
Entropy 2016, 18(8), 287; https://doi.org/10.3390/e18080287
Submission received: 7 June 2016 / Revised: 21 July 2016 / Accepted: 28 July 2016 / Published: 8 August 2016
(This article belongs to the Special Issue Entropy in Quantum Systems and Quantum Field Theory (QFT))

Abstract

:
The present work deals with the semi-classical tunnelling approach and the Hamilton–Jacobi method to study Hawking radiation from the dynamical horizon of both the homogeneous Friedmann–Robertson–Walker (FRW) model and the inhomogeneous Lemaitre–Tolman–Bondi (LTB) model of the Universe. In the tunnelling prescription, radial null geodesics are used to visualize particles from behind the trapping horizon and the Hawking-like temperature has been calculated. On the other hand, in the Hamilton–Jacobi formulation, quantum corrections have been incorporated by solving the Klein–Gordon wave equation. In both the approaches, the temperature agrees at the semiclassical level.
PACS:
04.70.-D; 04.50.-h; 04.62.+v

1. Introduction

The quantum description of a black hole (BH), namely the Hawking radiation (HR) is closely related to the existence of an event horizon to the BH. The derivation of Hawking that “BH evaporates particles” [1,2,3] was based on quantum field theory. Hartle and Hawking [4] subsequently derived the BH temperature at the semiclassical level using the Feynmann path integral. The mathematical complexity of the above procedures forces to develop semi-classical approaches [5,6,7,8,9,10,11] for studying BH radiation. However, these semi-classical techniques were classified into two approaches—the tunnelling approach of Parikh and Wilczek [5,6,9,10,11,12,13,14,15,16] and the standard Hamilton–Jacobi (HJ) method (known as complex path integral formalism) by Padmanabhan et al. [7,8].
The energy conservation in tunnelling of a thin shell from the hole is the main ingredient for the first approach and is often referred to as the radial null geodesic method. The imaginary part of the action from the s-wave emission is connected to the Boltzmann factor for emission to relate with Hawking temperature (HT). In the complex paths method, the action S(r,t) for a single scalar particle is obtained by solving the Klein–Gordon (KG) equation in gravitational background. Then, HT is obtained using “the principle of detailed balance” [7,8]. A remarkable behavior of the energy conservation in tunnelling is that the final spectrum is not strictly thermal [9,10,11], and this has important consequences on the black hole information paradox [17].
On the other hand, the discovery of HR [1,2] completes the cycle to describe BH as a thermodynamical object. The black body nature of BH reveals emission of thermal radiation with temperature proportional to its surface gravity. Before the discovery of HR, entropy was formulated by Bekenstein [18] as proportional to the horizon area of the BH. This new perspective leads to viewing general relativity (GR) from a completely different angle. The nice interrelationship between gravity and thermodynamics was first enlightened by Jacobson [19] replacing the usual definition of heat by flow of energy across a horizon. In fact, gravity might be due to the thermodynamics of the microstructure of spacetime. Then, this idea is extended to modified gravity theories [20,21,22,23]. At present, there is a general consensus that gravity might be originated by the thermodynamics of the unknown microstructure of spacetime. Thus, it is speculated that BH thermodynamics is playing the role of a bridge to put GR and quantum mechanics on the same platform—the challenging issue of quantum gravity. It is also an intuitive but general conviction that BHs result in highly excited states representing both the “hydrogen atom” and the “quasi-thermal emission” in quantum gravity. A recent approach [24,25,26] has shown that such an intuitive picture is more than a picture, discussing a model of quantum BH somewhat similar to the historical semi-classical model of the structure of a hydrogen atom introduced by Bohr in 1913 [27].
Normally, the global concept of event horizon is used to define the HT. However, difficulty arises in dynamical spacetime, where there is no existence of an event horizon even locally. In the recent past, Hayward et al. [28] formulated a locally defined HT for dynamical BH using the tunnelling idea of Parikh–Wilczek [5]. Then, Cai et al. [29] have shown HR from the locally defined apparent horizon of the FRW Universe, where HT is measured by an observer using the Kodama vector [30,31,32] inside the horizon. In the present work, we shall address this important issue in GR, namely the formulation of thermodynamics of dynamical spacetime. We shall derive HR and the corresponding HT using both the HJ method and the tunnelling approach. A comparison between these two formalism will be done for homogeneous (FRW) as well as an inhomogeneous (LTB) spherically symmetric spacetime model. The paper is organized as follows: Section 2 gives an overview of radial null geodesic method to derive HT. In Section 3, the HJ method with quantum prescription for the FRW model has been presented. Quantum corrected entropy has been formulated in Section 4. Section 5 and Section 6 deal with inhomogeneous LTB model, where HT has been evaluated using both the approaches. In the last Section we write a brief summary of the whole work. Finally, for the sake of completeness the paper ends with an Appendix concerning the calculation of Kodama Vector and Surface Gravity.

2. Radial Null Geodesic Approach: Hawking-Like Temperature

In semiclassical tunnelling analysis [5,6,9,10,11,12,13,33], the radial null geodesic method is a common way of evaluating Hawking radiation. This method is simple as compared to the HJ method, but it has some limitations, namely:
  • The method is applicable only for massless particles;
  • One has to use only Painleve type coordinates to avoid singularity at the horizon;
  • There is a discrepancy of factor two in this method;
  • There is no general method to include quantum effects.
The basic idea in this semiclassical approximation (Wentzel–Kramers–Brillouin (WKB) approximation) is that the emission rate for the s-wave emission of a massless particle can be related to the imaginary part of the action of a system. In the s-wave approximation, particles can be viewed as a massless shell moving along radial null geodesic. Now, compared to static BH cases, there is one basic difference, as in the present case, the metric coefficients depend on both radius and time. Hence, there is no longer time translation Killing vector field. We shall have to use the Kodama [30], which is time like inside the horizon and the associate energy.
The homogeneous and isotropic model of the Universe is described by the FRW metric as:
d s 2 = d t 2 + a 2 ( t ) 1 κ r 2 d r 2 + a 2 ( t ) r 2 d Ω 2 2 ,
where the FRW coordinates (t, r, θ, ϕ ) are orthogonal comoving coordinates with “t” the comoving time corresponding to a comoving observer, and κ is the intrinsic spatial curvature.
For the above FRW model of the Universe, if we make a change of the radial coordinate r R , where R = a r is known as the area radius, then the above standard FRW metric becomes a Painleve–Gullstrand like metric [34,35] as follows:
d s 2 = 1 R 2 R A 2 1 κ R 2 a 2 d t 2 2 H R 1 κ R 2 a 2 d t d R + d R 2 1 κ R 2 a 2 + R 2 d Ω 2 2 ,
where R A = 1 H 2 + κ a 2 is the radius of the apparent horizon. In the present case, for the metric (2), the Kodama vector and the corresponding energy are:
κ μ = 1 κ R 2 a 2 , 0 , 0 , 0 a n d ω = 1 κ R 2 a 2 S t ,
respectively. Thus, ω 1 κ R 2 a 2 is the energy of the particle as measured by an observer with the Kodama vector. The differential equation for the radial null geodesic (i.e., d s 2 = 0 = d Ω 2 2 ) has the form
d R d t = H R ± H 2 R 2 + 1 R 2 R A 2 ,
where ± sign are associated with the outgoing/ingoing null geodesic with the assumption that “t” increases towards future. As we are interested in the imaginary part of the action corresponding to the tunnelling process through a barrier (the classically forbidden region) as according to Parikh-Wilczek [5]:
I m S = I m R i n R o u t p R d R = I m R i n R o u t 0 p R d p R d R = I m R i n R o u t 0 E d H R ˙ d R ,
where we have used the Hamiltonian equation
R ˙ = H p R = d H d p R | R
.
Here, p R is the radial momentum, R i n and R o u t are positions very close to the horizon with R i n as the initial position and R o u t , a classical turning point. Using the value of R from Equation (4) into Equation (5), we get
I m S = I m R i n R o u t d R d H R ˙ = I m R i n R o u t d R R ˙ ω 1 + κ R 2 a 2 = ω R i n R o u t d R 1 + κ R 2 a 2 H 2 R 2 + 1 R 2 R A 2 H R = π R A ω ,
where integration over the Hamiltonian H gives the energy of the particle as ω 1 κ 2 a 2 as measured by an observer with Kodama vector. Now, comparing the tunnelling probability
Γ 2 I m S
with the Boltzmann factor e x p ω T , we have the temperature associated with the apparent horizon of the FRW Universe as [35]
T A = 2 π R A .
This is the semiclassical Hawking-like temperature of the FRW Universe for tunnelling of massless particles across the apparent horizon.
It should be noted that Im r i n r o u t p r d r is not canonically invariant and hence it is not a proper observable. The canonically invariant quantity Im p r d r over a closed path across the horizon cannot coincide with 2 Im r i n r o u t p r d r for the Painleve coordinates as a particle experiences barrier only from inside the horizon to outside, not the other way. However, for the invariant definition Γ e i Im p r d r , the Hawking temperature turns out to be twice the original temperature. This ambiguity in the factor of two has been discussed in [36,37,38,39,40,41,42].

3. Hamilton–Jacobi Method: Quantum Prescription

In this section, we shall deal with the tunnelling of massless particles beyond the semiclassical approximation by Hamilton–Jacobi (HJ) method. We shall start with a KG equation for a scalar field ϕ describing a massless scalar particle of the form
2 g ( g μ ν g ν ) ψ = 0 .
The explicit form of the KG equation in the background of the FRW metric (1) is given by
2 ψ t 2 ( 1 κ r 2 ) a 2 2 ψ r 2 + H ψ t + κ r a 2 ψ r = 0 .
It should be mentioned that here, due to the spherical symmetry of the FRW spacetime and consideration of the radial trajectories only, we have considered ( t r )—sector in the spacetime given by Equation (1), i.e., 2D hyperplane (t,r). Now, substituting the standard ansatz for the semiclassical wave function
ψ ( r , t ) = e x p i S ( r , t )
into the wave Equation (10), the differential equation for the action S becomes
S t 2 1 κ r 2 a 2 S r 2 i 2 S t 2 1 κ r 2 a 2 2 S r 2 + H S t + κ r a 2 S r = 0 .
As a first step to solve this partial differential equation (PDE), we expand the action in powers of Planck constant ℏ as [36]
S ( r , t ) = S 0 ( r , t ) + Σ k k S k ( r , t )
with k, a positive integer. Here, terms of the order of Planck’s constant, and its higher powers, are considered as quantum corrections over the semiclassical action S 0 . Now, substituting this ansatz for S in the differential Equation (12) and equating different powers of ℏ on both sides, we obtain the following set of PDEs:
0 : S 0 t 2 1 κ r 2 a 2 S 0 r 2 = 0 ,
1 : 2 S 0 t S 1 t 2 1 κ r 2 a 2 S 0 r S 1 t i 2 S 0 t 2 1 κ r 2 a 2 2 S 0 r 2 + H S 0 t + κ r a 2 S 0 r = 0 ,
2 : S 1 t 2 + 2 S 0 t S 2 t 1 κ r 2 a 2 S 1 r 2 + 2 S 0 r S 2 r i [ 2 S 1 t 2 1 κ r 2 a 2 2 S 1 r 2 + H S 1 t + κ r a 2 S 1 r ] = 0 ,
and so on. The above expansion of the action (see Equation (13)) to all orders in ℏ is due to [36]. However, Yale [43] criticized this expansion of the action in ℏ as [36] does not take into account the higher order corrections to the energy—only zero order energy is considered. Subsequently, Singleton et al. [44] has made an attempt to justify the above expansion of the action in [36] using some ancillary assumption, namely, self-similarity. Anyway, one can say that the expansion in Equation (13) may be conditionally valid.
Apparently, different order PDEs are very complicated, but, fortunately, there will be lot of simplifications if, in the PDE corresponding to k , all previous PDEs are used and finally we obtain identical PDE, namely
k : S k t 2 1 κ r 2 a 2 S k r 2 f o r k = 0 , 1 , 2 . . . .
We see that different order quantum corrections satisfy identical differential equations as the semiclassical action S 0 , thus the correction terms are not independent, rather proportional to S 0 (i.e., S k S 0 , for all k). To determine these proportionality constants, we shall use dimension analysis. Since S 0 has the dimension of ℏ, the proportionality constant d k for S k has the dimension of k . However, in standard units, namely G = c = k B = 1 , the Planck’s constant ℏ is of the order M p 2 ( M p is the Planck mass) and hence d k has the dimension of M 2 k , where M is identified as the mass of the Universe. Thus, the series expansion (13) can be written in terms of S 0 as
S ( r , t ) = S 0 ( r , t ) 1 + Σ k α k M 2 k ,
with α k as dimensionless constant parameters.
Now, a complete solution for S requires the solution of S 0 satisfying the PDE (14). Since the metric (1) is non-static, there is no time like Killing vector in the dynamical FRW spacetime. However, Kodama vector [28] has a similar role in FRW spacetime as the time like Killing vector does in the stationary BH spacetime. Inside the apparent horizon, Kodama vector is a time like vector, and hence there is a conserved energy of a particle moving in the time like Killing vector in the stationary BH spacetime.
The Kodama vector for the FRW metric is given by
K μ = ( 1 κ r 2 , H r 1 κ r 2 , 0 , 0 ) ,
with associated conserved energy of the particle as
ω = 1 κ r 2 S 0 t + H r 1 κ r 2 S 0 r .
Solving for S 0 t and S 0 r , using Equations (14) and (20), we get
S 0 = ω d t { 1 κ r 2 H a r } a ω d r 1 κ r 2 { 1 κ r 2 H a r } ,
where the ∓ sign corresponds to ingoing and outgoing scalar particle, respectively, for which the wave functions have the expressions [36]
ψ i n = e x p i 1 + Σ k α k M 2 k ω d t { 1 κ r 2 H a r } + a ω d r 1 κ r 2 { 1 κ r 2 H a r }
and
ψ o u t = e x p i 1 + Σ k α k M 2 k ω d t { 1 κ r 2 H a r } a ω d r 1 κ r 2 { 1 κ r 2 H a r } .
Now, across the horizon, the metric coefficients in the (r,t) sector alter their sign. Thus, the above time integration might have imaginary parts and make contributions to the probabilities for the ingoing and the outgoing particles. As a result, the probabilities are given by
P i n = | ψ i n | 2 = e x p [ 2 1 + Σ k α k M 2 k I m { ω d t { 1 κ r 2 H a r } + I m a ω d r 1 κ r 2 { 1 κ r 2 H a r } } ] ,
and
P o u t = | ψ o u t | 2 = e x p [ 2 1 + Σ k α k M 2 k I m { ω d t { 1 κ r 2 H a r } I m a ω d r 1 κ r 2 { 1 κ r 2 H a r } } ] .
In the present context for the cosmological spacetime, the outgoing probability has to be unity in the classical limit 0 , as there will be no observer and everything will go out [45]. Note that the situation is the opposite of what one has for the BH case, where P i n = 1 in the limit 0 . Hence, from Equation (25), in the limit 0 , we get,
I m ω d t { 1 κ r 2 H a r } = I m a ω d r 1 κ r 2 { 1 κ r 2 H a r } .
The above two integrals over t and r, respectively, are explicitly shown to be equal for the FRW model in [35]. Note that, in most tunnelling problems, one not only has a spatial integral contribution to the tunnelling rate but also a time integral contribution, which is critical to getting the correct thermodynamic properties (namely, temperature and entropy) for these spacetimes in the tunnelling picture (either null geodesic or HJ). This time contribution was first pointed out in [46,47,48]. Furthermore, this time contribution to the tunnelling amplitude resolves the “factor of 2” puzzle mentioned in Section 2. Moreover, in [46,47,48], it has been shown that Painleve types of coordinates are not the only coordinates to avoid singularity at the horizon, one may use Schwarzschild or Kruskal–Szekeres coordinates in dealing with tunnelling formulation, i.e., singularity at the horizon is not an issue if treated properly. Thus, P i n will have the simplified form
P i n = e x p 4 a ω 1 + Σ k α k M 2 k I m d r 1 κ r 2 { 1 κ r 2 H a r } = e x p 2 ω 1 + Σ k α k M 2 k π R a .
Hence, from “the principle of detailed balance”, i.e., [7,8]
P o u t = e x p ω T h P i n ,
we have
T k = 1 + Σ k α k M 2 k 1 1 2 π R A = 1 + Σ k α k M 2 k 1 T c .
This is the horizon temperature of the FRW model of the Universe after quantum corrections are taken into account. The quantum corrections to this modified temperature is very similar to that of Banerjee et al. [36] for general static BHs. In addition, the modified temperature has arbitrariness due to the choice of the parameters α k in the expression. Banerjee et al. [36] has shown that, for static BHs, different choices of α k ’s will lead to different physical interpretation. For future work, we shall attempt to study particles with non-zero mass and examine whether quantum corrections as well as Hawking-like temperature depend on the mass term.

4. Hamilton–Jacobi Method in the Lemaitre–Tolman–Bondi Model

The inhomogeneous spherically symmetric LTB spacetime model is given by the matric ansatz in a comoving frame as
d s 2 = d t 2 + R 2 1 + f ( r ) d r 2 + R 2 ( d θ 2 + s i n 2 θ d ϕ 2 ) ,
where R = R ( r , t ) is the area radius of the spherical surface and the curvature scalar f ( r ) classifies the spacetime as follows:
(a)
bounded: 1 < f ( r ) < 0 ,
(b)
marginally bounded: f ( r ) = 0 ,
(c)
unbounded: f ( r ) > 0 .
The Einstein field equations for the spacetime model can be written as: [49,50,51,52]
8 π G ρ = F ( r , t ) R 2 R a n d 8 π G p = F ˙ ( r , t ) R 2 R ˙ ,
and the evolution equation for R is given by
2 R R ¨ + R ˙ 2 + 8 π G p R 2 = f ( r ) .
Hence, the mass function
F ( r , t ) = R ( R ˙ 2 f ( r ) )
is related to the mass contained within the comoving radius “r”. The Universe is assumed to be filled with perfect fluid having energy-momentum tensor
T μ ν = ( ρ + p ) u μ u ν + p g μ ν ,
where the fluid 4-velocity u μ is normalized by u μ u ν = 1 and ρ, p are the energy density and pressure of the fluid, respectively. The energy-momentum conservation relation T μ ; ν ν = 0 gives
ρ ˙ + 3 H ( ρ + p ) = 0 , p = 0 ,
with H = 1 3 R ˙ R + 2 R ˙ R as the Hubble parameter; “dot” and “dash” over any quantity stand for differentiation with respect to “t” and “r”, respectively. Now, splitting the metric ansatz (30) on the surface of the two-sphere and on the 2D hypersurface normal to the two-sphere as
d s 2 = h a b d x a d x b + R 2 d Ω 2 2 ,
the dynamical apparent horizon is characterized by [25,26]
h a b a R b R = 0 ,
where h a b = d i a g [ 1 , R 2 1 + f ( r ) ] is the 2D metric normal to the two-sphere. Thus, the spherical surface of radius R = R A corresponding to the apparent horizon (i.e., marginally trapped surface) satisfies
R A = F ( r , t ) a n d R ˙ A 2 = 1 + f ( r ) .
For the sake of completeness, we recall that, in principle, the energy density, see Equation (35), could be phantom dark energy. In fact, the current observational data does not rule out this possibility [53,54]. For example, in the model in [53], which is founded on a brane–antibrane system, the equation of state parameter of the 4D Universe can change due to owing energy from an extra dimension and decrease from higher values of −1 (non-phantom phase) to lower values (phantom one). In general, we have some interesting properties about entropy and temperature concerning the phantom dark energy scenario. On one hand, as the internal phantom energy is negative, whereas the phantom temperature is definitely negative, and then hotter than any other sources in the Universe, its entropy is always positive, even though holding of the second law is not guaranteed by quantum-mechanical reasons [54]. On the other hand, concerning the observable matter in the Universe, the cosmic phantom field can be regarded as a cosmological source of information and negative entropy [54].
Trapping horizon ( R T ), a hypersurface foliated by marginal spheres, is characterized by [55]:
+ R T = 0 , i . e . , R ˙ T = 1 + f ( r ) ,
± = 2 t 1 + f ( r ) R r ,
are the null vectors to the two-sphere. Thus, both the horizons coincide for the LTB model as in the FRW spacetime.
The Misner–Sharp gravitational mass (in units of G = 1 ) is defined as: [55,56]
m ( r , t ) = R 2 ( 1 h a b a R b R ) .
One might note that this mass “m” is an invariant quantity on the 2D hypersurface normal to the two-sphere and m = m H = R A 2 on the horizon. Furthermore, one can introduce another invariant scalar associated with the horizon in the normal hypersurface known as dynamic surface gravity which is defined as [28]:
κ D = 1 2 R | H = 1 2 h a ( h h a b b R ) | H ,
which, for the given model, has the expression
κ D = 1 2 R t ( R ˙ R ) + 1 2 f ( 2 ) H .
For the LTB model with a decomposed metric (33), the Kodama vector κ has components
κ a ( r , t ) = 1 h ϵ a b b R , κ θ = 0 = κ ϕ ,
i.e.,
κ μ = 1 + f ( r ) , R ˙ R 1 + f ( r ) , 0 , 0 .
Thus, | | κ μ | | 2 = R ˙ 2 1 f ( r ) , i.e., the Kodama vector is time like, null or space like for inside, on the surface or outside the apparent (i.e., trapping) horizon, respectively. Note that the Kodama vector is very similar to the time like Killing vector for stationary BH spacetime and, consequently, an invariant energy associated with a particle is defined by the scalar on the normal space as [57]:
ω = κ a a I .
It should be noted that energy in general relativity is not an invariant. This is a consequence of the equivalence principle [58], which states that one can always choose a coordinate system, i.e., local Lorentz coordinate system, where the gravitational field is null. Thus, the gravitational energy cannot be localized [58]. One measures a particle’s energy, with respect to some natural timelike vector field. However, in Hawking’s black hole radiation, one should not always use the natural timelike Killing vector for a stationary black hole.
Here, the classical action I of the massless particle satisfies the HJ equation:
h a b a I b I = 0 , i . e . , I t 2 1 + f ( r ) R 2 I r 2 = 0 .
Hence, solving (47), we have
I r = ω R 1 + f ( r ) R ˙ 1 + f ( r ) a n d I t = ω R ˙ 1 + f ( r ) ,
where there is a pole at the horizon. The full classical action of an outgoing massless particle is
I = γ a I d x a ,
where γ is an oriented curve with positive orientation along the increasing values of x a =(t,r). As for massless particles, the radial motion is along a null direction, so from the metric, we have
0 = d s 2 = d t 2 + R 2 1 + f ( r ) d r 2 , i . e . , d t = ± R 1 + f ( r ) d r ,
for outgoing and ingoing particles, respectively. Thus, using Equation (48) in Equation (49), we have for outgoing particles
I = 2 d r ( r I ) = 2 γ ω R 1 + f ( r ) d r 1 R ˙ 1 + f ( r ) .
Now, expanding G ( r , t ) = 1 R ˙ 1 + f ( r ) in the neighbourhood of the horizon along a null direction, we obtain
G ( r , t ) R ¨ 1 + f Δ t R ˙ 1 + f R ˙ f ( r ) 2 ( 1 + f ) 3 2 Δ r H + · · · , = R ¨ R ( 1 + f ) R ˙ 1 + f + R ˙ f ( r ) 2 ( 1 + f ) 3 2 Δ r | H + · · · , = 2 R ( 1 + f ) κ D ( r r H ) + · · · .
Then, from Equation (51),
I = γ d r ω κ D ( r r H i 0 ) ,
which has a simple pole at r = r H . Using Feynmann’s i ϵ prescription, the imaginary part of the action (the real part has no physical consequence) can be written as [52]
I m I = π ω H κ D .
Hence, one might interpret T = κ D 2 π > 0 as the dynamical temperature associated with LTB spacetimes.
Moreover, we consider LTB tunnelling computation in the coordinate system ( r ˜ ,t,θ, ϕ ), where r ˜ = R . The metric ansatz (in γ gauge) becomes
d s 2 = A d t 2 2 B d r ˜ d t 2 + C ( d r ˜ ) 2 + r ˜ 2 d Ω 2 2 = q a b d x a d x b + r ˜ 2 d Ω 2 2 ,
where
A = 1 R ˙ 2 1 + f ( r ) = A ( r ˜ , t ) , B = R ˙ 1 + f ( r ) = B ( r ˜ , t ) , C = 1 1 + f ( r ) = C ( r ˜ , t )
.
The horizon is located at χ = 0 , where
χ = q a b a r ˜ b r ˜ = q r ˜ r ˜ = A A C + B 2 ,
i.e., A H = 0 gives the horizon. The Kodama vector is given by
κ μ = 1 B 2 + A C , 0 ^ ,
and the invariant energy has the expression
ω = t I B 2 + A C .
The dynamical surface gravity on the horizon is
κ D = 1 2 B 3 A B + 1 2 A ˙ C | H ,
where “dash” and “dot” represent partial derivatives with respect to r ˜ and t, respectively. The HJ equation for a massless particle along a radial trajectory has the explicit form
C ( t I ) 2 2 B ( t I ) ( r ˜ I ) + A ( r ˜ I ) 2 = 0 .
Now, integration over the temporal coordinate gives real contribution to the particle action (which has no physical significance). An imaginary contribution comes only from integration along the radial direction. Now, eliminating t I between Equations (58) and (60), we obtain
r ˜ I = ω ( B + B 2 + A C ) A B 2 + A C .
For the metric (55), by making a null expansion on the horizon, we obtain
0 = d s 2 = 2 d t d r ˜ + d R ˜ 2 1 + f ( r ) ,
and
1 R ˙ 2 1 + f ( r ) ) 2 R ˙ R ¨ 1 + f ( r ) d t 2 R ˙ R ˙ 1 + f ( r ) f ( r ) R ˙ 2 ( 1 + f ( r ) ) 2 d r ˜ = R ˙ R ¨ ( 1 + f ( r ) ) 3 2 + 2 R ˙ R ˙ 1 + f ( r ) f ( r ) R ˙ 2 ( 1 + f ( r ) ) 2 d r ˜ = 2 κ D H ( r ˜ r ˜ H ) .
Hence,
I m I = γ d r ˜ ω ( B + B 2 + A C ) 2 κ D H ( r ˜ r ˜ H i 0 ) = π ω H κ D H .
Thus, we have the same result as in Equation (54).

5. Determination of Entropy: Area Law

In this section, we shall examine the semiclassical Bekenstein–Hawking area law [1,2,18] for the non-static FRW spacetime at the apparent horizon. In addition, we shall calculate the corrections to the semiclassical entropy due to quantum effects of the Hawking-like temperature. Thus, we start with the thermodynamical law that expresses the energy conservation as
d M = T h d S A .
Here, T h is the temperature of the horizon (with quantum corrections), S A is the entropy of the horizon, and M is chosen to be the Misner–Sharp gravitational mass, defined as [59]
d M = d R A 2
on the horizon. Now, if we identify M in the quantum corrections (29) for the temperature as the above Misner–Sharp mass, then integrating (65), we have
S A = d M T h = 4 π M 1 + Σ α k M 2 k d M = 2 π M 2 + 4 π α 1 l n M 4 π α 2 M 2 + O ( 2 ) = 1 2 A 4 + 4 π α 1 l n M 16 π α 2 R A 2 + O ( 2 ) .
In the above equation, the first term on the right hand side is the semiclassical relation between area and entropy, which is the usual Bekenstein–Hawking area law with a discrepancy of factor 1 2 . The second term is the leading order quantum correction term and is the standard logarithmic correction term in black hole thermodynamics [45,60,61,62]. The higher order correction terms are in the inverse powers of the area of the horizon.

6. Tunnelling Approach

The basic idea in tunnelling method is that following the standard approach, i.e., in the semiclassical approximation (WKB approximation), the emission rate of tunnelling of a massless particle across the horizon can be related to the imaginary part of the action of the system. In the s-wave approximation, the particles are considered as massless shells, moving along a radial null geodesic.
Thus, for the metric (55), the radial null geodesic is characterized by
r ˜ ˙ = B ± B 2 + A C C ,
where as before ± sign indicates an outgoing or incoming null geodesic. Due to the tunnelling of the particles from the outside to the inside of the horizon, here we shall consider only an incoming geodesic. As we have seen that in the present context we need the imaginary part of the action produced by the tunnelling particles (the remaining part is always real), i.e., particles tunnelling through a barrier (the classically forbidden region), thus we obtain [27]:
I m S = I m r ˜ i r ˜ f p r ˜ d r ˜ = I m r ˜ c r ˜ f 0 p ˜ r d p r ˜ d r ˜ .
Here, the particle with radial momentum p r ˜ tunnels from the initial position r ˜ i just outside the horizon to the final point at r ˜ f which is a classical turning point, i.e., in the semiclassical analysis, the trajectory can represent a classically allowed motion
r ˜ ˙ = d H ˜ d p r ˜ | r ˜ ,
where the Hamiltonian H ˜ is the generator of the cosmic time “t”. Then, using Equations (68) and (70) in Equation (69), we write [27]
I m S = I m r ˜ i r ˜ f d r ˜ d H ˜ r ˜ ˙ = I m r ˜ i r ˜ f d r ˜ r ˜ ˙ ω B 2 + A C = ω I m r ˜ i r ˜ f C B 2 + A C { B 2 + A C B } = π r ˜ H ω .
One may note that in order to perform the integration over the Hamiltonian H, we get the energy as ω B 2 + A C , as measured by an observer with the Kodama vector. Then, using the emission rate, i.e., Γ e x p { 2 I m S } , we have
T = ω 2 I m S = 1 2 π r ˜ H ,
which is the HT as derived previously in Section 2.

7. Conclusions

In this work, we have studied Hawking-like radiation from the homogeneous FRW model and the inhomogeneous LTB model of the Universe using both the radial null geodesic method (tunnelling approach) and the HJ formalism approaches. In both of the methods, we have obtained identical Hawking-like temperature at the semiclassical level. “The factor of two” problem in the tunnelling approach has been overcome by considering the Kodama vector instead of time-like vectors and associated energy as the energy of the tunnelling particle. We have obtained quantum corrections to the Hawking-like temperature using the HJ method, and these corrections are similar to those for general static BHs. In addition, the quantum corrected entropy formula has been evaluated from the law of thermodynamics, and it is found that the usual entropy law with quantum corrections are the same as in BH thermodynamics with some discrepancies in the multiplicative factor. In the LTB model, the HT is measured by an observer with the Kodama vector (Equation (45)) inside the trapping (apparent) horizon. In both of the models, the trapping horizon of the Universe is related to the HT in contrast to the BH case, where HR is associated with an event horizon. The present work supports that of the FRW Universe [21] in the literature. In the LTB model, making a coordinate transformation (given by metric Equation (55)), we reach the same conclusion, indicating that HR is independent of coordinate choice. For future work, we shall attempt to find interpretation of the parameters involved in the quantum corrections. It is worthy to generalize the tunnelling approach for non-zero mass particles (i.e., time-like geodesics) as well as to incorporate the quantum correction terms (back reaction effects). Finally, we shall try to resolve the ambiguity in the multiplicative factors for the entropy-area formula.

Acknowledgments

The authors thank the unknown reviewers for useful comments, which permitted improvements to this paper. Christian Corda is funded by the Research Institute for Astronomy and Astrophysics of Maragha, Iran. Subenoy Chakraborty and Subhajit Saha are thankful to Inter-University Centre for Astronomy and Astrophysics, Pune, India for their warm hospitality and research facilities, as the work was completed there during a visit. In addition, Subenoy Chakraborty acknowledges the UGC-DRS Programme in the Department of Mathematics, Jadavpur University. Subhajit Saha has been supported by the Science and Engineering Research Board, Govt. of India under the National Post-doctoral Fellowship Scheme [File No. PDF/2015/000906].

Author Contributions

The Authors declare that each of them contributed in equal measure in realizing this paper. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix. Calculation of Kodama Vector and Surface Gravity

In a time-dependent spacetime, there is no (asymptotically timelike) Killing vector to define a preferred time coordinate. Kodama came forward with a divergence free vector field for any time-dependent spherically symmetric spacetime. The Kodama vector lies on ( 1 + 1 )-dimensional radial-temporal plane and is defined as
K a = ϵ a b b R ,
where ϵ a b is the ( 1 + 1 )-dimensional Levi–Civita tensor in the radial-temporal plane. It is easy to see that a K a = 0 , i.e., Kodama vector is free.
Furthermore, with the Kodama vector, there is an associated conserved current, namely,
J a = G a b K b ,
with a J a = 0 .
For the FRW model of spacetime with metrics given by Equation (1), the Kodama vector is given by
K b = a t b + H R r b .
Similarly, for the metric (2), the Kodama vector is given by Equation (3).
To find surface gravity, we write down the FRW metric (1) as
d s 2 = h a b ( x a ) d x a d x b + R 2 d Ω 2 2 ,
where a , b can take values 0 and 1. The two-dimensional metric
d h 2 = h a b ( x a ) d x a d x b ,
with h a b = diag { 1 , a 2 1 κ r 2 } is referred to as the normal metric with x 0 = t , x 1 = r . Let us consider the following scalar on the normal space
χ ( x ) = h a b a R b R = 1 H 2 + κ a 2 R 2 .
Then, the surface gravity on the apparent horizon is defined as
κ A = 1 2 χ R | R = R A = 1 R A .
Hence, the Hawking temperature on the apparent horizon is given by
T A = 1 2 π R A ,
as shown in Equation (8). Clearly, the surface gravity on the apparent horizon is constant.

References

  1. Hawking, S.W. Particle creation by black holes. Commun. Math. Phys. 1975, 43, 199–220. [Google Scholar] [CrossRef]
  2. Hawking, S.W. Black hole explosions? Nature 1974, 248, 30–31. [Google Scholar] [CrossRef]
  3. Wald, R.M. Quantum Field Theory in Curved Spacetime and Black Hole Thermodynamics; University of Chicago Press: Chicago, IL, USA, 1994. [Google Scholar]
  4. Hartle, J.B.; Hawking, S.W. Path-integral derivation of black-hole radiance. Phys. Rev. D 1976, 13, 2188. [Google Scholar] [CrossRef]
  5. Parikh, M.K.; Wilczek, F. Hawking Radiation as Tunneling. Phys. Rev. Lett. 2000, 85, 5024. [Google Scholar] [CrossRef] [PubMed]
  6. Parikh, M.K. Membrane Horizons: The Black Holes’ New Clothes. Ph.D. Thesis, Princeton University, Princeton, NJ, USA, 1998. [Google Scholar]
  7. Srinivasan, K.; Padmanabhan, T. Particle production and complex path analysis. Phys. Rev. D 1999, 60, 024007. [Google Scholar] [CrossRef]
  8. Sankarnarayanan, S.; Padmanabhan, T.; Srinivasan, K. Hawking radiation in different coordinate settings: Complex paths approach. Class. Quantum Grav. 2002, 19. [Google Scholar] [CrossRef]
  9. Corda, C. Non-strictly black body spectrum from the tunnelling mechanism. Ann. Phys. 2013, 337, 49–54. [Google Scholar] [CrossRef]
  10. Chakraborty, S.; Saha, S.; Corda, C. Quantum Corrected Non-Thermal Radiation Spectrum from the Tunnelling Mechanism. Galaxies 2015, 3, 103–112. [Google Scholar] [CrossRef]
  11. Corda, C.; Chakraborty, S.; Saha, S. Light from black holes and uncertainty in quantum gravity. EJTP 2015, 12, 107–116. [Google Scholar]
  12. Parikh, M.K. New coordinates for de Sitter space and de Sitter radiation. Phys. Lett. B 2002, 546, 189–195. [Google Scholar] [CrossRef]
  13. Parikh, M.K. A Secret Tunnel Through The Horizon. Gen. Relativ. Grav. 2004, 36, 2419–2422. [Google Scholar] [CrossRef]
  14. Parikh, M.K. Energy Conservation and Hawking Radiation. 2004; arXiv:hep-th/0402166. [Google Scholar]
  15. Kraus, P.; Wilczek, F. Self-Interaction Correction to Black Hole Radiance. Nucl. Phys. B 1995, 433, 403–420. [Google Scholar] [CrossRef]
  16. Kraus, P.; Keski-Vakkuri, E. Microcanonical D-Branes and Back Reaction. Nucl. Phys. B 1997, 491, 249–262. [Google Scholar]
  17. Corda, C. Time-dependent Schrodinger equation for black hole evaporation: No information loss. Ann. Phys. 2015, 353, 71–82. [Google Scholar] [CrossRef]
  18. Bekenstein, J.D. Black Holes and Entropy. Phys. Rev. D 1973, 7, 2333. [Google Scholar] [CrossRef]
  19. Jacobson, T. Thermodynamics of Spacetime: The Einstein Equation of State. Phys. Rev. Lett. 1995, 75, 1260. [Google Scholar] [CrossRef] [PubMed]
  20. Cai, R.G.; Kim, S.P. First Law of Thermodynamics and Friedmann Equations of Friedmann–Robertson– Walker Universe. J. High Energy Phys. 2005, 2005. [Google Scholar] [CrossRef]
  21. Akbar, M.; Cai, R.G. Friedmann equations of FRW universe in scalar-tensor gravity, f (R) gravity and first law of thermodynamics. Phys. Lett. B 2006, 635, 7–10. [Google Scholar] [CrossRef]
  22. Eling, C.; Guedens, R.; Jacobson, T. Nonequilibrium Thermodynamics of Spacetime. Phys. Rev. Lett. 2006, 96, 121301. [Google Scholar] [CrossRef] [PubMed]
  23. Akbar, M.; Cai, R.G. Thermodynamic behavior of field equations for f (R) gravity. Phys. Lett. B 2007, 648, 243–248. [Google Scholar] [CrossRef]
  24. Corda, C. Precise model of Hawking radiation from the tunnelling mechanism. Class. Quant. Grav. 2015, 32, 195007. [Google Scholar] [CrossRef]
  25. Corda, C. Black hole quantum spectrum. Eur. Phys. J. C 2013, 73. [Google Scholar] [CrossRef]
  26. Corda, C. Effective temperature, Hawking radiation and quasinormal modes. Int. J. Mod. Phys. D 2012, 21, 1242023. [Google Scholar] [CrossRef]
  27. Bohr, N. XXXVII. On the constitution of atoms and molecules. Philos. Mag. 1913, 26, 476–502. [Google Scholar] [CrossRef]
  28. Hayward, S.A.; Criscienzo, R.D.; Vanzo, L.; Nadalini, M.; Zerbini, S. Local Hawking temperature for dynamical black holes. Class. Quantum Grav. 2009, 26, 062001. [Google Scholar] [CrossRef]
  29. Cai, R.G.; Cao, L.M.; Hu, Y.P. Hawking radiation of an apparent horizon in a FRW universe. Class. Quantum Grav. 2009, 26, 155018. [Google Scholar] [CrossRef]
  30. Kodama, H. Conserved Energy Flux for the Spherically Symmetric System and the Backreaction Problem in the Black Hole Evaporation. Prog. Theor. Phys. 1980, 63, 1217–1228. [Google Scholar] [CrossRef]
  31. Minamitsuji, M.; Sasaki, M. Local conservation law and dark radiation in cosmological braneworld. Phys. Rev. D 2004, 70, 044021. [Google Scholar] [CrossRef]
  32. Racz, I. On the use of the Kodama vector field in spherically symmetric dynamical problems. Class. Quantum Grav. 2006, 23. [Google Scholar] [CrossRef]
  33. Medved, A.J.M. Radiation via tunneling from a de Sitter cosmological horizon. Phys. Rev. D 2002, 66, 124009. [Google Scholar] [CrossRef]
  34. Painlevé, P. La mécanique classique et la théorie de la relativité. C. R. Acad. Sci. 1921, 173, 677–680. (In French) [Google Scholar]
  35. Zhu, T. Hawking-like radiation as tunneling from the apparent horizon in an FRW universe. Int. J. Mod. Phys. D 2010, 19, 159–169. [Google Scholar] [CrossRef]
  36. Banerjee, R.; Majhi, B.R. Quantum Tunneling Beyond Semiclassical Approximation. J. High Energy Phys. 2008, 6. [Google Scholar] [CrossRef]
  37. Akhmedov, E.T.; Akhmedova, V.; Singleton, D. Hawking temperature in the tunneling picture. Phys. Lett. B 2006, 642, 124–128. [Google Scholar] [CrossRef]
  38. Akhmedov, E.T.; Akhmedova, V.; Pilling, T.; Singleton, D. Thermal radiation of various gravitational backgrounds. Int. J. Mod. Phys. A 2005, 22, 1705–1715. [Google Scholar] [CrossRef]
  39. Pilling, T. Black hole thermodynamics and the factor of 2 problem. Phys. Lett. B 2008, 660, 402–406. [Google Scholar] [CrossRef]
  40. Nakamura, T. Factor Two Discrepancy of Hawking Radiation Temperature. 2007; arXiv:0706.2916. [Google Scholar]
  41. Angheben, M.; Nadalini, M.; Vanzo, L.; Zerbini, S. Hawking Radiation as Tunneling for Extremal and Rotating Black Holes. J. High Energy Phys. 2005, 5. [Google Scholar] [CrossRef]
  42. Kerner, R.; Mann, R.B. Tunnelling, temperature, and Taub-NUT black holes. Phys. Rev. D 2006, 73, 104010. [Google Scholar] [CrossRef]
  43. Yale, A. There are no quantum corrections to the Hawking temperature via tunneling from a fixed background. Eur. Phys. J. C 2011, 71. [Google Scholar] [CrossRef]
  44. Singleton, D.; Vagenas, E.C.; Zhu, T. Self-similarity, conservation of entropy/bits and the black hole information puzzle. J. High Energy Phys. 2014, 1405. [Google Scholar] [CrossRef]
  45. Fursaev, D.N. Temperature and entropy of a quantum black hole and conformal anomaly. Phys. Rev. D 1995, 51. [Google Scholar] [CrossRef]
  46. Akhmedov, E.T.; Pilling, T.; Singleton, D. Subtleties in the quasi-classical calculation of Hawking radiation. Int. J. Mod. Phys. D 2008, 17, 2453–2458. [Google Scholar] [CrossRef]
  47. Akhmedova, V.; Pilling, T.; de Gill, A.; Singleton, D. Temporal contribution to gravitational WKB-like calculations. Phys. Lett. B 2008, 666, 269–271. [Google Scholar] [CrossRef]
  48. Akhmedova, V.; Pilling, T.; de Gill, A.; Singleton, D. Comments on anomaly versus WKB/tunneling methods for calculating Unruh radiation. Phys. Lett. B 2009, 673, 227–231. [Google Scholar] [CrossRef]
  49. Joshi, P.S. Global Aspects in Gravitation and Cosmology; Oxford University Press: Oxford, UK, 1979. [Google Scholar]
  50. Banerjee, A.; Debnath, U.; Chakraborty, S. Naked singularities in higher dimensional gravitational collapse. Int. J. Mod. Phys. D 2004, 12, 1255–1263. [Google Scholar] [CrossRef]
  51. Debnath, U.; Chakraborty, S. Gravitational collapse in higher dimensional space-time. Gen. Relativ. Grav. 2004, 36, 1243–1253. [Google Scholar] [CrossRef]
  52. Debnath, U.; Chakraborty, S. The study of gravitational collapse model in higher dimensional spacetime. Mod. Phys. Lett. A 2003, 18, 1265–1271. [Google Scholar] [CrossRef]
  53. Setare, M.R.; Sepehri, A.; Kamali, V. Constructing warm inflationary model in brane-antibrane system. Phys. Lett. B 2014, 735, 84–89. [Google Scholar] [CrossRef]
  54. González-Díaz, P.F.; Sigüenza, C.L. Phantom thermodynamics. Nucl. Phys. B 2004, 697, 363–386. [Google Scholar] [CrossRef]
  55. Cai, R.G.; Cao, L.M. Unified First Law and Thermodynamics of Apparent Horizon in FRW Universe. Phys. Rev. D 2007, 75, 064008. [Google Scholar] [CrossRef]
  56. Misner, C.W.; Sharp, D.H. Relativistic Equations for Adiabatic, Spherically Symmetric Gravitational Collapse. Phys. Rev. 1964, 136. [Google Scholar] [CrossRef]
  57. Hayward, S.A. Gravitational energy in spherical symmetry. Phys. Rev. D 1996, 53. [Google Scholar] [CrossRef]
  58. Misner, C.W.; Thorne, K.S.; Wheeler, J.A. Gravitation; W.H. Freeman and Company: London, UK, 1973. [Google Scholar]
  59. Criscienzo, R.D.; Hayward, S.A.; Nadalini, M.; Vanzo, L.; Zerbini, S. Hamilton-Jacobi tunneling method for dynamical horizons in different coordinate gauges. Class. Quantum Grav. 2010, 27, 015006. [Google Scholar] [CrossRef]
  60. Mann, R.B.; Solodukhin, S.N. Universality of Quantum Entropy for Extreme Black Holes. Nucl. Phys. B 1998, 523, 293–307. [Google Scholar] [CrossRef]
  61. Kaul, R.K.; Mazumder, P. Logarithmic Correction to the Bekenstein-Hawking Entropy. Phys. Rev. Lett. 2000, 84. [Google Scholar] [CrossRef] [PubMed]
  62. Ghosh, A.; Mitra, P. An improved estimate of black hole entropy in the quantum geometry approach. Phys. Lett. B 2005, 616, 114–117. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Chakraborty, S.; Saha, S.; Corda, C. Hawking-Like Radiation from the Trapping Horizon of Both Homogeneous and Inhomogeneous Spherically Symmetric Spacetime Model of the Universe. Entropy 2016, 18, 287. https://doi.org/10.3390/e18080287

AMA Style

Chakraborty S, Saha S, Corda C. Hawking-Like Radiation from the Trapping Horizon of Both Homogeneous and Inhomogeneous Spherically Symmetric Spacetime Model of the Universe. Entropy. 2016; 18(8):287. https://doi.org/10.3390/e18080287

Chicago/Turabian Style

Chakraborty, Subenoy, Subhajit Saha, and Christian Corda. 2016. "Hawking-Like Radiation from the Trapping Horizon of Both Homogeneous and Inhomogeneous Spherically Symmetric Spacetime Model of the Universe" Entropy 18, no. 8: 287. https://doi.org/10.3390/e18080287

APA Style

Chakraborty, S., Saha, S., & Corda, C. (2016). Hawking-Like Radiation from the Trapping Horizon of Both Homogeneous and Inhomogeneous Spherically Symmetric Spacetime Model of the Universe. Entropy, 18(8), 287. https://doi.org/10.3390/e18080287

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop