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 are not actually indivisible about a century later. An atom consists of three basic types of subatomic particles: electrons, protons, and neutrons, which occupy the tiny space in an atom. Protons and neutrons form the nucleus that contains most of an atom’s mass. Electrons are the lightest charged particles in nature and revolve around the nucleus of an atom. An electron is seemingly indivisible yet. Until today, we have not split an electron into two or more smaller particles. We only make the positive and negative electron annihilation. A free neutron is unstable, decaying into a proton, electron, and neutrino. However, a free proton is stable, and is composed of two up quarks and one down quark in the modern Standard Model. Furthermore, 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 will happen when light or electromagnetic, gravity waves and other waves make up the primary particles that constitute the fundamental elements of matter? The assumptions are derived then:

• Light or electromagnetic waves, weak interaction, gravity, and strong interaction are independent waves without rest masses. But their structures are different.
• 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)
• Weak and strong interaction 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)
• According to the unified electro-weak theory, light and weak interaction have the same speed c_L with spin number +1 or –1.
• Gravity and strong interaction have the same speed c_G without spin, and c_G is constant in a vacuum.
• The primary particles, which are electrons, electron neutrinos, and dark neutrinos in this paper, are made by the above four types of waves.

2. The Formation of Primary Particles

To suppose that the birth of primary particles may divide into the following two parts:

• The light and 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. The gravity and strong interaction couple together too (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.
• The E-W couple has no spin. The original spins of the coupled waves convert the electric or weak charge property.
• It makes a primary particle when the 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. The 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 with the same original spin.

Thus, in the spherical coordinate system, we have the uniform Laplace equation for the above 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)

whose general solution is

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

(6)

where A_j and B_j 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.

3. The Fields and Binding-Energy

We can now derive the electric, gravitational, weak interaction, and strong interaction fields Ê based on equation (6) and existing physical laws and data.

It is reasonable that equation (6) can be transformed into a pair of conjugate solutions.

$$\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}$$

(7)

Clearly, there is

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

(8)

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

And we let

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

(9)

The field Ê may be split into the macroscopic item ${\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.

Compared equation (7) 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}$$

(10)

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}$$

(11)

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.

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[1,2,3]. All of the above can help us to determine weak interaction field Ê_w from equation (7). Ê_w is so supposed to be equal to

$$\begin{array}{c}{\widehat{E}}_w=\{\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}\\ {\widehat{E}}_w^{*}=\{\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}\end{array}$$

(12)

where$\widehat{w}$is the mathematical weak charge, R_cw is the critical radius of 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 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,6]. Moreover, it is supposed that the fields of the four fundamental forces have a unified form in a very tiny range. Hence equation (7) can be translated into the strong interaction field Ê_S

$$\begin{array}{c}{\widehat{E}}_S=\{\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}\\ {\widehat{E}}_S^{*}=\{\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}\end{array}$$

(13)

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

Further, and m ≠ n.

Now we turn to determine the energy. The energy density of wave equation (1) is given by ${\left|\widehat{E}\right|}^2$, and the equations (1) and (2) are Lorentz invariance. So we hope the equation of energy is Lorentz invariance, too, and let the binding-energy E_a-b of the four 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|dxdydzdct}}$$

(14)

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

$${\displaystyle \underset{0}{\overset{2\pi }{\int }}{e}^{im\phi }{e}^{-in\phi }d\phi }=2\pi {\delta }_{mn}$$

(15)

and

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

(16)

When we use light as a measurement medium to determine the energy as usual, we can translate equation (14) 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}$$

(17)

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

The general energy expression can be calculated when equation (7) is substituted into equation (17).

$$\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 }}[\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 }}}}}\\ \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 }}]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 }}[\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 }}}}}\\ \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 }}]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}$$

(18)

Based on equation (18) and associated with equations (10) to (13), 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}}$$

(19)

$$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}}$$

(20)

$$E_{w-w}=\{\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}$$

(21)

and

$$E_{S-S}=\{\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}\hfill \\ ={\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}\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}$$

(22)

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}}$$

(23)

$$E_{e-w}=\{\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^4{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}$$

(24)

$$E_{e-S}=\{\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}$$

(25)

$$E_{G-w}=\{\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}$$

(26)

$$E_{G-S}=\{\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}$$

(27)

and

$$E_{w-S}=\{\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}$$

(28)

For the convenience of subsequent calculations, we require the results of our defined energy to be consistent with those of the conventional physical method. Compared equations (19) and (20) to electric and gravitational potential energy formulas, it can easily find that the relations of mathematical electric charge ê and mathematical mass $\widehat{m}$ to electric charge q and mass m are

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

(29)

and

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

(30)

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

4. The Structures of Primary Particles

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

• The whole biding-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 envelopes, which means the size of a primary particle, too, depends on the critical radius of weak or strong interaction.
• The two envelopes have the same inherent frequency ν_in, although this is not mathematically required.
• The degree of associated Legendre polynomials j is the same on the two envelopes.
• 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 biding-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}$$

