In the previous exploration of the fundamental interactions of the large-scale universe, we have discovered the existence of a $\gamma$ factor that affects Newton’s universal gravitation or the cosmic constant Einstein’s equations through the invariance of the gauge transformation of the principal and associated fiber bundles structure of the universe. It makes the calculated value of Newton’s universal gravitation in space-time larger than the galaxy scale multiplied by $1/\gamma$ to match the measured value, and makes the cosmological constant of Einstein’s equations no longer a constant but a function that increases with the increase of the scale of cosmic regions. The $\gamma$ factor can be used to explain that the so-called dark matter or dark energy phenomenon is actually caused by the deviation caused by the fact that Newton’s law of universal gravitation or Einstein’s equations are no longer completely valid in space-time scales larger than or equal to galaxies [1] . Dark matter and dark energy do not actually exist. The reason why the dark matter and dark energy hypothesis was proposed began because humans did not understand the role of the $\gamma$ factor [2] - [12] .

But what is the physical meaning of the $\gamma$ factor? The author of this article can first briefly answer this question: it reflects that the gravitational constant is not a fixed value but a variable! The gravitational constant is actually a function of the scale of cosmic regions! If we have such a conceptual shift, then the mystery of dark matter and dark energy can be easily explained.

The specific reasoning and analysis of the physical significance of the $\gamma$ factor and the conclusion that it leads to a non-constant gravitational constant is precisely the task of this article. Below, the author will elaborate on the description.

In references [1] [13] [14] , the author proposed a differential geometric representation of the cosmological landscape of the principal and associated bundles. A connection on the principal bundle $P\left(M,G\right)$ is a 1-form field ${\omega }_U$ that specifies a Lie algebra $\mathcal{G}$ value of $C^{\infty }$ on $U$ for a local trivial $T_U:{\pi }^{-1}\left[U\right]\to U\times G$, that is, a connection on $U\subset M$. In this case, if $T_V:{\pi }^{-1}\left[V\right]\to V\times G$ is another local trivial, that is, $U\cap V\ne \varnothing$, and the transition function from the local trivial $T_U$ to $T_V$ is $g_{UV}$, then the generalized Gauge equation holds [15] ,

${\omega }_V\left(Y\right)=\mathcal{A}{d}_{g_{UV}{\left(x\right)}^{-1}}{\omega }_U\left(Y\right)+L_{g_{UV}\left(x\right)\ast }^{-1}{g}_{UV\ast }\left(Y\right),\forall x\in U\cap V,Y\in T_{x}M$ (1)

where $L_{g_{UV}\left(x\right)}^{-1}$ is the inverse mapping of the left translation $L_{g_{UV}\left(x\right)}$ generated by $g_{UV}\left(x\right)\in G$, and $L_{g_{UV}\left(x\right)\ast }^{-1}\equiv {\left(L_{g_{UV}\left(x\right)}^{-1}\right)}_{\ast }$.

Furthermore, the principal bundle structure is composed of the structure group, the principal bundle manifold, and the bottom manifold has a principal bundle section as the gauge potential, a curvature as the gauge field strength, and satisfies the curvature gauge transformation relation between two different regions $U$ and $V$. For example, in general, if $g_{UV}:U\cap V\to G$ is a local trivial transition function from $T_U$ to $T_V$, ${\Omega }_V$ and ${\Omega }_U$ are the curvatures of the bottom manifold in the $V$ region or the $U$ region, respectively, then on the bottom manifold $U\cap V$ has:

${\Omega }_V=\mathcal{A}{d}_{g_{UV}^{-1}}{\Omega }_U$ (2)

where $\mathcal{A}{d}_{g_{UV}^{-1}}$ is defined as a linear transformation on the Lie algebra of the structure group $G$. If the structure group is a matrix group, then (2) can be further expressed as:

${\Omega }_V=g_{UV}^{-1}{\Omega }_U{g}_{UV}$ (3)

In fact, from Equations (2) and (3), we have:

${\Omega }_V=\mathcal{A}{d}_{g_{UV}^{-1}}{\Omega }_U=g_{UV}^{-1}{\Omega }_U{g}_{UV}$ (4)

To prove the above formula, we only need to prove that the tangent vectors on both sides of the above formula at point $e$ (corresponding to the unit element $e$ of the Lie group at $t=0$) in the tangent space (i.e. Lie algebra) are the same. To this goal, we can calculate the following formula to prove it (although more detailed proof can be found in ref. [1] ):

$\begin{array}{l}\text{Right tangent vector}={\frac{\text{d}}{\text{d}t}|}_{t=0}\left[g_{UV}^{-1}{\Omega }_U{g}_{UV}\right]={\frac{\text{d}}{\text{d}t}|}_{t=0}\left[I_{g_{UV}^{-1}}{\Omega }_U\right]\\ =I_{g_{UV}^{-1}{}^{*}}\left[{\frac{\text{d}}{\text{d}t}|}_{t=0}{\Omega }_U\right]=\mathcal{A}{d}_{g_{UV}^{-1}}{\Omega }_U=\text{Left tangent vector}\end{array}$ (5)

