Particle creation rate for dynamical black holes

We present the particle creation probability rate around a general black hole as an outcome of quantum fluctuations. Using the uncertainty principle for these fluctuation, we derive a new ultraviolet frequency cutoff for the radiation spectrum of a dynamical black hole. Using this frequency cutoff, we define the probability creation rate function for such black holes. We consider a dynamical Vaidya model and calculate the probability creation rate for this case when its horizon is in a slowly evolving phase. Our results show that one can expect the usual Hawking radiation emission process in the case of a dynamical black hole when it has a slowly evolving horizon. Moreover, calculating the probability rate for a dynamical black hole gives a measure of when Hawking radiation can be killed off by an incoming flux of matter or radiation. Our result strictly suggests that we have to revise the Hawking radiation expectation for primordial black holes that have grown substantially since they were created in the early universe. We also infer that this frequency cut off can be a parameter that shows the primordial black hole growth at the emission moment.


Introduction
Most calculations of the particle creation rate from black holes assume a stationary geometry. However, in an astrophysical context, matter and radiation will be falling into the black hole [1][2][3][4][5][6]. This has various consequences. Firstly, the usual black hole uniqueness theorems, which depend on a vacuum condition, will no longer hold, as recently pointed out by Hawking et al. [7]. Secondly, the usual calculations of black body radiation emission need to be amended to take the dynamic geometry into account (they almost all, explicitly or implicitly, rely on the Killing vectors in the vacuum, which are references for defining the mode functions). For the extended maximal geometry including a bifurcate Killing horizon [8], one can derive the Hawking radiation formula according to the Hawking-Hartle vacuum, which uses the Killing vector field to define positive and negative mode functions. These definitions implicitly rely on asymptotic flatness of spacetime, which does not hold in the case of a cosmological black hole [1][2][3][4][5].
As discussed in Refs. [9,10], Hawking radiation will be suppressed in dynamical black holes with matter or radiation accretion. This will be the case until very late times, when the dynamical horizon becomes a slowly evolving horizon. It has been claimed that this suppression effect will be negligible for dynamical black holes, but the argument has been made either on the basis of stationary models, or on the basis of scattering models which do not necessarily apply when infalling (positive density) radiation makes the apparent horizon spacelike [11]. We will model dynamic horizons in this paper, but the calculations are limited to the slowly evolving phase of a black hole. We show that the dynamical effect is negligible in this case (one might assume this is the case on physical grounds, but a calculation is needed to confirm this intuition). But nothing ensures that primordial black holes are slowly evolving, because the background thermal radiation provides the matter flux affecting this evolution and hence affecting Hawking radiation. We prove we can neglect Hawking Radiation in this case by calculating the probability creation rate. However, we do not estimate timescales.
To quantify the particle creation rate for a general black hole, we have to revise the analysis more fundamentally than using methods that rely on Killing vectors. It is well known that particle creation is a quantum effect which comes from quantum fluctuations. The basis of quantum fluctuations is Heisenberg's uncertainty principle, which somehow violates energy conservation at the quantum level (and for short times). The virtual created particle-antiparticle pairs can become real in special environments like strong electric fields (Schwinger effect), special boundaries for the field (Casimir effect), and a black hole event horizon (Hawking radiation). In this paper we derive the nature of particle creation around a general dynamical black hole according to the quantum fluctuations. Along the way, we calculate the probability rate of radiation for the general dynamical case, and give an upper bound for the energy spectrum of the radiation. For dynamical black holes, the spectrum range is smaller than in the non-dynamical case. We show that for a slowly evolving black hole the matter flux can be negligible and the black hole radiation luminosity is very close to that in the Schwarzschild (stationary) case. This is not surprising, but one needs a calculation to prove the black hole radiation is the same in these two cases.
Section 2 reviews the motivation for deriving a general picture of Hawking radiation. The main new results here are, in Sect. 3, we describe Hawking radiation according to particle creation in quantum field theory for a general black hole, generalizing the particle picture in a stationary black hole. Then in Sect. 4, using the Vaidya exact dynamical model, we calculate the probability rate for particle creation for some dynamical black holes with different matter fluxes. In Sects. 5 and 6 we apply this to galactic black holes and primordial black holes, respectively. In the latter case there can be a significant reduction of Hawking radiation.
We do not specifically investigate different Hawking radiation methods or contexts, but rather consider generic constraints on the radiation spectrum that apply in all cases. We illustrate them by some comments on the case of primordial black holes, when the matter flux is not negligible, and where the results might be very different from what is usually expected. To determine whether this is so or not will required detailed calculation of accretion of matter and radi-ation in the context of the early universe. The result is not obvious; we have to obtain clarity by investigating dynamical black holes, not the Schwarzschild case. We show how to approach it in the case of genuinely dynamical black holes.

