1. Introduction
Using classical general relativity, it was shown by Oppenheimer and Snyder [
1] that the ultimate fate of a spherically symmetric collapsing star must be a black hole. According to general relativity, classical black hole solutions have singularities arising during gravitational collapse. In particular, Penrose showed that even deviations from spherical symmetry cannot prevent spacetime singularities from arising [
2]. On the other hand, Hawking [
3,
4] used quantum field theory in strong gravitational fields and found that there must be a thermal flux of particle production, which is known as Hawking radiation. This means that a static observer located far away from the black hole should detect the temperature. Such temperature is very small and, as of today, has not been measured. It is widely believed that such spacetime singularities can be cured within a quantum theory of gravity. Regular black holes have attracted a lot of attention (see different regular black hole solutions [
5,
6,
7,
8,
9,
10,
11]), including a recent review [
12] and possible constraints to regular black holes with the Event Horizon Telescope image of Sagittarius A
[
13]. Different concerns have been raised about the stability of regular black holes. Specifically, it has been argued that regular black holes can be generically unstable because of the phenomenon known as the “mass inflation”, which can destabilize the inner horizon and the role of Hawking radiation to cure this instability [
14,
15], and the problem with such a claim see [
16]. In the present paper, we shall peruse a different scenario; namely, it was argued that ideas from T-duality can regularize the gravitational potential [
17,
18,
19,
20] and that this can play an important role to resolve the black hole singularity [
18]. Using T-duality, it is possible to show that the description of string theory below the length
is the same as its description above
. In this framework, the physics of four dimensions can be obtained by compactifying the other dimensions. For a single compact dimension with radius
R, one can use the boundary conditions by writing
where
w is known as the winding number. Furthermore, for the mass spectrum for such a system, we have
where
n is known as the Kaluza–Klein excitation level. The main idea behind T-duality is that the above spectrum does not change if we exchange the winding number
w and the Kaluza–Klein excitation level
n; namely, we can write [
18],
On physical grounds, this also suggests that the description of string theory below a certain length is equivalent to its description above it [
17]. Furthermore, by means of T-duality, one can show that Green’s function is invariant under
. In particular, for the Green’s function in the momentum space, the following was found [
17,
18,
19,
20],
in which
is a modified Bessel function of second kind. The zero-point length (
) is, therefore, produced by means of the compactified extra dimension of radius
R and cannot be probed below this length. Taking the limit
, we obtain the standard relation for the Green’s function, i.e.,
. In such a limit, the stringy effects are very small and can be totally neglected. Using the modified Green’s function, it was also shown that the point-like source distribution is replaced by a smearing matter density. Specifically, using the regularized potential due to the zero-point length
, and by solving the Poisson equation, one can obtain the energy density and the stress-energy tensor describing the smearing matter distribution [
18]. The static and spherically symmetric metric that solves the Einstein field equations with stringy effect is given by [
18]
where
M denotes the Komar mass. This is a very important solution because it is a non-perturbative solution that describes a static and spherically symmetric black hole geometry. For
, there exist two roots: the inner and outer horizon, which are
and
, respectively. We can also say that such a metric describes two possible phases of matter, the particle sector (
) and the black hole sector (
). For a very large mass, the solution is effectively the Schwarzschild black hole. The regularized electrostatic potential in T-duality was obtained in [
21]. Recently, such non-perturbative modifications were used to find a charged black hole solution in T-duality [
22], regular black holes in three dimensions [
23], entropic corrections to Friedmann’s equations [
24], regular black strings, and torus-like black holes [
25]. By considering Hawking evaporation, it was argued that black remnants should occur due to stringy theoretical effects [
26]. In the present work, we would like to study, in more detail, gravitational collapse and the final stage of the collapse using such stringy corrections.
This paper is outlined as follows. In
Section 2, we study the gravitational collapse of the interior of star. In
Section 3, we discuss the Planck star radius using an infalling observer point of view. In
Section 4, we discuss Hawking evaporation and the final Planck-size remnants. In
Section 5, we study the bouncing Planck star scenario. We comment on our results in
Section 6.
2. Gravitational Collapse and Planck Stars in T-Duality
Let us start by considering a spherically symmetric star composed of matter (dust) with vanishing pressure, which undergoes gravitational collapse. In general, the stress-energy tensor of the collapsing matter is given by
in which
is the energy density,
is the fluid 4-velocity, and
p is the pressure, which, in our case, vanishes. For such a fluid, we have to consider the energy conservation, i.e.,
, along with the Einstein equations. For the interior region, having a spherically symmetric star, we can write, in general [
27],
in which
is the area radius. In addition, we take
, which follows from the Einstein equations for the case of homogeneous dust [
27]. To study the gravitational collapse, we are going to use the well-known Tolman–Bondi spacetime by introducing the following function [
28,
29,
30,
31,
32,
33]
which leads to
For the exterior metric, we shall use the modified vacuum solution due to the stringy effects given by Equation (
5) and rewritten as
At this point, we utilize the Misner–Sharp mass function, which is defined by using the area radius, which, at fixed
, reads
from this equation, it follows that
A very important result that follows from the Tolman–Bondi spacetime is that one can obtain the Friedmann equations as a special case. To see this, we need to introduce the following relations
It can be easily seen how the FRW universe metric is obtained
Here,
k denotes the curvature of space, with
corresponding to flat, closed, and open universes, respectively. Using this equivalence, we can model and study the interior spherically symmetric homogeneous stars. Using Equations (12) and (13), we can find
It is important to note here that we shall neglect the backreaction effect of pre-Hawking radiation during the gravitational collapse. In particular, it has been shown that such an effect is small and cannot prevent the formation of the apparent horizon (see for example [
34]). The dynamical apparent horizon, a marginally trapped surface with vanishing expansion, is determined by the relation
where the two-dimensional metric reads
It is a simple calculation to find out the relation for the apparent horizon radius of the FRW universe
We are going to simplify the work; since
is a very small number, we can consider a series expansion around
via
Then, using the Friedmann Equation (
15), we obtain, in leading order terms,
The last equation can be further written as
where
is a constant defined as
This result is nothing but the corrected Friedmann equation reported recently in Ref. [
24] using a different approach (Verlinde’s entropic force scenario). In fact, it coincides with [
24] by taking
(dust). It is quite remarkable that we found a bridge between two different and competing directions in quantum gravity. On one hand, by considering string T-duality effects, we found the modified Friedmann Equation (
21), which coincides with the conclusions obtained from the loop quantum gravity approach [
35]
where
is the critical energy density
where
[
35] is the area gap that sets the discreteness scale of loop quantum gravity, and
is the Immirzi parameter. The correspondence is achieved by identifying
. A direct computation yields
. Using
[
18], we obtain
, which is in perfect agreement with the value proposed in loop quantum gravity
.
Once the gravitational collapse takes place, we can now use the modified Friedmann equations and explore the possibility that the collapse stops at some point due to the stringy corrections. To do so, we have to use the condition
along with
. From this condition, we can obtain the critical density; in fact, we obtain two branches of solution for the critical density
From this result, it follows that
which basically allows us to find the minimal quantity for the scale factor
Again, this is in perfect agreement with what has been found in Ref. [
24]. Such a critical density is, thus, inversely proportional to the minimal length
The above arguments show that during the gravitational collapsing phase, the singularity is never reached, and the interior solution of the black hole (black hole core) is a kind of very dense star. The possibility that a very dense star or Planck star exists inside the black hole was proposed in Ref. [
36]. The radius of such a star was conjectured to be proportional to the collapsed mass [
36]. It is very interesting that, as we shall see, such Planck stars hidden inside the stringy-corrected regular black holes can naturally appear in our analyses. In what follows, we are going to compute the radius of such a star. First, we need to rewrite the FRW metric in a simple form. Let us define the proper time
using
where
is the conformal time, along with radial coordinate defined as
From these equations, we obtain the FRW metric as
Choosing a surface
with fixed
, by matching the metrics, we can obtain the first equation
along with the second equation
From the last equation, it is not difficult to show that
Although we use a rather simple and idealized model of collapse, it highlights the main features of the interior dynamics of the interior of the star. Using the matching procedure of the interior and exterior metrics at the surface of the star, it is possible to study the motion of the star’s surface. In what follows, we shall show some important results; first, we are going to approximate the last equation as
where the plus/minus sign corresponds to the case of expansion or collapse. Because we are interested in the collapse, we chose the minus sign (
decreases with time), and by doing further simplifications, we obtain
where
We see that from the point of view of the outside observer, it takes an infinite amount of time
to see the formation of the black hole horizon
. The whole process is viewed in “very slow motion”. However, from the point of view from the inside, it takes a finite proper time for particles to reach the minimal distance. In the Oppenheimer–Snyder model [
1], the surface of a gravitationally collapsing spherically symmetric star made up of dust with radius
can be obtained via Equation (
33). At this point, let us define the following constant quantity
where
is the scale factor in the initial moment of collapse. We can see that this quantity is constant simply by taking
for the dust matter. At the initial time, we also have
, along with the radius of the star
. Furthermore, one can show that
From the corrected Friedmann equation (setting
), we obtain
or, in terms of
, we obtain
Solving the last two equations, exactly, is not an easy task. One simple guess is to try and generalize the parametric form
, which is a solution when
, then by using the following equation
we obtain
There are two branches in this solution that can describe the contraction and expansion. Again, finding an exact solution in closed form is outside the scope of the present work. Because the mass is conserved during the gravitational collapse (having in mind the Hawking radiation is very small), we must also have
once the Planck star is formed. At the surface of the star when the gravitational collapse stops, we also have
Note here that
. This means that we can obtain the proportionality
where we can identify
. Put in other words, during the gravitational collapse, the scale factor decreases, but the density per unit volume increases. Another way of stating this result is to say the mass of the collapsing matter is constant
When the gravitational collapse stops, we can find the radius of the Planck star using
; namely, we obtain
Using
, we estimate the radius of the Planck star as follows
The above value for the radius is in good agreement with Ref. [
36], having set
. The phenomenological aspects of Planck stars have been studied in Refs. [
37,
38,
39,
40]. For a stellar mass black hole with mass
and
m, we can obtain the radius of the Planck star (by restoring the constants
G and
c)
m. Although a small value, this shows that the radius of such a star is many orders of magnitudes greater compared to
. Such a star is hidden inside the event horizon of the black hole, with the geometry described by the metric (5).
3. A Free-Falling Observer and Planck Star Radius
Let us now study the whole process, as seen from a free-falling observer. To do this, we can use the Painlevé–Gullstrand coordinates through the definition of a new time coordinate as
for some arbitrary function
, along with the new metric
We see that there is no coordinate singularity at the horizon. The time coordinate of the Painlevé–Gullstrand metric is the same as the proper time of a freely falling observer who starts from infinity at zero velocity. We denote the Painlevé–Gullstrand coordinates as
and the Schwarzschild coordinates as
. One can use the Jacobian to relate these coordinates given by [
41]
along with the inverse of the transformation matrix
From the point of view of a static observer located far away from the black hole, the total energy momentum is the sum of the energy momentum of the black hole core or the Planck star energy density, and for the renormalized stress-energy tensor, we add the effect of Hawking radiation. For example, one can choose the Unruh vacuum stat (see [
41]). Hence, we can write
For a freely falling observer, we can write the components in Painlevé–Gullstrand coordinates by means of a coordinate transformation
As we did in the last section, we shall neglect here too the Hawking radiation effect as perceived by a freely falling observer and focus only on
. For simplicity, we work in
dimensions; this yields
with the components of the energy–momentum for the black hole core in Schwarzschild coordinates, given by
For a freely falling observer, the velocity in Painlevé–Gullstrand coordinates is given by
Using this velocity, we find that the energy density as measured by such an observer is given by
In other words, the energy density stays an invariant quantity. At this point, we use the condition
and after solving this equation, we obtain the minimal value at
This is in perfect agreement with the minimal scale factor found in the last section. We can now compute the total time measured by such a free-falling observer using
where we approximate
. After solving this integral, we obtain
The proper time of a particle is, therefore, finite. We can show that the time for light reaching the minimal distance, say, from the event horizon, is also finite. In this case, one can use
to find the radial equations, then using the integrating the equation, we obtain
where, again, we can use the approximation
. After solving this integral, we obtain a finite amount of time
Let us now use Equation (
51) to find the time
t measured by an observer located far away from the black hole
which yields
where
C is an integration constant. In the limit
, we obtain
, meaning that from the observer’s point of view, it takes an infinite amount of time to see the collapsing of matter. Due to the quantum gravitational effect, or the zero-point length effect, we found that the particles never reach the singularity, but this also implies the existence of Plank stars. This can be seen from Einstein field equations and using
, we must have
This shows that there is no singularity in the expression for the Ricci scalar, provided
. One can calculate one more scalar invariant, known as the Kretschmann scalar, given by the following result
To estimate the radius of the Planck star, here we recall that the Ricci scalar (
) for the above black hole given is found
From these two equations, we obtain
considering a series expansion around
, and by setting the radial coordinate to be the Planck star radius (we call it
), we obtain
Solving
, we obtain
which is in perfect agreement with Equation (
50) found in the last section in leading order terms.
5. Bouncing Planck Star: Black-Hole-to-White-Hole Transition
There is another possibility—perhaps a more interesting one—in which the Planck star bounces instead of decreasing its radius. This is due to the fact that the inner core (Planck star) solution may not be a stable state after all. Mathematically, bouncing at the critical point can be stated using the conditions:
,
, along with the condition
. This shows that there is a great level of similarity between the physics that describes the cosmic bounce and the possible bounce inside black holes. One can use the second modified Friedmann equation that describes the dynamical evolution reported in Ref. [
24]
where the minimal scale factor given by [
24]
We see that, in general, if we have matter with non-vanishing pressure, then the scale factor can be a function of
. Imposing the condition
, we obtain the interval
. On the other hand, if we use the equation of state via
, along with
, we obtain
which is further rewritten as
where
provided
. However, we must also have in mind that
; therefore, we are left with the allowed interval
. The general solution in this interval is given by
where we can take the interval
. At the initial moment
, one has
; hence, we can fix the constant
, which yields
. The interior metric now reads
As was argued in [
24], this metric can describe the bouncing universe. For reasons we elaborated above, we need a special form of matter with a specific interval for EoS parameter
in order to justify the bouncing effect. Coming back to our case, where we studied the collapsing of matter (dust) with zero pressure, i.e.,
, along with
, this means that the above bouncing condition is not satisfied. At this point a natural question arises: even if we have collapsing dust, which clearly does not satisfy the above bouncing condition, can we still say that the final state of the internal core of the black hole will be eternally stable? Of course, we do not know the answer to this question, but from a quantum mechanical point of view, we may speculate that the bouncing effect can also be a consequence of the black-hole-to-white-hole transition (BHWH). In other words, instead of the bouncing condition given by Equation (
84), which is a classical effect, we can have a purely quantum mechanical bounce due to the quantum tunneling effect. The idea behind the BHWH transition is not new; for example, in [
42], the authors tried to compute the probability amplitude between two configurations, say
and
, with the corresponding hypersurfaces
and
. In particular, the probability amplitude for the BHWH transition can be computed from [
42]
where
is a parameter measuring the width of the interpolating region. Furthermore, it was estimated for the BHWH transition probability an exponential decay law [
42]
with a mean lifetime
. There are other other arguments about the black-hole-to-white-hole transition. For instance, the probability increases with time if we take into account Hawking radiation (see [
38,
43]). This can be explained from the fact that as the mass of the Planck star decreases with time and the bouncing mass becomes smaller compared to the initial Planck mass, then, in accordance with the semiclassical standard tunneling factor
[
38,
43], the probability for the black-hole-to-white-hole transition increases as the mass decreases. Here, we note that
is the Euclidean action, where
. The tunneling probability per unit time can also be written (here, we restore
ℏ for a moment)
[
38,
43]. Here, let us consider again the Painlevé–Gullstrand coordinates that relate the time measured by an outside observer and the time measured by an outgoing observer, given by
along with the white hole metric
Due to the bouncing effect, the black hole becomes essentially a white hole with an explicit time-reversal symmetry. In that sense, a white hole is a solution in general relativity with a spacetime region that cannot be entered from the outside. From the point of view of an outside observer measuring in Schwarzschild coordinates, the time-reversed solution or the white hole geometry is the same as a black hole. As we saw, it takes a finite proper time to form a Planck star from the gravitational collapse, and yet, from the point of view of outside observer, due to the strong redshift effect, the gravitational collapse appears "frozen" in time due to the formation of the horizon. The same can be shown for the bouncing process. An outside observer sees the collapse/bouncing in “very slow motion”, and the entire process takes a long time. To see this, let us consider a white hole region with a time-reversed solution, i.e.,
in Equation (
89), then the time measured outside the white hole is given by
The first terms is finite proper time measured by the outgoing observer and can be computed via
and the result is similar to Equation (
64). For the time measured by the outside observer, we obtain
meaning that for a particle to reach the event horizon
, we need
. Put in other words, the bouncing effect of the star appears in “very slow motion” when observed from the outside. We can basically deduce the same conclusion using the matching of the interior and exterior metrics. To do such a computation, we need, of course, the explicit form of the scale factor. Let us take, just for fun, the exponential law, i.e.,
, then we obtain
where
; here,
C is some constant. In the initial time of expansion
, we must have
, i.e.,
should coincide with the radius of the Planck star. Now, assuming that during the expansion we reach the classical horizon radius with
, we obtain the total proper time
Because we have an expansion in this case, we need to take the plus sign in the right-hand side of Equation (
36), along with the time-reversed condition
. In doing so, we obtain
Solving this equation for the time leads to the following result
If we replace the expression for the proper time
, we finally obtain
The first term goes like
, and again, this confirms the fact that the time measured from the outside observer will be very large, i.e.,
, which is in agreement with Equation (
93). At this point, one can ask whether white holes can be stable remnants. Authors in [
43] argued that such a unitary process may not violate any known physics. This question is outside the scope of the present work, but if there is surrounding matter, most probably the white holes are unstable objects too and collapse again to black holes. According to [
43], there is a difference in the lifetime between black holes and white holes. The former is described by the law
and the latter by
. If such a spacetime bounce happens, there is a possibility that strings can increase the size. This is similar to the so-called fuzzball structure of the black hole, speculated in Ref. [
44].