When an electromagnetic wave propagate in a conductive medium, the scalar potential φ of the electromagnetic field is related to the distribution density ρ of free charges in the medium, and the magnetic vector potential $A$ has relation to the induced electric field. The definition formulas of φ and $A$ are as follow $-\nabla \phi -\partial A/\partial t=E$ and $\nabla \times A=B$. From Gauss’s law and Lorentz gauge [21] , it can be obtained that

${\nabla }^2\phi -\mu \epsilon \frac{{\partial }^2\phi }{\partial t^2}=-\rho /\epsilon$ (1)

which is called as d’Alembert’s equation of the scalar potential in the electromagnetic field [21] .

As shown in the Figure 1 , there is a single and static charge $q\left(r^{\prime },t\right)$ or some continuous, static and time-varying charges $\rho \left(r^{\prime },t\right)$ in the source region $V^{\prime }$ having a coordinate system $P^{\prime }\left(R^{\prime },{\theta }^{\prime },z^{\prime }\right)$ in the conductive medium. The field point region is V having a coordinate system $P\left(R,\theta ,z\right)$. O is the origin point of the coordinate system. The source point $P^{\prime }$ is in a stationary state.

In Figure 1 , $r$ and $r^{\prime }$ are the vectors from the origin point O to P and $P^{\prime }$; $R$ is the vector from $P^{\prime }$ to P, which is equal to $r-r^{\prime }$; $R$ is the magnitude of the vector $R$; $e_R$ is the unit vector of $R$, which is equal to $R/R$.

The solution of d’Alembert’s equation for φ generated by a single time-varying

charge $q\left(r^{\prime },t\right)$ is $\frac{q\left(r^{\prime },t-\frac{R}{C}\right)}{4\pi \epsilon R}$. By using the superposition principle [21] , the

scalar potential φ caused by $\rho \left(r^{\prime },t\right)$ can be expressed as

$\phi \left(R,t\right)=\frac{1}{4\pi \epsilon }{\displaystyle {\int }_{V^{\prime }}\frac{\rho \left(r^{\prime },t-\frac{R}{C}\right)}{R}\text{d}{V}^{\prime }}$, (2)

in which $t-R/C$ is the delay potential item [21] . Equation (2) is the general solution of Equation (1).

In Figure 1 the continuous and time-varying charges $\rho \left(r^{\prime },t\right)$ in $V^{\prime }$ are stationary and a time-varying and waving current source $J\left(r^{\prime },t\right)$ is excited by external electromagnetic field. At the same time assuming that the LEM wave travels along the positive z-axis in the field point region $V$ with a wave velocity $V_{ez}$ and a frequency ${\omega }_{\begin{array}{c}P\end{array}}$ gives that $\partial r^{\prime }/\partial t=0$, $\partial \rho \left(r^{\prime },t\right)/\partial t\ne 0$ and $\partial R/\partial t=\left(\partial R/\partial z\right)\cdot \left(\partial z/\partial t\right)=V_{ez}$. From Lorentz gauge and the law of charge conservation [21] , it is obtained that

$\begin{array}{c}\nabla \cdot A=-\mu \epsilon \frac{\partial \phi }{\partial t}=-\mu \epsilon \frac{1}{4\pi \epsilon }{\displaystyle {\int }_{V^{\prime }}\left[\frac{1}{R}\frac{\partial \rho \left(r^{\prime },t-\frac{R}{C}\right)}{\partial t}-\frac{1}{R^2}\frac{\partial R}{\partial t}\rho \right]\text{d}{V}^{\prime }}\\ =\frac{\mu }{4\pi }{\displaystyle {\int }_{V^{\prime }}\left[\frac{1}{R}{\nabla }^{\prime }\cdot J\left(r^{\prime },t-\frac{R}{C}\right)-\frac{1}{R^2}{V}_{ez}\rho \right]\text{d}{V}^{\prime }}\\ =\frac{\mu }{4\pi }{\displaystyle {\int }_{V^{\prime }}\left[{\nabla }^{\prime }\cdot \left(\frac{1}{R}J\right)-{\nabla }^{\prime }\left(\frac{1}{R}\right)\cdot J-\frac{1}{R^2}{J}_z\right]\text{d}{V}^{\prime }}\\ =\frac{\mu }{4\pi }{\displaystyle {\int }_{V^{\prime }}\left[\nabla \left(\frac{1}{R}\right)-\frac{1}{R^2}{J}_z\right]\text{d}{V}^{\prime }}=\frac{\mu }{4\pi }{\displaystyle {\int }_{V^{\prime }}\left(-\frac{1}{R^2}{e}_R\cdot J-\frac{1}{R^2}{J}_z\right)\text{d}{V}^{\prime }}\\ =\frac{\mu }{4\pi }{\displaystyle {\int }_{V^{\prime }}\left(-\frac{1}{R^2}{J}_t-\frac{1}{R^2}{J}_z\right)\text{d}{V}^{\prime }}\approx \frac{\mu }{4\pi }{\displaystyle {\int }_{V^{\prime }}-\frac{1}{R^2}{J}_t\text{d}{V}^{\prime }}\\ =\frac{\mu }{4\pi }{\displaystyle {\int }_{V^{\prime }}-\frac{1}{R^2}{e}_R\cdot J\text{d}{V}^{\prime }}\end{array}$