Black hole particle creation
Most Hawking radiation calculations apply only in the case of stationary geometries (Sect. 2.1), when timelike Killing vectors are present and the apparent and event horizons are identical. However, there are some methods that apply in the dynamic case (Sect. 2.2), when either the apparent horizon is timelike and lies outside the event horizon, which is the case of Hawking radiation back-reaction in an asymptotically flat context, or is spacelike and lies inside the event horizon, which is the case of cosmological black holes (with matter and infalling radiation) that is of interest to us.

Stationary case
The standard way to calculate the particle creation rate is by using the expectation value of the number operator, N = a † a, with an appropriate vacuum. The computation involves Bogoliubov transformations relating the past vacuum and future vacuum states. Quantizing the field, solutions of the classical equations can be written as a linear combination of positive-frequency and negative-frequency parts, corresponding, respectively, to particles and antiparticles.
To do this calculation we need to assume a stationary geometry so that we can define the positive and negative mode functions by using the timelike Killing vectors. In this way, the Euclidean action method based on a path integral [12] and gravitational anomaly methods [13] have been used to calculate Hawking radiation. They both are applicable in the case of a stationary geometry because a key role is played by the exponential relation between an affine parameter and a Killing vector parameter [14, (2.16)]. The observed radiation is related to the event horizon, because it must escape to infinity: we have to locate a boundary so that created particles can get to a distant observer, while trapped particles cannot. In the stationary case, the apparent horizon and event horizon are the same.

The non-stationary case
Consider two rather different non-stationary cases.
Hawking radiation back-reaction in an asymptotically flat context There is a scattering matrix method for determining the thermal spectrum that applies in this situation [15]. One can use this method to determine back-reaction effects on the Schwarzschild metric with fast moving particles. However, it seems to be hard to build the density matrix of cosmological black holes by this method in practice.