(31)

based on equation (17) E_pri is clearly equivalent to the rest mass.

Evidently, the total energy of a primary particle comprises the biding-energy or the rest mass and the energy in static fields.

4.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 weak interaction R_cw.

The charges in equations (10) and (12) are equal and minimal for an electron neutrino, which means ${ê}_{e\_\nu }={\widehat{w}}_{e\_\nu }$ when ${ê}_{e\_\nu }$ and ${\widehat{w}}_{e\_\nu }$ are the mathematical electric charge and the mathematical weak charge of an electron neutrino.

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

$${\widehat{E}}_{E,e\_\nu }\{\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}$$

(32)

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

Based on equation (31) and associated equations (19), (21), and (24), we can easily compute the biding-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}$$

(33)

Combining the envelopes’ characters b. with equations (10), (12), and (32), the fields around an electron neutrino Ê_{e_ν} can be directly written as

$${{\widehat{E}}_{e\_\nu }|}_{r>R_{cw}}\{\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}r>R_{cw}$$

(34)

4.2. 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 radius of strong interaction.
• Dark I and II have the same mathematical mass and mathematical strong charge.
• The mathematical mass ${\widehat{m}}_{D\_\nu }$ and the mathematical strong charge ${\widehat{s}}_{D\_\nu }$ are equal, i.e., ${\widehat{m}}_{D\_\nu }={\widehat{s}}_{D\_\nu }$. ${\widehat{s}}_{D\_\nu }$ is minimal.

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}\{\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}$$

(35)

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

Based on equation (31) and associated equations (20), (22), and (27), we can compute the biding-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}$$

(36)

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

$${\widehat{E}}_{E,D\_\nu II}\{\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}$$

(37)

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}$$

(38)

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

$${{\widehat{E}}_{D\_\nu I}|}_{r>R_{cS1}}\{\begin{array}{l}\begin{array}{c}\left[{\widehat{E}}_{G-D\_\nu I}+{\widehat{E}}_{S-D\_\nu I}\right]=\{\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}\end{array}\hfill \\ {\left[{\widehat{E}}_{G-D\_\nu I}+{\widehat{E}}_{S-D\_\nu I}\right]}^{*}=\{\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}\hfill \end{array}$$

(39)

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

$${{\widehat{E}}_{D\_\nu II}|}_{r>R_{cS2}}\{\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}r>R_{cS2}$$

(40)

Comparing equations (36) with (38) 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.

4.3. An Electron or A Positron

Electrons and positrons have the same structure. We will not distinguish significantly between electrons and positrons during 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 weak interaction R_cw, although there are three critical radii for 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$ are minimal, which means ${\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 biding-energy of an electron E_e.

$${\widehat{E}}_{E,e}\{\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}$$

(41)

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}$$

(42)

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 }$, 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}}$$

(43)

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{q_e^2}{c_L^2{R}_{cw}}}\left({\displaystyle \frac{R_{cS2}^2}{R_{cw}^2}}-1\right)$$

(44)

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

$${{\widehat{E}}_e|}_{r>R_{cw}}\{\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]}^{*}=\{\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}\hfill \end{array}$$

(45)

Next, the examination of equations (41) and (42) reveals that the 2^nd critical radius of the strong interaction R_cS2 should be the geometric characterization parameter of a 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}}$$

(46)

which shows that the gravity speed c_G is faster than the light speed c_L when the above equation compares with equation (42).

5. 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 for a distance of l. Particles I and II have the potential and kinetic energy in the separated state (Figure 2). At 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 in envelopes of primary particles.

Based on the law of conservation of energy, the initial energy of particles I and II is equal to the sum of the potential energy and the kinetic energy (including the magnetic energy for electrons) after their separation. Therefore, referring to equation (17) and taking the r_I coordinate system of the particle I as a reference (Figure 2), the potential energy of the two primary particles $E_{P_1^P_2}$ can be defined as

$$\begin{array}{c}{E}_{P_1^P_2}={\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}(r_I){\widehat{E}}_{P_2}^{*}(r_I)\right|r_I^2\sin {\theta }_{I}d{r}_{I}d{\theta }_{I}d{\phi }_{I}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}(r_I){\widehat{E}}_{P_2}^{*}(r_I)\right|r_I^2\sin {\theta }_{I}d{r}_{I}d{\theta }_{I}d{\phi }_{I}d{c}_{L}t}}}}\right)\\ ={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^k}}}\left({\displaystyle \underset{R_c}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left|{\widehat{E}}_{P_1}(r_I){\widehat{E}}_{P_2}^{*}(r_I)\right|r_I^2\sin {\theta }_{I}d{r}_{I}d{\theta }_{I}d{\phi }_I}}}-{\displaystyle \underset{R_c}{\overset{l}{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left|{\widehat{E}}_{P_1}(r_I){\widehat{E}}_{P_2}^{*}(r_{II})\right|r_I^2\sin {\theta }_{I}d{r}_{I}d{\theta }_{I}d{\phi }_I}}}\right)\end{array}$$