Therefore, under the cross-sectional transformation $\sigma$, on the bottom manifold (that is, the world we assume, with $U\subset V$, for example, $U$ represents the solar system region, and $V$ represents the region greater than or equal to the galaxy), having ${\omega }_U\to {\omega }_V$, that is, Equation (1) holds, and here $g_{UV}^{-1}$ and $g_{UV}$ are related to a gauge transformation.

Furthermore, it can be proved that Equation (3) is equivalent to the following similarity transformation Equation (6) [13] [14] :

${\Omega }_V=g_{UV}^{-1}{\Omega }_U{g}_{UV}\iff {{\hat{F}}^{\prime }}_{\mu \nu }=W{\hat{F}}_{\mu \nu }{W}^{-1}$ (6)

where ${{\hat{F}}^{\prime }}_{\mu \nu }$ is defined as the Lie algebraic matrix representation of the gravitational gauge field strength in region $V$ on the bottom manifold; ${\hat{F}}_{\mu \nu }$ is defined as the Lie algebraic matrix representation of the gravitational gauge field strength in region $U$ on the bottom manifold. $W^{-1}$ and $W$ can be regarded as the matrix expression of a certain gauge transformation, which can be called a gauge similarity transformation.

The author should emphasize that this gauge similarity transformation is affected by the regional scale, or depends on the change of regional scale. The scale of different regions $U$ or $V$ determines different gauge potentials or field strengths. The generalized gauge equation and the above gauge similarity transformation are the transformation relationship between the gauge potentials or field strengths of two different regions at the same space-time point of the bottom manifold (our universe). Not only that, after repeated consideration, the author also believes that: if the structure group is taken as the general linear group ( $GL\left(m,ℂ\right)$), then the generalized gauge transformation between the basic interactions corresponds to the gauge transformation between several large subgroups (for example, $U\left(1\right)$, $SU\left(2\right)$, $SU\left(3\right)$, $O\left(1,3\right)$), while the gauge transformation within a basic interaction (such as electromagnetic force) only corresponds to the gauge transformation within a “large” subgroup. For example, the gauge transformation within the electromagnetic interaction or weak interaction only involves the $U\left(1\right)$ or $SU\left(2\right)$ group [16] [17] , which are only subgroups of $GL\left(m,ℂ\right)$. These components are just a kind of “representation” of the original connection or curvature on the principal bundle in the real world (bottom manifold); in this sense, “quantization” or “classicization” is just a natural “representation” expression of the original connection or curvature selection in different regions of the real world, and it is not that they can replace each other, or that any physical field can be or needs to be quantized, especially gravity. It should be emphasized that the original connection (gauge potential) or curvature (gauge field strength) on the principal bundle is gauge invariant, and the gauge potential or field strength on the bottom manifold is just the projection component mapped (gauge selection) by the original connection or curvature of the principal bundle choosing different cross sections. This is the physical meaning of the unification of the basic interactions in the world [1] [13] [14] .

Based on the above assumption, this article believes that dark matter and dark energy are hypothetical hypotheses generated by the differences in gauge transformations of gravitational gauge potentials or gauge field strengths between different regional scales on the bottom manifold. There is no dark matter or dark energy, only changes in the gauge potential or field strength on the bottom manifold caused by gauge transformations. The original Einstein equations, Newton’s second law or the law of universal gravitation may undergo changes in form under such transformations, and may differ from their original form. There exists a difference in the $\gamma$ factor, which leads to the hypothesis of the existence of dark matter or dark energy [1] .

In fact, the derivation of the $\gamma$ factor can be achieved by constructing gauge transformations as follows:

$\begin{array}{c}\begin{array}{ccc}{{\hat{F}}^{\prime }}_{\mu \nu }& =& W{\hat{F}}_{\mu \nu }{W}^{-1}\\ {{\hat{F}}^{\prime }}_{\mu \nu }\equiv \text{gauge}\text{field}\text{strength}\text{in}V& \iff & {\hat{F}}_{\mu \nu }\equiv \text{gauge}\text{field}\text{strength}\text{in}U\end{array}\\ \downarrow \downarrow \\ {{\hat{F}}^{\prime }}_{\mu \nu }=\gamma \left(\frac{\left|a\right|}{a_0}\right){\hat{F}}_{\mu \nu }\end{array}$ (7)

