We next substitute the right hand side of Equation, (2-27), into the left hand side of relation, (2-28), since these are equivalent. We thereby obtain

${G}_{r}^{\left(4\right)}{M}_{r}^{\left(4\right)}{\rho}_{r}^{\left(4\right)}2{\text{\pi}}^{2}\text{d}r=-\text{d}{p}_{r}^{\left(4\right)}{A}_{r}^{\left(4\right)}-{p}_{r}^{\left(4\right)}6{\text{\pi}}^{2}{r}^{2}\text{d}r$ (2-29)

This we divide by, dr, and rearrange terms. The result is

$\text{d}{p}_{r}^{\left(4\right)}/\text{d}r{A}_{r}^{\left(4\right)}=-{G}_{r}^{\left(4\right)}{M}_{r}^{\left(4\right)}{\rho}_{r}^{\left(4\right)}2{\text{\pi}}^{2}-{p}_{r}^{\left(4\right)}6{\text{\pi}}^{2}{r}^{2}$

or, dividing by 4-D surface area,

$\text{d}{p}_{r}^{\left(4\right)}/\text{d}r=-{G}_{r}^{\left(4\right)}{M}_{r}^{\left(4\right)}{\rho}_{r}^{\left(4\right)}/{r}^{3}-3{p}_{r}^{\left(4\right)}/r$ (2-30)

We can compare this to our original ansatz, Equation, (2-14). Right away, we notice that the second term in Equation, (2-30), is new! Equation, (2-14), is missing this term, and cannot be correct because it does not take into account the change in surface area. Equation, (2-14), assumes that,
${A}_{r}^{\left(4\right)}={A}_{r+\text{d}r}^{\left(4\right)}$, which is patently false. If we retrace our steps, we recognize that
${p}_{r+\text{d}r}^{\left(4\right)}$ acts radially inwards, at radius,
$r+\text{d}r$, and it acts on surface area,
${A}_{r+\text{d}r}^{\left(4\right)}$. The radiative pressure,
${p}_{r}^{\left(4\right)}$, on the other hand, is directed radially out, and it acts on a different surface area,
${A}_{r}^{\left(4\right)}$. The difference in radiative forces in Equation, (2-28), requires us to take into account not only a change in pressure, but also, at the same time, the change in areas, upon which these pressures are acting. In the 4^{th} and 5^{th} lines of Equations, (2-28), we see this very explicitly. The change in area is what leads to the second term on the right hand side, in the last line of Equation, (2-28). In 3-D stellar hydrostatics, apparently, the second term must be negligible when compared to the first term, and for this reason, it is probably ignored. We do not do so here. That second term on the right hand side of Equation, (2-30), is crucial.

3. Choosing a Truncated Gaussian as a PDF for a Radiative Mass Distribution within the Black Hole; Determining the Shape Parameter

As was shown in Section 2, Equations, (2-6), and, (2-13), led to a naïve solution for determining the distribution of radiative mass within the black hole. If the distribution of radiative mass within the black hole were uniform, then a gravitational force within its interior cannot be accommodated. Nor would we be able to keep the inner surface temperature small with such a choice, and still pack the requisite mass within radius, R. We are thus led to consider other pdf’s which will bypass these problems. One good choice for our purposes is a truncated Gaussian [13] - [18] . In this section, we will fix the parameters associated with this distribution, such that we can model the interior for any size black hole. In this section, and the remaining sections, we will typically drop the superscript (4) on variables as a matter of convenience, as all quantities are defined in 4-D space.

First and foremost, we should recognize that other choices [19] [20] [21] [22] for a continuous truncated pdf are possible. They would lead, like the truncated Gaussian we are about to choose, to steep increases in radiative densities and pressures within the interior. Among the possibilities, we could list the truncated exponential distribution, the beta distribution, the Kumaraswamy distribution, etc. These all have support, $r\in \left[0,R\right]$, and satisfy Equations, (2-8), (2-9) and (2-12). Moreover, upon a judicious choice of parameters, these probability distributions can mimic the truncated Gaussian in many ways. The exponential involves only one parameter; the other two listed, require two parameters in order to specify the distribution. We will however choose the Gaussian or Normal distribution, ${f}_{r}=f\left(r,\mu ,\sigma ;a,b\right)=f\left(r,\mu ,\sigma ;0,R\right)$, where $\left(\mu ,\sigma \right)$ are the parameters to be determined, and $\left[a,b\right]=\left[0,R\right]$ are the lower and upper bounds for the r = radius value.

The truncated Gaussian, or truncated Normal, distribution with parameters, $\left(\mu ,\sigma \right)$, will serve our purposes well for the following reasons. We list them in bullet form.

1) The pdf is relatively simple, somewhat familiar, and non-trivial

2) Being a continuous distribution, this pdf can be made to smoothly increase, starting from a minimum value at, $r=R$, to a maximum value at, $r=0$.

3) This pdf becomes particularly simple if the mode is chosen at, $r=0$. In this situation, $\mu =0$, and our two parameters reduce effectively to one, the shape parameter, $\sigma $.

4) This pdf can be made to pack an incredible mass within a relatively small volume, especially in 4-D space, and still keep the inside surface temperature very low.

5) The pdf has maximum entropy of all the continuous distributions. In fact, random error follows this kind of distribution. It seems, in this regard, to be an efficient way to pack radiation energy. We have a distribution of blackbody energy based on temperature called Planck’s law. Superimposed on this is a distribution based on radius, our truncated Gaussian.

For these reasons, we pursue this choice. It may turn out, however, that other truncated pdf’s have more desirable characteristics and behavior. This can be worked out in other papers. We emphasize that all properties within the black hole in this paper are pdf-model dependent. This includes the calculation of entropy, gravitational constant, gravitational force, etc. and so the choice is very important. Great flexibility is present, however, for working out various scenarios, with this, and other probability distribution functions.

We will assume a-priori, that the mode is to be set at, $r=0$. As mentioned several times, this fixes the parameter $\mu $. Under this assumption, $\mu =0$, and our pdf, ${f}_{r}$, takes the following, much-simplified form,

${f}_{r}=f\left(r,0,\sigma ;0,R\right)=\sqrt{\left(2/\text{\pi}\right)}1/\sigma \text{exp}\left(-{r}^{2}/\left(2{\sigma}^{2}\right)\right)1/\text{erf}\left(R/\left(\sqrt{2}\sigma \right)\right)$ (3-1)

Thus, the maximum radiative energy density is at, $r=0$, the very center of the black hole. In this equation, $\text{erf}\left(x\right)$ is the error function, and $\sigma $ is the shape parameter. The corresponding cumulative distribution function, ${F}_{r}^{CDF}$, is found using the defining relation, (2-12). The result of this integration is,

${F}_{r}^{CDF}=F\left(r,0,\sigma ;0,R\right)=\text{erf}\left(r/\sqrt{2}\sigma \right)/\text{erf}\left(R/\sqrt{2}\sigma \right)$ (3-2)

This leads trivially to, ${F}_{R}^{CDF}=1$. See Equations, (2-12), and, (2-8). In the special case where, $\mu =0$, the mean is given by the following, much-simplified, function,

$Mean=\sqrt{\left(2/\text{\pi}\right)}\text{\hspace{0.05em}}\sigma \left[1-\text{exp}\left(-{R}^{2}/\left(2{\sigma}^{2}\right)\right)\right]1/\text{erf}\left(R/\sqrt{2}\sigma \right)$ (3-3)

The variance can also be calculated. The result is