When $\left|\frac{1}{R^2}{e}_R\cdot J\right|/\left|\frac{1}{R}\nabla \cdot J\right|=\frac{\left|V_{ez}{J}_t\right|}{\left|R{\omega }_{\begin{array}{c}P\end{array}}{J}_z\right|}\gg \frac{\left|V_{ez}\right|}{\left|R{\omega }_{\begin{array}{c}P\end{array}}\right|}\gg 1$ and $R\ll 3\times {10}^8\text{m}$ i.e., if it is ensured that $R\ll \frac{V_{ez}}{{\omega }_{\begin{array}{c}P\end{array}}}\approx \frac{C}{{\omega }_{\begin{array}{c}P\end{array}}}\ll C$, it can be considered that $\frac{R}{C}\ll 1\text{s}\approx 0\text{s}$, $t-\frac{R}{C}\approx t$ and

$\begin{array}{l}\frac{\mu }{4\pi }{\displaystyle {\int }_{V^{\prime }}\nabla \cdot \left[\frac{1}{R}J\left(r^{\prime },t-\frac{R}{C}\right)\right]\text{d}{V}^{\prime }}\approx \frac{\mu }{4\pi }{\displaystyle {\int }_{V^{\prime }}\nabla \cdot \left[\frac{1}{R}J\left(r^{\prime },t\right)\right]\text{d}{V}^{\prime }}\\ =\frac{\mu }{4\pi }{\displaystyle {\int }_{V^{\prime }}-\frac{1}{R^2}{e}_R\cdot J\text{d}{V}^{\prime }}+\frac{\mu }{4\pi }{\displaystyle {\int }_{V^{\prime }}\frac{1}{R}\nabla \cdot J\left(r^{\prime },t\right)\text{d}{V}^{\prime }}\approx \frac{\mu }{4\pi }{\displaystyle {\int }_{V^{\prime }}-\frac{1}{R^2}{e}_R\cdot J\text{d}{V}^{\prime }}\end{array}$, which gives

that

$A\left(R,t\right)=\frac{\mu }{4\pi }{\displaystyle {\int }_{V^{\prime }}\frac{1}{R}J\left(r^{\prime },t\right)\text{d}{V}^{\prime }}$, (3)

$\phi \left(R,t\right)=\frac{1}{4\pi \epsilon }{\displaystyle {\int }_{V^{\prime }}\frac{\rho \left(r^{\prime },t\right)}{R}\text{d}{V}^{\prime }}$. (4)

It is obvious that Equation (3) is not derived from Biot-Savart’s law [22] , and can be applied to the time-varying and waving current source $J\left(r^{\prime },t\right)$. Applying the second mean value theorem of integration [23] to Equation (3) and Equation (4) gives

$A\left(R,t\right)=\frac{\mu }{4\pi }\bar{J}\left(r^{\prime },t\right)F\left(R\right)$, (5)

$\phi \left(R,t\right)=\frac{1}{4\pi \epsilon }\bar{\rho }\left(r^{\prime },t\right)F\left(R\right)$, (6)

where $\bar{J}$and $\bar{\rho }$ are the space average values of $J$ and $\rho$ in $V^{\prime }$ while

$F\left(R\right)$ is defined as $F\left(R\right)={\displaystyle {\int }_{V^{\prime }}\frac{1}{R}\text{d}{V}^{\prime }}$.

Taking the limits of both sides of Equation (5) and Equation (6) gives