where, $\gamma \left(\frac{\left|a\right|}{a_0}\right)$ can be defined as a function of $\frac{\left|a\right|}{a_0}$, which the author suggests is called the gauge similarity transformation factor. $\left|a\right|$ is the absolute value of the acceleration related to the region, which can be defined as the acceleration of the particle relative to a basic frame determined by distant matter in the universe. $a_0$ is an acceleration constant related to our galaxy. Experimental findings [18] show that it is approximately $a_0~c{H}_0=2\times {10}^{-8}\text{cm}\cdot {\text{s}}^{-2}~5\times {10}^{-8}\text{cm}\cdot {\text{s}}^{-2}$. Surprisingly, the value of $a_0$ of 2 × 10⁻⁸ cm·s⁻² is very consistent with the measured value of the acceleration of the solar system in the gravitational field of the Milky Way (2.32 ± 0.16) × 10⁻⁸ cm·s⁻² [19] , indicating that the regional scale corresponding to $a_0$ is the scale of our solar system. Therefore, the solar system can be used as a reference system to determine the scope of application of relativity and Newtonian mechanics. Generally speaking, as the regional scale corresponding to $a_0$ increases, it should be more difficult for the overall regional matter to change its acceleration $\left|a\right|$ under the gravitational field of the galaxy. So, on average, if the acceleration $a_0$ is taken as the benchmark, the larger the galaxy scale, the smaller its overall acceleration $\left|a\right|$ should be. For example, if the regional scale of the solar system corresponding to $a_0$ is $U_{sun}$ and the galaxy is $V_{uni}$, there should be an inverse relationship,

$\frac{\left|a\right|}{a_0}~\frac{U_{sun}}{V_{uni}}$ (8)

Obviously, when the galaxy scale $V_{uni}$ approaches the universe scale, the above formula tends to 0. If the galaxy scale shrinks to the $U_{sun}$ scale, the above formula is 1. When the galaxy scale is very smaller than the reference system scale, the above formula tends to $\infty$. The corresponding acceleration shows that when $\left|a\right|\ll a_0$, $\gamma \left(\frac{\left|a\right|}{a_0}\right)\to \frac{\left|a\right|}{a_0}\to 0$, and when $\left|a\right|=a_0$, $\gamma \left(\frac{\left|a\right|}{a_0}\right)\to \gamma \left(1\right)$, and when $\left|a\right|\gg a_0$, $\gamma \left(\frac{\left|a\right|}{a_0}\right)\to 1$, so the gauge similarity transformation factor $\gamma$ varies between 0 and 1. In order to reflect this law, we speculatively choose the $\gamma \left(\frac{\left|a\right|}{a_0}\right)\equiv th\left(\frac{\left|a\right|}{a_0}\right)$ function to “simulate” this change. However, this does not take into account the fact that Einstein’s law of gravity and Newton’s second law or the law of universal gravitation is valid within the scale of the solar system. Therefore, the author generally believes that when $\left|a\right|\ge a_0$, $th\left(\beta \frac{\left|a\right|}{a_0}\right)\to 1$ should hold. At this time, the reference system should include the scale of our solar system, and a suitable coefficient $\beta$ should be introduced to make $\left|a\right|\ge a_0$,

$th\left(\beta \frac{\left|a\right|}{a_0}\right)\to 1$ (9)

An important question now is whether the guage transformation corresponding to the $\gamma$ coefficient constructed in this way really exists, that is, the transformation formula,

${{\hat{F}}^{\prime }}_{\mu \nu }=W{\hat{F}}_{\mu \nu }{W}^{-1}=\gamma \left(\frac{\left|a\right|}{a_0}\right){\hat{F}}_{\mu \nu }$ (10)

Do the guage transformations (matrices or operators) $W$ and $W^{-1}$ in the equation exist reasonably? To this end, by multiplying both sides of Equation (10) by $W$, one can obtain:

$W{\hat{F}}_{\mu \nu }=\gamma \left(\frac{\left|a\right|}{a_0}\right){\hat{F}}_{\mu \nu }W$ (11)

Which can deduce:

$\left[W,{\hat{F}}_{\mu \nu }\right]=W{\hat{F}}_{\mu \nu }-{\hat{F}}_{\mu \nu }W=\gamma \left(\frac{\left|a\right|}{a_0}\right){\hat{F}}_{\mu \nu }W-{\hat{F}}_{\mu \nu }W=\left(\gamma \left(\frac{\left|a\right|}{a_0}\right)-1\right){\hat{F}}_{\mu \nu }W$ (12)

Then, except for the limit point $\gamma \left(\frac{\left|a\right|}{a_0}\right)=1$, $\left(\gamma \left(\frac{\left|a\right|}{a_0}\right)-1\right)\ne 0$, and in general ${\hat{F}}_{\mu \nu }W\ne 0$, because if it is 0, then $W{\hat{F}}_{\mu \nu }=0$ is obtained from Equation (11), and $\left[W,{\hat{F}}_{\mu \nu }\right]$ is meaningless under this condition. Therefore, it can be deduced that Equation (12) is not commutative under the condition of Equation (11):

$\left[W,{\hat{F}}_{\mu \nu }\right]\ne 0$ (13)

This is very logical, because $W$ has the properties of an operator and ${\hat{F}}_{\mu \nu }$ belongs to the gravitational field strength. If they are commutative, it means that the gravitational field strength given by the gauge similarity transformation is unchanged. This is possible in the solar system, but on a space-time scale greater than or equal to the galaxy, we now know that this is usually not true, otherwise, there would be no hypothesis of dark matter.

Furthermore, using Equation (10), ${{\hat{F}}^{\prime }}_{\mu \nu }=W{\hat{F}}_{\mu \nu }{W}^{-1}=\gamma {\hat{F}}_{\mu \nu }$, we can reversely deduce the matrix representation of the gauge transformation (conversion function) as follows:

