Electron Mass Predicted From Substructure Stability in Electrodynamical Model

Modern physics has characterized spacetime, the interactions between particles, but not the nature of the particles themselves. Previous models of the electron have not specified its substance nor justified its cohesion. Here we present a relativistic electrodynamical model of the electron at rest, founded on natural interpretations of observables. Essentially intertwined positively and negatively charged subparticles revolve at light velocity in coplanar circular orbits, forming some coherent “envelope” and “nucleus”, possibly responsible for its wavelike and corpuscular behaviors, respectively. We show that the model can provide interpretations of fundamental constants, satisfy the Virial theorem, and exhibit cohesion and stability without invoking Poincaré stresses. Remarkably, the stability condition allows predicting electron mass, regarded as being a manifestation of its total (kinetic and potential) electromagnetic cohesion energy, and muon mass, directly from the substructure. Our study illustrates the possibility of constructing causal and objectively realist models of particles beneath the Compton scale. Finally, wave-corpuscle duality and the relation to quantum mechanics are discussed in the light of our electron model.

Introduction

Depending on the experiment, the most emblematic subatomic particle, the electron, has been found to interact as a point-like corpuscle in scattering experiments [1], or to behave as an extensible wave [2]. Elaborating on Bohr's interpretation of Quantum Mechanics [3], Heisenberg concluded that particles could neither be represented nor even apprehended by the human mind, and that only their abstract mathematical description existed [4]. For de Broglie however, “abstract presentations have no physical reality. Only the movement of elements localized in space, in the course of time, has physical reality” [5]. Hence, modern physics has identified with unprecedented precision the interactions and their underlying principles, has successfully described its environment, spacetime, but still lacks a characterization of the nature of its “objects,” the particles themselves.

Consequently, several kinds of electron models have been proposed: extended models [6], point-like models, and mixed models in which a point-like corpuscle follows an extended trajectory [7]. Early attempts included the spherical models of Abraham [8] and Lorentz [9], which led to theories of electromagnetic mass [10–13]. Spherical models soon evolved into the so-called ring models of Parson [14], Webster [15], Allen [16], and Compton [17], constituted of rotating infinitesimal charges and verifying the properties of classical magnetic moment and Compton scattering. Essential constraints however, such as electron cohesion and stability, could not be satisfied: new putative forces, denoted Poincaré stresses [18], were suggested to maintain the cohesion of the negatively charged electron. The abstract descriptions of quantum mechanical theories [19, 20] then successfully accounted for the wave-like behavior of the electron and probabilistically predicted [21] the values of most observables by considering a point-like particle, yet failed at interpreting fundamental constants or explaining how a point-like corpuscle could have spin or a finite energy density. Paradoxically, quantum mechanics revived geometric models when Schrödinger noticed within the Dirac equation itself a rapid oscillatory trembling motion, which he called Zitterbewegung (zbw) [22], exhibiting microcurrents arising at light velocity c. Surprisingly, the electron seemed to follow a helical trajectory of radius ƛ_c, the reduced Compton wavelength, surrounding the average travel direction (Figure 1A). Several such zbw models, identifying spin with orbital angular momentum, were interpreted classically [27–29]. Subsequent electrodynamical or hydrodynamical models involved fluids with spin [30], current loops of a certain thickness [31], Dirac-like Equations [32, 33], moving charged membranes [34], plasmoid fibers [35], or toroidal geometry [35, 36]. Wondering whether zbw could be a real phenomenon, Hestenes emphasized the need to investigate the electron substructure, suggested zbw could originate in the electron self-interaction [37], and showed zbw was compatible with the ring models [38].

FIGURE 1

Figure 1. Triolets and the helical trajectory. (A) In Schrödinger's Zitterbewegung (zbw) model, the wavefunction associated to the electron seems to revolve at light velocity along a helical trajectory of radius ƛ_c, the reduced Compton wavelength, surrounding the average travel direction, and to exhibit microcurrents. Quantum mechanics does not specify which forces could cause the electron, which is assumed to be point-like, to follow such a peculiar helical trajectory. (B) Triolets are colorless particles composed of three sparks, each bearing electric charge ±e/6 and a specific strong interaction color charge. Thus triolets bear electric charge ±e/6 or ±e/2 depending on their combination of sparks. They travel at light velocity c, possess angular momentum L_trlt, and triolets will thereafter be represented as upward or downward, filled or hollow triangles depending on their electric charge, as depicted here. (C) In our model, the electron is composed of triolets forming a nucleus and an envelope. It is conceivable that, in the absence of perturbation, the nucleus of the moving electron attracts envelope triolets and maintains them bound, thus explaining their helical trajectory. Conversely, envelope triolets would revolve at light velocity on an orbit of radius ƛ_c around the nucleus, exhibiting the zbw microcurrents, and guide the nucleus, as in pilot-wave theories [23–26], sensing the electromagnetic fields generated by the envelopes of other particles.

With the development of realist models of the electron emerged theories of electromagnetic mass. At first, the spherical models of Abraham [8], and Lorentz [9] seemed to fail to recover Einstein's relation E = mc² due to the appearance of a factor 4/3, but later proved to be compatible, once relativistic corrections were accounted for [12]. Stability of the sphere however still relied on Poincaré stresses or unknown surface tension [34], and electron mass could not be predicted from an objective criterion, but depended on the value taken by an arbitrary parameter, whose value is unconstrained, i.e., the radius of the sphere. Of note, the mass of subatomic particles is not predicted by quantum theories, and their values need to be inserted in calculations [12]. Most ring models [14–17] are prior to the discoveries of the spin, anomalous magnetic moment, and quantum mechanics. The ring model of Bergman and Wesley [31] exhibited cohesion and stability, but the expression for mass still involved an arbitrary parameter (i.e., width of current loop), and the substance constituting the electron remained indeterminate. More recently, Consa proposed a point-like electron following a toroidal trajectory [36], recovered mass independently of any arbitrary parameter, but did not specify how the trajectory developed nor demonstrated its stability. To our knowledge, the Virial theorem, which should be satisfied since the electron is a bound system, has not been considered in electron models. Potential energy is often equated to +mc², although cohesion potential energy should be negative for a bound system [as it is for the atom for instance [19]]. Kinetic energy is not usually accounted for, even though Lorentz [9], Hestenes [38] and others [e.g., [32]] noted the existence of a rotating motion and wondered whether kinetic energy did not contribute to rest mass. For Barut and Bracken, rest mass energy of the particle is the energy of the internal motion in the rest frame [29].

Hence, several issues remain to be addressed regarding the electron: for instance, which forces could cause the puzzling helical trajectory? What could be the nature of the substance constituting the electron? Could an electrodynamical description account for electron cohesion and stability? And could Lorentz' hypothesis advocating the electromagnetic origin of mass be simultaneously implemented from an objective criterion, instead of an arbitrary parameter?

In this study, we present a relativistic electrodynamical model of the electron at rest, in which charged subparticles follow definite trajectories. The model is based on two main hypotheses: (i) the existence of charged colorless subparticles called triolets, (ii) the assumption that triolets revolve at light velocity on coplanar circular orbits, constituting an envelope and nucleus, depending on their electromagnetic charges. As the electron is coherent, it is assumed that the model satisfies the Virial theorem. Constraints capturing the measured values of several observables (classical and anomalous magnetic moments, spin, Compton wavelength, kinetic energy) are formulated. Using Liénard-Wichert potentials, we then determine the specific kinds and numbers of triolets satisfying envelope and nucleus stability. Remarkably, we find that these kinds and numbers are precisely those that allow predicting electron mass and muon mass electromagnetically directly from the substructure, thus implementing Lorentz' hypothesis. Electron mass is effectively derived from an expression of substructure stability, which constitutes an objective criterion in our view. Our system also illustrates the possibility of constructing causal, local, objective, and realist models of particles beneath the Compton scale. Finally, we discuss novel perspectives suggested by the model, relative to the understanding of wave-corpuscle duality and to its relation to quantum theory.

Description of the Model and Hypotheses

In a previous study, we proposed that just six kinds of indestructible elementary subparticles denoted sparks, bearing electric charge ±e/6 and a specific strong interaction color charge, are necessary and sufficient to reconstruct all subatomic particles, so that sparks are conserved and reorganized across particle decays and annihilations [Avner, Boillot, Richard, submitted]. Since sparks are subject to both the strong and electromagnetic interactions, with the former dominating at short distances [20], groups of three sparks could presumably assemble beforehand to form composite colorless particles, thereafter called triolets, bearing charge +e/6, –e/6, +e/2, or –e/2 (Figure 1B). Henceforth, we shall suppose that the electron is exclusively composed of triolets, which travel at light velocity [7], exhibit some intrinsic angular momentum L_trlt, and being colorless, are submitted to electromagnetic and centrifugal forces only (hypothesis A).

Following de Broglie's proposition, we aim at constructing a plausible electrodynamical model of the electron at rest, in which positive and negative triolets form an electromagnetically bound system, exhibit the zbw microcurrents, and account for all experimentally measured observables. The electron is considered here as a particle of a certain extension, composed of revolving charged subparticles, the triolets, thereby exhibiting magnetic moment and intrinsic angular momentum (its spin) sensed by other particles. We know that the measured value of the electron magnetic moment is the sum of Bohr magneton μ_B = –eℏ/2m, predicted by both classical physics and quantum mechanics, where ℏ is reduced Planck constant, m the electron mass, and e the elementary charge, and an anomalous magnetic moment [39], which accounts for a small fraction a_anml ≃ 0.001159 of the previous and is only predicted by quantum electrodynamics [20]. Remarkably, the value of the classical magnetic moment of the electron can be derived by considering a charge (–e) revolving on a circular orbit of radius ƛ_c [19]. Hence, we reckoned the classical and anomalous magnetic moments could, respectively, be produced by two different components of the electron, namely a negatively charged envelope and a neutrally charged nucleus, also possibly responsible for the electron's wavelike and corpuscular behaviors, respectively. The peculiar helical trajectory of the electron predicted by zbw could then be naturally apprehended by considering that zbw describes the dynamics of envelope triolets, which are attracted and bound to the nucleus (Figure 1C). Electron spin could correspond to the sum of angular momenta of envelope triolets. Moreover, we shall regard electron mass as being a manifestation of the total electromagnetic cohesion energy E of the particle, as Lorentz hypothesized [9], through Einstein's formula m = E/c². The latter interpretation of the mass is naturally suggested by the observation that the muon possesses a mass ~206.77 times bigger than that of the electron, while its Compton wavelength is ~206.77 times smaller, as would be the case for a mass of electromagnetic origin, presenting a potential proportional to inverse distance.

The net electromagnetic forces acting on any particular envelope triolet should mostly depend on its surrounding triolets. The envelope could be organized into a complex structure, with triolets irregularly distributed along the orbits, or revolving at various radii, or experiencing fluctuations. To facilitate calculations however, we chose to make approximations and consider triolets at radial equilibrium rotating in the same direction on four coplanar circular orbits of different radii depending on their four different electromagnetic charges (hypothesis B, Figure 2). In our model, positive and negative nucleus triolets are intertwined to maintain their cohesion and could rotate along two close yet separate orbits due to the charged envelope. This could cause in turn a similar arrangement in the envelope, which would exhibit predominantly intertwined triolets, in spite of the excess of negative triolets. We are aware our model is only an approximation, even if we reckon that a collection of fluctuating ±e/6 and ±e/2 triolets traveling at light velocity could possibly converge toward such a configuration. Because of their stronger charges, ±e/2 triolets could be more tightly bound and form a condensed nucleus, while ±e/6 triolets would be bound more loosely and constitute the envelope.

FIGURE 2

Figure 2. Model of the electron at rest. In our simplified model, triolets rotate at light velocity in the same direction along four different coplanar circular orbits depending on their electric charge, constituting an envelope, and nucleus. Negative triolets are more numerous at the envelope, while the nucleus is neutrally charged. Due to the charged envelope, nucleus triolets are separated into two close orbits depending on their charge. Envelope triolets similarly revolve on separated orbits whose radii are close to the reduced Compton wavelength. Possible triolet configurations (triolet kinds and numbers, angular distributions, orbital radii) must fulfill constraints expressing radial stability and the measured values of charge, spin, magnetic moments, and mass. Due to consecutive negative triolets in the envelope, intertwined envelope triolets assemble into stretches separated by a distance d_env.

In addition, as the electron is a bound system whose inner potentials allegedly depend on position coordinates only and not velocities (justification is given below), the Virial theorem should be verified [40]: for inverse square law electromagnetic interactions, one typically has E = U/2 and E = –T, where T is the total internal kinetic energy and U the internal potential energy. Therefore, T and U should, respectively, amount to +mc² and −2mc², resulting in total internal energy E = T + U = –mc² corresponding to electron mass, the minus sign being indicative of a bound system. Finally, we shall admit that, for the electron at rest, envelope triolets approximately follow a circular trajectory of radius ƛ_c = ℏ/mc, as suggested by the classical derivation of Bohr's magneton, and by zbw-like models. Interpretations of fundamental constants associated to the electron, such as reduced Planck constant ℏ and fine-structure constant α = e²/4πε₀ℏc (where ε₀ designates vacuum permittivity), should also emerge from the model.

Formulation of the Model

