Journal of Modern Physics
Vol.3 No.9A(2012), Article ID:23081,9 pages DOI:10.4236/jmp.2012.329154

A Guess Model of Black Holes and the Evolution of Universe

Zhenhua Mei1, Shuyu Mei2

1College of Chemistry and Molecular Engineering, Qingdao University of Science and Technology, Qingdao, China

2Medical College, Qingdao University, Qingdao, China


Received June 10, 2012; revised July 16, 2012; accepted July 23, 2012

Keywords: Black Hole; Evolution; Universe; Gravitation; Entropy


Based on the gravitational theory, fundamental data, and comprehensible suppositions, an evolution model of the universe was proposed. The universe exists in explosion and constringency mobile equilibrium state. The critical sizes of celestial bodies were calculated in their evolution process.

1. Introduction

As the theoretical ratiocination (Einstein’s theory of relativity) and the cumulating of the more and more obtained data or phenomena from astronomic observation of the celestial bodies, people began to have a general impression of the universe, such as the clangorous words expressed: big bang, expanding universe, and black hole et al. (Ginsburg 1985). From a point of view that the microcosm decides the macrocosm, more information of the universe can be predicted from present knowledge of the elementary particles.

In this paper, some parameters of evolution of the universe were calculated. The main thoughts were the Newton’s universal gravitation and the structures of matter.

2. Preparation of Research

2.1. Fundamental Data and Suppositions

The gravitational constant, G, is 6.670 × 10–11 N m2·kg−2. The radius of hydrogen atom, a0, is 5.29 × 10−11 m. The mass density of atomic nucleus, ρn, was determined (Erdei 1976) to be 2 ×1017 kg·m−3. The mass of neutron, mn, is 1.67495 × 10−27 kg, radius rn, 1.2598 × 10−15 m. The combining energy of nucleus is 8 MeV. The gravitational acceleration on earth surface is 9.81 m·s−2. The sun has the mass of 1.989 × 1030 kg, density of 1.409 × 103 kg·m−3, radius of 6.960 × 108 m.

Suppositions: the combining energy of nucleus (strong interaction in nucleus) came from gravitation, and the nucleus consisted of gravitons; the combinative form between nucleus was its overlaps one another.

2.2. Formula Derivation A—Gravitational Acceleration of a Particle from Solid Sphere

According to Newton’s law, the gravitational force between two distant particles has the form:

, and. Hence, the gravitational acceleration,.

However, to a neighboring large body (has a mass of M, radius of R), the radius of the body cannot be neglected; its accurate result can be obtained by processing a mathematical integral.

As Figure 1 expressed, a solid sphere has a radius of R, a homogeneous density of ρ, the particle has a mass, m, and has a distance, nR, from the center, the gravitational acceleration, a, of the particle was derived as follows:

Figure 1. Sketch map of a particle (m), which with a distance, nR, from the center of the solid sphere.


where, “M” is the arbitrary point in the solid sphere, “h” the point corresponding to “M” horizontally. “r” is the distance of the point “M” from original point “o”, ρ′ is ideally homogeneous mass density of a solid sphere, ρ the average one. k stands for the calibration factor of average mass density.

For a general integral


For integral


For general integral

Let f =






Substituted (3) with (4)


As the results of (2) and (5) expressed, thus, (1) becomes

Let, then,


F(n) is a function of n. The, F – n, diagram was showed in Figure 2.

For the earth, the average mass density ρ is 5.5153 × 103 kg·m−3, the radius R is 6.356078 × 106 m. The gravitational acceleration on earth surface was calculated (n = 1) to be 12.68 m·s−2. However, the real value of gravitational acceleration on earth is 9.81 m·s−2. It is 0.7736 times than the calculated one. It is because of the earth that does not have an ideally homogeneous density used in above formula derivation process. The ratio can be used as calibration factor, k, of average density of a solid sphere (ρ′ = kρ. For earth, k = 0.7736; and it was supposed to be fit for any other celestial bodies).

2.3. Formula Derivation B—Self-Gravitational Pressure in Center of Solid Sphere

The gravitational pressure in center from solid sphere can be calculated according to Equation (6).

, (where,)


The integral, , can be calculated diagrammatically. The diagram showed in Figure 3.

The area beneath the curve calculated by computer gave the result of 0.62129. Thus, p = 0.62129 × 4Gkρ2R2/3. Supposing the considered celestial body has the same value of calibration factor k of earth, then we obtained

p = 0.6408Gρ2R2.(7.1)

2.4. Formula Derivation C—Gravitational Potential Energy from Solid Sphere

The gravitational potential energy, Eg, from a solid sphere can be calculated.


According to relation (6), the integral, , can be calculated diagrammatically. The diagram showed in Figure 4.