First, consider that there are many matrices expressing the gravitational strength that is diagonal matrix of the 2nd-order metric tensor [1] [13] [14] , that is,

$\left\{{\hat{F}}_{\mu \nu }\right\}=\left\{g_{jj}\right\}\equiv \left(\begin{array}{c}{g}_{00}\\ 0\\ 0\\ 0\end{array}\begin{array}{c}0\\ g_{11}\\ 0\\ 0\end{array}\begin{array}{c}0\\ 0\\ g_{22}\\ 0\end{array}\begin{array}{c}0\\ 0\\ 0\\ g_{33}\end{array}\right)$ (14)

and the matrix expression of the gauge transformed gravitational strength is:

$\left\{{{\hat{F}}^{\prime }}_{\mu \nu }\right\}\equiv \gamma \left(\begin{array}{c}{g}_{00}\\ 0\\ 0\\ 0\end{array}\begin{array}{c}0\\ g_{11}\\ 0\\ 0\end{array}\begin{array}{c}0\\ 0\\ g_{22}\\ 0\end{array}\begin{array}{c}0\\ 0\\ 0\\ g_{33}\end{array}\right)$ (15)

This means that the gravitational strength is a diagonal matrix expression of the second-order metric tensor corresponding to a diagonal matrix $\left\{g_{jj}\right\}$ of $\left\{{\hat{F}}_{\mu \nu }\left(x\right)\right\}$. So we can always find such a matrix similarity transformation $\left\{W\right\}$, which transforms the diagonalization matrix of ${\hat{F}}_{\mu \nu }\left(x\right)$ into $\gamma \left\{g_{jj}\right\}$, and the similarity transformation expression is deduced as follows:

$\begin{array}{l}\left(\begin{array}{c}{w}_{00}\\ w_{10}\\ w_{20}\\ w_{30}\end{array}\begin{array}{c}{w}_{01}\\ w_{11}\\ w_{21}\\ w_{31}\end{array}\begin{array}{c}{w}_{02}\\ w_{12}\\ w_{22}\\ w_{32}\end{array}\begin{array}{c}{w}_{03}\\ w_{13}\\ w_{23}\\ w_{33}\end{array}\right)\left(\begin{array}{c}{g}_{00}\\ 0\\ 0\\ 0\end{array}\begin{array}{c}0\\ g_{11}\\ 0\\ 0\end{array}\begin{array}{c}0\\ 0\\ g_{22}\\ 0\end{array}\begin{array}{c}0\\ 0\\ 0\\ g_{33}\end{array}\right){\left(\begin{array}{c}{w}_{00}\\ w_{10}\\ w_{20}\\ w_{30}\end{array}\begin{array}{c}{w}_{01}\\ w_{11}\\ w_{21}\\ w_{31}\end{array}\begin{array}{c}{w}_{02}\\ w_{12}\\ w_{22}\\ w_{32}\end{array}\begin{array}{c}{w}_{03}\\ w_{13}\\ w_{23}\\ w_{33}\end{array}\right)}^{-1}\\ =\gamma \left(\begin{array}{c}{g}_{00}\\ 0\\ 0\\ 0\end{array}\begin{array}{c}0\\ g_{11}\\ 0\\ 0\end{array}\begin{array}{c}0\\ 0\\ g_{22}\\ 0\end{array}\begin{array}{c}0\\ 0\\ 0\\ g_{33}\end{array}\right)\end{array}$ (16)

where defining

$\left\{W\right\}\equiv \left(\begin{array}{c}{w}_{00}\\ w_{10}\\ w_{20}\\ w_{30}\end{array}\begin{array}{c}{w}_{01}\\ w_{11}\\ w_{21}\\ w_{31}\end{array}\begin{array}{c}{w}_{02}\\ w_{12}\\ w_{22}\\ w_{32}\end{array}\begin{array}{c}{w}_{03}\\ w_{13}\\ w_{23}\\ w_{33}\end{array}\right)=\left(W_0{W}_1{W}_2{W}_3\right)$ (17)

Then, one can get:

$\left(\begin{array}{c}{w}_{00}{g}_{00}\\ w_{10}{g}_{00}\\ w_{20}{g}_{00}\\ w_{30}{g}_{00}\end{array}\begin{array}{c}{w}_{01}{g}_{11}\\ w_{11}{g}_{11}\\ w_{21}{g}_{11}\\ w_{31}{g}_{11}\end{array}\begin{array}{c}{w}_{02}{g}_{22}\\ w_{12}{g}_{22}\\ w_{22}{g}_{22}\\ w_{32}{g}_{22}\end{array}\begin{array}{c}{w}_{03}{g}_{33}\\ w_{13}{g}_{33}\\ w_{23}{g}_{33}\\ w_{33}{g}_{33}\end{array}\right)=\gamma \left(\begin{array}{c}{g}_{00}{w}_{00}\\ g_{11}{w}_{10}\\ g_{22}{w}_{20}\\ g_{33}{w}_{30}\end{array}\begin{array}{c}{g}_{00}{w}_{01}\\ g_{11}{w}_{11}\\ g_{22}{w}_{21}\\ g_{33}{w}_{31}\end{array}\begin{array}{c}{g}_{00}{w}_{02}\\ g_{11}{w}_{12}\\ g_{22}{w}_{22}\\ g_{33}{w}_{32}\end{array}\begin{array}{c}{g}_{00}{w}_{03}\\ g_{11}{w}_{13}\\ g_{22}{w}_{23}\\ g_{33}{w}_{33}\end{array}\right)$ (18)