$\underset{V^{\prime }\to 0}{\lim }\phi \left(R,t\right)=\phi \left(R,t\right)=\underset{V^{\prime }\to 0}{\text{lim}}\frac{1}{4\pi \epsilon }\bar{\rho }\left(r^{\prime },t\right)F\left(R\right)=\frac{1}{4\pi \epsilon }\rho \left(r^{\prime },t\right)\underset{V^{\prime }\to 0}{\text{lim}}F\left(R\right)$, (7)

$\underset{V^{\prime }\to 0}{\lim }A\left(R,t\right)=A\left(R,t\right)=\underset{V^{\prime }\to 0}{\text{lim}}\frac{\mu }{4\pi }\bar{J}\left(r^{\prime },t\right)F\left(R\right)=\frac{\mu }{4\pi }J\left(r^{\prime },t\right)\underset{V^{\prime }\to 0}{\text{lim}}F\left(R\right)$. (8)

Since $\underset{V^{\prime }\to 0}{\text{lim}}F\left(R\right)$ is independent of the $\nabla$ operator and time $t$, and the source charges $\rho \left(r^{\prime },t\right)$ are continuously distributed and stationary, it follows

that ${\nabla }^2\phi =\frac{1}{4\pi \epsilon }{\nabla }^2\rho \underset{V^{\prime }\to 0}{\text{lim}}F\left(R\right)$, $\frac{{\partial }^2\phi }{\partial t^2}=\frac{1}{4\pi \epsilon }\frac{{\partial }^2\rho }{\partial t^2}\underset{V^{\prime }\to 0}{\text{lim}}F\left(R\right)$. Substituting the

above two equations into Equation (1), it is obtained that

${\nabla }^2\rho -\mu \epsilon \frac{{\partial }^2\rho }{\partial t^2}=-\frac{\rho }{\epsilon }4\pi \epsilon \frac{１}{\underset{V^{\prime }\to 0}{\text{lim}}F\left(R\right)}$. While the ρ in it is canceled, a new operator

$\square$ is defined as

$\square ={\nabla }^2-\mu \epsilon \frac{{\partial }^2}{\partial t^2}=-4\pi \frac{１}{\underset{V^{\prime }\to 0}{\text{lim}}F\left(R\right)}$. (9)

And then applying the new operator $\square$ to both sides of Equation (8) yields

${\nabla }^{2}A-\mu \epsilon \frac{{\partial }^{2}A}{\partial t^2}=-\mu J\left(r^{\prime },t\right)$. (10)

Equation (2) and Equation (12) are well known as d’ Alembert’s equations of the electromagnetic field [21] in the conductive medium.

According to Equation (28), the wave speed of $B_z$ in a vacuum is equal to the speed of light C. According to Equation (23), the frequency ${\omega }_{\begin{array}{c}P\end{array}}$ of $B_z$ can be

expressed as $\left|{\omega }_{\begin{array}{c}P\end{array}}\right|=\left|-\text{j}\frac{1}{{\tau }_2}\right|=\frac{{\omega }^2{\epsilon }_0}{{\sigma }_0}={10}^{29}\text{Hz}$, where ${\epsilon }_0~{10}^{-11}\text{F}/\text{m}$, is the vacuum

dielectric constant [21] , ${\sigma }_0~{10}^{-14}\text{S}/\text{m}$ is the atmospheric electrical conductivity [26] , which is approximately the vacuum electrical conductivity, $\omega$ taken as 10¹³ Hz (the frequency of light waves) is the source light waves frequency. The

typical wavelength of $B_z$ is calculated as ${\lambda }_z=\frac{C}{{\omega }_P}~{10}^{-21}\text{m}$. $B_{z0}$ equaling $B_{t0}\underset{V^{\prime }\to 0}{\lim }\left|\frac{\nabla \times {\stackrel{\rightharpoonup }{A}}_{\theta }}{\nabla \times {\stackrel{\rightharpoonup }{A}}_t}\right|=B_{t0}\underset{V^{\prime }\to 0}{\lim }\left|\frac{{\displaystyle {\int }_{V^{\prime }}-r^{-2}\left(e_r\times J_{\theta }\right)\text{d}{V}^{\prime }}}{{\displaystyle {\int }_{V^{\prime }}-r^{-2}\left(e_r\times J_t\right)\text{d}{V}^{\prime }}}\right|=B_{t0}\underset{V^{\prime }\to 0}{\lim }\left|\frac{{\bar{J}}_{\theta }}{{\bar{J}}_t}\right|=B_{t0}\frac{J_{\theta }}{J_t}$ is far

