The Models of Primary Particles

1. Introduction

We know intimately the term "atom" which comes from ancient Greek and means "uncuttable" or translated as "indivisible." In the early 19th century, the scientist John Dalton introduced the modern definition of an atom to characterize chemical elements. It was discovered that Dalton's atoms were not actually indivisible about a century later. An atom is made up of three basic types of subatomic particles: electrons, protons, and neutrons, which occupy the tiny space inside an atom. Protons and neutrons form the nucleus, which contains most of the mass of an atom. Electrons are the lightest charged particles in nature and orbit the nucleus of an atom. An electron is seemingly indivisible yet. So far, we have not split an electron into two or more smaller particles. We just annihilate the positive and negative electrons. A free neutron is unstable and decays into a proton, an electron and an electron neutrino. A free proton, on the other hand, is stable, and consists of two up quarks and one down quark in the modern Standard Model. Another question is whether a quark can be cut into smaller parts or whether the matter is infinitely divisible.

This paper tries to answer the above questions from a different perspective. What happens when light or electromagnetic, gravitational and other waves make up the primary particles that are the fundamental elements of matter? The assumptions are then derived:

a.

Light or electromagnetic waves, the weak interaction, gravity, and the strong interaction are independent waves without rest mass. However, their structures are different.

b.

Light and gravity can be described by the wave Equation with the field strength Ê and speed c.

$${\nabla }^2Ê={\displaystyle \frac{1}{c^2}}{\displaystyle \frac{{\partial }^2Ê}{\partial t^2}}$$

(1)

• c. The weak and strong interactions can be described by the 4-dimensional Laplace Equation with field strength Ê' and speed c'.

$${\nabla }^2{Ê}^{\text{'}}=-{\displaystyle \frac{1}{c{\text{'}}^2}}{\displaystyle \frac{{\partial }^2{Ê}^{\text{'}}}{\partial t^2}}$$

(2)

• d. According to the electroweak theory, light and the weak interaction have the same speed c_L with spin number +1 or –1.
• e. Gravity and the strong interaction have the same speed c_G without spin, and c_G is constant in a vacuum.
• f. The primary particles, which are electrons, electron neutrinos, and dark neutrinos in this paper, are made of the four waves mentioned above.

The Formation of Primary Particles

Suppose that the birth of primary particles can be divided into the following two parts:

• Light and the weak interaction couple together (hereafter referred to as the E-W couple) when they have the same spin number and the second-order partial derivatives of their fields with respect to time ${\displaystyle \frac{{\partial }^2Ê}{\partial t^2}}$ are equal. Gravity and the strong interaction also couple together (hereafter referred to as the G-S couple) when the second-order partial derivatives of their fields with respect to time ${\displaystyle \frac{{\partial }^2{Ê}^{\text{'}}}{\partial t^2}}$ are equal. So we have

$$c_L^2{\nabla }^2({Ê}_e+{Ê}_w)=0$$

(3)

and

$$c_G^2{\nabla }^2({Ê}_G+{Ê}_S)=0$$

(4)

Where Ê_e, Ê_w, Ê_G, and Ê_S are electric, weak interaction, gravitational, and strong interaction fields.

• b. The original spins of the light and the weak interaction convert the polarity of the electric and weak charges when an E-W couple is formed.
• c. It makes a primary particle when two coupled waves attract each other and shrink to a tiny sphere. One E-W couple and one G-S couple produce an electron or a positron whose charge property depends on the original spin of the E-W couple. Dark neutrinos are composed of two G-S couples. Two E-W couples with different original spin compress themselves into an electron neutrino. But they cannot attract each other if they have the same original spin.

Thus, in the spherical coordinate system, we have the uniform Laplace Equation for Equations (3), (4)

$$c^2\left[{\displaystyle \frac{1}{r^2}}{\displaystyle \frac{\partial }{\partial r}}\left(r^2{\displaystyle \frac{\partial Ê}{\partial r}}\right)+{\displaystyle \frac{1}{r^2\sin \theta }}{\displaystyle \frac{\partial }{\partial \theta }}\left(\sin \theta {\displaystyle \frac{\partial Ê}{\partial \theta }}\right)+{\displaystyle \frac{1}{r^2{\sin }^2\theta }}{\displaystyle \frac{{\partial }^2Ê}{\partial {\phi }^2}}\right]=0$$

(5)

2. The Fields and Binding Energies

We can now derive the electric, gravitational, weak interaction, and strong interaction fields Ê based on Equation (5) and existing physical laws and data. Since Equation (5) has infinite pairs of conjugate particular solutions[1]

$$\begin{array}{c}Ê={\displaystyle \frac{1}{c^2}}\left(A_j{r}^j+{\displaystyle \frac{B_j}{r^{l+1}}}\right){\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )e^{-ik\phi }}\\ {Ê}^{*}={\displaystyle \frac{1}{c^2}}\left(A_j{r}^j+{\displaystyle \frac{B_j}{r^{j+1}}}\right){\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )e^{ik\phi }}\end{array}$$

(6)

where Aj and Bj are constants, ${\mathrm{P}}_j^k(\cos \theta )$ are associated Legendre polynomials, j and k are integers, j = 0, 1, 2, 3, …, k ≤ j, and j is called the degree of associated Legendre polynomials.

Clearly, there are

$${[{\widehat{E}}_a+{\widehat{E}}_b]}^{*}={\widehat{E}}_a^{*}+{\widehat{E}}_b^{*}$$

(7)

and

$${({\widehat{E}}^{*})}^{*}=\widehat{E}$$

(8)

where the subscripts a and b denote electric, gravitational, weak interaction, or strong interaction.

The field Ê may be split into the macroscopic term ${\displaystyle \frac{1}{c^2}}\left(A_j{r}^j+{\displaystyle \frac{B_j}{r^{j+1}}}\right)$ and the quantum factors ${\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )e^{-ik\phi }}$, ${\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )e^{ik\phi }}$. The degree of associated Legendre polynomials j rules properties of the field Ê because it is not only an exponent of r in the macroscopic item but also impacts forms of the quantum factors.

Comparing Equation (6) to Gauss's law of electrostatics and Newton's law of gravity, the electric field Ê_e is

$$\begin{array}{c}{\widehat{E}}_e=\pm {\displaystyle \frac{1}{c_L^2}}{\displaystyle \frac{B_1}{r^2}}{\displaystyle \sum _{k=0}^{-1}{\mathrm{P}}_1^k(\cos \theta )e^{-ik\phi }}=\pm {\displaystyle \frac{\widehat{q}}{c_L^2{r}^2}}(\cos \theta +\sin \theta e^{-i\phi })\\ {\widehat{E}}_e^{*}=\pm {\displaystyle \frac{1}{c_L^2}}{\displaystyle \frac{B_1}{r^2}}{\displaystyle \sum _{k=0}^1{\mathrm{P}}_1^k(\cos \theta )e^{ik\phi }}=\pm {\displaystyle \frac{\widehat{q}}{c_L^2{r}^2}}(\cos \theta +\sin \theta e^{i\phi })\end{array}$$

(9)

and the gravitational field Ê_G is

$$\begin{array}{c}{\widehat{E}}_G=-{\displaystyle \frac{1}{c_G^2}}{\displaystyle \frac{B_1}{r^2}}{\displaystyle \sum _{k=0}^1{\mathrm{P}}_1^k(\cos \theta )e^{-ik\phi }}=-{\displaystyle \frac{\widehat{m}}{c_G^2{r}^2}}(\cos \theta +\sin \theta e^{-i\phi })\\ {\widehat{E}}_G^{*}=-{\displaystyle \frac{1}{c_G^2}}{\displaystyle \frac{B_1}{r^2}}{\displaystyle \sum _{k=0}^1{\mathrm{P}}_1^k(\cos \theta )e^{ik\phi }}=-{\displaystyle \frac{\widehat{m}}{c_G^2{r}^2}}(\cos \theta +\sin \theta e^{i\phi })\end{array}$$

(10)

where '–' means attractive interaction, and '+' means repulsive interaction, as usual in this paper, $\widehat{q}$is the mathematical electric charge, and $\widehat{m}$is the mathematical mass or the gravitational charge.

The weak interaction has an intensity of a similar magnitude to the electromagnetic force at very short distances (around 10^-18 meters), but this starts to decrease exponentially with increasing distance. Its effective range is about 10^-17 to 10^-16 meters^{[2, 3]}. All of the above can help us to determine the weak interaction field Ê_w from Equation (6). Ê_w is therefore supposed to be equal to

$$\begin{array}{c}{\widehat{E}}_w=\left\{\begin{array}{l}\pm {\displaystyle \frac{1}{c_L^2}}{\displaystyle \frac{B_1}{r^2}}{\displaystyle \sum _{k=0}^1{\mathrm{P}}_1^k(\cos \theta )e^{-ik\phi }}=\pm {\displaystyle \frac{\widehat{w}}{c_L^2{r}^2}}(\cos \theta +\sin \theta e^{-i\phi })r\le R_{cw}\hfill \\ \pm {\displaystyle \frac{1}{c_L^2}}{\displaystyle \frac{B_m}{r^{m+1}}}{\displaystyle \sum _{k=0}^m{\mathrm{P}}_m^k(\cos \theta )e^{-ik\phi }}=\pm {\displaystyle \frac{\widehat{w}{R}_{cw}^{m-1}}{c_L^2{r}^{m+1}}}{\displaystyle \sum _{k=0}^m{\mathrm{P}}_m^k(\cos \theta )e^{-ik\phi }}r>R_{cw}\hfill \end{array}\right.\\ {\widehat{E}}_w^{*}=\left\{\begin{array}{l}\pm {\displaystyle \frac{1}{c_L^2}}{\displaystyle \frac{B_1}{r^2}}{\displaystyle \sum _{k=0}^1{\mathrm{P}}_1^k(\cos \theta )e^{ik\phi }}=\pm {\displaystyle \frac{\widehat{w}}{c_L^2{r}^2}}(\cos \theta +\sin \theta e^{i\phi })r\le R_{cw}\hfill \\ \pm {\displaystyle \frac{1}{c_L^2}}{\displaystyle \frac{B_m}{r^{m+1}}}{\displaystyle \sum _{k=0}^m{\mathrm{P}}_m^k(\cos \theta )e^{ik\phi }}=\pm {\displaystyle \frac{\widehat{w}{R}_{cw}^{m-1}}{c_L^2{r}^{m+1}}}{\displaystyle \sum _{k=0}^m{\mathrm{P}}_m^k(\cos \theta )e^{ik\phi }}r>R_{cw}\hfill \end{array}\right.\end{array}$$

(11)

where$\widehat{w}$is the mathematical weak charge, R_cw is the critical radius of the weak interaction, and m is an integer greater than 1.

The strong force is a short-range interaction (around 10^-15 meters) similar to the weak force. But its range is more complex than that of the weak force. At distances comparable to the diameter of a proton, it is approximately 100 times as strong as the electromagnetic force. At smaller distances, however, it becomes weaker. In particle physics, this effect is known as asymptotic freedom^{[4, 5]}. Furthermore, it is supposed that the fields of the four fundamental forces have a unified form in a very small region. Hence, Equation (6) can be translated into the strong interaction field Ê_S

$$\begin{array}{c}{\widehat{E}}_S=\left\{\begin{array}{l}-{\displaystyle \frac{B_1}{c_G^2{r}^2}}{\displaystyle \sum _{k=0}^1{\mathrm{P}}_1^k(\cos \theta )e^{-ik\phi }}=-{\displaystyle \frac{\widehat{s}}{c_G^2{r}^2}}(\cos \theta +\sin \theta e^{-i\phi })r\le R_{cS1}\hfill \\ -{\displaystyle \frac{A_{1}r}{c_G^2}}{\displaystyle \sum _{k=0}^1{\mathrm{P}}_1^k(\cos \theta )e^{-ik\phi }}=-{\displaystyle \frac{\widehat{s}r}{c_G^2{R}_{cS1}^3}}(\cos \theta +\sin \theta e^{-i\phi })R_{cS1}<r\le R_{cS2}\hfill \\ -{\displaystyle \frac{1}{c_G^2}}{\displaystyle \frac{B_n}{r^{n+1}}}{\displaystyle \sum _{k=0}^n{\mathrm{P}}_n^k(\cos \theta )e^{-ik\phi }}=-{\displaystyle \frac{\widehat{s}{R}_{cS2}^{n+2}}{c_G^2{R}_{cS1}^3{r}^{n+1}}}{\displaystyle \sum _{k=0}^n{\mathrm{P}}_n^k(\cos \theta )e^{-ik\phi }}r>R_{cS2}\hfill \end{array}\right.\\ {\widehat{E}}_S^{*}=\left\{\begin{array}{l}-{\displaystyle \frac{B_1}{c_G^2{r}^2}}{\displaystyle \sum _{k=0}^1{\mathrm{P}}_1^k(\cos \theta )e^{ik\phi }}=-{\displaystyle \frac{\widehat{s}}{c_G^2{r}^2}}(\cos \theta +\sin \theta e^{i\phi })r\le R_{cS1}\hfill \\ -{\displaystyle \frac{A_{1}r}{c_G^2}}{\displaystyle \sum _{k=0}^1{\mathrm{P}}_1^k(\cos \theta )e^{ik\phi }}=-{\displaystyle \frac{\widehat{s}r}{c_G^2{R}_{cS1}^3}}(\cos \theta +\sin \theta e^{i\phi })R_{cS1}<r\le R_{cS2}\hfill \\ -{\displaystyle \frac{1}{c_G^2}}{\displaystyle \frac{B_n}{r^{n+1}}}{\displaystyle \sum _{k=0}^n{\mathrm{P}}_n^k(\cos \theta )e^{ik\phi }}=-{\displaystyle \frac{\widehat{s}{R}_{cS2}^{n+2}}{c_G^2{R}_{cS1}^3{r}^{n+1}}}{\displaystyle \sum _{k=0}^n{\mathrm{P}}_n^k(\cos \theta )e^{ik\phi }}r>R_{cS2}\hfill \end{array}\right.\end{array}$$

(12)

where$\widehat{s}$is the mathematical strong charge, R_cS1 and R_cS2 are the 1^st and 2^nd critical radii of the strong interaction, and n is an integer greater than 1.

Further, it is assumed that R_cS1 << R_cw < R_cS2 and m ≠ n.

Now we turn to determine the energy. The energy density of the wave Equation (1) is given by ${\left|\widehat{E}\right|}^2$, and the Equations (1) and (2) are Lorentz invariant. So we hope that the Equation of energy is also Lorentz invariant, and let the binding energy E_a^b of the four fundamental fields be

$$E_{a^b}\propto {\displaystyle \frac{1}{t_2-t_1}}{\displaystyle \underset{V}{\iiint }{\displaystyle \underset{c{t}_1}{\overset{c{t}_2}{\int }}\left|{\widehat{E}}_a{\widehat{E}}_b^{*}\right|dxdydzd(ct)}}$$

(13)

Note that ${\widehat{E}}_a$ and ${\widehat{E}}_b^{*}$ are independent of time, and there are

$${\displaystyle \underset{0}{\overset{2\pi }{\int }}{e}^{ik\phi }{e}^{-i{k}^{\prime }\phi }d\phi }=2\pi {\delta }_{k{k}^{\prime }}$$

(14)

and

$${\displaystyle \underset{0}{\overset{\pi }{\int }}{\mathrm{P}}_j^k(\cos \theta ){\mathrm{P}}_{j^{\prime }}^k(\cos \theta )sin\theta d\theta }={\displaystyle \frac{(j+k)!}{(j-k)!}}{\displaystyle \frac{2}{2j+1}}{\delta }_{j{j}^{\prime }}={\mathrm{N}}_j^k{\delta }_{j{j}^{\prime }}$$

(15)

When we use light as a measurement medium to determine the energy as usual, we can translate Equation (13) into the spherical coordinate system with the optical medium, which is

$$\begin{array}{c}{E}_{a^b}={\displaystyle \frac{1}{(t_2-t_1)2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{R}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}{\displaystyle \underset{c_L{t}_1}{\overset{c_L{t}_2}{\int }}\left|{\widehat{E}}_a{\widehat{E}}_b^{*}\right|r^2\sin \theta drd\theta d\phi d(c_{L}t)}}}}\\ ={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{R}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left|{\widehat{E}}_a{\widehat{E}}_b^{*}\right|r^2\sin \theta drd\theta d\phi }}}\end{array}$$

(16)

where R is the radius at which two fields begin to interact with each other.

The general energy expression can be calculated when Equation (6) is substituted into Equation (16).

$$\begin{array}{c}{E}_{a^b}={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{R}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left|{\widehat{E}}_a{\widehat{E}}_b^{*}\right|r^2\sin \theta drd\theta d\phi }}}\\ ={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{R}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left[\frac{1}{c_a^2}\left|A_{aj}{r}^j+\frac{B_{aj}}{r^{j+1}}\right|{\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )e^{-ik\phi }}\right.}}}\\ \left.\times {\displaystyle \frac{1}{c_b^2}}\left|A_{bj}{r}^j+{\displaystyle \frac{B_{bj}}{r^{j+1}}}\right|{\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )e^{ik\phi }}\right]r^2\sin \theta drd\theta d\phi \\ ={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{R}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left[\frac{1}{c_b^2}\left|A_{bj}{r}^j+\frac{B_{bj}}{r^{j+1}}\right|{\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )e^{-ik\phi }}\right.}}}\\ \left.\times {\displaystyle \frac{1}{c_a^2}}\left|A_{aj}{r}^j+{\displaystyle \frac{B_{aj}}{r^{j+1}}}\right|{\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )e^{ik\phi }}\right]r^2\sin \theta drd\theta d\phi \\ ={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{R}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left|{\widehat{E}}_b{\widehat{E}}_a^{*}\right|r^2\sin \theta drd\theta d\phi }}}\\ =E_{b^a}\\ ={\displaystyle \frac{c_L^2}{c_a^2{c}_b^2}}{\displaystyle \underset{R}{\overset{\infty }{\int }}\left|A_{aj}{r}^j+\frac{B_{aj}}{r^{j+1}}\right|\left|A_{bj}{r}^j+\frac{B_{bj}}{r^{j+1}}\right|r^{2}dr}\end{array}$$