The area beneath the curve calculated by computer

Figure 2. Diagram of function F (n).

Figure 3. Partial diagram of function F (n). The variable n changes from zero to one.

Figure 4. Partial diagram of function F(n). The variable n changes from 1 to ∞.

gave the result of 1.7193. Therefore, we have

Eg = 1.7193 × 4GkρmR2 = 5.320GρmR2.(8.1)

3. Target Calculation

3.1. Maximal Radius of Celestial Bodies

Atoms have similar energy demand when one electron is compressed in nucleus forming a new atom (which has lower atomic number). Taking atom zinc as an example, when one electron of Zn is compressed in nucleus, the isotopic atom Cu is then formed. The energy demand, ΔEa, can be calculated. According to quantum theory, the energy of an outer electron of a atom has the form E ≈ −13.6Z2/n2 (in ground state). Thus,




For an extreme hot electron of atom of hydrogenΔEa = ΔE1 + 13.6 = 63.6 eV = 1.02 × 10−17 J As the gravitational force increasing, many radioactive isotopes formed with the atomic number decreasing. Considering the atom of hydrogen, that has the maximum compression radius a0, we guessed that the hydrogen atoms would be the last element to be compressed to nucleus. Then, the maximum pressure, pmax,a, for compressing a hydrogen atom to neutron could be derived as:

ΔEa = FΔr = pSΔr = pSa0.


where, Δr was replaced with the radius of hydrogen atom a0.

The pressure comes from the gravitation of celestial body. Combining Equations (7.1) and (9), we have 0.6408Gρ2R2 = 1.64 × 1013then,


If a solid celestial body has the same average mass density of earth (5.5153), then

To a hot gas celestial body, which has the same average mass density of sun (1.409), then, we have


3.2. Minimum Neutron Star

A model of neutron star was showed in Figure 5. A neutron body with the radius of R is in the center of the star. A general matter layer with the thickness of (n − 1)R is around the neutron body.

The gravitational pressure on the surface of neutron body can be calculated. It contains two parts.

Equation (6) can express linearly in a short segment.

F(n) = 1.355571n, (n = 0 → 1).

According to Equation (7),

where, ρn is the nuclear mass density, ρ the mass density layer matter. Substituting the parameters, we have


Figure 5. A model of neutron star.

Combining Equations (9) and (10), we have


According to Equations (6) and (11), we obtained following results:

R = 11.2 m, nR = 1121 m R = 16 m, nR = 32 m R = 261 m, nR = 261.6 m R = 15257 m, nR = 15257 m.

Combining Equations (7.1) and (13), we have

3.3. Minimal Mass of Black Hole

When the radius of a neutron star achieved its maximal value, as the mass continuing accumulating, a graviton body began growing, the radius of neutron star then decreased correspondingly. When the mass of graviton body grew large enough, the neutron star began to have the ability to draw back a photon; we call it becoming a body of black hole. The maximum radius of black hole can be calculated by imitating the process of what done in paragraph 3.2. Just replacing ρ and ρn with the value of ρn and ρg (showed in Equation (20)), Equation (11) became Equation (12). Where ρg is stands for the mass density of graviton.


According to Equations (12) and (13), following relations were obtained:

R = 0.971 m, nR = 9714 m, M = 7.679 ×1029 kg

R = 9.714 m, nR = 9714 m, M = 7.679 ×1029 kg

R = 97.10 m, nR = 9710 m, M = 7.671 ×1029 kg

R = 1154 m, nR = 9279 m, M = 6.722 ×1029 kg

R = 4078 m, nR = 7340 m, M = 4.537 ×1029 kg

R = 10793 m, nR = 11873 m, M = 3.672 ×1030 kg

R = 14090 m, nR = 14908 m, M = 7.825 ×1030 kg

R = 42004 m, nR = 42273 m, M = 1.971 ×1032 kg

R = 74934 m, nR = 75084 m, M = 1.114×1033 kg.

3.4. Calculation of Graviton

Nuclear has the combining energy of about 8 MeV (Wichmann 1971). Deducting the influence of proton, the maximum combining energy of nuclear would be 9.15 MeV, i.e., 1.47 × 10−12 J. The parameter of the radius of a graviton, r0, is the first needs to obtain. Depending on above suppositions described in Section 2.1, and the normal gravitational potential energy formula, E = Gm1m2/r, substituting m1, m2 with the mass of neutron, E with the combining energy in nucleus, then we have,

Supposing the energy demand for destroy a neutron, ΔEn, is equal to its maximum combining energy, then we have

ΔEn = 1.47 ×10−12 J.

Following the derivate process in Section 3.1, but replaced ΔEa with ΔEn, a0 with rn, Equation (9) becomes