Hawking radiation for black holes in a cosmological
context Recently, the tunnelling method was developed to express particle creation events with minimal conditions [19][20][21][22]. In the particle picture of radiation for a stationary black hole, the tunnelling method clarifies that both particle and antiparticle tunnelling contributes to the rate for the Hawking process, where the particle will be outside the horizon and the antiparticle will be inside the horizon. In terms of a classical interpretation, this means the black hole tidal force near the horizon is so strong that it does not let the created particle pair annihilate each other. But this picture is not complete because a quantum particle can be in the classically forbidden region. Also most of the Hawking radiation has wavelengths of order m, and it is only well away from the black hole that one can identify outgoing Hawking "particles" unambiguously. The semi-classical energy-momentum tensor looks like a flux of "real" positive energy particles only at r > 3 m at best, and near the horizon has an ingoing flow of negative energy with no simple particle interpretation. For the most part creation of the Hawking radiation is a highly nonlocal process that cannot be described as taking place in a locally flat region around a single observer. However, one can obtain valid local relations by determining the expectation value of the matter stress tensor at each point.
The best way to describe the quantum particle is by the wave function amplitude, not the tidal force work to create them. This description was used in the tunnelling method calculations of [19][20][21]. Quantitatively, one can describe quantum particle creation by using the matter flux F and energy density E of radiation determined from the stress tensor expectation value T μν [10,35]: This shows pure created matter properties relative to an observer with four-velocity u μ , where n μ is a spacelike normal direction to this velocity, n μ u μ = 0. These quantities are observer dependent and for both an infalling observer and static observer are positive [35] so that it is meaningless to distinguish a place where the Hawking radiation originates [34]. This method shows that the apparent horizon plays a key role in the particle creation process [16][17][18]23]. Therefore, this method inherently can be extended to the case of dynamical black holes. However, it requires an important assumption: that the essential adiabatic condition (eikonal approximation) must be satisfied for the wave equation around the apparent horizon, which is related to redshifting and amplification of waves.
Reference [9] showed the Shankaranarayanan, Padmanabhan and Visser method (wave interpretation) and the Parikh and Wilczek method (particle method) are the same in this context, and are in fact both aspects of the tunnelling method. It was shown that this adiabatic condition gives important limitations on which cosmological black holes can emit Hawking radiation, and described the range of the resulting mass spectrum, depending on the flux. Considering the created matter flux and densities, Padmanabhan [35] confirms our analysis that the matter flux and density is positive outside the horizon for both static and infalling observers. Moreover. we can see created particles with all different wavelengths outside the horizon according to the tunnelling method, but their observed rates are different.
There is another approach which shows by using the energy conservation principle that a dynamical black hole cannot emit radiation, as described in [10]. An observer with four-velocity u μ can measure the particle energy as where p μ is the particle four-momentum. As depicted in the left hand side of Fig. 1, the particle and antiparticle energy, E = −p μ u μ , will both have the same sign at the point B and so energy conservation will be violated around a dynamical black hole if we have pair creation around the dynamical apparent horizon. But a created particle around the point A can have a different sign of energy from the antiparticle (which preserves energy conservation) and so can become a long-lived actual particle. This energy calculation is a classical result. However, we know that energy conservation will be broken in the quantum physics case on small enough scales. From the uncertainty principle (UP) we know that which is the covariant form of the UP presented in Appendix A (Eq. (A.5)). The particle/antiparticle picture can be applied in the dynamic case using this uncertainty principle, which shows how particles can be created around a black hole by vacuum fluctuations. In the next section we will how this can be used to characterise the particle creation phenomenon around a dynamical black hole. We do not claim that created particle must have a special wavelength, but rather give a constraint on the frequency (wavelength) range that is derived from the dynamical nature of the process. We do not use any specific one of the different methods often used, and particularly the tunnelling method, to get our result: we just use the quantum uncertainty principle to give a constraint on the spectrum of Hawking radiation that will hold for all methods.