smaller than $B_{t0}$, which imply that it is possible that conventional equipment for detecting TEM waves may fail to detect LEM waves.

The ratio of the energy density of LEM wave to TEM waves can be expressed as

$E_P/E_T=\hslash {\omega }_{\begin{array}{c}P\end{array}}/\hslash \omega ={10}^{29}\text{Hz}/{10}^{13}\text{Hz}={10}^{16}$, (36)

which means that apart from $E_T=\hslash \omega$, $B_z$ must acquire huge additionally free energy from zero-point vacuum field [5] [27] . The mechanism for the free energy gained by $B_z$ is likely associated with its induced vortex-magnetic field $E_{\theta }$, which has the capability to extract zero-point vacuum energy through the entangled interaction between its produced torsion field and zero-point vacuum energy field. As per Equation (32), the total energy density of $B_z$ can be represented as $E_P=-\text{d}W/\text{d}t+E_T$, where $-\text{d}W/\text{d}t$ denotes work performed by the zero-point vacuum energy field on $B_z$. Because $E_P\gg E_T$, it follows that $E_P\approx -\text{d}W/\text{d}t$; thus, indicating that most of $B_z$’s energy originates from its interaction with “zero-point vacuum energy field”.

When two source light waves are superimposed in a vacuum, the resultant electromagnetic fields seemly both vanish. But the resultant magnetic vector potential

$A_{\theta }$ still persist and produce $E_{\theta }=-\partial A_{\theta }/\partial t=\sqrt{2}C{B}_0\exp \left[\text{j}\left(\frac{{\omega }_P}{C}z–{\omega }_{P}t\right)\right]$ and $B_z=\nabla \times A_{\theta }=\sqrt{2}{B}_0\exp \left[\text{j}\left(\frac{{\omega }_P}{C}z-{\omega }_{P}t\right)\right]$ called as “vacuum artificial free energy

field” with huge energy ( $\hslash {\omega }_{\begin{array}{c}P\end{array}}\gg \hslash \omega$).

This mechanism for LEM waves absorbing free energy is validated by the experiments of Meyl [10] [11] who found that LEM wave exhibits distinctive energetic properties and could increase ATP (adenosine triphosphate) content 40% in Epomoea chrysanthemum mitochondria, thereby extending the flowering period by nearly 10%. We can utilize this way to generate substantially free energy that might replace petrochemical and electrical energy and be the main energy in the future. Moreover, this free energy can be transmitted to thousands of households through a wireless network via the wireless energy transmission mode of LEM waves. This is an enticing new energy application scenario.

Owing to the skin effect, TEM waves are unable to penetrate conductive medium and can only propagate within a thin layer $\delta$ on its surface [28] . $B_z$ with a frequency 10²⁹ Hz possessess “high-frequency transparency” [29] [30] and can freely traverse through metallic media and human body (weak electrolyte). Particularly, it can non-destructively penetrate the blood barrier of the human brain, convey the drug information it carries to the brain tissue, and treat related diseases. What’s more, neutrinos are assumed as the propagating particles of $B_z$. The extremely small cross-section of 10⁻⁴³ cm² for neutrino to interact with atomic nuclei results in a minuscule probability of its capture by an atomic nucleus within a square centimeter area, thus demonstrating a harmlessly penetrative capability of $B_z$ to organism cells [31] .

The hydrogen bond energy between water molecule-groups is relatively weak and susceptible to be disrupted. Moreover, the vortex movement of LEM wave is also coherent with torsion field to extract the zero point energy, producing high-energy rays, neutrons and high-energy particles in its vortex center, accompanied by highly directional cold nuclear fusion [32] - [34] . On the one hand, the huge energy released by this process first breaks the fragile hydrogen bonds between water molecules (bond energy 23 kJ/mol, 0.24 eV) and produces smaller molecule groups of water. At the same time, it enhances the HO-H covalent bond energy of water molecules (bond energy 492 kJ/mol, 5.10 eV) [35] . On the other hand, it makes the bound electrons of water molecules be extremely active and become free electrons [35] . The increase of free electron concentration will weaken the ionization of water molecules, and finally form an ionization equilibrium state, which can increase the conductivity of water, and finally form high conductivity, weak alkaline-small molecule group-high energy state water. In this way, the energy of LEM waves is capable of being absorbed and stably stored by water for a very long time [36] .

