Metrological analysis showed that, under the law of mass-energy equivalence (LME) [1] , physical space must be finite and time dilation shall occur to unit of time defined on atomic clock (atomic time, AT) in motion therein [2] . It was also shown that, under the law of Planck on photon energy [3] , in conjunction with the law of energy conservation, Planck constant is a state invariant (SIT) in centripetal force field and time dilation shall occur in such field regardless of choice/definition of energy/time [4] . Static electric field (SEF) is a centripetal force field, therefore, presence of SEF shall alter unit of AT therein. Therefore, duration of unit of AT in SEF is not SIT but state variant of the field, i.e., function of location of clock defining particle at rest in the field. By definition, speed of light in vacuo (SLV) is inversely proportional to duration of unit of AT [5] . Metrological analysis showed that length of unit of length in space is SIT even under field environment [2] [4] . Therefore, SLV is not and cannot be SIT in SEF but varies in the field. However, constancy, hence state invariance, of SLV was assumed a prior in Schrödinger equation [6] and general category of wave mechanic equations (WME). It is therefore necessary to revise such in order to comply with LME and state-dependent SLV.

The original law of Coulomb on electrostatics (LCE) [7] for charge particles at rest in space is expressed as, under infinite space approximation (ISA),

$F_{k,j}=\frac{q_j{q}_k}{4\text{π}{\epsilon }_0}{f}_{k,j},f_{k,j}\equiv \frac{r_k-r_j}{{\Vert r_k-r_j\Vert }^3}$.(1)

$F_{k,j}$: Electrostatic force experienced by charge particle $k$ at rest in field of charge particle $j$ at rest. ${\epsilon }_0$: Electric constant. $q$: Rest charge of particle. $r$: Location vector of charge particle.

Retaining symmetry of the original law on the ground of interaction, LCE is rewritten as

$F_{k,j,x}=\frac{q_{j,x}{q}_{k,x}}{\sqrt{4\text{π}{\epsilon }_{0,j,x}}\sqrt{4\text{π}{\epsilon }_{0,k,x}}}{f}_{k,j},x=f,0,e$.(2)

$x$: State indicator, perspective at location of associated entity at rest in SEF.

By definition [8] ,

$\begin{array}{l}\frac{1}{4\text{π}{\epsilon }_0}=\frac{\alpha \hslash c}{e^2}\to \frac{1}{4\text{π}{\epsilon }_{0,x}}=\frac{{\alpha }_x{\hslash }_x{c}_x}{e_x^2}\to \\ F_{k,j,x}=\sqrt{{\alpha }_{j,x}{\hslash }_{j,x}{c}_{j,x}}\sqrt{{\alpha }_{k,x}{\hslash }_{k,x}{c}_{k,x}}\frac{q_{j,x}}{e_{j,x}}\frac{q_{k,x}}{e_{k,x}}{f}_{k,j}.\end{array}$(3)

$a$: Fine structure constant. $\hslash$: Planck constant, $\hslash \equiv h/2\text{π}$. $c$: SLV. $e$: Elementary charge.

By definition of charge,

$q\equiv n_{q}e\to q_x\equiv n_{q,x}{e}_x$.(4)

$n_q$: Numeral aspect of charge.

By the rule of metrology [4] ,

$n_{q,x}=\frac{q_x}{e_x}\subset \text{SIT}\to n_{q,x}=n_q\to q_x=n_q{e}_x$.(5)

That is, numeral aspect of charge of charge particle, referred to as reduced charge, is SIT if charge of an entity and unit of the charge are in identical state. Therefore, counted in reduced form, amount of charges of any charge system is SIT. Therefore, electric neutrality of electric neutral system is SIT, i.e., such system shall remain electrically neutral regardless of state of the system if charges of the system are counted in reduced form. Accordingly,

$F_{k,j,x}=c_i\sqrt{{\alpha }_{j,x}{\hslash }_{j,x}}\sqrt{{\alpha }_{k,x}{\hslash }_{k,x}}{\beta }_{e,j,x}^{1/2}{n}_{q,j}{\beta }_{e,k,x}^{1/2}{n}_{q,k}{f}_{k,j},{\beta }_{e,y,x}\equiv \frac{c_{y,x}}{c_{y,i}},y\in \left\{j,k\right\}$.(6)

${\beta }_e$: Referred to as Coulomb Factor. $c_i$: SLV as measured/defined in Rest State (RS) [2] . $c_{j,x}$: SLV as measured by particle $j$ at rest in field of others.

Since Planck constant is SIT in SEF,

${\hslash }_x={\hslash }_i\equiv \hslash \to F_{k,j,x}=\hslash c_{\text{i}}\sqrt{{\alpha }_{j,x}}\sqrt{{\alpha }_{k,x}}{\beta }_{e,j,x}^{1/2}{n}_{q,j}{\beta }_{e,k,x}^{1/2}{n}_{q,k}{f}_{k,j}$.(7)

Assumption 3: static fine structure constant is SIT,

${\alpha }_x={\alpha }_i\equiv \alpha \to F_{k,j,x}=\alpha \hslash c_i{\beta }_{e,j,x}^{1/2}{n}_{q,j}{\beta }_{e,k,x}^{1/2}{n}_{q,k}{f}_{k,j}$. (8)

That is, it is assumed that outcome of measurement of static fine structure constant shall not be affected by where in a field setup for such measurement is located regardless of strength of the field as long as the field is static in Rest Frame (RF) [2] and setup for the measurement is at rest in the field. Therefore, under Assumption 3, static fine structure constant is identical/one and same/indifferent whether measurement of the entity is conducted in far field (field strength approaching zero) or vicinity of electric Schwarzschild Sphere (ESS) as long as field and setup for measurement are both at rest with respect to RF.

Denote

$q^{\ast }\equiv {\beta }_{e,q}^{1/2}{n}_q\to F_{k,j,x}=\alpha \hslash c_i{q}_{j,x}^{\ast }{q}_{k,x}^{\ast }{f}_{k,j}$.(9)

$q^{\ast }$: Rest charge parameter of charge particle in reduced form at rest in field of others.

This is the LCE in infinite space under Assumption 3.

Consider a charge particle in RS in S³ (three dimensional finite space) with its location therein assigned as origin of RF. Then, exists SEF in association with the charge at rest in S³ with center of the field at the origin. Such field can be probed by test charge (effect of probe charge to field being probed is ignored/disregarded by definition of test charge). Thus, from Equation (2),

$F_{p,x}=\frac{q_{f,x}{q}_{p,x}{f}_s}{\sqrt{4\text{π}{\epsilon }_{0,f,x}}\sqrt{4\text{π}{\epsilon }_{0,p,x}}}{e}_{\phi }$.(10)