Particle creation rate for general black holes
In quantum physics, quantum fluctuations allow the creation of particle-antiparticle pairs of virtual particles. These pairs exist for an extremely short time, and then reciprocally annihilate. In some cases, however, it is possible to boost the pair apart using external energy, therefore they avoid annihilation and become actual (long-lived) particles. It was shown in [10] that to have long-lived particle creation (which respects energy conservation) around a black hole, it is necessary that one particle falls into the black hole (that is, inside the apparent horizon) and the second remains outside the apparent horizon. In addition we know that energy conservation can be violated for short times according to the UP. Now we want to see when this principle allows that a virtual particle pair becomes an actual pair around a dynamical black hole.
Consider the right hand side of Fig. 1. At first the black holes radius grows, but finally it becomes isolated [28] at point A. Note that here the event horizon is the causal boundary (a null surface) for the region from which the particle can escape from falling into the black hole. The only case when particle creation can occur is when the antiparticle falls into the apparent horizon so as to have negative energy, and the particle gets to infinity (or to a distant observer) from outside the event horizon. According to the UP (3), if a created particle near the horizon (in the simple case, a photon which moves on a null geodesic at the point B) is supposed to be seen at future infinity and the related antiparticle falls inside the apparent horizon, the particle-antiparticle distance must satisfy where x 0 is the distance between event and apparent horizon seen by a free-falling observer released on the far distance. A dfree observer is defined as a geodesically moving observer who starts from rest at a large distance, which is outside the apparent horizon because we are considering the dynamical case with infalling positive density radiation. Here an observer is defined to be 'at rest' if they are such that the area coordinate R is instantaneously constant for them at the time considered: This is possible when the surfaces R = const are timelike, which will be true at the event horizon or large distance which is outside the apparent horizon in the dynamic case with infalling matter. This observer is in fact a Kodama observer [41]. It is shown in Appendix B that this vector is timelike outside the apparent horizon, even though the apparent horizon itself is spacelike if matter is falling in. Then τ = x 0 is the time that a created particle takes to move from the event horizon to the apparent horizon, in the dfree frame. As it is outside the event horizon, this particle can be seen from outside. Here we assume that the time scale for a free-falling observer to pass between the two horizons is shorter than the time scale of the apparent horizon evolution.
According to the UP, the allowed range of momentum for these created particles is In term of the energy-time UP, we get On defining w := E/h, the relation determines the allowed energy spectrum of radiation from a dynamical black hole in terms of the massless particle frequency [the energy must be positive because of (2)]. To have creation of a real particle as seen by the dfree observer the energy should be in this interval, where all quantities are measured relative to the dfree observer.
There is a similar result for the Schwinger effect [24]. To have actual particles arising from a virtual pair, the probability to find one particle at the distance x 0 from another one is proportional to exp(− x 0 λ ), where λ =h mc is the Compton length for a particle of mass m. In the case that the particles are photons, we have photons created if λ ν > x 0 or w < 1 x 0 , which is similar to Eq. (7). There are two further constraints on the black hole radiation spectrum [25]. The first one comes from the adiabatic condition where κ is the surface gravity defined by [22] in the case of an evolving black hole, which says that the spacetime geometry must be slowly evolving on the time scale set by the frequency of the Hawking photon. The second one is related to a phase space effect: where m extremal is the extremal mass for charged or rotating black holes. This ultraviolet cutoff says that the emitted photon energy can never be greater than the available mass energy.
To calculate these spectral effects on the total luminosity of the black hole, we consider a Planck-distributed flux of Hawking-like radiation which satisfies the adiabatic condition given in Eq. (6.26) of the first reference in [26,27]: namely where κ = −p (u) p(u) with U = p(u) which U and u are affine parameters at future and past infinity, respectively. 1 From [25], when we do not take the uncertainty principle into account, the total luminosity of the black hole radiation can be written in the form with f (w) coming from a grey-body factor which includes backscattering, and the temperature is T =h κ 2π (13) where κ is the surface gravity. Considering the allowed range of the energy spectrum (7), the radiation luminosity around the dynamical horizon becomes where w max = Minimum{w up , w ps }.
Note that Eqs. (9)- (14) all hold in the dfree frame. We can define another useful quantity in this frame: We call P the probability creation rate for a dynamical black hole. This quantity says by what fraction the luminosity will be decreased when considering a dynamical black hole relative to a stationary black hole with the same mass. The effect of the UP (2) is via the quantity w max , which gives an ultraviolet cutoff of the radiation spectrum. It is clear from (10) and (15) Another feature of this creation rate function is that it is proportional toh, which indicates that this is a quantum effect. In the classical limith → 0 there is no particle creation. This quantity is calculated relative to a far distant observer. For this observer, the thermal form of the luminosity is kept. Note that for calculating the quantity w up , we need to calculate the distance between two horizons as seen by a dfree observer on the horizon. For highly dynamical black holes, the distance between the apparent and event horizons becomes greater and w up will be smaller, hence we have less radiation flux: τ 1 ⇒ P 1. In contrast if the radiation influx into the black hole becomes smaller, the distance between the two horizons becomes smaller and w up will be greater. Then τ 0 ⇒ P 1.
Note that a far distant freely falling observer dfree is defined for well-defined backgrounds such as a black hole embedded in asymptotically Minkowski or FLRW spacetimes. If R is again the area distance coordinate and τ is the proper time of the dfree observer, the reference frame dfree is defined as that for which This is uniquely defined because unlike Minkowski spacetime, these spacetimes are not invariant under the Lorentz group (although they are asymptotically flat). One could use Gullstrand-Painlevé coordinates to describe the free-falling observers at large distance. For a more generic backgrounds there is no unique method for defining the reference frame in the curved spacetime.
We fix the frequency relative to this observer dfree by w up = w dfree as the observed frequency. The standard method of deriving Hawking radiation uses the same reference for the thermal radiation frequency. In the Bogoliubov transformation method [14] we choose the observer vacuum at a large distance from Schwarzschild black hole and read the thermal radiation frequency from it. For reading the thermal radiation from the quantum stress tensor we use the same distant observer to read the luminosity of thermal radiation [10]. Since this observer will not see the gravitational time dilation and velocity Doppler shift effects when the particle is released, the intrinsic frequency of the created particle is equal to the frequency measured by this observer at that time.
Now we want to see how the particle goes to infinity. An arbitrary observer near the horizon will find that the created particle frequency near the horizon is given by where k μ is the null particle four vector, H means the horizon position, and u μ is the observer four vector. This frequency will be changed by the gravitational time dilation and velocity Doppler shift relative to the far distant observer. As shown in Fig. 1, different freely falling observers will see different frequencies for the received radiation frequency.
In the spherically symmetric case where the metric take the form let us transform the metric to the areal coordinate, R(t, r ). With this transformation, d R = R dr +Ṙdt, so We know that a free-falling observer does not detect anything abnormal when falling into the black hole. We can transform the time of this metric to the proper time of a free-falling observer, i.e. τ = t + C(t, R). Moreover, we demand that the constant time slice be flat (this means that we choose C so that the coefficient of d R 2 becomes 1). As a result, we get the following form of the metric relative to this free-fall observer: In the Schwarzschild case where a 2 = b −2 = (1−2m/R), the dfree observer has proper time Consequently, we get the Painlevé metric As is clear from this form of the metric, the areal coordinate R is the proper length relative to the proper time coordinate in the radial direction. This quantity is invariant for any time and spatial (radial) coordinate transformation, so different radial observers will see the same areal distance between the two points. The backreacted energy due to Hawking radiation will reduce the black hole mass, determined as the Misner-Sharp mass m = 1 2 R H [21], which quantifies the black hole radius. For the dfree observer, the distance between two horizons can be expressed as The dτ = constant hypersurface in (21) is flat, which leads to (24). In the case that the black hole mass grows from m = m 0 to m iso and becomes isolated, the free-falling observer will see We choose the dfree observer's frame as a standard reference frame for measuring radiation frequencies. But a different observer at different radial point will measure different frequencies. To see this change in the frequencies, we know that the gravitational time dilation and Doppler effect can be calculated by the redshift formula. As shown in Fig. 1 the redshift of the frequency from a point C to point O where the dfree observer is located is Hence for an observer at point C the radiation spectrum is blueshifted with frequency w = (1+z)w dfree . This is like the Schwarzschild case in which the standard Hawking thermal flux is measure by a distant observer at rest. Their result will be the same as that of the dfree obsrver, because their velocities are the same at a large distance and the redshift relation (26) depends only on the four-velocity, not the acceleration (which will be different for them). Any other observer who is nearer will see a blueshifted spectrum [35,36] in the case that the observer starts off from rest.