Lee Sichen has demonstrated through nuclear magnetic resonance experiments that the size of water molecular groups can be reduced after absorbing LEM waves [37] . The results of the experiments show that after the LEM waves generated by LEM wave generator was directly irradiated to water for 100 minutes, there exhibited a significant increase in the half height and width change amplitude of the nuclear magnetic resonance signal of oxygen isotope O17 in water, reaching up to +10%. This LEM wave’s energy can remain stable in water for approximately one years [36] . Smaller water molecular groups are more adept at transporting oxygen and nutrient molecules into human cells, thereby enhancing blood oxygen concentration and nutritional absorption capacity. This is why Guopu from the Jin Dynasty of China called that “Qi” (LEM waves) “stops when it meets water and disperses when it encounters wind” [37] - [39] .

The particle model of magnetic P-wave $B_z$ in Figure 4 representing a high energy state consists of two neutrinos with same rotation and spin directions [18] [36] [40] . Assuming the neutrinos’ revolution velocity $V_R$ and radius $R$ are equal to light speed C and the wavelength ${\lambda }_z$ of $B_z$ in Figure 5 , the centripetal acceleration $a_R$ of the neutrinos’ revolution around the center core of $B_z$ is shown as

$a_R=-\frac{C^2}{{\lambda }_z}{e}_R=-\frac{C}{{\tau }_2}{e}_R$. (37)

When $R\ll C/{\omega }_P=3\times {10}^{-21}\text{m}$, and the source point region $V^{\prime }$ approaches zero that results in that $R=r$. In the field point coordinate system $P\left(r,\theta ,z\right)$

shown in Figure 1 , supposing the field point space $V$ expands along the positive r-axis at a velocity $V_r$, we can gain

$\begin{array}{c}{B}_z=\nabla \times A_{\theta }=\frac{{\mu }_0}{4\pi }{\displaystyle {\int }_{V^{\prime }}\nabla \times \left(\frac{1}{r}{J}_{\theta }\right)\text{d}{V}^{\prime }}\\ =\frac{{\mu }_0}{4\pi }{\displaystyle {\int }_{V^{\prime }}\left(-r^{-2}{J}_{\theta }{e}_z\right)\text{d}{V}^{\prime }}+jk\frac{{\mu }_0}{4\pi }{\displaystyle {\int }_{V^{\prime }}\left(-r^{-1}{J}_{\theta }{e}_r\right)\text{d}{V}^{\prime }}\\ \approx \frac{{\mu }_0}{4\pi }{\displaystyle {\int }_{V^{\prime }}\left(-r^{-2}{J}_{\theta }{e}_z\right)\text{d}{V}^{\prime }},\end{array}$ (38)

$E_{\theta }=-\nabla {\phi }_{\theta }-\frac{\partial A_{\theta }}{\partial t}=\frac{{\mu }_0}{4\pi }{\displaystyle {\int }_{V^{\prime }}\left(-r^{-2}{V}_r{J}_{\theta }-r^{-1}{V}_{r}j{\omega }_P{J}_{\theta }\right)\text{d}{V}^{\prime }}.$ (39)

From Equation (21), we are able to obtain

$\nabla \times B_z=\frac{1}{C^2}\frac{\partial E_{\theta }}{\partial t}=\frac{1}{C^2}{E}_{\theta }/{\tau }_2=\frac{{\mu }_0}{4\pi }{\displaystyle {\int }_{V^{\prime }}{r}^{-3}{J}_{\theta }\text{d}{V}^{\prime }}.$ (40)

From Equation (39), it is given that

$\begin{array}{c}\frac{\partial E_{\theta }}{\partial t}=\frac{\partial }{\partial t}\left(-\frac{\partial A_{\theta }}{\partial t}\right)\\ =a_r\frac{{\mu }_0}{4\pi }{\displaystyle {\int }_{V^{\prime }}-r^{-2}{J}_{\theta }\text{d}{V}^{\prime }}+V_r^2\frac{{\mu }_0}{4\pi }{\displaystyle {\int }_{V^{\prime }}2r^{-3}{J}_{\theta }\text{d}{V}^{\prime }}\\ +V_{r}j{\omega }_P\frac{{\mu }_0}{4\pi }{\displaystyle {\int }_{V^{\prime }}2r^{-2}{J}_{\theta }\text{d}{V}^{\prime }}-{\omega }_P^2{A}_{\theta }\\ =B_z\times a_r+V_r^2\nabla \times B_z+\left(2\frac{V_r}{C}-1\right)A_{\theta }/{\tau }_2^2\\ =B_z\times a_r+\frac{V_r^2}{C^2}{E}_{\theta }/{\tau }_2+\left(2\frac{V_r}{C}-1\right)A_{\theta }/{\tau }_2^2=E_{\theta }/{\tau }_2.\end{array}$ (41)