That is:

$\left\{{\hat{F}}_{\mu \nu }\right\}\left(W_0{W}_1{W}_2{W}_3\right)=\left(g_{00}{W}_0{g}_{11}{W}_1{g}_{22}{W}_2{g}_{33}{W}_3\right)$ (19)

The above formula can be simplified as:

$\left\{{\hat{F}}_{\mu \nu }\right\}{W}_i=g_{ii}{W}_i,i=0,1,2,3$ (20)

Now, we can associate $\left\{{\hat{F}}_{\mu \nu }\right\}$ with the eigenvalues and eigenvectors of its matrix, i.e. $g_{ii}$ is the i-th eigenvalue of the matrix $\left\{{\hat{F}}_{\mu \nu }\right\}$, and $W_i$ is the eigenvector corresponding to $g_{ii}$. Because $\left\{W\right\}$ must be a reversible matrix, the eigenvectors of $\left\{{\hat{F}}_{\mu \nu }\right\}$ need to be linearly independent, that is, the necessary and sufficient condition for the 4-order square matrix $\left\{{\hat{F}}_{\mu \nu }\right\}$ to be similar to the diagonal matrix $\left\{g_{jj}\right\}$ is that the eigenvectors $W_0,W_1,W_2$ and $W_3$ of $\left\{{\hat{F}}_{\mu \nu }\right\}$ are linearly independent of each other. Of course, under the condition of the Equation (16) this is satisfied, hence the linearly independent eigenvectors $W_0,W_1,W_2$ and $W_3$ of $\left\{{\hat{F}}_{\mu \nu }\right\}$ can be constructed. For example, from Equation (18), by selecting:

1) $w_{01}=1=w_{02}=w_{03}$, one can obtain:

$w_{00}=\frac{-g_{11}-g_{22}-g_{33}+3\gamma g_{00}}{\left(1-\gamma \right)g_{00}}$ (21)

2) $w_{10}=1=w_{12}=w_{13}$, one can obtain:

$w_{11}=\frac{-g_{00}-g_{22}-g_{33}+3\gamma g_{11}}{\left(1-\gamma \right)g_{11}}$ (22)

3) $w_{20}=1=w_{21}=w_{23}$, one can get:

$w_{22}=\frac{-g_{00}-g_{11}-g_{33}+3\gamma g_{22}}{\left(1-\gamma \right)g_{22}}$ (23)

4) $w_{30}=1=w_{31}=w_{32}$, one can get:

$w_{33}=\frac{-g_{00}-g_{11}-g_{22}+3\gamma g_{33}}{\left(1-\gamma \right)g_{33}}$ (24)

Finally, one can arrive at:

$W_0=\left(\begin{array}{c}\frac{-g_{11}-g_{22}-g_{33}+3\gamma g_{00}}{\left(1-\gamma \right)g_{00}}\\ 1\\ 1\\ 1\end{array}\right)$ (25)

$W_1=\left(\begin{array}{c}1\\ \frac{-g_{00}-g_{22}-g_{33}+3\gamma g_{11}}{\left(1-\gamma \right)g_{11}}\\ 1\\ 1\end{array}\right)$ (26)

$W_2=\left(\begin{array}{c}1\\ 1\\ \frac{-g_{00}-g_{11}-g_{33}+3\gamma g_{22}}{\left(1-\gamma \right)g_{22}}\\ 1\end{array}\right)$ (27)

$W_3=\left(\begin{array}{c}1\\ 1\\ 1\\ \frac{-g_{00}-g_{11}-g_{22}+3\gamma g_{33}}{\left(1-\gamma \right)g_{33}}\end{array}\right)$ (28)

The above results allow us to obtain a 4th-order gauge transformation matrix expression: $\left\{W\right\}=\left(W_0{W}_1{W}_2{W}_3\right)$, and this expression is non-degenerate, so it is reasonable, indicating that this gauge transformation that extends with the scale of the universe can reasonably exist, and is actually determined or selected by the intrinsic properties and structure of the principal and associated bundles of the universe. The $\gamma$ factor derived from this gauge transformation invariance is actually a function of the space-time scale. If the solar system is defined as a standard unit space-time scale, then $\gamma$ decreases as the space-time scale increases relative to the scale of the universal region in units of this standard scale, which can be simulated by the hyperbolic tangent, that is, the gauge similarity transformation equation can now be expressed as:

${{\hat{F}}^{\prime }}_{\mu \nu }=W{\hat{F}}_{\mu \nu }{W}^{-1}=th\left(\frac{\left|a\right|}{a_0}\right){\hat{F}}_{\mu \nu }$ (29)