(47)

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 }}\left|{\widehat{E}}_{P_1}(r_I){\widehat{E}}_{P_2}^{*}(r_{II})\right|r_I^2\sin {\theta }_{I}d{r}_{I}d{\theta }_{I}d{\phi }_I}}}$ converts to kinetic energy (including the magnetic energy for electrons).

It is therefore assumed that the Self-Conjugate Mechanism remains between two interacting primary particles. According to the existing knowledge of physics, 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 particles have achieved the self-conjugacy after rotating one cycle with angular velocity ω. There are self-conjugate lines in pairs, which are further presumed at $I:(r_c,\theta ,\omega t_c)$ for particle I and $II:(r_c,\theta ,\omega t_c+\pi )$ for particle II in their respective spherical coordinate systems, i.e., these lines are in the opposite position while r_I = r_II = r, θ_I = θ_II = θ (Figure 2). Thus, at the self-conjugate lines and using equation (7), equation (47) could 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}}}\left({\displaystyle \underset{R_c}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left|{\widehat{E}}_{P_1}(r_I){\widehat{E}}_{P_2}^{*}(r_I)\right|r_I^2\sin {\theta }_{I}d{r}_{I}d{\theta }_{I}d{\phi }_I}}}-{\displaystyle \underset{R_c}{\overset{l}{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left|{\widehat{E}}_{P_1}(r_I){\widehat{E}}_{P_2}^{*}(r_I)\right|r_I^2\sin {\theta }_{I}d{r}_{I}d{\theta }_{I}d{\phi }_I}}}\right)\\ ={\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}(r){\widehat{E}}_{P_2}^{*}(r)\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{l}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}{\displaystyle \underset{0}{\overset{2\pi }{\int }}\left[{\overline{\widehat{E}}}_{P_1}(r){\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )e^{-ik\phi }}\right]\left[{\overline{\widehat{E}}}_{P_2}^{*}(r){\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{l}{\overset{\infty }{\int }}{\displaystyle \underset{0}{\overset{\pi }{\int }}\left[{\overline{\widehat{E}}}_{P_1}(r){\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )}\right]\left[{\overline{\widehat{E}}}_{P_2}^{*}(r){\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )}\right]r^2\sin \theta drd\theta {\displaystyle \underset{2\eta \pi }{\overset{2\eta \pi +2\pi }{\int }}{e}^{im(\omega t_c)}{e}^{-in(\omega t_c+\pi )}d(\omega 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 }}\left[{\overline{\widehat{E}}}_{P_1}(r){\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )}\right]\left[{\overline{\widehat{E}}}_{P_2}^{*}(r){\displaystyle \sum _{k=0}^j{\mathrm{P}}_j^k(\cos \theta )}\right]r^2\sin \theta drd\theta \times }}\\ \{\begin{array}{l}{\displaystyle \underset{2\eta \pi }{\overset{2\eta \pi +2\pi }{\int }}{e}^{im(\omega t_c)}{e}^{-in(\omega t_c+\pi )}d(\omega t_c)}=-2\pi {\delta }_{mn}\mathrm{Same}\mathrm{rotation}\mathrm{direction}\mathrm{of}\omega \hfill \\ {\displaystyle \underset{2\eta \pi }{\overset{2\eta \pi +2\pi }{\int }}{e}^{im(\omega t_c)}{e}^{-in(-\omega t_c+\pi )}d(\omega t_c)}=0\mathrm{Opposite}\mathrm{rotation}\mathrm{direction}\mathrm{of}\omega \hfill \end{array}\\ =\{\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{Same}\mathrm{rotation}\mathrm{direction}\mathrm{of}\omega \hfill \\ 0\mathrm{Opposite}\mathrm{rotation}\mathrm{direction}\mathrm{of}\omega \hfill \end{array}\end{array}$$

(48)

where t_c is the time of conjugation of two primary particles, and η = 0, 1, 2, …, η. Here, the minus sign only indicates that the energy spreads out in this process. In other words, equation (48) should have the plus sign when two primary particles move in the opposite mode of Figure 2, i.e., two rested on each other particles are attracted at the initial distance l and approach each other, since the r_II coordinate system is reversed in Figure 2, which means that the energy collects. Later, we will continue to use equation (48) and ignore the minus sign, without distinguishing whether the energy is dissipated or accumulated in the motion of two primary particles.

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 forms between the two particles because each of the two particles consists of two coupled waves. Therefore, equation (48) can be translated into

$$\begin{array}{c}{E}_{P_1^P_2}={\displaystyle \frac{c_L^2}{2\pi {\displaystyle \sum _{k=0}^j{\mathrm{N}}_j^{\left|k\right|}}}}{\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^{\left|k\right|}}}}\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}$$