(17)

Based on Equation (17) and associated with Equations (9) to (12), we first compute the self-binding energy of the four fields. The self-binding energy of the four fields E_e^e, E_G^G, E_w^w, and E_S^S are

$$E_{e^e}=c_L^2{\displaystyle \underset{r}{\overset{\infty }{\int }}\frac{{\widehat{q}}^2}{c_L^4{r}^2}dr}={\displaystyle \frac{{\widehat{q}}^2}{c_L^{2}r}}$$

(18)

$$E_{G^G}=c_L^2{\displaystyle \underset{r}{\overset{\infty }{\int }}\frac{{\widehat{m}}^2}{c_G^4{r}^2}dr}={\displaystyle \frac{c_L^2{\widehat{m}}^2}{c_G^{4}r}}$$

(19)

$$E_{w^w}=\left\{\begin{array}{l}{c}_L^2{\displaystyle \underset{r}{\overset{R_{cw}}{\int }}\frac{{\widehat{w}}^2}{c_L^4{r}^2}dr}+c_L^2{\displaystyle \underset{R_{cw}}{\overset{\infty }{\int }}\frac{{\widehat{w}}^2{R}_{cw}^{2m-2}}{c_L^4{r}^{2m}}dr}={\displaystyle \frac{{\widehat{w}}^2}{c_L^2}}\left[{\displaystyle \frac{1}{r}}-{\displaystyle \frac{2m-2}{(2m-1)R_{cw}}}\right]r\le R_{cw}\hfill \\ c_L^2{\displaystyle \underset{r}{\overset{\infty }{\int }}\frac{{\widehat{w}}^2{R}_{cw}^{2m-2}}{c_L^4{r}^{2m}}dr}={\displaystyle \frac{{\widehat{w}}^2{R}_{cw}^{2m-2}}{c_L^2(2m-1)r^{2m-1}}}r>R_{cw}\hfill \end{array}\right.$$

(20)

and

$$E_{S^S}=\left\{\begin{array}{l}\begin{array}{l}{c}_L^2{\displaystyle \underset{r}{\overset{R_{cS1}}{\int }}\frac{{\widehat{s}}^2}{c_G^4{r}^2}dr}+c_L^2{\displaystyle \underset{R_{cS1}}{\overset{R_{cS2}}{\int }}\frac{{\widehat{s}}^2{r}^4}{c_G^4{R}_{cS1}^6}dr}+c_L^2{\displaystyle \underset{R_{cS2}}{\overset{\infty }{\int }}\frac{{\widehat{s}}^2{R}_{cS2}^{2n+4}}{c_G^4{R}_{cS1}^6{r}^{2n}}dr}\\ ={\displaystyle \frac{c_L^2{\widehat{s}}^2}{c_G^4}}\left[{\displaystyle \frac{1}{r}}-{\displaystyle \frac{6}{5R_{cS1}}}+{\displaystyle \frac{2(n+2)R_{cS2}^5}{5(2n-1)R_{cS1}^6}}\right]r\le R_{cS1}\end{array}\hfill \\ \begin{array}{l}{c}_L^2{\displaystyle \underset{r}{\overset{R_{cS2}}{\int }}\frac{{\widehat{s}}^2{r}^4}{c_G^4{R}_{cS1}^6}dr}+c_L^2{\displaystyle \underset{R_{cS2}}{\overset{\infty }{\int }}\frac{{\widehat{s}}^2{R}_{cS2}^{2n+4}}{c_G^4{R}_{cS1}^6{r}^{2n}}dr}\\ ={\displaystyle \frac{c_L^2{\widehat{s}}^2}{c_G^4}}\left[{\displaystyle \frac{1}{5R_{cS1}^6}}\left(R_{cS2}^5-r^5\right)+{\displaystyle \frac{R_{cS2}^5}{(2n-1)R_{cS1}^6}}\right]R_{cS1}<r\le R_{cS2}\end{array}\hfill \\ c_L^2{\displaystyle \underset{r}{\overset{\infty }{\int }}\frac{{\widehat{s}}^2{R}_{cS2}^{2n+4}}{c_G^4{R}_{cS1}^6{r}^{2n}}dr}={\displaystyle \frac{c_L^2{\widehat{s}}^2}{c_G^4}}\times {\displaystyle \frac{R_{cS2}^{2n+4}}{(2n-1)R_{cS1}^6{r}^{2n-1}}}r>R_{cS2}\hfill \end{array}\right.$$

(21)

The binding energy of the four fields, such as E_e^G, E_e^w, E_e^S, etc. are

$$E_{e^G}=c_L^2{\displaystyle \underset{r}{\overset{\infty }{\int }}\frac{\widehat{q}}{c_L^2{r}^2}\frac{\widehat{m}}{c_G^2{r}^2}{r}^{2}dr}={\displaystyle \frac{\widehat{q}\widehat{m}}{c_G^{2}r}}$$

(22)

$$E_{e^w}=\left\{\begin{array}{l}{c}_L^2{\displaystyle \underset{r}{\overset{R_{cw}}{\int }}\frac{\widehat{q}}{c_L^2{r}^2}\frac{\widehat{w}}{c_L^2{r}^2}{r}^{2}dr}+0={\displaystyle \frac{\widehat{q}\widehat{w}}{c_L^2}}\left({\displaystyle \frac{1}{r}}-{\displaystyle \frac{1}{R_{cw}}}\right)r\le R_{cw}\hfill \\ 0r>R_{cw}\hfill \end{array}\right.$$

(23)

$$E_{e^S}=\left\{\begin{array}{l}{c}_L^2{\displaystyle \underset{r}{\overset{R_{cS1}}{\int }}\frac{\widehat{q}}{c_L^2{r}^2}\frac{\widehat{s}}{c_G^2{r}^2}{r}^{2}dr}+c_L^2{\displaystyle \underset{R_{cS1}}{\overset{R_{cS2}}{\int }}\frac{\widehat{q}}{c_L^2{r}^2}\frac{\widehat{s}r}{c_G^2{R}_{cS1}^3}{r}^{2}dr}+0={\displaystyle \frac{\widehat{q}\widehat{s}}{c_G^2}}\left({\displaystyle \frac{1}{r}}-{\displaystyle \frac{3}{2R_{cS1}}}+{\displaystyle \frac{R_{cS2}^2}{2R_{cS1}^3}}\right)r\le R_{cS1}\hfill \\ c_L^2{\displaystyle \underset{r}{\overset{R_{cS2}}{\int }}\frac{\widehat{q}}{c_L^2{r}^2}\frac{\widehat{s}r}{c_G^2{R}_{cS1}^3}{r}^{2}dr}+0={\displaystyle \frac{\widehat{q}\widehat{s}}{2c_G^2{R}_{cS1}^3}}\left(R_{cS2}^2-r^2\right)R_{cS1}<r\le R_{cS2}\hfill \\ 0r>R_{cS2}\hfill \end{array}\right.$$

(24)

$$E_{G^w}=\left\{\begin{array}{l}{c}_L^2{\displaystyle \underset{r}{\overset{R_{cw}}{\int }}\frac{\widehat{m}}{c_G^2{r}^2}\frac{\widehat{w}}{c_L^2{r}^2}{r}^{2}dr}+0={\displaystyle \frac{\widehat{m}\widehat{w}}{c_G^2}}\left({\displaystyle \frac{1}{r}}-{\displaystyle \frac{1}{R_{cw}}}\right)r\le R_{cw}\hfill \\ 0r>R_{cw}\hfill \end{array}\right.$$

(25)

$$E_{G^S}=\left\{\begin{array}{l}{c}_L^2{\displaystyle \underset{r}{\overset{R_{cS1}}{\int }}\frac{\widehat{m}}{c_G^2{r}^2}\frac{\widehat{s}}{c_G^2{r}^2}{r}^{2}dr}+c_L^2{\displaystyle \underset{R_{cS1}}{\overset{R_{cS2}}{\int }}\frac{\widehat{m}}{c_G^2{r}^2}\frac{\widehat{s}r}{c_G^2{R}_{cS1}^3}{r}^{2}dr}+0={\displaystyle \frac{c_L^2\widehat{m}\widehat{s}}{c_G^4}}\left({\displaystyle \frac{1}{r}}-{\displaystyle \frac{3}{2R_{cS1}}}+{\displaystyle \frac{R_{cS2}^2}{2R_{cS1}^3}}\right)r\le R_{cS1}\hfill \\ c_L^2{\displaystyle \underset{r}{\overset{R_{cS2}}{\int }}\frac{\widehat{m}}{c_G^2{r}^2}\frac{\widehat{s}r}{c_G^2{R}_{cS1}^3}{r}^{2}dr}+0={\displaystyle \frac{c_L^2\widehat{m}\widehat{s}}{2c_G^4{R}_{cS1}^3}}\left(R_{cS2}^2-r^2\right)R_{cS1}<r\le R_{cS2}\hfill \\ 0r>R_{cS2}\hfill \end{array}\right.$$

(26)

and

$$E_{w^S}=\left\{\begin{array}{l}{c}_L^2{\displaystyle \underset{r}{\overset{R_{cS1}}{\int }}\frac{\widehat{w}}{c_L^2{r}^2}\frac{\widehat{s}}{c_G^2{r}^2}{r}^{2}dr}+c_L^2{\displaystyle \underset{R_{cS1}}{\overset{R_{cw}}{\int }}\frac{\widehat{w}}{c_L^2{r}^2}\frac{\widehat{s}r}{c_G^2{R}_{cS1}^3}{r}^{2}dr}+0={\displaystyle \frac{\widehat{w}\widehat{s}}{c_G^2}}\left({\displaystyle \frac{1}{r}}-{\displaystyle \frac{3}{2R_{cS1}}}+{\displaystyle \frac{R_{cw}^2}{2R_{cS1}^3}}\right)r\le R_{cS1}\hfill \\ c_L^2{\displaystyle \underset{r}{\overset{R_{cw}}{\int }}\frac{\widehat{w}}{c_L^2{r}^2}\frac{\widehat{s}r}{c_G^2{R}_{cS1}^3}{r}^{2}dr}+0={\displaystyle \frac{\widehat{w}\widehat{s}}{2c_G^2{R}_{cS1}^3}}\left(R_{cw}^2-r^2\right)R_{cS1}<r\le R_{cw}\hfill \\ 0r>R_{cw}\hfill \end{array}\right.$$

(27)

For the convenience of subsequent calculations, we require the results of our defined energy to be consistent with those of the conventional physical method. Comparing Equations (18) and (19) with the electric and gravitational potential energy formulas, it is easy to see that the relations of the mathematical electric charge $\widehat{q}$ and the mathematical mass $\widehat{m}$ to the electric charge q and the mass m are as follows

$$\widehat{q}=c_L\sqrt{k}q$$

(28)

and

$$\widehat{m}={\displaystyle \frac{c_G^2\sqrt{G}m}{c_L}}$$

(29)

where k is the Coulomb constant, and G is the gravitational constant.

The Structures of Primary Particles

A primary particle looks like a tiny spheroidal balloon with two envelopes (Figure 1). Each of them is formed by an E-W couple or a G-S couple. The envelopes can characterize as:

• The whole binding energy of the coupled waves concentrates on the envelopes.
• The macroscopic items of combined field strengths of the two coupled waves are equal on the envelopes. Outside the envelopes, the coupled waves become two independent static fields. But there are no fields inside the envelopes.
• The size of the envelope, that is the size of a primary particle, depends on the critical radius of the weak or strong interaction.
• The two envelopes have the same inherent frequency ν_in, although this is not mathematically required.
• The degree of the associated Legendre polynomials j is the same on the two envelopes.
• The behaviors of the two envelopes obey the Self-Conjugate Mechanism, which requires that one occupies the surface of ${\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )e^{-ik\phi }}$ and the other must take up ${\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )e^{ik\phi }}$, or they are conjugate to each other.

Hence, the binding energy of a primary particle E_pri can be generally described as

$$\begin{array}{c}{E}_{pri}={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{R}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left|[{\widehat{E}}_a+{\widehat{E}}_b]{[{\widehat{E}}_c+{\widehat{E}}_d]}^{*}\right|r^2\sin \theta drd\theta d\phi }}}\\ ={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{R}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left|{\widehat{E}}_a{\widehat{E}}_c^{*}+{\widehat{E}}_a{\widehat{E}}_d^{*}+{\widehat{E}}_b{\widehat{E}}_c^{*}+{\widehat{E}}_b{\widehat{E}}_d^{*}\right|r^2\sin \theta drd\theta d\phi }}}\\ =E_{a^c}+E_{a^d}+E_{b^c}+E_{b^d}\end{array}$$

(30)

based on Equation (16), where the subscripts a to d denote the electric, gravitational, weak interaction, or strong interaction.

E_pri is clearly equivalent to the rest mass. The total energy of a primary particle comprises the binding energy or the rest mass and the energy in static fields.

2.1. An Electron Neutrino

An electron neutrino is composed of two E-W couples with different original spin. In order to explore its structure, these assumptions should be adopted:

• Its radius r_{e_ν} is equal to the critical radius of the weak interaction R_cw.
• The charges in Equations (9) and (11) are equal and minimal for an electron neutrino, i.e., ${\widehat{q}}_{e\_\nu }={\widehat{w}}_{e\_\nu }$, if ${\widehat{q}}_{e\_\nu }$ and ${\widehat{w}}_{e\_\nu }$ are the mathematical electric charge and the mathematical weak charge of an electron neutrino.

Integrating Equations (9), (11), and the above characters, we have two field Equations on envelopes of an electron neutrino Ê_{E, e_ν}

$${\widehat{E}}_{E,e\_\nu }\left\{\begin{array}{l}{\left[{\widehat{E}}_e+{\widehat{E}}_w\right]}_{r=R_{cw}}={\displaystyle \frac{{\widehat{q}}_{e\_\nu }+{\widehat{w}}_{e\_\nu }}{c_L^2{R}_{cw}^2}}(\cos \theta +\sin \theta e^{-i\phi })\hfill \\ -{\left[{\widehat{E}}_e+{\widehat{E}}_w\right]}_{r=R_{cw}}^{*}=-{\displaystyle \frac{{\widehat{q}}_{e\_\nu }+{\widehat{w}}_{e\_\nu }}{c_L^2{R}_{cw}^2}}(\cos \theta +\sin \theta e^{i\phi })\hfill \end{array}\right.$$

(31)

where '–' only indicates that two E-W couples are attracted to each other on the envelopes.

Based on Equation (30) and associated Equations (18), (20), and (23), we can easily compute the binding energy of an electron neutrino E_{e_ν}.

$$\begin{array}{c}{E}_{e\_\nu }={\left[E_{e^e}+2E_{e^w}+E_{w^w}\right]}_{r=R_{cw}}\\ ={\left[{\displaystyle \frac{{\widehat{q}}_{e\_\nu }^2}{c_L^{2}r}}+{\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{w}}_{e\_\nu }}{c_L^2}}\left({\displaystyle \frac{1}{r}}-{\displaystyle \frac{1}{R_{cw}}}\right)+{\displaystyle \frac{{\widehat{w}}_{e\_\nu }^2}{c_L^{2}r}}\left({\displaystyle \frac{1}{r}}-{\displaystyle \frac{2m-2}{(2m-1)R_{cw}}}\right)\right]}_{r=R_{cw}}\\ ={\displaystyle \frac{1}{c_L^2{R}_{cw}}}\left({\widehat{q}}_{e\_\nu }^2+{\displaystyle \frac{{\widehat{w}}_{e\_\nu }^2}{2m-1}}\right)={\displaystyle \frac{2m}{(2m-1)R_{cw}}}{\displaystyle \frac{{\widehat{q}}_{e\_\nu }^2}{c_L^2}}={\displaystyle \frac{2m}{(2m-1)R_{cw}}}{\displaystyle \frac{{\widehat{w}}_{e\_\nu }^2}{c_L^2}}\end{array}$$