$\begin{array}{c}Var={\sigma}^{2}\left[1-R/\sigma \left[1+\text{erf}\left(R/\sqrt{2}\sigma \right)\right]/\text{erf}\left(R/\sqrt{2}\sigma \right)\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{2}{\text{\pi}}{\left[1-\text{exp}\left(-{R}^{2}/\left(2{\sigma}^{2}\right)\right)/\text{erf}\left(R/\sqrt{2}\sigma \right)\right]}^{2}\end{array}$ (3-4)

Other higher moments can be calculated. We note that for the truncated Gaussian, the mean is not simply equal to, $\mu $. Neither is the standard deviation, equal to, $\sigma $. The mean is calculated using Equation, (3-3). We see that its value depends on $\sigma $, the shape parameter, and radius, R. The standard deviation is given by the expression, $STD=\sqrt{Var}$, where the variance is written out in Equation, (3-4).

We know that Equation, (2-7), is a requirement. We can specialize this equation to zero radius, and, $r=R$. In the first instance, we obtain, with the help of Equation, (3-1),

${\rho}_{0}{A}_{0}={M}_{R}{f}_{0}={M}_{R}\sqrt{\left(2/\text{\pi}\right)}1/\sigma 1/\text{erf}\left(R/\sqrt{2}\sigma \right)$ (3-5)

In the second instance, we find, again with the help of Equation, (3-1),

${\rho}_{R}{A}_{R}={M}_{R}{f}_{R}={M}_{R}\sqrt{\left(2/\text{\pi}\right)}1/\sigma \text{exp}\left(-{R}^{2}/2{\sigma}^{2}\right)/\text{erf}\left(R/\sqrt{2}\sigma \right)$ (3-6)

The superscripts have been dropped as we are dealing here exclusively with 4-D quantities.

Equation, (3-5), is to be understood, in the following sense. In the limit where, $r\to 0$, this equality holds. A density function such as, ${\rho}_{0}$, requires a finite volume, and so, on the left hand side of Equation, (3-5), we require a minimum volume. The right hand side is definitely not volume dependent, and unequal to zero. For left and right hand sides to match, we imagine a limiting process where the right hand side is always finite, but the left hand side approaches, but never reaches, $r=0$. It is in this limit that we are to understand Equation, (3-5). We also notice that the right hand side of Equation, (3-6), is unequal to zero (unless $\sigma =\infty $ ), and thus finite in value.

We next divide the left hand side of Equation, (3-5), by the left hand side of Equation, (3-6). We do the same thing on the right hand side. Performing this operation lets us write

${\rho}_{0}{A}_{0}/\left({\rho}_{R}{A}_{R}\right)={f}_{0}/{f}_{R}=\text{exp}\left({R}^{2}/2{\sigma}^{2}\right)$ (3-7)

We shall see, shortly, that the quantity, ${R}^{2}/2{\sigma}^{2}$, is, in practice, already large when applied to black holes. Thus the ratio, ${f}_{0}/{f}_{R}$, is dramatically large.

More generally, following the steps above, if the radius lies within the range, $0\le r\le R$, we can use Equations, (2-7), and, (3-1), to demonstrate that

${\rho}_{r}{A}_{r}/\left({\rho}_{R}{A}_{R}\right)={f}_{r}/{f}_{R}=\text{exp}\left[\left({R}^{2}-{r}^{2}\right)/2{\sigma}^{2}\right]$ (3-8)

With this equation, we see that the product, ${\rho}_{r}{A}_{r}$, is monotonically decreasing with increasing radius. We have a maximum value for linear radiative mass density, ${\rho}_{r}{A}_{r}$, at, $r=0$, and this steadily decreases, until we achieve a minimum at $r=R$. Not only is radiative mass density, ${\rho}_{r}{A}_{r}$, decreasing with increasing r, but so too are, the internal energy density, ${u}_{r}{A}_{r}$, the radiative force, ${p}_{r}{A}_{r}={F}_{r}={F}_{r}^{\left(4\right)}$, and the radiative heat density, ${q}_{r}{A}_{r}$. Refer to Equations, (2-21), (2-22), (2-23), and, (2-24), where these specific proportionalities are established. These quantities, must therefore, also decrease according to Equation, (3-8), for increasing radius. It remains to determine the parameter, $\sigma $, and this we do next.

One way to find the shape parameter, $\sigma $, is to start with Equation, (2-25). We rewrite it here in the form,

${f}_{r}{M}_{R}=10{F}_{r}/{c}^{2}$ (3-9)

On the left hand side, we have the pdf, whereas on the right hand side we have the 4-D radiative force, at radius, r, within the black hole. We specialize Equation, (3-9), to radius, $r=R$. So we are just inside the event horizon, on the 4-D side. This gives

${f}_{R}{M}_{R}=10{F}_{R}/{c}^{2}$ (3-10)

We next recognize that for an isolated, static black hole (ISBH), with no heat inflow,

$5{F}_{R}^{\left(4\right)}=4{F}_{R}^{\left(3\right)}$ (ISBH) (3-11)

This was proven in reference [5] , using the generalized Stefan-Boltzmann law at the 3-D/4-D interface. Substituting this expression into the right hand side of Equation, (3-10), renders

${f}_{R}{M}_{R}=8{F}_{R}^{\left(3\right)}/{c}^{2}$ (ISBH) (3-12)

The quantity, ${F}_{R}^{\left(3\right)}$, is the radiative force right outside the event horizon, on the 3-D side. For an isolated, static black hole, it was assumed that the radiative temperature just outside the event horizon is equal to, 2.725 K, the CMB temperature. Thus, ${F}_{R}^{\left(3\right)}$, can easily be determined numerically. The result is

$\begin{array}{c}{F}_{R}^{\left(3\right)}={p}_{R}^{\left(3\right)}{A}_{R}^{\left(3\right)}=1/3{u}_{R}^{\left(3\right)}4\text{\pi}{R}^{2}\\ =1/3{a}^{\left(3\right)}{\left(2.725\right)}^{4}4\text{\pi}{R}^{2}\\ =1/3\left(7.5657\times {10}^{-16}\right){\left(2.725\right)}^{4}4\text{\pi}{R}^{2}\\ =1.747\times {10}^{-13}{R}^{2}\end{array}$ (ISBH) (3-13)

In the first line, we made use of the first equality in Equations, (2-16). For the 2^{nd} and 3^{rd} lines, it is well-known that, in 3-D space, the internal energy density can be expressed as,
${u}_{R}^{\left(3\right)}={a}^{\left(3\right)}{T}^{4}=\left(7.5657\times {10}^{-16}\right){T}^{4}$, in MKS units.

We next substitute the right hand side of Equation, (3-13), into the right hand side of Equation, (3-12), in order to eliminate 3-D radiative force, ${F}_{R}^{\left(3\right)}$. Bringing over the factor of total radiative mass, ${M}_{R}$, to the other side, we obtain

${f}_{R}=8\left(1.747\times {10}^{-13}\right){R}^{2}/\left({M}_{R}{c}^{2}\right)$ (ISBH) (3-14)

Next, because of the Schwarzschild relation, $R=2G{M}_{R}/{c}^{2}$, we can relate the radius of the black hole to its mass, and vice versa. This allows us to rewrite Equation, (3-14), entirely in terms of the mass of the black hole, or entirely, in terms of its radius. In terms of its mass, we find,

${f}_{R}=\left(3.416\times {10}^{-83}\right){M}_{R}$ (ISBH) (3-15)

In terms, of its radius, we obtain,

${f}_{R}=\left(2.303\times {10}^{-56}\right)R$ (ISBH) (3-16)

Either way, for a given size, or mass, we now have an expression for, ${f}_{R}$, the probability density function, just inside the 4-D event horizon.

Once we know, ${f}_{R}$, we can determine, $\sigma $, the shape parameter. Using Equation, (3-1), it is clear that

${f}_{R}=\sqrt{\left(2/\text{\pi}\right)}1/\sigma \text{exp}\left(-{R}^{2}/2{\sigma}^{2}\right)/\text{erf}\left(R/\sqrt{2}\sigma \right)$ (3-17)

Unfortunately, the right hand side of this equation is a complicated function of $\sigma $, as $\sigma $ occurs in three separate places. Thus, an analytical solution does not seem possible. We therefore resort to a numerical solution, using a trial and error approach. We step through various $\sigma $ values and look for a match between the right hand side of Equation, (3-17), and either one of Equations, (3-15), or (3-16). This is best achieved by using a spreadsheet, and we will give specific solutions for various black holes.

We will consider three black holes. The first will have the mass which equals that of the sun (not realistic), which we call black hole, A, abbreviated, $B{H}_{A}$. The second will have a mass ten times that of the sun. This we refer to as black hole, B, or $B{H}_{B}$. The third black hole is supermassive black hole with a mass, a million times that of the sun. These kinds of black holes are found at galactic centers. This is black hole, C, or $B{H}_{C}$, for short. Using these masses, and the Schwarzschild relation, we find

$\left({M}_{A},{M}_{B},{M}_{C}\right)=\left(1.989\times {10}^{30},1.989\times {10}^{31},1.989\times {10}^{36}\right)\text{kg}$

$\to \left({R}_{A},{R}_{B},{R}_{C}\right)=\left(2.954\times {10}^{3},2.954\times {10}^{4},2.954\times {10}^{9}\right)\text{meters}$ (3-18)

We next determine the probability distribution function at radius, R, which is, ${f}_{R}$, for all three black holes using either Equation, (3-15), or Equation, (3-16), whichever is convenient. We find

$\left({f}_{R,A},{f}_{R,B},{f}_{R,C}\right)=\left(6.8292\times {10}^{-53},6.8292\times {10}^{-52},6.8292\times {10}^{-47}\right){\text{meters}}^{-1}$ (3-19)

Using these values, we next determine the associated $\sigma $ values. We use Equation, (3-17), and our trial and error approach utilizing a spread sheet. We obtain a match for the following shape parameter values, which apply for each of the black holes considered.

$\left({\sigma}_{A},{\sigma}_{B},{\sigma}_{C}\right)=\left(1.951\times {10}^{2},1.992\times {10}^{3},2.241\times {10}^{8}\right)\text{meters}$ (3-20)

In this way, we can fix the parameters associated with our choice of probability distribution function. Fixing the $\mu $ values is easy, if we assume a maximum in ${f}_{r}$ at the origin. If the mode is at, $r=0$, then necessarily, we must set, $\left({\mu}_{A},{\mu}_{B},{\mu}_{C}\right)=\left(0,0,0\right)$. Finding the shape parameter, $\sigma $, however, is much more complicated. There does not appear to be an easier way of determining the $\sigma $ value, other than this trial and error approach, and looking for a match in solution.

Once more, we have focused on an isolated, static black hole (ISBH) because our goal in this paper is to illustrate a method, and not work out all possibilities. If we have net radiative inflow, then necessarily, $\text{d}Q/\text{d}t>0$. For this situation, we would have an evolving inside surface temperature at the event horizon, one that is continuously changing with increasing event horizon radius, $R=R\left(t\right)$. The inside surface temperature, ${T}_{1}$, decreases as R increases. The net radiative heat inflow can be determined from the outside blackbody temperature, and the size (or mass) of the black hole. This changing inside surface temperature, on the 4-D side, would translate into other values for, ${f}_{R}={f}_{R}\left(t\right)$, and, consequently, $\sigma =\sigma \left(t\right)$, which hypothetically can be worked out. We would obtain a different distribution of radiative mass density, or better said, a changing distribution, within the black hole because of Equation, (3-8), and the different $\sigma $, and, ${f}_{R}$ values. The point is that we would know how to proceed for this more involved situation. In addition, in the interests of keeping the discussion manageable, we have only considered three massive black holes. A whole catalog can be created for differing masses, and different heat inflows.

We close this section by giving the dimensionless ratio of

$\left({\sigma}_{A}/{R}_{A},{\sigma}_{B}/{R}_{B},{\sigma}_{C}/{R}_{C}\right)=\left(0.0660,0.0674,0.0759\right)$ (3-21)

Even though $\sigma $ is not the standard deviation for a truncated Normal distribution (see Equation, (3-4)), the Pearson ratio is a measure of dispersion. Due to the relatively small values indicated by Equations, (3-21), we can expect a high concentration of radiative mass, radiative pressure, radiative heat density, etc. located close to the center of the black hole. Towards the surface these values will drop down to almost insignificant amounts. This will be made apparent when we construct tables, in Sections 4 and 5, for the three black holes under consideration. There, we will see very clearly how the values are concentrated within a relatively narrow radius, close to the very core of the black hole. It is interesting to notice that this dispersion does not increase significantly, even for a very massive black hole.

4. Determining Temperature, Radiative Pressure, Internal Energy Densities, and Radiative Forces within the Black Hole

We have spent some time determining the shape parameter associated with our choice of probability distribution. This will enable us to determine the temperature as a function of radius, within the black hole, itself. Once we have the temperature distribution, we can also find other key thermodynamic quantities of interest within the black hole, and specify these as a function of 4-D radius. We are thinking of the radiative pressure, the internal energy density, the radiative heat density, and the total (radiative mass) energy density. These quantities depend on the temperature and temperature only. Moreover, if temperature is a function of radius, so too, are these quantities.

The radiative force within the black hole will also depend on temperature, but, in addition, on the surface area. The surface area, in turn, depends on the 4-D radius. By definition, according to Equation, (2-23), ${F}_{r}\equiv {p}_{r}{A}_{r}={p}_{r}2{\text{\pi}}^{2}{r}^{3}$, where we have dropped the superscript on, ${F}_{r},{p}_{r},{A}_{r}$, for ease of writing. Unless otherwise indicated, all quantities in this section, and the next, are 4-dimensional quantities. Therefore, we will dispense with the superscripts, unless we have to distinguish between 3-D and 4-D quantities. At radius, r, the force, ${F}_{r}$, is a radiative force pushing the surface area, ${A}_{r}=2{\text{\pi}}^{2}{r}^{3}$, out. At radius, $r+\text{d}r$, we have another radiative force, ${F}_{r+\text{d}r}$, which pushes the surface area, ${A}_{r+\text{d}r}$, at $r+\text{d}r$, from the outside in, because it acts from the top down. Together with gravity, which also pulls the layer between radii, r, and, $r+\text{d}r$, in, these forces guarantee for hydrostatic equilibrium, according to Equation, (2-26). Layer by layer, the black hole is thus kept in stable equilibrium. For a given radius, r, the forces, ${F}_{r}$, and ${F}_{r+\text{d}r}$, depend only on temperature, ${T}_{r}=T\left(r\right)$. And because, ${T}_{r}=T\left(r\right)>{T}_{r+\text{d}r}=T\left(r+\text{d}r\right)$, we have a greater radiative force pushing out at radius, r, than pushing in at radius, $r+\text{d}r$. In other words, ${F}_{r}>{F}_{r+\text{d}r}$. Refer to Equation, (2-26).

From Section 2, we know that Equation, (2-17), holds. With this equation we can calculate the internal energy density at any temperature, ${T}_{r}$. Once we have the internal energy density, it is easy to calculate the radiative pressure, the heat density, and the entropy density. We use Equations, (2-18). To find the total energy “density”, and the radiative mass density, we employ Equations, (2-20), and, (2-21). Thus, the temperature, ${T}_{r}$, is needed for all these quantities. But, ${T}_{r}$, can be determined. We start with Equation, (2-25). This equation can be rewritten in the forms,

${\rho}_{r}{A}_{r}={f}_{r}{M}_{R}$

$\left({e}_{r}/{c}^{2}\right){A}_{r}={f}_{r}{M}_{R}$

$\left(10/4\right)\left({u}_{r}/{c}^{2}\right){A}_{r}={f}_{r}{M}_{R}$

$\left(10/4\right)\left({a}^{\left(4\right)}{T}_{r}^{5}\right)/{c}^{2}{A}_{r}={f}_{r}{M}_{R}$

$\left(10/4\right)\left({a}^{\left(4\right)}{T}_{r}^{5}\right)/{c}^{2}2{\text{\pi}}^{2}{r}^{3}=\sqrt{\left(2/\text{\pi}\right)}1/\sigma \text{exp}\left(-{r}^{2}/\left(2{\sigma}^{2}\right)\right)1/\text{erf}\left(R/\left(\sqrt{2}\sigma \right)\right){M}_{R}$ (4-1)

For the 2^{nd} line, we substituted the first line in Equation, (2-21). In the 3^{rd} line, we employed the last equality in relations, (2-22). The 1^{st} line in Equation, (2-17), was invoked for the 4^{th} line in Equations, (4-1). And finally for the 5^{th} line, we substituted, Equation, (3-1), on the right hand side, and our, by now familiar, expression for the 4-dimensional surface area,
${A}_{r}=2{\text{\pi}}^{2}{r}^{3}$, on the left hand side. Upon inspection of the right hand side in the last of the series of Equations, in, (4-1), we see that this depends on the variables, r, R, and,
$\sigma $. The left hand side depends on the temperature,
${T}_{r}$, and the radius, r. Therefore, for a given r value, we can use this equation to determine the temperature as a function of radius. All that we need is the shape parameter and the mass, or size, of the black hole. Those values were determined in the previous section for the various specific black holes under consideration.

Equation, (4-1), however, is a rather complicated function. Another way to calculate the temperature, ${T}_{r}$, within the black hole is less complicated, and offers additional insight. We shall follow this approach, even though it is less direct, as it will allow us to derive the other quantities in a more systematic manner. We proceed in a series of steps, and use the three black holes introduced in the previous section, to illustrate the results. In this way numerical comparisons can be made, and highlighted. The three massive black holes under consideration, are those specified in Equations, (3-18), where we have specified both their masses, and their associated radii.

The first step is to determine the temperature just inside the event horizon, ${T}_{1}$. For this we use our generalized Stefan-Boltzmann law, derived in reference [5] . The equation reads,

$\left(1/c\right)\text{d}Q/\text{d}t=\left(4/3\right){a}^{\left(3\right)}{T}_{2}^{4}4\text{\pi}{R}^{2}-\left(5/4\right){a}^{\left(4\right)}{T}_{1}^{5}2{\text{\pi}}^{2}{R}^{3}$ (4-2)

In this equation, the constants, ${a}^{\left(3\right)}$, and ${a}^{\left(4\right)}$, have the numerical values in MKS units,

${a}^{\left(4\right)}\equiv \left(\text{3\pi}/2c\right){\sigma}^{\left(4\right)}=4.7481\times {10}^{-13}$ (4-3)

${a}^{\left(3\right)}\equiv \left(4/c\right){\sigma}^{\left(3\right)}=7.5657\times {10}^{-16}$ (4-4)

The quantity, ${\sigma}^{\left(3\right)}=5.67\times {10}^{-8}\text{Watts}/{\text{m}}^{\text{2}}\cdot {\text{K}}^{\text{4}}$, is the Stefan-Boltzmann constant in three dimensional space. In 4-D space, its extension [8] is given by, ${\sigma}^{\left(4\right)}=3.021\times {10}^{-5}\text{Watts}/{\text{m}}^{\text{3}}\cdot {\text{K}}^{\text{5}}$. In Equation, (4-2), the quantity, $\text{d}Q/\text{d}t$, is the net radiative heat entering the black hole from the 3-D side. We have shown in a previous work [5] that, $\text{d}Q/\text{d}t\ge 0$, is a direct consequence of our 4-D black hole radiation model, as we are ignoring/discounting radiative evaporation, for example, due to Hawking radiation. These are second order effects, at best, if they exist, as shown in reference [5] . The temperatures, ${T}_{2}$, and ${T}_{1}$, in Equation, (4-2), are the temperatures just outside, and just inside, the event horizon, respectively. The temperature, ${T}_{2}$, is on the 3-dimensional side, whereas, the temperature, ${T}_{1}$, is on the 4-D side within the black hole. The, R, is, of course, the radius of the black hole event horizon, proportional to its mass, because of the Schwarzschild relation.

We notice that Equation, (4-2), contains two terms on the right hand side. The first term on the right hand side, is the radiative heat leaving 3-D space, and entering the 4-D black hole. We assume that the event horizon is infinitely thin. The second term on the right hand side is the radiative heat exiting the 4-D black hole, and entering 3-D space. Both propagate at velocity, c. Another way of writing Equation, (4-2), is thus

$1/c\text{d}Q/\text{d}t=1/c\text{d}{Q}^{\left(3\right)}/\text{d}t-1/c\text{d}{Q}^{\left(4\right)}/\text{d}t$ (4-5)

In this expression, $\text{d}{Q}^{\left(3\right)}/\text{d}t$ is the radiative heat streaming into the black hole, and, $\text{d}{Q}^{\left(4\right)}/\text{d}t$, is the radiative heat flowing out. This expression leads to a necessarily positive radiative surface tension, and was a primary motivation for the construction of a black hole as a 4-D object.

There are other ways to write Equation, (4-2), the generalized Stefan-Boltzmann law. A particularly straight-forward way is in terms of radiative pressure, and radiative force. In this formulation,

$\left(1/c\right)\text{d}Q/\text{d}t=4{p}_{2}^{\left(3\right)}{A}_{R}^{\left(3\right)}-5{p}_{1}^{\left(4\right)}{A}_{R}^{\left(4\right)}$ (4-6)

here, the surface areas, ${A}_{R}^{\left(3\right)}=4\text{\pi}{R}^{2}$, and, ${A}_{R}^{\left(4\right)}=2{\text{\pi}}^{2}{R}^{3}$, are the 3-dimensional and 4-dimensional surface areas, respectively, at the event horizon. Also, it is to be noticed that the 3-D radiative pressure, ${p}_{2}^{\left(3\right)}=1/3{u}_{2}^{\left(3\right)}=1/3{a}^{\left(3\right)}{T}_{2}^{4}$. Thus, ${p}_{2}^{\left(3\right)}$, is defined in terms of the temperature right outside the event horizon. By contrast, the 4-D radiative pressure is defined in terms of a different temperature, ${T}_{1}$. There, we have, ${p}_{1}^{\left(4\right)}=1/4{u}_{1}^{\left(4\right)}=1/4{a}^{\left(4\right)}{T}_{1}^{5}$, where, ${T}_{1}$, is the temperature right inside the event horizon, within the 4-D black hole.

Equation, (4-6), for net radiative heat inflow, can also be rewritten as,

$\left(1/c\right)\text{d}Q/\text{d}t=4{F}_{R}^{\left(3\right)}-5{F}_{R}^{\left(4\right)}$ (4-7)

where,

${F}_{R}^{\left(3\right)}\equiv {p}_{2}^{\left(3\right)}{A}_{R}^{\left(3\right)}$ and ${F}_{R}^{\left(4\right)}\equiv {p}_{1}^{\left(4\right)}{A}_{R}^{\left(4\right)}$ (4-8a,b)

The particular factors of, 4, and, 5, in Equations, (4-6), and, (4-7), are imposed upon us by the dimensionality of the differing spaces. The radiative force, ${F}_{R}^{\left(3\right)}$, acts on the event horizon and pulls the event horizon membrane radially in. The radiative force, ${F}_{R}^{\left(4\right)}$, by contrast, acts on the infinitely thin event horizon, and pushes the membrane radially out. Because $\text{d}Q/\text{d}t\ge 0$, we see that, $4{F}_{R}^{\left(3\right)}\ge 5{F}_{R}^{\left(4\right)}$. The radiative surface tension, defined as, ${F}_{\text{ST}\text{\hspace{0.17em}}\text{radiative}}\equiv {F}_{R}^{\left(3\right)}-{F}_{R}^{\left(4\right)}$, was thus shown to be always positive.

We will next assume that there is no net inflow. In other words, $\text{d}Q/\text{d}t=0$, and we have what we refer to as an isolated, static black hole (ISBH). Net inflow will be treated in another paper and it will be easy to extend these ideas once we understand the process for an isolated static black hole. For an ISBH, we assume that the outside temperature is, ${T}_{2}=2.725\text{\hspace{0.17em}}\text{K}$, a working hypothesis. This cannot be proven, just surmised, because, at this small temperature, net inflows will not lead to appreciable expansion. In other words, expansion, if it existed, cannot be observed as it would be so slight. As such, at ${T}_{2}=2.725\text{\hspace{0.17em}}\text{K}$, we may, or, may not, have inflow. We assume no inflow. Under this assumption, the CMB photons surrounding the black hole will lead to a temperature, called the equilibrium temperature, just inside the event horizon equal to,

${T}_{1}^{5}R=0.0597$ or ${T}_{1}=0.569{R}^{-1/5}$ (ISBH) (4-9a,b)

This is obtained by solving Equation, (4-2), after we set the left hand side equal to zero. Remember that, for, ${T}_{2}=2.725\text{\hspace{0.17em}}\text{K}$, we assume we have an isolated, static black hole.

Even though Equations, (4-9a,b), were derived for an isolated, static black hole, they tell us how to proceed for a dynamic black hole (DBH), with $\text{d}Q/\text{d}t>0$. For a different outside surface temperature, ${T}_{2}>2.725\text{\hspace{0.17em}}\text{K}$, we will have net inflow, and, $R=R\left(t\right)$. The relation to be used here for ${T}_{1}$ is the same as before, except that now, R is evolving, and getting larger. As a consequence, ${T}_{1}$, will start to decrease by the equations above. The radius will start to increase with net heat inflow, but until it has achieved a new value, we must assume that the inside temperature hasn’t caught up yet. In other words it is still the temperature as determined by Equation (4-9b). We can argue as follows. The inside surface temperature determines, through the shape parameter $\sigma $, the inner mass distribution within the black hole. This, in turn, determines its total mass. The total mass determines the total radius. And the total radius relates back to the inside temperature by use of Equations, (4-9a,b). Until such time that the radius increases, ${T}_{1}$, will be fixed by the above equation. In a sense, a dynamic black hole is an evolving set of isolated, static black holes, each with a different radius. The energy absorbed will determine the final radius, irrespective of how much time it takes. The rate of inflow will determine that.

If we go back in cosmological time, the CMB temperature was higher. Thus, Equations, (4-9a,b), would have to be modified as we assumed that the CMB temperature was 2.725 K to obtain the factors of 0.0597 and 0.569 in these expressions. For higher (or lower) CMB temperatures, we have to scale both factors as follows:

$\left(0.0597\right)\to {\left({T}_{CMB}/2.725\right)}^{4}\left(0.0597\right)$, and,

$\left(0.569\right)\to {\left({T}_{CMB}/2.725\right)}^{4/5}\left(0.569\right)$,

to take into account different CMB temperatures. If we go back in cosmological time, for the same radius, the equilibrium inside surface temperature was higher. That means for the same outside temperature load, ${T}_{2}>{T}_{CMB}$, inflow must have been less. In future epoch, it will be greater under the same conditions. As, ${T}_{CMB}\to 0$, we obtain (reproduce) the permanent inflow model discussed in the introduction.

From Equations, (4-9a), we see that the mass (size) of a black hole determines the inside surface temperature. The greater the radius, the less of an inside surface temperature, we will have. For the three black holes considered in relations, (3-18), upon using Equations, (4-9a), we find that

$\left({T}_{1,A},{T}_{1,B},{T}_{1,C}\right)=\left(0.115\text{\hspace{0.17em}}\text{K},0.0726\text{\hspace{0.17em}}\text{K},0.00726\text{\hspace{0.17em}}\text{K}\right)$ (4-10)

These temperatures are many orders of magnitude greater than the corresponding Hawking temperatures, and yet, they are below the minimum outside surface temperature of, 2.725 K. For example, the predicted Hawking temperature for a black hole having a mass equal to that of the sun is only, 62 nK, which is much less than 0.115 K, shown in Equation (4-10). But then, the space considered was 3-dimensional, and a totally different specific process was considered. We also remark that the variation in temperature between the various massive black holes is not as dramatic as might be expected. This is due to the ${R}^{-1/5}$ factor in Equation, (4-9b), which in effect, is equivalent to ${M}_{R}^{-1/5}$ because of the Schwarzschild proportionality between mass and radius. The Hawking temperature, by contrast, equals, ${T}_{H}=\hslash {c}^{3}/\left(8\text{\pi}GM{k}_{B}\right)$, and therefore is proportional to ${M}_{R}^{-1}$. His temperature, therefore, falls off much more dramatically than ours, as one increases black hole mass. This first step of finding the inside 4-D surface temperature, really consists of applying Equation, (4-9a), or, Equation, (4-9b), to find the inside surface temperature.

The second step consists of strict calculations. In this step, we determine the internal energy density, the radiative pressure, the radiative heat density, and the total entropy density, defined by Equations, (2-17), and, (2-18). These quantities are all to be calculated just inside the event horizon, where we know the temperature, ${T}_{1}$. For the internal energy density, ${u}_{R}={u}_{1}={u}_{1}^{\left(4\right)}$, we use Equation, (2-17). The respective temperatures are listed in Equations, (4-10). The results for the three black holes considered are thus,

$\left({u}_{1,A},{u}_{1,B},{u}_{1,C}\right)=\left(9.60\times {10}^{-18},9.60\times {10}^{-19},9.60\times {10}^{-24}\right)\text{J}/{\text{m}}^{\text{4}}$ (4-11)

These values are very small, but then the temperature, just within the event horizon is very small. We can now use Equations, (2-18), to find the radiative pressures, the heat densities, the entropy densities. The results, with appropriate units, are:

$\left({p}_{1,A},{p}_{1,B},{p}_{1,C}\right)=\left(2.40\times {10}^{-18},2.40\times {10}^{-19},2.40\times {10}^{-24}\right)\text{N}/{\text{m}}^{\text{3}}$ (4-12)

$\left({q}_{1,A},{q}_{1,B},{q}_{1,C}\right)=\left(1.20\times {10}^{-17},1.20\times {10}^{-18},1.20\times {10}^{-23}\right)\text{J}/{\text{m}}^{\text{4}}$ (4-13)

$\left({s}_{1,A},{s}_{1,B},{s}_{1,C}\right)=\left(1.04\times {10}^{-16},1.65\times {10}^{-17},1.65\times {10}^{-21}\right)\text{J}/\left({\text{m}}^{\text{4}}\cdot \text{K}\right)$ (4-14)

Equations, (4-11), (4-12), (4-13), and (4-14), hold for an isolated, static black hole (ISBH), and until the radius has changed, for a dynamic black hole (DBH), as well.

In step three, we calculate the total (radiative mass) energy “density”, the 4-D mass density, and the radiative force, all at, $r=R$, just inside the 4-D event horizon. For the total energy density, we use Equation, (2-20), in conjunction with Equations, (4-11), (4-12), and, (4-13). The results are:

$\left({e}_{1,A},{e}_{1,B},{e}_{1,C}\right)=\left(2.40\times {10}^{-17},2.40\times {10}^{-18},2.40\times {10}^{-23}\right)$ (4-15)

These values, numerically, are twice the radiative heat density, which confirms the second equality in the Equations, (2-22). In fact, all the equalities in Equation, (2-22), are found to hold true. For the mass density, we can use the first line in Equation, (2-21), together with Equation, (4-15). Here, we obtain

$\left({\rho}_{1,A},{\rho}_{1,B},{\rho}_{1,C}\right)=\left(2.67\times {10}^{-34},2.67\times {10}^{-35},2.67\times {10}^{-40}\right)\text{kg}/{\text{m}}^{\text{4}}$ (4-16)

These values are exceedingly low, but then we are in 4-D space. And, for the radiative force, Equation, (2-23), can be utilized where we use the radiative pressure values indicated in Equations, (4-12), together with the radii specified in Equations, (2-18). The radii are needed to calculate the 4-D surface areas using the formula, ${A}_{R}=2{\text{\pi}}^{2}{R}^{3}$. The net results of these calculations give,

$\left({F}_{1,A},{F}_{1,B},{F}_{1,C}\right)=\left(1.22\times {10}^{-6},1.22\times {10}^{-4},1.22\times {10}^{6}\right)\text{Newtons}$ (4-17)

These 4-D radiative forces, pushing out against the event horizon, seem incredibly small, considering the fact that we are dealing with very massive objects. They, in fact, are. We should keep in mind, however, that they push against a greater 4-D surface area. Radiative forces, as it will turn out, are insignificant when compared to gravitational forces, except in the most extreme of circumstances, close to the event horizon of a black hole. This will be seen shortly. Within the interior however, radiative forces will have to be comparable to gravitational forces for hydrostatic equilibrium. See Equation, (2-26). On the surface we do not have hydrostatic equilibrium; we have, instead, an intrinsic surface tension which is non-zero due to the changing dimensionality of space. Equations, (4-15), (4-16), and (4-17), hold only at the event horizon, just inside the 4-D black hole, because it is there that we have established a temperature. As we shall see shortly, all the values indicated by Equations, (4-11), through to, (4-17), will increase dramatically as we enter the interior of the black hole. The subscript, “1” in all these relations can be interchanged with, “R”, because the temperature, ${T}_{1}$, holds at radius R, inside the black hole. As examples, we can write ${p}_{1,A}={p}_{R,A},{p}_{1,B}={p}_{R,B},{q}_{1,C}={q}_{R,C}$, etc.

Thus far, all the pressures, densities, and forces are given at radius, R, because at that radius, we know the temperature. We now look into the interior, where $r<R$. In step 4, we will calculate the mass density as a function of radius, ${\rho}_{r}=\rho \left(r\right)$. For that, we will employ Equation, (3-8). This can be written in the form,

$\left({\rho}_{r}{r}^{3}\right)/\left({\rho}_{R}{R}^{3}\right)={f}_{r}/{f}_{R}=\text{exp}\left[\left({R}^{2}-{r}^{2}\right)/2{\sigma}^{2}\right]$ (4-18)

We notice that in the limit where, $r\to 0$, the right hand side is non-zero. Therefore the left hand side must also be finite in this limit. This can only mean one thing for the left hand side, namely, that ${\rho}_{r}$ increases more rapidly, than ${r}^{3}$ decreases, in this limit. Also remember that ${\rho}_{r}$ is really a density. A density, a-priori, requires a finite volume, however small. The pdf chosen is seen to give us a finite result for the right hand side given this limit. And so, without further analysis, we must accept that the left hand side is also finite in the limit where, $r\to 0$.

For each of the black holes considered, we know the shape parameter,
$\sigma $, and the total radii, R. Therefore, for a specified
$r<R$ value, we can evaluate the right hand side of Equation, (4-18). We also know the
${\rho}_{R}$ values on the left hand side, which hold at the surface. These are given by Equations, (4-16). We can therefore determine the inner radiative mass density,
${\rho}_{r}$, as a function of radius. These results have to be given in table form, because we are considering separate, r/R values. The tables are presented in Tables 1-3. Table 1, applies to a black hole having a mass equal to that of the sun. Table 2, holds for a black hole, having a mass ten times that of the sun. And Table 3 refers to the 3^{rd} black hole under consideration, one having a mass a million times that of the sun. All tables will increase the radius in steps indicated. And so, in the first column, we have
$r=0.03R,0.06R,0.09R,\cdots $. Notice that we do not start with,
$r=0$, as it is meaningless to define a density, or for that matter, a temperature, for zero volume. We therefore start with an almost infinitesimally small, but finite radius,
$r={10}^{-6}R$. The associated,
${\rho}_{r}$, values are calculated under column 3 in these tables. Column 2 is reserved for,
${f}_{r}/{f}_{R}=\text{exp}\left[\left({R}^{2}-{r}^{2}\right)/2{\sigma}^{2}\right]$, which is specific to our choice of pdf.

The fifth step is about calculating the total energy “density”, using the first line in relation, (2-21). And then, by using Equations, (2-22), we evaluate the radiative pressure, the radiative heat density, and the internal energy density. Since the 4-D mass density, ${\rho}_{r}$, within the black hole is a function of radius, r, so too are, ${e}_{r},{p}_{r},{q}_{r}$ and, ${u}_{r}$. The total energy density is given by the simple relation, ${e}_{r}={\rho}_{r}{c}^{2}$, where c is the speed of light. Again, because we are dealing with functions that depend on, r/R, we have to give the results in table form. These are presented in Tables 1-3, where each table refers to a black hole of a specific mass. The values for, ${u}_{r},{p}_{r}$ and, ${q}_{r}$, are highlighted under columns, 4, 5, and 6, respectively. The values for total (radiative mass) energy densities, ${e}_{r}$, are not

Table 1. For a blackhole having 1 solar mass.

Table 2. For a Blackhole having 10 solar masses.

Table 3. For a Blackhole having 10^{6} solar masses.

tabulated, as this is simply the sum of columns, 4, 5 and 6. Or, what is equivalent, we could take the entries under column, 3, for, ${\rho}_{r}$, in the tables, and, multiply those by ${c}^{2}$ to obtain ${e}_{r}$.

In step 6, we finally evaluate the temperature, ${T}_{r}$. This step is relatively easy. Instead of using Equation, (4-1), we will use a much simpler expression, Equation (2-17). Solving Equation, (2-17), for temperature, ${T}_{r}$, gives

${T}_{r}={\left({u}_{r}/{a}^{\left(4\right)}\right)}^{1/5}$ (4-19)

In this equation, ${a}^{\left(4\right)}=4.7481\times {10}^{-13}$, in MKS units, as specified in Equation, (2-17). Using the ${u}_{r}$ values specified in the tables, given under column, 4, we can thus evaluate the corresponding temperatures by use of Equation, (4-19). These temperatures are listed under column, 7, in each of our three tables. When comparing the temperatures at the very core for the three black holes considered, we see something very interesting. The temperature actually goes down at the core as the black hole increases its mass. But then the surface temperature just inside the event horizon, ${T}_{1}$, also decreases its temperature as radius, R, increases in value. This would mean that the smallest black holes have the highest temperatures at their core. This fact could have important ramifications for high energy physics, and perhaps even modeling massive particles as mini black holes.

It also shows that, for a dynamic black hole, the inside surface temperatures will decrease as the radius expands. Assume that a black hole expands from a mass of one solar radius to a mass of ten times that value. Depending on the rate of inflow, this will occur within a certain period of time. Due to the Schwarzschild condition, the radius also increases ten-fold. We are in effect going from Table 1 and Table 2. Now look under the temperature column. We see that not only has the surface temperature just inside the event horizon decreased, but also the core temperature. Think of a dynamic black hole as a series of snapshots of isolated, static black holes, each with an increased radius.

Retracing our steps, we should recognize that, ${u}_{r},{q}_{r},{p}_{r},{e}_{r}$, are all proportional to one another, and thus to ${\rho}_{r}$. Using the proportionalities, we can claim that,

${u}_{r}/{u}_{R}={q}_{r}/{q}_{R}={p}_{r}/{p}_{R}={e}_{r}/{e}_{R}={\rho}_{r}/{\rho}_{R}={\left({T}_{r}/{T}_{R}\right)}^{5}$ (4-20)

By setting up a ratio, the constants relating these individual quantities factor out. For the final equality, in Equations, (4-20), we employed the relation, (2-17). Moreover, the mass density, ${\rho}_{r}$, in turn, is related to ${\rho}_{R}$, through Equation, (3-18). Because of Equation, (3-18), we can further write,

$\begin{array}{c}{u}_{r}{A}_{r}/\left({u}_{R}{A}_{R}\right)={q}_{r}{A}_{r}/\left({q}_{R}{A}_{R}\right)={p}_{r}{A}_{r}/\left({p}_{R}{A}_{R}\right)={e}_{r}{A}_{r}/\left({e}_{R}{A}_{R}\right)\\ ={\rho}_{r}{A}_{r}/\left({\rho}_{R}{A}_{R}\right)={f}_{r}/{f}_{R}=\text{exp}\left[\left({R}^{2}-{r}^{2}\right)/2{\sigma}^{2}\right]\end{array}$ (4-21)

An equivalent way of writing this equation is in terms of radii. Because, ${A}_{r}=2{\text{\pi}}^{2}{r}^{3}$, and ${A}_{R}=2{\text{\pi}}^{2}{R}^{3}$, it is possible to write,

$\begin{array}{c}{u}_{r}{r}^{3}/\left({u}_{R}{R}^{3}\right)={q}_{r}{r}^{3}/\left({q}_{R}{R}^{3}\right)={p}_{r}{r}^{3}/\left({p}_{R}{R}^{3}\right)={e}_{r}{r}^{3}/\left({e}_{R}{R}^{3}\right)\\ ={\rho}_{r}{r}^{3}/\left({\rho}_{R}{R}^{3}\right)={f}_{r}/{f}_{R}=\text{exp}\left[\left({R}^{2}-{r}^{2}\right)/2{\sigma}^{2}\right]\end{array}$ (4-22)

Moreover, we recognize that, by Equation, (2-23), the radiative force at radius, r, is given by,
${F}_{r}={p}_{r}{A}_{r}$. From this, it follows that,
${F}_{R}={p}_{R}{A}_{R}$. The 3^{rd} equality in Equation, (4-21), can therefore be re-expressed as,

${F}_{r}/{F}_{R}={f}_{r}/{f}_{R}=\text{exp}\left[\left({R}^{2}-{r}^{2}\right)/2{\sigma}^{2}\right]$ (4-23)

We see from Equations, (4-21), (4-22), and (4-23), that all quantities on the left hand side, are, in effect, determined by the ratio,
${f}_{r}/{f}_{R}$. The numerical values for
${f}_{r}/{f}_{R}$ are specified in the 2^{nd} column in each of the tables. Again, this specific ratio is fixed by our choice of pdf.

The temperature is also related to this ratio, ${f}_{r}/{f}_{R}$. The temperature profile within the black hole, is found using the internal energy density, by Equations, (2-17), or, what is equivalent, (4-19). Using the version indicated by Equation, (4-19), and combining this with the first equality in Equations, (4-20), we are justified in writing,

${T}_{r}/{T}_{R}={\left({u}_{r}/{u}_{R}\right)}^{1/5}={\left({\rho}_{r}/{\rho}_{R}\right)}^{1/5}$ (4-24)

Moreover, multiplying both left and right hand sides by, ${\left(r/R\right)}^{3/5}$, we find that,

$\left({T}_{r}/{T}_{R}\right){\left(r/R\right)}^{3/5}={\left[{\rho}_{r}{r}^{3}/\left({\rho}_{R}{R}^{3}\right)\right]}^{1/5}={\left({f}_{r}/{f}_{R}\right)}^{1/5}=\text{exp}\left[\left({R}^{2}-{r}^{2}\right)/10{\sigma}^{2}\right]$ (4-25)

The last two equalities in Equation, (2-22), were utilized in obtaining the results in Equations, (4-25). We notice from relation, (4-25), that

$\left({T}_{r}/{T}_{R}\right)={\left(R/r\right)}^{3/5}{\left({f}_{r}/{f}_{R}\right)}^{1/5}$ (4-26)

Equations, (4-24), (4-25), and, (4-26), are some other ways to calculate the temperature, ${T}_{r}$, within the black hole. However, Equation, (4-19), seems to us to be the most direct way to evaluate these temperatures.

The seventh step is to evaluate the 4-D entropy density function, ${s}_{r}=s\left(r\right)$. This is considered a state variable, i.e., it will depend on the state of the system, defined as a particular layer of the black hole, at a specific temperature. The easiest way to find, ${s}_{r}$, is to use the last equation, in relations, (2-18). We know the heat density, ${q}_{r}$, as well as temperature within a specific layer, ${T}_{r}$. These values were tabulated under columns, 6, and 7, respectively. Thus we can find, ${s}_{r}$, by simply dividing the entry in one column, column 6, by the corresponding entry in the other, column 7. The results for, ${s}_{r}$, are listed under column 8, in Tables 1-3.

An interesting formulation for entropy density is to use a ratio. We can compare the entropy density within the interior of the black hole, ${s}_{r}$, to its value just inside the event horizon, where ${s}_{r}={s}_{R}$. In step 2, we evaluated, $\left({s}_{R,A},{s}_{R,B},{s}_{R,C}\right)=\left({s}_{1,A},{s}_{1,B},{s}_{1,C}\right)$. See Equation, (4-14). We find, upon constructing the ratio,

$\begin{array}{c}{s}_{r}/{s}_{R}=\left({q}_{r}/{q}_{R}\right)\left({T}_{R}/{T}_{r}\right)={\left(R/r\right)}^{3}\left({f}_{r}/{f}_{R}\right)\left({T}_{R}/{T}_{r}\right)\\ ={\left(R/r\right)}^{3}\left({f}_{r}/{f}_{R}\right){\left(R/r\right)}^{-3/5}{\left({f}_{r}/{f}_{R}\right)}^{-1/5}\\ ={\left(R/r\right)}^{12/5}{\left({f}_{r}/{f}_{R}\right)}^{4/5}\\ ={\left(R/r\right)}^{12/5}\text{exp}\left[0.4\left({R}^{2}-{r}^{2}\right)/{\sigma}^{2}\right]\end{array}$ (4-27)

For the 1^{st} line, we used the last relation in Equations, (2-18), and the 2^{nd} equality in Equations, (4-22). For the second line, Equation, (4-26), was employed. The 3^{rd} line involves simplification. And, as for the 4^{th} line, the last equality, in Equation, (4-21), was utilized. From the last line in relations, (4-27), we see that the entropy density increases rather dramatically if one decreases the radius, r. As stated, the values for,
${s}_{r}$, are specified in the tables. They are to be found under column, 8.

In step 9, we calculate the 4-D surface area the 4-D radiative force, ${F}_{r}={p}_{r}{A}_{r}$. These values depend on the radius, and also, on the temperature at radius, r. The evaluation of surface area is easy; we simply use ${A}_{r}=2{\text{\pi}}^{2}{r}^{3}=2{\text{\pi}}^{2}{R}^{3}{\left(r/R\right)}^{3}$. These values are not tabulated. For, ${F}_{r}={p}_{r}{A}_{r}$, we multiply the radiative pressure, entered under column 5, by the corresponding surface area. Those results are tabulated under column, 9, in the three tables.

As was done for entropy density, we can construct a ratio, which will tell us how the radiative force increases within the interior of the black hole. But this was already done by way of Equation, (4-23). Again, the factor,
${f}_{r}/{f}_{R}$, comes into play, and this factor is listed under the 2^{nd} column in all three tables. Therefore, because,
${F}_{r}/{F}_{R}={f}_{r}/{f}_{R}$, we need only look under the 2^{nd} column in the three tables, to see how radiative force increases as one penetrates into the interior. Needless to say, the radiative force increases dramatically with decreasing r value. Although ridiculously weak at the surface of the black hole, its value within the interior is strong. We can think of the black hole as a multi-layered, 4-D capacitor, where energy densities, radiative forces, and radiative pressures increase substantially with decreasing radii.

In step 10, we determine the temperature gradient, $\text{d}{T}_{r}/\text{d}r$, within the 4-D black hole. For this, we will utilize Equation, (4-25), the first and the last entry. Written out somewhat differently, we have

${T}_{r}^{5}{r}^{3}={T}_{R}^{5}{R}^{3}\text{exp}\left[\left({R}^{2}-{r}^{2}\right)/10{\sigma}^{2}\right]$ (4-28)

We differentiate both the left, and the right, hand sides with respect to radius, r. The result is

$5{T}_{r}^{4}\text{d}{T}_{r}/\text{d}r{r}^{3}+{T}_{r}^{5}3{r}^{2}={T}_{R}^{5}{R}^{3}\text{exp}\left[\left({R}^{2}-{r}^{2}\right)/10{\sigma}^{2}\right]\left(-r/{\sigma}^{2}\right)$ (4-29)

We next divide the left hand side of Equation, (4-29), by the left hand side of Equation, (4-28). We do the same with the right hand side; divide the right hand side of Equation, (4-29), by the right hand side of Equation, (4-28). In this way we obtain,

$\left(5/{T}_{r}\right)\text{d}{T}_{r}/\text{d}r+3/r=\left(-r/{\sigma}^{2}\right)$ (4-30)

Hence, it follows that,

$\text{d}{T}_{r}/\text{d}r=-\left(r/{\sigma}^{2}+3/r\right){T}_{r}/5$ (4-31)

The negative sign tells us that we will have a decrease in temperature, in the amount of, $\text{d}{T}_{r}$, for a corresponding increase in radius, from r to $r+\text{d}r$. This decrease depends, of course, on the shape parameter, $\sigma $, as well as on the radius, r, and temperature, ${T}_{r}$.

We know the values on the right hand side of Equation, (4-31). The radius is specified under column, 1, in all tables. We use, $r=\left(r/R\right)R$. For the temperatures, ${T}_{r}$, we utilize the corresponding entries under column 7 in all three tables. And the shape parameters are listed in Equations, (3-20). Thus, the right hand side of Equation, (4-31), can be evaluated. We list, $\text{d}{T}_{r}/\text{d}r$, under column, 10, in Tables 1-3.

For step, 11, we calculate the 4-D blackbody specific heat, ${C}_{r}\equiv \text{d}{u}_{r}/\text{d}{T}_{r}$, at radius, r. Using Equation, (2-17), we see that,

${C}_{r}\equiv \text{d}{u}_{r}/\text{d}{T}_{r}=5{a}^{\left(4\right)}{T}_{r}^{4}$ (4-32)

However, we also know from the last two equalities in relations, (2-18), that it is permissible to write,

${s}_{r}=5/4{u}_{r}/{T}_{r}=5/4{a}^{\left(4\right)}{T}_{r}^{4}$ (4-33)

We utilized Equation, (2-17), for the last equality. Upon comparing Equations, (4-32), with Equation, (4-33), it is obvious that,

${C}_{r}=4{s}_{r}$ (4-34)

This very simple result makes it easy to calculate the specific heat, or, $\text{d}{u}_{r}/\text{d}{T}_{r}$. We simply take the entries for, ${s}_{r}$, specified under column 8 in the tables, and multiply them by a factor of 4. Those values are listed under column 11 in all tables.

Our final step, in this section, is step, 12. Here, we calculate the internal energy density gradient, $\text{d}{u}_{r}/\text{d}r$. We’ll use the fact that,

$\text{d}{u}_{r}/\text{d}r=\left(\text{d}{u}_{r}/\text{d}{T}_{r}\right)\left(\text{d}{T}_{r}/\text{d}r\right)={C}_{r}\left(\text{d}{T}_{r}/\text{d}r\right)$ (4-35)

The entries for, ${C}_{r}$, are given under column 11 in Tables 1-3. We also know the corresponding values for, $\text{d}{T}_{r}/\text{d}r$, as these are listed under column 10. Thus, we are able to evaluate the right hand side of Equation, (4-35). We display those values, in our final column, column 12, in all tables. We remark that because, $\text{d}{T}_{r}/\text{d}r$, is negative, $\text{d}{u}_{r}/\text{d}r$, will also be negative. Therefore, for an increase in 4-D radius, from radius, r, to radius, $r+\text{d}r$, we expect a decrease in internal energy density, $\text{d}{u}_{r}$, at a specific radius, r. Not only that, but the radiative pressure, the heat density, and the radiative mass density, will also decrease. Refer to Equations, (4-20), where this can be seen directly.

5. Determining Radiative Mass, Work Done against Radiative Pressure, Gravitational Forces, and Entropy within the Black Hole

In this section, we will use the two stellar-like equations, which have been generalized to four spatial dimensions, in Section 2. We wish to calculate some important global attributes which hold within, and on the surface, of the black hole. By doing so, we will appreciate the role that the chosen pdf plays in the determination of these physical characteristics. We are looking specifically, at the radiative mass, the work needed to assemble the black hole against radiative pressure, the gravitational forces both within, and on the surface of the black hole, and the entropy. To keep the discussion simple, we will consider only an isolated, static black hole, where there is no net radiative inflow. Net inflow can be considered in a follow up paper. We will also assume that the radiative temperature just outside the event horizon, on the 3-D side, is the CMB temperature, 2.725 K.

As in the previous section, we will proceed in steps, as it will help us to organize our thoughts in a coherent way. Also, we will dispense with the superscripts as all quantities, because, unless otherwise indicated, the quantities are 4-D. The three black holes under consideration are the same as before. Black hole, $B{H}_{A}$, will have a mass equal to that of the sun. Black holes, $B{H}_{B}$ and, $B{H}_{C}$, have masses, ten and a million times that of the sun, respectively. Finally, because these global quantities or attributes are dependent on radius, we will have to present our results in table form. We will have three tables; Table 4 refers to black hole, $B{H}_{A}$, Table 5 to black hole, $B{H}_{B}$, and Table 6 relates to black hole, $B{H}_{C}$.

We start with the cumulative distribution function (cdf), given by Equation, (3-2). The defining relation, for, ${F}_{r}^{CDF}$, is Equation, (2-12). Written out once more, we have

${F}_{r}^{CDF}\equiv {\displaystyle {\int}_{0}^{r}{f}_{r}\text{d}r}=\text{erf}\left(r/\sqrt{2}\sigma \right)/\text{erf}\left(R/\sqrt{2}\sigma \right)$ (5-1)

For each of the black holes under consideration, black holes $B{H}_{A},B{H}_{B}$, and, $B{H}_{C}$, we have a specific R, and, a specific $\sigma $ value. These are indicated in Equations, (3-18) and (3-20). Therefore, it is possible to evaluate, numerically, the right hand side of Equation, (5-1), for various $r=\left(r/R\right)R$ values. These ${F}_{r}^{CDF}$ values are listed under column 2, in Tables 4-6. Under column 1, we list the internal radii, in steps of, $r=0,0.03R,0.06R,0.09R$, etc.. This is step 1. All tables are presented at the end of the paper.

Step 2 will determine the radiative mass, ${M}_{r}$, starting from $r=0$, up to, and including radius, r. From Equation, (2-11), we know that,

${M}_{r}={M}_{R}{F}_{r}^{CDF}$ (5-2)

Table 4. For a blackhole having 1 solar mass.

Table 5. For a blackhole having 10 solar masses.

Table 6. For a blackhole having 10^{6} solar masses.

Thus, for a determination of,
${M}_{r}$, all that is needed is for us to multiply the total mass of the black hole,
${M}_{R}$, by the cdf, specified in the^{ }2^{nd} column, in Tables 4-6. These values are entered in the third column in the corresponding tables. It is to be noticed, upon inspection of the results, that most of the radiative mass is concentrated relatively close to the center of the black hole. This holds for all three black holes considered. This is a consequence of the pdf chosen, a truncated Gaussian, and the fact that the values indicated in Equations, (3-21), are relatively low in value. In other words, the dispersion is low.

In step 3, we seek to find the total radiative energy, ${E}_{r}$, starting from, $r=0$, up to, and including radius, r. We are also interested in finding the total internal energy, ${U}_{r}$, the total heat trapped, ${Q}_{r}$, and the total work done against radiative pressure in assembling the black hole, ${W}_{r}$. The quantities are to be evaluated, starting from, $r=0$, up to, and including radius, r. We start with, ${E}_{r}$. To find the total radiative energy, we integrate over the energy density function, but 4-dimensionally, from, $r=0$, to a final value, r. Within an infinitesimally small, 4-D volume element, $\text{d}{V}_{r}={A}_{r}\text{d}r=2{\text{\pi}}^{2}{r}^{3}\text{d}r$, we will have trapped electromagnetic energy in the amount, $\text{d}{E}_{r}={e}_{r}\text{d}{V}_{r}={\rho}_{r}{c}^{2}\text{d}{V}_{r}$. This is stored energy in a 4-D layer, between radius, r, and radius, $r+\text{d}r$. We integrate to find the total energy contained up to radius, r. The result is,

${E}_{r}={\displaystyle {\int}_{0}^{r}\text{d}{E}_{r}}={\displaystyle {\int}_{0}^{r}{e}_{r}\text{d}{V}_{r}}={\displaystyle {\int}_{0}^{r}{\rho}_{r}{c}^{2}\text{d}{V}_{r}}={M}_{r}{c}^{2}$ (5-3)

The mass of the black hole is thought to be made up of radiative energy, in all its forms, internal, plus pressure, plus heat. For the last equality, we made use of Equation, (2-4). If we specialize Equation, (5-3), to $r=R$, we obtain the full mass, ${M}_{R}$, for the mass of the black hole on the right hand side. See Equation, (2-5). This is the same mass that we observe three-dimensionally.

The next quantity to be determined is the internal energy, ${U}_{r}$. We follow the same steps as before. Within an infinitesimal volume element, $\text{d}{V}_{r}={A}_{r}\text{d}r$, we must have internal energy in the amount, $\text{d}{U}_{r}={u}_{r}{A}_{r}\text{d}r$. We integrate to find

${U}_{r}={\displaystyle {\int}_{0}^{r}\text{d}{U}_{r}}={\displaystyle {\int}_{0}^{r}{u}_{r}\text{d}{V}_{r}}={\displaystyle {\int}_{0}^{r}\frac{4}{10}{e}_{r}\text{d}{V}_{r}}=\frac{4}{10}{M}_{r}{c}^{2}$ (5-4)

For the 3^{rd} equality, we recognized that,
${u}_{r}=\frac{4}{10}{e}_{r}$. See the last equality in Equations, (2-22). We proceed similarly for the heat trapped, and the work done against radiative pressure. For the former, we find,

${Q}_{r}={\displaystyle {\int}_{0}^{r}\text{d}{Q}_{r}}={\displaystyle {\int}_{0}^{r}{q}_{r}\text{d}{V}_{r}}={\displaystyle {\int}_{0}^{r}\frac{1}{2}{e}_{r}\text{d}{V}_{r}}=\frac{1}{2}{M}_{r}{c}^{2}$ (5-5)

And for the latter, we obtain,

${W}_{r}={\displaystyle {\int}_{0}^{r}\text{d}{W}_{r}}={\displaystyle {\int}_{0}^{r}{p}_{r}\text{d}{V}_{r}}={\displaystyle {\int}_{0}^{r}\frac{1}{10}{e}_{r}\text{d}{V}_{r}}=\frac{1}{10}{M}_{r}{c}^{2}$ (5-6)

We made use of the relations in Equations, (2-22), in order to relate heat density, and radiative pressure, to, ${e}_{r}$. If we specialize the Equations, (5-4), (5-5), and (5-6), to $r=R$, we see that,

${U}_{R}=\frac{4}{10}{M}_{R}{c}^{2},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{Q}_{R}=\frac{1}{2}{M}_{R}{c}^{2},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{W}_{R}=\frac{1}{10}{M}_{R}{c}^{2}$ (5-7)

The factors of $\frac{4}{10},\frac{5}{10}$, and, $\frac{1}{10}$, in all these relations remain the same,

whether we consider, a layer, an enclosing volume, or the total volume. The factors are imposed upon us by the dimensionality of space and the fact that we are dealing with blackbody radiation, period. What is obvious with these relations is the fact that the internal energy makes up 40% of the total energy, the trapped heat represents 50% of the total energy, and the work done against pressure makes up 10% of the total. In 3-D space, the proportions would be different. In

three dimensions, the corresponding proportions would be, $\frac{3}{8}=37.5\%$, $\frac{4}{8}=50\%$, and, $\frac{1}{8}=12.5\%$.

Whenever we integrate pressure over volume, we obtain work done. Equation, (5-6), is particularly noteworthy from that point of view. We interpret this equation as the work done in assembling the black hole against radiative pressure, which will want to blow the black hole out. What keeps this radiative mass held together (contained) must be gravity. Our hydrostatic equation shows that. Layer by layer, gravity holds the radiative forces in check. See Equation, (2-26). In a sense, a black hole is something like a 4-D version of ball lightning, but because of the 4-D space, it is nonluminous. See the waterfall model described in reference [4] . We can think of it as an exotic type of 4-D capacitor, filled with blackbody radiative energy in a multi-layered fashion.

Coming back to our tables, we list ${U}_{r},{Q}_{r}$, and, ${W}_{r}$, under columns 4, 5 and 6, respectively. In each of our Tables 4-6, we have calculated these values. The total (radiative mass) energy is not tabulated. This can be found most easily, by multiplying, the radiative mass, entered under column 3, by ${c}^{2}$. Another way is to sum columns 4, 5, and 6, which gives the same results.

The next quantity to be considered is the product, ${G}_{r}{M}_{r}={G}_{r}^{\left(4\right)}{M}_{r}^{\left(4\right)}$, which brings us to step 4. In this step, we introduce the 4-D gravitational “constant”, ${G}_{r}={G}_{r}^{\left(4\right)}$, which is not really a constant as it will depend on the 4-D radius. It is model dependent, unlike Newton’s constant, ${G}_{N}$. In other words, its value depends on the pdf chosen. To determine the product of, ${G}_{r}{M}_{r}$, we will use Equation, (2-30). This can be rewritten in the form,

${G}_{r}{M}_{r}{\rho}_{r}=-\left(\text{d}{p}_{r}/\text{d}r\right){r}^{3}-3{p}_{r}{r}^{2}$ (5-8)

On the left hand side of this equation, we express the radiative mass density, ${\rho}_{r}$, in terms of the radiative pressure, ${p}_{r}$. From Equation, (2-21), it can be seen that,

${\rho}_{r}{A}_{r}=10{p}_{r}/{c}^{2}{A}_{r}$ (5-9)

Moreover, using Equation, (4-21), it is also true that,

${p}_{r}{A}_{r}/\left({p}_{R}{A}_{R}\right)=\text{exp}\left[\left({R}^{2}-{r}^{2}\right)/2{\sigma}^{2}\right]$ (5-10)

We differentiate both left and right hand sides of this equation with respect to radius, r. The result is

$\left[\left(\text{d}{p}_{r}/\text{d}r\right){A}_{r}+{p}_{r}\left(\text{d}{A}_{r}/\text{d}r\right)\right]/\left({p}_{R}{A}_{R}\right)=\text{exp}\left[\left({R}^{2}-{r}^{2}\right)/2{\sigma}^{2}\right]\left(-r/{\sigma}^{2}\right)$ (5-11)

Next, we divide the left hand side of Equation, (5-11), by the left hand side of Equation, (5-10), and we do the same on the right hand side. We obtain in this way, the new result

$\left(1/{p}_{r}\right)\left(\text{d}{p}_{r}/\text{d}r\right)+\left(1/{A}_{r}\right)\left(\text{d}{A}_{r}/\text{d}r\right)=\left(-r/{\sigma}^{2}\right)$ (5-12)

However, $\left(1/{A}_{r}\right)\left(\text{d}{A}_{r}/\text{d}r\right)=3/r$. And therefore, Equation, (5-12), can be rewritten as,

$\left(\text{d}{p}_{r}/\text{d}r\right)=\left(-{p}_{r}\right)\left(3/r+r/{\sigma}^{2}\right)$ (5-13)

We next substitute both Equations, (5-9), and, (5-13), into Equation, (5-8). The relation, (5-9), is used to eliminate, ${\rho}_{r}$, on the left side of Equation, (5-8), whereas relation,(5-13), will eliminate the $\left(\text{d}{p}_{r}/\text{d}r\right)$ term on the right hand side. In this way, we obtain,

${G}_{r}{M}_{r}10{p}_{r}/{c}^{2}={p}_{r}\left(3/r+r/{\sigma}^{2}\right){r}^{3}-3{p}_{r}{r}^{2}$ (5-14)

Two terms cancel on the right hand side, and we divide the remaining expression by, ${p}_{r}$. After rearranging terms, we find the much simpler result,

${G}_{r}{M}_{r}=0.1{c}^{2}\left({r}^{4}/{\sigma}^{2}\right)$ (5-15)

This is a very interesting expression because it depends on radius, and the shape parameter, $\sigma $.

Now, for the three black holes considered we know the shape parameters. These are given in Equations, (3-20). Therefore, we can evaluate, ${G}_{r}{M}_{r}$, as a function of radius, using the above Equation, (5-15). The results are present under column, 7, in each of the three Tables 4-6.

This brings us to step 5. In this step, we calculate, ${G}_{r}={G}_{r}^{\left(4\right)}$. From Equation, (2-11), we saw that that the radiative mass, ${M}_{r}$, at radius, r, is proportional to the cumulative distribution function, ${F}_{r}^{CDF}$, by means of the relation, ${M}_{r}={F}_{r}^{CDF}{M}_{R}$. Moreover, ${M}_{R}$, being the total mass of the black hole, can be related to its total radius, R, by way of the Schwarzschild relation, $R=2G{M}_{R}/{c}^{2}$. We can therefore write in place of ${M}_{r}$, the equation, ${M}_{r}={F}_{r}^{CDF}{M}_{R}={F}_{r}^{CDF}\left(R{c}^{2}/2G\right)$. We substitute this for, ${M}_{r}$ on the left hand side of Equation, (5-15). This gives us,

${G}_{r}{F}_{r}^{CDF}\left(R{c}^{2}/2G\right)=0.1{c}^{2}\left({r}^{4}/{\sigma}^{2}\right)$ (5-16)

Simplifying, and rearranging terms, it follows that,

${G}_{r}{F}_{r}^{CDF}=0.2\left(G/R\right)\left({r}^{4}/{\sigma}^{2}\right)$ (5-17)

This equation shows us the relation between the 4-D gravitational “strength”, ${G}_{r}$, and Newton’s constant, G. If we consider the special case where, $r=R$, this reduces to,

${G}_{R}=0.2G\left({R}^{3}/{\sigma}^{2}\right)$ (5-18)

At, $r=R$, we have a sharp discontinuity in gravitational strength due to the fact that $0.2\left({R}^{3}/{\sigma}^{2}\right)\ne 1$. In other words the gravitational strength does not match at the event horizon, when we compare its value on the 3-D side versus on the 4-D side. As one enters the 4-D black hole from 3-D space, the gravitational strength will increase abruptly in value, from, G, to ${G}_{R}={G}_{R}^{\left(4\right)}$. However, as argued in the previous paper [5] , this is to be expected. Many physical quantities will change their value abruptly due to the discontinuity of space. Gravity is no exception. We referred to this state of affairs in our previous work as our “waterfall model”. We emphasize, once more, that the results indicated by Equations, (5-15), and, (5-18), are pdf dependent. If we had introduced a different pdf, other than the truncated Gaussian, to model the interior, radiative mass, distribution of the black hole, then we would have obtained different expressions for the strength of gravity both within, and on the surface, of the black hole.

We also notice, from Equation, (5-17), that in the limit where, $r\to 0$, both the left hand side and right hand side vanish. In this limit, ${F}_{r}^{CDF}\to 0$, as can be seen from Equation, (2-12). At the very center of the black hole, it would appear as if there is no gravitational strength. In actual fact, it can take on any value in this singular limit, as Equation, (5-17), would still be satisfied. However, as soon as the radius becomes finite, not any value is allowed; Equation, (5-17), determines a specific value for, ${G}_{r}$. For a given size black hole, where we have a definitive $\sigma $ value, the gravitational strength will increase as, ${G}_{r}~{r}^{4}$, as indicated by relation, (5-17). For each of the black hole masses considered, we have the corresponding R and $\sigma $ values. See Equations, (3-18), and, (3-20). These values can be substituted in Equation, (5-17), to find ${G}_{r}$, for various r values. We also need the ${F}_{r}^{CDF}$ values, but these are specified under column 2, in each of the tables. The calculated ${G}_{r}$ values, using Equation, (5-17), are listed under column 8, in each of the Tables 4-6.

It would have also been possible to obtain, ${G}_{r}$, using the product, ${G}_{r}{M}_{r}$, specified under column 7, and dividing out the corresponding ${M}_{r}$ terms, which are indicated under column 3. However, this approach would not have given us explicit mathematical expressions for, ${G}_{r}$.

In step 6, we determine the gravitational acceleration, ${g}_{r}={g}_{r}^{\left(4\right)}$, within the 4-D black hole. This is most easily accomplished using a 4-D version of Gauss’s law. In 4-D space, we generalize Gauss’s law to read,

$\oint {g}_{r}\cdot \text{d}{A}_{r}}=-2{\text{\pi}}^{2}{r}^{3}{g}_{r}={M}_{r}/{\epsilon}_{r$ (5-19)

The ${M}_{r}$ stands for the radiative mass which is enclosed within radius, r, and we have assumed spherical symmetry. The negative sign is due to ${g}_{r}$ pointing radially inwards. We have introduced a “gravitational permittivity”, ${\epsilon}_{r}$, on the right hand side. As we shall see, this quantity will be radius dependent, and hence the subscript, r. This is the gravitational permittivity to be associated with the photon-packed blackbody medium. The permittivity really refers to the permittivity of a material substance, which, in this instance, is blackbody radiation at a specific temperature. Hence, it is radius dependent.

We rewrite Equation, (5-19), as

${g}_{r}={g}_{r}^{\left(4\right)}=-\left(\frac{1}{2{\text{\pi}}^{2}{\epsilon}_{r}}\right){M}_{r}/{r}^{3}=-{G}_{r}{M}_{r}/{r}^{3}$ (5-20)

where,

${G}_{r}={G}_{r}^{\left(4\right)}\equiv 1/\left(2{\text{\pi}}^{2}{\epsilon}_{r}\right)$ (5-21)

We can think of Equation, (5-21), as the defining equation for ${\epsilon}_{r}$. Because, ${G}_{r}$, is a function of radius, r, so too, is, ${\epsilon}_{r}$. This makes sense as an energetically denser medium should be reflected in the gravitational permittivity of the medium. Moreover, Equation, (5-20), seems a logical 4-D extension of the familiar 3-D equation, ${g}^{\left(3\right)}=-G{M}_{r}/{r}^{2}$.

We next substitute Equation, (5-15), into the right hand side of Equation, (5-20), and cancel terms. We thereby obtain a relatively simple result,

${g}_{r}={g}_{r}^{\left(4\right)}=-0.1{c}^{2}r/{\sigma}^{2}=-0.2G{M}_{R}/R\left(r/{\sigma}^{2}\right)$ (5-22)

For the second equality, we made use of the Schwarzschild equation. Equation, (5-22), is a very elegant result. It tells us that the gravitational acceleration increases linearly with radius, r, within the 4-D black hole of a specified mass, or radius. This is so in spite of the fact that the medium does not have uniform density. Also remember that the mass, or radius, of the black hole will determine its shape parameter, $\sigma $, as was seen in Section 3. As such, $\sigma $, is, effectively, a constant. Equation, (5-22), is specific to the pdf chosen, which, in our situation, is a truncated Normal distribution. At the event horizon where the radius, $r=R$, we see that, ${g}_{R}=-0.1{c}^{2}R/{\sigma}^{2}$. Equations, (5-20), and, (5-22), are to be noted for their simplicity. A different pdf would have led to different results.

We have calculated the gravitational acceleration, using Equation, (5-22), for the three black holes we are analyzing. The results are presented in Tables 4-6, under column 9. The gravitational acceleration will have its greatest value at the event horizon, much like the Earth.

For step 7, we calculate the gravitational potential, ${\phi}_{r}$, as a function of radius, r, which also will include, $r=R$. We know that the gravitational acceleration is negative the gradient of the gravitational potential. If we have spherical symmetry, then, ${g}_{r}=-\text{d}{\phi}_{r}/\text{d}r$, where ${g}_{r}$ is directed radially in. Therefore, $\text{d}{\phi}_{r}=-{g}_{r}\text{d}r$, and we can integrate. For the gravitational acceleration, ${g}_{r}$, we will use the first equality in Equation, (5-22). Thus, upon substitution,

${\int}_{0}^{r}\text{d}{\phi}_{r}}={\displaystyle {\int}_{0}^{r}\left(0.1{c}^{2}\right)\left(r/{\sigma}^{2}\right)\text{d}r$

${\phi}_{r}-{\phi}_{0}=\left(0.1{c}^{2}\right)\left({r}^{2}/2{\sigma}^{2}\right)$ (5-23)

This can also be rewritten as,

${\phi}_{r}-{\phi}_{0}=\left(0.2G{M}_{R}/R\right)\left({r}^{2}/2{\sigma}^{2}\right)$ (5-24)

because by the Schwarzschild condition, ${c}^{2}=2G{M}_{R}/R$. At the particular radius, $r=R$, Equations, (5-23), and, (5-24), become,

${\phi}_{R}={\phi}_{0}+\left(0.1{c}^{2}\right)\left({R}^{2}/2{\sigma}^{2}\right)={\phi}_{0}+\left(0.1G{M}_{R}R/{\sigma}^{2}\right)$ (5-25)

The value, ${\phi}_{0}$, is our constant of integration.

We want the 4-D gravitational potential to match the 3-D gravitational potential at, $r=R$. This will guarantee an equipotential surface for the event horizon, which has exactly the same value in both 3-D and 4-D space. For this purpose, we will demand that, ${\phi}_{R}^{\left(4\right)}={\phi}_{R}^{\left(3\right)}$, at radius, R. Using the second expression in Equation, (5-25), we obtain

${\phi}_{0}+\left(0.1G{M}_{R}R/{\sigma}^{2}\right)=-G{M}_{R}/R$ (5-26)

This fixes the constant of integration to equal,

${\phi}_{0}=-G{M}_{R}/R\left[1+0.1{R}^{2}/{\sigma}^{2}\right]$ (5-27)

Coming back to our original expression, (5-24), we see that this can now be rewritten as,

${\phi}_{r}=\left(-G{M}_{R}/R\right)\left[1+0.1\left({R}^{2}-{r}^{2}\right)/{\sigma}^{2}\right]$ (5-28)

This is our expression for the gravitational potential, within the black hole, and on its surface. We note that when, $r=R$, we retrieve our familiar, 3-D expression for the gravitational potential.

Utilizing Equation, (5-28), we can find ${\phi}_{r}$ for any value of radius, r. The right hand side can easily be evaluated for a specific black hole mass. Given the black hole radius, or mass, we can also find, $\sigma $, using the methods illustrated in Section 3. We have found the shape parameter for three black holes; they are indicated in Equations, (3-20). For these very specific black holes, we have calculated the gravitational potential as a function of r. The values are listed under column 10, in Tables 4-6. The gravitational potential is most negative at the very interior of the black hole.

The next step is step 8. In this step, we will evaluate the gravitational force. This calculation is very significant is that it will allow us to determine the total surface tension, which is needed to find the total work done in expanding the event horizon. See reference [5] . To determine the gravitational force, we will make use of the gravitational acceleration. If we consider a layer of thickness, dr, between radii, r, and, $r+\text{d}r$, the gravitational force acting on that layer, must be,

$\begin{array}{c}\text{d}{F}_{G,r}=\left({g}_{r}\right)\left(\text{d}{M}_{r}\right)=\left({g}_{r}\right)\left({\rho}_{r}\text{d}{V}_{r}\right)\\ =\left({g}_{r}\right)\left({\rho}_{r}{A}_{r}\text{d}r\right)=\left(-0.1{c}^{2}r/{\sigma}^{2}\right)\left({M}_{R}{f}_{r}\text{d}r\right)\end{array}$ (5-29)

For the 2^{nd} equality, we have made use of Equations, (5-22), and, (2-25). To find the total gravitational force, at a specific radius r within the black hole, we integrate Equation, (5-29). We obtain as a result,

$\begin{array}{c}{F}_{G,r}={\displaystyle {\int}_{0}^{r}\left({g}_{r}\right)\left({\rho}_{r}{A}_{r}\text{d}r\right)}={\displaystyle {\int}_{0}^{r}\left(-0.1{c}^{2}r/{\sigma}^{2}\right)\left({M}_{R}{f}_{r}\text{d}r\right)}\\ =\left(-0.1{M}_{R}{c}^{2}\right){\displaystyle {\int}_{0}^{r}\left(r/{\sigma}^{2}\right)\left({f}_{r}\text{d}r\right)}\\ =\left(-0.1{M}_{R}{c}^{2}\right){\displaystyle {\int}_{0}^{r}\left(r/{\sigma}^{2}\right)\left({f}_{R}\mathrm{exp}\left[\frac{{R}^{2}-{r}^{2}}{2{\sigma}^{2}}\right]\text{d}r\right)}\\ =\left(-0.1{M}_{R}{c}^{2}{f}_{R}\right){\displaystyle {\int}_{0}^{r}\left(r/{\sigma}^{2}\right)\left(\mathrm{exp}\left[\frac{{R}^{2}-{r}^{2}}{2{\sigma}^{2}}\right]\text{d}r\right)}\\ =-{F}_{R}\left[{\text{e}}^{{R}^{2}/2{\sigma}^{2}}-{\text{e}}^{\left({R}^{2}-{r}^{2}\right)/2{\sigma}^{2}}\right]\end{array}$ (5-30)

For the 3^{rd} line, we have substituted Equation, (4-21). And for the 5^{th} line, we have used Equation, (3-10), to simplify the constant factor in front of the integral.

We notice that Equation, (5-30), equals zero when r = 0, as it must. We also notice that, if $r=R$, then we obtain,

${F}_{G,R}=-{F}_{R}\left[{\text{e}}^{{R}^{2}/2{\sigma}^{2}}-1\right]$ (5-31)

Both Equations, (5-30), and (5-31), are rather simple. As always, they are a direct consequence of the pdf chosen. We can use these equations directly to find the 4-D gravitational force as a function of radius. We know, the, $\sigma $, as well as the, R, values for the three black holes under consideration. See Equations, (3-20), and, (3-18). We also know the value for the radiative force, ${F}_{R}={F}_{R}^{\left(4\right)}$, just inside the event horizon. These are given by Equations, (4-17), for the three black holes. The values for the gravitational force are tabulated under column 11 in each of the Tables 4-6.

In the first paper [5] , we spent a considerable of time talking about the significance of the surface tension. In fact, we noted that it consisted of two components, a radiative as well as a gravitational component,

${F}_{\text{ST}}=\left({F}_{\text{ST}\text{\hspace{0.17em}}\text{radiative}}+{F}_{G,R}^{\left(4\right)}\right)={F}_{R}^{\left(3\right)}-{F}_{R}^{\left(4\right)}+{F}_{G,R}^{\left(4\right)}$ (5-32)

See Equations, (3-10), or, (3-13), in reference [5] . The radiative force, ${F}_{R}^{\left(3\right)}$, is the 3-D radiative force pushing the event horizon in, whereas, ${F}_{R}^{\left(4\right)}$, is the 4-D radiative force pushing the event horizon out. The surface tension is important in determining the work needed in expanding the black hole. Very simply, to increase the radius of a black hole by an amount, dR, from surface radius, R, to, $R+\text{d}R$, it takes the following amount of work, $\text{d}W={F}_{\text{ST}}\text{d}R$. We compare the radiative surface tension in Equation, (5-32), to the gravitational surface tension within the same equation. From reference [5] , we obtained

${M}_{BH}=\left({M}_{\text{sun}},10{M}_{\text{sun}},{10}^{6}{M}_{\text{sun}}\right)$

$\to {F}_{\text{ST}\text{\hspace{0.17em}}\text{radiative}}=\left(3.05\times {10}^{-7},3.05\times {10}^{-5},3.05\times {10}^{+5}\right)$ (5-33)

Refer to Equation, (3-19), specifically, in the work cited. These values hold for an isolated, static black hole, where there is virtually no net radiative heat coming in. We compare these values to the gravitational force values, indicated in Tables 4-6. There we obtain, at the surface,

${M}_{BH}=\left({M}_{\text{sun}},10{M}_{\text{sun}},{10}^{6}{M}_{\text{sun}}\right)$

$\to {F}_{\text{ST}\text{\hspace{0.17em}}\text{gravity}}={F}_{G,R}=\left(7.30\times {10}^{43},7.16\times {10}^{45},6.36\times {10}^{43}\right)$ (5-34)

In both equations, forces are measured in Newtons. We see that there is no comparison. Radiative forces cannot compete against gravitational forces when expanding a black hole in this situation, where there is virtually no net radiative heat inflow.

If there is substantial net heat inflow, the radiative surface tension contribution may not negligible. Consider an external, outside temperature of, ${T}_{2}={10}^{9}\text{\hspace{0.17em}}\text{K}$, which would correspond to hard X-rays, having energies in the 100 keV range. We can work out the radiative forces, in terms of radius, R. We obtain using our, by now, familiar relations,

${F}_{R}^{\left(3\right)}={p}_{2}^{\left(3\right)}{A}_{R}^{\left(3\right)}=\frac{1}{3}{a}^{\left(3\right)}{\left({10}^{9}\right)}^{4}4\text{\pi}{R}^{2}=3.17\times {10}^{21}{R}^{2}$ (5-35)

${F}_{R}^{\left(4\right)}={p}_{1}^{\left(4\right)}{A}_{R}^{\left(4\right)}=\frac{1}{4}{a}^{\left(4\right)}{\left(0.0597\times {R}^{-1}\right)}^{5}2{\text{\pi}}^{2}{R}^{3}=1.40\times {10}^{-13}{R}^{2}$ (5-36)

We have used Equations, (4-8a,b), (2-16a), (2-17) and, (4-9a), to obtain the above. The 3-D internal energy density can be expressed as, ${u}_{R}^{\left(3\right)}={a}^{\left(3\right)}{T}^{4}=\left(7.5657\times {10}^{-16}\right){T}^{4}$, in MKS units, a well-known result. The radiative surface tension is defined as, ${F}_{\text{ST}\text{\hspace{0.17em}}\text{radiative}}\equiv {F}_{R}^{\left(3\right)}-{F}_{R}^{\left(4\right)}$. For the three black holes, we therefore find, using the definition for radiative surface tension, and Equations, (5-35), with, (5-36),

${M}_{BH}=\left({M}_{\text{sun}},10{M}_{\text{sun}},{10}^{6}{M}_{\text{sun}}\right)$

$\to {F}_{\text{ST}\text{\hspace{0.17em}}\text{radiative}}=\left(2.77\times {10}^{28},2.77\times {10}^{30},2.77\times {10}^{40}\right)$ (5-37)

We have used the radii listed in Equations, (3-18). These values in Equations, (5-37), are still much less than those indicated by Equation, (5-34). However, for very massive black holes ( ${M}_{BH}>{10}^{6}{M}_{\text{sun}}$ ), with external temperatures very high ( ${T}_{2}>{10}^{9}\text{\hspace{0.17em}}\text{K}$ ), we can see that we are approaching a situation where the radiative contributions to surface tension can no longer simply be ignored.

Our final step is step 9. In this step, we evaluate the entropy of a black hole, ${S}_{r}$, as a function of radius. We treat the entropy as an extrinsic quantity, dependent on volume. Our volume, however, is a 4-dimensional one. In a volume element, $\text{d}{V}_{r}$, between r and $r+\text{d}r$, we can expect a value for entropy in the amount, $\text{d}{S}_{r}={s}_{r}\text{d}{V}_{r}$, where ${s}_{r}$ is the entropy density in layer, r and $r+\text{d}r$. The entropy up to radius, r, can be found by integrating this expression. We find that,

${S}_{r}={S}_{r}^{\left(4\right)}={\displaystyle {\int}_{0}^{r}{s}_{r}\text{d}{V}_{r}}={\displaystyle {\int}_{0}^{r}{s}_{r}{A}_{r}\text{d}r}$ (5-38)

The, ${A}_{r}=2{\text{\pi}}^{2}{r}^{3}$, is the 4-D surface area at radius, r. On the left hand side, we have fixed the constant of integration, ${S}_{0}=0$, as we believe that there can be no entropy if the volume is zero.

We know, however, that the entropy density can be related to radiative pressure, when multiplied by a volume or surface area. From Equations, (2-18), we can write,

${s}_{r}{A}_{r}=5/4{u}_{r}{A}_{r}/{T}_{r}=5{p}_{r}{A}_{r}/{T}_{r}$ (5-39)

This expression will be substituted into Equation, (5-29), and we obtain,

${S}_{r}={\displaystyle {\int}_{0}^{r}{s}_{r}{A}_{r}\text{d}r}={\displaystyle {\int}_{0}^{r}5{p}_{r}{A}_{r}/{T}_{r}\text{d}r}=5{\displaystyle {\int}_{0}^{r}{F}_{r}/{T}_{r}\text{d}r}$ (5-40)

For the last equality, we employed Equation, (2-23), which relates radiative force to radiative pressure.

Next, we utilize, Equation, (4-23). Substituting this for, ${F}_{r}$, allows us to re-express Equation, (5-40), as follows.

${S}_{r}=5{\displaystyle {\int}_{0}^{r}{F}_{r}/{T}_{r}\text{d}r}=5{F}_{R}{\displaystyle {\int}_{0}^{r}\text{exp}\left[\left({R}^{2}-{r}^{2}\right)/2{\sigma}^{2}\right]/{T}_{r}\text{d}r}$ (5-41)

We will also make use of relation, (4-25), which we will rewrite as,

$1/{T}_{r}=1/{T}_{R}{\left(r/R\right)}^{3/5}\text{exp}\left[-\left({R}^{2}-{r}^{2}\right)/10{\sigma}^{2}\right]$ (5-42)

We substitute this Equation, (5-42), into the integrand of Equation, (5-41), in order to eliminate the $1/{T}_{r}$ term. Upon simplifying, we find that

${S}_{r}=\left(5{F}_{R}/{T}_{R}\right)\left(1/{R}^{3/5}\right){\displaystyle {\int}_{0}^{r}{r}^{3/5}\text{exp}\left[0.4\left({R}^{2}-{r}^{2}\right)/{\sigma}^{2}\right]\text{d}r}$ (5-43)

This is an expression, which can be integrated. The result is

${S}_{r}=\left(-{F}_{R}/{T}_{R}\right){R}^{-\frac{3}{5}}{\left(5/2\right)}^{\frac{9}{5}}\left[\Gamma \left(\frac{4}{5},\frac{2{r}^{2}}{5{\sigma}^{2}}\right)-\Gamma \left(\frac{4}{5}\right)\right]{\sigma}^{8/5}{\text{e}}^{\frac{2{R}^{2}}{5{\sigma}^{2}}}$ (5-44)

In this equation, $\Gamma \left(s\right)$, is the gamma function and, $\Gamma \left(s,x\right)$, is the upper incomplete gamma function, defined in the Appendix. Using the identities given in the Appendix, we can re-express Equation, (5-44), in terms of the lower incomplete gamma function, $\gamma \left(s,x\right)$. The result is,

${S}_{r}=\left(+{F}_{R}/{T}_{R}\right){R}^{-\frac{3}{5}}{\left(5/2\right)}^{\frac{9}{5}}\gamma \left(\frac{4}{5},\frac{2{r}^{2}}{5{\sigma}^{2}}\right){\sigma}^{8/5}{\text{e}}^{\frac{2{R}^{2}}{5{\sigma}^{2}}}$ (5-45)

This gives the entropy as a function of radius within the black hole. In hindsight, Equation, (5-45), could have been obtained by using Equation, (4-27). This would have been a little more direct, but then we would not have seen the connection between temperature and radius, as specified in Equation, (5-42). And the quantity, ${s}_{R}$, would also have to be worked out.

The lower incomplete gamma function, $\Gamma \left(s,x\right)$, in Equation, (5-45), can be evaluated for various radii within the black hole since we know the $\sigma $ value. We also know the radiative force, ${F}_{R}$, just inside the black hole, as well as its radius. Finally, the temperature, ${T}_{R}={T}_{1}$, is also known. Therefore, the entropy, ${S}_{r}$, can be calculated for each of the black holes considered. The results are presented under column 12 in the Tables 4-6. For the various radii, Equations, (3-18), were utilized. For the shape parameters, Equations, (3-20), were employed. The corresponding 4-D inside surface temperatures were specified by Equations, (4-10). And in relations, (4-17), we have the respective radiative forces, just inside the event horizons for the black holes being analyzed. All the quantities in Equation, (5-36), are 4-D quantities.

We note that because of the mathematical inequality, $\Gamma \left(s\right)>\Gamma \left(s,x\right)$, the entropy, as given by expression, (5-35), will always be positive. We also give the Excel function for the lower incomplete gamma function, as this is somewhat unusual. To calculate, $\gamma \left(s,x\right)$, we may use

$\gamma \left(s,x\right)=EXP\left(GAMMALN\left(s\right)\right)\ast GAMMA.DIST\left(x,s,1,TRUE\right)$ (5-46)

Utilizing this expression, one can verify these results.

We are now in a position to discuss our results. We first notice upon inspection of Tables 4-6, that the radiative mass in each instance is very much concentrated within the core of the black holes. In fact, within a radius of 15% of the total radius, we have the following percentages of total mass contained, for the three black holes being considered:

$\left({M}_{r=0.15R}^{A},{M}_{r=0.15R}^{B},{M}_{r=0.15R}^{C}\right)=\left(97.7\%{M}_{R}^{A},97.4\%{M}_{R}^{B},95.2\%{M}_{R}^{C}\right)$ (5-47)

These percentages indicate that the radiative mass is very much localized within a relatively small core volume. The cumulative probability distribution function, ${F}_{r}^{CDF}$, listed under column 2, represents the actual percentage of total mass contained within a given radius. See Equation, (2-11).

The same percentages hold for the total radiative energy, ${E}_{r=0.15R}$, at $r=0.15R$, the internal energy, ${U}_{r=0.15R}$, at $r=0.15R$, and the heat energy, ${Q}_{r=0.15R}$, at $r=0.15R$. All these quantities are proportional to, ${M}_{r=0.15R}$, as can be seen by Equations, (4-22), and, (4-23). Once we know how radiative mass is distributed within the black hole, we effectively know how energy, in all its forms, is likewise distributed. We have chosen a specific pdf to model the mass distribution within the interior of the black hole. This fixes column 2 in all these tables. And from column 2, we derive the other quantities. Again, we can think of the black hole as an exotic type of capacitor in that it can store radiative energy. How specifically that energy is packed will depend on the pdf, which is chosen.

A third realization is that gravity varies within the black hole. Surprisingly, it is radius dependent, a consequence of the pdf chosen. Thus, in 4-D space, the 4-D gravitational “constant”,
${G}_{r}={G}_{r}^{\left(4\right)}$, is not a constant. Also at the event horizon, its value changes abruptly. Upon entering the black hole from 3-D space, its value increases dramatically from,
$G={G}_{N}$, which is Newton’s constant, to
${G}_{R}={G}_{R}^{\left(4\right)}$. We can set up a ratio of
${G}_{R}/G$. This ratio is not dimensionless as the gravitational “constant” in 4-D space has units of, Newtons m^{3}/kg^{2}. We find for the three black holes under investigation, the following increases upon entering the black hole form 3-D space:

$\left({G}_{R}^{A}/G,{G}_{R}^{B}/G,{G}_{R}^{C}/G\right)=\left(1.35\times {10}^{5},1.30\times {10}^{6},1.03\times {10}^{11}\right)\text{meters}$ (5-48)

The larger the black hole, the greater the increase upon entering 4-D space. For supermassive black holes, this increase is particularly impressive.

A fourth surprise is entropy. Contrary to current thinking, the entropy in our model suggests that a black hole may be a highly ordered state. Our entropy value, in relative terms, is small, very much smaller than the Bekenstein result. Also, our entropy scales with increasing radius, comparable to ${R}^{2}$, but not with a clean power law. Our expression for entropy is given by Equation, (5-45). At radius, $r=R$, which characterizes the event horizon, this quantity has been calculated in each of the tables. For the three black holes being analyzed, we find

$\left({S}_{R}^{A},{S}_{R}^{B},{S}_{R}^{C}\right)=\left(1.63\times {10}^{37},6.59\times {10}^{40},7.24\times {10}^{46}\right)$ (5-49)

The units for 4-D entropy are, Joule/K, in the MKS system. According to Bekenstein, the sun has an entropy of roughly, 10^{35} J/K. For a black hole having the mass of the sun, he calculates a value equal to, 1.5 × 10^{54} J/K. Our value for a black hole having the mass of the sun is different, 1.63 × 10^{37} J/K, which is about one hundred times larger than that of the sun, but much, much lower than his estimate. See the first equality in Equations, (5-49). Moreover, his entropy scales as,
${S}_{\text{Bekenstein}}~{R}^{2}$. Our entropy also scales, but the scaling behavior is complicated. Upon inspection of the values indicated by Equations, (5-49), there is much variation between the various masses, and corresponding radii. If our result is to be believed, most of the entropy in the universe is, thus, not in the form of black holes. As with other calculated global quantities, it needs to be emphasized that our formula for entropy depended on our choice for distribution function. If we had chosen a different pdf, other than a truncated Gaussian, the results for entropy would be entirely different.

As with radiative mass, most of the entropy, is to be found within the very interior of the black hole. Resorting to the tables, we can calculate how much entropy we have, as a percentage of the total, at a radius of, $r=0.15R$. We find,

$\left({S}_{r=0.15R}^{A},{S}_{r=0.15R}^{B},{S}_{r=0.15R}^{C}\right)=\left(91.2\%{S}_{R}^{A},90.4\%{S}_{R}^{B},84.9\%{S}_{R}^{C}\right)$ (5-50)

We see that we have accounted for 85%+ of the total entropy, within an inner radius, 15% of the total.

6. Summary and Conclusions

In this work, we analyzed what the internal structure of a black hole could look like, given our 4-D blackbody radiation model for a black hole. We treat the black hole as a self-contained, spherically symmetric, steady state, 4-D spatial ball filled with radiation. The blackbody radiation produces the radiative mass and we chose a particular pdf, the truncated Gaussian, to pack the mass a certain way within its interior. Any truncated pdf with support, $r\in \left[0,R\right]$, hypothetically, could have been chosen, but the truncated Gaussian is relatively simple, somewhat familiar, and leads to acceptable results. We fixed the parameters by assuming that we have maximum radiative energy at the center, which allowed us to set, $\mu =0$. We were also able to determine the shape parameter, $\sigma $, by using the temperature just inside the event horizon. This inner surface temperature was calculated using a generalized version of the Stefan-Boltzmann law, derived in a previous paper, and which allows for radiative transfers between spaces having different dimensions. In a very real sense, we can think of a black hole as an exotic type of capacitor, but one which stores energy in 4-D space, in a non-uniform way, and, in blackbody form. The distribution of radiative energy is not constant; rather it follows a pdf, which has to be specified and fulfill certain conditions. We treat the black hole as a 4-D construct, for reasons given in a previous work. If our model is correct, the rip or tear in the space-time continuum occurs not at the center of the black hole, but rather at its surface, because it is there that we are transitioning between 3-D and 4-D space. The event horizon is assumed infinitely thin.

In Section 2, we considered two stellar-like equations, conservation of mass and hydrostatic equilibrium. We generalized these equations to 4-D space. We discovered that the radiative mass density cannot be uniformly distributed. First, the temperature at the surface of the black hole, on the inside, would be too high, if we want to accommodate the required mass. Second, and perhaps more importantly, we would have no gravitational forces within its interior. For that to happen, a temperature gradient is necessary. The only gravity which would then exist would be at the surface. We therefore introduced a generalized Gaussian distribution, but truncated so that only values of radii within the black hole are admissible. The 4-D generalization of the hydrostatic equilibrium Equation, (2-14), was also considered. We found that this had to be modified because it is not just pressure which is changing when one changes the radius. We also have to change the surface area upon which this pressure is acting. Our modification to this equation is relation, (2-30). In Section 2, we argued that the internal energy density, the heat density and the radiative pressure all contributed towards the total energy/radiative mass density. A temperature change will carry over to all these variables, and not selectively to a subset. All these variables are defined in terms of a specific temperature, and if the temperature changes, so, too, do all these variables change, simultaneously.

In Section 3, we looked at the specifics of the truncated Gaussian distribution. We specified the pdf, and the corresponding cdf, for the case where the mode is assumed to be at the center, $r=0$. This fixed our first parameter, $\mu $, associated with this distribution. In this case where we have maximum radiative mass density at the center, $\mu =0$. The second parameter, $\sigma $, called the shape parameter, was determined using the inside surface temperature of the black hole. We determined three values for sigma, for three specific black holes. The results are given by Equations, (3-20). The black hole masses considered were highlighted in Equations, (2-18). For this choice of distribution, Equations, (3-6), and, (3-8), tell us how the radiative mass is distributed (packed) within the black hole. In this paper, we considered, by and large, only black holes for which there is no radiative heat inflow. This was to keep the discussion simple. We did indicate, however, how to generalize these results for the situation where there is net heat inflow. Equations, (2-14), (2-15), and, (2-16) apply if there is no net heat inflow. This allowed us to find the pdf value, ${f}_{R}$, at $r=R$. We then used Equation, (2-18), to fix the $\sigma $ value for each black hole having a different mass, or alternatively, radius.

Section 4 concerned itself with determining the temperature profile, the various densities, the radiative pressures, and the radiative forces within the black hole. All these quantities are functions of the 4-D radius, r. After a review of the generalized Stefan-Boltzmann law for a 3-D/4-D interface, we found, ${T}_{1}$, which is the temperature just inside the event horizon, on the 4-D side. For the three black holes considered in this paper, these temperatures were evaluated. They are given by Equations, (4-10). The calculated temperatures, specified by Equations, (4-10), are much higher than the corresponding Hawking temperatures, but then again, they are defined within a different spatial dimension, and are obtained by different means. The three black holes we focused on, numerically, were those having masses, and radii, as specified under Equations, (3-18).

Once we have the temperatures at the surface, we can use our probability distribution function, to determine the temperature profiles, ${T}_{r}=T\left(r\right)$ within the interior. This will also allow us to calculate, within the interior, the internal energy density, ${u}_{r}=u\left(r\right)$, the radiative pressure, ${p}_{r}=p\left(r\right)$ the heat density, ${q}_{r}=q\left(r\right)$, the total radiative mass density, ${\rho}_{r}=\rho \left(r\right)$, the entropy density, ${s}_{r}=s\left(r\right)$ and the radiative force, ${F}_{r}=F\left(r\right)$. We presented our results in table form, one for each black hole under consideration. Table 1, our first table, gave the results for black hole, $B{H}_{A}$, which has a mass equal to that of the sun. Table 2, holds for a black hole, $B{H}_{B}$, which has a mass, which is ten times that of the sun. And, Table 3 corresponds to a supermassive black hole, $B{H}_{C}$, which is a million times more massive than that of the sun. The equations that we used to calculate these values were given in the text.

In Section 5, we determined some other attributes associated with the three black holes, which were under investigation. Specifically, we focused on evaluating the total enclosed radiative mass, ${M}_{r}=M\left(r\right)$, up to and including radius, r. We also looked at the enclosed radiative internal energy, ${U}_{r}=U\left(r\right)$, the trapped heat, ${Q}_{r}=Q\left(r\right)$ and the work done against radiative pressure, ${W}_{r}=W\left(r\right)$, up to, and including, radius, r. These were given by Equations, (5-4), (5-5), and (5-6), respectively. The total radiative energy, ${E}_{r}=E\left(r\right)$, is the sum of these three components, and this is what produces the radiative mass in our model. For their individual evaluation, integration over 4-D space was necessary. Their values, as well as those for, ${M}_{r}$, are listed in Tables 4-6, for various radii within the black holes. Table 4 holds for black hole, $B{H}_{A}$, which has one solar mass. Table 5 and Table 6, apply for black holes, $B{H}_{B}$, and, $B{H}_{C}$, respectively, where black hole, $B{H}_{B}$, has a mass ten times that of the sun. And black hole, $B{H}_{C}$, refers to a black hole which is a million times more massive than the sun. We saw that these values were essentially proportional to the cumulative distribution function, ${F}_{r}^{CDF}$, which is indicated under column 2 in all these tables. Most of the radiative mass (85%+), and therefore energy, is to be found within a core radius that is 15% of the total radius.

In this section, we also evaluated/calculated the gravitational acceleration, ${g}_{r}={g}_{r}^{\left(4\right)}$, the gravitational strength, ${G}_{r}={G}_{r}^{\left(4\right)}$, the gravitational force, ${F}_{G,r}={F}_{G,r}^{\left(4\right)}$, and the gravitational potential, ${\phi}_{r}=\phi \left(r\right)$, within the black hole, as functions of radius, r. All these quantities were given by relatively simple expressions, a result of the pdf chosen, and they are also listed in the tables indicated above. We saw that upon entering the 4-D black hole, the gravitational “constant” is no longer a constant, but varies with respect to radius. Not only that, but its value increases abruptly upon entering 4-D space. The increases in gravitational “constant” over the Newtonian value is given by Equations, (5-48), which are the increases, for the three black holes under consideration.

Finally, in this section, we evaluated the entropy, ${S}_{r}=S\left(r\right)$, from, $r=0$, up to and including radius, r. This was given by a rather complicated function, Equation, (5-45), which involves a lower incomplete gamma function. This result is a consequence of the pdf chosen, a truncated Normal distribution. From Equation, (5-45), we notice that the entropy depends on the full radius, R, the inside surface temperature of the black hole, ${T}_{R}={T}_{1}$, the radiative force just inside the event horizon, ${F}_{R}={p}_{1}{A}_{R}$, and the shape parameter, $\sigma $. We have each of these quantities for the three black holes being analyzed. Thus, we can tabulate the entropy, ${S}_{r}$, as a function of radius. The results are presented in each of the Tables 4-6, under column 12. The total entropy, ${S}_{R}$, at radius, R, is especially interesting. We find that it is very much less than the corresponding Bekenstein result, especially for very large black holes. See Equations, (5-49), where specific values are calculated. Our results also scale with the radius of the black hole, but in a somewhat complicated manner. See Equation, (5-45), where the R dependency is indicated. According to Bekenstein, the entropy should scale as, ${S}_{\text{Bekenstein}}~{R}^{2}$, making it an intrinsic variable in 3-D space. In sharp contrast to current thinking, our conclusion is that the black hole may be a highly ordered state, if our model is to be taken seriously. Moreover, in Equations, (5-50), we observe that the entropy is pretty much determined within a radius, which is 15% of the full radius. Almost the entire entropy is contained within the core.

Our model suggests that black holes are, in fact, highly ordered states. In this context, it has been noted that black holes may be well-defined quantum mechanical systems, having ordered discrete quantum spectra [23] [24] [25] . Within this framework, black holes consist of quasi-normal modes (QNM’s), which are to be considered normal quantum levels within the black hole. The external structure of a black hole can thus be approximated as a quantized system analogous to the excited electron states within the Bohr model. As a consequence, evaporative emissions in the form of Hawking radiation would consist of pure states. This interesting model supports the notion that black holes are, in reality, highly ordered quantized states in a semi-classical approximation. We considered here in this paper a classical picture for black holes, but a quantum mechanical treatment would certainly be an interesting extension of this work.

We focused in this paper on one pdf, and one scenario, an isolated, static black hole, with no heat inflow. Obvious extensions would be to consider net radiative heat inflows, where $\text{d}Q/\text{d}t>0$. There, we would have to specify the specific amount of temperature on the outside of the event horizon. Another extension would be to consider other truncated pdf’s with the same support, $r\in \left[0,R\right]$, to model the interior radiative mass distribution inside the black hole. These may have more favorable attributes or characteristics. Our goal in this work was to develop a framework, for considering what could happen within the interior of a black hole. It was not, obviously, to play out all scenarios.

Another generalization of this work would be to add other forms of radiation, not just blackbody, and consider their effect within the interior of the black hole. We could also incorporate massive particles, obeying Fermi-Dirac statistics. If the temperature is above the so-called “freeze-out” temperature, these particles are in the form of radiation. Below the freeze-out temperature, we expect distinct massive particles. We have seen that the temperature within the black hole assume a distribution of temperatures. Thus, there could be regions within the black hole where we have particles, and other regions, where we have radiations. For a one solar mass black hole (see Table 1), the temperature decreases from, $T=4.14\times {10}^{12}\text{\hspace{0.17em}}\text{K}$, at a radius, $r={10}^{-6}R$, down to, $T=1.09\times {10}^{9}\text{\hspace{0.17em}}\text{K}$, at a radius, $r=0.2R$. These are high temperatures, which correspond to energies of, 357 MeV, down to, 94 keV, respectively. Above, 214 MeV, most of the particles in the standard model [1] [2] [3] [4] are in thermodynamic equilibrium with the blackbody photons. To take into account both fermionic and bosonic relativistic species, we must specify a scale factor, ${g}_{\text{TOTAL}}^{*}\left(T\right)$, which is temperature dependent. At about, 214 MeV, the scale factor, ${g}_{\text{TOTAL}}^{*}=62.5$. We can multiply the internal energy density of blackbody photons by a factor, such as this, to estimate what the black hole is capable of holding energy-wise, if we wish to take into account other forms of radiation. This is another area for further study.

Above background energies of, 1 GeV, which corresponds to a temperature of, 1.16 × 10^{13} K, the scale factor becomes,
${g}_{\text{TOTAL}}^{*}=75.5$. And at background energies of approximately, 1 TeV, which corresponds to a temperature of, 1.16 × 10^{16} K, the scale factor increases to,
${g}_{\text{TOTAL}}^{*}=106.75$. Above, 1 TeV, all particles in the standard model are relativistic. Interestingly such temperatures are only exceeded for incredibly small black holes. We saw that the smaller the black hole, the greater the temperature within its very core. As mentioned in the text, in passing, this may have ramifications for particle physics, and the modeling of elementary particles as mini black holes. This is a further area of investigation. One can consider the modeling of elementary particles as mini-black holes.

It is difficult to see, at present, how this model can be tested observationally. We can estimate outside surface temperatures by the electromagnetic emissions being given off next to a black hole. We can also estimate the black hole mass, and size, by the trajectory of stars in its vicinity. This should allow us to estimate inside surface temperature. Net heat flow can then be calculated using our generalized Stefan-Boltzmann equation. We have assumed that the heat flowing in is a direct function of the surface temperature on the outside, versus that on the inside. For a temperature on the inside, we assumed a-priori, that, for a specific radius, it can be calculated as if the black hole were an isolated, static black hole. This is a working assumption, which may or may not hold true. With all these assumptions, we might observe how a black hole expands with time, and see if our relations make sense in predicting radiative heat inflow with corresponding expansion. We believe that 4-D space is probably beyond our observational reach. Probing this space to verify the radial properties within the black hole may well prove to be the ultimate challenge.

Finally, relating to probing 4-D space, Hawking radiation could be considered. Hawking radiation was largely ignored in this work as it was argued that evaporative processes were 2^{nd} order effects. Our inside surface black hole temperatures were shown to be much higher than the corresponding Hawking temperatures, and thus evaporative processes were largely ignored/discounted. We calculated and compared temperatures in the text for various massive black holes. We also showed that the inside temperature is much, much smaller than the outside 3-D temperatures, with or without radiative inflow. (In brane theory, one has a similar result… the temperatures are higher in higher dimensional space.) However, Hawking radiation and evaporative processes should be looked at, and studied, more carefully in light of the peculiarities of our 4-D spatial model. If the inside surface temperatures are higher, this actually may have important ramifications for leakage.

As is known [26] , the spectrum of particle emissions in terms of energy flux is dependent on the type of black hole considered as well as type of particle emitted. For example, a non-rotating black hole emits spin zero particles pretty much like a perfect blackbody spectrum. A rotating black hole, on the other hand, emits a series of peaks, which is dependent on the rotation parameter and angular velocity. See, for example, Figure 3-3 and Figure 3-4, versus Figure 3-5, in the reference cited, where the latter figure refers to a spinning black hole. Furthermore, qualitatively, higher spin particles follow the same pattern. By considering evaporative processes, one might conceivably infer what is happing within the 4-D black hole. This is obviously a further line of inquiry, and research, which can, and should be pursued.

Acknowledgements

With these two papers, the author would like to thank his children, Kira, Juergen, and Henry, for their love and support. They have been a constant source of encouragement, inspiration, joy and goodwill, throughout my entire life. They are, first and foremost, my proudest achievement, and will always remain so.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

Cite this paper

Pilot, C. (2019) On the Internal Structure of a Black Hole Utilizing a 4-D Spatial Blackbody Radiation Model. Journal of High Energy Physics, Gravitation and Cosmology, 5, 719-772. https://doi.org/10.4236/jhepgc.2019.53039

References

- 1. Kolb, E.W. and Turner, M.S. (1989) The Early Universe. Addison-Wesley, Reading.
- 2. Baumann, D.D. (2015) Lecture Notes on Cosmology. http://theory.uchicago.edu/~liantaow/my-teaching/dark-matter-472/lectures.pdf
- 3. Mather, J.C., et al. (1999) Calibrator Design for the COBE Far-Infrared Absolute Spectrophotometer (FIRAS). The Astrophysical Journal, 512, 511-520. https://doi.org/10.1086/306805
- 4. Husdal, L. (2016) On Effective Degrees of Freedom in the Early Universe. Galaxies, 4, 78. https://doi.org/10.3390/galaxies4040078
- 5. Christopher, P. (2019) Are Black Holes 4-D Spatial Spheres Filled with Black Body Radiation? Generalization of the Stefan-Boltzmann Law and Young-Laplace Relation for Spatial Radiative Transfers. Journal of High Energy Physics, Gravitation, and Cosmology, 5, 638-682. https://doi.org/10.4236/jhepgc.2019.53036
- 6. Zeilik, M.A. and Gregory, S.A. (1998) Introductory Astronomy & Astrophysics. 4th Edition, Saunders College Publishing, Philadelphia.
- 7. Hansen, C.J., Kawaler, S.D. and Trimble, V. (2004) Stellar Interiors. 2nd Edition, Springer, Berlin.
- 8. Landsberg, P.T. and De Vos, A. (1989) The Stefan Boltzmann Constant in an N-Dimensional Space. Journal of Physics A: Mathematics and General, 22, 1073-1084. https://doi.org/10.1088/0305-4470/22/8/021
- 9. Menon, V.J. and Agrawal, D.C. (1998) Comment on “The Stefan-Boltzmann Constant in N-Dimensional Space”. Journal of Physics A: Mathematics and General, 31, 1109-1110. https://doi.org/10.1088/0305-4470/31/3/021
- 10. Barrow, J.D. and Hawthorne, W.S. (1990) Equilibrium Matter Fields in the Early Universe. Monthly Notices of the Royal Astronomical Society, 243, 608-609.
- 11. Gonzalez-Ayala, J., Perez-Oregon, J., Cordero, R. and Angulo-Brown, F. (2015) A Possible Cosmological Application of Some Thermodynamic Properties of the Black Body Radiation in N-Dimensional Euclidean Spaces. Entropy, 17, 4563-4581. https://doi.org/10.3390/e17074563
- 12. Gonzalez-Ayala, J., Cordero, R. and Angulo-Brown, F. (2016) Is the (3+1)-d Nature of the Universe a Thermodynamic Necessity? EPL (Europhysics Letters), 113, 4006. https://doi.org/10.1209/0295-5075/113/40006
- 13. Johnson, N.L., Kotz, S. and Balakrishnan, N. (1994) Continuous Univariate Distributions. Volume 1, Wiley, Hoboken.
- 14. Johnson, N.L. and Kotz, S. (1970) Continuous Univariate Distributions-1. Chapter 13, John Wiley & Sons, Hoboken.
- 15. Barr, D.R. and Sherrill, E.T. (1999) Mean and Variance of Truncated Normal Distributions. The American Statistician, 53, 357-361. https://doi.org/10.1080/00031305.1999.10474490
- 16. Burkardt, J. (2014) The Truncated Normal Distribution. Department of Scientific Computing Website, Florida State University, Tallahassee.
- 17. Lecture 4: Selection. web.ist.utl.pt. Instituto Superior Técnico, Lisboa, 1.
- 18. Fernandez-de-Cossio-Diaz, J. (2017) Truncated Normal: Compute Mean and Variance of the Univariate Truncated Normal Distribution (Works Far from the Peak).
- 19. Johnson, N.L., Kotz, S. and Balakrishnan, N. (1995) Beta Distributions. In: Johnson, N.L., Kotz, S. and Balakrishnan, N., Eds., Continuous Univariate Distributions, Vol. 2, 2nd Edition, Wiley, Hoboken, Chapter 21.
- 20. Jones, M.C. (2009) Kumaraswamy’s Distribution: A Beta-Type Distribution with Some Tractability Advantages. Statistical Methodology, 6, 70-81. https://doi.org/10.1016/j.stamet.2008.04.001
- 21. Forge Project. Handbook on Probability Distributions: R-Forge Distributions Core Team University Year 2009-2010. https://r-forge.r-project.org/scm/viewvc.php/*checkout*/pkg/inst/doc/probdistrmain.pdf?revision=24&root=distributions&pathrev=24
- 22. Another Useful Handbook of Probability Distribution Functions. http://www.stat.rice.edu/~dobelman/textfiles/DistributionsHandbook.pdf
- 23. Corda, C., Hendi, S.H., Katebi, R. and Schmidt, N.O. (2013) Bohr-Like Black Holes. Journal of High Energy Physics, 06, 008.
- 24. Corda, C. (2015) Time Dependent Schroedinger Equation for Black Hole Evaporation; No Information Loss. Annals of Physics, 353, 71-82. https://doi.org/10.1016/j.aop.2014.11.002
- 25. Corda, C. (2015) Precise Model of Hawking Radiation from the Tunneling Mechanism. Classical and Quantum Gravity, 32, Article ID: 195007. https://doi.org/10.1088/0264-9381/32/19/195007
- 26. Calmet, X., Carr, B. and Winstanley, E. (2014) Quantum Black Holes. Springer Briefs in Physics, Springer Verlag, New York.

Appendix

Here, we present certain formulae, as it relates to the truncated Gaussian distribution, and identities relating to the calculation of the entropy. These are given for reference purposes. The symbols are explained in the text.

The truncated Gaussian distribution function is given by,

${f}_{r}=f\left(r,0,\sigma ;0,R\right)=\sqrt{\left(2/\text{\pi}\right)}1/\sigma \text{exp}\left(-{r}^{2}/\left(2{\sigma}^{2}\right)\right)1/\text{erf}\left(R/\left(\sqrt{2}\sigma \right)\right)$ (3-1)

The corresponding cumulative distribution function reads,

${F}_{r}^{CDF}=F\left(r,0,\sigma ;0,R\right)=\text{erf}\left(r/\sqrt{2}\sigma \right)/\text{erf}\left(R/\sqrt{2}\sigma \right)$ (3-2)

We have assumed that the radius has support, $r\in \left[a,b\right]=\left[0,R\right]$, and that the mode equals zero. Thus, $\mu =0$. The shape parameter, $\sigma $, is determined by the inside surface radiative force, ${F}_{R}={F}_{R}^{\left(4\right)}$, within the 4-D black hole. The function, $\text{erf}\left(x\right)$, is the error function.

$\sqrt{\left(\text{\pi}/2\right)}\text{\hspace{0.05em}}\sigma {\text{e}}^{{R}^{2}/\left(2{\sigma}^{2}\right)}\text{erf}\left(r/\sqrt{2}\sigma \right)={\text{e}}^{\left({R}^{2}-{r}^{2}\right)/\left(2{\sigma}^{2}\right)}{F}_{r}^{CDF}/{f}_{r}$ (A-1)

${F}_{R}^{\left(4\right)}/{f}_{R}={p}_{1}^{\left(4\right)}{A}_{R}^{\left(4\right)}/\left({\rho}_{1}^{\left(4\right)}{A}_{R}^{\left(4\right)}\right){M}_{R}=0.1{M}_{R}{c}^{2}$ (A-2)

${F}_{r}^{\left(4\right)}/{f}_{r}={p}_{r}^{\left(4\right)}{A}_{r}^{\left(4\right)}/\left({\rho}_{r}^{\left(4\right)}{A}_{r}^{\left(4\right)}\right){M}_{R}=0.1{e}_{r}/{e}_{r}{M}_{R}{c}^{2}=0.1{M}_{R}{c}^{2}$ (A-3)

${\int}_{0}^{r}{\text{e}}^{\left({R}^{2}-{r}^{2}\right)/\left(2{\sigma}^{2}\right)}\text{d}r}={\displaystyle {\int}_{0}^{r}{f}_{r}/{f}_{R}\text{d}r}=1/{f}_{R}{\displaystyle {\int}_{0}^{r}{f}_{r}\text{d}r}={F}_{r}^{CDF}/{f}_{R$ (A-4)

${F}_{r}^{CDF}/{f}_{R}={M}_{R}{F}_{r}^{CDF}/\left({M}_{R}{f}_{R}\right)={M}_{r}/\left({\rho}_{R}{A}_{R}\right)={M}_{r}{c}^{2}/\left({e}_{R}{A}_{R}\right)=0.1{M}_{r}{c}^{2}/{F}_{R}$ (A-5)

Some identities relating to the calculation of entropy follow. Let,

$\gamma \left(s,x\right)=\text{lower incomplete gamma function}={\displaystyle {\int}_{0}^{x}{t}^{s-1}{\text{e}}^{-t}\text{d}t}$ (A-6)

$\Gamma \left(s,x\right)=\text{upper incomplete gamma function}={\displaystyle {\int}_{x}^{\infty}{t}^{s-1}{\text{e}}^{-t}\text{d}t}$ (A-7)

Then,

$\gamma \left(s,x\right)+\Gamma \left(s,x\right)=\Gamma \left(s\right)$ (A-8)

where, $\Gamma \left(s\right)$, is the gamma function. From this, it follows that

$\gamma \left(s,x\right)=\Gamma \left(s\right)-\Gamma \left(s,x\right)>0$ (A-9)

This inequality will always lead to a positive value for entropy.

$\begin{array}{c}{\displaystyle {\int}_{0}^{r}{\text{e}}^{0.4\left({R}^{2}-{r}^{2}\right)/{\sigma}^{2}}{r}^{3/5}\text{d}r}=-\frac{{5}^{\frac{4}{5}}}{{2}^{\frac{9}{5}}}\left[\Gamma \left(\frac{4}{5},\frac{2{r}^{2}}{5{\sigma}^{2}}\right)-\Gamma \left(\frac{4}{5}\right)\right]{\sigma}^{8/5}{\text{e}}^{0.4{R}^{2}/{\sigma}^{2}}\\ =+\frac{{5}^{\frac{4}{5}}}{{2}^{\frac{9}{5}}}\gamma \left(\frac{4}{5},\frac{2{r}^{2}}{5{\sigma}^{2}}\right){\sigma}^{8/5}{\text{e}}^{0.4{R}^{2}/{\sigma}^{2}}\end{array}$ (A-10)

The symbol, r, refers to the 4-D radius within the black hole. Other identities are derived in the text.