(49)

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 the same as the potential energy when primary particles I and II move relative to each other. Reversing the Newtonian mechanics definition of work, i.e., $F=\nabla W=\nabla E_{P_1^P_2}$, 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}}$$

(50)

Two self-conjugate primary particles have no initial phase difference in zenith and azimuthal angle under the requirements of the self-conjugate lines assumption. From equation (48), they have potential energy or force when they rotate in the same direction, while they have zero potential energy or rest when they do in opposite directions. Therefore, the force mode has two forms, all right-handed or spin-up and all left-handed or spin-down, which shows that each primary particle has spin values of $\pm {\displaystyle \frac{1}{2}}$.

5.1. Two Particles of the Same Type

Start by computing the interaction between two electron neutrinos. Combining equations (19), (21), and (49), 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})}}+{\displaystyle \frac{R_{cw}^{2m-2}{\widehat{w}}_{e\_\nu }^2}{c_L^2(2m-1){(l+R_{cw})}^{2m-1}}}\right]\\ \approx {\displaystyle \frac{2{\widehat{q}}_{e\_\nu }^2}{c_L^2(l+R_{cw})}}={\displaystyle \frac{2{\widehat{w}}_{e\_\nu }^2}{c_L^2(l+R_{cw})}}\mathrm{when}l\ge R_{cw}\end{array}$$

(51)

Since two equations (34) have attractive and repulsive states under the Self-Conjugate Mechanism, from equation (50), 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-2}}}\right]\\ \approx \pm {\displaystyle \frac{2{\widehat{q}}_{e\_\nu }^2}{c_L^2{(l+R_{cw})}^2}}=\pm {\displaystyle \frac{2{\widehat{w}}_{e\_\nu }^2}{c_L^2{(l+R_{cw})}^2}}\mathrm{when}l\ge R_{cw}\end{array}$$

(52)

where the sign ‘–’ or ‘+’ depends on the self-conjugate forms between two electron neutrinos, and the “–” or “+” should be random.

Association equation (49) with equations (20), (22), (27), and (39), 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}}\\ =\{\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]\\ \approx {\displaystyle \frac{2c_L^2{\widehat{s}}_{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]\\ ={\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}}}+{\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{s}}_{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]\\ ={\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}\end{array}$$

(53)

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

$$F_{D\_\nu I^D\_\nu I}=\{\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]\\ \approx -{\displaystyle \frac{2c_L^2{\widehat{s}}_{D\_\nu }^2}{c_G^4{R}_{cS1}^6}}{(l+R_{cS1})}^4=-{\displaystyle \frac{2c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{R}_{cS1}^6}}{(l+R_{cS1})}^{4}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}}+{\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{s}}_{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]\\ =-{\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}$$

(54)

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}}}+{\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{s}}_{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]\\ ={\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}$$

(55)

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{s}}_{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]\\ =-{\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}$$

(56)

The potential energy between two electrons E_e^e has two forms because an electron is composed of one E-W couple and one G-S couple. Same as electron neutrinos, the two forms should be random and rely on the self-conjugate forms between two electrons. Combining equations (19) to (28) and (49), we can compute 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.

$$\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}}\\ =\{\begin{array}{l}\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{\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}\\ \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})}}\mathrm{when}l\mathrm{is}\mathrm{large}\end{array}$$

(57)

Another is E_e^eII when the two E-W couples are conjugate to the two G-S couples.

$$\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}}\\ =\{\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\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}\end{array}$$

(58)

Derivation of the last two equations can yield the forces between two electrons in both forms that are

$$\begin{array}{c}{F}_{e^eI}=\{\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}\\ \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}}\mathrm{when}l\mathrm{is}\mathrm{large}\end{array}$$

(59)

and

$$F_{e^eII}=\{\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}$$

(60)

plus ${\widehat{s}}_e={\widehat{s}}_{D\_\nu }$ in equations (57) to (60).

Comparing equations (57) and (58), (59) and (60) shows that the potential energy and force between two electrons are very different in the two self-conjugate forms, with the smaller one close to zero.

5.2. Two Particles of the Different Type

Similar to the last section, this section still starts by computing the interaction of an electron neutrino with another primary particle. Combining equations (23) to (28), (49), and (50), 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}}\\ =\{\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\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}\end{array}$$

(61)

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

$$E_{e\_\nu ^D\_\nu I}=\{\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}$$

(62)

plus ${\widehat{m}}_{D\_\nu }={\widehat{s}}_{D\_\nu }$ in above two equations.

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})}}={\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{s}}_{D\_\nu }}{c_G^2(l+R_{cS2})}}\end{array}$$

(63)

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}}=-{\displaystyle \frac{2{\widehat{q}}_{e\_\nu }{\widehat{s}}_{D\_\nu }}{c_G^2{(l+R_{cS2})}^2}}$$