$F_{p,x}$: Electrostatic force probed by probe at rest in field. $q_{f,x}$: Rest charge of particle causing the field. ${\epsilon }_{0,f,x}$: Electric constant as measured at rest at location of the particle causing the field. $q_{p,x}$: Rest charge of probe at rest in the field. ${\epsilon }_{0,p,x}$: Electric constant as measured at rest at location of probe in the field. s: Internal distance between charge particle and probe. $f_s$: Function of $s$. $e_{\phi }$: Unit vector of $s$ at location of probe in field in S³. $x$: State indicator.

By definition of electric field,

$E\equiv \frac{F_p}{q_p}\to E_{f,x}\equiv \frac{F_{p,x}}{q_{p,x}}=\frac{q_{f,x}}{\sqrt{4\text{π}{\epsilon }_{0,f,x}}\sqrt{4\text{π}{\epsilon }_{0,p,x}}}{f}_s{e}_{\phi }$.(11)

$E$: Electric field. $q_p$: Probe charge. $E_f$: SEF associated with charge causing the field.

In classical field theories, field of an object is an intrinsic property of the object alone that is unaltered regardless of presence/absence and/or action/inaction of others. Accordingly, plurality of fields is vector additive in space. However, due to LME and state invariance of Planck constant in field, presence of charge in field of charge particle shall cause alteration of the electric constant associated with the particle hence rest charge parameter of same, therefore alteration of the field of the charge particle. Therefore, classical field is incompatible with LME, and SEF as expressed in Expression (11) also depends on others, hence plurality of the fields is not simple vector additive but interactive.

By flux conservation theorem of Gauss for SEF,

${\displaystyle \underset{S}{∯}{\epsilon }_{0}E\cdot dS}=q\to {\displaystyle \underset{\text{SPS}}{∯}\sqrt{{\epsilon }_{0,f,x}}\sqrt{{\epsilon }_{0,p,x}}{E}_{f,x}\cdot dS}=q_{f,x}$.(12)

$S$: Any nonlocal simultaneous enclosure of simple/single connectivity. $dS$: Areal element of the enclosure. $q$: Charge enclosed by the enclosure. SPS: Spherical probing sphere.

In spherical coordinates of S³ [2] ,

$d{S}_{\phi }=e_{\phi }{R}_e^2{\sin }^2\phi \sin \vartheta d\vartheta d\theta \to f_s=\frac{1}{R_e^2{\sin }^2\left[s/R_e\right]}$.(13)

$R_e$: External radius of S³. $d{S}_{\phi }$: Areal element of SPS.

Therefore, LCE in S³ is expressed as, for SEF centered at internal origin,

$F_e=\frac{\alpha \hslash c_i{q}_s^{\ast }{q}_e^{\ast }}{R_e^2{\sin }^2\left[s/R_e\right]}{e}_{\phi }$.(14)

$F_e$: Electrostatic force experienced by charge particle at rest in SEF. $q_s^{\ast }$: Rest charge parameter of charge causing the field. $q_e^{\ast }$: Rest charge parameter of charge particle at rest in the field. $s$: Internal distance between field causing charge and charge in the field.

In vicinity of the origin,

$s\ll R_e\to F_e\simeq \frac{\alpha \hslash c_i{q}_s^{*}{q}_e^{*}}{s^2}\hat{s}$.(15)

That is, in near field of origin (in comparison with external radius of finite space), LCE is of the familiar form, as an approximation.

Consider an electron at rest in field of a proton in RS at origin of RF. From Equation (15), electrostatic force experienced by the electron is, under ISA,

$\begin{array}{c}{\beta }_{e,p,x}=1\\ {\beta }_{e,e,x}\equiv {\beta }_e\end{array},\begin{array}{c}{n}_{q,p}=+1\\ n_{q,e}=-1\end{array}\to F_{e,x}=-\frac{\alpha \hslash c_i{\beta }_e^{1/2}}{r^2}\hat{r}$.(16)

${\beta }_{e,p,x}$: Coulomb Factor at rest at location of proton in field of others. ${\beta }_{e,e,x}$: Coulomb Factor at rest at location of electron in SEF of proton. $n_{q,p}$: Reduced charge of proton. $n_{q,e}$: Reduced charge of electron. $F_{e,x}$: Electrostatic force experienced by electron at rest in SEF of proton centered at origin. $r$: Distance between electron and proton. $\hat{r}$: Unit vector of $r$ at location of electron in SEF of proton.

Denote

$r_e\equiv \frac{\alpha \hslash c_i}{E_{e,i}}\text{,}\rho \equiv \frac{r}{r_e}\to F_{e,x}=-\frac{E_{e,i}}{r_e}\frac{{\beta }_e^{1/2}}{{\rho }^2}\hat{\rho }$.(17)

$r_e$: Characteristic length of field (CLF) of elementary charge in RS. $E_{e,i}$: Self energy of electron in RS. $\rho$: Distance between electron and proton, in reduced unit.

The CLF above is also known as classical electron radius [9] .

Consider relocating the electron in SEF of proton in infinitesimal displacement/velocity. Work done to the electron is

$d{w}_{e,x}=-F_{e,x}\cdot dr\to \frac{d{w}_{e,x}}{d\rho }=E_{e,i}\frac{{\beta }_e^{1/2}}{{\rho }^2}$.(18)

$d{w}_{e,x}$: Work done for infinitesimal displacement of electron in infinitesimal velocity in SEF of proton.

By the law of energy conservation, work done to the electron in such manner is gained by the electron as self energy increment,

$d{E}_{e,x}=d{w}_{e,x}\to \frac{d{E}_{e,x}}{d\rho }=E_{e,i}\frac{{\beta }_e^{1/2}}{{\rho }^2}$.(19)

$E_{e,x}$: Self energy of electron at rest in SEF of proton.

It has been shown that static Planck constant in centripetal force field is field invariant [4] and

$E_{s,0,f}={\beta }_f{E}_{s,i},{\beta }_f\equiv \frac{c_{s,0,f}}{c_i}$.(20)

$E_{s,0,f}$: Self energy of a particle at rest in static centripetal force field. $c_{s,0,f}$: SLV as measured at rest at location of the particle in the field. $E_{s,i}$: Self energy of the particle in RS.

SEF is centripetal force field, therefore, with Equation (19),

$\frac{d{\beta }_e}{d\rho }=\frac{{\beta }_e^{1/2}}{{\rho }^2}\to {\beta }_e^{1/2}=1-\frac{1}{2\rho }\to \begin{array}{c}{\rho }_{\text{ESS}}\equiv {\rho |}_{{\beta }_e=0}=\frac{1}{2}\\ {\rho }_{\text{ESS}}\le \rho <\infty \end{array}$.(21)

${\rho }_{\text{ESS}}$: Radius of ESS in reduced unit.