(32)

Combining the envelopes' characters b. with Equations (9), (11), and (31), the fields around an electron neutrino Ê_{e_ν} can be directly written as

$${\left.{\widehat{E}}_{e\_\nu }\right|}_{r>R_{cw}}\left\{\begin{array}{l}\begin{array}{l}\left[{\widehat{E}}_{e-e\_\nu }+{\widehat{E}}_{w-e\_\nu }\right]=\pm {\displaystyle \frac{{\widehat{q}}_{e\_\nu }}{c_L^2{r}^2}}(\cos \theta +\sin \theta e^{-i\phi })\pm {\displaystyle \frac{{\widehat{w}}_{e\_\nu }{R}_{cw}^{m-1}}{c_L^2{r}^{m+1}}}{\displaystyle \sum _{k=0}^m{\mathrm{P}}_m^k(\cos \theta )e^{-ik\phi }}\\ {\left[{\widehat{E}}_{e-e\_\nu }+{\widehat{E}}_{w-e\_\nu }\right]}^{*}=\mp {\displaystyle \frac{{\widehat{q}}_{e\_\nu }}{c_L^2{r}^2}}(\cos \theta +\sin \theta e^{i\phi })\mp {\displaystyle \frac{{\widehat{w}}_{e\_\nu }{R}_{cw}^{m-1}}{c_L^2{r}^{m+1}}}{\displaystyle \sum _{k=0}^m{\mathrm{P}}_m^k(\cos \theta )e^{ik\phi }}\end{array}\hfill \end{array}\right.r>R_{cw}$$

(33)

2.1. Dark Neutrinos

Two G-S couples make a dark neutrino. However, the strong interaction field has two critical radii, so there are two types of dark neutrinos, and they are named Dark I and Dark II. Similar to Section 4.1, it is assumed that:

The sizes of Dark I and II are equal to the 1^st and 2^nd critical radii of the strong interaction.

Dark I and II have the same mathematical mass and the same mathematical strong charge.

The mathematical strong charge ${\widehat{s}}_{D\_\nu }$ is minimal and equal to the mathematical mass ${\widehat{m}}_{D\_\nu }$, i.e., ${\widehat{m}}_{D\_\nu }={\widehat{s}}_{D\_\nu }$.

Replicating the process of the previous section, we have the fields of a Dark I on the envelopes Ê_{E, D_νI}

$${\widehat{E}}_{E,D\_\nu I}\left\{\begin{array}{l}{\left[{\widehat{E}}_G+{\widehat{E}}_S\right]}_{r=R_{cS1}}={\displaystyle \frac{{\widehat{m}}_{D\_\nu }+{\widehat{s}}_{D\_\nu }}{c_G^2{R}_{cS1}^2}}(\cos \theta +\sin \theta e^{-i\phi })\hfill \\ -{\left[{\widehat{E}}_G+{\widehat{E}}_S\right]}_{r=R_{cS1}}^{*}=-{\displaystyle \frac{{\widehat{m}}_{D\_\nu }+{\widehat{s}}_{D\_\nu }}{c_G^2{R}_{cS1}^2}}(\cos \theta +\sin \theta e^{i\phi })\hfill \end{array}\right.$$

(34)

where '–' only means that two G-S couples are attracted to each other on the envelopes.

Based on Equation (30) and associated Equations (19), (21), and (26), we can compute the binding energy of a Dark I E_{D_νI}

$$\begin{array}{c}{E}_{D\_\nu I}={\left[E_{G^G}+2E_{G^S}+E_{S^S}\right]}_{r=R_{cS1}}\\ ={\displaystyle \frac{c_L^2}{c_G^4}}{\left[{\displaystyle \frac{{\widehat{m}}_{D\_\nu }^2}{r}}+{\widehat{m}}_{D\_\nu }{\widehat{s}}_{D\_\nu }\left({\displaystyle \frac{1}{r}}-{\displaystyle \frac{3}{2R_{cS1}}}+{\displaystyle \frac{R_{cS2}^2}{2R_{cS1}^3}}\right)+{\widehat{s}}_{D\_\nu }^2\left({\displaystyle \frac{1}{r}}-{\displaystyle \frac{6}{5R_{cS1}}}+{\displaystyle \frac{2(n+2)R_{cS2}^5}{5(2n-1)R_{cS1}^6}}\right)\right]}_{r=R_{cS1}}\\ \approx {\displaystyle \frac{2(n+2)R_{cS2}^5}{5(2n-1)R_{cS1}^6}}{\displaystyle \frac{c_L^2{\widehat{s}}_{D\_\nu }^2}{c_G^4}}={\displaystyle \frac{2(n+2)R_{cS2}^5}{5(2n-1)R_{cS1}^6}}{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4}}\end{array}$$

(35)

To get the fields of a Dark II on the envelopes Ê_{E, D_νII} and the binding energy of a Dark II E_{D_νII}, we imitate the last process and have

$${\widehat{E}}_{E,D\_\nu II}\left\{\begin{array}{l}{\left[{\widehat{E}}_G+{\widehat{E}}_S\right]}_{r=R_{cS2}}={\displaystyle \frac{1}{c_G^2}}\left({\displaystyle \frac{{\widehat{m}}_{D\_\nu }}{R_{cS2}^2}}+{\displaystyle \frac{R_{cS2}}{R_{cS1}^3}}{\widehat{s}}_{D\_\nu }\right)(\cos \theta +\sin \theta e^{-i\phi })\hfill \\ -{\left[{\widehat{E}}_G+{\widehat{E}}_S\right]}_{r=R_{cS2}}^{*}=-{\displaystyle \frac{1}{c_G^2}}\left({\displaystyle \frac{{\widehat{m}}_{D\_\nu }}{R_{cS2}^2}}+{\displaystyle \frac{R_{cS2}}{R_{cS1}^3}}{\widehat{s}}_{D\_\nu }\right)(\cos \theta +\sin \theta e^{i\phi })\hfill \end{array}\right.$$

(36)

and

$$\begin{array}{c}{E}_{D\_\nu II}={\left[E_{G-G}+2E_{G-S}+E_{S-S}\right]}_{r=R_{cS2}}\\ ={\displaystyle \frac{c_L^2}{c_G^4}}{\left[{\displaystyle \frac{{\widehat{m}}_{D\_\nu }^2}{r}}+{\displaystyle \frac{{\widehat{s}}_{D\_\nu }^2}{5R_{cS1}^6}}\left(R_{cS2}^5-r^5\right)+{\displaystyle \frac{R_{cS2}^5{\widehat{s}}_{D\_\nu }^2}{(2n-1)R_{cS1}^6}}\right]}_{r=R_{cS2}}\\ \approx {\displaystyle \frac{R_{cS2}^5}{(2n-1)R_{cS1}^6}}{\displaystyle \frac{c_L^2{\widehat{s}}_{D\_\nu }^2}{c_G^4}}={\displaystyle \frac{R_{cS2}^5}{(2n-1)R_{cS1}^6}}{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4}}\\ \approx {\displaystyle \frac{5}{2(n+2)}}{E}_{D\_\nu I}\end{array}$$

(37)

Following the computations of the particle external field in the previous section, we can obtain the fields around a Dark I Ê_{D_νI}

$${\left.{\widehat{E}}_{D\_\nu I}\right|}_{r>R_{cS1}}\left\{\begin{array}{l}\begin{array}{c}\left[{\widehat{E}}_{G-D\_\nu I}+{\widehat{E}}_{S-D\_\nu I}\right]=\left\{\begin{array}{l}{\widehat{E}}_{G-D\_\nu I}-{\displaystyle \frac{{\widehat{s}}_{D\_\nu }r}{c_G^2{R}_{cS1}^3}}(\cos \theta +\sin \theta e^{-i\phi })R_{cS1}<r\le R_{cS2}\hfill \\ {\widehat{E}}_{G-D\_\nu I}-{\displaystyle \frac{{\widehat{s}}_{D\_\nu }{R}_{cS2}^{n+2}}{c_G^2{R}_{cS1}^3{r}^{n+1}}}{\displaystyle \sum _{k=0}^n{\mathrm{P}}_n^k(\cos \theta )e^{-ik\phi }}r>R_{cS2}\hfill \\ \mathrm{where}{\widehat{E}}_{G-D\_\nu I}=-{\displaystyle \frac{{\widehat{m}}_{D\_\nu }}{c_G^2{r}^2}}(\cos \theta +\sin \theta e^{-i\phi })\hfill \end{array}\right.\end{array}\hfill \\ {\left[{\widehat{E}}_{G-D\_\nu I}+{\widehat{E}}_{S-D\_\nu I}\right]}^{*}=\left\{\begin{array}{l}{\widehat{E}}_{G-D\_\nu I}^{*}-{\displaystyle \frac{{\widehat{s}}_{D\_\nu }r}{c_G^2{R}_{cS1}^3}}(\cos \theta +\sin \theta e^{i\phi })R_{cS1}<r\le R_{cS2}\hfill \\ {\widehat{E}}_{G-D\_\nu I}^{*}-{\displaystyle \frac{{\widehat{s}}_{D\_\nu }{R}_{cS2}^{n+2}}{c_G^2{R}_{cS1}^3{r}^{n+1}}}{\displaystyle \sum _{k=0}^n{\mathrm{P}}_n^k(\cos \theta )e^{ik\phi }}r>R_{cS2}\hfill \\ \mathrm{where}{\widehat{E}}_{G-D\_\nu I}^{*}=-{\displaystyle \frac{{\widehat{m}}_{D\_\nu }}{c_G^2{R}_{cS1}^2}}(\cos \theta +\sin \theta e^{i\phi })\hfill \end{array}\right.\hfill \end{array}\right.$$

(38)

and the fields around a Dark II Ê_{D_νII}

$${\left.{\widehat{E}}_{D\_\nu II}\right|}_{r>R_{cS2}}\left\{\begin{array}{l}\begin{array}{c}\left[{\widehat{E}}_{G-D\_\nu II}+{\widehat{E}}_{S-D\_\nu II}\right]=-{\displaystyle \frac{{\widehat{m}}_{D\_\nu }}{c_G^2{r}^2}}(\cos \theta +\sin \theta e^{-i\phi })\\ -{\displaystyle \frac{{\widehat{s}}_{D\_\nu }{R}_{cS2}^{n+2}}{c_G^2{R}_{cS1}^3{r}^{n+1}}}{\displaystyle \sum _{k=0}^n{\mathrm{P}}_n^k(\cos \theta )e^{-ik\phi }}\\ {\left[{\widehat{E}}_{G-D\_\nu II}+{\widehat{E}}_{S-D\_\nu II}\right]}^{*}=-{\displaystyle \frac{{\widehat{m}}_{D\_\nu }}{c_G^2{R}_{cS1}^2}}(\cos \theta +\sin \theta e^{i\phi })\\ -{\displaystyle \frac{{\widehat{s}}_{D\_\nu }{R}_{cS2}^{n+2}}{c_G^2{R}_{cS1}^3{r}^{n+1}}}{\displaystyle \sum _{k=0}^n{\mathrm{P}}_n^k(\cos \theta )e^{ik\phi }}\end{array}\hfill \end{array}\right.r>R_{cS2}$$

(39)

Comparing Equations (35) with (37) reveals that, as far as measurements, namely energy, are concerned, there is little difference between Dark I and Dark II. However, their volumes are significant differences in the microscopic domain. There should only be Dark IIs in most cases following the principle of energy minimization.

2.3. An Electron or A Positron

Electrons and positrons have the same structure. We will not distinguish significantly between electrons and positrons in the subsequent descriptions and computations. One E-W couple and one G-S couple attract each other to form an electron or a positron, so its structure is the most complex in primary particles. Following the assumptions about dark neutrinos, it is supposed that:

The radius of an electron r_e equals the critical radius of the weak interaction R_cw, although there are three critical radii for the weak and strong interactions.

The mathematical electric charge ${\widehat{q}}_e$ and the mathematical weak charge ${\widehat{w}}_e$ are equal, i.e., ${\widehat{q}}_e={\widehat{w}}_e$.

The mathematical strong charge ${\widehat{s}}_e$ is minimal, i.e., ${\widehat{s}}_e={\widehat{s}}_{D\_\nu }$.

Referring to the way we did in the previous sections, we can obtain the field of an electron on the envelopes Ê_{E, e} and the binding energy of an electron E_e.

$${\widehat{E}}_{E,e}\left\{\begin{array}{l}{\left[{\widehat{E}}_e+{\widehat{E}}_w\right]}_{r=R_{cw}}={\displaystyle \frac{{\widehat{q}}_e+{\widehat{w}}_e}{c_L^2{R}_{cw}^2}}(\cos \theta +\sin \theta e^{-i\phi })\hfill \\ -{\left[{\widehat{E}}_G+{\widehat{E}}_S\right]}_{r=R_{cw}}^{*}=-{\displaystyle \frac{1}{c_G^2}}\left({\displaystyle \frac{{\widehat{m}}_e}{R_{cw}^2}}+{\displaystyle \frac{{\widehat{s}}_e{R}_{cw}}{R_{cS1}^3}}\right)(\cos \theta +\sin \theta e^{i\phi })\hfill \\ {\displaystyle \frac{{\widehat{q}}_e+{\widehat{w}}_e}{c_L^2{R}_{cw}^2}}={\displaystyle \frac{1}{c_G^2}}\left({\displaystyle \frac{{\widehat{m}}_e}{R_{cw}^2}}+{\displaystyle \frac{{\widehat{s}}_e{R}_{cw}}{R_{cS1}^3}}\right)\hfill \end{array}\right.$$

(40)

and

$$\begin{array}{c}{E}_e={\left[E_{e^G}+E_{e^S}+E_{G^w}+E_{w^S}\right]}_{r=R_{cw}}\\ ={\left[{\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_e}{c_G^{2}r}}+{\displaystyle \frac{{\widehat{q}}_e{\widehat{s}}_e}{2c_G^2{R}_{cS1}^3}}\left(R_{cS2}^2-r^2\right)\right]}_{r=R_{cw}}\\ ={\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_e}{c_G^2{R}_{cw}}}+{\displaystyle \frac{{\widehat{q}}_e{\widehat{s}}_e}{2c_G^2{R}_{cS1}^3}}\left(R_{cS2}^2-R_{cw}^2\right)\end{array}$$

(41)

According to the previous assumptions R_cS1 << R_cw < R_cS2, ${\widehat{q}}_e={\widehat{w}}_e$, and ${\widehat{s}}_e={\widehat{s}}_{D\_\nu }={\widehat{m}}_{D\_\nu }$, we have

$${\displaystyle \frac{2{\widehat{q}}_e}{c_L^2{R}_{cw}^2}}={\displaystyle \frac{2{\widehat{w}}_e}{c_L^2{R}_{cw}^2}}\approx {\displaystyle \frac{{\widehat{s}}_{D\_\nu }{R}_{cw}}{c_G^2{R}_{cS1}^3}}={\displaystyle \frac{{\widehat{m}}_{D\_\nu }{R}_{cw}}{c_G^2{R}_{cS1}^3}}$$

(42)

and

$$E_e\approx {\displaystyle \frac{{\widehat{q}}_e{\widehat{s}}_{D\_\nu }}{2c_G^2{R}_{cS1}^3}}\left(R_{cS2}^2-R_{cw}^2\right)={\displaystyle \frac{{\widehat{q}}_e^2}{c_L^2{R}_{cw}}}\left({\displaystyle \frac{R_{cS2}^2}{R_{cw}^2}}-1\right)$$

(43)

Now we directly give the result for fields around an electron Ê_e.

$${\left.{\widehat{E}}_e\right|}_{r>R_{cw}}\left\{\begin{array}{l}\begin{array}{l}\left[{\widehat{E}}_{e-e}+{\widehat{E}}_{w-e}\right]=\pm {\displaystyle \frac{{\widehat{q}}_e}{c_L^2{r}^2}}(\cos \theta +\sin \theta e^{-i\varphi })\pm {\displaystyle \frac{{\widehat{w}}_e}{c_L^2{r}^2}}(\cos \theta +\sin \theta e^{-i\phi })r>R_{cw}\hfill \end{array}\hfill \\ {\left[{\widehat{E}}_{G-e}+{\widehat{E}}_{S-e}\right]}^{*}=\left\{\begin{array}{l}-{\displaystyle \frac{{\widehat{m}}_e}{c_G^2{r}^2}}(\cos \theta +\sin \theta e^{i\varphi })-{\displaystyle \frac{{\widehat{s}}_{e}r}{c_G^2{R}_{cS1}^3}}(\cos \theta +\sin \theta e^{i\phi })R_{cw}<r\le R_{cS2}\hfill \\ -{\displaystyle \frac{{\widehat{m}}_e}{c_G^2{r}^2}}(\cos \theta +\sin \theta e^{i\varphi })-{\displaystyle \frac{{\widehat{s}}_{D\_\nu }{R}_{cS2}^{n+2}}{c_G^2{R}_{cS1}^3{r}^{n+1}}}{\displaystyle \sum _{k=0}^n{\mathrm{P}}_n^k(\cos \theta )e^{ik\phi }}r>R_{cS2}\hfill \end{array}\right.\hfill \end{array}\right.$$

(44)

Next, the examination of Equations (40) and (41) reveals that the 2^nd critical radius of the strong interaction R_cS2 should be the geometric characterization parameter of the G-S envelope rather than the 1^st critical radius of the strong interaction R_cS1. It is further assumed that R_cw and R_cS2 are proportional to the wavelengths of the E-W and the G-S couples, respectively, i.e., $R_{cw}=\xi {\lambda }_{E-\mathrm{W}}$, $R_{cS2}=\xi {\lambda }_{G-S}$. Thus, we can write directly with the envelopes characters d.

$${\displaystyle \frac{R_{cS2}}{R_{cw}}}={\displaystyle \frac{{\lambda }_{G-S}}{{\lambda }_{E-W}}}={\displaystyle \frac{c_G}{c_L}}$$

(45)

which shows that the speed of gravity c_G is faster than the speed of light c_L when the above Equation compares with Equation (41).

3. The Interactions Between Two Primary Particles

Imagine that two static primary particles are initially rested on each other and then separated by a repulsive force over a distance of l. Particles I and II have potential and kinetic energy in the separated state (Figure 2). In the initial state, particles I and II are equivalent in that they share a common center of the sphere (Figure 2) because there are no fields inside the envelopes of the primary particles.

According to the law of conservation of energy, the initial external energy of particles I and II, or the total energy in their external fields, is equal to the sum of potential energy and non-potential energy (including kinetic energy, magnetic energy for electrons, etc.) after their separation. Referring to Equation (16), the potential energy of the two primary particles $E_{P_1^P_2}$ should be associated with the following Equation

$$\begin{array}{c}{E}_{P_1^P_2}\propto {\displaystyle \frac{1}{(t_2-t_1)2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}\times \\ \left({\displaystyle \underset{R_c}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}{\displaystyle \underset{c_L{t}_1}{\overset{c_L{t}_2}{\int }}\left|{\widehat{E}}_{P_1}{\widehat{E}}_{P_2}^{*}\right|r^2\sin \theta drd\theta d\phi d(c_{L}t)}}}}-{\displaystyle \underset{R_c}{\overset{l}{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}{\displaystyle \underset{c_L{t}_1}{\overset{c_L{t}_2}{\int }}\left|{\widehat{E}}_{P_1}{\widehat{E}}_{P_2}^{*}\right|r^2\sin \theta drd\theta d\phi d(c_{L}t)}}}}\right)\\ ={\displaystyle \frac{1}{(t_2-t_1)2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{l}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}{\displaystyle \underset{c_L{t}_1}{\overset{c_L{t}_2}{\int }}\left|{\widehat{E}}_{P_1}{\widehat{E}}_{P_2}^{*}\right|r^2\sin \theta drd\theta d\phi d(c_{L}t)}}}}\\ ={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{l}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left|{\widehat{E}}_{P_1}{\widehat{E}}_{P_2}^{*}\right|r^2\sin \theta drd\theta d\phi d(c_{L}t)}}}l\mathrm{is}\mathrm{constant}\end{array}$$