Following the derivate process in Section 3.2, but replaced ρ and ρn with ρn, and ρg, pmax, a with pmax, n, Equation (11) became Equation (14).



A Smallest Known Black Hole (named XTE J1650- 500) has the mass, Mb, of 7.558 × 1030 kg (Robert & Rob 2008). Perhaps it is not the minimal black hole theoretically; however, it was considered to be true here temporarily. According to the data and Equation (8), we have



Combining Equations (15), (16) and (17), we have Equations (18) and (19).



Combining Equations (18) and (19) diagrammatically, we obtained n = 1.0190.

R = 14090 m.

nR = 14358 m.

ρg = 6.3096 × 1017 kg·m−3.(20)


where mg is stands for the mass of a graviton.

4. Evolvement of the Universe

The evolutive process of the universe was divided by us into six parts; and the process goes a circle (a → f →a).

a) Growing of normal celestial body. Preponderant normal celestial bodies will attract substances from its surroundings. It gradually grows as it accumulates the mass. It would have the maximum mass of 5.06 × 1029 kg, radius, 4.41 × 108 km.

b) Formation of neutron star. As the continuously accumulating of matter, the normal molecular celestial body begins to compress accompanied with neutron body engendered in the center. The compressing process looks like a process of phase variation. In the process, it will release large quantity of energy (lose of masses), become an extremely hot body and shine. A maximum neutron star has a maximum mass of 2.98 ×1030 kg, maximum radius of 15.3 km.

c) Growing of graviton star. With continuously compressing, a graviton body would begin to grow in the center. The graviton star has a critical mass of 7.558 × 1030 kg, radius of 14.9 km.

d) Formation of black hole. When graviton star grows to its critical value, it becomes a minimum black hole in the mean time. It can capture photons just near its surface. Afterward, it gradually grows in mass with the accumulating of matter. As the mass increasing, far more distant photons can be hold in. The compressing energy that released absorbed by itself. As an isolated system, the entropy of the black hole decreased spontaneously in its growing process.

e) Explosion of black hole. When the density of black hole increased large enough, the distance between gravitons going over a critical small value, the gravitational mechanism destroyed. The gravitation becomes repulsive force. Then, explosion occurred; the entropy increased sharply. Neutrons, protons, electrons, etc. particles and afterward atoms were produced subsequently in the explosion process. After explosion, this part of universe exists in expansion state in a unabiding period.

f) Constringency of partial universe. In the late period of explosion, the running matters slowed down, spread in wide space, mixed, and became part of cosmic dusts. Then, partial universe would exist in a contractive state in a long run. In intermediate constringency, the celestial bodies formed.

5. Discussion

When a celestial body has abundant material source in its surroundings, it will grow quickly and has a larger mass than we calculated. If the celestial body exists in process b, it would be a fixed star. The larger mass it has, the more quickly it will be compressed, and more efficiency the energy would release. i.e., the larger mass it owns, the higher brightness it will has. In this way, a shining star could burn with no existence of special fuels of nuclear fusion, but merely existence of normal atoms.

According to the calculated results in this paper, it is evident that every galaxy will have a big black hole in its center.

The universe consists of vast cosmic dusts and galaxies. They go along circle evolutive processes described above. Occasionally a limit big black hole explodes, causing the expansion in a partial space. However, in a long time scale, most space contracted slowly. The universe exists in a mobile equilibrium state of explosion and constringency (Thomas & Hermann 1948). Those observed stars, which leave away at acceleratory velocity (Hubble 1929), might result in other reasons. It is by no means that the whole universe is expanding now. Reversely, our surrounding universe showed us its contractive views. We cannot imagine that the shape of our Milky Galaxy is formed in an expansion process.

6. Tags

A roughly information about graviton was derived in this paper. It is an exiting result to the researchers in fundamental physics. Future observation of existence of new smaller black hole is apparently meaningful.


  1. V. L. Ginsburg, “On Physics and Astrophysics,” Nauka Press, 1985.
  2. G. Erdei, “Principles of Structure of Matter,” Mir Press, 1971.
  3. E. H. Wichmann, “Quantum Physics,” Berkeley Physics Course, McGraw Hill, 1971.
  4. N. Robert and G. Rob, “NASA Scientists Identify Smallest Known Black Hole,” NASA’s Goddard Space Flight Center, 2008.
  5. G. Thomas and B. Hermann, “The Steady-State Theory of the Expanding Universe,” Royal Astronomical Society, Vol. 108, 1948, p. 252.
  6. E. Hubble, “A Relation between Distance and Radial Velocity among Extra-Galactic Nebulae,” Proceedings of the National Academy of Sciences, Vol. 15, 1929, p. 16.