Figure 1 plots out self energy of electron in SEF of proton, in reduced unit. Note that left bound of the domain is at radius of ESS. Left side of the left bound is forbidden zone since, by definition, ${\beta }_e^{1/2}$ must be a real entity.

For S³,

${\beta }_e^{1/2}=1-\frac{\cot \phi }{2R_e},\phi \equiv \frac{\rho }{R_e}\to {\phi }_{\text{ESS}}=\arctan \frac{1}{2R_e}<\phi \le \frac{\text{π}}{2}$.(22)

$R_e$: External radius of S³ in unit of CLF.

No point in finite space is free of field of any kind. Therefore, reference state for the field is chosen at internal radius of the space, whereat, strength of the field centered at internal origin is minimal hence an approximation to true RS.

From Equation (17), under ISA,

$\begin{array}{c}{ℱ}_e\equiv \frac{r_e}{E_{e,i}}{F}_e\to {ℱ}_{e,x}=-\left(1-\frac{1}{2\rho }\right)\frac{\hat{\rho }}{{\rho }^2}\text{,}\frac{1}{2}\le \rho <\infty \to \\ {ℱ}_{e,x,\max }=\frac{16}{27},{\rho }_{{ℱ}_{e,x},\max }=\frac{3}{4},ℱ\equiv \Vert ℱ\Vert \end{array}$.(23)

${ℱ}_{e,x}$: Electrostatic force experienced by an electron at rest in SEF of proton, in reduced unit.

That is, strength of electrostatic attraction between electron and proton is not a monotonic function of the distance between the parties but shall reach maxima at 3/4 of the distance unit and then approach zero towards ESS and become so at ESS. Therefore, if electron is sufficiently close to proton then electrostatic interaction between the parties shall become diminishing. On the other hand, reduced charges of the parties are SIT hence shall remain intact.

From Equation (20), under ISA,

$E_{e,s,0,e}=E_{e,x}={\beta }_e{E}_{e,i}={\left(1-\frac{1}{2\rho }\right)}^2{E}_{e,i}\to {E_{e,x}|}_{{\rho }_{\text{ESS}}\le \rho <\infty }<E_{e,i}$.(24)

$E_{e,s,0,e}$: Self energy of electron at rest in SEF of proton. $E_{e,x}$: Rest energy of electron at rest in SEF of proton.

That is, self energy of an electron shall approach zero towards ESS and become so at ESS. Therefore, if an electron is at rest in vicinity of ESS then self energy of the electron hence difference thereof shall be approaching zero. Further, at ESS, electron can neither emit nor absorb photon of any kind. For same reason, if an electron is in motion in close proximity of ESS then emission/absorption of photon shall not be an effective channel for the electron to lose/gain energy. Therefore, close proximity of ESS is zone of optical silence.

By definition of Coulomb Factor, Expression (20),

$c_{e,x}={\beta }_e{c}_i\to {c_{e,x}|}_{{\rho }_{\text{ESS}}\le \rho <\infty }<c_i$.(25)

That is, SLV in SEF shall approach zero towards ESS and become so at ESS. Therefore, if an electron is at rest in vicinity of ESS then SLV as measured by the electron shall be slower than that as measured away from ESS. Therefore, if an electron is in motion in vicinity of ESS then velocity of the electron shall be slower than that in region away from ESS even if it were approaching local SLV in the field. Further, if an electron were to land at ESS then magnitude of its touchdown velocity cannot exceed local SLV thereat, which is zero.

By definition of SLV [10] , with state invariance of unit of length and Planck constant in SEF,

${\mathbb{U}}_{t,x}=\frac{{\mathbb{U}}_{t,i}}{{\beta }_e}\to {{\mathbb{U}}_{t,x}|}_{{\rho }_{\text{ESS}}\le \rho <\infty }>{\mathbb{U}}_{t,i}$.(26)

${\mathbb{U}}_{t,x}$: Unit of self time of electron at rest in SEF of proton. ${\mathbb{U}}_{t,i}$: Unit of Rest Time [2] .

That is, time dilation of unit of self time of an electron in SEF shall occur regardless of choice/definition of the self time. If the electron is to approach ESS then the clock defining self time of the electron shall cease to work. On the other hand, lacking of particular set of recurring events, by itself, shall have no effect to any other event of the electron, i.e., such by itself shall not prevent other event of the electron from happening, since events happen or not happen regardless of status of assigned set of recurring events [2] .

By LME, from Equation (20),

$m_{e,s,0,e}=m_{e,x}=\frac{m_{e,i}}{{\beta }_e}\to {m_{e,s,0,e}|}_{{\rho }_{\text{ESS}}\le \rho <\infty }>m_{e,i}$.(27)

$m_{e,s,0,e}$: Self mass of electron at rest in SEF of proton. $m_{e,x}$: Rest mass of electron at rest in SEF of proton. $m_{e,i}$: Self/rest mass of electron in RS.

That is, rest mass of an electron at rest in SEF of proton shall approach infinity towards ESS and become so at ESS. Therefore, if an electron is at rest in vicinity of ESS then rest mass of the electron as measured by itself and/or at rest at location of the electron in the field shall be heavier than otherwise. Therefore, according to the law of gravitation of Newton, gravitation interaction between the parties shall become non-negligible. Therefore, the electron shall be gravitationally attracted towards or repelled from the proton pending on type of rest mass of electron with respect to that of proton.

In particle model ignoring spin of electron, self field interaction and magnetic field induced by motion of electron, other types of forces, etc., equation of motion of an electron in SEF of proton is, under Assumption 3 and ISA, from Equation (17),

$\frac{d{P}_{e,f,u,e}}{d{t}_{e,s,u,e}}=F_{e,f,0,e}=-\frac{m_{e,i}{c}_i^2}{r_e}\frac{{\beta }_e^{1/2}}{{\rho }^2}\hat{\rho },{\beta }_e^{1/2}=1-\frac{1}{2\rho },\frac{1}{2}\le \rho <\infty$.(28)

$P_{e,f,u,e}$: Field momentum of an electron in motion in SEF of proton at origin, as perceived at rest at location of the electron in the field. $t_{e,s,u,e}$: Self time of the electron in motion in the field [4] .

According to the second law of Newtonian mechanics, alteration of momentum of a particle is caused by and equal to net force experienced by the particle. However, momentum is an attribute of the object it is in association with that is perspective dependent, and there are two but only two perspectives in particle-field system, i.e., perspective of particle and that of field. From perspective of a particle, momentum of the particle is and is always zero regardless of the force it is experienced since a particle is and is always at rest with respect to itself, regardless. Therefore, the momentum referred to by the law has to be the momentum perceived from perspective of field, i.e., field momentum. There also have two but only two types of time available in particle-field system, i.e., self time of particle and local time in field. However, due to time dilation in centripetal force field, Equation (26), there can be no such thing, in general, as locally defined common time of the field. In other words, local time in such field is generally non differentiable. Therefore, the time differential referred to by the law has to be self time of particle.