(46)

where ${\widehat{E}}_{P_1}$ and ${\widehat{E}}_{P_2}$ are the external fields of primary particles I and II, l is the distance between the two particles (Figure 2), R_c is the larger of the two particles’ radii, and ${\displaystyle \underset{R_c}{\overset{l}{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}{\displaystyle \underset{c_L{t}_1}{\overset{c_L{t}_2}{\int }}\left|{\widehat{E}}_{P_1}{\widehat{E}}_{P_2}^{*}\right|r^2\sin \theta drd\theta d\phi d(c_{L}t)}}}}$ converts to the non-potential energy.

From Equations (18) and (19), the integration result of Equation (46) in the macroscopic world is consistent with the expressions for electric potential energy and gravitational potential energy. Then Equation (46) must still hold in the microscopic world. It is therefore assumed that the Self-Conjugate Mechanism between two interacting primary particles remains everywhere. The Self-Conjugate Mechanism makes two sets of fields around one particle conjugate to two sets of fields around another particle, depending on the rotation of two particles. In other words, two primary particles have achieved the self-conjugacy after rotating one cycle with their respective angular velocities ω₁ and ω₂ (Figure 2), and ${\omega }_1=\pm {\omega }_2=\omega$. There are self-conjugate lines in pairs, and it is further presumed that in their respective spherical coordinate systems. If two primary particles are repulsive to each other (Figure 2), one is at $I:(r_c,\theta ,\pi +{\omega }_1{t}_c)$ for particle I and $II:(r_c,\theta ,\pi -{\omega }_2{t}_c)$ for particle II, and another is at $I:(r_c,\theta ,\pi -{\omega }_1{t}_c)$ for particle I and $II:(r_c,\theta ,\pi +{\omega }_2{t}_c)$ for particle II. If two primary particles are attractive, then r_I and r_II in Figure 2 are reversed, so that the upper (left) line in Figure 2 becomes $I:(r_c,\theta ,{\omega }_1{t}_c)$ for particle I and $II:(r_c,\theta ,2\pi -{\omega }_2{t}_c)$ for particle II, and the lower (right) line becomes $I:(r_c,\theta ,2\pi -{\omega }_1{t}_c)$ for particle I and $II:(r_c,\theta ,{\omega }_2{t}_c)$ for particle II. Obviously, these lines are in the opposite position meanwhile r_I = r_II = r, θ_I = θ_II = θ (Figure 2).

Notice that there are