Our system captures the measured values of charge, magnetic moments, spin and kinetic energy, and will be validated by showing that cohesion and stability can be satisfied, and potential energy (and thus electron mass) can be recovered. Let us here mathematically formulate the constraints: (i) a charge –e carried by N_env = N_env+ + N_env− triolets of charge ±e/n_env at the envelope; (ii) a classical magnetic moment μ_B generated by envelope triolets rotating at radii ρ_env+ = η_env+ƛ_c, ρ_env− = η_env−ƛ_c, and producing currents I_env+, I_env−; (iii) an anomalous magnetic moment a_anml·μ_B generated by N_nuc nucleus triolets (N_nuc+ = N_nuc−) of charge ±e/n_nuc rotating in the same direction as envelope triolets at radii ρ_nuc+ = η_nuc+ƛ_c, ρ_nuc− = η_nuc−ƛ_c, with momentum p_nuc+ ≃ p_nuc− = p_nuc and producing currents I_nuc+, I_nuc−; (iv) an internal kinetic energy T = ∑_ip_ic = +mc²; (v) a spin S_env = +ℏ/2 generated by envelope triolets of momentum p_env+, p_env−:

$\begin{array}{l}-e=e[\left(\frac{N_{env+}}{n_{env}}\right)-\left(\frac{N_{env-}}{n_{env}}\right)+\left(\frac{N_{nuc+}}{n_{nuc}}\right)\\ -\left(\frac{N_{nuc-}}{n_{nuc}}\right)],& (1)\end{array}$

$\begin{array}{ll}\frac{-e\hslash }{2m}=I_{env+}{\pi \rho }_{env+}^2+I_{env-}{\pi \rho }_{env-}^2,& (2)\end{array}$

$\begin{array}{ll}\frac{-a_{anml}e\hslash }{2m}=I_{nuc+}{\pi \rho }_{nuc+}^2+I_{nuc-}{\pi \rho }_{nuc-}^2,& (3)\end{array}$

$\begin{array}{l}{\displaystyle \sum _i{p}_{i}c}=(N_{env+}{p}_{env+}+N_{env-}{p}_{env-}\\ +N_{nuc}{p}_{nuc})c,& (4)\end{array}$

$\begin{array}{ll}\frac{\hslash }{2}=N_{env+}{\rho }_{env+}{p}_{env+}+N_{env-}{\rho }_{env-}{p}_{env-}.& (5)\end{array}$

The fact that the muon has same spin as the electron, despite possessing a smaller Compton wavelength and the same number of triolets according to our chemical theory [Avner, Boillot, Richard, submitted], suggests that the angular momentum of envelope triolets could be a constant ρ_env+p_env+ ≃ ρ_env−p_env− ≡ L_trlt,env, yielding from (5):

$\begin{array}{ll}\hslash =2N_{env}{L}_{trlt,env}.& (6)\end{array}$

L_trlt,env is possibly determined by the triangular substructure of envelope triolets made of three strongly interacting sparks, and could be at the basis of Planck's constant. Further constraints are also deduced (see Values of Observables) from Equations (1–5):

$\begin{array}{ll}{n}_{env}=N_{env-}-N_{env+},& (7)\end{array}$

$\begin{array}{ll}{n}_{env}=(N_{env-}{\eta }_{env-}-N_{env+}{\eta }_{env+}),& (8)\end{array}$

$\begin{array}{ll}{a_{anml}n}_{nuc}=N_{nuc+}({\eta }_{nuc-}-{\eta }_{nuc+}),& (9)\end{array}$

$\begin{array}{l}\frac{T_{env}+T_{nuc}}{m{c}^2}=1\simeq \frac{1}{b_{env}}\left(\frac{N_{env+}}{{\eta }_{env+}}+\frac{N_{env-}}{{\eta }_{env-}}\right)\\ \begin{array}{l}+\frac{N_{nuc}}{b_{nuc}{\eta }_{nuc}},\end{array}& (10)\end{array}$

$\begin{array}{ll}{b}_{env}=2N_{env},& (11)\end{array}$

where b_env and b_nuc are dimensionless numbers. Assuming that η_env+ ≃ η_env− ≡ η_env ≃ 1, we deduce (in Values of Observables) from Equations (10–11) that the kinetic energies of the envelope and nucleus are approximately equal T_env+ ≃ T_nuc ≃ +mc²/2, leading to relation:

$\begin{array}{ll}{b}_{nuc}{\eta }_{nuc}=2N_{nuc}.& (12)\end{array}$

Furthermore, the electromagnetic force acting on a nucleus triolet due to the envelope charge and current and the electromagnetic force exerted on an envelope triolet due to the net nucleus magnetic moment were derived but found to be negligible when compared to intra-component interactions. This suggests that each component is only loosely bound to the other, almost constituting an independent system, and thus verifies the Virial theorem independently (see Values of Observables), yielding for potential energies U_env ≃ U_nuc ≃ –mc² and total energies E_env ≃ E_nuc ≃ –mc²/2. The system of Equations (9–11) further allows to determine η_nuc+ (Values of Observables) for each value of N_nuc+:

$\begin{array}{ll}{\eta }_{nuc+}=\frac{N_{nuc+}}{b_{nuc}}\left[2-\frac{a_{anml}{b}_{nuc}{n}_{nuc}}{2N_{nuc+}^2}+\sqrt{4+{\left(\frac{a_{anml}{b}_{nuc}{n}_{nuc}}{2N_{nuc+}^2}\right)}^2}\right],& (13)\end{array}$

while η_nuc− is then given by Equation (9).

System cohesion and stability can be formulated by ensuring triolets are at radial equilibrium. As triolets are electrically charged and travel at light velocity, we use Liénard-Wichert potentials from relativistic electrodynamics [41] to express the radial components of electric field E_ij⊥ and magnetic field B_ij emitted by triolet T_j of charge q_j at retarded time t', radius ρ_j and retarded angle θ'_j, and sensed at distance R_ij—electromagnetic fields traveling at light velocity in vacuum—by triolet T_i arriving at the vertical (angle 0), radius ρ_i, at time t (Figure 3A). From known electrodynamical expressions [41] for these fields, using cylindrical unit vectors and coordinates, and Figure 3B, we deduce (Forces and Potentials):

$\begin{array}{ll}{\text{E}}_{\text{ij}\perp }=\frac{q_j\text{sin }{\gamma }_j}{4\pi {\epsilon }_0{R}_{ij}{\rho }_i{\left(1+\text{sin }{\gamma }_j\right)}^2}\widehat{\rho },& (14)\end{array}$

$\begin{array}{ll}{\text{B}}_{\text{i}\text{j}}=\frac{-q_j}{4\pi {\epsilon }_{0}c{R}_{ij}{\rho }_j{(1+\text{sin }{\gamma }_j)}^2}\widehat{\text{z}},& (15)\end{array}$

where R_ij and γ_j are defined by:

$\begin{array}{ll}{R}_{ij}^2={\rho }_i^2+{\rho }_j^2-2{\rho }_i{\rho }_j\text{cos }{\theta }_j^{\prime },& (16)\end{array}$

$\begin{array}{ll}\text{sin }{\gamma }_j=\frac{{\rho }_i}{R_{ij}}\text{sin }{\theta }_j^{\prime }.& (17)\end{array}$

Note that these fields depend on position coordinates only, not velocities, thereby justifying the use of the Virial theorem. We then derive expressions (Forces and Potentials) for the net radial Lorentz force F_ij⊥ due to triolet T_j exerted on triolet T_i belonging to the same component, and for the centrifugal force F_ctfg,i experienced by triolet T_i:

$\begin{array}{ll}{\text{F}}_{\text{ij}\perp }=\frac{q_i{q}_j}{4\pi {\epsilon }_0{R}_{ij}{\left(1+\text{sin }{\gamma }_j\right)}^2}\left[\frac{\text{sin }{\gamma }_j}{{\rho }_i}+\frac{1}{{\rho }_j}\right]\widehat{\rho },& (18)\end{array}$

$\begin{array}{ll}{\vec{\text{F}}}_{ctfg,i}=\frac{hc}{b_i{\rho }_i^2}\widehat{\rho },& (19)\end{array}$

where b_i stands for b_env (respectively, b_nuc) when T_i belongs to the envelope (resp. the nucleus). In the electron at rest, assuming triolets remain at radial equilibrium, the net radial component of the Lorentz force exerted by other triolets should compensate the centrifugal force. Neglecting the small contribution of the envelope onto the nucleus and vice-versa, and expressing equilibrium for triolet T_i along the radial direction and rearranging (Triolets at Radial Equilibrium), we obtain for the envelope and nucleus:

$\begin{array}{l}\frac{1}{\alpha }\simeq \frac{-b_{env}}{n_{env}^2}{\displaystyle \sum _{j\in env}^{N_{env}-1}}\frac{{\rho }_i^2\text{sgn}\left(i.j\right)}{R_{ij}{\left(1+\text{sin }{\gamma }_j\right)}^2}\left(\frac{\text{sin }{\gamma }_j}{{\rho }_i}+\frac{1}{{\rho }_j}\right)\\ \begin{array}{l}\equiv G_{i\in env}({\eta }_i),\end{array}& (20)\end{array}$

$\begin{array}{l}\frac{1}{\alpha }\simeq \frac{-b_{nuc}}{n_{nuc}^2}{\displaystyle \sum _{j\in nuc}^{N_{nuc}-1}}\frac{{\rho }_i^2\text{sgn}\left(i.j\right)}{R_{ij}{\left(1+\text{sin }{\gamma }_j\right)}^2}\left(\frac{\text{sin }{\gamma }_j}{{\rho }_i}+\frac{1}{{\rho }_j}\right)\\ \begin{array}{l}\equiv G_{i\in nuc}({\eta }_i),\end{array}& (21)\end{array}$

where sgn(i·j) is the sign of the product of the charges of triolets T_i and T_j, and α is the fine-structure constant, which is found to be related to the ratio between the centrifugal force and the net radial electromagnetic force experienced by any single triolet inside the electron. We assume positive and negative triolets are intertwined and uniformly distributed along the orbits except—as negative triolets are more numerous at the envelope—consecutive negative envelope triolets, which presumably repel to produce stretches of alternatively charged triolets separated by empty space (Figure 2). Let d_env designate the distance (using the number of missing triolets as units) between the stretches. The expressions under the sums in Equations (20–21) can be calculated by first considering the non-retarded angular positions θ_j of triolets distributed along the circular orbit, then by determining the corresponding retarded angles θ'_j, as illustrated in Retarted Angles, using Newton's recursion method for instance onto transcendental equation:

$\begin{array}{ll}{\left({\theta }_j-{\theta }_j^{\prime }\right)}^2=1-2\frac{{\rho }_i}{{\rho }_j}\text{cos }{\theta }_j^{\prime }+{\left(\frac{{\rho }_i}{{\rho }_j}\right)}^2,& (22)\end{array}$

and then deriving γ_j from Equation (17). Equations (20–21) will help us derive adequate values for N_env, N_nuc, n_env, n_nuc, b_env, b_nuc.

FIGURE 3

Figure 3. Geometric diagrams. (A) The influence of electromagnetic fields due to triolet T_j onto T_i: let triolets T_i and T_j belong to the same component (envelope or nucleus). Triolet T_j rotates at light velocity along circular orbit of radius ρ_j and arrives at angle θ_j at time t, but was at position T_j' at retarded angle θ_j' and time t' when it emitted electromagnetic fields that reached triolet T_i revolving along coplanar circular orbit of radius ρ_i and arriving at angle 0 (vertical y axis) at time t. The retarded electromagnetic fields can be expressed using Liénard-Wichert potentials. This figure applies to all envelope and nucleus triolets. (B) Diagram showing vectors and angles involved in the demonstration of the expressions of electromagnetic fields and potentials. (C) Diagram depicting the case ρ_j= ρ_i.

The potential energy due to interactions between the nucleus and envelope being negligible, the total potential energy of our system is approximately U_tot ≃ U_env + U_nuc, where U_env, and U_nuc are, respectively, the envelope and nucleus potential energies, which are evaluated in Potential Energy:

$\begin{array}{ll}{U}_{env}\simeq \frac{2\alpha m{c}^2}{n_{env}^2}{\displaystyle \sum _{i\in env}^{N_{env}}}{\displaystyle \sum _{j\ne i}^{N_{env}-1}}\frac{\text{sgn}(i.j)}{H_{ij}(1+\text{sin }{\gamma }_j)},& (23)\end{array}$

$\begin{array}{ll}{U}_{nuc}\simeq \frac{2\alpha m{c}^2}{n_{nuc}^2}{\displaystyle \sum _{i\in nuc}^{N_{nuc}}}{\displaystyle \sum _{j\ne i}^{N_{nuc}-1}}\frac{\text{sgn}(i.j)}{H_{ij}(1+\text{sin }{\gamma }_j)},& (24)\end{array}$

where H_ij = R_ij/ƛ_c. Assuming η_env+ ≃ η_env− ≃ 1, we demonstrate (Potential Energy) from b_env=2N_env (11) and Equation (20), which expresses the radial stability of every envelope triolet, that Equation (23) yields U_env≃ –mc². Likewise, assuming η_nuc+ ≃ η_nuc−, we demonstrate (Potential Energy) from b_nucη_nuc = 2N_nuc (12) and Equation (21), which expresses the radial stability of every nucleus triolet, that Equation (24) yields U_nuc ≃ –mc². Hence, we find: U_tot ≃ U_env + U_nuc ≃ −2mc², as expected from the Virial theorem, and recover electron mass. As substructure stability implies radial equilibrium for all envelope and nucleus triolets (20–21), it allows predicting electron mass. We find it remarkable that the same number of triolets allows to recover both substructure stability and electron mass.

Determination of Suitable Configurations

The problem then reduces to determining triolet configurations, i.e., sets of values for {n_env, n_nuc, N_env+, N_env−, N_nuc, b_env, b_nuc, η_env+, η_env−, η_nuc+, η_nuc−, d_env}, that verify radial equilibrium for every triolet and correctly predict the total energy. We shall estimate the stability and total energy in three different models of the envelope successively, each lying at a different level of approximation. The three models are: the one-orbit model, where all envelope triolets rotate on the same orbit η_env+≃η_env−≡η_env; the two-orbits model, where positively-charged envelope triolets revolve on orbit of radius η_env+ and negative triolets at radius η_env−; the n-orbits model where every envelope triolet i rotates on a circular orbit of specific but fixed radius η_i.