By definition of field momentum, under LME,

$P_{e,f,u,e}\equiv m_{e,f,u,e}{v}_{e,f,u,e}=\frac{m_{e,f,0,e}}{{\beta }_u}{c}_{e,f,u,e}u=\frac{m_{e,f,0,e}}{{\beta }_u}{c}_{e,f,0,e}u=m_{e,i}{c}_i\frac{u}{{\beta }_u}$.(29)

$m_{e,f,u,e}$: Mass of an electron in motion in SEF of proton as perceived at rest at location of the electron in the field. $v_{e,f,u,e}$: Velocity of an electron in motion in SEF of proton as perceived at rest at location of the electron in the field. ${\beta }_u$: Lorentz Factor of electron in absolute motion in finite space/bound state [2] . $u$: Reduced velocity of electron in motion.

By definition of self time, under LME,

$d{t}_{e,s,u,e}=\frac{d{t}_i}{{\beta }_u{\beta }_e}\to {\beta }_u\frac{d}{d\tau }\frac{u}{{\beta }_u}=-\frac{1}{{\beta }_e^{1/2}}\frac{\hat{\rho }}{{\rho }^2},d\tau \equiv \frac{c_i}{r_e}d{t}_i$.(30)

$t_i$: Rest Time [2] [4] .

Thus, in reduced units, equation of motion of an electron in SEF of proton is

$a-f_{u}u+f_e\hat{\rho }=0,a\equiv \frac{du}{d\tau }\text{,}{f}_u\equiv \frac{d\ln {\beta }_u}{d\tau }\text{,}{\beta }_u\equiv \sqrt{1-u^2}\text{,}{f}_e\equiv \frac{1}{{\beta }_e^{1/2}{\rho }^2}$.(31)

Therefore, motion of charge particle in SEF is identical in form to that of mass particle in static gravitation field (SGF) [4] .

With cylindrical coordination system, Equation (31) is decomposed to

$\begin{array}{l}\begin{array}{c}0=\rho \psi +2\omega q-f_u\rho \omega \\ 0=\chi -\rho {\omega }^2-f_{u}q+f_e\end{array},\begin{array}{c}q\equiv {\rho }^{\prime }\\ \chi \equiv q^{\prime }\end{array},\begin{array}{c}\omega \equiv {\theta }^{\prime }\\ \psi \equiv {\omega }^{\prime }\end{array},\text{<msup> <mo> ′ </mo> </msup>}\equiv \frac{d}{d\tau }\to \\ \begin{array}{c}\psi =q\omega \left(f_e-\frac{2}{\rho }\right)\\ \chi =\rho {\omega }^2-\left(1-q^2\right)f_e\end{array}\to \begin{array}{c}d\ln \frac{{\beta }_e}{{\beta }_u}=0\\ d\ln \frac{{\rho }^2\omega }{{\beta }_u}=0\end{array}\to \begin{array}{c}\frac{{\beta }_e}{{\beta }_u}=ϵ\\ \frac{{\rho }^2\omega }{{\beta }_u}=\gamma \end{array}\end{array}$.(32)

$ϵ$: Integration constant, total energy of an electron in SEF of proton in reduced unit. γ: Integration constant, angular momentum of an electron in SEF of proton in reduced unit.

Thus,

$\begin{array}{c}{\omega }^2=\frac{{\gamma }^2{\beta }_e^2}{{ϵ}^2{\rho }^4},q^2=1-\frac{{\beta }_e^2}{{ϵ}^2}\left(\text{1+}\frac{{\gamma }^{\text{2}}}{{\rho }^2}\right),ϵ\ne 0\to \\ \underset{{\beta }_e\to 0}{\lim }{\beta }_u=0,\underset{{\beta }_u\to 0}{\lim }\rho \omega =0\to \underset{{\beta }_e\to 0}{\lim }{u}^2=1,\underset{{\beta }_u\to 0}{\lim }{q}^2=1\end{array}$.(33)

That is, if an electron were to rendezvous onto ESS of proton then, at landing site, transverse velocity of the electron shall approach zero and impact velocity approach local SLV, which is zero along surface normal at landing site, while total energy and angular momentum of the electron shall remain intact during such rendezvous.

For an electron in orbital motion outside ESS of SEF of proton, at apsis distances,

${q|}_{\rho =a,b}=0,\frac{1}{2}<b\le a\to \begin{array}{c}{ϵ}^2=\frac{\left(a^2-b^2\right){\beta }_{e,a}^2{\beta }_{e,b}^2}{a^2{\beta }_{e,b}^2-b^2{\beta }_{e,a}^2}\\ {\gamma }^2=\frac{a^2{b}^2\left({\beta }_{e,a}^2-{\beta }_{e,b}^2\right)}{a^2{\beta }_{e,b}^2-b^2{\beta }_{e,a}^2}\end{array}$.(34)

$a$ and $b$: Maximum and minimum distance between orbiting electron and SEF center in reduced unit.

Therefore, orbital energy of electron in SEF of proton is completely determined by orbital apsides, and circular orbit is of relatively lower energy. Further, exists minima in orbital energy of electron in circular motion in SEF of proton, referred to herein as ground state,

$\underset{b\to a}{\lim }{ϵ}_{\text{o}}=\frac{{\left(2a-1\right)}^{5/2}}{4a^2\sqrt{2a-3}}\to {a>\frac{3}{2}\text{,}{ϵ}_{\text{o},\min }={ϵ}_{\text{o}}|}_{a=3}=\frac{1}{4}{\left(\frac{5}{3}\right)}^{5/2},u_{\text{o,min}}={\left(\frac{2}{5}\right)}^{1/2}$.(35)

${ϵ}_{\text{o}}$: Total energy of an electron in circular orbit in SEF of proton, in reduced unit.

If an electron is orbiting in vicinity of ESS of proton closer than that of the ground state then orbital energy of the electron shall become higher for such maneuvering, which shall approach infinity at orbital radius ${\rho }_{\text{ESS}}+1$, i.e., one $r_e$ away from the ESS. Therefore, even under tremendous electrostatic attraction in vicinity of ESS, an electron in circular orbital motion cannot maneuver onto the ESS. On the other hand, even in the ground state, total energy of an electron is marginally lower than that in RS, i.e., $~0.987E_{e,i}$.

By definition of apsidal precession, from Equation (33),

$\begin{array}{l}{T}_{AP}=\frac{\text{π}{\left(a+b\right)}^{3/2}}{\sqrt{2}}\left(1+\frac{2}{a+b}+\cdots \right),\\ {\theta }_{AP}=\frac{5\text{π}}{4}\left(\frac{1}{a}+\frac{1}{b}\right)+\frac{95\text{π}}{64}{\left(\frac{1}{a}+\frac{1}{b}\right)}^2+\cdots .\end{array}$(36)