In Equation (41), the first term presumably represents the portion where the vacuum LEM wave $E_{\theta }$ generates the acceleration $a_r$ (gravitational field), the second term might be the part that gives rise to the weak nuclear force field, and the third term undoubtedly is the component that generates the electromagnetic field. Supposing that the portions of the vacuum LEM waves $B_z$ and $E_{\theta }$ that generates the acceleration $a_r$ still compliy with the wave equations of the LEM waves shown as Equation (24) and Equation (25), i.e.,

$\frac{1}{{\mu }_0{\epsilon }_0}{B}_z=C{e}_r\times E_{\theta }$ and $E_{\theta }=-C{e}_r\times B_z$, (42)

thus, we can obtain

$a_r=E_{\theta }\times \frac{B_z}{B_z^2{\tau }_2}=-\frac{C}{{\tau }_2}{e}_r=-\frac{C^2}{{\lambda }_z}{e}_r=a_R,$ (43)

which means that the acceleration produced by $E_{\theta }$ and $B_z$ is precisely the acceleration generated by the centripetal force field of the neutrino constituting the LEM wave and revolving around its center core in Figure 5 . To sum up, it should be inferred that vacuum LEM waves are capable of generating gravitational fields. This conclusion still requires a substantial amount of experimentation for validation. If the aforementioned theory holds true, humanity can utilize LEM waves to achieve anti-gravity fields. When the mass of an object becomes zero, in accordance with special relativity [41] , the object can reach the speed of light. Furthermore, based on the relativistic length contraction effect and time dilation effect [42] [43] , the space at the speed of light can be reduced to zero, and time comes to a halt, meaning that humanity can achieve interstellar teleportation and explore any location in the universe within a finite life cycle. Based on Equation (41) and Equation (43), it is obtained that

$\frac{1}{{\mu }_0{\epsilon }_0}\frac{\partial B_z}{\partial t}=a_r\times E_{\theta },$ (44)

$\frac{\partial E_{\theta }}{\partial t}=-a_r\times B_z,$ (45)

which are the fundamental formulas called as the Maxwell’s equations of the gravitational field generated by vacuum LEM waves. Analyzing (41) on the Earth’s

surface, we can gain $a_r=\frac{V_0^2}{r_D}=\frac{7.9\times 7.9\times {10}^6}{6.371\times {10}^6}=9.8\text{m}/{\text{s}}^2$ while $V_r$ presenting

space expanding speed of our vacuum is equal to $r_{D}H=6.371\times {10}^6\times 67.8\times$ ${10}^3/\left(3.08\times {10}^{22}\right)=1.4\times {10}^{-11}\text{m}/\text{s}\ll C$. Here, $V_0$ is the first cosmic velocity, $r_D$ is the radius of the Earth and $H$ is Hubble constant [44] , which gives that

$\frac{V_r}{C}~{10}^{-19}\ll 1$, $\left|B_z\times a_r\right|=\left(9.8/3\right)\times {10}^{-8}{\omega }_P{A}_{\theta }$, $\left|\frac{V_r^2}{C^2}{E}_{\theta }/{\tau }_2\right|={10}^{10}{\omega }_P{A}_{\theta }$ and

$\left|\left(2\frac{V_r}{C}-1\right){\omega }_P^2{A}_{\theta }\right|={10}^{29}{\omega }_P{A}_{\theta }$. Evidently, at this circumstance, the electromagnetic force: the weak nuclear force: the gravitational force =10³⁷:10¹⁸:1. When $V_r=C$, $\left|B_z\times a_r\right|=\left|E_{\theta }/{\tau }_2\right|=\left|\left(2\frac{V_r}{C}-1\right){\omega }_P^2{A}_{\theta }\right|={10}^{29}{\omega }_P{A}_{\theta }$, that is, when an object moves at

the speed of light, its electromagnetic field = weak nuclear force field = the gravitational field. However, since the object’s mass is zero at this circumstance, the gravitational force acting on the object is zero while its electromagnetic force remains unchanged, but the weak nuclear force is increased to equal the electromagnetic force. From the viewpoint of LEM waves, the weak nuclear field, gravitational field and electromagnetic field seem to be unified.