Here, the region $U$ where ${\hat{F}}_{\mu \nu }$ is located is smaller or much smaller than the region $V$ where ${{\hat{F}}^{\prime }}_{\mu \nu }$ is located, or physically speaking, $U$ ~ Solar System, $V$ ≥ “Galaxy”. Therefore, it can be understood that the invariance of the gauge transformation caused by the choice of this kind of internal gauge in the universe finally determines $\gamma$, which is essentially determined by the internal space-time structure of the principal and associated fiber bundles of the universe. The most powerful evidence supporting this “theory of the determination of the intrinsic properties of the universe” is that the universe is indeed observed to be expanding at an accelerated rate in experiments, from which it can be obtained that a $\gamma$ factor does have an impact on Newton’s gravitational formula [18] [19] .

In order to ensure that Einstein’s law of gravitational attraction remains unchanged and Newton’s law of universal gravitation (and second law) remains unchanged within the scale of the solar system, the gauge similarity transformation should be given by:

${{\hat{F}}^{\prime }}_{\mu \nu }=W{\hat{F}}_{\mu \nu }{W}^{-1}={\hat{F}}_{\mu \nu }$ (30)

In the case of a galaxy greater than or equal to galaxies, for example, $\left|a\right|<a_0$, then $0<\gamma \left(\frac{\left|a\right|}{a_0}\right)<1$, and the gauge transformation above changes to:

${{\hat{F}}^{\prime }}_{\mu \nu }=\gamma \left(\frac{\left|a\right|}{a_0}\right){\hat{F}}_{\mu \nu }<{\hat{F}}_{\mu \nu }$ (31)

Thus, Newton’s second law from the perspective of our solar system becomes:

$m\gamma \left(\frac{\left|a\right|}{a_0}\right)a<ma=F$ (32)

Newton’s second law needs to be revised, which is consistent with MOND’s view [18] [20] - [22] .

In fact, from the point of view that the gauge field strength is the curvature of the manifold at the bottom of the principal fiber bundle, considering the linear Einstein equations under the Newtonian approximation, we can have [23] [24] :

$m\frac{{\text{d}}^2{x}^i}{{\left(\text{d}{x}^0\right)}^2}=-m{\Gamma }_{00}^i$ (33)

where, ${\Gamma }_{00}^i$ is the Christoffel symbol. According to Newton’s approximation, $g_{ab}={\eta }_{ab}+h_{ab}$, etc., we can get:

$m\frac{{\text{d}}^2{x}^i}{{\left(\text{d}{x}^0\right)}^2}=-m\frac{1}{2}{h}_{00,i}$ (34)

and considering in the classical gravitational field, Newton’s second law can be written as:

$m_I\frac{{\text{d}}^2{x}^i}{\text{d}{t}^2}=m_g\left(-\frac{\partial \phi }{\partial x^i}\right)$ (35)

where $\phi$ is the Newtonian gravitational potential. Comparing the above two equations, we get:

$\phi =\frac{1}{2}{h}_{00}+\text{const}$ (36)

Assume that the gravitational field at infinity disappears as $\phi =0$, then the metric should return to the Minkowski metric, $h_{00}=0$, so the constant in the above formula is 0, which allows one to obtain:

$g_{00}=-1+h_{00}=-\left(1+\frac{2\phi }{c^2}\right)$ (37)

However, considering the gauge transformation and Newton approximation:

${\hat{F}}_{\mu \nu }\to \gamma {\hat{F}}_{\mu \nu }\to \gamma {\hat{F}}_{\mu \nu }\to \gamma {\Gamma }_{00}^i\to \gamma g_{00}\to \gamma h_{00}\to \gamma \phi$ (38)

in the classical gravitational field above, Newton’s second law becomes:

$m\left(-\frac{\partial \gamma \phi }{\partial x^i}\right)=m\gamma a<ma=F$ (39)

This once again shows that the difference in gravity is reflected in the $\gamma$ factor, which is the origin of the dark matter hypothesis.

In addition, if the scale of cosmic space-time region $V_{uni}$ $\gg$ “reference system scale, or $\gg$ galaxy scale”, then $\frac{\left|a\right|}{a_0}\ll 1$, $th\left(\beta \frac{\left|a\right|}{a_0}\right)\to 0$, its physical meaning is that the gravitational gauge field strength is ${{\hat{F}}^{\prime }}_{\mu \nu }=th\left(\beta \frac{\left|a\right|}{a_0}\right){\hat{F}}_{\mu \nu }\to 0$, which corresponds to the gravitational gauge field strength of the manifold region on the cosmic scale is 0. Equation (39) means that Newton’s second law becomes:

$\underset{\frac{\left|a\right|}{a_0}\to 0}{\lim }\left[mth\left(\beta \frac{\left|a\right|}{a_0}\right)a\right]\to 0$ (40)

This means that the curvature of the entire region of the universe is zero, and the attraction disappears. The author speculates that at this scale, the universe may be driven by the energy conservation law of the initial kinetic energy of the Big Bang and the disappearance of the attraction, and it may also produce repulsive force. This part of the conversion energy of the disappearance of gravity plus the initial kinetic energy left by the original Big Bang may be the origin of the so-called dark energy [1] .