$T_{AP}$: Apsidal orbital period, duration from a moment of an apside configuration to next of same, in reduced Rest Time. ${\theta }_{AP}$: Apsidal precession of an electron orbiting in SEF of proton per apsidal orbital period.

That is, electron orbiting in SEF of proton shall have apsidal precession, which is always positive, i.e., beyond 2per apsidal orbital period.

In contrast to particle model, object in wave mechanics (WM) [11] is perceived as extended object (EO) or in association with such. Due to time dilation in centripetal force field, self time of an EO in such field is generally indefinable hence the concept meaningless. On the other hand, total energy of an electron in such field is invariant to location, time, state of motion of the electron, etc., ignoring spin, gravitation, interaction of self field and magnetic field associated with motion of the electron, and assuming stativity of the field. Such invariance of such attribute can be utilized to decipher energy state of such EO in such field by WM.

From Expression (33), under Assumption 3 and ISA,

${ϵ}_e\equiv \frac{E_{e,f,u,e}}{E_{e,i}}=\frac{{\beta }_e}{{\beta }_u}\text{,}{ϵ}_e\ne 0\to \frac{{ϵ}_e^2}{{\beta }_e^2}=\frac{1}{1-u^2},{\beta }_e\ne 0$.(37)

${ϵ}_e$: Total energy of an electron in motion in SEF as measured at rest in the field, in reduced unit. $E_{e,f,u,e}$: Total energy of an electron in motion in SEF as measured at rest in the field. $E_{e,i}$: Self energy of electron in RS.

From Expression (29),

$\begin{array}{c}{P}_{e,f,u,e}=m_{e,i}{c}_i\frac{u}{{\beta }_u}\to p_e\equiv \frac{P_{e,f,u,e}}{m_{e,i}{c}_i}=\frac{u}{{\beta }_u}\to \\ 1+p_e^2=\frac{1}{1-u^2}\to -p_e^2+\frac{{ϵ}_e^2}{{\beta }_e^2}=1.\end{array}$(38)

$P_{e,f,u,e}$: Field momentum of an electron in motion in SEF as perceived at rest in the field. $p_e$: Field momentum of an electron in motion in SEF as perceived at rest in the field, in reduced unit.

According to WM [6] , exists correspondence relationship between field momentum of an object and spacial gradient of wave function (WF) of same, expressed as

$P_e\sim i\hslash {\nabla }_r,i\equiv \sqrt{-1}\to p_e\sim \frac{1}{m_{e,i}{c}_i}\frac{i{\hslash }_i}{r_e}{\nabla }_{\rho }=\frac{i}{{\alpha }_e}{\nabla }_{\rho }$.(39)

${\nabla }_r$: Gradient operator in unit of length of space. ${\nabla }_{\rho }$: Gradient operator in unit of CLF. ${\alpha }_e$: Fine structure constant.

Thus, under Assumption 3,

$\begin{array}{c}\frac{{\Delta }_{\rho }}{{\alpha }_e^2}+\frac{{ϵ}_e^2}{{\beta }_e^2}=\hat{I},{\beta }_e^{1/2}=1-\frac{1}{2\rho },\frac{1}{2}<\rho <\infty \to \\ \left(\frac{{\Delta }_{\rho }}{{\alpha }_e^2}+\frac{{ϵ}_e^2}{{\beta }_e^2}-1\right)\Psi =0\text{,}\begin{array}{c}\underset{\rho \to 1/2}{\lim }\Psi =0\\ \underset{\rho \to \infty }{\lim }\Psi =0\end{array}.\end{array}$(40)

${\Delta }_{\rho }$: Laplace operator in unit of CLF. $\hat{I}$: Identity operator. Ψ: WF of EO associated with an electron in motion in SEF of proton.

This is the WME complying with LME and finite space/bound state for electron in SEF of proton under the approximations aforementioned. In such, attributes are expressed in reduced form, i.e., in units of the corresponding ones measured/defined in RS, e.g., energy in unit of $m_{e,i}{c}_i^2$, momentum in $m_{e,i}{c}_i$, length in $r_e$, etc. Note that it is the domain boundaries of the Coulomb Factor that determines the boundaries of ρ in the WME, even though the equation itself may appear not being bounded by such. The boundary conditions of Ψ are to ensure convergence of the WF in entire domain that is valid physically.

In far field,

$\rho \gg \frac{1}{2}\to \frac{1}{{\beta }_e^2}\simeq 1+\frac{2}{\rho }\to \left(\frac{{\Delta }_{\rho }}{2{\alpha }_e^2}+\frac{1}{\rho }\right)\Psi \simeq -{\epsilon }_e\Psi ,\begin{array}{l}{ϵ}_e\equiv 1-{\epsilon }_e\hfill \\ 0\le {\epsilon }_e\ll 1\hfill \end{array}$.(41)

Therefore, the classical equation is but an approximation of WME (40) under the far field condition, whiles form of the approximated WME may lead to over interpretation of the corresponding WF in near field.

Note also that any WF has or should have unit in association since WF is an attribute of the object it is in association with, and attribute must have unit by definition of attribute [2] . Assuming unit associated with WF is field invariant, then such unit can be eliminated from WME due to form of WME. In addition, if a WF is a solution of a WME then any WF in constant proportion to the WF is also a solution of the WME, due to form of WME. Therefore, standardization of WF is necessary, via, e.g., normalization.

Formal solution of WME (40) involves D-finite function, which is not transparent for analysis. To solve the WME for analysis, take the first order far field approximation while retaining the domain boundaries of $\rho$,

$\begin{array}{c}\rho \gg \frac{1}{2}\to \frac{1}{{\beta }_e^2}\simeq 1+\frac{2}{\rho -1/2}\equiv f_1,\frac{1}{{\rho }^2}\simeq \frac{1}{{\left(\rho -1/2\right)}^2}\to \\ \frac{1}{X}\frac{d^{2}X}{d{x}^2}-\frac{l\left(l+1\right)}{x^2}-\frac{1}{4}+\frac{\delta }{x}\simeq 0,\begin{array}{c}\Psi \equiv R{Y}_l^m\\ X\equiv \rho R\end{array},\\ x\equiv 2\nu \left(\rho -1/2\right),\nu \equiv {\alpha }_e\sqrt{1-{ϵ}_e^2},\delta \equiv \frac{{\alpha }_e{ϵ}_e^2}{\sqrt{1-{ϵ}_e^2}}\end{array}$.(42)

$Y_l^m$: Normalized spherical harmonic function.

Solution of the approximated WME is

