A neutron star with mass greater that the Tolman-Oppenheimer-Volkoff (TOV) limit MTOV~3 solar mass MSUN becomes a black hole in which the star core collapses to a gravitational mass singularity in general relativity [ p172], [ black hole/§2.3]. The quark structure of the neutrons has been neglected. But in these dense stellar objects, gravitational forces are strong and their gradients large so that such a neglect cannot be carried out without appropriate justification. Further, such a singularity contradicts quantum mechanics.
Dark matter needed to account for many observation, e. g., the galaxy rotation curve, and the dark energy required to drive the rapidly expanding universe remain hypothetical as they cannot be observed.
Recently, it was shown  that the “hidden”, unobservable relative energy between the diquark and quark in nucleons can interact with their ambient gravitational potential and possibly play the roles of dark matter and dark energy in the universe. The underlying theory is the scalar strong interaction hadron theory SSI  .
In this paper, it is shown that, as a neutron on the surface of a neutron star with mass ~MTOV falls toward the star center, not only gravitational energy is released, but a negative relative energy also emerges simultaneously. These two energies can cancel each other and prevent the collapse of the star and formation of a mass singularity. The negative relative energy generated in an expanding galaxy can play the role of dark matter and account for the galaxy rotation curve. The positive relative energy generated in outer regions of the universe can play the role of dark energy, leads to accelerating expansion of the universe and may eventually give rise to cosmic rays.
The relevant parts of SSI are outlined in Sections 2-5. Creation of negative relative energy applied to neutron star collapse is considered in Sections 6-7. In Section 8, such negative relative energy reinforces the existing gravitational potential in an expanding galaxy to account for the galaxy rotation curve qualitatively. The accelerating expansion of the universe by means of the positive relative energy created in the outer parts of the universe is qualitatively outlined in Section 9. Some scenarios of the outer regions of the universe including cosmic ray generation are shown in Section 10. The so-obtained scenario of the universe is summarized in Section 11.
2. Outline of Construction of Baryon Wave Equations in SSI  
The starting point is the Dirac equations for three quarks A, B and C with masses mA, mB and mC located at xI, xII and xIII, respectively. The quarks interact with each other via scalar strong potentials VAB(xI), VAC(xI), in [ (2.1-3) or   (9.1.1-3)]. These equations are in van der Waerden’s spinor form which is manifestly Lorentz covariant contrary to the conventional bispinor Dirac equation which is not. For quark B, for instance, the equations read
where , . The spinor indices run from 1 to 2 and the quark spinors χB and ψB are linear combinations of the four components in the conventional Dirac bispinor ψ according to [  (C11)]. mB is quark mass, about 2/3 of the proton mass [  Table 5.2]. The scalar strong potentials VBA and VBC generated by quark A and C, respectively, are governed by [ (2.4) or   (9.1.4)]. There is no gauge fields for electromagnetic interactions between the quarks. Here, the scalar gravitational potential VGB from ambient matter and energy, neglected in the references, has been included.
To construct baryon wave functions, the three pairs of quark equations, one of them being (2.1), together with the cited strong potential equations, are multiplied together and the products of the quark spinors and those of the strong potentials are generalized to baryon wave functions and baryon potentials nonseparable in xI, xII and xIII according to [  Sec. 9.2]. For ground state baryons, quark C at xIII has been merged into quark A at xI to form a diquark at xI. The nearly intractable three body problem is thus reduced to a manageable two body problem.
In [  Sec. 9.3], each quark mass above has been generalized to a quark mass operator operatng on a flavour function in a fictitious, complex internal or flavour space zp, where p runs from 1 to 6 representing the six quark flavours. The quark mass is recovered as eigenvalue of such an operator. The three flavour functions, one for each quark, are analogously multiplied together and generalized to a nonseparable baryon flavour function. The product of the three mass operators are similarly generalized to a nonseparable baryon mass operator in the flavour space and leads to a mass Mb [  (9.3.19)] representing the quark mass contribution to the baryon mass.
The baryon flavour function can now be removed leaving behind the ground state baryon wave equation [ (2.9)], [  (9.3.16)]
The interbaryon strong potential Φb is a triple product of three strong interaction quark potentials of the type of VBA and VBC in (2.1) and is governed by [  (9.2.13b)]. The six component χ and ψ can be decomposed into a quartet or spin 3/2 part, which can be separated off, leaving behind a doublet or spin 1/2 part [  (10.0.6)]:
where χ0 and ψ0 the wave functions of the doublet baryons.
3. Laboratory and Relative Spaces
Since quarks cannot be observed, their coordinate spaces are converted into an observable laboratory space Xμ for the baryon and a relative space x between the diquark and the quark via the linear transformation given above [ (6.2)] or by [ (5.1)], [  (3.1.3a)]
For observable particles, am is often determined by that Xμ is the center of mass of these particles. If these particles have equal mass, am = 1/2. Such kind of determination cannot be done here because quarks are not observable individually. am has been taken to be an arbitrary real constant and represents a new degree of freedom. The relative space is “hidden” [  (3.1.3)] and cannot be observed. If it were observable, then (3.1), with a given am, leads to that both the diquark at xI and quark at xII become observables, contrary to experience. This can also be seen directly in the first of (3.1) in which the right side members cannot be measured; hence xμ is also “hidden”, independent of am. It turns out that the bulk of hadron physics in SSI lies in such “hidden” spaces.
The baryon wave functions in (2.1) have been factorized into the form of [  (10.1.1)],
[  (3.1.6)] = (3.3)
where EK is the energy of the baryon and K its momentum. x0 is the relative time and −ωK the associated relative energy in the “hidden” relative space and can also not be observed within SSI.
In spherical coordinates, [  (3.1.7b)], the doublet wave functions in (3.2) with total angular momentum j = 1/2 and orbital angular momentum l = 0 read [  (10.2.3)]
where the Y’s are the usual spherical harmonics. is found by changing the signs of f0(r) in (3.4).
4. Radial Wave Equations in Relative Space, Solutions and Results
Consider baryons at rest, K = 0. am = 1/2 is set as in the meson case [ (3.5.7)], [ (5.7.2)]. Similarly, the “hidden” relative energy −ω0 = 0 is also set following the meson case [ (3.5.6)], [ (5.7.1)]. Insertion of (3.1-4) into (2.3) using [  (3.1.4)] leads to the radial wave equations [ (6.9)], [  (10.2.12)]
where the subscript d denotes doublet and [  (10.2.2a)] gives the diquark-quark strong interaction potential
here, the nonlinear potential Φcd(x) [  (10.2.2)] vanishes for large normalization volume Ωcb®∞ in [  (10.3.14)]. The db’s are unknown integration constants.
The two coupled third order Equations (4.1) have been converted into six first order equations [ (10.7.5)], [ (10.4.5)] which have six linearly independent solutions each associated with its own λ+ values in [  (10.2 8a)] found from the sixth order indicial equation [ (10.4.4)]. These eventually led to two free parameters w(1) and w(2) in [  (11.1.1)]. It turned out that only db4 = 0 led to results in near agreement with data [ §11.1.3]. Thus, there are four db and two w totalling six constants that can be varied to satisfy the six boundary conditions posed by the six first order system.
Due to the large number of unknown constants, (4.1-2) could not be solved as a conventional eigenvalue problem. A less ambitious approach has been adopted. The known mass of the neutron is used as input for the eigenvalue E0 and the quark masses obtained from meson spectra given in [  Table 5.2] are used as input for Mb according to (2.2c) where mA = mC is the d quarks mass and mB the u quark mass. The six first order equations [ (10.7.5)], [ (10.4.5)] have been solved on a computer. It turned out that given a db2 value, the remaining five constants can be freely chosen such that g0(r) and f0(r) converge at large r. This is allowed since this takes place in the “hidden” relative space.
Here, the word “finite” on the first line of the second paragraph on [ p239] turned out to be incorrect in subsequent numerical calculations and should be changed to “infinite”; db2 can be varied continuously to yield converging g0(r) and f0(r). One so-obtained neutron wave functions are plotted in [ Figure 12.1], [ Figure 11.1b] reproduced in Figure 1.
These wave functions have led to the nearly correct predictions of the neutron life and the electron asymmetry parameter A or the neutrino asymmetry parameter B [ Table 12.1] in its beta decay.
5. Relative Energy and Gravitation
If am = 1/2 and −ω0 = 0 used in (4.1) were not assumed, a derivative in (2.3) operated on (3.2) will, in terms of (3.1) [see   (3.1.4)], be changed as follows
The expressions on both sides of the arrow will be equal if [ (6.6) or   (3.1.10a)]
which keeps (4.1) invariant. The relative energy −ω0 is “hidden” with respect to strong and electromagnetic but not gravitational interactions . This is indicated by VBG in (2.1), which has been neglected in SSI because it is generally small next to the strong VBA, VBC there. In neutron stars and black holes, however, VBG can be large and may not be neglected. Since am, hence also ω0, can take on any real value, from −∞ to +∞, the relative energy can be important is astrophysics.
Figure 1. g00(r) and f00(r) are g0(r) and f0(r) in (4.1) normalized according to [  (10.3.14)] for the db2 = −0.3202 case in [  Table 11.1] where the four db and two w constants are given. r is the diquark-quark distance. This wave function has been used to evaluate neutron decay parameters in [ §12.6.3].
In , −ω0 > 0 is identified as dark energy and −ω0 < 0 plays the role of dark matter. Both E0 and ω0 are observables with respect to gravitational forces. Lorentz invariance leads to that the associated baryon momentum K in (3.3) and the relative momentum represented by in (5.1) are also observables. Thus, the conjugate variables X, the coordinate of the baryon, and x, the relative coordinate in (3.1), are also observables, remembering that x is still a “hidden” variable with respect to strong forces. By (3.1), both the diquark coordinate xI and quark coordinate xII are now observables with respect to gravitation indicated by VBG in (2.1).
6. Neutron Star Collapse Scenario
Consider the following idealized scenario. A neutron star with a mass MNS = MNSa, a critical mass, such that the quantum degeneracy pressure of the neutrons and the strong neutron-neutron repulsion precisely balance off the gravitational pressure. Let an external neutron arrive at the surface of this star. The star now becomes too heavy and a gravitational collapse is anticipated to start. Conventionally, this collapse gives rise to a black hole and continues until the star ends up in a mass singularity at its center; the quark structure of the neutron is ignored. Actually, in neutron stars and black holes, the gravitational forces are strong and their gradients large that the diquark-quark structure of the neutron may need be taken into account. This scenario is illustrated in Figure 2.
In Figure 2, an external neutron arrives on the surface of the neutron star at position a on the right part of the figure. The interquark distance is fm according to Figure 1. Xa lies in the middle corresponding to am = 1/2 in (3.1) and (4.1) and hence ω0 = 0 by (5.2). This neutron falls from Xa to Xb. To begin with, the three positions xI,a, Xa and xII,a are simply
Figure 2. Illustration of the above idealized scenario. The horizontal line represents the radius R inside the neutron star. The center of the star is at R = 0. An external neutron arrives at Xa on the surface of this star at R = Ra = xI,a in position a and falls to Xb closer to the star center in position b. xI,a and xII,a are the positions of the diquark and quark of the arriving neutron, respectively. When this neutron has moved to Xb, the diquark and quark positions become at first xI,b and xII,b, respectively; the diquark-quark distance is the same in both positions. Due to gravitational collapse, however, the neutron star becomes smaller at position Xb so that the diquark-quark distance in the arriving neutron is also reduced via the replacement of xII,b by xII,bG. Rb = xI,b is radius of the collapsed star when the neutron has moved to Xb.
shifted to the left to become xI,b, Xb and xII,b, respectively. The diquark-quark distance remains unchanged. In the same fall, however, the neutron star radius at position b is reduced to Rb < Rq due to gravitational collapse. The incoming neutron is thus similarly compressed to become smaller by a factor of Rb/Ra so that rb is also reduced to a smaller value rbG by the same factor;
As was mentioned at the end of Section 5, gravity interacts directly with the quark at xII in (2.1), hence also with that in Figure 2. Phenomenologically, it also interacts with the diquark at xI in Figure 2. An equation of motion for the diquark on par with (2.1) for the quark has not been constructed. In an eventual such equation, the diquark mass is expected to be about twice the mass of the quark in (2.1) or at least greater than it. Assuming this, gravity will pull the diquark closer to the star center relative to the quark, as is indicated in Figure 2. Note that the quark and the diquark are point particles in quantum theory and hence cannot be “spaghettified” here. The radius Rb of the star at position b is taken to be xI,b
If Xb were taken to be the star radius Rb at position b, then the diquark coordinate xI,b in Figure 2 would be negative at the star center Xb = 0, contrary to the requirement that a radius ≥ 0. Practically, Xb and Rb are of the magnitude km while their difference is of magnitude fm; both are nearly the same.
The radial wave Equations (4.1-2) hold for the arriving neutron in position a. At position b, (6.1) shows that the diquark-quark distance or radius r in (4.1) is reduced by a factor of α. Equations (4.1-2) remain invariant under the transformations , , , , . The neutron mass E0 and the quark masses in Mb at position b are unchanged. So does also the confining constant db2 = −0.3202 in Figure 1; this value is tied up by neutron decay [ §12.6.3]. The remaining three db and two w constants mentioned above Figure 1 can be varied to produce converging g0(αr) and f0(αr) analogous to the computations that led to Figure 1. These wave functions will obviously be much narrower than those in Figure 1 having an average width of αra~4α fm.
7. Removal of Mass Singularity and “Black” Neutron Star
Application of (3.1), (5.2) and (6.1) to position b in Figure 1 yields
where −ω0 is the relative energy gained when the external neutron arriving at position a falls to position b. In this fall, the gravitational energy EG released is given by
where G is the gravitational constant and RSa is the Schwarzschild radius of this collapsing star with radius Ra in Figure 1 These two energies cancel each other when the star’s radius
i.e., when the neutron star in Figure 1 with the critical mass MNSa has a radius Ra equal to its Schwarzschild radius RSa. At the star center Rb = 0, both energies become infinite but these both singularities cancel; no net mass singularity is created when (7.3) holds. This is consistent with that SSI is a quantum theory not allowing such a singularity. Note that when approaching the star center and Rb becomes of the size of the neutron of a few fm. For Ra ~10 km; am = Ra/2Rb is a huge number ~1018.
In this case, the fall of this external neutron actually gains no energy because EG − ω0 = 0 for all Rb. This neutron therefore remains at position a and becomes “weightless”. The anticipated collapse of the neutron star does not start and no mass singularity is created.
The radius Ra of the initial neutron star has been estimated by equating gravitational pressure to the pressure of the degenerate neutrons [ p161]
which together with (7.2-3) yields Ra = RSa = 10.8 km and MNSa = 3.6MSUN. This mass is consistent with the similarly estimated maximum mass 4.3MSUN [ p162] of a neutron star. It is of somewhat higher than MTOV~3MSUN. cited in Section 1. However, estimates of MTOV are uncertain and vary widely due to that the equations of state for extremely dense matter are not well known and to uncertain effects of neutron-neutron repulsion and of pressure in general relativity [ p162]. Thus, the so-estimated star mass may be used on equal footing with those in the literature.
A possible scenario is as follows. As additional neutrons arrive at the star surface, a thin shell of weightless neutrons is added to the star at first. Let the mass of this shell be ΔMNSa, the new Schwarzschild radius will be and the new star radius will be , assuming that the neutron density in the shell is the same as that in the star.
Since RSaΔ > RaΔ, a black hole is created. As a new neutron arrive at the new star surface with radius RaΔ, it will tend to fall inwards. This time, the gravitational energy (7.2) will slightly exceed the relative energy loss (7.1) and this new neutron tends to fall slowly. However, the pressure of the degenerate neutron gas in the shell is now unopposed by the weightless neutrons and this gas will therefore expand accompanied by a reduced density. If the expansion reaches a radius > RSaΔ, the back hole is lost and this enlarged star will be a visible neutron star similar to the initial star in Figure 2. If the expansion stops just inside RSaΔ, the created black hole persists and one arrives at a situation nearly the same as that on the right part of Figure 2 with Ra replaced by RaΔ and (7.3) holds. This process may be repeated and more shells are added. It can in principle lead to a super-massive black hole. How the neutrons, part of them “weightless”, are distributed inside such a black hole is unknown and needs not be known, inasmuch as they are inside its Schwarzschild radius and cannot be observed any way; there is no conflict with observation.
In this scenario, the neutrons in such a black hole fill it up to its Schwarzschild radius. Such a black hole may more suitably be called a “black” neutron star. The absence of mass singularity here is consistent with the conjecture that an eventual future quantum gravity theory will not contain any singularity.
8. Galaxy Rotation Curve without Dark Matter
There is a large body of data that require the presence of dark matter. The first one is the galaxy rotation curve [ Dark Matter §3.1]. However, nothing is known about this hypothetical dark matter. In , this matter has been identified as the negative relative energy −ω0 < 0 between the diquark and quark in nucleon. Here, the role of the dark matter in accounting for this rotation curve will be played by this −ω0.
Consider again an idealized scenario as follows. A spiral galaxy with an appreciable part of its mass consisting of hydrogen gas was in its earlier stage of development. In that stage, this galaxy was smaller, denser and hotter according to the big bang model and the gravitational potential resisting such an expansion was insufficient. It therefore expanded. The situation is analogous to a violation of Jeans criterion for a gas cloud.
To illustrate the mechanism, let us turn off the gravity for a moment. Follow now the movement of a proton, denoted by c, in a hydrogen molecule. In this thermal expansion, this molecule will collide with other molecules, exchange energy and momentum with them and end up in a new position, labeled d, farther out from the galaxy center. In this journey, the forces involved are all Coulomb forces, between the proton and the orbiting electron and between the orbiting electrons in other molecules. The protons get dragged along; their quark structure is not involved in the expansion.
The situation is illustrated in Figure 3.
In Figure 3, the coordinate of proton c is Xc which lies in the middle between the diquark at x1,c and quark at xI1,c. This corresponds to am = 1/2 and ω0 = 0 in (7.1) for a free baryon at rest, similar to those for the neutron in position a in Figure 2. In this thermal expansion, these three coordinates are simply shifted to the right with . The new xI,c and xII,c are on both sides of Xd with the same separation distance.
Turn now on gravity, the situation is reversed. Gravitational effects on electrons are negligible due to their small mass. Gravitational forces now act on this proton at Xd and try to pull it back together with other protons acted upon. However, they turn out to be too small to account for the galaxy rotation curve and prevent the escape of the outer stars. Here, such escapes are prevented by including the negative energies generated by differentiated gravitational pull on quarks and proton.
Just like the neutron star case mentioned below (6.1), gravity also acts directly on the quarks and tends to pull them towards the galaxy center, with the heavier diquark at xI,d closer to this center than does the lighter quark at xI1,d, as is shown in Figure 3. Here, the diquark-quark distance fm as was mentioned below Figure 2. This distance is fixed by the strong interquark potential (4.2) and is unaffected by gravity here; the core collapse situation in position b of Figure 2 is absent here.
Applying (3.1) the gravity pulled position d in Figure 3 yields the relative energy −ω0 given by
Figure 3. Illustration of the expanding galaxy scenario. The horizontal line represents the radius Rg inside the galaxy. The center of the galaxy is at Rg = 0. A proton in position c is moved to position d in the expansion. Xc denotes the position of the proton at position c, xI,c that of the diquark in this proton, and xII,c that of the quark. In the expansion, this proton c moves from Xc to Xd in position d. When gravity is turned on, the quark coordinates become xI,d and xII,d, respectively.
Figure 4. Illustration of the expanding universe scenario. The horizontal line represents the distance Ru from some unspecified inner region of the universe. Xe denotes the coordinate of a proton in a hydrogen molecule at position e, xI,e that of the diquark in this proton and x1I,e that of the quark. Xe lies in the middle between them. A second hydrogen molecule containing a second proton with its diquark and quark at the same xI,e and xII,e is also present. The coordinate of this second proton may be Xe− or Xe+ close to Xe. The second proton at Xe− will generate a positive relative energy that pushes its both quarks outwards. The coordinate of this proton lags behind. This process continues and this second proton arrives at position f with the coordinates xI,f and xII,f for the quarks and Xe−f for this second proton.
This energy has the same negative sign as the gravitational potential energy produced by matter inside Xd and hence reinforces it to become large enough to keep the outer stars of the galaxy from escaping. It is due to the lag of the “hidden” quark coordinates xI,d, x1I,d behind the observable proton coordinate Xd in Figure 3, caused by gravitational pull, in the expansion. This lag is caused by that the mass of the diquark and the quark is greater than the proton mass E0 itself, as can be seen below (2.1) and (6.1); the unobservable quarks experience greater gravitational pull than does the proton containing them. The difference between these two masses is taken up by the strong, confining potential Fbd(r) of (4.2) in the “hidden”, relative space.
The required energy is -w0 » -5.5E0 per proton, when averaged over the whole universe, and is generated for “free” at “no cost”. This −ω0 value leads to am ≈ 6 (  §6.1) so that the “whole” size of this gravitationally polarized proton in Figure 3 becomes Xd − xI,d ≈ 24 fm by (8.1). Further, the proton in a hydrogen atom moves in an orbit around the center of gravity of the electron-proton system. This orbit has a diameter ≈ 60 fm so that Xd is uncertain by this amount which is greater than the above 24 fm in Figure 3; there will be no observational difference.
This scenario, characterized by Figure 3, is also expected to have dominated in the earlier stages of the universe when it was smaller, hotter, denser, and fast expanding.
9. Accelerating Expansion of the Universe without Dark Energy
The observed accelerating expansion of the universe is currently considered to be due to assumed dark energy in the outer regions of the universe. This hypothetical energy may here be replaced by the positive relative energy −ω0 > 0 corresponding to am < 1/2 in (8.1).
As in Sections 6 and 8, consider the following idealized scenario. The expansion mechanism in Section 8 applied to a galaxy can analogously be used in some later stage of the development of the universe. In the outer part of the universe, its expansion leads to that the hydrogen gas density decreases and the gas temperature drops there. The expansion nearly comes to a halt. In this region, the gas is tenuous, cold and experiences very weak gravitational force.
Consider a proton in a hydrogen molecule of this gas. The configuration of this proton in position d of Figure 3, with the galaxy replaced by the universe, has now largely returned to its original form in position c of Figure 3. This proton is also similar to the neutron in position a of Figure 2 with am » 1/2 and the pulling energy -w0 < 0 approaching -w0 » 0. The situation is illustrated in position e of Figure 4.
Let a second hydrogen molecule arrive at position e simultaneously. It contains a second proton with its diquark and quark at the same xI,e and xII,e. Due to differences in the paths of these both molecules before they reach position e, this second proton may end up at nearby or where .
If its position is Xe+, then (8.1) yields a negative relative energy . Analogous to the situation below (8.1), this small energy tends to pull the quarks toward the left in Figure 4 into the inner region of the universe. This tendency is countered by the slightly higher gas pressure inside Xe. No appreciable net movement of the quarks is expected; the situation is stable.
If its position is Xe−, then (8.1) yields a positive relative energy . Contrary to the above Xe+ case, this small energy tends to push the quarks outward, towards the right part of Figure 4. This tendency is reinforced by the slightly higher gas pressure inside Xe; Δ increases and am decreases. The quarks move toward the outer region of the universe but the observable proton coordinate Xe− lags behind; this situation is opposite to that for Xe+ above or to that in position d in Figure 3 mentioned below (8.1). After some time, the unobservable quarks arrive at xI,f and xII,f in Figure 4. Application of (8.1) to the coordinates in position f indicates that am has decreased from 1/2−Δ at position e to <0 at position f so that the positive relative energy −ω0 by (8.1) has grown appreciably to >E0/2. This new, larger positive relative energy further pushes the quarks outwards at a greater pace and tends to increase the lag of the proton coordinate Xe-f behind the quark coordinates xI,f, xII,f in Figure 4. This by (8.1) generates still greater −ω0 > 0 which gives still greater pushing force on the quarks. This is an unstable, “run away” situation which can in principle continue forever. The driving relative energy is “free”, “costs nothing” and is not conserved. If this mechanism is applied to all protons in the outer regions of the universe, an acceleratingly expanding universe emerges.
Current data show that dark energy exceeds the energy of ordinary matter by a factor of ≈14 when averaged over the universe. Identifying this energy with −ω0, (8.1) leads to am = −13.5 (  §6.1). This value and (8.1) give the “whole” size of this second proton at position f in Figure 4 xII,f − Xe−f = 58 fm which is close to but still within the error margin ≈ 60 fm for Xe−f mentioned at the end of Section 8. This margin is small relative to the size of hydrogen atom and does not alter its behavior as a hydrogen atom.
This scenario, characterized by Figure 4, also dominates in the later stages of the universe when it has become large.
10. Scenarios in the Outer Regions of the Uiniverse
10.1. Possible Plasma Creation
As the above “run way” situation continues, the “whole” size of this proton, xII,f − Xe−f in Figure 4 or the proton-quark separation gets larger and eventually reaches 0.55Å, the radius of hydrogen atom. Since the observable proton coordinate Xe-f is located at the center of the atom, the quarks will in the same orbit as the electron. These two objects may come to be close to each other and interact via some so far unspecified mechanism. Since the electron is much lighter than a quark, which has a mass of about 2/3 of the proton mass mentioned below (2.1), it will be assumed that the electron gets ejected, noting that the ionization energy of hydrogen atom is negligible relative to the energy scale involved here. This hydrogen atom becomes ionized.
The relative energy corresponding to this case is by (8.1) −ω0 ≈ 1.29 × 1013 eV. In this scenario, hydrogen atoms in the expanding outer regions of the universe having this energy or greater turn into plasma.
10.2. Generation of Cosmic Rays
According to 10.1., the hydrogen gas in outer regions of the universe can expand until an average atom acquires an energy ~1.29 × 1013 eV, beyond which the atom becomes ionized. For atoms carrying higher energies, the gas becomes a tenuous plasma in intergalactic space. In this plasma, the protons continue to gain energy by the increasing −ω0 but the ejected electrons lag behind, inasmuch as gravity is absent in Dirac’s equation for an electron in a Coulomb field. This leads to a proton current flowing radially outwards which in its turn generates transverse magnetic field. As there are inhomogeneities in gas distribution in this region, the magnetic field will also vary in space.
As the “run away” expansion in Section 9 continues, some protons become very energetic. The magnetic field will cause them to move perpendicular to their path and to the direction of the magnetic field. Some of these high energy protons may be moved by magnetic fields such that they return to the inner part of the universe. Such protons can be a source of cosmic rays with energy > 1013 eV,
For a cosmic ray proton having a high energy of, say, 1020 eV, The “whole” size of this proton xII,f − Xe−f in Figure 4 given by (8.1) is still rather small, ~0.4 mm (unobservable), which is hardly detectable even if it were observable.
10.3. Diquark-Quark “Flip” and Magnetic Curtailment
In Figures 2-4, the diquark is closer to the center of the star, galaxy or universe; xI < xII. This is due to that the heavier diquark experiences a stronger pull towards the centers than does the lighter quark. This pulling force becomes very weak far out in the universe where the positive relative energy increases by the above run away instability pushing the quarks outwards with ever increasing proton energy and greater proton-quark separation or the “whole” size of the proton, x1I,f - Xe-f in Figure 4.
This may in principle go on forever. This pushing force is stronger on the heavier diquark than that on the quark and may eventually push the diquark past the quark farther out and ”flip” their positions from xI,f < xII,f in Figure 4 to xI > xII farther away. In this scenario, (8.1) shows that am, hence also −ω0, changes sign; the expansion force switches to contracting force. The expansion will slow down and change to contraction. If such processes continue, an “oscillating” scenario of the outer regions of the universe may emerge.
But even before this scenario, another one may take place. The protons with energy > 1013 eV are part of a plasma with transverse magnetic field in §10.2. Some of them will eventually move in a direction perpendicular to the direction of Ru in Figure 4. In that motion, the acceleration mechanism in Section 9 is no longer active. In this scenario, infinite proton energy and infinite expansion of the universe’s outer regions may be curtailed by the magnetic fields.
11. Summary-Scenario of the Universe
The present results lead to the following scenario for the universe.
1) The universe contains neither dark matter nor dark energy.
2) A heavy neutron star with mass = MTOV and radius = its Schwarzschild radius does not collapse into a mass singularity and may be called a “black” neutron star. If every neutron star gets heavier, becomes a black hole and passes through this stage in its development, there will be no mass singularity in the universe, in agreement with that quantum mechanics does not allow such a singularity.
3) The run away instability in Section 9 provides a mechanism for an accelerating expansion of the universe. The driving positive relative energy is “free” and “costs nothing”. The density and temperature of the universe will decrease with time. This scenario may in principle go on forever.
4) In outer parts of the universe, fast expanding hydrogen gas may turn into plasma and part of it may become cosmic rays. There may also exist a scenario in which the above free expansion can halt and eventually reverts to contraction.
The above results are derived phenomenologically by joining SSI to aspects of general relativity. A formal integration of these both theories is beyond reach; no quantum gravity theory exists presently.