The theory of general relativity GR  has been proved to be accurate in most cases, any attempt to complement it must be taken carefully.
Despite its immense success, GR theory lacks in explaining the behavior of the gravity inside black holes . Besides, the calculation of the acceleration of gravity outside a black hole or any other massive celestial object having a strong gravitational pull is the subject of numerically solving complex nonlinear equations . Moreover, finally, there is a bridge that has yet to be built between GR and quantum mechanics .
The present theory has the merit and the novelty of helping to explain or solve the three points above, besides offering a simple computational way to calculate gravity inside and outside black holes, or in any other massive object, using Fourier transforms, such as explained later.
As a start point for the theory of black holes, the author considers the existence of two kinds of energy sources for the gravitational field: The rest mass, which converts into energy by E = Mc2 and the energy due the motion of particles.
For calculations are used the Fourier transform, widely used in quantum mechanics for representing the wave function in the position or momentum space, among many other uses .
Regarding dark matter and dark energy, both have been the subject of discussion for a long time , with some works done in an attempt to relate the two effects .
Based on the work discussed above, it is possible to consider, by hypothesis, the possibility of the quantization of the spatial frequency k, that is, only occurring at fixed intervals, such as km, 2km, 3km, where km stands for “minimum spatial frequency”.
The proposed quantization, when applied in the inverse Fourier transform, reveals areas of extra gravitational force spaced at 2π/km, which fall as inscribed in an inverse distance envelope (1/r). Therefore, choosing the correct value of 2π/km in the range of kiloparsecs (kpc) is sufficient to explain the practically flat rotation curves of most galaxies . Also, the appearance of adjacent regions of negative gravitational force due to the proposed quantization is also verified. Areas of negative gravitational force cancel, on average, the effect of areas of positive gravitational force, which means that an object traveling between galaxies experiences, on average, a null gravitational effect.
The null average effect of gravity between galaxies allows us to demonstrate that the mean luminosity density of the universe is a sufficient factor to explain the repulsive force responsible for the accelerating expansion of the universe.
By now, it is not clear what causes the universe to expand, puzzling the scientists even more since the discovery that the expansion is in fact accelerating  .
The repulsive force of light as responsible for the accelerating expansion of the universe is mostly unconsidered due to its small magnitude compared to gravitational effects, being necessary the elaboration of alternative concepts such as negative vacuum pressure or dark energy, among others .
However, considering the average annulment of the force of gravity between galaxies, it is clear that luminosity could be a predominant factor and the only one necessary to explain the expansion of the universe. The physical model developed in the last part of the present paper has yielded results that agree very well with values that are within the current acceptable limits for the Hubble constant ≈ 67 - 73 km∙s−1∙Mpc−1, obtained respectively by  and .
The present model excludes the existence of non-baryonic matter, or the existence of dark matter or dark energy, further concludes that the expansion of the universe follows a coasting universe model, expanding forever, never contracting.
2. Decomposition of Matter in k-Space
The k-space or reciprocal space is the spatial frequency domain of a spatial Fourier transform, as shown in Equation (1). It is possible to take a distribution of mass in position space, and through the use of a spatial Fourier transform to decompose it in different spatial frequencies.
As our object of study are mostly radially symmetric spherical bodies, such as stars, black holes, and so on, it is showed, using spherical coordinates, that the Fourier transform for such a case can be simplified as in Equation (2) .
where r is the radial coordinate, and k is the spatial frequency in k-space, which unit is m−1. For a spherical object of mass M, radius R, and volume V, the radial density D(r) is as shown in Equation (3).
Substituting Equation (3) in Equation (2) and resolving for D(k), we get
where R is the radius of the spherical object, k is the spatial frequency and M is the mass of the object. The function D(k) has unit kg and can be understood as being the decomposition of the mass M in different spatial frequencies k’s. Equation (4) is plotted in Figure 1.
From Figure 1, we see that the decomposition of an object of mass M, with uniform density and radius R in k-space, has a maximum value of M and a “bandwidth” inversely proportional to the radius R. However, we can observe the existence of negative values, which does not necessarily mean negative mass, which has no physical significance. Instead negatives values in k-space are understood as spatial frequencies shifted 180˚ in phase, being necessary for the correct reconstruction
Figure 1. Decomposition in k-space of an object with mass M, radius R, and uniform density.
of the object after the application of the inverse Fourier transform, as will be explained below.
Reconstruction of the Gravitational Potential
To facilitate the understanding of topics to be discussed later, we will normalize the function D(k), that is, divide it by M and introduce the function F(k) which is simply D(k)/M.
The inverse Fourier transform associated with the Fourier transform of Equation (2) is defined as
Upon inspection, it is verified that the right transformation for the classical Newtonian gravitational potential is given by Equation (6).
Which gives the correct value for the classical Newtonian gravitational potential, of the form GM/r, where G is the gravitational constant.
Plotting the Equation (6) along with its negative derivative with respect to r ( ), which is the force per unit mass i.e., gravitational acceleration, we obtain as shown in Figure 2 the correct value for the gravitational potential and force per unit mass (gravitational acceleration) in the classical Newtonian form. As a further verification of the presented method, in Figure 3
Figure 2. Gravitational potential (Solid line) and force per unit mass (Dashed line) obtained from Equation (6) for a solid sphere.
Figure 3. Gravitational potential (Solid line) and force per unit mass (Dashed line) for a hollow sphere.
is plotted the value of the gravitational potential for a hollow sphere using the pair of direct and inverse Fourier transforms given respectively by the Equations (2) and (6) and using Equation (3) considering the volume V as the volume of the outer shell. Where the integration limits for the Equation (2) are taken between 0.9R and R. The respective value of its negative derivative (acceleration) is also showed in the dashed line.
Such previous examples prove the validity of the adopted method. Also, other validated cases are not shown here.
The same concept could be expanded to objects that do not necessarily have radial symmetry, for that it is enough to use the correct pairs of Fourier transforms for such cases. These results indicate that any radial distribution of matter can be decomposed in k-space, using Equation (2), and by using Equation (6) to reconstruct the corresponding classical Newtonian gravitational potential, which by the negative value of its differentiation in space results in force per unit mass (acceleration).
3. Motion Energy with Increasing in Mass Density
In the relativity theory , the total energy E of an object is
With Mc2 being the rest mass energy (from now on the mass M is also considered the rest mass) and its “motion energy” given by the product of its momentum and the speed of light c. Based on that, the theory presented here assumes that the use of Equation (6) only considers the rest mass M or in a similar way, the “rest energy” Mc2 to correctly calculate the classical Newtonian gravitational potential.
However, the motion energy of an object is not only due its translational motion, being also contained in the form of quantum motion of its constituent particles.
Such phenomenon is described in the well-known “Uncertainty Principle” , shown in Equation (8) and which essentially states that the greater the certainty in determining the momentum of the particle, the less certainty in its position and vice versa.
where ħ is the reduced Planck constant.
This principle can be used to estimate the average values for the motion energy of a quantum particle, which is inversely proportional to the dimension occupied by the particle ∆x.
It is assumed here that the increase in mass density will act in reducing the radius of the particles (more particles packed together) and consequently increasing the motion energy of its constituent particles. This can be clearly seen in Figure 1, where the “bandwidth” is proportional to 1∕R, meaning that as the radius of the spherical object decreases the amount of energy in the “momentum space” k spreads out, having consequently more energy in the form of quantum motion.
Gravitational Potential due the Motion Energy of Quantum Particles
The spatial frequency k is directly related to the momentum of a particle through the De Broglie relation hk, where h is the Planck constant. As the Fourier transform indicated in Equation (2) is a 3D transform, the same is valid for the inverse transform. So, integrating the function all over the 3D space of the spatial frequency k, we obtain the expression in Equation (9), which indicates the gravitational potential due to the motion energy.
where physical constants and extra terms were adjusted in order to give the same dimension as in Equation (6), since the unit for k is m−1.
As the Equation (7) indicates that the rest mass energy and the momentum energy are perpendicular to each other we can suppose that the angle total energy E is shifted in arctan (Momentum energy/Total Energy) in relation to the rest mass. This all can be summarized in the Equation (10).
where UT is the complete solution for the gravitational potential. The terms in the square root are respectively proportional to the rest mass energy and motion energy and the term inside the arctan function represents the ratio between them. The values 2 multiplying the motion energy were chosen in order to adjust the equation to the GR theory as seen later. The Equation (10) is reduced to the classical gravitational potential as in Equation (6) in the limit , where RS is the Schwarzschild radius.
4. Non-Rotating Black Holes
Making use of the Equation (10), we are able to plot for the condition for the formation of a black hole , being .
For comparison we use the relativistic gravitational potential as obtained by  and indicated in Equation (11).
By analyzing Figure 4, which is the plotting of Equation (18), we realize that the gravitational force only goes to infinity at the border of a black hole.
Which gives a pretty good agreement and besides that offers an understanding of how the gravitational field behaves internally to the massive body, showing that the gravitational potential only tends to infinity in the border of the black hole.
Such procedure can be used for calculation of the gravitational field near massive bodies approaching the Schwarzschild limit.
An important question that can be raised throughout this topic is how the mass contained in an object can deflect the space-time fabric to exert a force at the distance on some other object? As an answer, we can infer that any object that has mass has constituent particles that are always in motion, and that movement would be transmitted to the space-time fabric in a way to combine with the other vibrations originated by other particles building up a resulting total deflection of the space-time fabric in all internal and external regions of such object.
The only difference between the deflection caused by rest mass and the motion energy of these particles would be the way they "move", causing different effects on the space-time fabric that would later be transmitted to all regions of the space.
Up to here, we have worked with gravitational effects more evident over small
Figure 4. Gravitational potential field for a black hole. Numerical solution using Equation (10) (Circle), and the relativistic gravitational potential using Equation (11) (Solid line) for comparison.
distances on cosmological scales 1 kpc, contained in smaller dimensions than those of our solar system. Since the extra gravitational effect of black holes seems more obvious in their surroundings, it became asymptotically more similar to the classic Newtonian model for greater distances, since the extra effect due to the motion energy of the quantum particles decays more abruptly as distance increases.
However, for larger cosmological scales in the order of kpc, we verify the existence of the phenomenon known as dark matter, and on even larger scales (intergalactic medium) we verify the existence of the phenomenon of dark energy.
In the next topic, we will demonstrate that both dark matter and dark energy can be caused by the existence of the quantization of the spatial frequency k, which can be understood as a limitation of the matter moving at all values of spatial frequency k, or a limitation of the space-time field of being excited in all spatial frequencies, with such excitations only occurring in steps.
We must clarify here that the quantization of the spatial frequency does not affect the results obtained in the present topic, since it only has effects on a scale much larger than that of our solar system.
5. Quantization of the Spatial Frequency k
We here postulate the existence of a minimum value for the spatial frequency k, being the other values multiples of this minimum value. Therefore, the spatial frequency can only acquire discrete values, such as km, 2km, 3km, and so on.
This leads to a discrete representation of the integral in Equation (6), with an increment of km, as indicated in Equation (12).
The value of km should be considered as very small, in such a way that the value of 2π/km is on the galaxies size scale (kpc). Therefore, the celestial objects considered in this work are galaxies in general, where the effects of the quantization of the spatial frequency are evident.
5.1. Dark Matter
To demonstrate how the quantization of the spatial frequency k allows us to deduce the effect of dark matter, we will take the example of the Milky Way.
The distribution of mass in the bulge of the Milky Way Db(r) can be modeled according to Equation (13) .
With and . This equation can be approximated using a bell curve, as indicated in Equation (14).
In Figure 5, we can see that the two curves are similar, and the use of one or the other is indifferent to our goal of showing the effect of the quantization of the spatial frequency k in the behavior of the gravitational force. However, we will choose to use the bell curve because it is easier to solve analytically.
After the application of Equation (14) in Equation (2), with , we obtain the mass distribution in k-space, also a bell curve, in terms of the spatial frequency k, as shown in Equation (15).
Applying this value in Equation 12, with and comparing the effect with the classic Newtonian force, we get Figure 6. From Figure 6, we can see that the effect of the quantization of the spatial frequency creates regions of extra gravitational force that decrease inscribed in a 1/r envelope. This phenomenon is like that found in signal analysis, more specifically in the sampling theory, where the discretization of the signal in one domain, produces extra spectrums in the other domain .
Figure 5. Equation (13) (Solid line) and Equation (14) (dashed line) for comparison.
Figure 6. Force per unit mass (acceleration) in the classical Newtonian (solid line) and new proposed model (dashed line) with . The point line shows that the resulting force is inscribed in a 1/r envelope (points).
Being therefore sufficient the occurrence of these regions to explain the fact that most spiral galaxies have flat rotational velocity curves, since such galaxies have a high concentration of mass in a volume in their center (bulge) as modeled in the present case. Besides, there are adjacent regions of negative gravitational force, which effect is to expel objects orbiting the galaxy towards regions of positive gravitational force, being therefore most likely responsible for the creation of voids in the galaxies, where there is a lower concentration of matter.
The present model is perfectly symmetrical, requiring more realistic models for the precise determination of where the matter will accumulate. More precisely, it is necessary and suggested by the author to run complex computational simulations to analyze the evolution of spiral galaxies, considering the quantization of the spatial frequencies.
Now applying the relationship , where is the rotational velocity, a(r) is the acceleration, and r is the radial distance, we obtain Figure 7 where the negative values (gravitational repulsion) are not considered.
From Figure 7, fringes of extra gravitational force are created around the central bulge, providing areas where the velocity of rotation is practically constant. These fringes are also probably responsible for the creation of galactic arms, where areas of lower matter density are those with repulsive gravitational force (negative velocity). Regardless, it is suggested again as confirmation of the present theory the execution of computer simulations of the evolution of galaxies, considering the
Figure 7. Rotational velocity using the relation in the classical Newtonian way (solid line) and the new proposed model (dashed line) with .
quantization of the spatial frequency k. The quantization of the spatial frequency used in Figure 6 and Figure 7 is , which is not necessarily the exact value, but apparently close to the real value. Computer simulations must be performed in order to determine a more accurate value.
5.2. Dark Energy
Knowing that the acceleration acc, in terms of the displacement variable s, can be expressed as indicated in Equation (16).
Integrating both sides and solving for the velocity v, we obtain
Therefore, we have an expression that relates the variations of acceleration in space s, with velocity v acquired by a mass. Consequently, we notice that if the integral of the acceleration in space is limited ( ), the velocity v will have a constant value, causing a zero acceleration of a mass placed in that field. Now considering a(r) as this acceleration, we have the necessary conditions to verify if the acceleration of gravity with the effect of the quantization of the spatial frequency causes a net acceleration or not. The force per unit mass a(r) (acceleration) due to the quantization varies positively and negatively inscribed in an envelope of 1/r (Figure 6). Thus, it can be proved by integral convergence rules, that the velocity acquired by a mass placed in this gravitational field does not exhibit any net acceleration.
According to what was discussed above, we must assume that, on average, the net gravitational force between galaxies is canceled. It hence requires a force of repulsion much smaller than what is currently used to explain the accelerating expansion of the universe.
The next topic goes on to explain that the average luminous density of the universe is a sufficient factor to explain the commonly accepted rates of acceleration of the universe.
6. Luminosity as the Factor for Universe Expansion
In this topic, the influence of luminosity on the expansion of the universe is evaluated considering the net gravitational effect of the attraction between galaxies as null, as previously discussed.
Luminosity is defined as the total electromagnetic power emitted by a source, and for astronomy, such sources can be stars, galaxies, supernovae, or any other astronomical object . The unity of luminosity is W, and this power is due to photons that individually carry a discrete amount of energy given by hf. Typically, an astronomical object such as a star emits its luminosity in a broad range of frequencies (bolometric luminosity). Each photon is associated with a momentum given by hf/c, where c is the speed of light in the vacuum. It means that a stream of photons emitted by a source can exert a force over an object, being this force proportional to the cross-section area S of the object which lies perpendicular to the direction of the stream of photons. For an object apart from the light source, the power per square meter (W/m2), commonly named “apparent brightness,” is proportional to the inverse of the square of the distance (1/r2). This relation is quite useful for our study, since we can borrow some concepts used in classical gravity, which also follow an inverse square law.
Considering that the universe is homogeneous and isotropic when viewed at large scales, we can make use of the shell theorem, which allows significant simplifications for objects subjected to gravitational forces inside or outside a shell or sphere.
We realize, in Figure 8, with the application of the shell theorem, that the entire universe outside the sphere does not exert any gravitational force over the astronomical object A; and all the gravitational effect of the objects inside the sphere over the object A can be replaced by an equivalent object in the center of the sphere with a mass equivalent to the sum of all masses within the sphere.
Since apparent brightness (W/m2) also follows the inverse square law, just like the force of gravity, we can, by analogy, use the same simplification as for the luminosity.
Figure 8. Shell theorem for gravity applied for a homogeneous and isotropic universe (a), and the equivalent gravitational effect in (b). This simplification is also valid for the “luminosity push”.
Using again Figure 8, we can also state that the equivalent “luminosity push” felt by an object A in Figure 8(a) is derived not taking into account all the luminosities sources outside the sphere, but replacing all the luminosity sources inside the sphere by an equivalent one at the center with a luminosity equivalent to the sum of all luminosities inside the sphere (Figure 8(b)).
This simplification, again, is only possible because, just like gravity, the apparent brightness also follows the inverse square law (1/r2). However, this simplification has some constraints; among them, we can list:
1) The presence of intergalactic dust causes an attenuation of the order 0.01 mag∙h∙Gpc−1 , being negligible for a radius R < 1 Gpc.
2) The mean free path M.F.P for a photon traveling through the local universe is given by Equation (18), Where Rg is the mean galaxy radius, and n is mean density of galaxies per volume. n can be estimated as 0.014 h3∙Mpc−3 , and Rg as 0.015 Mpc. Replacing these values into Equation (18) we obtain M.P.F = 101.05 h−1 Gpc.
So, assuming the two constraints above, we can fairly say the universe is transparent for our analysis, being this way, the simplification given by Figure 8 is valid.
6.1. Derivation of the Hubble Parameter
Since the universe is expanding, the mean luminosity density given by varies in time. It leads to a total luminosity L inside a sphere of a comoving radius R(t) given by:
Generating a luminosity pressure P
According to , the Tolman relation, which expresses the redshift dependence of the surface’s bolometric luminosities of the galaxies B, is
where z is the redshift, and B0 is the surface’s bolometric luminosity at z = 0. The Tolman relation implies that galaxies at redshift z will appear fainter by a factor of (1 + z)4. This relation is valid for an observer placed at z = 0. If the observer is moved to a position at the point A, at a redshift z, in relation to the center of a sphere, as indicated in Figure 8, we must reverse this relation, indicating that in the past, the surface brightness of the galaxies was brighter. Knowing that, and that the scale factor a and the redshift are related by , the luminosity pressure, as indicated in Equation (20), should be modified by
A galaxy with a face-on surface S, as indicated in Figure 9, opacity k and mass M, is subject to an acceleration given by
Replacing L of Equation (19) in Equation (22), produces
Applying the scale factor , we get the following relationship for acceleration
Figure 9. Galaxy with Area S,opacity K and mass M in the surface of a sphere of radius R.
The Equation (25) is derived in .
Where t1 is the present time and is the luminosity density of the present time. Replacing Equations (24) and (25) in Equation (23) yield:
The Equation (26) is correct for galaxies, which are partially opaque . In general, its opacity depends on the wavelength and topology of the galaxies. Most of the opacity of galaxies is related to absorption and scattering of photons, and it is directly related to its dust content, occurring mostly for higher frequencies in the UV-Blue range of wavelengths. Extensive studies showing such relations are found in  .
Studies conducted by  placed the opacity K of spiral galaxies in the range 0.48 ± 0.15. Integrating Equation (26) in a and dividing by a2, produces
Which is an expression for the Hubble constant since , where B is a term to be determined. The term B can be determined using initial conditions, assuming that the gravitational force between galaxies is on average, nullified. This means that for galaxies, placed, as shown in Figure 10, the luminosity pressure is dominant, where the scale factor is ≈2Rg, with Rg being the average galactic radius. We can call the process, where the galaxies start to move away as “luminosity dominance,” and the time of its occurrence as tD, during this time we have
where is the mean distance between galaxies at the “luminosity dominance” time td, and is the mean distance between galaxies at the present
Figure 10. Galaxies apart in 2.Rg at time tD, when the pressure of light begins to dominate over the gravitational force.
The mean distance between galaxies at time tD, as seen in Figure 10, is approximately twice as the galactic radius, being the average galaxy radius commonly assumed to be 10 kpc  . Then, , and the present time average distance between galaxies can be derived using .
where n is the present time density number of galaxies. Another important consideration to be made is assuming that at the time of the “luminosity dominance” tD the recessional velocity is zero since it is the time when the luminosity repulsive force starts to overcome the gravitational force. So, consequently
Applying Equations (28) and (30) in Equation (27), we obtain
Solving for B yields
Replacing this in Equation (27) we obtain
Rearranging for the Hubble constant
6.2. Derivation of the Scale Factor A
OBS: The factor scale factor a here presented is not the acceleration.
Considering the Equation derived for the Hubble constant (Equation (34)), we can verify that the second term inside the square root is equal to the first term when t =tD. For t >tD the first term increasingly dominates as the time increases. So, for we can approximate the Equation (34) as
Rearranging we obtain
Equation (36) can be solved for the scale factor a. We can “guess” a solution form for the scale factor a as being
where t is time, D and U are constants to be determined, replacing Equation (37) in Equation (36) and rearranging we get
Since the right term is a constant, we need to set U to force the left term to be also a constant, this is done doing
Replacing U = 1 in Equation (38) and solving for D yields
Therefore, the approximated scale factor for is
The Equation (41) is a good approximation for values of but is not the whole answer for a(t), since if we derive twice the value present in Equation (41) it will become zero, what is not valid, since the acceleration value for is given in Equation (37). Therefore, we can refine the value of a(t) doing the same procedure as before for the acceleration term in Equation (37). Calling the acceleration term as aacc(t), and guessing the answer to be in the form , we get
Using the same procedure as in Equation (39) leads to U = 5/2, and solving for D produces
Adding this term to Equation (41), yields
Considering our initial time as tD= 0, we have the value of a(tD) in such a case. Adding this to Equation (44) we obtain the complete solution for the scale factor as a function of time t expressed in Equation (45).
6.3. Results of Acceleration for the Milky Way
According to , the luminosity density at the present time is estimated in and the number density of galaxies is estimated to be . Using the Equation (29) with , we obtain . Therefore, considering , we get
 have reported that the mass of the Milky Way galaxy is 1.5 trillion solar masses, being about 90% of this mass dark matter. So, the visible, baryonic mass is about 0.15 trillion solar masses. The present work uses only the visible matter for the calculation, since as reported in the present work, the extra dark matter is only consequence of the effect of the quantization of the spatial frequency, not being accounted as inertial mass.
 also reported a Milky Way radius of about 39 kpc, over twice as much as previously reported, being most of this extra part due to the globular cluster out in the galactic halo. Therefore, for the calculation, only the disc radius is considered, being about 15 kpc , resulting in a face-on area of about 148 kpc2. For the opacity, k is considered the values obtained by  being about 0.48 for spiral galaxies. So, in resume, we have the following parameters, for the Milky Way galaxy, in SI base units:
Applying those values for Equation (34), for the present time t1, with a(t1) = 1, we obtain
which agrees very well within the interval of acceptable values for the Hubble constant ≈ 67 - 73 km∙s−1∙Mpc−1, obtained respectively by  .
Using Equation 34, we can also plot the normalized value of the Hubble parameter (H(a)/H(t1)) against the scale factor a, with the result indicated in Figure 11.
Using the Equations (34) and (45) we can also plot the normalized value of the Hubble parameter (H(t)/H(t1)) and scale factor a against time in years (Yr), with the result indicated in Figure 12.
Figure 11. Normalized value of the Hubble parameter H(a)/H(t1) plotted against the scale factor a.
Figure 12. Normalized value of the Hubble parameter (H(t)/H(t1)) in solid line and scale factor (point line) plotted against time in years.
Looking at Figure 12, we realize that the universe is coasting type, and thus will expand steadily and never contracts. Analyzing the Equation (35), we realize that the only term that varies from one galaxy to another is the expression S.K/M. As the Hubble constant is approximately the same for a wide range of galaxy types, not depending on their inclination too, we are led to conclude that the term S.K/M is approximately constant for all galaxies and remains approximately constant no matter the inclination. Analyzing in a simple way, we conclude that this could be the case since for inclined galaxies, the cross-section area S is reduced, and since the photons have to cross more matter, intuitively, the opacity K will grow at the same proportion, keeping the term S.K approximately constant.
The developed theory proved to be entirely satisfactory in developing a model of the gravitational force for massive objects, such as non-rotating black holes, in line with Einstein’s theory of relativity (GR). As merit, the present theory is able to determine the behavior of the acceleration of gravity within a black hole. Besides, the theory offers a straightforward and simple method to numerically calculate the effect of gravity near massive objects. It also provides a “bridge” between gravity and quantum physics, since the deflection of the space-time fabric can be understood as originated in the quantum scale.
The developed theory also proved to be quite satisfactory in the development of a model for the dark matter, besides demonstrating the average absence of gravitational force between galaxies. From the absence of average force between galaxies, it was possible to demonstrate that the average luminosity density of the universe is enough to explain the currently accepted rates of expansion of the universe. As a suggestion for future studies, the author recommends simulations of evolution of the galaxies, taking into consideration the effect of quantization of the spatial frequency k.
 Brustein, R. and Medved, A.J.M. (2019) Resisting Collapse: How Matter inside a Black Hole Can Withstand Gravity. Physical Review D, 99, Article ID: 064019.
 Wang, B., Abdalla, E., Atrio-Barandela, F. and Pavón, D. (2016) Dark Matter and Dark Energy Interactions: Theoretical Challenges, Cosmological Implications and Observational Signatures. Reports on Progress in Physics, 79, Article ID: 096901.
 Riess, A.G., Filippenko, A.V., Challis, P., Clocchiatti, A., Diercks, A., Garnavich, P.M., Gilliland, R.L., Hogan, C.J., Jha, S., Kirshner, R.P. and Leibundgut, B.R.U.N.O. (1998) Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant. The Astronomical Journal, 116, 1009.
 Perlmutter, S., Aldering, G., Goldhaber, G., Knop, R.A., Nugent, P., Castro, P.G., Deustua, S., Fabbro, S., Goobar, A., Groom, D.E. and Hook, I.M. (1999) Measurements of Ω and Λ from 42 High-Redshift Supernovae. The Astrophysical Journal, 517, 565.
 Freedman, W.L., Madore, B.F., Gibson, B.K., Ferrarese, L., Kelson, D.D., Sakai, S., Mould, J.R., Kennicutt Jr., R.C., Ford, H.C., Graham, J.A. and Huchra, J.P. (2001) Final Results from the Hubble Space Telescope Key Project to Measure the Hubble Constant. The Astrophysical Journal, 553, 47.
 Blanton, M.R., Dalcanton, J., Eisenstein, D., Loveday, J., Strauss, M.A., SubbaRao, M., Weinberg, D.H., Anderson Jr, J.E., Annis, J., Bahcall, N.A. and Bernardi, M. (2001) The Luminosity Function of Galaxies in SDSS Commissioning Data. The Astronomical Journal, 121, 2358-2380.
 Reynolds, A.P., Parmar, A.N., Hakala, P.J., Pollock, A.M.T., Williams, O.R., Peacock, A. and Taylor, B.G. (1999) The EXOSAT Medium-Energy Slew Survey Catalog. Astronomy and Astrophysics Supplement Series, 134, 287-300.
 White III, R.E., Keel, W.C. and Conselice, C.J. (2000) Seeing Galaxies through Thick and Thin. I. Optical Opacity Measures in Overlapping Galaxies. The Astrophysical Journal, 542, 761-778.
 Calzetti, D., Armus, L., Bohlin, R.C., Kinney, A.L., Koornneef, J. and Storchi-Bergmann, T. (2000) The Dust Content and Opacity of Actively Star-Forming Galaxies. The Astrophysical Journal, 533, 682-695.