$$\begin{array}{c}{\displaystyle \underset{2\eta \pi }{\overset{2\eta \pi +2\pi }{\int }}{e}^{ik(\pi +{\omega }_1{t}_c)}{e}^{-i{k}^{\prime }(\pi -{\omega }_2{t}_c)}d({\omega }_1{t}_c)}=\left\{\begin{array}{l}{\displaystyle \underset{2\eta \pi }{\overset{2\eta \pi +2\pi }{\int }}{e}^{i(k-k^{\prime })\pi }{e}^{i(k+k^{\prime }){\omega }_1{t}_c}d({\omega }_1{t}_c)}{\omega }_2={\omega }_1\hfill \\ {\displaystyle \underset{2\eta \pi }{\overset{2\eta \pi +2\pi }{\int }}{e}^{i(k-k^{\prime })\pi }{e}^{i(k-k^{\prime }){\omega }_1{t}_c}d({\omega }_1{t}_c)}{\omega }_2=-{\omega }_1\hfill \end{array}\right.\\ =\left\{\begin{array}{l}\begin{array}{l}0\mathrm{Same}\mathrm{rotational}\mathrm{direction}\\ \end{array}\hfill \\ 2\pi {\delta }_{k{k}^{\prime }}\mathrm{Opposite}\mathrm{rotational}\mathrm{directions}\hfill \end{array}\right.\end{array}$$

(47)

and

$$\begin{array}{c}{\displaystyle \underset{2\eta \pi }{\overset{2\eta \pi +2\pi }{\int }}{e}^{ik({\omega }_1{t}_c)}{e}^{-i{k}^{\prime }(2\pi -{\omega }_2{t}_c)}d({\omega }_1{t}_c)}={\displaystyle \underset{2\eta \pi }{\overset{2\eta \pi +2\pi }{\int }}{e}^{ik({\omega }_1{t}_c)}{e}^{i{k}^{\prime }({\omega }_2{t}_c)}d({\omega }_1{t}_c)}\\ =\left\{\begin{array}{l}\begin{array}{l}0\mathrm{Same}\mathrm{rotational}\mathrm{direction}\\ \end{array}\hfill \\ 2\pi {\delta }_{k{k}^{\prime }}\mathrm{Opposite}\mathrm{rotational}\mathrm{directions}\hfill \end{array}\right.\end{array}$$

(48)

where η = 0, 1, 2, …, η and t_c is the time it takes two primary particles to complete one rotation.

Thus, at a self-conjugate line, taking the r_I coordinate system of the particle I as a reference (Figure 2), and using Equations (6), (47) and (48), Equation (46) can be rewritten as

$$\begin{array}{c}{E}_{P_1^P_2}={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{l}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{2\eta \pi }{\overset{2\eta \pi +2\pi }{\int }}\left|{\widehat{E}}_{P_1}(r,\theta ,{\xi }_1){\widehat{E}}_{P_2}^{*}(r,\theta ,{\xi }_2)\right|r^2\sin \theta drd\theta d({\omega }_1{t}_c)}}}\\ ={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{l}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{2\eta \pi }{\overset{2\eta \pi +2\pi }{\int }}\left[{\widehat{E}}_{P_1}(r){\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )e^{ik{\xi }_1}}\right]\left[{\widehat{E}}_{P_2}(r){\displaystyle \sum _{k^{\prime }=0}^{j^{\prime }}{\mathrm{P}}_{j^{\prime }}^{k^{\prime }}(\cos \theta )e^{i{k}^{\prime }{\xi }_2}}\right]r^2\sin \theta drd\theta d({\omega }_1{t}_c)}}}\\ =\left\{\begin{array}{l}{\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{l}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left|{\widehat{E}}_{P_1}{\widehat{E}}_{P_2}^{*}\right|r^2\sin \theta drd\theta d\phi }}}\mathrm{Opposite}\mathrm{rotational}\mathrm{directions}\hfill \\ 0\mathrm{Same}\mathrm{rotational}\mathrm{direction}\hfill \end{array}\right.\end{array}$$

(49)

where $\left({\xi }_1,{\xi }_2\right)=\left(\pi +{\omega }_1{t}_c,\pi -{\omega }_2{t}_c\right)$ or $\left({\xi }_1,{\xi }_2\right)=\left({\omega }_1{t}_c,2\pi -{\omega }_2{t}_c\right)$.

Thus we can follow the results from the previous chapters when there is the potential energy between two primary particles, and the distance of l remains constant. There are only two conjugate patterns between the two particles because each of the two particles consists of two coupled waves. Therefore, Equation (49) can be translated to

$$\begin{array}{c}{E}_{P_1^P_2}={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}{\displaystyle \underset{l}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left|{\widehat{E}}_{P_1}{\widehat{E}}_{P_2}^{*}\right|r^2\sin \theta drd\theta d\phi }}}={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}\times \\ {\displaystyle \underset{l}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left|\pm [{\widehat{E}}_{a-P_1}+{\widehat{E}}_{b-P_1}]{[{\widehat{E}}_{c-P_2}+{\widehat{E}}_{d-P_2}]}^{*}\pm [{\widehat{E}}_{w-P_1}+{\widehat{E}}_{x-P_1}]{[{\widehat{E}}_{y-P_2}+{\widehat{E}}_{z-P_2}]}^{*}\right|r^2\sin \theta drd\theta d\phi }}}\\ ={\left[E_{a-P_1^c-P_2}+E_{a-P_1^d-P_2}+E_{b-P_1^c-P_2}+E_{b-P_1^d-P_2}\right]}_{r=l+R_c}\\ +{\left[E_{w-P_1^y-P_2}+E_{w-P_1^z-P_2}+E_{x-P_1^y-P_2}+E_{x-P_1^z-P_2}\right]}_{r=l+R_c}\end{array}$$

(50)

where the subscripts a to d, and w to z denote the electric, gravitational, weak interaction, or strong interaction.

According to Newtonian mechanics, the work done is equal to the potential energy when the primary particles I and II move relative to each other. Reversing the Newtonian definition of work, i.e., $F=\nabla W=\nabla E_{P_1^P_2}$, and retaining the habit of using the negative sign as the sign of attraction, the force between two primary particles $F_{P_1^P_2}$ is therefore

$$F_{P_1^P_2}=\pm {\displaystyle \frac{d{E}_{P_1^P_2}}{dl}}$$

(51)

From Equation (49) it can be seen that two self-conjugate primary particles have potential energy or force when they rotate in opposite directions, while they have zero potential energy or rest when they rotate in the same direction. Clearly there are two modes of opposite rotation for two self-conjugate primary particles, which shows that each primary particle has spin values of $\pm {\displaystyle \frac{1}{2}}$.

The result of Equation (50) correlates the conjugate forms of all four fundamental waves, making subsequent calculations based on it lengthy. Therefore, in the following Equations , only the terms associated with the mathematical electric charge $\widehat{q}$ and the mathematical mass $\widehat{m}$ are shown to make the results concise, since the previous assumption was that ${\widehat{q}}_{e\_\nu }={\widehat{w}}_{e\_\nu }$, ${\widehat{q}}_e={\widehat{w}}_e$, ${\widehat{s}}_e={\widehat{s}}_{D\_\nu }={\widehat{m}}_{D\_\nu }$. The assumption R_cS1 << R_cw < R_cS2 is also used to further simplify the results of the later Equations .

3.1. Two Particles of the Same Type

Start by computing the interaction between two electron neutrinos. Combining Equations (18), (20), and (50), we have the potential energy E_{e_ν^e_ν}

$$\begin{array}{c}{E}_{e\_\nu ^e\_\nu }=2{\left[E_{e-e\_\nu ^e-e\_\nu }+2E_{e-e\_\nu ^w-e\_\nu }+E_{w-e\_\nu ^w-e\_\nu }\right]}_{r=l+R_{cw}}\\ =2\left[{\displaystyle \frac{{\widehat{q}}_{e\_\nu }^2}{c_L^2(l+R_{cw})}}+0+{\displaystyle \frac{R_{cw}^{2m-2}{\widehat{w}}_{e\_\nu }^2}{c_L^2(2m-1){(l+R_{cw})}^{2m-1}}}\right]\\ ={\displaystyle \frac{2{\widehat{q}}_{e\_\nu }^2}{c_L^2(l+R_{cw})}}\left[1+{\displaystyle \frac{R_{cw}^{2m-2}}{(2m-1){(l+R_{cw})}^{2m-2}}}\right]\\ \approx {\displaystyle \frac{2{\widehat{q}}_{e\_\nu }^2}{c_L^2(l+R_{cw})}}l\ge R_{cw}\end{array}$$

(52)

Since two Equations (33) have attractive and repulsive states under the Self-Conjugate Mechanism, from Equation (51), the force of two electron neutrinos F_{e_ν^e_ν} is

$$\begin{array}{c}{F}_{e\_\nu ^e\_\nu }=\pm 2\left[{\displaystyle \frac{{\widehat{q}}_{e\_\nu }^2}{c_L^2{(l+R_{cw})}^2}}+{\displaystyle \frac{R_{cw}^{2m-2}{\widehat{w}}_{e\_\nu }^2}{c_L^2{(l+R_{cw})}^{2m}}}\right]\\ =\pm {\displaystyle \frac{2{\widehat{q}}_{e\_\nu }^2}{c_L^2{(l+R_{cw})}^2}}\left[1+{\displaystyle \frac{R_{cw}^{2m-2}}{c_L^2{(l+R_{cw})}^{2m-2}}}\right]\\ \approx \pm {\displaystyle \frac{2{\widehat{q}}_{e\_\nu }^2}{c_L^2{(l+R_{cw})}^2}}l\ge R_{cw}\end{array}$$

(53)

where the sign '–' or '+' depends on the self-conjugate patterns between two electron neutrinos, and the "–" or "+" should be random.

Association Equation (50) with Equations (19), (21), (26), and (38), we can compute the potential energy between two Dark Is E_{D_νI^D_νI}

$$\begin{array}{c}{E}_{D\_\nu I^D\_\nu I}=2{\left[E_{G-D\_\nu I^G-D\_\nu I}+2E_{G-D\_\nu I^S-D\_\nu I}+E_{S-D\_\nu I^S-D\_\nu I}\right]}_{r=l+R_{cS1}}\\ =\left\{\begin{array}{l}\begin{array}{c}{\displaystyle \frac{2c_L^2}{c_G^4}}\left[{\displaystyle \frac{{\widehat{m}}_{D\_\nu }^2}{l+R_{cS1}}}+{\displaystyle \frac{{\widehat{m}}_{D\_\nu }{\widehat{s}}_{D\_\nu }}{R_{cS1}^3}}\left(R_{cS2}^2-{(l+R_{cS1})}^2\right)+{\displaystyle \frac{{\widehat{s}}_{D\_\nu }^2}{R_{cS1}^6}}\left({\displaystyle \frac{R_{cS2}^5-{(l+R_{cS1})}^5}{5}}+{\displaystyle \frac{R_{cS2}^5}{2n-1}}\right)\right]\\ ={\displaystyle \frac{2c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{R}_{cS1}^6}}\left[{\displaystyle \frac{R_{cS1}^6}{l+R_{cS1}}}+R_{cS1}^3\left(R_{cS2}^2-{(l+R_{cS1})}^2\right)+{\displaystyle \frac{R_{cS2}^5-{(l+R_{cS1})}^5}{5}}+{\displaystyle \frac{R_{cS2}^5}{2n-1}}\right]\\ \approx {\displaystyle \frac{2c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{R}_{cS1}^6}}\left[{\displaystyle \frac{R_{cS2}^5-{(l+R_{cS1})}^5}{5}}+{\displaystyle \frac{R_{cS2}^5}{2n-1}}\right]0\le l\le R_{cS2}-R_{cS1}\\ \end{array}\hfill \\ \begin{array}{l}{\displaystyle \frac{2c_L^2}{c_G^4}}\left[{\displaystyle \frac{{\widehat{m}}_{D\_\nu }^2}{l+R_{cS1}}}+0+{\displaystyle \frac{R_{cS2}^{2n+4}}{(2n-1)R_{cS1}^6}}{\displaystyle \frac{{\widehat{s}}_{D\_\nu }^2}{{(l+R_{cS1})}^{2n-1}}}\right]\\ ={\displaystyle \frac{2c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4}}\left[{\displaystyle \frac{1}{l+R_{cS1}}}+{\displaystyle \frac{R_{cS2}^{2n+4}}{(2n-1)R_{cS1}^6{(l+R_{cS1})}^{2n-1}}}\right]\\ ={\displaystyle \frac{2c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4(l+R_{cS1})}}\left[1+{\displaystyle \frac{R_{cS2}^{2n+4}}{(2n-1)R_{cS1}^6{(l+R_{cS1})}^{2n-2}}}\right]l>R_{cS2}-R_{cS1}\end{array}\hfill \end{array}\right.\end{array}$$

(54)

and the force between two Dark Is F_{D_νI^D_νI}

$$F_{D\_\nu I^D\_\nu I}=\left\{\begin{array}{l}\begin{array}{l}-{\displaystyle \frac{2c_L^2}{c_G^4}}\left[{\displaystyle \frac{{\widehat{m}}_{D\_\nu }^2}{{(l+R_{cS1})}^2}}+{\displaystyle \frac{2{\widehat{m}}_{D\_\nu }{\widehat{s}}_{D\_\nu }}{R_{cS1}^3}}(l+R_{cS1})+{\displaystyle \frac{{\widehat{s}}_{D\_\nu }^2}{R_{cS1}^6}}{(l+R_{cS1})}^4\right]\\ =-{\displaystyle \frac{2c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{(l+R_{cS1})}^2}}\left[1+{\displaystyle \frac{2{(l+R_{cS1})}^3}{R_{cS1}^3}}+{\displaystyle \frac{{(l+R_{cS1})}^6}{R_{cS1}^6}}\right]0\le l\le R_{cS2}-R_{cS1}\\ \end{array}\hfill \\ \begin{array}{l}-{\displaystyle \frac{2c_L^2}{c_G^4}}\left[{\displaystyle \frac{{\widehat{m}}_{D\_\nu }^2}{{(l+R_{cS1})}^2}}+0+{\displaystyle \frac{R_{cS2}^{2n+4}}{R_{cS1}^6}}{\displaystyle \frac{{\widehat{s}}_{D\_\nu }^2}{{(l+R_{cS1})}^{2n}}}\right]\\ =-{\displaystyle \frac{2c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{(l+R_{cS1})}^2}}\left[1+{\displaystyle \frac{R_{cS2}^{2n+4}}{R_{cS1}^6{(l+R_{cS1})}^{2n-2}}}\right]l>R_{cS2}-R_{cS1}\end{array}\hfill \end{array}\right.$$

(55)

Duplicating the last process yields the potential energy between two Dark IIs E_{D_νII^D_νII}

$$\begin{array}{c}{E}_{D\_\nu II^D\_\nu II}=2{\left[E_{G-D\_\nu II^G-D\_\nu II}+2E_{G-D\_\nu II^S-D\_\nu II}+E_{S-D\_\nu II^S-D\_\nu II}\right]}_{r=l+R_{cS2}}\\ ={\displaystyle \frac{2c_L^2}{c_G^4}}\left[{\displaystyle \frac{{\widehat{m}}_{D\_\nu }^2}{l+R_{cS2}}}+0+{\displaystyle \frac{R_{cS2}^{2n+4}}{(2n-1)R_{cS1}^6}}{\displaystyle \frac{{\widehat{s}}_{D\_\nu }^2}{{(l+R_{cS2})}^{2n-1}}}\right]\\ ={\displaystyle \frac{2c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4(l+R_{cS2})}}\left[1+{\displaystyle \frac{R_{cS2}^{2n+4}}{(2n-1)R_{cS1}^6{(l+R_{cS2})}^{2n-2}}}\right]\end{array}$$

(56)

and the force between two Dark IIs F_{D_νII^D_νII}

$$\begin{array}{c}{F}_{D\_\nu II^D\_\nu II}=-{\displaystyle \frac{2c_L^2}{c_G^4}}\left[{\displaystyle \frac{{\widehat{m}}_{D\_\nu }^2}{{(l+R_{cS2})}^2}}+{\displaystyle \frac{R_{cS2}^{2n+4}}{R_{cS1}^6}}{\displaystyle \frac{{\widehat{s}}_{D\_\nu }^2}{{(l+R_{cS2})}^{2n}}}\right]\\ =-{\displaystyle \frac{2c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{(l+R_{cS2})}^2}}\left[1+{\displaystyle \frac{R_{cS2}^{2n+4}}{R_{cS1}^6{(l+R_{cS2})}^{2n-2}}}\right]\end{array}$$

(57)

The potential energy between two electrons (or two positrons) and between an electron and a positron E_e^e has two forms because an electron (or a positron) consists of an E-W couple and a G-S couple. As with electron neutrinos, the two forms should be random and based on the self-conjugate patterns between two electrons. Combining Equations (18) to (27) and (50), we can calculate the two forms of the potential energy. One is E_e^eI when the two E-W couples are conjugate, and the two G-S couples are conjugate (hereafter referred to as the E-P I (Electric Pattern I)).

$$\begin{array}{c}{E}_{e^eI}={\left[E_{e-e^e-e}+2E_{e-e^w-e}+E_{w-e^w-e}+E_{G-e^G-e}+2E_{G-e^S-e}+E_{S-e^S-e}\right]}_{r=l+R_{cw}}\\ =\left\{\begin{array}{l}\begin{array}{l}{\displaystyle \frac{{\widehat{q}}_e^2}{c_L^2(l+R_{cw})}}+0+{\displaystyle \frac{R_{cw}^{2m-2}{\widehat{w}}_e^2}{c_L^2(2m-1){(l+R_{cw})}^{2m-1}}}+{\displaystyle \frac{c_L^2{\widehat{m}}_e^2}{c_G^4(l+R_{cw})}}+{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{s}}_e}{c_G^4{R}_{cS1}^3}}\left[R_{cS2}^2-{(l+R_{cw})}^2\right]\\ +{\displaystyle \frac{c_L^2{\widehat{s}}_e^2}{c_G^4}}\left[{\displaystyle \frac{1}{5R_{cS1}^6}}\left(R_{cS2}^5-{(l+R_{cw})}^5\right)+{\displaystyle \frac{R_{cS2}^5}{(2n-1)R_{cS1}^6}}\right]0\le l\le R_{cS2}-R_{cw}\\ \end{array}\hfill \\ \begin{array}{l}{\displaystyle \frac{{\widehat{q}}_e^2}{c_L^2(l+R_{cw})}}+{\displaystyle \frac{R_{cw}^{2m-2}{\widehat{w}}_e^2}{c_L^2(2m-1){(l+R_{cw})}^{2m-1}}}+{\displaystyle \frac{c_L^2{\widehat{m}}_e^2}{c_G^4(l+R_{cw})}}\\ +{\displaystyle \frac{c_L^2{R}_{cS2}^{2n+4}}{c_G^4(2n-1)R_{cS1}^6}}{\displaystyle \frac{{\widehat{s}}_e^2}{{(l+R_{cw})}^{2n-1}}}l>R_{cS2}-R_{cw}\end{array}\hfill \end{array}\right.\\ =\left\{\begin{array}{l}\begin{array}{l}{\displaystyle \frac{{\widehat{q}}_e^2}{c_L^2(l+R_{cw})}}\left[1+{\displaystyle \frac{R_{cw}^{2m-2}}{c_L^2(2m-1){(l+R_{cw})}^{2m-2}}}\right]+{\displaystyle \frac{c_L^2{\widehat{m}}_e^2}{c_G^4(l+R_{cw})}}+{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4{R}_{cS1}^3}}\left[R_{cS2}^2-{(l+R_{cw})}^2\right]\\ +{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4}}\left[{\displaystyle \frac{1}{5R_{cS1}^6}}\left(R_{cS2}^5-{(l+R_{cw})}^5\right)+{\displaystyle \frac{R_{cS2}^5}{(2n-1)R_{cS1}^6}}\right]0\le l\le R_{cS2}-R_{cw}\\ \end{array}\hfill \\ \begin{array}{l}{\displaystyle \frac{{\widehat{q}}_e^2}{c_L^2(l+R_{cw})}}\left[1+{\displaystyle \frac{R_{cw}^{2m-2}}{c_L^2(2m-1){(l+R_{cw})}^{2m-2}}}\right]+{\displaystyle \frac{c_L^2{\widehat{m}}_e^2}{c_G^4(l+R_{cw})}}\\ +{\displaystyle \frac{c_L^2{R}_{cS2}^{2n+4}}{c_G^4(2n-1)R_{cS1}^6}}{\displaystyle \frac{{\widehat{m}}_{D\_\nu }^2}{{(l+R_{cw})}^{2n-1}}}l>R_{cS2}-R_{cw}\end{array}\hfill \end{array}\right.\\ \approx {\displaystyle \frac{{\widehat{q}}_e^2}{c_L^2(l+R_{cw})}}+{\displaystyle \frac{c_L^2{\widehat{m}}_e^2}{c_G^4(l+R_{cw})}}l\mathrm{is}\mathrm{large}\end{array}$$