$\Psi =\frac{1}{\rho }\left(c_1{M}_{\delta }^{l+1/2}\left[2\nu \left(\rho -\frac{1}{2}\right)\right]+c_2{W}_{\delta }^{l+1/2}\left[2\nu \left(\rho -\frac{1}{2}\right)\right]\right)Y_l^m$.(43)

$c_{1\&2}$: Linear combination coefficients. $M_a^b$: Whittaker M function. $W_a^b$: Whittaker W function.

The solution above is generally divergent at the domain boundaries pending on values of δ and $l$. However, solution of WME represents physical aspect of the object in association with and physical entity cannot have aspect, attribute, or property being/becoming divergent, even at domain boundary. Therefore, Expression (43) cannot be the WF representing the object. On the other hand, the expression is a general solution of Equation (42) hence the solution of the WME according to the uniqueness theorem of solution of differential equation [12] .

To obtain WF for WME (40), let

$\delta =n,n=1,2,\cdots \to {ϵ}_n^2=\frac{1}{2}\frac{n^2}{{\alpha }_e^2}\left(\sqrt{1+4\frac{{\alpha }_e^2}{n^2}}-1\right)=1-\frac{{\alpha }_e^2}{n^2}+2\frac{{\alpha }_e^4}{n^4}+\cdots$.(44)

Define

$\begin{array}{c}\begin{array}{c}{\Phi }_{n,l,m}\equiv R_n^l{Y}_l^m\\ n>l\ge \left|m\right|\ge 0\end{array},R_n^l\equiv \frac{{\alpha }_e{ϵ}_n{\rho }^{-1}}{n\sqrt{\left(n-l-1\right)!\left(n+l\right)!}}{W}_n^{l+\frac{1}{2}}\left[\frac{2{\alpha }_e^2{ϵ}_n^2}{n}\left(\rho -\frac{1}{2}\right)\right]\to \\ \underset{\rho \to 1/2}{\lim }{\Phi }_{n,l,m}=0,\underset{\rho \to \infty }{\lim }{\Phi }_{n,l,m}=0\end{array}$.(45)

Φ: Base function for WME (40).

All the base functionsabove are real, finite, continuous, differentiable and integrable in the domain meeting the boundary conditions that form a base function set that is complete and irreducible. Any real, finite, and continuous function in real domain with specified boundary conditions can be expressed as a linear combination of members of complete and irreducible set of real, finite, continuous, and differentiable functions in same domain under same boundary conditions, if the function is decomposable with respect to such set [13] . WF is real, finite, and continuous function meeting boundary conditions specified, by definition of WF. Therefore, WF of WME (40) can be expressed as

$\begin{array}{c}\Psi ={\displaystyle \underset{n,l,m}{\sum }{c}_{n,l,m}{\Phi }_{n,l,m}},{\displaystyle \underset{n,l,m}{\sum }{c}_{n,l,m}^2}\ne 0\to \\ {\displaystyle \underset{n,l,m}{\sum }{c}_{n,l,m}\left({\Delta }_{\rho }+\frac{{\alpha }_e^2{ϵ}_e^2}{{\beta }_e^2}-{\alpha }_e^2\right){\Phi }_{n,l,m}}=0\end{array}$.(46)

$c_{n,l,m}$: Linear combination coefficient, real numeral constant determined by function decomposition procedure.

Therefore,

${\displaystyle \underset{n,l,m}{\sum }{c}_{n,l,m}\left(\frac{{ϵ}_e^2}{{\beta }_e^2}-f_1{ϵ}_n^2-\frac{l\left(l+1\right)}{{\alpha }_e^2}\left(\frac{1}{{\rho }^2}-\frac{1}{{\left(\rho -1/2\right)}^2}\right)\right)|R_n^l{Y}_l^m\rangle }=0$.(47)

$|\rangle$: Dirac notation, right bracket [14] .

For $l\ne 0$,

${\displaystyle \underset{n,n^{\prime }>l>0}{\sum }{c}_{n,l}\langle R_{n^{\prime }}^l|\frac{{ϵ}_e^2}{{\beta }_e^2}-f_1{ϵ}_n^2-\frac{l\left(l+1\right)}{{\alpha }_e^2}\left(\frac{1}{{\rho }^2}-\frac{1}{{\left(\rho -1/2\right)}^2}\right)|R_n^l\rangle }=0$.(48)

$\langle |$: Dirac notation, left bracket [14] .

The equations are linearly independent with respect to each other and organized per the variable $l$.

In matrix form,

$\begin{array}{c}\left(\eta {ϵ}_e^2-\nu -l\left(l+1\right)\mu \right)c=0\\ \begin{array}{c}{\eta }_{j,k}\equiv \langle R_j^l|\frac{1}{{\beta }_e^2}|R_k^l\rangle \text{,}{\mu }_{j,k}\equiv \frac{1}{{\alpha }_e^2}\langle R_j^l|\frac{1}{{\rho }^2}-\frac{1}{{\left(\rho -1/2\right)}^2}|R_k^l\rangle \\ {\nu }_{j,k}\equiv \langle R_j^l|f_1|R_k^l\rangle {ϵ}_k^2,j,k\in \left\{n>l>0\right\}.\end{array}\end{array}$.(49)

$c$: Real coefficient vector to be determined. 0: Column vector of zeros.

Since the base function set is complete and irreducible, matrix $\eta$ is invertible. Therefore,

$\left(I{ϵ}_e^2-\lambda \right)c=0,\lambda \equiv {\eta }^{-1}\left(\nu +l\left(l+1\right)\mu \right)\to \left|I{ϵ}_e^2-\lambda \right|\left|c\right|=\left|0\right|=0$.(50)

I: Identity matrix. $\left|\right|$: For vector, square root of dot product of the vector; for square matrix, determinant of the matrix.

Since the linear coefficient vector $c$ is nonzero, determinant of the matrix above, which is but a polynomial of ${ϵ}_e^2$, has to be zero,

$\left|c\right|\ne 0\to \left|I{ϵ}_e^2-\lambda \right|=0\to {ϵ}_{e,n,l}^2={\lambda }_{n,l}\to {ϵ}_{e,n,l}=\sqrt{{\lambda }_{n,l}}$.(51)

${\lambda }_{n,l}$: The nth root of the polynomial, i.e., eigenvalue of matrix $\lambda$, which is $l$-dependent.

That is, total energy of an electron in SEF of proton is discrete by WM. In comparison with that in the classical case, the energy levels are $l$-dependent.

For $l=0$, element of matrix $\eta$ is divergent due to insufficient rank of zero point of the base functions at left boundary of $\rho$. Therefore, take far field approximation of ${\beta }_e^{-2}$ to the second order,