(64)

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}}\\ =\{\begin{array}{l}\begin{array}{l}{\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{q}}_e}{c_L^2(l+R_{cw})}}+{\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\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})}}+{\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})}}l>R_{cS2}-R_{cw}\hfill \end{array}\\ \approx \{\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{s}}_{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}\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}\\ \approx {\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{q}}_e}{c_L^2(l+R_{cw})}}={\displaystyle \frac{{\widehat{w}}_{e\_\nu }{\widehat{w}}_e}{c_L^2(l+R_{cw})}}\mathrm{when}l\ge R_{cw}\end{array}$$

(65)

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

$$\begin{array}{c}{F}_{e\_\nu ^e}=\{\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}\\ \approx \{\begin{array}{l}{\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{s}}_{D\_\nu }}{c_G^2{R}_{cS1}^3}}(l+R_{cw})0\le l\le R_{cS2}-R_{cw}\hfill \\ {\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}\\ \approx \pm {\displaystyle \frac{{\widehat{q}}_{e\_\nu }{\widehat{q}}_e}{c_L^2{(l+R_{cw})}^2}}=\pm {\displaystyle \frac{{\widehat{w}}_{e\_\nu }{\widehat{w}}_e}{c_L^2{(l+R_{cw})}^2}}\mathrm{when}l\ge R_{cw}\end{array}$$

(66)

where ${\widehat{q}}_{e\_\nu }={\widehat{w}}_{e\_\nu }$, ${\widehat{q}}_e={\widehat{w}}_e$, and ${\widehat{s}}_e={\widehat{s}}_{D\_\nu }$ in equations (65) and (66), and the sign ‘–’ or ‘+’ randomizes in equation (66).

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 (20), (22), (27), (49), and (50), 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})}}+{\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{s}}_{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]={\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}$$

(67)

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

$$F_{D\_\nu I^D\_\nu II}=-{\displaystyle \frac{2c_L^2{\widehat{s}}_{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]=-{\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]$$

(68)

Based on equations (23) to (28), (43), (49), and (50), 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}}\\ =\{\begin{array}{l}\begin{array}{l}{\displaystyle \frac{{\widehat{q}}_e{\widehat{m}}_{D\_\nu }}{c_G^2(l+R_{cw})}}+{\displaystyle \frac{{\widehat{q}}_e{\widehat{s}}_{D\_\nu }}{2c_G^2{R}_{cS1}^3}}\left(R_{cS2}^2-{(l+R_{cw})}^2\right)+{\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})}}+{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4(l+R_{cw})}}+{\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}\\ \approx \{\begin{array}{l}\begin{array}{l}{\displaystyle \frac{c_L^2{\widehat{s}}_{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]\\ ={\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}\end{array}\hfill \\ \begin{array}{l}{\displaystyle \frac{c_L^2{\widehat{s}}_{D\_\nu }^2}{c_G^4(l+R_{cw})}}\left[{\displaystyle \frac{R_{cw}^3}{2R_{cS1}^3}}+{\displaystyle \frac{R_{cS2}^{2n+4}}{(2n-1)R_{cS1}^6{(l+R_{cw})}^{2n-2}}}\right]\\ ={\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4(l+R_{cw})}}\left[{\displaystyle \frac{R_{cw}^3}{2R_{cS1}^3}}+{\displaystyle \frac{R_{cS2}^{2n+4}}{(2n-1)R_{cS1}^6{(l+R_{cw})}^{2n-2}}}\right]l>R_{cS2}-R_{cw}\end{array}\hfill \end{array}\end{array}$$

(69)

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

$$\begin{array}{c}{F}_{D\_\nu I^e}=\{\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}\\ \approx \{\begin{array}{l}-{\displaystyle \frac{c_L^2{\widehat{s}}_{D\_\nu }^2}{c_G^4{R}_{cS1}^6}}{(l+R_{cw})}^4={\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 \\ \begin{array}{l}-{\displaystyle \frac{c_L^2{\widehat{s}}_{D\_\nu }^2}{c_G^4{(l+R_{cw})}^2}}\left[{\displaystyle \frac{R_{cw}^3}{2R_{cS1}^3}}+{\displaystyle \frac{R_{cS2}^{2n+4}}{R_{cS1}^6{(l+R_{cw})}^{2n-2}}}\right]\\ =-{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{(l+R_{cw})}^2}}\left[{\displaystyle \frac{R_{cw}^3}{2R_{cS1}^3}}+{\displaystyle \frac{R_{cS2}^{2n+4}}{R_{cS1}^6{(l+R_{cw})}^{2n-2}}}\right]l>R_{cS2}-R_{cw}\end{array}\hfill \end{array}\end{array}$$

(70)

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})}}+{\displaystyle \frac{c_L^2{\widehat{m}}_e{\widehat{m}}_{D\_\nu }}{c_G^4(l+R_{cS2})}}+{\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}}}\\ \approx {\displaystyle \frac{c_L^2{\widehat{s}}_{D\_\nu }^2}{c_G^4(l+R_{cS2})}}\left[{\displaystyle \frac{R_{cw}^3}{2R_{cS1}^3}}+{\displaystyle \frac{R_{cS2}^{2n+4}}{(2n-1)R_{cS1}^6{(l+R_{cS2})}^{2n-2}}}\right]\\ ={\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4(l+R_{cS2})}}\left[{\displaystyle \frac{R_{cw}^3}{2R_{cS1}^3}}+{\displaystyle \frac{R_{cS2}^{2n+4}}{(2n-1)R_{cS1}^6{(l+R_{cS2})}^{2n-2}}}\right]\end{array}$$