An example: the Vaidya black hole
To estimate the creation rate function for a dynamical black hole, we apply this formalism to a dynamical Vaidya metric [33] (an exact spherically symmetric solution of the Einstein equations for infalling radiation). In principle, these calculations can be extended to non-spherically symmetric cases. Note that this geometry, which exactly represents infalling radiation, will apply only until Hawking radiation starts and the outgoing flux becomes dominant.

Ingoing Vaidya metric
The Vaidya spacetime metric is (27) and from the Einstein equations has stress-energy tensor where m(v) is the mass and d 2 = dθ 2 + sin 2 θ dφ 2 . The outgoing, , and ingoing null geodesic, n, which are cross normalised, .n = −1, are Their expansions are If v μ = μ − Cn μ is the tangent vector to the apparent horizon, where C is the evolution parameter, we find the following equation for Vaidya black holes: This shows that the parameter C is proportional to the matter flux which falls into the black hole. Assume that the evolution is characterised by a time scale so that where V = v and m i is the initial mass. To describe the near-equilibrium limit which we expect (see [32]) to occur at large (in this case "large" means large relative to m i ), we set For the perturbative calculations (1/L) is a useful expansion parameter.
Defining a new radial coordinate R by the Vaidya metric becomes 4.2 The radiation formula: slowly evolving case As shown in [32], up to second order in (1/L), the distance between the event horizon and apparent horizon can be obtained from This expansion ensures that the black hole horizon is a slowly evolving horizon (SEH) [29]. Here the dot indicates a derivative with respect to V . The distance between the event horizon and apparent horizon for the dfree observer is From (7) the allowed energy spectrum for this black hole is (38) and the ultraviolet frequency for L up is Therefore, by determining this frequency we can determine the frequency spectrum range and then the particle creation rate from (14). Note that Eq. (37) (which characterises the distance between the two horizons in the slowly evolving horizon case) does not work for dynamical black holes. It was derived on the assumption thatṁ 1.