$\begin{array}{c}\rho \gg \frac{1}{2}\to \frac{1}{{\beta }_e^2}\simeq 1+\frac{2}{\rho -1/2}+\frac{3}{{\left(\rho -1/2\right)}^2}\equiv f_2\to \\ {\displaystyle \underset{n,n^{\prime }>0}{\sum }{c}_{n,0}\langle R_{n^{\prime }}^0|f_2{ϵ}_e^2-f_1{ϵ}_n^2|R_n^0\rangle }\simeq 0.\end{array}$(52)

In matrix form,

$\left(I{ϵ}_e^2-\lambda \right)c\simeq 0,\lambda \equiv {\eta }^{-1}\nu \text{,}\begin{array}{c}{\eta }_{j,k}\equiv \langle R_j^0|f_2|R_k^0\rangle \\ {\nu }_{j,k}\equiv \langle R_j^0|f_1|R_k^0\rangle {ϵ}_k^2\end{array}\to {ϵ}_{e,n,0}\simeq \sqrt{{\lambda }_{n,0}}$.(53)

Due to the far field approximation taken, the result above is not exact with respect to WME (40). On the other hand, retention of the boundary condition of ρ assures correctness of qualitative aspect of profiles of the WFs in near field, especially that of the ground state, which is distinctly different from that in the classical case. Figure 2 plots out some of the base functions for an electron in SEF of proton, wherein, abscissa of the plot is in unit of $r_e/{\alpha }_e^2$. It can be seen from the plot that profiles of the base functions are similar to that in the classical case except the one corresponding to the ground state, which is zero at left bound of the domain hence different from that in the classical case. In general, base functions are not eigenfunctions of WME and linear combinations of the base functions are, that is also different with respect to that in the classical case.

Regard hydrogen atom as a binary particle system (BPS) at rest in RF. From Equation (8), under ISA,

$\begin{array}{l}{F}_{e,x}=-\frac{E_{e,i}{\beta }_{e,x}^{1/2}{\beta }_{p,x}^{1/2}}{r_e}\frac{{\rho }_e-{\rho }_p}{{\Vert {\rho }_e-{\rho }_p\Vert }^3},F_{p,x}=-\frac{E_{e,i}{\beta }_{p,x}^{1/2}{\beta }_{e,x}^{1/2}}{r_e}\frac{{\rho }_p-{\rho }_e}{{\Vert {\rho }_p-{\rho }_e\Vert }^3}\to \\ \begin{array}{c}\frac{d{E}_{e,x}}{E_{e,i}}=d{\beta }_{e,x}=-\frac{r_e{F}_{e,x}\cdot d{\rho }_e}{E_{e,i}}=\frac{{\beta }_{e,x}^{1/2}{\beta }_{p,x}^{1/2}}{{\Vert {\rho }_e-{\rho }_p\Vert }^3}\left({\rho }_e-{\rho }_p\right)\cdot d{\rho }_e\\ \frac{d{E}_{p,x}}{E_{p,i}}=d{\beta }_{p,x}=-\frac{r_e{F}_{p,x}\cdot d{\rho }_p}{E_{p,i}}=\frac{\mu {\beta }_{p,x}^{1/2}{\beta }_{e,x}^{1/2}}{{\Vert {\rho }_p-{\rho }_e\Vert }^3}\left({\rho }_p-{\rho }_e\right)\cdot d{\rho }_p\end{array},\begin{array}{c}x\equiv s,0,e\\ \mu \equiv \frac{m_{e,i}}{m_{p,i}}\end{array}\end{array}$.(54)

$F_{e\left|\right|p}$: Electrostatic force experienced by electron or proton in field of the other. $E_{e\left|\right|p,i}$: Self/rest energy of electron or proton in RS. $r_e$: CLF, classical electron radius. ${\beta }_{e\left|\right|p}$: Coulomb Factor of electron or proton. ${\rho }_{e\left|\right|p}$: Location vector of electron or proton, in reduced unit. $E_{e\left|\right|p,x}$: Self/rest energy of electron or proton at rest in field of the other. $m_{e\left|\right|p,i}$: Self/rest mass of electron or proton in RS. $x$: State indicator.

Set origin of RF to RS rest mass center of the BPS and approximate that

$\begin{array}{c}\begin{array}{c}s\equiv {\rho }_e-{\rho }_p\\ {\rho }_e\simeq +\frac{1}{1+\mu }s\\ {\rho }_p\simeq -\frac{\mu }{1+\mu }s\end{array}\to \begin{array}{c}\frac{d{\beta }_{e,x}}{ds}=\frac{{\beta }_{e,x}^{1/2}{\beta }_{p,x}^{1/2}}{\left(1+\mu \right)s^2}\\ \frac{d{\beta }_{p,x}}{ds}=\frac{{\mu }^2{\beta }_{p,x}^{1/2}{\beta }_{e,x}^{1/2}}{\left(1+\mu \right)s^2}\end{array}\to \frac{d}{ds}\left({\mu }^2{\beta }_{e,x}-{\beta }_{p,x}\right)=0\to \\ {\beta }_{p,x}=1-{\mu }^2\left(1-{\beta }_{e,x}\right)\to \frac{d{\beta }_{e,x}^{1/2}}{ds}=\frac{\sqrt{1-{\mu }^2\left(1-{\beta }_{e,x}\right)}}{2\left(1+\mu \right)s^2}\end{array}$.(55)

$s$: Distance vector between particles of the BPS in reduced unit.

Therefore,

$\begin{array}{c}{\beta }_{e,x}^{1/2}=\cosh \left[\frac{\mu }{1+\mu }\frac{1}{2s}\right]-\frac{1}{\mu }\sinh \left[\frac{\mu }{1+\mu }\frac{1}{2s}\right]\to \\ s_{\text{ESS}}=\frac{\mu }{2}\frac{\text{arccoth}{\mu }^{-1}}{1+\mu }\equiv s_0\end{array}$.(56)

$s_{\text{ESS}}$: Electric Schwarzschild distance between particles of the BPS, in reduced unit.

That is, minimal distance between parties of the BPS is a bit shorter of ½ of the CLF. Although expression for self energy of the electron appears different in comparison to that of Expression (21), the profile is actually quite similar to that of electron in SEF of proton, as can be seen in Figure 1 .

From Equation (54), with Equation (55),

$\begin{array}{c}\frac{d{P}_{e,f,u,e}}{d{t}_{e,s,u,e}}=-\frac{E_{e,i}}{r_e}{\beta }_{e,x}^{1/2}{\beta }_{p,x}^{1/2}\frac{\hat{s}}{s^2}\to \\ a_e-f_u{u}_e+f_e\hat{s}=0,f_u\equiv \frac{d\ln {\beta }_{u,e}}{d\tau },f_e\equiv \frac{{\beta }_{p,x}^{1/2}}{{\beta }_{e,x}^{1/2}{s}^2}\to \\ d\ln \frac{{\beta }_{e,x}}{{\beta }_{e,u}}=0\to \frac{{\beta }_{e,x}}{{\beta }_{e,u}}={ϵ}_e.\end{array}$(57)