(58)

Another is E_e^eII when the two E-W couples are conjugate to the two G-S couples (hereafter referred to as the E-P II).

$$\begin{array}{c}{E}_{e^eII}=2{\left[E_{e-e^G-e}+E_{e-e^S-e}+E_{w-e^G-e}+E_{w-e^S-e}\right]}_{r=l+R_{cw}}\\ =\left\{\begin{array}{l}{\displaystyle \frac{2{\widehat{q}}_e{\widehat{m}}_e}{c_G^2(l+R_{cw})}}+{\displaystyle \frac{{\widehat{q}}_e{\widehat{s}}_e}{c_G^2{R}_{cS1}^3}}\left[R_{cS2}^2-{(l+R_{cw})}^2\right]+0+0\hspace{1em}0\le l\le R_{cS2}-R_{cw}\hfill \\ {\displaystyle \frac{2{\widehat{q}}_e{\widehat{m}}_e}{c_G^2(l+R_{cw})}}+0+0+0l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ =\left\{\begin{array}{l}{\displaystyle \frac{2{\widehat{q}}_e{\widehat{m}}_e}{c_G^2(l+R_{cw})}}+{\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2{R}_{cS1}^3}}\left[R_{cS2}^2-{(l+R_{cw})}^2\right]0\le l\le R_{cS2}-R_{cw}\hfill \\ {\displaystyle \frac{2{\widehat{q}}_e{\widehat{m}}_e}{c_G^2(l+R_{cw})}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\end{array}$$

(59)

Derivation of the last two Equations can yield the forces between two electrons (or two positrons) in both forms that are

$$\begin{array}{c}{F}_{e^eI}=\left\{\begin{array}{l}\begin{array}{l}\pm {\displaystyle \frac{{\widehat{q}}_e^2}{c_L^2{(l+R_{cw})}^2}}\pm {\displaystyle \frac{R_{cw}^{2m-2}{\widehat{w}}_e^2}{c_L^2{(l+R_{cw})}^{2m}}}-{\displaystyle \frac{c_L^2{\widehat{m}}_e^2}{c_G^4{(l+R_{cw})}^2}}-{\displaystyle \frac{2c_L^2{\widehat{m}}_e{\widehat{s}}_e}{c_G^4{R}_{cS1}^3}}(l+R_{cw})\\ -{\displaystyle \frac{c_L^2{\widehat{s}}_e^2}{c_G^4{R}_{cS1}^6}}(l+R_{cw})4 0\le l\le R_{cS2}-R_{cw}\end{array}\hfill \\ \pm {\displaystyle \frac{{\widehat{q}}_e^2}{c_L^2{(l+R_{cw})}^2}}\pm {\displaystyle \frac{R_{cw}^{2m-2}{\widehat{w}}_e^2}{c_L^2{(l+R_{cw})}^{2m}}}-{\displaystyle \frac{c_L^2{\widehat{m}}_e^2}{c_G^4{(l+R_{cw})}^2}}-{\displaystyle \frac{c_L^2{R}_{cS2}^{2n+4}}{c_G^4{R}_{cS1}^6}}{\displaystyle \frac{{\widehat{s}}_e^2}{{(l+R_{cw})}^{2n}}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ =\left\{\begin{array}{l}\begin{array}{l}\pm {\displaystyle \frac{{\widehat{q}}_e^2}{c_L^2{(l+R_{cw})}^2}}\left[1+{\displaystyle \frac{R_{cw}^{2m-2}}{{(l+R_{cw})}^{2m-2}}}\right]-{\displaystyle \frac{c_L^2{\widehat{m}}_e^2}{c_G^4{(l+R_{cw})}^2}}-{\displaystyle \frac{2c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4{R}_{cS1}^3}}(l+R_{cw})\\ -{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{R}_{cS1}^6}}{(l+R_{cw})}^{4}0\le l\le R_{cS2}-R_{cw}\end{array}\hfill \\ \pm {\displaystyle \frac{{\widehat{q}}_e^2}{c_L^2{(l+R_{cw})}^2}}\left[1+{\displaystyle \frac{R_{cw}^{2m-2}}{{(l+R_{cw})}^{2m-2}}}\right]-{\displaystyle \frac{c_L^2{\widehat{m}}_e^2}{c_G^4{(l+R_{cw})}^2}}-{\displaystyle \frac{c_L^2{R}_{cS2}^{2n+4}}{c_G^4{R}_{cS1}^6}}{\displaystyle \frac{{\widehat{m}}_{D\_\nu }^2}{{(l+R_{cw})}^{2n}}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ \approx \pm {\displaystyle \frac{{\widehat{q}}_e^2}{c_L^2{(l+R_{cw})}^2}}-{\displaystyle \frac{c_L^2{\widehat{m}}_e^2}{c_G^4{(l+R_{cw})}^2}}l\mathrm{is}\mathrm{large}\end{array}$$

(60)

and

$$\begin{array}{c}{F}_{e^eII}=\left\{\begin{array}{l}-{\displaystyle \frac{2{\widehat{q}}_e{\widehat{m}}_e}{c_G^2{(l+R_{cw})}^2}}-{\displaystyle \frac{2{\widehat{q}}_e{\widehat{s}}_e}{c_G^2{R}_{cS1}^3}}(l+R_{cw})0\le l\le R_{cS2}-R_{cw}\hfill \\ -{\displaystyle \frac{2{\widehat{q}}_e{\widehat{m}}_e}{c_G^2{(l+R_{cw})}^2}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ =\left\{\begin{array}{l}-{\displaystyle \frac{2{\widehat{q}}_e{\widehat{m}}_e}{c_G^2{(l+R_{cw})}^2}}-{\displaystyle \frac{2{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2{R}_{cS1}^3}}(l+R_{cw})0\le l\le R_{cS2}-R_{cw}\hfill \\ -{\displaystyle \frac{2{\widehat{q}}_e{\widehat{m}}_e}{c_G^2{(l+R_{cw})}^2}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\end{array}$$

(61)

A comparison of Equations (60) and (61) shows that the force between two electrons are very different in the two self-conjugate patterns. Equation (60) has the square terms of the mathematical electric charge ${\widehat{q}}_e$ and the mathematical weak charge ${\widehat{w}}_e$, which makes the force properties obey the charges of the two electrons. However, there are only the product terms of the mathematical electric charge ${\widehat{q}}_e$ and the mathematical mass ${\widehat{m}}_e$, and the mathematical electric charge ${\widehat{q}}_e$ and the mathematical strong charge ${\widehat{s}}_e$ in Equation (61), which means that the force is attractive. This indicates that even two electrons (or two positrons) can come together with small potential energies and forces if they interact with each other in the form E-P II.

3.2. Two Particles of the Different Type

Similar to the last section, this one starts by computing the interaction of an electron neutrino with another primary particle. Combining Equations (22) to (27), (50), and (51), the potential energy between an electron neutrino and a Dark I E_{e_ν^D_νI} is

$$\begin{array}{c}{E}_{e\_\nu ^D\_\nu I}=2{\left[E_{e-e\_\nu ^G-D\_\nu I}+E_{e-e\_\nu ^S-D\_\nu I}+E_{w-e\_\nu ^G-D\_\nu I}+E_{w-e\_\nu ^S-D\_\nu I}\right]}_{r=l+R_{cw}}\\ =\left\{\begin{array}{l}{\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cw})}}+{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{s}}_{D\_\nu }}{c_G^2{R}_{cS1}^3}}\left(R_{cS2}^2-{(l+R_{cw})}^2\right)+0+0 0\le l\le R_{cS2}-R_{cw}\hfill \\ {\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cw})}}+0+0+0l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ =\left\{\begin{array}{l}{\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cw})}}\left[1+{\displaystyle \frac{R_{cS2}^2(l+R_{cw})-{(l+R_{cw})}^3}{2R_{cS1}^3}}\right]0\le l\le R_{cS2}-R_{cw}\hfill \\ {\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cw})}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\end{array}$$

(62)

and the force between an electron neutrino and a Dark I F_{e_ν^D_νI} is

$$\begin{array}{c}{F}_{e\_\nu ^D\_\nu I}=\left\{\begin{array}{l}-{\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{m}}_{D\_\nu }}{c_G^2{(l+R_{cw})}^2}}-{\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{s}}_{D\_\nu }}{c_G^2{R}_{cS1}^3}}(l+R_{cw})0\le l\le R_{cS2}-R_{cw}\hfill \\ -{\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{m}}_{D\_\nu }}{c_G^2{(l+R_{cw})}^2}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ =\left\{\begin{array}{l}-{\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{m}}_{D\_\nu }}{c_G^2{(l+R_{cw})}^2}}\left[1-{\displaystyle \frac{{(l+R_{cw})}^3}{R_{cS1}^3}}\right]0\le l\le R_{cS2}-R_{cw}\hfill \\ -{\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{m}}_{D\_\nu }}{c_G^2{(l+R_{cw})}^2}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\end{array}$$

(63)