The connection or curvature of large-scale spacetime can certainly affect the dynamical equations. In fact, the physical essence reflected by the gauge similarity transformation Formula (4) can be realized by introducing a factor that varies with the space-time scale into the cosmological term of the Einstein equations, for example:

$G_{ab}+\Lambda \left(\frac{\left|a\right|}{a_0}\right)g_{ab}=8\pi T_{ab}$ (41)

where $G_{ab}$ is the Einstein tensor, $g_{ab}$ is the metric tensor, $T_{ab}$ is the energy-momentum tensor, Λ was originally the cosmological constant, and is now related to $\gamma \left(\frac{\left|a\right|}{a_0}\right)$, that is,

$-\Lambda g_{ab}=8\pi \left(\frac{1-\gamma }{\gamma }\right)T_{ab}$ (42)

For a detailed proof, please refer to [1] . In fact, according to the similarity gauge transformation and Einstein equation, let $\gamma =th\left(\beta \frac{\left|a\right|}{a_0}\right)$, we have:

$R_{ab}-\frac{1}{2}R{g}_{ab}=\frac{8\pi }{\gamma }{T}_{ab}$ (43)

Further, consider the Einstein equation with the cosmological constant Λ:

$R_{ab}-\frac{1}{2}R{g}_{ab}+\Lambda g_{ab}=8\pi T_{ab}$ (44)

can be transformed into:

$R_{ab}=8\pi \left(T_{ab}-\frac{1}{2}T{g}_{ab}\right)+\Lambda g_{ab}$ (45)

with

$R=-8\pi T+4\Lambda$ (46)

then by substituting three Equations (44)-(46) into Equation (43), one can obtain:

$\frac{8\pi }{\gamma }{T}_{ab}=R_{ab}-\frac{1}{2}R{g}_{ab}=8\pi T_{ab}-\Lambda g_{ab}$ (47)

So, we have Equation (42), and under this condition, Equation (43) and Einstein’s equation with cosmological constant (41) can be derived from each other.

The gauge similarity transformation $R_{ab}\to \gamma R_{ab}$ (note that according to the principal bundle theory, curvature is equivalent to field strength) corresponding to condition (42) is equivalent to adding the cosmological constant Λ to the Einstein equations. This shows that in a sense, adding a certain cosmological constant to the Einstein equations is equivalent to applying a certain gauge transformation to its curvature. The cosmological constant is no longer a constant and is related to the coefficient $\gamma$ of the gauge field strength after the gauge similarity transformation. This is the physical meaning of the cosmological constant. Moreover,Λ is affected by $\gamma \left(\frac{\left|a\right|}{a_0}\right)$, and Λ is also a function of the space-time scale, that is, it becomes very large as the space-time scale increases, which naturally solves the “cosmological constant problem of physicists or astronomers”, that is, the cosmological constant increases to a very large extent as the universal region scale increases [15] ; when the space-time scale tends to the universe scale, $\gamma \left(\frac{\left|a\right|}{a_0}\right)\to \frac{\left|a\right|}{a_0}\to 0$, $\left|\Lambda \right|\to \infty$, the gauge gravitational field strength is 0, and the Einstein equation enters the chaotic region; when the space-time scale is smaller than the solar system scale, $\gamma \left(\frac{\left|a\right|}{a_0}\right)\to 1$, $\Lambda \to 0$, the gravitational field strength remains unchanged after the gauge transformation, the original Einstein equation is held, and Newton’s second law is held. The conditions that must be satisfied for the specific change of Λ are determined by Equation (42).

In short, from the above analysis of the cosmological constant of the Einstein equation, the author found that since the cosmological constant now changes with the factor $\gamma \left(\frac{\left|a\right|}{a_0}\right)$ and is no longer a constant, but a function that changes with the scale of the cosmic region, or is related to the regional acceleration ratio $\frac{\left|a\right|}{a_0}$, the Einstein equation will change with the size of the cosmic region and present a different form, so that the Newton’s law or law of universal gravitation observed at the same point in the bottom manifold (that is, our universe) in the $U$ region may not be the same as the Newton’s law or law of universal gravitation observed at the same point in the $V$ region. The transformation between them satisfies the similar gauge transformation (or generalized gauge equation); and it is the difference between both that creates the fictitious origin of dark matter and dark energy.

In fact, from Equation (41), the Einstein equation containing the cosmological term becomes:

$R_{ab}-\frac{1}{2}R{g}_{ab}=8\pi T_{ab}-\Lambda g_{ab}=8\pi \frac{T_{ab}}{\gamma \left(\frac{\left|a\right|}{a_0}\right)}$ (48)

that is, it can be abbreviated as:

$G_{ab}=8\pi \frac{T_{ab}}{\gamma \left(\frac{\left|a\right|}{a_0}\right)}$ (49)

Here, $G_{ab}$ certainly is the Einstein tensor. So, the Einstein equations for the universe on a large scale are suggested to be given by:

$G_{ab}=R_{ab}-\frac{1}{2}R{g}_{ab}+\Lambda \left(\frac{\left|a\right|}{a_0}\right)g_{ab}=8\pi T_{ab}$ (50)

Therefore, the evolution law of Einstein’s equations in large-scale cosmic regions is affected by $\Lambda \left(\frac{\left|a\right|}{a_0}\right)$ or $\gamma \left(\frac{\left|a\right|}{a_0}\right)$ and is expressed as follows:

$\{\begin{array}{l}{G}_{ab}=8\pi T_{ab},\left|a\right|\ll a_0,\gamma \to 1,\text{Einstein}\text{equation},\text{Newton}\text{gravity}\text{valid}\text{region}\\ G_{ab}=8\pi \frac{T_{ab}}{\gamma },0<\gamma <1,\text{Einstein},\text{Newton}\text{gravity}\text{deviation}\text{region,}\text{vertual}\text{dark}\text{matter}\uparrow \\ -\Lambda g_{ab}=8\pi \left(\frac{1-\gamma }{\gamma }\right)T_{ab}\to \infty ,\left|a\right|\gg a_0,\gamma \to 0,\text{Chaotic}\text{region}\text{of}\text{gravity}\downarrow \text{or}\text{repulsion}\uparrow \end{array}$(51)

In detail, the mechanism of the above picture at the two extreme points of curvature 0 and 1 can be explained as follows:

1) The 0 curvature interval is a fluctuation zone or a chaotic zone where gravity disappears completely or repulsion increases. At the 0 limit point, Einstein’s equations and Newton’s gravitational law do not hold. The repulsion energy converted from the disappearance of gravity is the source of dark energy.

2) Whether the universe will continue to evolve depends on which part of the force disappears more. According to the principle of energy conservation, the author speculates that: from the first branch of the simulation function $\gamma =th\left(\beta \frac{\left|a\right|}{a_0}\right)$ to point 0, the gravitational force disappears more, the repulsive force will be generated more, and the universe may continue to evolve towards the second branch (see details in ref. [1] ). On the contrary, it will rebound and shrink towards the first branch with a high probability. It is not impossible that the universe will repeatedly explode and shrink in the first branch with a certain probability to oscillate and evolve. The driving force of evolution ultimately comes from the principle of energy conservation. In the process of the first branch, universe curvature tending to 1, as $\left|a\right|\ll a_0$, the universe appears as a closed 4-dimensional spherical surface, and its space-time scale is gradually compressed until it becomes a singularity. Then the Big Bang begins a new curvature evolution from $1\to 0$; at the limit point +1, Einstein’s equation and Newton’s law of universal gravitation are both valid. During the first branch: $0<\gamma <1$, Einstein’s equation and Newton’s law of universal gravitation are only “partially” valid compared with that in the solar system. The difference between the solar system and beyond galaxy on the Einstein equation and Newton’s law of universal gravitation is just the source of the illusory dark matter. At this time, Newton’s law of universal gravitation is corrected to: $F=G\frac{mM}{\gamma r^2}$, and the Einstein field equation is corrected to: $G_{ab}=8\pi \frac{T_{ab}}{\gamma }$.

Obviously, $\gamma$ makes Newton’s gravitational constant no longer a constant. Can we say that $\gamma$ is ultimately determined by dark matter or dark energy? This statement seems to be divorced from the $\gamma$ factor, the gauge gravitational potential or the gauge field strength, which is essentially a reflection of the principal bundle connection or curvature. At present, it has not been truly measured whether dark matter and dark energy really exist, let alone that both dark matter and dark energy are some kinds of gauge field. In fact, the existence of $\gamma$ can be simply explained by the fact that the gravitational constant is no longer a constant but a function of the evolution of the scale of the universal region. This function becomes smaller as the scale of the universal region increases in units of the solar system, that is,

$F=\frac{G}{\gamma }\frac{mM}{r^2}=G^{\prime }\frac{mM}{r^2}\to G^{\prime }\equiv \frac{G}{\gamma }$ (52)

Therefore, the above formula once again clearly draws the most important conclusion of this article: the gravitational constant $G\to G^{\prime }$ is not a constant but a function that changes with the size of the universal region! The produced reason for the dark matter and dark energy hypothesis may be because the gravitational constant is a variable!

Based on the cosmic evolution picture constructed by the author using the principal and associated fiber bundles, we systematically proposed the $\gamma$ factor theory that the Newton’s law of universal gravitation and the cosmological constant of Einstein’s equation must be corrected in the large-scale space-time structure of the universe with the help of gauge invariance. From this, the cyclic hypothesis of cosmic evolution is proposed, and further finds that the gravitational constant that people think is natural is not a constant but a function of the scale of the universal region! Therefore, the produced reason for the dark matter and dark energy hypothesis may be that the gravitational constant is a variable. The existence of dark matter and dark energy may be just an illusory hypothesis, and their origin comes from the understanding that the gravitational constant is constant.

This study is supported by the Shenzhen Xinzhou Biological Information Technology Co., Ltd.