${ϵ}_e$: Integration constant, total energy of electron in hydrogen atom in RS.

Therefore,

$\left(\frac{{\Delta }_s}{{\alpha }^2}+\frac{{ϵ}_e^2}{{\beta }_{e,x}^2}-1\right)\Psi =0,\alpha \equiv \frac{{\alpha }_e}{1+\mu },\begin{array}{c}\underset{s\to s_0}{\lim }\Psi =0\\ \underset{s\to \infty }{\lim }\Psi =0\end{array}$.(58)

${\Delta }_s$: Laplace operator with respect to $s$.

Expand the Coulomb Factor of the electron,

$s\gg s_0\to \begin{array}{c}\frac{1}{{\beta }_{e,x}^2}\simeq 1+\frac{2}{\left(1+\mu \right)\left(s-s_0\right)}+\frac{\gamma }{\left(1+\mu \right){\left(s-s_0\right)}^2}\equiv f_2\\ \gamma \equiv \frac{5-2\mu /\text{arctanh}\mu -{\mu }^2}{2\left(1+\mu \right)}.\end{array}$(59)

Define the base function set as

${\Phi }_{n,l,m}\equiv R_n^l{Y}_l^m,R_n^l\propto \frac{1}{s}{W}_n^{{\iota }_{n,l}}\left[2{\nu }_n\left(s-s_0\right)\right],\begin{array}{c}n>l\ge \left|m\right|\ge 0\\ n,l,m\in \text{Integers}\end{array}.$ (60)

Φ: Base function for WME (58).

Wherein,

$\begin{array}{l}{\nu }_n\equiv \alpha \sqrt{1-{ϵ}_n^2},{\iota }_{n,l}\equiv \sqrt{{\left(l+\frac{1}{2}\right)}^2-\frac{{\alpha }^2\gamma {ϵ}_n^2}{1+\mu }},\\ {ϵ}_n^2=\frac{{\left(1+\mu \right)}^2{n}^2}{2{\alpha }^2}\left(\sqrt{1+\frac{4{\alpha }^2}{{\left(1+\mu \right)}^2{n}^2}}-1\right)\end{array}$.(61)

WME (58) becomes

$\begin{array}{l}{\displaystyle \underset{n,l,m}{\sum }{c}_{n,l,m}\left(\frac{{ϵ}_e^2}{{\beta }_{e,x}^2}-f_2{ϵ}_n^2-\frac{l\left(l+1\right)}{{\alpha }^2}\left(\frac{1}{s^2}-\frac{1}{{\left(s-s_0\right)}^2}\right)\right)}|R_n^l{Y}_l^m\rangle =0\to \\ {\displaystyle \underset{n,n^{\prime }>l}{\sum }{c}_{n,l}\langle R_{n^{\prime }}^l|\frac{{ϵ}_e^2}{{\beta }_{e,x}^2}-f_2{ϵ}_n^2-\frac{l\left(l+1\right)}{{\alpha }^2}\left(\frac{1}{s^2}-\frac{1}{{\left(s-s_0\right)}^2}\right)}|R_n^l\rangle =0\end{array}$.(62)

In matrix form, for $l\ne 0$,

$\begin{array}{c}\left(\eta {ϵ}_e^2-\nu -l\left(l+1\right)\mu \right)c=0\\ \begin{array}{c}{\eta }_{j,k}\equiv \langle R_j^l|\frac{1}{{\beta }_{e,x}^2}|R_k^l\rangle \text{,}{\mu }_{j,k}\equiv \frac{1}{{\alpha }^2}\langle R_j^l|\frac{1}{s^2}-\frac{1}{{\left(s-s_0\right)}^2}|R_k^l\rangle \\ {\nu }_{j,k}\equiv \langle R_j^l|f_2|R_k^l\rangle {ϵ}_k^2,j,k\in \left\{n>l>0\right\}.\end{array}\end{array}$(63)

Therefore,

$\left(I{ϵ}_e^2-\lambda \right)c=0,\lambda \equiv {\eta }^{-1}\left(\nu +l\left(l+1\right)\mu \right)\to {ϵ}_{e,n,l}=\sqrt{{\lambda }_{n,l}}$.(64)

For $l=0$, with far field approximation of the second order,

$s\gg s_0\to \eta ϵc=0\text{,}\begin{array}{c}{\eta }_{j,k}\equiv \langle R_j^l|f_2|R_k^l\rangle \\ {ϵ}_{j,k}\equiv {\delta }_{j,k}\left({ϵ}_e^2-{ϵ}_k^2\right)\end{array}\to {ϵ}_{e,n,0}={ϵ}_n$.(65)

${\delta }_{j,k}$: Kronecker Delta function.

The result is not exact due to the far field approximations as well as the assumption on stativity of the rest mass center of the BPS.

Similar WME for energy state of a mass particle in SGF is obtainable if the particle is perceived as an EO or in association with such. Total energy of a mass particle in SGF is [4]

${ϵ}_g\equiv \frac{E_{f,u,g}}{E_i}=\frac{{\beta }_g}{{\beta }_u}=\frac{1}{\sqrt{1-u^2}}{\left(1-\frac{4}{\rho }\right)}^{1/4}\to u^2+\frac{{\beta }_g^2}{{ϵ}_g^2}=1\to -p_g^2+\frac{{ϵ}_g^2}{{\beta }_g^2}=1$. (66)

$E_{f,u,g}$: Total energy of a mass particle in motion in SGF as measured at rest at location of the particle in the field. $E_i$: Self energy of the mass particle in RS. ${ϵ}_g$: Total energy of a mass particle in motion in SGF as measured at rest at location of the particle in the field in reduced unit. ${\beta }_g$: Schwarzschild Factor for mass particle in SGF. $p_g$: Field momentum of a mass particle in reduced unit.

The equation is for total energy of a mass particle in motion in SGF under ISA, ignoring spin and self field interaction associated with motion of the particle and assuming stativity of the field. Self energy of mass particle in SGF is also shown in Figure 1 . Note that left bound of the domain is at radius of Schwarzschild Sphere (SS) of the field, which is four length units in terms of $r_g$ as defined below.

According to WM,

$\begin{array}{c}{P}_g\sim i\hslash {\nabla }_r\to p_g\sim \frac{1}{m_i{c}_i}\frac{i\hslash }{r_g}{\nabla }_{\rho },r_g\equiv \frac{G_i{M}_i}{c_i^2}\to \\ \frac{{\Delta }_{\rho }}{{\alpha }_g^2}+\frac{{ϵ}_g^2}{{\beta }_g^2}=\hat{I},{\alpha }_g\equiv \frac{G_i{M}_i{m}_i}{\hslash c_i}.\end{array}$(67)