Repeating the previous processes gives the potential energy between an electron neutrino and a Dark II E_{e_ν^D_νII}

$$\begin{array}{c}{E}_{e\_\nu ^D\_\nu II}=2{\left[E_{e-e\_\nu ^G-D\_\nu II}+E_{e-e\_\nu ^S-D\_\nu II}+E_{w-e\_\nu ^G-D\_\nu II}+E_{w-e\_\nu ^S-D\_\nu II}\right]}_{r=l+R_{cS2}}\\ ={\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cS2})}}+0+0+0\\ ={\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cS2})}}\end{array}$$

(64)

the force between an electron neutrino and a Dark I F_{e_ν^D_νII}

$$F_{e\_\nu ^D\_\nu II}=-{\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{m}}_{D\_\nu }}{c_G^2{(l+R_{cS2})}^2}}$$

(65)

the potential energy between an electron neutrino and an electron E_{e_ν^e}

$$\begin{array}{c}{E}_{e\_\nu ^e}={\left[E_{e-e\_\nu ^e-e}+E_{e-e\_\nu ^w-e}+E_{w-e\_\nu ^e-e}+E_{w-e\_\nu ^w-e}\right]}_{r=l+R_{cw}}\\ +{\left[E_{e-e\_\nu ^G-e}+E_{e-e\_\nu ^S-e}+E_{w-e\_\nu ^G-e}+E_{w-e\_\nu ^S-e}\right]}_{r=l+R_{cw}}\\ =\left\{\begin{array}{l}\begin{array}{l}{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{q}}_e}{c_L^2(l+R_{cw})}}+0+0+{\displaystyle \frac{R_{cw}^{2m-2}{\widehat{w}}_{e\_\nu }{\widehat{w}}_e}{c_L^2(2m-1){(l+R_{cw})}^{2m-1}}}\\ +{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{m}}_e}{c_G^2(l+R_{cw})}}+{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{s}}_e}{2c_G^2{R}_{cS1}^3}}\left(R_{cS2}^2-{(l+R_{cw})}^2\right)+0+0\hspace{1em}0\le l\le R_{cS2}-R_{cw}\\ \end{array}\hfill \\ {\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{q}}_e}{c_L^2(l+R_{cw})}}+0+0+{\displaystyle \frac{R_{cw}^{2m-2}{\widehat{w}}_{e\_\nu }{\widehat{w}}_e}{c_L^2(2m-1){(l+R_{cw})}^{2m-1}}}+{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{m}}_e}{c_G^2(l+R_{cw})}}+0+0+0l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ =\left\{\begin{array}{l}\begin{array}{l}{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{q}}_e}{c_L^2(l+R_{cw})}}\left[1+{\displaystyle \frac{R_{cw}^{2m-2}}{(2m-1){(l+R_{cw})}^{2m-2}}}\right]+{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{m}}_e}{c_G^2(l+R_{cw})}}\\ +{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{m}}_{D\_\nu }}{2c_G^2{R}_{cS1}^3}}\left(R_{cS2}^2-{(l+R_{cw})}^2\right)0\le l\le R_{cS2}-R_{cw}\\ \end{array}\hfill \\ {\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{q}}_e}{c_L^2(l+R_{cw})}}\left[1+{\displaystyle \frac{R_{cw}^{2m-2}}{(2m-1){(l+R_{cw})}^{2m-2}}}\right]l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ \approx {\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{q}}_e}{c_L^2(l+R_{cw})}}l\ge R_{cw}\end{array}$$

(66)

and the force between an electron neutrino and an electron F_{e_ν^e}

$$\begin{array}{c}{F}_{e\_\nu ^e}=\left\{\begin{array}{l}\pm {\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{q}}_e}{c_L^2{(l+R_{cw})}^2}}\pm {\displaystyle \frac{R_{cw}^{2m-2}{\widehat{w}}_{e\_\nu }{\widehat{w}}_e}{c_L^2{(l+R_{cw})}^{2m}}}-{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{m}}_e}{c_G^2{(l+R_{cw})}^2}}-{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{s}}_e}{c_G^2{R}_{cS1}^3}}(l+R_{cw})0\le l\le R_{cS2}-R_{cw}\hfill \\ \pm {\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{q}}_e}{c_L^2{(l+R_{cw})}^2}}\pm {\displaystyle \frac{R_{cw}^{2m-2}{\widehat{w}}_{e\_\nu }{\widehat{w}}_e}{c_L^2{(l+R_{cw})}^{2m}}}-{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{m}}_e}{c_G^2{(l+R_{cw})}^2}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ =\left\{\begin{array}{l}\pm {\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{q}}_e}{c_L^2{(l+R_{cw})}^2}}\left[1+{\displaystyle \frac{R_{cw}^{2m-2}}{{(l+R_{cw})}^{2m-2}}}\right]-{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{m}}_e}{c_G^2{(l+R_{cw})}^2}}-{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{m}}_{D\_\nu }}{c_G^2{R}_{cS1}^3}}(l+R_{cw})0\le l\le R_{cS2}-R_{cw}\hfill \\ \pm {\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{q}}_e}{c_L^2{(l+R_{cw})}^2}}\left[1+{\displaystyle \frac{R_{cw}^{2m-2}}{{(l+R_{cw})}^{2m-2}}}\right]l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ \approx \pm {\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{q}}_e}{c_L^2{(l+R_{cw})}^2}}l\ge R_{cw}\end{array}$$

(67)

where the sign '–' or '+' randomizes in Equations (67).

It is next computed that the potential energies and forces between the two types of dark neutrinos and between each of them and an electron. Based on Equations (19), (21), (26), (50), and (51), the potential energy between a Dark I and a Dark II E_{D_νI^D_νII} is

$$\begin{array}{c}{E}_{D\_\nu I^D\_\nu II}=2{\left[E_{G-D\_\nu I^G-D\_\nu II}+E_{G-D\_\nu I^S-D\_\nu II}+E_{S-D\_\nu I^G-D\_\nu II}+E_{S-D\_\nu I^S-D\_\nu II}\right]}_{r=l+R_{cS2}}\\ =2\left[{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4(l+R_{cS2})}}+0+0+{\displaystyle \frac{c_L^2{R}_{cS2}^{2n+4}}{c_G^4(2n-1)R_{cS1}^6}}{\displaystyle \frac{{\widehat{s}}_{D\_\nu }^2}{{(l+R_{cS2})}^{2n-1}}}\right]\\ ={\displaystyle \frac{2c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4(l+R_{cS2})}}\left[1+{\displaystyle \frac{R_{cS2}^{2n+4}}{(2n-1)R_{cS1}^6{(l+R_{cS2})}^{2n-2}}}\right]\end{array}$$

(68)

and the force between a Dark I and a Dark II F_{D_νI^D_νII} is

$$\begin{array}{c}{F}_{D\_\nu I^D\_\nu II}=-{\displaystyle \frac{2c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{(l+R_{cS2})}^2}}-{\displaystyle \frac{2c_L^2{R}_{cS2}^{2n+4}}{c_G^4{R}_{cS1}^6}}{\displaystyle \frac{{\widehat{s}}_{D\_\nu }^2}{{(l+R_{cS2})}^{2n}}}\\ =-{\displaystyle \frac{2c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{(l+R_{cS2})}^2}}\left[1+{\displaystyle \frac{R_{cS2}^{2n+4}}{R_{cS1}^6{(l+R_{cS2})}^{2n-2}}}\right]\end{array}$$

(69)

Based on Equations (22) to (27), (42), (50), and (51), the potential energy between a Dark I and an electron E_{D_νI^e} is

$$\begin{array}{c}{E}_{D\_\nu I^e}={\left[E_{G-D\_\nu I^e-e}+E_{G-D\_\nu I^w-e}+E_{S-D\_\nu I^e-e}+E_{S-D\_\nu I^w-e}\right]}_{r=l+R_{cw}}\\ +{\left[E_{G-D\_\nu I^G-e}+E_{G-D\_\nu I^S-e}+E_{S-D\_\nu I^G-e}+E_{S-D\_\nu I^S-e}\right]}_{r=l+R_{cw}}\\ =\left\{\begin{array}{l}\begin{array}{l}{\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cw})}}+0+{\displaystyle \frac{{\widehat{q}}_e{\widehat{s}}_{D\_\nu }}{2c_G^2{R}_{cS1}^3}}\left(R_{cS2}^2-{(l+R_{cw})}^2\right)+0+{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4(l+R_{cw})}}\\ +{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }{\widehat{s}}_e}{2c_G^4{R}_{cS1}^3}}\left(R_{cS2}^2-{(l+R_{cw})}^2\right)+{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{s}}_{D\_\nu }}{2c_G^4{R}_{cS1}^3}}\left(R_{cS2}^2-{(l+R_{cw})}^2\right)\\ +{\displaystyle \frac{c_L^2{\widehat{s}}_{D\_\nu }{\widehat{s}}_e}{c_G^4}}\left[{\displaystyle \frac{1}{5R_{cS1}^6}}\left(R_{cS2}^5-{(l+R_{cw})}^5\right)+{\displaystyle \frac{R_{cS2}^5}{(2n-1)R_{cS1}^6}}\right]0\ge l\ge R_{cS2}-R_{cw}\\ \end{array}\hfill \\ {\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cw})}}+0+0+0+{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4(l+R_{cw})}}+0+0+{\displaystyle \frac{c_L^2{\widehat{s}}_{D\_\nu }{\widehat{s}}_e{R}_{cS2}^{2n+4}}{c_G^4(2n-1)R_{cS1}^6{(l+R_{cw})}^{2n-1}}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ =\left\{\begin{array}{l}\begin{array}{l}{\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2}}\left[{\displaystyle \frac{1}{l+R_{cw}}}+{\displaystyle \frac{R_{cS2}^2-{(l+R_{cw})}^2}{2R_{cS1}^3}}\right]+{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4}}\left[{\displaystyle \frac{1}{l+R_{cw}}}+{\displaystyle \frac{R_{cS2}^2-{(l+R_{cw})}^2}{2R_{cS1}^3}}\right]\\ +{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{R}_{cS1}^3}}\left[{\displaystyle \frac{R_{cS2}^2-{(l+R_{cw})}^2}{2}}+{\displaystyle \frac{R_{cS2}^5-{(l+R_{cw})}^5}{5R_{cS1}^3}}+{\displaystyle \frac{R_{cS2}^5}{(2n-1)R_{cS1}^3}}\right]0\ge l\ge R_{cS2}-R_{cw}\\ \end{array}\hfill \\ {\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cw})}}+{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4(l+R_{cw})}}+{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2{R}_{cS2}^{2n+4}}{c_G^4(2n-1)R_{cS1}^6{(l+R_{cw})}^{2n-1}}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ \approx \left\{\begin{array}{l}{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{R}_{cS1}^6}}\left[{\displaystyle \frac{R_{cS2}^5-{(l+R_{cw})}^5}{5}}+{\displaystyle \frac{R_{cS2}^5}{2n-1}}\right]0\ge l\ge R_{cS2}-R_{cw}\hfill \\ {\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cw})}}+{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2{R}_{cS2}^{2n+4}}{c_G^4(2n-1)R_{cS1}^6{(l+R_{cw})}^{2n-1}}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\end{array}$$

(70)

the force between a Dark I and an electron F_{D_νI^e} is

$$\begin{array}{c}{F}_{D\_\nu I^e}=\left\{\begin{array}{l}\begin{array}{l}-{\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2{(l+R_{cw})}^2}}-{\displaystyle \frac{{\widehat{q}}_e{\widehat{s}}_{D\_\nu }}{c_G^2{R}_{cS1}^3}}(l+R_{cw})-{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4{(l+R_{cw})}^2}}\\ -{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }{\widehat{s}}_e}{c_G^4{R}_{cS1}^3}}(l+R_{cw})-{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{s}}_{D\_\nu }}{c_G^4{R}_{cS1}^3}}(l+R_{cw})-{\displaystyle \frac{c_L^2{\widehat{s}}_{D\_\nu }{\widehat{s}}_e}{c_G^4{R}_{cS1}^6}}{(l+R_{cw})}^{4}0\ge l\ge R_{cS2}-R_{cw}\\ \end{array}\hfill \\ -{\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2{(l+R_{cw})}^2}}-{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4{(l+R_{cw})}^2}}-{\displaystyle \frac{c_L^2{\widehat{s}}_{D\_\nu }{\widehat{s}}_e{R}_{cS2}^{2n+4}}{c_G^4{R}_{cS1}^6{(l+R_{cw})}^{2n}}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ =\left\{\begin{array}{l}\begin{array}{l}-{\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2{(l+R_{cw})}^2}}\left[1+{\displaystyle \frac{{(l+R_{cw})}^3}{R_{cS1}^3}}\right]-{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4{(l+R_{cw})}^2}}\left[1+{\displaystyle \frac{{(l+R_{cw})}^3}{R_{cS1}^3}}\right]\\ -{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{R}_{cS1}^3}}(l+R_{cw})\left[1+{\displaystyle \frac{{(l+R_{cw})}^3}{R_{cS1}^3}}\right]0\ge l\ge R_{cS2}-R_{cw}\\ \end{array}\hfill \\ -{\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2{(l+R_{cw})}^2}}-{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4{(l+R_{cw})}^2}}-{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2{R}_{cS2}^{2n+4}}{c_G^4{R}_{cS1}^6{(l+R_{cw})}^{2n}}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\\ \approx \left\{\begin{array}{l}-{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{R}_{cS1}^6}}{(l+R_{cw})}^{4}0\ge l\ge R_{cS2}-R_{cw}\hfill \\ -{\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2{(l+R_{cw})}^2}}-{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4{(l+R_{cw})}^2}}-{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2{R}_{cS2}^{2n+4}}{c_G^4{R}_{cS1}^6{(l+R_{cw})}^{2n}}}l>R_{cS2}-R_{cw}\hfill \end{array}\right.\end{array}$$

(71)

the potential energy between a Dark II and an electron E_{D_νII^e} is