(71)

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

$$\begin{array}{c}{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}}}\\ \approx -{\displaystyle \frac{c_L^2{\widehat{s}}_{D\_\nu }^2}{c_G^4{(l+R_{cS2})}^2}}\left[{\displaystyle \frac{R_{cw}^3}{2R_{cS1}^3}}+{\displaystyle \frac{R_{cS2}^{2n+4}}{R_{cS1}^6{(l+R_{cS2})}^{2n-2}}}\right]\\ =-{\displaystyle \frac{c_L^2{\widehat{m}}_{D\_\nu }^2}{c_G^4{(l+R_{cS2})}^2}}\left[{\displaystyle \frac{R_{cw}^3}{2R_{cS1}^3}}+{\displaystyle \frac{R_{cS2}^{2n+4}}{R_{cS1}^6{(l+R_{cS2})}^{2n-2}}}\right]\end{array}$$

(72)

6. The Structure Values of Primary Particles

In this chapter, we attempt to import the existing physical data into the computational results of the previous chapters 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 macro and micro levels. We can be confident that the available measurements reflect the characteristics of the electric charge and no other factors. Thus, comparing equation (45) to Gauss’s law of electrostatics and combined equation (29), it can be easily found that the relation of the mathematical electric charge ${\widehat{q}}_e$ of an electron to an electron charge q_e in present physical data.

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

(73)

It is somewhat difficult to obtain the charges of gravity and the strong force because gravity is the weakest force in the four fundamental interactions, and the strong force is a short-range interaction. In the models of this paper, although the gravitational and inertial mass are different, it can be hypothesized that Newton’s law of gravity is accurate to a large extent because only electron neutrinos in the model do not have gravitational waves and because electron neutrinos have much less biding-energy or rest mass. We focus on protons ---- stable, heavy subatomic particles ---- to explore these charges and start with the components of a proton. A combination of two attractive primary particles is easily separated by external action since there is no third particle to hold it back. However, the aggregation of multiple strongly attracted primary particles can cause the annihilation of these particles. Therefore, a combination of three high biding-energy, mutually attractive primary particles should be an extremely stable particle. Similar to the Quark model[7], it is supposed that a proton is composed of two Dark IIs and one positron instead of two up quarks and one down quark and has the structure of Figure 3 that will not delved into.

From equations (55), (56), (71), and (72), a proton might have the structure of Figure 3, and the mutual distances between the three primary particles are small, i.e., $l\approx 0$, in a proton. It is further supposed that the kinetic energies of the three primary particles’ mutual motion (including magnetic energy) are negligibly small compared to their binding-energies and potential energies. Hence, the energy of a proton E_p equals

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

(74)

when equations (38), (55), and (71) are subsequently substituted in 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,8]. Furthermore,

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

(75)

can be obtained from equations (30) and (40), where m_p is the mass of a proton.

Next, by combining equations (43), (44), and (73) to (75), we have

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

(76)

that we can solve and yield [8]