Numerical values: galactic black holes
For a galactic black hole with mass 10 6 M with accretion rate of about dm dv = 10 3 M year and assuming m d 2 m dv 2 dm dv , in SI units so we have In SI units, and Since the w ps > w up , the ultraviolet frequency spectrum will be w min w w max = w up (44) with w min given by (42) and w up by (41). With this frequency spectrum, the numerical value for the probability creation rate (16) of this black hole is This is not strange because the matter fluxṁ 1. To present the horizon phase in terms of geometrical parameters, we need to define an invariant parameter. The C parameter is not invariant under changing the null geodesic parameter, so we have to make an invariant quantity from this function and calculate it for the black hole apparent horizon. This new quantity is and we get for the Vaidya model. The numeric value ofC for a galactic black hole is which verifies that the black hole horizon is a slowly evolving horizon.
As a result, we consider Hawking radiation for dynamical black holes that have slowly evolving horizons. For all values of the matter flux such thatC ∼ṁ 1 and we have a slowly evolving horizon, for the probability creation rate (16) we get It is inconceivable to have cosmological black holes in such an isolated phase. As a result, in practice for cosmological black holes we have The final conclusion is that the usual galactic mass black holes in a cosmological context will be slowly evolving and will emit thermal radiation, but the temperature is below the CMB (cosmic microwave background) thermal temperature so it will not be visible. As an aside, the black hole heat capacity is negative. This means that the greater the black hole mass, the less the Hawking temperature. Therefore we cannot apply the classical thermodynamic scenario to them. Namely, if heat comes from the hotter CMB radiation to the black holes, the CMB flux grows the black hole mass [9] and the black hole temperature will decrease. Hence, the black hole becomes colder. There is not a classical thermodynamic interpretation as in classical general relativity. Black hole radiation connects geometrically to other heat fluxes like the CMB. The assumed accretion rate is much greater than that normally considered astrophysically reasonable, and it is vastly super-Eddington. Therefore, the fact that the dynamic effects are small in this case is evidence that deviations from the stationary black hole Hawking radiation are negligible for any astrophysical black hole in the present universe. We show here that the Hawking radiation luminosity for astrophysical black hole is very close to the stationary black hole luminosity.

Dynamic black holes: primordial black holes
To have a good intuition for the spherical dynamical black hole particle creation rate in the case of dynamic black holes, consider a dynamical black hole such that its initial apparent horizon radius is r i = 2m i where m i is the initial Misner-Sharp mass [31] and it is radius grows (due to the infalling matter flux) till it becomes isolated and its final radius becomes As shown in Eq. (25) the distance between the two horizons relative to the dfree observer is x 0 = 2δm. Hence, the frequency w up can be obtained from Eq. (8).
We define the new parameter n as which shows how much the black hole horizon grows relative to its initial radius. To get the order of magnitude for the probability creation rate (16) we approximate the spectrum as a Planckian spectrum, as is also proposed in the tunnelling method [37] in calculating the luminosity order of magnitude.
Neglecting w min , we find the probability creation rate (16) for different numbers n as shown in Table 1. This table shows an interesting result: that for a dynamical black hole such that its final radius becomes 25 times greater than its initial radius, the luminosity of the particle creation rate dies off by a factor P = 0.18. This result is compatible with the intuition from Fig. 1 that bigger n means bigger growth of the horizon radius, and less created particles around r i at point B can escape to outside of r f at point A.
This result changes the expected primordial black hole radiation when they are located in the hot soup of the very early radiation dominated era. It is probable that primordial black holes that were created in the early universe eat this much matter flux, and their initial radius grows by this amount. The UP frequency visibility, w up , is another observational effect in our analysis. For a Primordial black hole with mass m = 10 15 gr if n = 0.1 or x 0 = 0.1r i the upper bound frequency becomes w up = 10 11 H z which can be observable. Therefore, observing the radiating frequency cutoff one can verify that this black hole was in a dynamical phase at the emission moment.