$$\begin{array}{c}{E}_{D\_\nu II^e}={\left[E_{G-D\_\nu II^e-e}+E_{G-D\_\nu II^w-e}+E_{S-D\_\nu II^e-e}+E_{S-D\_\nu II^w-e}\right]}_{r=l+R_{cS2}}\\ +{\left[E_{G-D\_\nu II^G-e}+E_{G-D\_\nu II^S-e}+E_{S-D\_\nu II^G-e}+E_{S-D\_\nu II^S-e}\right]}_{r=l+R_{cS2}}\\ ={\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cS2})}}+0+0+0+{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4(l+R_{cS2})}}+0+0+{\displaystyle \frac{c_L^2{\widehat{s}}_{D\_\nu }{\widehat{s}}_e{R}_{cS2}^{2n+4}}{c_G^4(2n-1)R_{cS1}^6{(l+R_{cS2})}^{2n-1}}}\\ ={\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cS2})}}+{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4(l+R_{cS2})}}+{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2{R}_{cS2}^{2n+4}}{c_G^4(2n-1)R_{cS1}^6{(l+R_{cS2})}^{2n-1}}}\end{array}$$

(72)

and the force between a Dark II and an electron F_{D_νII^e} is

$$F_{D\_\nu II^e}=-{\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2{(l+R_{cS2})}^2}}-{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4{(l+R_{cS2})}^2}}-{\displaystyle \frac{c_L^2{\widehat{s}}_{D\_\nu }{\widehat{s}}_e{R}_{cS2}^{2n+4}}{c_G^4{R}_{cS1}^6{(l+R_{cS2})}^{2n}}}$$

(73)

4. The Structure Values of Primary Particles

In this chapter, we attempt to incorporate the existing physical data into the computational results of the previous chapters in order to obtain the structural values of primary particles. Nowadays, we are fully aware of the characteristics of electricity, such as the charge, the potential energy, and the field, at both the macroscopic and microscopic levels. We can be confident that the available measurements reflect the characteristics of the electric charge and no other factors. Thus, comparing Equation (44) with Gauss's law of electrostatics and combining Equation (28), it is easy to see that the relationship between the mathematical electric charge of an electron ${\widehat{q}}_e$ and an electron charge q_e in the present physical data.

$${\widehat{q}}_e=c_L\sqrt{k}{q}_e$$

(74)

It is somewhat difficult to obtain the charges of gravity and the strong force because gravity is the weakest of the four fundamental interactions, and the strong force is a short-range interaction. In the models of this paper, although the gravitational and inertial masses are different, it can be hypothesized that Newton's law of gravity is largely accurate because only electron neutrinos in the model do not have gravitational waves and because electron neutrinos have much less binding energy or rest mass. We focus on protons ---- stable, heavy subatomic particles ---- to explore the charges of dark neutrinos. A proton is supposed to consist of one Dark II and one positron with the opposite spin directions, because a combination of three or more primary particles must have one or more pairs of primary particles rotating in the same direction or in the rest state. In other words, the stability of the combination decreases as the number of primary particles increases. From Equations (72) and (73), a proton could have the structure shown in Figure 3, and the two primary particles are close together, i.e., l = 0.

It is further supposed that the non-potential energies of the two primary particles are negligibly small compared with their binding energies and potential energy. Hence, the energy of a proton E_p equals

$$\begin{array}{c}{E}_p=E_{D\_\nu II}+E_e+{\left.E_{D\_\nu II^e}\right|}_{l=0}\\ \approx {\displaystyle \frac{c_L^2{R}_{cS2}^5{\widehat{m}}_{D\_\nu }^2}{(2n-1)c_G^4{R}_{cS1}^6}}+E_e+{\left[{\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cS2})}}+{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4(l+R_{cS2})}}+{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2{R}_{cS2}^{2n+4}}{c_G^4(2n-1)R_{cS1}^6{(l+R_{cS2})}^{2n-1}}}\right]}_{l=0}\\ \approx {\displaystyle \frac{c_L^2{R}_{cS2}^5{\widehat{m}}_{D\_\nu }^2}{(2n-1)c_G^4{R}_{cS1}^6}}+E_e+{\displaystyle \frac{c_L^2{R}_{cS2}^5{\widehat{m}}_{D\_\nu }^2}{(2n-1)c_G^4{R}_{cS1}^6}}\\ \approx {\displaystyle \frac{2c_L^2{R}_{cS2}^5{\widehat{m}}_{D\_\nu }^2}{(2n-1)c_G^4{R}_{cS1}^6}}\approx 2E_{D\_\nu II}\end{array}$$

(75)

when Equations (37) and (72) are subsequently substituted into the above Equation, and the energy of an electron E_e is neglected in the above Equation because E_p is 1836 times E_e (E_p = 938MeV, E_e = 0.511MeV)[7]. Furthermore,

$${\widehat{m}}_{D\_\nu }={\displaystyle \frac{c_G^2\sqrt{G}}{2c_L}}{m}_p$$

(76)

can be obtained from Equations (29) and (39), where m_p is the mass of a proton.

Next, by combining Equations (42), (43), and (74) to (76), we have

$$\left\{\begin{array}{l}{\displaystyle \frac{R_{cw}^3}{R_{cS1}^3}}={\displaystyle \frac{4\sqrt{k}{q}_e}{\sqrt{G}{m}_p}}\hfill \\ {\displaystyle \frac{\sqrt{k}{q}_e\sqrt{G}{m}_p}{4}}{\displaystyle \frac{R_{cS2}^2-R_{cw}^2}{R_{cS1}^3}}=E_e\hfill \\ {\displaystyle \frac{G{m}_p^2}{2(2n-1)}}{\displaystyle \frac{R_{cS2}^5}{R_{cS1}^6}}=E_p\hfill \end{array}\right.$$

(77)

that we can solve and yield

$$\begin{array}{l}{\left[{\displaystyle \frac{2(2n-1)E_p}{G{m}_p^2}}\right]}^{{\displaystyle \frac{2}{5}}}{R}_{cS1}^{{\displaystyle \frac{2}{5}}}-{\displaystyle \frac{4E_e}{\sqrt{k}{q}_e\sqrt{G}{m}_p}}{R}_{cS1}={\left[{\displaystyle \frac{4\sqrt{k}{q}_e}{\sqrt{G}{m}_p}}\right]}^{{\displaystyle \frac{2}{3}}}\\ {\left[{\displaystyle \frac{2(2n-1)\times 938.272\times {10}^6\times 1.60218\times {10}^{-19}}{6.67430\times {10}^{-11}\times {(1.67262\times {10}^{-27})}^2}}\right]}^{{\displaystyle \frac{2}{5}}}{R}_{cS1}^{{\displaystyle \frac{2}{5}}}\\ -{\displaystyle \frac{4\times 510.999\times {10}^3\times 1.60218\times {10}^{-19}}{\sqrt{8.98755\times {10}^9}\times 1.60218\times {10}^{-19}\times \sqrt{6.67430\times {10}^{-11}}\times 1.67262\times {10}^{-27}}}{R}_{cS1}\\ ={\left[{\displaystyle \frac{4\sqrt{8.98755\times {10}^9}\times 1.60218\times {10}^{-19}}{\sqrt{6.67430\times {10}^{-11}}\times 1.67262\times {10}^{-27}}}\right]}^{{\displaystyle \frac{2}{3}}}\\ {\left[(2n-1)R_{cS1}\right]}^{{\displaystyle \frac{2}{5}}}-3.2758\times {10}^{11}{R}_{cS1}=5.6137\times {10}^{-10}\end{array}$$

(78)

Evidently, R_cS1 is very small from the above Equation. Hence

$$\begin{array}{l}{\left[(2n-1)R_{cS1}\right]}^{{\displaystyle \frac{2}{5}}}\approx 5.6137\times {10}^{-10}\\ R_{cS1}={\displaystyle \frac{7.47\times {10}^{-24}}{2n-1}}\left(\mathrm{m}\right)={\displaystyle \frac{7.47\times {10}^{-9}}{2n-1}}\left(\mathrm{fm}\right)\end{array}$$

(79)

Next

$$\begin{array}{c}{R}_{cw}={\left({\displaystyle \frac{4\sqrt{k}{q}_e}{\sqrt{G}{m}_p}}\right)}^{{\displaystyle \frac{1}{3}}}{R}_{cS1}={\left({\displaystyle \frac{4\sqrt{8.98755\times {10}^9}\times 1.60218\times {10}^{-19}}{\sqrt{6.67430\times {10}^{-11}}\times 1.67262\times {10}^{-27}}}\right)}^{{\displaystyle \frac{1}{3}}}{R}_{cS1}\\ =1.6444\times {10}^6{R}_{cS1}={\displaystyle \frac{1.23\times {10}^{-17}}{2n-1}}\left(\mathrm{m}\right)={\displaystyle \frac{1.23\times {10}^{-2}}{2n-1}}\left(\mathrm{fm}\right)\end{array}$$

(80)

${\left[{\displaystyle \frac{2(2n-1)E_p}{3G{m}_p^2}}\right]}^{{\displaystyle \frac{2}{5}}}{R}_{cS1}^{{\displaystyle \frac{2}{5}}}={\left[{\displaystyle \frac{4\sqrt{k}{q}_e}{\sqrt{G}{m}_p}}\right]}^{{\displaystyle \frac{2}{3}}}$${\displaystyle \frac{2(2n-1)E_p}{G{m}_p^2}}{R}_{cS1}={\left[{\displaystyle \frac{4\sqrt{k}{q}_e}{\sqrt{G}{m}_p}}\right]}^{{\displaystyle \frac{5}{3}}}$${\displaystyle \frac{R_{cS2}^5}{R_{cS1}^5}}={\left[{\displaystyle \frac{4\sqrt{k}{q}_e}{\sqrt{G}{m}_p}}\right]}^{{\displaystyle \frac{5}{3}}}$${\displaystyle \frac{G{m}_p^2}{2(2n-1)}}{\displaystyle \frac{R_{cS2}^5}{R_{cS1}^6}}=E_p$${\displaystyle \frac{\sqrt{k}{q}_e\sqrt{G}{m}_p}{4}}{\displaystyle \frac{R_{cS2}^2-R_{cw}^2}{R_{cS1}^3}}=E_e$

$$\begin{array}{c}{R}_{cS2}=\sqrt{{\displaystyle \frac{4E_e}{\sqrt{k}{q}_e\sqrt{G}{m}_p}}{R}_{cS1}^3+R_{cw}^2}\approx R_{cw}\left[1+{\displaystyle \frac{2E_e}{\sqrt{k}{q}_e\sqrt{G}{m}_p}}\times {\displaystyle \frac{R_{cS1}^2}{R_{cw}^2}}\times R_{cS1}\right]\\ =R_{cw}\left[1+{\displaystyle \frac{2E_e}{\sqrt{k}{q}_e\sqrt{G}{m}_p}}\times {\left({\displaystyle \frac{\sqrt{G}{m}_p}{4\sqrt{k}{q}_e}}\right)}^{{\displaystyle \frac{2}{3}}}\times {\displaystyle \frac{G{m}_p^2}{2(2n-1)E_p}}{\left({\displaystyle \frac{4\sqrt{k}{q}_e}{\sqrt{G}{m}_p}}\right)}^{{\displaystyle \frac{5}{3}}}\right]\\ =R_{cw}\left[1+{\displaystyle \frac{1}{2n-1}}{\displaystyle \frac{4E_e}{E_p}}\right]=R_{cw}\left[1+{\displaystyle \frac{1}{2n-1}}{\displaystyle \frac{4\times 510.999\times {10}^3}{938.272\times {10}^6}}\right]\\ =\left(1+{\displaystyle \frac{3.20\times {10}^{-3}}{2n-1}}\right)R_{cw}\\ =\left({\displaystyle \frac{1.23\times {10}^{-2}}{2n-1}}\right)\left(1+{\displaystyle \frac{3.20\times {10}^{-3}}{2n-1}}\right)\left(\mathrm{fm}\right)\end{array}$$

(81)

Based on Equation (75), the binding energy of a Dark II E_{D_νII} is

$$E_{D\_\nu II}\approx {\displaystyle \frac{E_p}{2}}={\displaystyle \frac{938.272}{2}}=469\left(\mathrm{MeV}\right)$$

(82)

From Equation (37), the binding energy of a Dark I E_{D_νI} is

$$\begin{array}{c}{E}_{D\_\nu I}\approx {\displaystyle \frac{2(n+2)}{5}}{E}_{D\_\nu II}={\displaystyle \frac{(n+2)E_p}{5}}=(n+2)\times {\displaystyle \frac{938.272}{5}}\\ =(n+2)\times 188(\mathrm{MeV})\end{array}$$

(83)

However, I have not found a way to calculate the energy of an electron neutrino E_{e_ν}, which also leads to an inability to determine the tiny difference between gravitational mass and inertial mass, especially when ${\widehat{m}}_e={\displaystyle \frac{c_G^2\sqrt{G}{m}_e}{c_L}}$.

It can be determined that m > n since R_cS2 ≈ R_cw and the effective range of the weak interaction is smaller than the one of the strong interaction. I would further suggest that n = 3 and m = 5 because they seem to fit the current knowledge of the strong and weak interactions, and 1, 3, 5 is a tiny, pretty odd series. Of course, the exact number of m, n can only be found experimentally.

The speed of gravity c_G is easily obtained from Equations (45) and (81)

$$c_G={\displaystyle \frac{R_{cS2}}{R_{cw}}}{c}_L=\left(1+{\displaystyle \frac{3.20\times {10}^{-3}}{2n-1}}\right)c_L$$

(84)

Moreover, the energy definition Equations (13) and (16) reveal that regardless of the speed of light c_L or the speed of gravity c_G, there is an invariant that is the rest mass m₀. Therefore, we can obtain the relation of momentum and energy between the two measurement media of light and gravitational waves

$${\displaystyle \frac{1}{c^2}}\left(p_x^2+p_y^2+p_z^2-{\displaystyle \frac{E^2}{c^2}}\right)={\displaystyle \frac{1}{{c^{\prime }}^2}}\left({p^{\prime }}_x^2+{p^{\prime }}_y^2+{p^{\prime }}_z^2-{\displaystyle \frac{{E^{\prime }}^2}{{c^{\prime }}^2}}\right)$$

(85)

from the relationship between momentum and energy in Special Relativity $p_x^2+p_y^2+p_z^2-{\displaystyle \frac{E^2}{c^2}}=-m_0^2{c}^2$. Equation (85) directly derives the relation of space-time between the two measurement media

$${\displaystyle \frac{1}{c^2}}\left(d{x}^2+d{y}^2+d{z}^2-c^{2}d{t}^2\right)={\displaystyle \frac{1}{{c^{\prime }}^2}}\left(d{x^{\prime }}^2+d{y^{\prime }}^2+d{z^{\prime }}^2-c^{2}d{t^{\prime }}^2\right)$$

(86)

So the proper time τ is the same in all reference frames with the two measurement media.