We shall first estimate the number of triolets N_env present in the envelope by considering the one-orbit model. Assuming η_env+ ≃ η_env− and η_nuc+ ≃ η_nuc−, we have R_ij ≃ 2ρ_icosγ_j (Figure 3C) both at the envelope and nucleus, and Equations (20–21) can be approximated to:

$\begin{array}{ll}\frac{1}{\alpha }\simeq \frac{-b_{env}}{2n_{env}^2}{\displaystyle \sum _{j\in env}^{N_{env}-1}}\frac{\text{sgn}\left(i.j\right)}{\text{cos }{\gamma }_j\left(1+\text{sin }{\gamma }_j\right)}\simeq G_{i\in env}({\eta }_i),& (25)\end{array}$

$\begin{array}{ll}\frac{1}{\alpha }\simeq \frac{-b_{nuc}}{2n_{nuc}^2}{\displaystyle \sum _{j\in nuc}^{N_{nuc}-1}}\frac{\text{sgn}\left(i.j\right)}{\text{cos }{\gamma }_j\left(1+\text{sin }{\gamma }_j\right)}\simeq G_{i\in nuc}({\eta }_i).& (26)\end{array}$

Recalling that b_env is related to N_env via b_env = 2N_env (11), and setting values for input parameters {n_env, d_env}, the iteration over N_env values in Equation (25) enabled us to determine values for b_env and N_env approximately verifying Equations (25) and (11) simultaneously. Due to the asymmetry in the arrangement of envelope triolets, we found these Equations were satisfied for different values of N_env depending on the triolet T_i under consideration. In the case n_env = 6, d_env = 0 for instance, we found positive triolets approximatively satisfied these conditions for N_env ≃ 108, while negative triolets did so for N_env ≃ 144, thus justifying the necessity of considering two distinct orbits in the envelope. Although these figures should be regarded as merely indicative, cases d_env = 1 and d_env = 2 also pointed at average value N_env = 126, corresponding to N_env+ = 60 and N_env−= 66, and we shall be considering only this case in the remainder of our analysis. For the nucleus, in the absence of a constraint like Equation (11), values for b_nuc and η_nuc satisfying Equations (26) and (12) simultaneously were determined for every iterated value of N_nuc. However, when accounting for the correction due to envelope current (first two terms, Triolets at Radial Equilibrium):

$\begin{array}{ll}{G}_{env>i\in nuc}\approx \frac{b_{nuc}\text{sgn}(i)}{n_{nuc}}\left[\frac{{\eta }_{nuc}^3}{2}+\frac{3{\eta }_{nuc}^4}{8}\right],& (27)\end{array}$

the best estimate appeared to be N_nuc = 18 (Table 1). Note that input values other than n_env = 6, n_nuc = 2 did not yield any possible solutions.

TABLE 1

Table 1. Stability and energy of various nucleus configurations.

Now, considering n_env = 6, n_nuc = 2, d_env = 2, η_env+ ≃ η_env− ≃ 1 (one-orbit model), and putting in the value obtained above for N_env, we evaluated potential energies U_env, U_nuc using Equations (13, 23–24) and found U_env = −0.997·mc² (Table 2), U_nuc ≃ −1.000·mc² (Table 1). The total potential energy therefore amounts to U_tot ≃ −1.997·mc², close to our expected result. Hence, recalling that kinetic energies satisfy T_env ≃ T_nuc ≃ +mc²/2, then T_tot ≃ T_env + T_nuc ≃ +mc², E_tot ≃ T_tot+ U_tot ≃ –mc², and the mass of the electron is deduced directly from our model substructure. Likewise, since the muon is seen as an excited state of the electron [6] according to our chemical theory [Avner, Boillot, Richard, submitted], presumably displaying a similar arrangement of triolets albeit on a smaller scale, muon mass can also be successfully calculated by replacing m by muon mass m_μ in expressions (23–24), or equivalently ƛ_c by the reduced muonic Compton Wavelength ƛ_muon.

TABLE 2

Table 2. Stability and energy of various envelope models.

We next evaluated the cohesion and stability of individual triolets. For the symmetric nucleus, we computed the right-hand side G_nuc of Equation (21) for every triolet; for n_nuc = 2, N_nuc = 18 for instance, accounting for the correction due to the envelope current, we obtained G_nuc(η_nuc+) ≃ 136.83, G_nuc(η_nuc−) ≃ 137.24 (Table 1). For the asymmetric envelope, which can be divided into six identical stretches of 21 triolets in the case N_env = 126, we computed the right-hand sides G_env of Equation (20) for every triolet belonging to the first stretch and compared the results with the left-hand side 1/α ≃ 137.036, which they should yield if triolets were truly at radial equilibrium. For the one-orbit model, setting n_env = 6, values of G_env disagreed with the expected value for all values of d_env (the case d_env = 2 is given in Table 3). Clearly, in these conditions at least, the centrifugal and net electromagnetic forces fail to compensate and to ensure radial equilibrium, one dominating over the other, and triolets would be moving radially as well as azimuthally. Therefore, we considered the two-orbits model with η_env+≃ 0.977, η_env− ≃ 1.023, for which we obtained an acceptable energy value (Table 2). Once again, we found that radial equilibrium was not verified for many envelope triolets, especially for consecutive negative triolets or those adjacent to them (Table 3). Hence, we decided to complicate our model again and considered envelope triolets orbiting at various but fixed radii ρ_i (n-orbits model) instead of the probably too general ρ_env+ and ρ_env−. We heuristically determined fixed radii exhibiting reasonable stability for all envelope triolets, then used an optimization algorithm, described in Optimization Algorithm, to make every triolet tend toward radial equilibrium, minimizing criterion K, the average absolute deviation from 1/α per triolet:

$\begin{array}{ll}K=\frac{1}{N_{env}}{\displaystyle \sum _i^{N_{env}}}|G_{i\in env}({\eta }_i)-\frac{1}{\alpha }|,& (28)\end{array}$

which effectively constitutes a measure of global stability of envelope triolets. Our algorithm converged toward a solution yielding acceptable energy and global stability (Table 2). The stability values G_env(η_i) of individual envelope triolets belonging to the first stretch in the n-orbits model are shown in Table 3: most values appeared to be close to 1/α. We found that our optimization algorithm nicely converged toward stable solutions. However, the latter were highly dependent on initial conditions, and a thorough optimization study is needed to ensure local minima are avoided.

TABLE 3

Table 3. Stability of individual envelope triolets.

Discussion

In this study, we presented a relativistic electrodynamical model of the electron based on natural interpretations of its associated observables. Our electron model is composed of triolets that revolve along coplanar circular orbits constituting an envelope and nucleus, which could be responsible for its wavelike and corpuscular behaviors, respectively. These two components would thus constitute a natural solution to wave-corpuscle duality. Capturing the values of charge, spin, magnetic moments, Compton wavelength and kinetic energy, we created a triolet-based configuration that verified cohesion and stability without invoking Poincaré stresses, and predicted electron and muon mass, defined as electromagnetic cohesion energy, directly from substructure stability. Importantly, our model accounts for kinetic energy and presents a negative cohesion potential energy, in agreement with the Virial theorem. In our model, the numbers of triolets in the envelope and nucleus are the adjusting parameters, and the same numbers are found to account both for substructure stability and electron mass. Notably, electron mass can be derived directly from an expression of substructure stability. Our study therefore implements Lorentz' hypothesis, which advocates the electromagnetic origin of mass, from an objective criterion, even if satisfaction of the criterion itself relies on two parameters, i.e., the numbers of triolets in the envelope and nucleus. Noteworthy, these parameters are not arbitrary, but instead are strongly constrained by several relations (11, 12, 20, 21, 27) that fix their values in our model. Altogether, we believe our study establishes that deterministic electrodynamical models of subatomic particles can be constructed beneath the Compton scale, in agreement with an objectively realist conception of physics.

Envelope triolets could also fluctuate radially or otherwise in time, possibly constituting a periodic wave that revolves at light velocity. This system has not been investigated here, but is of interest because this periodic wave could correspond to the wave associated to the electron, first imagined by de Broglie and later represented by wavefunction |ψ> in Schrödinger's wave mechanics or Dirac's quantum mechanics. It is conceivable that a wave made of envelope triolets, if it exists, attracts and drives the nucleus in the manner of the de Broglie-Bohm guiding wave [23, 24], sensing the electromagnetic fields generated by the envelopes belonging to other particles. Hence, envelope triolets could undulate and incarnate wavefunction |ψ>, whose concrete existence has recently been reconsidered [42]. Note further that nucleus triolets could also form a wave, reminiscent of the second wave described in de Broglie's double solution theory [23]. Specifically, triolets could propagate in a highly dynamical manner and experience irregular fluctuations, as in the hydrodynamical model of Bohm and Vigier [25]. Importantly, it has been suggested that solutions of this type could account both for quantum phenomena [26] and for quantum principles [37]. Bell also wrote that such solutions were compatible with the predictions of quantum mechanics [43]. Further, it is conceivable that such a complex envelope can exhibit several stable states, much like modes for a vibrating rope. These could correspond to the eigenstates of quantum mechanics. In the general case, the envelope would be in an unstable state, but could converge toward one of its eigenstates upon measurement, which could be conceived as the sum of interactions between system subparticles and apparatus subparticles. Such propositions constitute an interpretation of von Neumann's reduction of the wave packet [44], and would provide a possible solution to the measurement problem of quantum mechanics [45].

These considerations suggest that quantum theories, which encompass all subatomic phenomena and whose standard interpretation states that everything is intrinsically probabilistic, could eventually emerge [46, 47] from a relativistic electrodynamical description in agreement with the deterministic paradigm, which supports the causality principle, objective reality, and governs macroscopic physics. In this perspective, Schrödinger and Dirac Equations would constitute high-level descriptions of the dynamics of envelope triolets. Our study therefore provides new insight regarding the unification of the two apparently irreconcilable paradigms in physics: the deterministic and quantum paradigms.

Now, how exactly does the electron appear to be point-like in corpuscular interactions? How does our model relate to the observation that the electron seems spherical [48], or that its spin, charge and orbital components seem to be separable [49–51]? How would the moving electron, which exhibits a wave satisfying de Broglie relation p=h/λ, be described? Could our description be regarded as an attempt to create a corpuscular counterpart to wave mechanics? Could analogous electrodynamical models be similarly constructed for other subatomic particles [52]? Could our extended model of the electron bring insight to the nature of molecular bonding, or to the arrangement of electrons inside atoms? And finally, what would be the implications for the interpretation of quantum mechanics [45]? How would quantum properties, such as the existence of eigenstates, the measurement problem or entanglement, and quantum phenomena, such as the two-slits experiment or the one-dimensional potential well, be understood in the light of our model? We believe the aforementioned questions should stimulate discussion and foster novel investigations.

Methods

Values of Observables

Charge

The charge of the electron is given by the N_env+ triolets of charge (+e/n_env), the N_env− triolets of charge (–e/n_env), the N_nuc+ triolets of charge (+e/n_nuc) and the N_nuc− triolets of charge (–e/n_nuc):

$\begin{array}{l}-e=e\left[\left(\frac{-N_{env-}}{n_{env}}\right)+\left(\frac{N_{env+}}{n_{env}}\right)+\left(\frac{-N_{nuc-}}{n_{nuc}}\right)+\left(\frac{N_{nuc+}}{n_{nuc}}\right)\right].\end{array}$

Assuming the nucleus is neutrally charged (hypothesis B), implying N_nuc+= N_nuc−, we deduce:

$\begin{array}{ll}{n}_{env}=N_{env-}-N_{env+}.& (\mathrm{A1})\end{array}$

Nucleus and Envelope Orbits

Let us suppose triolets of charges (+e/n_env), (–e/n_env), (+e/n_nuc), (–e/n_nuc) revolve along four coplanar circular orbits of radii:

$\begin{array}{ll}\{\begin{array}{c}\begin{array}{c}{\rho }_{env+}={\eta }_{env+}{ƛ}_C\\ {\rho }_{env-}={\eta }_{env-}{ƛ}_C\end{array}\\ \begin{array}{c}{\rho }_{nuc+}={\eta }_{nuc+}{ƛ}_C\\ {\rho }_{nuc-}={\eta }_{nuc-}{ƛ}_C\end{array}\end{array}& (\mathrm{A2})\end{array}$

where ƛ_c = ℏ/mc is the reduced Compton wavelength, and η's are dimensionless real numbers.

Classical and Anomalous Magnetic Moments

Let us express the classical magnetic moment μ_B= –eℏ/2m = ∑_iI_iA_i= ∑_iQ_iA_i/t_i, where I_i is the current generated by triolet T_i, Q_i its charge, t_i = c/2πρ_i the time taken to go through a full orbit at light velocity c, and A_i the area formed by this orbit. The magnetic moment is due to a net charge (–e) made of N_env = N_env++ N_env− triolets revolving in the same direction along envelope orbits of radii ρ_env+ and ρ_env−:

$\begin{array}{l}{\mu }_B=\frac{-e\hslash }{2m}=\frac{Q_{env+}{A}_{env+}}{t_{env+}}+\frac{Q_{env-}{A}_{env-}}{t_{env-}},\\ \frac{-e\hslash }{2m}=\frac{N_{env+}e}{n_{env}}\frac{c\pi {\rho }_{env+}^2}{2\pi {\rho }_{env+}}+\frac{N_{env-}(-e)}{n_{env}}\frac{c\pi {\rho }_{env-}^2}{2\pi {\rho }_{env-}},\\ \frac{-e\hslash }{2m}=\frac{ec}{2n_{env}}(N_{env+}{\eta }_{env+}-N_{env-}{\eta }_{env-})\frac{\hslash }{mc},\\ \begin{array}{l}(N_{env-}{\eta }_{env-}-N_{env+}{\eta }_{env+})=n_{env}.\end{array}& (\mathrm{A3})\end{array}$

As the anomalous magnetic moment μ_nuc = –a_anml(eℏ/2m), with a_anml ≃ 0.001159, is relatively small, let us assume it is produced by an equal number N_nuc+ = N_nuc− of positive and negative triolets of charge (±e/n_nuc) revolving in the same direction as envelope triolets along nucleus orbits of slightly different radii due to the net envelope charge:

$\begin{array}{l}{\mu }_{nuc}=-a_{anml}\frac{e\hslash }{2m}=\frac{Q_{nuc+}{A}_{nuc+}}{t_{nuc+}}+\frac{Q_{nuc-}{A}_{nuc-}}{t_{nuc-}},\\ \begin{array}{l}{N}_{nuc+}({\eta }_{nuc-}-{\eta }_{nuc+})=a_{anml}{n}_{nuc}.\end{array}& (\mathrm{A4})\end{array}$

Virial Theorem

The virial theorem states that if a system remains bound, and if its inner potentials do not depend on velocities but only on positions, then the kinetic and potential energies take on definite shares in the total energy, depending on the degree n of the forces that apply. As the electron is a bound system, and as in our system the magnetic force will be found to depend on position coordinates ρ and γ only, the theorem applies and, for electromagnetic interactions in r⁻², it stipulates that:

$\begin{array}{ll}\{\begin{array}{c}T=m{c}^2\\ U=-2m{c}^2\\ E=T+U=-m{c}^2\end{array}& (\mathrm{A5})\end{array}$

where T, U, and E, respectively, designate the internal kinetic energy, internal potential energy, and total internal energy of the system. Note that the potential and total energies are negative, as they should be for a bound system.

Kinetic Energy

The kinetic energy is given by:

$\begin{array}{l}T=m{c}^2={\displaystyle \sum _i}{p}_{i}c=N_{env+}{p}_{env+}c+N_{env-}{p}_{env-}c\\ \begin{array}{l}+N_{nuc+}{p}_{nuc+}c+N_{nuc-}{p}_{nuc-}c,\end{array}& (\mathrm{A6})\end{array}$

suggesting:

$\begin{array}{ll}\{\begin{array}{c}\begin{array}{c}\begin{array}{c}{p}_{env+}=mc/K_{env+}\\ p_{env-}=mc/K_{env-}\end{array}\\ p_{nuc+}=mc/K_{nuc+}\end{array}\\ p_{nuc-}=mc/K_{nuc-}\end{array},& (\mathrm{A7})\end{array}$

where the K's remain to be determined, thus yielding from Equation (A6):

$\begin{array}{ll}1=\frac{N_{env+}}{K_{env+}}+\frac{N_{env-}}{K_{env-}}+\frac{N_{nuc+}}{K_{nuc-}}+\frac{N_{nuc+}}{K_{nuc-}}.& (\mathrm{A8})\end{array}$

Note that we may assume that nucleus triolets possess comparable momentum p_nuc+ ≃ p_nuc− = p_nuc, and that their orbit radius is approximately ρ_nuc+ ≃ ρ_nuc− = ρ_nuc, since (ρ_nuc+– ρ_nuc−) is very small according to Equation (A4).

Spin

Since particles as different as quarks and leptons (which possess different numbers of sparks according to our chemical model [Avner, Boillot, Richard, submitted]) share same spin, the latter can be interpreted as being the total angular momentum the particle conveys to the objects it encounters, i.e., the sum of the angular momenta of its envelope triolets. For the electron, assuming all triolets revolve in the same positive direction, it is written using Equations (A2, A6):

$\begin{array}{ll}S=+\frac{\hslash }{2}={\displaystyle \sum _i}{\rho }_i{p}_i=N_{env+}{\rho }_{env+}{p}_{env+}+N_{env-}{\rho }_{env-}{p}_{env-},& (\mathrm{A9})\end{array}$

$\begin{array}{ll}\frac{1}{2}=\frac{N_{env+}{\eta }_{env+}}{K_{env+}}+\frac{N_{env-}{\eta }_{env-}}{K_{env-}},& (\mathrm{A10})\end{array}$

Further, as the muon is composed of the same number of triolets as the electron according to our chemical model and exhibits a Compton length much smaller than that of the electron [Avner, Boillot, Richard, submitted], spin ℏ/2 is thus independent of the radii of triolets' orbits. A necessary and sufficient condition is then that variables K's be proportional to η's:

$\begin{array}{ll}\{\begin{array}{c}\begin{array}{c}\begin{array}{c}{K}_{env+}=b_{env+}{\eta }_{env+}\\ K_{env-}=b_{env-}{\eta }_{env-}\end{array}\\ K_{nuc+}=b_{nuc}{\eta }_{nuc+}\end{array}\\ K_{nuc-}=b_{nuc}{\eta }_{nuc-}\end{array}& (\mathrm{A11})\end{array}$

where b_env+, b_env−, b_nuc+, b_nuc− are values independent of radii, in order that the η's cancel out in Equation (A10), yielding:

$\begin{array}{ll}\frac{1}{2}=\frac{N_{env+}}{b_{env+}}+\frac{N_{env-}}{b_{env-}}.& (\mathrm{A12})\end{array}$

The angular momentum of triolet i is given by:

$\begin{array}{ll}{L}_i=p_i{\rho }_i=\frac{mc}{b_i{\eta }_i}.{\eta }_i\frac{\hslash }{mc}=\frac{\hslash }{b_i},& (\mathrm{A13})\end{array}$

implying for spin and kinetic energy:

$\begin{array}{ll}\frac{b_{env}}{2}=N_{env+}+N_{env-},& (\mathrm{A14})\end{array}$

$\begin{array}{ll}1=\frac{N_{env+}}{b_{env}{\eta }_{env+}}+\frac{N_{env-}}{b_{env}{\eta }_{env-}}+\frac{N_{nuc+}}{b_{nuc}{\eta }_{nuc+}}+\frac{N_{nuc-}}{b_{nuc}{\eta }_{nuc-}}.& (\mathrm{A15})\end{array}$

Definition of Planck's Constant

Supposing angular momentum L_trlt,env is a constant common to every envelope triolet, the expression for the spin, from Equation (A9), due to the envelope is:

$\begin{array}{l}\frac{\hslash }{2}=N_{env+}{L}_{trlt,env}+N_{env-}{L}_{trlt,env}=N_{env}{L}_{trlt,env},\end{array}$

and thus:

$\begin{array}{ll}\hslash =2N_{env}{L}_{trlt,env},& (\mathrm{A16})\end{array}$

meaning that the constant angular momentum L_trl,env common to every envelope triolet could be at the basis of Planck's constant.

Kinetic Energy of the Nucleus and Envelope

From Equations (A6, A7, A11), the kinetic energy of the nucleus is given by:

$\begin{array}{l}{T}_{nuc}=N_{nuc+}{p}_{nuc+}c+N_{nuc-}{p}_{nuc-}c\\ \begin{array}{l}=\frac{m{c}^2{N}_{nuc+}}{b_{nuc}}\left(\frac{1}{{\eta }_{nuc+}}+\frac{1}{{\eta }_{nuc-}}\right).\end{array}& (\mathrm{A17})\end{array}$

Likewise, the kinetic energy of the envelope is:

$\begin{array}{l}{T}_{env}=N_{env+}{p}_{env+}c+N_{env-}{p}_{env-}c,\\ \begin{array}{l}{T}_{env}=\frac{m{c}^2}{b_{env}}\left(\frac{N_{env+}}{{\eta }_{env+}}+\frac{N_{env-}}{{\eta }_{env-}}\right).\end{array}& (\mathrm{A18})\end{array}$

Now, assuming η_env+≃ η_env− ≃ 1 according to Schrödinger's Zitterbewegung, T_env becomes, using Equation (A14):

$\begin{array}{ll}{T}_{env}\simeq \frac{m{c}^2}{b_{env}}(N_{env-}+N_{env+})\simeq \frac{1}{2}m{c}^2,& (\mathrm{A19})\end{array}$

and thus:

$\begin{array}{ll}{T}_{nuc}=T-T_{env}\simeq \frac{1}{2}m{c}^2.& (\mathrm{A20})\end{array}$

The forthcoming study of the interactions between the nucleus and envelope will show that they are negligible compared to intra-component forces (nucleus onto itself, envelope onto itself). The two components therefore almost behave as two bound independent systems, and thus presumably obey the Virial theorem separately. Hence, since we have T_nuc ≃ T_env≃ mc²/2, we should also obtain U_nuc ≃ U_env≃ –mc² so that the total energies amount to: E_nuc ≃ E_env ≃ –mc²/2 and E_tot≃ –mc².

Determination of η_nuc+ and η_nuc-

In order to determine ${\eta }_{nu{c}^{+}}$ and ${\eta }_{nu{c}^{-}}$, considering Equations (A17) and (A20), we have:

$\begin{array}{ll}\left(\frac{1}{{\eta }_{nuc+}}+\frac{1}{{\eta }_{nuc-}}\right)\simeq \frac{b_{nuc}}{2N_{nuc+}}.& (\mathrm{A21})\end{array}$

The latter expression, together with Equation (A4), can allow us to determine η_nuc+ and η_nuc− in terms of N_nuc+, a_anml, n_nuc, and b_nuc:

$\begin{array}{l}\frac{b_{nuc}}{2N_{nuc+}}\simeq \frac{1}{{\eta }_{nuc+}}+\frac{1}{\left({\eta }_{nuc+}+\frac{a_{anml}{n}_{nuc}}{N_{nuc+}}\right)},\\ \frac{1}{{\eta }_{nuc+}}\left(\frac{b_{nuc}{\eta }_{nuc+}}{2N_{nuc+}}-1\right)=\frac{1}{{\eta }_{nuc+}}\left(\frac{1}{1+\frac{a_{anml}{n}_{nuc}}{N_{nuc+}{\eta }_{nuc+}}}\right),\\ 1=\left(1+\frac{a_{anml}{n}_{nuc}}{N_{nuc+}{\eta }_{nuc+}}\right)\left(\frac{b_{nuc}{\eta }_{nuc+}}{2N_{nuc+}}-1\right),\\ {\eta }_{nuc+}^2\left(\frac{b_{nuc}}{2N_{nuc+}}\right)+{\eta }_{nuc+}\left(\frac{a_{anml}{b}_{nuc}{n}_{nuc}}{2N_{nuc+}^2}-2\right)\\ -\left(\frac{a_{anml}{n}_{nuc}}{N_{nuc+}}\right)=0,\\ \Delta =4+{\left(\frac{a_{anml}{b}_{nuc}{n}_{nuc}}{2N_{nuc+}^2}\right)}^2,\end{array}$

and taking the positive solution, we find:

$\begin{array}{ll}{\eta }_{nuc+}=\frac{N_{nuc+}}{b_{nuc}}\left(2-\frac{a_{anml}{b}_{nuc}{n}_{nuc}}{2N_{nuc+}^2}+\sqrt{4+{\left(\frac{a_{anml}{b}_{nuc}{n}_{nuc}}{2N_{nuc+}^2}\right)}^2}\right),& (\mathrm{A22})\end{array}$

and η_nuc− can then be derived from Equation (A4).

Forces and Potentials

Centrifugal Force of a Triolet

Assuming triolets travel at light velocity, the centrifugal force [16] of triolet T_i, revolving along orbit of radius ρ_i = η_iƛ_c, is in cylindrical coordinates:

$\begin{array}{ll}{F}_{ctf,i}=\frac{p_i{v}_i}{{\rho }_i}=\frac{mc.c}{b_i{\eta }_i{\rho }_i}=\frac{\hslash c}{b_i{\rho }_i^2},& (\mathrm{B1})\end{array}$

where b_i stands for b_env (respectively, b_nuc) when T_i belongs to the envelope (resp. the nucleus). This expression applies both to nucleus and envelope triolets.

Electromagnetic Force Exerted on Nucleus Triolet i Due to Current at Envelope

The electromagnetic force exerted onto nucleus triolet i is given by the Lorentz force written using scalar potential V and vector potential A:

$\begin{array}{l}{\overrightarrow{\text{F}}}_{env\pm >i}=\frac{\text{sgn}(i).e}{n_{nuc}}[-\overrightarrow{\nabla }{V}_{env\pm >i}-\frac{\partial }{\partial t}{\overrightarrow{\text{A}}}_{env\pm >i}+c\widehat{\theta }\\ \times \left(\overrightarrow{\nabla }\times {\overrightarrow{\text{A}}}_{env\pm >i}\right)],& (\mathrm{B2})\end{array}$

if all triolets revolve in the same positive direction. The expressions for the scalar and vector potentials and their derivatives must be determined.

As a net charge (–e) circulates around the envelope, the scalar potential and vector potential, for r_nuc < r_env and cos θ = 0 (since the orbit is in the plane z = 0), are given [41] by:

$\begin{array}{ll}{V}_{env\pm >i\in nuc}=\frac{Q_{env\pm }}{4\pi {\epsilon }_0{r}_{env\pm }}{\displaystyle \sum _{l=0,2,4\dots }^{\infty }}{\left[P_l(0)\right]}^2{\left(\frac{{\rho }_i}{{\rho }_{env\pm }}\right)}^l,& (\mathrm{B3})\end{array}$

$\begin{array}{ll}{A}_{env\pm >i\in nuc}=\frac{{\mu }_0{I}_{env\pm }}{2}{\displaystyle \sum _{l=1,3,5\dots }^{\infty }}\frac{{\left[P_l^1(0)\right]}^2}{l\left(l+1\right)}{\left(\frac{{\rho }_i}{{\rho }_{env\pm }}\right)}^l,& (\mathrm{B4})\end{array}$

where the P_l(x) and P_l¹(x), respectively, designate the Legendre polynomials and associated Legendre polynomials, yielding:

$\begin{array}{ll}{V}_{env\pm >i\in nuc}\simeq \frac{Q_{env\pm }}{4\pi {\epsilon }_0}\left[\frac{1}{{\rho }_{env\pm }}+\frac{1}{4}\frac{{\rho }_i^2}{{\rho }_{env\pm }^3}+\frac{9}{64}\frac{{\rho }_i^4}{{\rho }_{env\pm }^5}\right],& (\mathrm{B5})\end{array}$

$\begin{array}{ll}\frac{\partial V_{env\pm >i\in nuc}}{\partial {\rho }_{nuc\pm }}\simeq \frac{Q_{env\pm }}{4\pi {\epsilon }_0}\left[\frac{1}{2}\frac{{\rho }_i}{{\rho }_{env\pm }^3}+\frac{9}{16}\frac{{\rho }_i^3}{{\rho }_{env\pm }^5}\right].& (\mathrm{B6})\end{array}$

Recalling μ₀= 1/(ε₀c²), v = c and Q_env± = ±N_env±e/n_env:

$\begin{array}{l}\frac{{\mu }_0{I}_{env\pm }}{2}=\frac{{\mu }_0{Q}_{env\pm }}{{2t}_{env\pm }}\\ \begin{array}{l}\simeq \frac{1}{2{\epsilon }_0{c}^2}\left(\frac{\pm N_{env\pm }e}{n_{env}}\right)\left(\frac{c}{2\pi {\rho }_{env\pm }}\right),\end{array}& (\mathrm{B7})\end{array}$

$\begin{array}{ll}{A}_{env\pm >i\in nuc}\simeq \frac{\pm N_{env\pm }e}{4\pi {\epsilon }_{0}c{n}_{env}}\left[\frac{1}{2}\frac{{\rho }_i}{{\rho }_{env\pm }^2}+\frac{3}{16}\frac{{\rho }_i^3}{{\rho }_{env\pm }^4}\right],& (\mathrm{B8})\end{array}$

$\begin{array}{l}\frac{\partial \left(A_{env\pm >i\in nuc}\widehat{\theta }\right)}{\partial t_{nuc}}=-A_{env\pm >i\in nuc}\frac{c}{{\rho }_{nuc}}\widehat{\rho }\\ \begin{array}{l}\simeq \frac{-(\pm N_{env\pm }e)}{4\pi {\epsilon }_0{n}_{env}}\left[\frac{1}{2{\rho }_{env\pm }^2}+\frac{3}{16}\frac{{\rho }_i^2}{{\rho }_{env\pm }^4}\right]\widehat{\rho },\end{array}& (\mathrm{B9})\end{array}$

$\begin{array}{l}\vec{\nabla }\times \vec{\text{A}}=\left|\begin{array}{ccc}\widehat{\rho }& \widehat{\theta }& \widehat{k}\\ \frac{\partial }{\partial \rho }& \frac{\partial }{\rho \partial \theta }& \frac{\partial }{\partial z}\\ 0& A_{env\pm >i\in nuc}& 0\end{array}\right|=\frac{\partial (A_{env\pm >i\in nuc})}{\partial {\rho }_{nuc}}\widehat{k}\\ \begin{array}{l}\simeq \frac{(\pm N_{env\pm }e)}{4\pi {\epsilon }_{0}c{n}_{env}}\left[\frac{1}{2{\rho }_{env\pm }^2}+\frac{9}{16}\frac{{\rho }_i^2}{{\rho }_{env\pm }^4}\right]\widehat{\text{k}}.\end{array}& (\mathrm{B10})\end{array}$

The electromagnetic force (B2) exerted on a nucleus triolet T_i by the envelope is then given by:

$\begin{array}{l}{\overrightarrow{\text{F}}}_{env\pm >i\in nuc}\simeq \frac{-\text{sgn}\left(i\right)}{n_{nuc}}\frac{\left(\pm N_{env\pm }{e}^2\right)}{4\pi {\epsilon }_0{n}_{env}{\rho }_{env\pm }^2}[\frac{1}{2}\frac{{\rho }_i}{{\rho }_{env\pm }}+\frac{3}{8}\frac{{\rho }_i^2}{{\rho }_{env\pm }^2}\\ +\frac{9}{16}\frac{{\rho }_i^3}{{\rho }_{env\pm }^3}]\widehat{\rho }.& (\mathrm{B11})\end{array}$

Electromagnetic Force Exerted on Envelope Triolet i Due to Current Flowing at Nucleus

According to Equation (A3), the magnetic moment due to the nucleus is:

$\begin{array}{ll}{\mu }_{nuc}=\frac{-a_{anml}e\hslash }{2m}=\frac{N_{nuc+}ec}{2n_{nuc}}({\rho }_{nuc+}-{\rho }_{nuc-}).& (\mathrm{B12})\end{array}$

The vector potential and its derivatives are given [41] by:

$\begin{array}{ll}{\vec{\text{A}}}_{nuc>i\in env}\simeq \frac{{\mu }_0}{4\pi }\frac{{\mu }_{nuc}}{{\rho }_i^2}\widehat{\theta },& (\mathrm{B13})\end{array}$

$\begin{array}{ll}{\vec{\text{A}}}_{nuc>i\in env}\simeq \frac{1}{8\pi {\epsilon }_{0}c}\frac{N_{nuc+}e}{n_{nuc}}\frac{({\rho }_{nuc+}-{\rho }_{nuc-})}{{\rho }_i^2}\widehat{\theta },& (\mathrm{B14})\end{array}$

$\begin{array}{l}\frac{\partial \left(A_{nuc>i\in env}\widehat{\theta }\right)}{\partial t_i}=-A_{nuc>i\in env}\frac{c}{{\rho }_i}\widehat{\rho }\\ \begin{array}{l}=\frac{-1}{8\pi {\epsilon }_0}\frac{N_{nuc+}e}{n_{nuc}}\frac{({\rho }_{nuc+}-{\rho }_{nuc-})}{{\rho }_i^3}\widehat{\rho }\end{array}& (\mathrm{B15})\end{array}$

$\begin{array}{l}\vec{\nabla }\times \vec{\text{A}}=\left|\begin{array}{ccc}\widehat{\rho }& \widehat{\theta }& \widehat{k}\\ \frac{\partial }{\partial \rho }& \frac{\partial }{\rho \partial \theta }& \frac{\partial }{\partial z}\\ 0& A_{nuc>i\in env}& 0\end{array}\right|=\frac{\partial (A_{nuc>i\in env})}{\partial {\rho }_i}\widehat{k}\\ \begin{array}{l}=\frac{(-1)}{4\pi {\epsilon }_{0}c}\frac{N_{nuc+}e}{n_{nuc}}\frac{({\rho }_{nuc+}-{\rho }_{nuc-})}{{\rho }_i^3}\widehat{\text{k}}.\end{array}& (\mathrm{B16})\end{array}$

As the net nucleus charge is zero, and using Equation (A4), the force is defined by:

$\begin{array}{l}{\vec{\text{F}}}_{nuc>i\in env}=\frac{\text{sgn}(i).e}{n_{env}}\left[-\frac{\partial }{\partial t}{\vec{\text{A}}}_{nuc>i}+c{\widehat{\theta }}_j\left({\vec{\nabla }\times \vec{\text{A}}}_{nuc>i}\right)\right],\\ {\vec{\text{F}}}_{nuc>i\in env}=\frac{3}{8\pi {\epsilon }_0}\frac{\text{sgn}(i)e^2{N}_{nuc+}}{n_{nuc}{n}_{env}}\frac{({\eta }_{nuc+}-{\eta }_{nuc-}){ƛ}_C}{{\rho }_i^3}\widehat{\rho },\\ \begin{array}{l}{\vec{\text{F}}}_{nuc>i\in env}=\frac{-3}{8\pi {\epsilon }_0}\frac{\text{sgn}(i)e^2}{n_{env}}\frac{a_{anml}{ƛ}_C}{{\rho }_i^3}\widehat{\rho }.\end{array}& (\mathrm{B17})\end{array}$

Electromagnetic Force Exerted on Triolet i at Radius ρ_i Due to Triolet j at Radius ρ_j

Every triolet experiences the fields emitted by all other triolets belonging to the same or adjacent orbit in the same component. Here we estimate the electromagnetic field and force exerted by a single triolet revolving on the same or adjacent orbit.

Let triolet T_j' (ρ_jsin θ'_j, ρ_jcos θ'_j) of charge q_j, revolving at light velocity on circular orbit of radius ρ_j, be positioned at angle θ'_j at retarded time t', and emitting an electromagnetic field received at time t by triolet T_i(0, ρ_i) of charge q_i revolving at light velocity on circular orbit of radius ρ_i, and arriving at angle θ_i = 0 on vertical axis y (Figure 3A). We have:

$\begin{array}{l}\overrightarrow{{\text{T}}_j{\text{T}}_i}\left(\begin{array}{c}-{\rho }_j\text{ sin }{\theta }_j^{\prime }\\ {\rho }_i-{\rho }_j\text{ cos }{\theta }_j^{\prime }\end{array}\right),\\ \begin{array}{l}{T_j{T}_i}^2={\rho }_i^2+{\rho }_j^2-2{\rho }_i{\rho }_j\text{cos }{\theta }_j^{\prime }\equiv R_{ij}^2,\end{array}& (\mathrm{B18})\end{array}$

$\begin{array}{ll}{\widehat{\text{n}}}_{ji}=\frac{\overrightarrow{{\text{T}}_{\text{j}}{\text{T}}_{\text{i}}}}{T_j{T}_i}\left(\begin{array}{c}-\frac{{\rho }_j}{R_{ij}}\text{sin }{\theta }_j^{\prime }\\ \frac{{\rho }_i-{\rho }_j\text{ cos }{\theta }_j^{\prime }}{R_{ij}}\end{array}\right),& (\mathrm{B19})\end{array}$

$\begin{array}{ll}{R}_{ij}=\sqrt{{\rho }_i^2+{\rho }_j^2-2{\rho }_i{\rho }_j\text{cos }{\theta }_j^{\prime }}.& (\mathrm{B20})\end{array}$

The trajectory, velocity and acceleration of triolet T_j are, respectively, given by:

$\begin{array}{ll}{\vec{\text{w}}}_j(t^{\prime })={\rho }_j\left(\text{sin }\omega t^{\prime }\widehat{\text{x}}+\text{cos }\omega t^{\prime }\widehat{\text{y}}\right),& (\mathrm{B21})\end{array}$

$\begin{array}{ll}{\vec{\text{v}}}_j(t^{\prime })={\rho }_j\omega \left(\text{cos }\omega t^{\prime }\widehat{\text{x}}-\text{sin }\omega t^{\prime }\widehat{\text{y}}\right),& (\mathrm{B22})\end{array}$

$\begin{array}{ll}{\vec{\text{a}}}_j(t^{\prime })={-\rho }_j{\omega }^2\left(\text{sin }\omega t^{\prime }\widehat{\text{x}}+\text{cos }\omega t^{\prime }\widehat{\text{y}}\right),& (\mathrm{B23})\end{array}$

with ω being the angular velocity, satisfying relations c = ρω and θ' = ωt'. Since v = c, β = v/c = 1, we also have:

$\begin{array}{l}{\widehat{\rho }}_{\text{j}}\left(\begin{array}{c}\text{sin }{\theta }_j^{\prime }\\ \text{cos }{\theta }_j^{\prime }\\ 0\end{array}\right),{\beta }_{\text{j}}\left(\begin{array}{c}\text{cos }{\theta }_j^{\prime }\\ -\text{sin }{\theta }_j^{\prime }\\ 0\end{array}\right),{\stackrel{.}{\beta }}_{\text{j}}\left(\begin{array}{c}-c\text{ sin }{\theta }_j^{\prime }/{\rho }_j\\ -c\text{ cos }{\theta }_j^{\prime }/{\rho }_j\\ 0\end{array}\right)\\ \begin{array}{l}=\frac{-c}{{\rho }_j}{\widehat{\rho }}_{\text{j}}\end{array}& (\mathrm{B24})\end{array}$

$\begin{array}{ll}g=1-{\beta }_{\text{j}}.{\widehat{\text{n}}}_{ji}=1-\text{cos }\left(\frac{\pi }{2}+{\gamma }_j\right)=1+\text{sin }{\gamma }_j.& (\mathrm{B25})\end{array}$

The electric and magnetic fields emitted by T_j and received by T_i are given [41] by:

$\begin{array}{ll}{E}_{\text{j}}=\frac{q_j}{4\pi {\epsilon }_0}\left[\frac{\left({\widehat{\text{n}}}_{ji}-{\beta }_{\text{j}}\right)\left(1-{\beta }^2\right)}{g^3{R}_{ij}^2}+\frac{{\widehat{\text{n}}}_{ji}\times \left[\left({\widehat{\text{n}}}_{ji}-{\beta }_{\text{j}}\right)\times {\stackrel{.}{\beta }}_{\text{j}}\right]}{c{g}^3{R}_{ij}}\right],& (\mathrm{B26})\end{array}$

$\begin{array}{ll}{B}_{\text{j}}=\frac{{\mu }_0{q}_j}{4\pi }\left[\frac{\left({\text{v}}_{\text{j}}\times {\widehat{\text{n}}}_{ji}\right)\left(1-{\beta }^2\right)}{g^3{R}_{ij}^2}+\frac{\left({\beta }_{\text{j}}\times {\widehat{\text{n}}}_{ji}\right)\left({\stackrel{.}{\beta }}_{\text{j}}.{\widehat{\text{n}}}_{ji}\right)+g{\stackrel{.}{\beta }}_{\text{j}}\times {\widehat{\text{n}}}_{ji}}{g^3{R}_{ij}}\right].& (\mathrm{B27})\end{array}$

From Figure 3B, it can be seen that:

$\begin{array}{ll}{\stackrel{.}{\beta }}_{\text{j}}.{\widehat{\text{n}}}_{ji}=\frac{c}{{\rho }_j}\text{cos }{\gamma }_j,& (\mathrm{B28})\end{array}$

$\begin{array}{ll}\text{cos}\frac{{\overline{\gamma }}_j}{2}|{\widehat{\text{n}}}_{ji}-{\beta }_{\text{j}}|,& (\mathrm{B29})\end{array}$

$\begin{array}{ll}\text{cos}\frac{{\overline{\gamma }}_j}{2}=\sqrt{\frac{1}{2}(1+\text{cos }{\overline{\gamma }}_j)}=\sqrt{\frac{1}{2}(1+\text{sin }{\gamma }_j)}.& (\mathrm{B30})\end{array}$

And thus:

$\begin{array}{ll}{\widehat{\text{n}}}_{ji}.\left({\widehat{\text{n}}}_{ji}-{\beta }_{\text{j}}\right)=|{\widehat{\text{n}}}_{ji}-{\beta }_{\text{j}}|.\text{cos}\frac{{\overline{\gamma }}_j}{2}=2{\text{ cos}}^2\frac{{\overline{\gamma }}_j}{2}=1+\text{sin }{\gamma }_j.& (\mathrm{B31})\end{array}$

From Equations (B25, B28) and identity: a×(b×c) = (a·c)b–(a·b)c, we deduce:

$\begin{array}{l}{\widehat{n}}_{ji}\times \left[\left({\widehat{n}}_{ji}-{\beta }_{\text{j}}\right)\times {\stackrel{.}{\beta }}_j\right]=\left({\stackrel{.}{\beta }}_j.{\widehat{n}}_{ji}\right)\left({\widehat{n}}_{ji}-{\beta }_j\right)-\left[{\widehat{n}}_{ji}.\left({\widehat{n}}_{ji}-{\beta }_j\right)\right]{\stackrel{.}{\beta }}_{\text{j}},\\ \begin{array}{l}{\widehat{n}}_{ji}\times \left[\left({\widehat{n}}_{ji}-{\beta }_j\right)\times {\stackrel{.}{\beta }}_{\text{j}}\right]=\frac{c}{{\rho }_j}\text{cos }{\gamma }_j\left({\widehat{n}}_{ji}-{\beta }_j\right)-\left(1+\text{sin }{\gamma }_j\right){\stackrel{.}{\beta }}_{\text{j}},\end{array}& (\mathrm{B32})\end{array}$

implying, since 1–β² = 0 and using Equations (B24, B25, B32):

$\begin{array}{l}{E}_{\text{j}}=\frac{q_j}{4\pi {\epsilon }_0}\left[\frac{{\widehat{\text{n}}}_{ji}\times \left[\left({\widehat{\text{n}}}_{ji}-{\beta }_j\right)\times {\stackrel{.}{\beta }}_j\right]}{c{g}^3{R}_{ij}}\right],\\ \begin{array}{l}{E}_{\text{j}}=\frac{q_j}{4\pi {\epsilon }_0{R}_{ij}{\rho }_j}\left[\left({\widehat{\text{n}}}_{ji}-{\beta }_j\right)\frac{\text{cos }{\gamma }_j}{{\left(1+\text{sin }{\gamma }_j\right)}^3}+{\widehat{\rho }}_j\frac{1}{{\left(1+\text{sin }{\gamma }_j\right)}^2}\right].\end{array}& (\mathrm{B33})\end{array}$

From Figure 3B, we also have:

$\begin{array}{ll}{\beta }_j\times {\widehat{\text{n}}}_{ji}=-\text{sin }\left(\frac{\pi }{2}+{\gamma }_j\right)\widehat{\text{z}}=-\text{cos }{\gamma }_j\widehat{\text{z}},& (\mathrm{B34})\end{array}$

$\begin{array}{ll}{\stackrel{.}{\beta }}_j\times {\widehat{\text{n}}}_{ji}=-\frac{c{\widehat{\rho }}_{\text{j}}}{{\rho }_j}\times {\widehat{\text{n}}}_{ji}=-\frac{c}{{\rho }_j}\text{sin }\left(\pi -{\gamma }_j\right)\widehat{\text{z}}=-\frac{c}{{\rho }_j}\text{sin }{\gamma }_j\widehat{\text{z}},& (\mathrm{B35})\end{array}$

yielding, using Equations (B27, B28) and μ₀ = 1/(ε₀c²):

$\begin{array}{l}{\text{B}}_{\text{j}}=\frac{{\mu }_0{q}_j}{4\pi }\left[\frac{-\frac{c}{{\rho }_j}\text{cos }{\gamma }_j\text{cos }{\gamma }_j-\frac{c}{{\rho }_j}(1+\text{sin }{\gamma }_j)\text{sin }{\gamma }_j}{{(1+\text{sin }{\gamma }_j)}^3{R}_{ij}}\right]\widehat{\text{z}},\\ \begin{array}{l}{\text{B}}_{\text{j}}=\frac{-q_j}{4\pi {\epsilon }_{0}c{R}_{ij}{\rho }_j{(1+\text{sin }{\gamma }_j)}^2}\widehat{\text{z}}.\end{array}& (\mathrm{B36})\end{array}$

The magnetic force is directed along ρ_i since B_j is along z. But to express the equilibrium we need to find the component of E_j along ρ_i, and thus we need:

$\begin{array}{ll}{\widehat{\text{n}}}_{ji}.{\widehat{\rho }}_{\text{i}}=\frac{1}{R_{ij}}\left({\rho }_i-{\rho }_j\text{ cos }{\theta }_j^{\prime }\right),& (\mathrm{B37})\end{array}$

$\begin{array}{ll}{\widehat{\rho }}_{\text{i}}.\left(-{\beta }_{\text{j}}\right)=\text{cos }\left({\theta }_j^{\prime }-\frac{\pi }{2}\right)=\text{sin }{\theta }_j^{\prime },& (\mathrm{B38})\end{array}$

$\begin{array}{ll}{\widehat{\rho }}_{\text{i}}.{\widehat{\rho }}_{\text{j}}=\text{ cos }{\theta }_j^{\prime },& (\mathrm{B39})\end{array}$

yielding from Equation (B33):

$\begin{array}{l}{E}_{ji\perp }=\frac{q_j}{4\pi {\epsilon }_0{R}_{ij}{\rho }_j}[\frac{1}{R_{ij}}\left({\rho }_i-{\rho }_j\text{cos }{\theta }_j^{\prime }\right)\frac{\text{cos }{\gamma }_j}{{\left(1+\text{sin }{\gamma }_j\right)}^3}\\ +\frac{\text{sin }{\theta }_j^{\prime }\text{cos }{\gamma }_j}{{\left(1+\text{sin }{\gamma }_j\right)}^3}+\frac{\text{cos }{\theta }_j^{\prime }}{{\left(1+\text{sin }{\gamma }_j\right)}^2}]\widehat{\rho }.& (\mathrm{B40})\end{array}$

This can be rearranged by expressing θ'_j as a function of γ_j and vice versa. From Equations (B19, B24):

$\begin{array}{l}\text{cos }{\gamma }_j=-{\widehat{\rho }}_{\text{j}}.{\widehat{\text{n}}}_{ji}=-\text{sin }{\theta }_j^{\prime }\left(-\frac{{\rho }_j}{R_{ij}}\text{sin }{\theta }_j^{\prime }\right)\\ -\text{cos }{\theta }_j^{\prime }\left(\frac{{\rho }_i-{\rho }_j\text{ cos }{\theta }_j^{\prime }}{R_{ij}}\right),\\ \begin{array}{l}\text{cos }{\gamma }_j=\frac{1}{R_{ij}}\left({\rho }_j-{\rho }_i\text{ cos }{\theta }_j^{\prime }\right).\end{array}& (\mathrm{B41})\end{array}$

Similarly, from Equation (B28):

$\begin{array}{l}\text{sin }{\gamma }_j\widehat{\text{z}}={\widehat{\rho }}_{\text{j}}\times {\widehat{\text{n}}}_{ji}=\left|\begin{array}{ccc}\widehat{x}& \widehat{y}& \widehat{z}\\ \text{sin }{\theta }_j^{\prime }& \text{cos }{\theta }_j^{\prime }& 0\\ -\frac{{\rho }_j}{R_{ij}}\text{sin }{\theta }_j^{\prime }& \frac{{\rho }_i-{\rho }_j\text{cos }{\theta }_j^{\prime }}{R_{ij}}& 0\end{array}\right|,\\ \begin{array}{l}\text{ sin }{\gamma }_j=\frac{{\rho }_i}{R_{ij}}\text{sin }{\theta }_j^{\prime }.\end{array}& (\mathrm{B42})\end{array}$

Relations (B41) and (B42) may be reversed:

$\begin{array}{ll}\text{sin }{\theta }_j^{\prime }=\frac{R_{ij}}{{\rho }_i}\text{sin }{\gamma }_j,& (\mathrm{B43})\end{array}$

$\begin{array}{ll}\text{cos }{\theta }_j^{\prime }=\frac{1}{{\rho }_i}({\rho }_j-R_{ij}\text{cos }{\gamma }_j).& (\mathrm{B44})\end{array}$

Then, using these to rearrange Equation (B40) and developing:

$\begin{array}{l}\frac{1}{R_{ij}^2}\left({\rho }_i-{\rho }_j\text{ cos }{\theta }_j^{\prime }\right)\left({\rho }_j-{\rho }_i\text{ cos }{\theta }_j^{\prime }\right)\\ =\frac{1}{R_{ij}^2}\left[{\rho }_i{\rho }_j{\text{sin}}^2{\theta }_j^{\prime }-R_{ij}^2\text{cos }{\theta }_j^{\prime }\right],\\ \frac{\text{sin }{\theta }_j^{\prime }}{R_{ij}}\left({\rho }_j-{\rho }_i\text{ cos }{\theta }_j^{\prime }\right)+\text{cos }{\theta }_j^{\prime }\left(1+\frac{{\rho }_i}{R_{ij}}\text{sin }{\theta }_j^{\prime }\right)\\ =\text{cos }{\theta }_j^{\prime }+\frac{{\rho }_j}{R_{ij}}\text{sin }{\theta }_j^{\prime },\end{array}$

we obtain using Equation (B42):

$\begin{array}{l}{E}_{ji\perp }=\frac{q_j\text{sin }{\theta }_j^{\prime }}{4\pi {\epsilon }_0{R}_{ij}^2{\left(1+\text{sin }{\gamma }_j\right)}^3}\left[\frac{{\rho }_i}{R_{ij}}\text{sin }{\theta }_j^{\prime }+1\right]\widehat{\rho },\\ E_{ji\perp }=\frac{q_j\text{sin }{\gamma }_j}{4\pi {\epsilon }_0{R}_{ij}{\rho }_i{\left(1+\text{sin }{\gamma }_j\right)}^2}\widehat{\rho }.& (\mathrm{B45})\end{array}$

The Lorentz force is then:

$\begin{array}{l}{F}_{\text{ij}\perp }=q_i\left(E_{ij\perp }+c{\widehat{\theta }}_i\times B_{ij}\right),\\ F_{\text{ij}\perp }=\frac{q_i{q}_j}{4\pi {\epsilon }_0{R}_{ij}{\rho }_i}\left[\frac{\text{sin }{\gamma }_j}{{\left(1+\text{sin }{\gamma }_j\right)}^2}\right]\widehat{\rho }\\ +\frac{q_i{q}_j}{4\pi {\epsilon }_0{R}_{ij}{\rho }_j}\left[\frac{1}{{\left(1+\text{sin }{\gamma }_j\right)}^2}\right]\widehat{\rho },\\ F_{\text{ij}\perp }=\frac{q_i{q}_j}{4\pi {\epsilon }_0{R}_{ij}{\left(1+\text{sin }{\gamma }_j\right)}^2}\left[\frac{\text{sin }{\gamma }_j}{{\rho }_i}+\frac{1}{{\rho }_j}\right]\widehat{\rho }.& (\mathrm{B46})\end{array}$

The scalar and vector Liénard-Wichert retarded electromagnetic potentials [41] are:

$\begin{array}{ll}{V}_{ij}=\frac{q_j}{4\pi {\epsilon }_0{\left(R_{ij}-{\beta }_{\text{j}}.{\text{R}}_{\text{i}\text{j}}\right)}_{rtrd}}=\frac{q_j}{4\pi {\epsilon }_0{R}_{ij}\left(1+\text{sin }{\gamma }_j\right)},& (\mathrm{B47})\end{array}$

$\begin{array}{ll}{\text{A}}_{\text{i}\text{j}}=\frac{{\mu }_0}{4\pi }{\left(\frac{q_j{v}_j{\widehat{\theta }}_j}{R_{ij}-{\beta }_{\text{j}}.{\text{R}}_{\text{i}\text{j}}}\right)}_{rtrd}=\frac{q_j{\widehat{\theta }}_j}{4\pi {\epsilon }_{0}c{R}_{ij}\left(1+\text{sin }{\gamma }_j\right)}.& (\mathrm{B48})\end{array}$

Approximation ρ_i=ρ_j. When making this approximation (one-orbit model), from Figure 3C, R_ij becomes:

$\begin{array}{ll}{R}_{ij}=2{\rho }_i\text{cos }{\gamma }_j.& (\mathrm{B49})\end{array}$

Note that if ρ_i= ρ_j, Equation (B46) then becomes:

$\begin{array}{ll}{\text{F}}_{\text{ij}\perp }=\frac{q_i{q}_j}{8\pi {\epsilon }_0{\rho }_i^2\text{cos }{\gamma }_j(1+\text{sin }{\gamma }_j)}\widehat{\rho }.& (\mathrm{B50})\end{array}$

Triolets at Radial Equilibrium

Equilibrium of Envelope Triolets

Envelope triolets are submitted to the centrifugal force (B1), the magnetic force due to the net nucleus magnetic moment (B17), and the net electromagnetic force due to the other envelope triolets (B46). Equilibrium for env– triolets can be written:

$\begin{array}{l}0=\frac{\hslash c}{b_{env}{\rho }_{env-}^2}+\frac{(-e)}{n_{env}}{\displaystyle \sum _j^{N_{env}-1}}\frac{e}{n_{env}}\frac{1}{4\pi {\epsilon }_0}\frac{\text{sgn}\left(j\right)}{R_{ij}{(1+\text{sin }{\gamma }_j)}^2}\\ \times \left[\frac{\text{sin }\gamma }{{\rho }_{env-}}+\frac{1}{{\rho }_j}\right]+\frac{3}{8\pi {\epsilon }_0}\frac{e^2}{n_{env}}\frac{a_{anml}{ƛ}_C}{{\rho }_{env-}^3}.\end{array}$

And rearranging to isolate the fine-structure constant:

$\begin{array}{l}\frac{4\pi {\epsilon }_0\hslash c}{e^2}=\frac{1}{\alpha }=\frac{b_{env}{\rho }_{env-}^2}{n_{env}}[{\displaystyle \sum _j^{N_{env}-1}\frac{1}{n_{env}}}\frac{\text{sgn}\left(j\right)}{R_{ij}{(1+\text{sin }{\gamma }_j)}^2}\\ \times \left(\frac{\text{sin }\gamma }{{\rho }_{env-}}+\frac{1}{{\rho }_j}\right)-\frac{3a_{anml}{ƛ}_C}{2{\rho }_{env-}^3}].& (\mathrm{C1})\end{array}$

Likewise, equilibrium for env+ triolets can be written:

$\begin{array}{l}\frac{1}{\alpha }=\frac{b_{env}{\rho }_{env+}^2}{n_{env}}[-{\displaystyle \sum _j^{N_{env}-1}\frac{1}{n_{env}}\frac{\text{sgn}\left(j\right)}{R_{ij}{(1+\text{sin }{\gamma }_j)}^2}\left(\frac{\text{sin }\gamma }{{\rho }_{env+}}+\frac{1}{{\rho }_j}\right)}\\ +\frac{3a_{anml}{ƛ}_C}{2{\rho }_{env+}^3}].& (\mathrm{C2})\end{array}$

Neglecting the term due to the nucleus magnetic moment, the equations become:

$\begin{array}{ll}\frac{1}{\alpha }=\frac{-b_{env}}{n_{env}^2}\left[{\displaystyle \sum _j^{N_{env}-1}}\frac{{\rho }_i^2\text{sgn}\left(i.j\right)}{R_{ij}{(1+\text{sin }{\gamma }_j)}^2}\left(\frac{\text{sin }\gamma }{{\rho }_i}+\frac{1}{{\rho }_j}\right)\right]\equiv G_{env}.& (\mathrm{C3})\end{array}$

The fine structure constant therefore appears to be naturally related to the ratio between the centrifugal force and the net electromagnetic force experienced by a single triolet. Making the ρ_i=ρ_j approximation (B49), we obtain:

$\begin{array}{ll}\frac{1}{\alpha }=\frac{-b_{env}}{2n_{env}^2}\left[{\displaystyle \sum _j^{N_{env}-1}}\frac{\text{sgn}\left(i.j\right)}{\text{cos }{\gamma }_j(1+\text{sin }{\gamma }_j)}\right].& (\mathrm{C4})\end{array}$

Equilibrium of Nucleus Triolets

Nucleus triolets are submitted to the centrifugal force (B1), the electromagnetic force due to the envelope (B11), and the net electromagnetic force due to the other nucleus triolets (B46). Equilibrium for nuc– triolets is thus written:

$\begin{array}{l}\frac{1}{\alpha }\simeq \frac{b{\rho }_{nuc-}^2}{n_{nuc}}[\frac{{\rho }_{nuc-}}{2n_{env}}\left(\frac{N_{env+}}{{\rho }_{env+}^3}-\frac{N_{env-}}{{\rho }_{env-}^3}\right)\\ +{\displaystyle \sum _j^{N_{nuc}-1}\frac{1}{n_{nuc}}\frac{\text{sgn}\left(j\right)}{R_{ij}{\left(1+\text{sin }{\gamma }_j\right)}^2}\left(\frac{\text{sin }\gamma }{{\rho }_{nuc-}}+\frac{1}{{\rho }_j}\right)}].& (\mathrm{C5})\end{array}$

Similarly we have for the nuc+ triolets:

$\begin{array}{l}\frac{1}{\alpha }\simeq \frac{b{\rho }_{nuc+}^2}{n_{nuc}}[\frac{{\rho }_{nuc+}}{2n_{env}}\left(\frac{N_{env+}}{{\rho }_{env+}^3}-\frac{N_{env-}}{{\rho }_{env-}^3}\right)\\ -{\displaystyle \sum _j^{N_{nuc}-1}\frac{1}{n_{nuc}}\frac{\text{sgn}\left(j\right)}{R_{ij}{\left(1+\text{sin }{\gamma }_j\right)}^2}\left(\frac{\text{sin }\gamma }{{\rho }_{nuc+}}+\frac{1}{{\rho }_j}\right)}].& (\mathrm{C6})\end{array}$

Neglecting the term due to the envelope current, the equations become:

$\begin{array}{ll}\frac{1}{\alpha }=\frac{{-b}_{nuc}}{n_{nuc}^2}\left[{\displaystyle \sum _j^{N_{nuc}-1}}\frac{{\rho }_i^2\text{sgn}\left(i.j\right)}{R_{ij}{(1+\text{sin }{\gamma }_j)}^2}\left(\frac{\text{sin }\gamma }{{\rho }_i}+\frac{1}{{\rho }_j}\right)\right]\equiv G_{nuc}.& (\mathrm{C7})\end{array}$

Making the ρ_i=ρ_j approximation (B49), we obtain:

$\begin{array}{ll}\frac{1}{\alpha }=\frac{-b_{nuc}}{2n_{nuc}^2}\left[{\displaystyle \sum _j^{N_{nuc}-1}}\frac{\text{sgn}\left(i.j\right)}{\text{cos }{\gamma }_j(1+\text{sin }{\gamma }_j)}\right].& (\mathrm{C8})\end{array}$

Also, the correction due to envelope current (first two terms) is:

$\begin{array}{l}{G}_{env>i\in nuc}\approx \frac{-b_{nuc}\text{sgn}(i)(-n_{env})}{n_{nuc}{n}_{env}}\left[\frac{{\rho }_{nuc}^3}{2{\rho }_{env}^3}+\frac{3{\rho }_{nuc}^4}{8{\rho }_{env}^4}\right]\\ \begin{array}{l}\approx \frac{b_{nuc}\text{sgn}(i)}{n_{nuc}}\left[\frac{{\eta }_{nuc}^3}{2}+\frac{3{\eta }_{nuc}^4}{8}\right].\end{array}& (\mathrm{C9})\end{array}$

Retarted Angles

Evaluating the Values of Retarded Angle θ_j' From Non-retarded Angle θ_j

If we suppose triolets are uniformly distributed along the circular orbits (this is certainly true of the nucleus since we have N_nuc+ = N_nuc−, but is an approximation in the case of the envelope, as there are more negative than positive triolets), then angle θ_j (expressed in radians) determining the position of the j^th triolet (starting at 1) at non-retarded time t on the orbit is defined by:

$\begin{array}{ll}{\theta }_{j\in nuc}=2\pi \frac{j}{N_{nuc}}.& (\mathrm{D1})\end{array}$

Note that, for the envelope, we also need to account for the empty space of length d_env (using the number of missing triolets as units) separating the n_env stretches of triolets, yielding for triolets T_j belonging to the first stretch:

$\begin{array}{ll}{\theta }_{j\in 1stStretch}=2\pi \frac{j}{\left(N_{env}+n_{env}{d}_{env}\right)}.& (\mathrm{D2})\end{array}$

To evaluate θ_j' determining the angular position T_j' at retarded time t' when the electromagnetic field was emitted toward triolet T_i, which arrives at angle 0 (vertical y axis) at time t to receive the field, we use the following relation, derived from Figure 3A:

$\begin{array}{ll}{R}_{ij}={\rho }_j\delta {\theta }_j={\rho }_j\left({\theta }_j-{\theta }_j^{\prime }\right).& (\mathrm{D3})\end{array}$

Then squaring Equations (B20) and (D3) and equating, we obtain:

$\begin{array}{ll}{({\theta }_j-{\theta }_j^{\prime })}^2=1-2\left(\frac{{\rho }_i}{{\rho }_j}\right)\text{cos }{\theta }_j^{\prime }+{\left(\frac{{\rho }_i}{{\rho }_j}\right)}^2.& (\mathrm{D4})\end{array}$

Given ρ_i, ρ_j, and θ_j, the retarded angles θ'_j may be numerically determined by recurrence, using a computer program that implements Newton method for instance, to resolve transcendental Equation (D4) for all triolets of angular position θ_j expressed in radians. The corresponding values of γ_j are then estimated using Equation (B42).

Potential Energy

Electric Potential Energy

By definition, the electric potential energies at the envelope and nucleus are defined by:

$\begin{array}{ll}{U}_{elec,env}={\displaystyle \sum _i^{N_{env}}}{\displaystyle \sum _{j\ne i}^{N_{env}-1}}{q}_i{V}_{ij}={\displaystyle \sum _{i\in env}}{\displaystyle \sum _{j\ne i}}\frac{q_i{q}_j}{4\pi {\epsilon }_0{R}_{ij}\left(1+\text{sin }{\gamma }_j\right)},& (\mathrm{E1})\end{array}$

$\begin{array}{ll}{U}_{elec,env}=\frac{\alpha m{c}^2}{n_{env}^2}{\displaystyle \sum _{i\in env}}{\displaystyle \sum _{j\ne i}}\frac{\text{sgn}(i.j)}{H_{ij}\left(1+\text{sin }{\gamma }_j\right)},& (\mathrm{E2})\end{array}$

where H_ij = R_ij/ƛ_c. Likewise, we have:

$\begin{array}{ll}{U}_{elec,nuc}=\frac{\alpha m{c}^2}{n_{nuc}^2}{\displaystyle \sum _{i\in nuc}}{\displaystyle \sum _{j\ne i}}\frac{\text{sgn}(i.j)}{H_{ij}\left(1+\text{sin }{\gamma }_j\right)}.& (\mathrm{E3})\end{array}$

Making the ρ_i = ρ_j approximation (B49), we obtain:

$\begin{array}{ll}{U}_{elec,env}=\frac{\alpha m{c}^2}{2n_{env}^2}{\displaystyle \sum _{i\in env}}{\displaystyle \sum _{j\ne i}}\frac{\text{sgn}(i.j)}{{\eta }_j\text{cos }{\gamma }_j\left(1+\text{sin }{\gamma }_j\right)},& (\mathrm{E4})\end{array}$

$\begin{array}{ll}{U}_{elec,nuc}=\frac{\alpha m{c}^2}{2n_{nuc}^2}{\displaystyle \sum _{i\in nuc}}{\displaystyle \sum _{j\ne i}}\frac{\text{sgn}(i.j)}{{\eta }_j\text{cos }{\gamma }_j\left(1+\text{sin }{\gamma }_j\right)}.& (\mathrm{E5})\end{array}$

Magnetic Potential Energy

The magnetic potential energy U_mag and electric potential energy U_elec are, respectively, the opposite of the magnetic work and electric work [41] given by:

$\begin{array}{ll}{W}_{mag}=\frac{1}{2{\mu }_0}{\displaystyle {\int }_{allspace}}{B}^{2}d\tau =-U_{mag},& (\mathrm{E6})\end{array}$

$\begin{array}{ll}{W}_{elec}=\frac{{\epsilon }_0}{2}{\displaystyle {\int }_{allspace}}{E}^{2}d\tau =-U_{elec}.& (\mathrm{E7})\end{array}$

Now, the vector expression relating the magnetic field to the electric field:

$\begin{array}{ll}\vec{B}=\frac{1}{c}\widehat{n}\times \vec{E}& (\mathrm{E8})\end{array}$

holds in relativistic electrodynamics with particles going at light velocity, yielding:

$\begin{array}{ll}{W}_{mag}=\frac{1}{2{\mu }_0{c}^2}{\displaystyle {\int }_{allspace}}{E}^{2}d\tau ,& (\mathrm{E9})\end{array}$

and since we know that c² = 1/ε₀μ₀, we have:

$\begin{array}{ll}{W}_{mag}=\frac{{\epsilon }_0}{2}{\displaystyle {\int }_{allspace}}{E}^{2}d\tau =W_{elec}.& (\mathrm{E10})\end{array}$

Therefore:

$\begin{array}{ll}{U}_{mag}=U_{elec}.& (\mathrm{E11})\end{array}$

Total Potential Energy

Neglecting the potential energy of the envelope acting on the nucleus U_env>nuc, and the potential energy of the nucleus acting on the envelope U_nuc>env, the electron potential energy is approximately:

$\begin{array}{ll}{U}_{tot}\simeq U_{env}+U_{nuc},& (\mathrm{E12})\end{array}$

where U_env is the envelope potential energy and U_nuc the nucleus potential energy. Using Equations (E2, E11), we obtain:

$\begin{array}{ll}{U}_{env}=U_{env,mag}+U_{env,elec}=2U_{env,elec},& (\mathrm{E13})\end{array}$

$\begin{array}{ll}{U}_{env}=\frac{2\alpha m{c}^2}{n_{env}^2}{\displaystyle \sum _{i\in env}}{\displaystyle \sum _{j\ne i}}\frac{\text{sgn}(i.j)}{H_{ij}\left(1+\text{sin }{\gamma }_j\right)},& (\mathrm{E14})\end{array}$

where H_ij= R_ij/ƛ_c. Likewise, using Equation (E.3) we have:

$\begin{array}{ll}{U}_{nuc}=\frac{2\alpha m{c}^2}{n_{nuc}^2}{\displaystyle \sum _{i\in nuc}}{\displaystyle \sum _{j\ne i}}\frac{\text{sgn}(i.j)}{H_{ij}\left(1+\text{sin }{\gamma }_j\right)}.& (\mathrm{E15})\end{array}$

Compatibility Between Potential Energies and Radial Equilibrium Equations

It can be shown that Equations (E14, E15) are compatible with Equations (C4, C8) if we assume η_nuc+≃ η_nuc− and η_env+ ≃ η_env− ≃ 1. Indeed, Equation (C4) becomes:

$\begin{array}{l}\left[{\displaystyle \sum _j^{N_{env}-1}}\frac{\text{sgn}(i.j)}{2\text{ cos }{\gamma }_j(1+\text{sin }{\gamma }_j)}\right]\simeq \frac{-n_{env}^2}{\alpha b_{env}}.\end{array}$

Then, by replacing the relation above into Equation (98), since η_env+≃η_env−≃1, we obtain:

$\begin{array}{l}{U}_{env}=\frac{2\alpha m{c}^2}{n_{env}^2}{\displaystyle \sum _{i\in env}}{\displaystyle \sum _{j\ne i}}\frac{\text{sgn}(i.j)}{2{\eta }_i\text{ cos }{\gamma }_j(1+\text{sin }{\gamma }_j)}\\ \simeq \frac{-2\alpha m{c}^2}{n_{env}^2}{N}_{env}\frac{n_{env}^2}{\alpha b_{env}},\\ U_{env}\simeq \frac{-2N_{env}m{c}^2}{b_{env}}.\end{array}$

Since b_env = 2N_env, this yields: U_env ≃ –mc² as expected. Likewise, Equation (C8) becomes:

$\begin{array}{l}\left[{\displaystyle \sum _j^{N_{nuc}-1}}\frac{\text{sgn}(i.j)}{2\text{ cos }{\gamma }_j(1+\text{sin }{\gamma }_j)}\right]\simeq \frac{-n_{nuc}^2}{\alpha b_{nuc}}.\end{array}$

Then, by replacing the relation above into Equation (E15), we obtain:

$\begin{array}{l}{U}_{nuc}=\frac{2\alpha m{c}^2}{n_{nuc}^2}{\displaystyle \sum _{i\in nuc}}{\displaystyle \sum _{j\ne i}}\frac{\text{sgn}(i.j)}{2{\eta }_i\text{ cos }{\gamma }_j(1+\text{sin }{\gamma }_j)}\simeq \frac{-2\alpha m{c}^2}{n_{nuc}^2}\frac{N_{nuc}}{{\eta }_{nuc}}\frac{n_{nuc}^2}{\alpha b_{nuc}},\\ U_{nuc}\simeq \frac{-2N_{nuc}m{c}^2}{b_{nuc}{\eta }_{nuc}}.\end{array}$

Since 2N_nuc = b_nuc η _nuc (A21), we obtain: U_nuc ≃ –mc² as expected.

Optimization Algorithm

An optimization algorithm has been devised and implemented to determine a set of optimum orbital radii for envelope triolets by minimizing average absolute deviation K, in the n-orbits model where each triolet possesses a specific radii η_i at the envelope. An approximate solution is determined heuristically before applying this algorithm. The algorithm next considers in turn every envelope triolet belonging to the first stretch, tries five different radii surrounding the current radius, and computes for each the stability of all envelope triolets. The radius yielding best overall stability is then attributed to the corresponding triolets in all six stretches. Once the procedure has been applied to all triolets of the first stretch, it is run again, considering five closer radii this time (thus slowly reducing the noise), until convergence toward an optimum solution is reached. The corresponding pseudocode is shown below. The algorithm was applied with the following values: delta = 0.00201, step = 0.00005, n_env= 6, N_env+ = 60, N_env− = 66, d_env = 2.

Data Availability Statement

All datasets generated for this study are included in the article/supplementary material.

Author Contributions

SA conceived the study, formed the hypotheses, constructed the model, wrote down and solved the equations, implemented the computations, and wrote the manuscript. FB helped solve the equations, devised the optimization algorithm, independently implemented the computations, and reviewed the manuscript. All authors approved the submitted version.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors wish to thank Patrick Richard (IFSTTAR, University Gustave Eiffel) for helpful advice and Gilles Salbert (IGDR, University of Rennes 1) for support and reading the manuscript.

References

1. Bourilkov D. Hint for axial-vector contact interactions in the data on e⁺e⁻ → e⁺e⁻ (γ) at center-of-mass energies 192-208 GeV. Phys Rev D. (2001) 64:071701. doi: 10.1103/PhysRevD.64.071701

CrossRef Full Text | Google Scholar

2. Davisson CJ, Germer LH. Diffractions of electrons by a crystal of Nickel. Phys. Rev. (1927) 30:705.

Google Scholar

3. Bohr N. The quantum postulate and the recent development of atomic theory. Nature. (1928) 121:580–90.

Google Scholar

4. Heisenberg W. Das naturgesetz und die struktur der materie, belser-presse, 1967. trans. Natural Law and The Structure of Matter. Belser: Warm Wind Books (1981).

Google Scholar

5. de Broglie L. Treize remarques sur divers sujets de physique théorique. Ann. Fond. L. de Broglie. (1976) 1:116–28.

Google Scholar

6. Dirac P. An extensible model of the electron. Proc. Royal Soc. Lond A. (1962) 268:57–67.

Google Scholar

7. Bunge M. A picture of the electron. Nuovo Cimento. (1955) 1:977–85.

Google Scholar

8. Abraham M. Prinzipien der dynamik des elektrons. Ann Phys. (1903) 1:105–79.

Google Scholar

9. Lorentz HA. The Theory of Electrons. 2nd ed. Leipzig: Teubner (1916).

Google Scholar

10. Fermi E. Über einen Widerspruch zwischen der elektrodynamischen und relativistischen theorie der elektromagnetischen masse. Phys Zeitschrift. (1922) 23:341–4.

11. Dirac P. Classical theory of radiating electrons. Proc Royal Soc Lond A. (1938) 167:148.

Google Scholar

12. Rohrlich F. Self-energy and stability of the classical electron. Am J Phys. (1960) 28:639.

Google Scholar

13. Schwinger J. Electromagnetic mass revisited. Found Phys. (1983) 13:373–83.

Google Scholar

14. Parson AL. A magneton theory of the structure of the atom smithson. Misc Collect. (1915) 65:1–80.

Google Scholar

15. Webster LW. The theory of electromagnetic mass of the parson magneton and other non-spherical systems. Phys Rev. (1917) 9:484–99.

Google Scholar

16. Allen HS. The case for a ring electron. Proc Phys Soc Lond. (1918) 31:49.

Google Scholar

17. Compton AH. The size and shape of the electron. Phys Rev. (1919) 14:247.

Google Scholar

18. Poincaré H. Sur la dynamique de l'électron. Rend Circolo Mat Palermo. (1906) 21:129–76.

Google Scholar

19. Eisberg R, Resnick R. Quantum Physics of Atoms, Molecules, Solids, Nuclei and Particles, 2nd ed. Hamilton, ON: John Wiley and Sons (1985).

Google Scholar

20. Griffiths DJ. Introduction to Elementary Particles, 2nd ed. Weinheim: Wiley-VCH. (2008). doi: 10.1002/9783527618460.index

CrossRef Full Text | Google Scholar

21. Born M. Quantenmechanik der Stoßvorgänge. Zeitschrift für Physik. (1926) 38:803–27.

Google Scholar

22. Schrödinger E. Uber die kraftefreie Bewegung in der relativistischen Quantum-mechanik, Sitz. Preuss. Akad. Wiss. Physik-Math. (1930)24:418.

Google Scholar

23. de Broglie L. La mécanique ondulatoire et la structure atomique de la matière et du rayonnement. J Phys Radium. (1927) 8:225–41.

Google Scholar

24. Bohm D. A suggested interpretation of the quantum theory in terms of hidden variables. Phys Rev. (1952) 85:166. doi: 10.1103/PhysRev.85.166

CrossRef Full Text | Google Scholar

25. Bohm D, Vigier JP. Model of the causal interpretation of quantum theory in terms of a fluid with irregular fluctuations. Phys Rev. (1954) 96:208–16.

Google Scholar

26. Bohm D, Dewdney C, Hiley BH. A quantum approach to the Wheeler delayed-choice experiment. Nature. (1985) 315:294–7.

Google Scholar

27. Weyssenhoff J. On two relativistic models of dirac's electron. Acta Phys Pol. (1947) 9:47–53.

PubMed Abstract | Google Scholar

28. Huang K. On the zitterbewegung of the dirac electron. Am J Phys. (1952) 20:479–84.

Google Scholar

29. Barut AO, Bracken AJ. Zitterbewegung and the internal geometry of the electron. Phys Rev D. (1981) 23:10.

Google Scholar

30. Barut AO, Zanghi N. Classical model of the dirac electron. Phys Rev Lett. (1984) 52:2009.

Google Scholar

31. Bergman DL, Wesley JP. Spinning charge ring model of electron yielding anomalous magnetic moment. Galilean Electrodynamics. (1990) 1:63–7.

Google Scholar

32. Pavšič M, Recami E, Rodrigues WA Jr, Maccarrone GD, Raciti F, Salesi G. Spin and electron structure. Phys Lett B. (1993)318:481–8.

Google Scholar

33. Salesi G, Recami E. Hydrodynamical reformulation and quantum limit of the barut-zanghi theory, found. Phys Lett. (1997) 10:533–46.

Google Scholar

34. Barut AO, Pavšič M. Radiation reaction and the electromagnetic energy-momentum of moving relativistic charged membranes. Phys Lett B. (1994) 331:45–50.

Google Scholar

35. Bostick WH. Mass, charge and current: the essence and morphology. Physics Essays. (1991) 4:45.

Google Scholar

36. Consa O. Helical solenoid model of the electron. Progr Phys. (2018) 14:80–89.

Google Scholar

37. Hestenes D. Quantum mechanics from self-interaction. Found Phys. (1985) 15:63–87.

Google Scholar

38. Hestenes D. The zitterbewegung interpretation of quantum mechanics. Found Phys. (1990) 20:1213–32.

Google Scholar

39. Odom B, Hanneke D, D'Urso B, Gabrielse G. New measurement of the electron magnetic moment using a one-electron quantum cyclotron. Phys Rev Lett. (2006) 97:030801.

Google Scholar

40. Goldstein H, Poole C, Safko J. Classical Mechanics, 3rd ed. Longman: Addison Wesley (2002). doi: 10.1119/1.1484149

CrossRef Full Text | Google Scholar

41. Zangwill A. Modern Electrodynamics. Cambridge: Cambridge University Press (2013).

Google Scholar

42. Pusey MF, Barrett J, Rudolph T. On the reality of the quantum state. Nat Phys. (2012) 8:475–8.

Google Scholar

43. Bell JS. On the impossible pilot wave. Found Phys. (1982) 12:989–99.

Google Scholar

44. von Neumann J. Mathematical Foundations of Quantum Mechanics. Princeton, NJ: Princeton Univ. Press (1955).

Google Scholar

45. Wheeler JA, Zurek WH. Quantum Theory and Measurement. Princeton, NJ: Princeton University Press (1983).

PubMed Abstract | Google Scholar

46. Einstein A. Physics and reality. J Franklin Institute. (1936) 221:349–82.

Google Scholar

47. Dürr D, Goldstein S, Zanghi N. Quantum equilibrium and the origin of absolute uncertainty. J Stat Phys. (1992) 67:843–907. doi: 10.1007/BF01049004

CrossRef Full Text | Google Scholar

48. Hudson JJ, Kara DM, Smallman IJ, Sauer BE, Tarbutt MR, Hinds EA. Improved measurement of the shape of the electron. Nature. (2011) 473:493–6. doi: 10.1038/nature10104

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Kim BJ, Koh H, Rotenberg E, Oh SJ, Eisaki H, Motoyama N, et al. Distinct spinon and holon dispersions in photoemission spectral functions from one-dimensional SrCuO₂. Nat Phys. (2006) 2:397–401. doi: 10.1038/nphys316

CrossRef Full Text | Google Scholar

50. Jompol Y, Ford CJB, Griffiths JP, Farrer I, Jones GAC, Anderson D, et al. Probing spin-charge separation in a tomonaga-luttinger liquid. Science. (2009) 325:597–601. doi: 10.1126/science.1171769

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Schlappa J, Wohlfeld K, Zhou KJ, Mourigal M, Haverkort MW, Strocov VN, et al. Spin-orbital separation in the quasi-one-dimensional Mott insulator Sr₂CuO₃. Nature. (2012) 485:82–5. doi: 10.1038/nature10974

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Eidelman S, Hayes KG, Olive KA, Aguilar-Benitez M, Amsler C, Asner D. (Particle Data Group). Review of particle physics. Phys Lett B. (2004) 592:1–5.

Google Scholar