$$\begin{array}{l}{\left[{\displaystyle \frac{8(2n-1)E_p}{3G{m}_p^2}}\right]}^{{\displaystyle \frac{2}{5}}}{R}_{cS1}^{{\displaystyle \frac{2}{5}}}-{\displaystyle \frac{8E_e}{\sqrt{k}{q}_e\sqrt{G}{m}_p}}{R}_{cS1}={\left[{\displaystyle \frac{8\sqrt{k}{q}_e}{\sqrt{G}{m}_p}}\right]}^{{\displaystyle \frac{2}{3}}}\\ {\left[{\displaystyle \frac{8(2n-1)\times 938.272\times {10}^6\times 1.60218\times {10}^{-19}}{3\times 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{8\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{8\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}}}-5.8394\times {10}^{11}{R}_{cS1}=7.9425\times {10}^{-10}\end{array}$$

(77)

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 7.9425\times {10}^{-10}\\ R_{cS1}={\displaystyle \frac{1.78\times {10}^{-23}}{2n-1}}\left(\mathrm{m}\right)={\displaystyle \frac{1.78\times {10}^{-8}}{2n-1}}\left(\mathrm{fm}\right)\end{array}$$

(78)

Next

$$\begin{array}{c}{R}_{cw}={\left({\displaystyle \frac{8\sqrt{k}{q}_e}{\sqrt{G}{m}_p}}\right)}^{{\displaystyle \frac{1}{3}}}{R}_{cS1}={\left({\displaystyle \frac{8\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}\\ =2.0718\times {10}^6{R}_{cS1}={\displaystyle \frac{3.68\times {10}^{-17}}{2n-1}}\left(\mathrm{m}\right)={\displaystyle \frac{3.68\times {10}^{-2}}{2n-1}}\left(\mathrm{fm}\right)\end{array}$$

(79)

Since in the solution of R_cS1 from equation (76), we actually make ${\left[{\displaystyle \frac{8(2n-1)E_p}{3G{m}_p^2}}\right]}^{{\displaystyle \frac{2}{5}}}{R}_{cS1}^{{\displaystyle \frac{2}{5}}}={\left[{\displaystyle \frac{8\sqrt{k}{q}_e}{\sqrt{G}{m}_p}}\right]}^{{\displaystyle \frac{2}{3}}}$, i.e., ${\displaystyle \frac{8(2n-1)E_p}{3G{m}_p^2}}{R}_{cS1}={\left[{\displaystyle \frac{8\sqrt{k}{q}_e}{\sqrt{G}{m}_p}}\right]}^{{\displaystyle \frac{5}{3}}}$, which can derive ${\displaystyle \frac{R_{cS2}^5}{R_{cS1}^5}}={\left[{\displaystyle \frac{8\sqrt{k}{q}_e}{\sqrt{G}{m}_p}}\right]}^{{\displaystyle \frac{5}{3}}}$ from ${\displaystyle \frac{3G{m}_p^2}{8(2n-1)}}{\displaystyle \frac{R_{cS2}^5}{R_{cS1}^6}}=E_p$ and tell us that R_cS2 and R_cw are almost equal. This computational procedure is not precise enough, so we use${\displaystyle \frac{\sqrt{k}{q}_e\sqrt{G}{m}_p}{8}}{\displaystyle \frac{R_{cS2}^2-R_{cw}^2}{R_{cS1}^3}}=E_e$ in equation (76) to solve R_cS2

$$\begin{array}{c}{R}_{cS2}=\sqrt{{\displaystyle \frac{8E_e}{\sqrt{k}{q}_e\sqrt{G}{m}_p}}{R}_{cS1}^3+R_{cw}^2}\approx R_{cw}\left[1+{\displaystyle \frac{4E_e}{\sqrt{k}{q}_e\sqrt{G}{m}_p}}{\displaystyle \frac{R_{cS1}^2}{R_{cw}^2}}{R}_{cS1}\right]\\ =R_{cw}\left[1+{\displaystyle \frac{4E_e}{\sqrt{k}{q}_e\sqrt{G}{m}_p}}{\left({\displaystyle \frac{\sqrt{G}{m}_p}{8\sqrt{k}{q}_e}}\right)}^{{\displaystyle \frac{2}{3}}}{\displaystyle \frac{3G{m}_p^2}{8(2n-1)E_p}}{\left({\displaystyle \frac{8\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{12E_e}{E_p}}\right]=R_{cw}\left[1+{\displaystyle \frac{1}{2n-1}}{\displaystyle \frac{12\times 510.999\times {10}^3}{938.272\times {10}^6}}\right]\\ =\left(1+{\displaystyle \frac{6.54\times {10}^{-3}}{2n-1}}\right)R_{cw}\\ =\left({\displaystyle \frac{3.68\times {10}^{-2}}{2n-1}}\right)\left(1+{\displaystyle \frac{6.54\times {10}^{-3}}{2n-1}}\right)\left(\mathrm{fm}\right)\end{array}$$

(80)

Based on equation (74) the biding-energy of a Dark II E_{D_νII} is

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

(81)

From equation (38), the biding-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}{15}}=(n+2)\times {\displaystyle \frac{938.272}{15}}\\ =(n+2)\times 62.6(\mathrm{MeV})\end{array}$$

(82)

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 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 strong interaction. I would suggest further that n = 3 and m = 5 because they seem to fit the current knowledge of strong and weak interactions, and 1, 3, 5 is a tiny, pretty odd series. Of course, exactly how many m, n can only be obtained experimentally.

The gravity speed c_G can easily get from equations (46) and (80)

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

(83)

In addition, the energy definition equations (14) and (17) reveal that regardless of the light speed c_L or the gravity speed 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)$$

(84)

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 (84) 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)$$

(85)

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

7. Conclusions and Discussion

In this paper, field and energy equations for strong and weak interactions are derived based on a set of fundamental assumptions that are consistent with existing knowledge. With the help of these equations, the structures of the primary particles are established, and the interactions between two primary particles are analyzed. The characteristic parameters of the primary particles (Table 1) are computed from the present physical data.

Moreover, we can derive the matrix $\left(\begin{array}{cc}\cos \theta & \sin \theta e^{-i\phi }\\ \sin \theta e^{i\phi }& -\cos \theta \end{array}\right)$ when we rewrite the quantum factors in equations (32), (35), (37), and (41) as $\left(\begin{array}{cc}\cos \theta & \sin \theta e^{-i\phi }\\ \sin \theta e^{i\phi }& \cos \theta \end{array}\right)$ and multiply it by the Pauli matrix $\left(\begin{array}{cc}1& 0\\ 0& -1\end{array}\right)$. $\left(\begin{array}{cc}\cos \theta & \sin \theta e^{-i\phi }\\ \sin \theta e^{i\phi }& -\cos \theta \end{array}\right)$ has well-known results in quantum mechanics with eigenvalues of ±1 and eigenvectors of $\left(\begin{array}{c}\cos {\displaystyle \frac{\theta }{2}}\\ \sin {\displaystyle \frac{\theta }{2}}{e}^{i\phi }\end{array}\right)$, $\left(\begin{array}{c}\sin {\displaystyle \frac{\theta }{2}}{e}^{-i\phi }\\ -\cos {\displaystyle \frac{\theta }{2}}\end{array}\right)$, i.e., spin values $\pm {\displaystyle \frac{1}{2}}$. Thus the conclusions of this paper are compatible with existing physical-mathematical methods.

Furthermore, the interactive analysis reveals that:

• The Self-Conjugation condition is that two primary particles have no initial phase difference in zenith and azimuthal angle. Two Self-Conjugation primary particles have potential energy or force when they rotate in the same direction. However, they have zero potential energy or rest when they spin in opposite directions. It is one of the foundations of the Pauli exclusion principle.
• Dark Is have the asymptotic freedom characteristic, but following the principle of energy minimization, there should only be Dark IIs in most cases.
• The force between two electrons has three values, one large, one small, and one zero.
• Whether two electron neutrinos or an electron neutrino and an electron attract or repel each other is randomized. Because of this, electron neutrinos are a weak destabilizer in the nucleus, and even though the binding-energy of electron neutrinos is the smallest, no signs of neutrino destruction have been found so far.
• Primary particles behave perfect tiny spheres in terms of energy and interactions, but they also look like uneven minuscule spheres in external fields. Which is the reality of a primary particle? Observation or mathematics? The answer should be that „the Moon is always there, doesn’t matter we see it or not”, however, the Moon is changed when we see it.

In the interactive analysis, the Self-Conjugate Mechanism plays a significant role in the microscopic domain. Consider the case when external actions cause a combination of two primary particles to rotate from the same to different directions until they are opposite and separated. So the two primary particles have the opposite direction of spin, or in other words, they entangle themselves after the separation. Evidently, the distance of quantum entanglement is finite, not infinite like currently believed, because the strength of fields around primary particles decreases with distance. The transmission speed of quantum entanglement is the speed of light or gravity. Einstein is right in this case.

In addition, the gravity speed c_G is slightly larger than the light speed c_L $\left({\displaystyle \frac{c_G}{c_L}}=1+{\displaystyle \frac{6.54\times {10}^{-3}}{2n-1}}\right)$, and the rest mass m₀, the proper time τ are the invariants with the two measurement media of light and gravity. Gravity might have similar quantum characteristics, such as $E_G=h_G{\nu }_G$. What are the quantum characteristics of weak and strong interactions?

It can be further hypothesized that a dark matter is formed when the positron in a proton is replaced by one Dark II. Therefore, the energy of a dark matter E_dm is

$$\begin{array}{c}{E}_{dm}=3E_{D\_\nu II}+3{\left[E_{D\_\nu II^D\_\nu II}\right]}_{l=0}\\ ={\displaystyle \frac{3c_L^2{R}_{cS2}^5{\widehat{m}}_{D\_\nu }^2}{(2n-1)c_G^4{R}_{cS1}^6}}+{\displaystyle \frac{6c_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]}_{l=0}\\ \approx {\displaystyle \frac{9c_L^2{R}_{cS2}^5{\widehat{m}}_{D\_\nu }^2}{(2n-1)c_G^4{R}_{cS1}^6}}\approx 9E_{D\_\nu II}={\displaystyle \frac{3}{2}}{E}_p\end{array}$$

(86)

and the mathematical mass of a dark matter ${\widehat{m}}_{dm}$ is

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

(87)

from equations (38) and (55). Equations (86) and (87) remind us again that the difference between gravitational and inertial mass is almost negligible. Certainly, it is also reasonable to assume that a combination of three Dark IIs was annihilated by itself, and there are only free-state Dark IIs as dark matters in our cosmic, since no particle heavier than a proton is as stable as a proton. Or, are we listening to a concerto for solo, triple Dark IIs, and conventional matters in our universe today?

Let us imagine now what will happen when all primary particles in our universe will gather, be crushed, till be destroyed at a point. This consequence should be similar to electron annihilation. There will only be the four fundamental waves with vast high energy in our universe at that moment, namely the scene of the Big Bang.

It is worth noting that this paper cannot answer:

• What is the angular velocity ω of a primary particle, and is ω stable or variable?
• Why positrons and negative protons are so rare in our universe, implying that Dark IIs confine only positrons.