Conclusion
As this is a quantum effect, to understand the nature of particle creation near a general black hole we need to adopt a radiation scenario that takes into account quantum fluctuations.
The essence of quantum fluctuations is based on Heisenberg's uncertainty principle (UP) at the quantum level. We have used this principle to study the nature of Hawking radiation near a dynamical black hole. Applying the UP allows one to set limits on the energy spectrum of the radiation. The radiation spectrum has a dependency on the distance x 0 between the apparent and event horizons which is measured by the dfree observer. The new proposed spectrum (14) will reduce to the standard spectrum when the geometry becomes stationary. We also defined the probability creation rate (16) to show by how much black hole radiation will be decreased due to the dynamical nature of a black hole. Note that we did not claim that particle creation takes place just on the horizon. We determined the minimum distance between the horizons allowing pairs to be real in the dynamical case. From that we derived an upper bound for the radiation frequency cut off.
To apply this to known black hole models, first we took the Vaidya model where we can explicitly calculate the distance between the two horizons. But it is in an adiabatic phase, and has a slowly evolving horizon. We have shown that the probability creation rate for such a Vaidya black hole which has a slowly evolving horizon is P ∼ 1. This verifies our last work [10], where we inferred that black holes with slowly evolving horizons emit Hawking radiation. This result also is compatible with the thermodynamics principle that says temperature is meaningful only in the case of equilibrium or near equilibrium systems.
To study dynamical models, we need to make numerical simulations. It is suggested by this investigation that a black hole starts emitting radiation when it becomes slowly evolving. But our spectrum discussion concerns general black holes and will apply in that case as well as the dynamic case. In addition, to get some intuition for the probability creation rate of a dynamical black hole, we calculated this quantity for a spherical black hole such that its apparent horizon grew several times due to accretion. As shown in Table 1, the probability creation rate will not be negligible in the case that its horizon grows 10 times greater. This seems to be an important result for primordial black hole radiation. We have inferred that this frequency cut off, w up , can be observable for dynamical primordial black holes, and we can consider it to be a parameter which shows the primordial black hole growth at the emission moment. Since these black holes were created in the radiation dominated era when they were surrounded by a hot soup, these black holes can grow substantially. This paper confirms our last paper's results [9] that we have to revise the primordial black hole radiation formula to take this into account.
The case of a primordial black hole considered here does exhibit a breakdown of the quasi-stationary approximation, but one might assume that Hawking radiation is totally negligible compared with the thermal radiation of the early universe and classical gravitational radiation associated with the black hole formation process. There may be no observational consequences of the result. However, our result challenges the dynamical primordial black hole radiation result which is mentioned in the literature [39,40]. To check this result one needs to quantify the Hawking radiation flux and compare it with the incoming flux of thermal radiation of the early universe. The scope of this paper is not to consider in detail primordial or astrophysical BHs, it is to consider general principles at work. Table 1 is just an example, but is not applied to any specific simulation.
with the weak energy condition where w is constant. The Einstein equations give and M, the Misner-Sharp mass, is defined by It can be shown [31] that if a perfect fluid collapses the apparent horizon will form at R = 2M. Consider a two-sphere with R = constant. We can define timelike and spacelike normal vectors to this compact surface. Consider the radial vector It can thus easily be seen that the Kodama vector outside the apparent horizon (R > 2M) is timelike; inside the apparent horizon (R < 2M) it is spacelike; and on the apparent horizon (R = 2M) it is null. On the other hand, we know that the event horizon is located outside the apparent horizon. As a result, the Kodama vector is timelike on the event horizon.