arXiv:1409.8271v13 [physics.class-ph] 4 Jan 2019 eoiyi h ia hoyi h eoiyo ih,and Schr¨odinger light, motion. of mean velocity the the around is rapidly theory instantaneous oscillates elec- electron Dirac the rapidly-varying the Also, in a velocity [4]. has moment electron unlike dipole the theory, tric Dirac loop, the current in associated loop, as a current dipole classsical magnetic a ideal the elec- an with in as the not modeled coupling Although is is -orbit tron loop the electron. current of the classical analysis of semiclassical a representation that good clear a been has it [3], oc sacneuneo lsia electrodynamics. quantum classical the of explaining consequence a for as basis force a the provide if may expected correct forces magnetic present that the propose of to purpose the is to the paper compared Nonetheless, physics atomic force. gener- in Coulomb are weak which as remain, regarded forces interac- ally magnetic Coulomb only the by then If for tion, accounted electrodynamics. are classical forces in electric obvious not is decay.mechanics not does state ground radiate, the not proton. so will the and the the electron cancels the towards theory acceleration, exactly electron no Bohm the Undergoing state of - ground attraction the Broglie Coulomb in de force the with quantum motion in orbital decay, the Kepler-like radiative to obtains embodying mechanics atom classical as hydrogen of viewed na¨ıve mechanics. application a quantum be and Where can classical between which of difference existence the force, implies quantum [1] a mechanics quantum of mulation I LSIA PNIGPRIL MODEL PARTICLE SPINNING CLASSICAL II. ic ulcto fteDrcter fteelectron the of theory Dirac the of publication Since luil nlgeo h unu oc fBohmian of force quantum the of analogue plausible A for- Bohm - Broglie de the of potential quantum The h pe flgt upsto fsc oini h ai of basis the is motion such cha of point Supposition par a light. by of th of shown consisting canceled speed is each effectively the it particles be paper similar can this two In force between Coulomb proton. the force the quantum to predicts the electron because the atom, of hydrogen tion the of state ground omdlto ftemgei oc ihsailpro equa period spatial with force magnetic the of non-relativisti that modulation Co attributa shown to the be further may is opposite mechanics It but Bohmian equal electrodynamics. of circul is force their quantum force and the magnetic that mass the equal of ph have part the particles phase, by the modulated When sinusoidally force, motions. Coulomb radi a the contains luminall to motions two equal circulatory between their to force compared magnetic large The mechanics. quantum nted rgi omfruaino unu ehnc,th mechanics, quantum of formulation Bohm - Broglie de the In nepeaino unu ehnc 2 is [2] mechanics quantum of interpretation .INTRODUCTION I. osBh’ unu oc aeaCasclOrigin? Classical a Have Force Quantum Bohm’s Does Dtd uy3,2019) 30, July (Dated: [email protected] ai .Lush C. David rudsaeqatmfrei yrgncatoms. hydrogenic in the force to quantum similar state remarkably ground is that in- circumstances particles spin-spin certain mass) magnetic under equal (i.e., the similar for between teraction expectation representa- classical leads shown, charges a point be to circulating will as particles As charged of of tion factor function. wave phase the complex Dirac provide the of the to for phase interpretation [6] The physical proposed a elec- been the has [5]. to motion moment rise zitterbewegung around magnetic giving light spin, intrinsic of electron speed tron the the of re- at direction motion be the- the circular may Dirac a zitterbewegung as the The garded of motion “zitterbewegung.” oscillatory ory electron the called uial icltn hre oee,asmlrap- simpler a However, of charge. consisting circulating each particles luminally interac- mass a magnetic equal the two is between consider Rivas tion to by slightly charge only circulating extended luminally the of structure has approach. [8] this particles on interac- Dirac based The dipole obtained 1/2 been ideal spin magneton. also two an Bohr for of the Lagrangian field half tion magnetic of circulating time the moment the the to with of distance equal sufficient is field charge at acceleration defined, wavelength magnetic Compton is averaged reduced it the ra- where half for and of shows, motion also of Rivas off dius falling distance. magnitude, with equal the squarely of with inverse charge the magnitude stationary in a from consistent mo- of is far is field the sufficently field motion, and, of the electric undefined) of plane center average is the it the in where then (except tion, motion isotropic light, circulatory radial, of the purely of Rivas speed velocity the charge. the is if point that classical shown circulating has consisting luminally electron a an of of structure [7] Rivas electromagnetic classical by the in determination of origin the on an relies have electrodynamics, plausibly may theory Bohm’s nti ae h nlsso h lcrmgei field electromagnetic the of analysis the paper this In of force quantum the that paper, this of thesis The eaiemto ewe h atce leads particles the between motion relative c livresur a atwt magnitude with part law square inverse al g oigwt icltr oinat motion circulatory with moving rge oted rgi wavelength. Broglie de the to l s ieec ewe h circulatory the between difference ase tcascleetoyaissimilarly electrodynamics classical at lm oc.Ti assapossibility a raises This force. ulomb tr oin r lge u u of out but aligned are motions atory the ftemgei oc htoccurs that force magnetic the of t xcl acl h olm attrac- Coulomb the cancels exactly l otemgei oc fclassical of force magnetic the to ble -icltn hre o separation for charges y-circulating Zitterbewegung lcrni ttoayi the in stationary is electron e nepeainof interpretation 2 proach than Lagrangian mechanics will be used to char- field of a charge of equal value to the field source charge, acterize the interaction. Rather than attempting to ap- the magnetic interaction has nominally the same strength ply the interaction Lagrangian, it will simply be assessed as the Coulomb interaction. However, it is modulated by what is the magnetic force which one luminally circulat- both the relative orientation and the relative phase of ing charge exerts upon another charge undergoing similar the circular motions. Furthermore, the relevant phase motion, while neglecting the influence of the second cir- difference varies with interparticle separation, because it culating charge upon the first. One of the particles is thus is between the test particle present-time phase and the taken as the source of an electromagnetic field acting on field source particle retarded-time phase. a second, test, particle. As it has already been shown by Rivas that the average of the electric acceleration field is Coulomb-like in the far field, except in the plane of the motion where it is undefined, it is a small step to see III. RETARDED ELECTROMAGNETIC FIELD OF A RELATIVISTICALLY CIRCULATING that the magnetic force on a second luminally circulat- CHARGE ing charge will also follow an inverse-square depedence on interparticle separation in the far field. The magnetic field strength is equal (in Gaussian units) to the elec- In this section the electromagnetic field of a luminally- tric field strength, while the magnetic force magnitude circulating point charge is evaluated explicitly, so that is bounded by the magnetic field strength times the test in the next section the magnetic force on the luminally- particle velocity ratio to the speed of light. For luminally circulating test charge can be evaluated. Sections III A circulating charges, the magnetic force in the far field is and III B here reproduce Rivas’s in [7] with notation thus similar in distance dependence to the Coulomb force tailored to support the derivation of the magnetic force of between static charges. one luminally-circulating charge upon another of Section For further simplification it will also be assumed that IV. the radii of the circular motions of the two particles are The electromagnetic field due to a point charge in equal, and that the particles are widely separated com- arbitrary motion can be determined from the Li´enard- pared to the radius of the motions. Also, the angular Wiechert potentials. The electric and magnetic fields, E momenta of the circulatory motions of the two particles, and B, at a field point r and time t for an arbitrarily representing the spin vectors, will be assumed aligned. moving field-source charge qs obtained from the retarded Because the test particle is moving at the speed of light in Li´enard-Wiechert potentials may be expressed (in Gaus- a field that is of the same mean magnitude as a Coulomb sian units) as [9]

n × (n − β) × β˙ n − β qs E(r,t)= qs 3 + 3 , (1) "γ2 (1 − β · n) R2 # c  (1− β · n) R  ret ret  

(1), also undefined when β = n, can be referred to as the electric acceleration field, or as the radiative field 2 B(r,t) = [n × E]ret , (2) owing to its inverse R rather than inverse R dependence. However, the acceleration field becomes a solution to the where n is a unit vector in the direction of the field electromagnetic wave equation only in the far field. It point from the position of the field-source charge at the can also give rise to more rapidly decaying field terms retarded time tr = t − R/c, with R the magnitude of the that are dynamically relevant in the near field. displacement from the charge position at the retarded time to the field point r at time t. With the field-source charge velocity v, β = v/c, and γ ≡ (1 − β2)−1/2. The subscript “ret” indicates that quantities in the brackets A. Electric Field of a Classical Zitterbewegung are evaluated at the retarded time. The overdot repre- Particle sents differentiation with respect to t, so β˙ = a/c where a = aaˆ is the acceleration. To determine the magnetic force between two rela- The first term on the right hand side of Eq. (1), often tivistically circulating charges, the Li´enard-Wiechert re- referred to as the electric velocity field, vanishes in the tarded electric acceleration field is first evaluated for the limit of β approaching unity, except if β = n, when it is source charge undergoing circular motion. (To avoid the undefined due to vanishing of the inverse factor of (1 − singular condition where β = n, it is sufficient assume β · n)3. The second term on the right hand side of Eq. that the test charge center of circular motion is away from 3 the plane of circular motion of the source charge. This length. Also, for the stationary zitter particle, the accel- 2 also ensures that the electric velocity field vanishes, for eration magnitude a = c|β˙ | = c /λ0. Defining the plane the luminal field-source charge.) Next, the magnetic field of the motion as the xy or z = 0 plane of a right-handed can be evaluated using Eq. (2), and the magnetic force Cartesian coordinate system (x,y,z), then for the case on the luminally-circulating test charge follows straight- of counterclockwise zitter motion, kˆ = [cos ωt, sin ωt, 0] forwardly. with ω = c/λ0. For brevity, and after [10], it will be convenient to refer If R = r − k = rrˆ − λ0kˆ is the displacement from to the luminal circular motion of the point charge as the the instantaneous charge position to the field point, then “zitter” motion, and to a charge undergoing such motion 2 2 1/2 2 1/2 R ≡ |R| = (r − 2λ0r · kˆ + λ0 ) = r(1 − 2ǫrˆ · kˆ + ǫ ) as a zitter particle. In spite of the luminal circular motion with ǫ = λ /r, and n ≡ R/R becomes of its charge, the zitter particle is considered stationary 0 when its center of motion is fixed. For the stationary zitter particle, it is useful to sup- rˆ − ǫkˆ pose that the center of the zitter motion is the origin of n = . (3) ˆ 2 1/2 a Cartesian coordinate system, so that r of Eqs. (1) and (1 − 2ǫrˆ · k + ǫ ) (2) can represent the vector displacement from the cen- ter of the zitter motion of the source charge to the field The field-source charge acceleration is toward the cen- ˆ point. Then if k is the vector displacement from the cen- ter of the circular motion, so aˆ = −k, and for counter- clockwise circular motion β × aˆ = zˆ, and so ter to the charge instantaneous position, k = λ0kˆ where kˆ is a unit vector and λ0 is the radius of the circular mo- tion. With angular momentum due to the zitter motion rˆ × kˆ given as L = λ0γmv with v = c, taking γm = me where (n − β) × aˆ = − − zˆ. (4) ˆ 2 1/2 me is the observed electron mass, and λ0 = ~/2mec, (1 − 2ǫrˆ · k + ǫ ) obtains that the electron intrinsic angular momentum is ~/2. The diameter of the circulatory motion that has an- Also for counter-clockwise circular motion, kˆ×zˆ = −β, gular momentum ~/2 is thus one reduced Compton wave- and cross-multiplying Eq. (4) by Eq. (3) on the left yields

ˆ ˆ ˆ rˆ · k rˆ − ǫk − k − ǫrˆ rˆ × zˆ + ǫβ n × ((n − β) × aˆ)= − − , (5)  1− 2ǫrˆ · kˆ+ ǫ2 (1 − 2ǫrˆ · kˆ + ǫ2)1/2

where the vector identity a × (b × c) = (a · c) b − the electric acceleration field term from Eq. (1) can now (a · b) c has been used. be written in terms of the distance r from the zitter parti- 2 Recalling that β˙ = aaˆ/c, with a = c /λ0 for the sta- cle center of motion (rather than in terms of the distance tionary zitter particle and that R = r(1 − 2ǫrˆ · kˆ + ǫ2)1/2, R from the luminally-circulating charge position) as

rˆ · kˆ rˆ − ǫkˆ − kˆ − ǫrˆ −qs rˆ × zˆ + ǫβ Eacc(r,t)= + , (6) 3 2 3/2 2 λ0 (1 − β · n) r  (1 −2ǫrˆ · kˆ + ǫ ) 1 − 2ǫrˆ · kˆ + ǫ   

where all quantities on the right hand side are evalu- of the time-averaged electric acceleration field intensity, ated at the retarded time. it is useful to consider the behavior of the field on the z- Although the electric acceleration field due to the zit- axis, i.e., perpendicular to the plane of the zitter motion ˆ ter particle has an explicit inverse dependence on the of the source charge. When rˆ = zˆ, then rˆ·k = rˆ ×zˆ = 0, β n distance r, the factors of ǫ = λ0/r give dependencies of · = 0, and so inverse higher powers of r, including, but not limited to, inverse r2. To see clearly the inverse square of distance character 4

compared to the zitter radius at all points off the z = 0 plane. On the z = 0 plane, (1 − β · n) can vanish and q kˆ + ǫzˆ ǫβ E r z s the electric field average cannot be calculated. acc( = zˆ,t)= 2 3/2 − 2 . (7) λ0z "(1 + ǫ ) 1+ ǫ #

The unit vectors kˆ and β average to zero over a cycle B. Magnetic Field of a Classical Zitterbewegung of the zitter motion, and on the z-axis ǫ = λ0/z. The Particle average electric acceleration field on the z-axis is thus The magnetic field due to the zitter particle charge acceleration is evaluated using Eqs. (2) and (6), but it qs zˆ < Eacc(r = zzˆ) >= . (8) is convenient to use the intermediate result of Eq. (5). z2 (1 + ǫ2)3/2   Taking the cross product of Eq. (5) with n (as given by Eq. (3)) on the left obtains that As z becomes large, 1 + ǫ2 approaches unity, and so the long distance behavior of the time-averaged electric acceleration field of the fixed zitter particle, perpendicu- (1 − ǫ2)rˆ × kˆ lar to the plane of the zitter motion, is indistinguishable n × [n × ((n − β) × aˆ)] = ˆ 2 3/2 from the electric field of a static charge of the same mag- (1 − 2ǫrˆ · k + ǫ ) nitude. (Eq. (8) here is equivalent to Eq. (6.3) of [11].) rˆ − ǫkˆ × (rˆ × zˆ)+ ǫ rˆ − ǫkˆ × β Rivas [11] has numerically integrated the electric ac- − , (9) celeration field of the zitter particle to find the average   (1 − 2ǫrˆ · kˆ +ǫ2)  at points off the z-axis, and at near as well as distant separation from center of the zitter motion. The aver- and with the vector identity a × (b × c) = (a · c) b − age field is purely radial and isotropic at distances large (a · b) c,

ˆ ˆ ˆ (1 − ǫ2)rˆ × kˆ (rˆ · zˆ) rˆ − zˆ − ǫ k · zˆ rˆ − k · rˆ zˆ + ǫ rˆ − ǫk × β n × [n × ((n − β) × aˆ)] = − . (10) (1 − 2ǫrˆ · kˆ + ǫ2)3/2  (1 −2ǫrˆ ·kˆ + ǫ2)   

2 2 1/2 With a = c /λ0, β˙ = aaˆ/c, and R = r(1 − 2ǫrˆ · kˆ + ǫ ) , the magnetic acceleration field of the zitterbewegung particle can now be found from Eqs. (2) and (6) to be

2 (rˆ · zˆ) rˆ − zˆ − ǫ kˆ · zˆ rˆ − kˆ · rˆ zˆ + ǫ rˆ − ǫkˆ × β qs (1 − ǫ )rˆ × kˆ Bacc = − , (11) 3 2 2    2 3/2   λ0 (1 − β · n) r "(1 − 2ǫrˆ · kˆ + ǫ ) (1 − 2ǫrˆ · kˆ + ǫ )  

where the right-hand side is evaluated at the retarded on a test particle in arbitrary motion, due to the accel- time. eration field of a source particle undergoing the circular Since for luminal charge velocity the electric velocity zitter motion of radius λ0 can be evaluated by taking a field vanishes everywhere that it is defined, and the ac- cross product of Eq. (11) with βt, and using again the celeration field is undefined under the same conditions as vector identity a×(b × c) = (a · c) b−(a · b) c, to obtain the velocity field, Eq. (11), where it is well defined, de- that scribes the complete magnetic field of the zitter particle.

IV. RADIALLY-DIRECTED INVERSE-SQUARE LAW MAGNETIC FORCE BETWEEN DIRAC PARTICLES

The force on a test charge qt moving with velocity βt in a magnetic field B is F = qtβt × B. The magnetic force 5

2 ˆ ˆ q q (1 − ǫ ) βt · k rˆ − (βt · rˆ) k F = t s 3 2 2 λ0 (1 − β · n) r  (1h− 2ǫrˆ ·kˆ + ǫ ) i

ˆ  2 ˆ ˆ (rˆ · zˆ) βt × rˆ − βt × zˆ + ǫ k · rˆ βt × zˆ + ǫ [(βt · β) rˆ − (βt · rˆ) β] ǫ (βt · β) k − βt · k β − + , (12) (1 − 2ǫrˆ · kˆ + ǫ2)3/2 h(1 − 2ǫrˆ · kˆ + ǫ2)3/2 i 

where the right-hand side is evaluated at the retarded ǫ (βt · β) rˆ in the numerator of the second fraction in the time, except βt. It is apparent that the magnetic large square brackets of Eq. (12)), force will have a radially-directed term ( i.e., the term

qtqs ǫβrˆ qtqsrˆ β F r = − βt · 3 = − 2 βt · 3 , (13) λ0 β n rˆ kˆ 2 3/2 r β n rˆ kˆ 2 3/2 "(1 − · ) r(1 − 2ǫ · + ǫ ) #ret " (1 − · ) (1 − 2ǫ · + ǫ ) #ret

where in the second equality rˆ and r have been taken moving with the test particle, for interparticle separation outside the brackets indicating retardation, because the large compared to λ0, i.e., large compared to half a re- center of motion of the zitterbewegung field-source par- duced Compton wavelength, it is reasonable to disregard ticle is assumed to be stationary, and it has been used the change in position of the test particle as it moves in again that ǫ = λ0/r. For large enough interparticle sep- its zitter orbit, and consider only its velocity direction aration, ǫ << 1 and so time dependence. Then,

β β rˆ −qtqs t · s,ret βt·βs,ret = [− sin ωt, cos ωt, 0]·[− sin(ωt+δφ), cos(ωt+δφ), 0] F r ≈ , (14) (1 − β · n)3 r2 (15) where δφ is the phase difference of the test particle where in the numerator β is rewritten as βs,ret to em- present-time zitter motion to the retarded-time source phasize its association with the source particle and at particle zitter motion. Carrying out the multiplication the retarded time. The term given by Eq. (14), that obtains that contributes to the magnetic force caused by a stationary zitter particle on another zitter particle, differs from the Coulomb force caused by a stationary charge on another βt · βs,ret = cos δφ. (16) β β classical point charge by the factor − t · s,ret and the in- Thus, depending on the zitter phase difference between 3 3 verse factor of (1−β·n) . The inverse factor of (1−β·n) the particles, for aligned spins, the magnetic force may modulates the magnitude of the radial magnetic force at effectively either cancel or double the Coulomb force be- the frequency of the zitter motion, but averages to unity tween the particles. over a zitter cycle. In the case of the test particle circulatory motion op- To make a preliminary assessment of how the magnetic posite to that of the field source particle, force radial term given by Eq. (14) will affect the motion of a zitterbewegung test particle, the test particle may 2 2 be supposed to be executing a circular motion with the βt · βs,ret = −(sin ωt − cos ωt)cos δφ − sin 2ωt sin δφ, same radius λ0 as that of the field source particle, and (17) also at the speed of light so that βt = 1. (The inverse which averages to zero over a zitter cycle. However, (1 − β · n)3 factor will be disregarded here, supposing the magnetic force inverse square law radial component that its physical effects average out or that the test par- may be found to nonetheless have an influence on the ticle is located close enough to the z-axis so that β · n motion, when the average is more properly carried out is negligibly small.) For simplicity, suppose both parti- including the change in position of the test charge over cles have their zitter motion in the x-y plane and in the its zitter orbit. Also, there are other inverse-square law same direction. Also, while strictly rˆ is directed from magnetic force terms in Eq. (12) which have yet to be the source particle center of motion toward a field point considered. 6

For aligned zitter motions, it is straightforward to as- tum force, as shown in the appendix, does not undergo sess how the change in the source particle zitter phase such sign reversals in the ground state. induced by retardation might influence the test particle motion under the combination of electric and magnetic forces. The phase of the test particle zitter motion can be expressed as φt = φt0 + ωt, where φt0 ≡ φt(t = 0) is the V. TIME DILATION AND RELATIONSHIP TO test particle zitter phase at t = 0. Similarly, for aligned THE DE BROGLIE WAVELENGTH spins, the phase of the source particle zitter motion can be expressed as φs = φs0 + ωt, where φs0 is the source Up to this point the inverse-square-law magnetic force particle zitter phase at t = 0. Then the source particle between zitter particles has been calculated assuming phase at the retarded time is that the centers of zitter motion of the source and test particles are relatively stationary. Since it was noted that the magnetic force between stationary zitter par- φs(tr) ≡ φs,ret = φs0 + ω (t − rret/c)= φs(t) − rret/λ0, ticles varies sinusoidally with the separation between the (18) particles, it is naturally of interest to consider what is where rret ≡ r(tr). Thus, if δφ0 ≡ φt0 − φs0 is the the effect of relative motion between the zitter particles. phase difference between the zitter motions neglecting The effect on the relative zitter phase due to relative mo- retardation, then the phase difference δφret ≡ φt − φs,ret tion can be straightforwardly taken into account. The accounting for retardation is effect of relativistic time dilation on the zitter frequency of a moving zitter particle must also be taken into ac- count. For the present any physical consequences are yet rret(t) δφret(t)= δφ0 + . (19) to be worked out, but it seems worth noting that the λ 0 wavelength associated with the beat frequency between Therefore, over the course of an interparticle distance the rest frame and time dilated zitter frequencies can be change by an amount equal to π times the reduced Comp- equated with the de Broglie wavelength in the limit of ton wavelength 2λ0 of the particles, the magnetic force small relative velocity between the zitter particles. radial term of Eq. (13) will undergo a full cycle of si- Suppose that the test zitter particle is uniformly trans- ¯ nusoidal modulation. Also, because the radial magnetic lating with speed v = c|βt|, where the overbar indicates force can have a magnitude up to as large as the Coulomb an average over a test particle zitter period, in the iner- force between two static charges of the same magnitude, tial reference frame where the field-source zitter particle it can be expected to significantly influence the motion. is stationary. Then, by the Lorentz transformation, the However, it’s difficult to envision how the magnetic force phase of the test particle zitter motion as observed in derived here can be directly linked to the Bohmian quan- the source zitter particle rest frame can be written as tum force, since a change in separation of only one-half [12] ωτ = γω(t−υ ·r/c2), where τ is the time coordinate the Compton wavelength will invert the sign of the force. in the test zitter particle rest frame, γ = (1−(v/c)2)−1/2, Since the Bohr radius can be expressed as the Comp- and r here is the displacement from the source zitter par- ton wavelength divided by the fine-structure constant, ticle to the test zitter particle. The modulating factor on this corresponds to a change in separation of only about the inverse-square radial magnetic force in the numerator π/(2(137)) times the Bohr radius. The Bohmian quan- of the right-hand side of Eq. (14) becomes

2 2 βt · βs,ret = [− sin γω(t − υ · r/c ), cos γω(t − υ · r/c ), 0] · [− sin(ωtr + δφ), cos(ωtr + δφ0), 0], (20)

where δφ0 here is a constant to account for the phase difference between the source and test particle zitter mo- tions at (source particle rest frame) time t = 0. Carry- ing out the multiplication and applying the trigonometric βt · βs,ret = 2 identity sin α sin β + cos α cos β = cos(α − β) obtains cos (γ − 1)ωt − γωυ · r/c + ω(rret/c) − δφ0 . (22)  

2 Since the source zitter particle is stationary, rret(t) = βt · βs,ret = cos γω(t − υ · r/c ) − ωtr − δφ0 . (21) r(t) = r0 − vrt, where r0 ≡ r(t = 0) and vr = cβr is   the component of the test zitter particle average velocity Substituting for tr = t − rret/c, directly toward the source zitter particle, 7

where p ≡ γmev is the test zitter particle momentum. This differs from the de Broglie wavelength λde Broglie = −1 βt · βs,ret = h/p by the factor (1 + γ )/2, i.e., 2 cos (γ − 1)ωt − γωυ · r/c − ωβrt + φr , (23)

where φr ≡ ωr0/c − δφ0 is a constant.  For the case of circular motion of the test zitter particle around the source zitter particle, υ · r and βr are zero and so 1+ γ−1 vT = λ . (28) ret,nonrad 2 de Broglie

βt · βs,ret = cos((γ − 1)ωt + φr). (24)

2 Substituting for ω = c/λ0 = 2c me/~ obtains that 2 (γ − 1)ω =2c me(γ − 1)/~. The frequency νret,nonrad of the sinusoidal modulation of βt · βs,ret, for the retarded −1 field acting on a nonradially-moving test zitter particle Thus, in the low velocity limit where (1 + γ )/2 ≈ 1, is then the de Broglie wavelength can be equated with the spatial period of the beat due to the difference of the source and test particles’ zitter frequencies caused by time dilation. ω(γ − 1) 2c2m (γ − 1) Therefore, the de Broglie wavelength in the low-velocity ν = = e . (25) ret,nonrad 2π h limit can be obtained as a modulation of the magnetic force caused by the zitterbewegung of one Dirac particle The period Tret,nonrad of the sinusoidal modulation of acting upon another of equal mass. At large velocity, −1 βt · βs,ret is thus (1 + γ )/2 ≈ 1/2, and the modulation spatial period of the magnetic force acting on a (non-radially here) moving 1 h Dirac particle approaches half the de Broglie wavelength. Tret,nonrad = = 2 , (26) νret,nonrad 2c me(γ − 1) Therefore, in the low velocity limit, where (1 + and the distance traveled by the test zitter particle in γ−1)/2 ≈ 1, the de Broglie wavelength can be equated − one period of the modulation is (using γ 2 =1 − β2), with the spatial period of the beat due to the difference of the source and test particle zitter frequencies caused by time dilation. Over the course of travel of one de Broglie hv wavelength, for non-radial low-velocity relative motion, vTret,nonrad = 2 2c me(γ − 1) the magnetic force due to the zitter motions correspond- hv ing to aligned spins will undergo one cycle of sinusoidal = 2 −1 2c γme(1 − γ ) modulation. hv(1 + γ−1) = 2 −2 2c γme(1 − γ ) For βr nonzero, with the trigonometric identity cos(α+ h(1 + γ−1) β) = cos α cos β − sin α sin β, Eq. (23) can be rewritten = , (27) as 2p

2 2 βt · βs,ret = cos((γ − 1)ωt − γωυ · r/c ) cos(ωβrt − φr) − sin((γ − 1)ωt + γωυ · r/c ) sin(ωβrt − φr), (29)

illustrating that for non-zero radial motion it is not motion. Therefore modulation of the inverse-square mag- possible to factor βt · βs,ret such that one factor is a si- netic force due to radial relative motion is not obviously nusoidal modulation of frequency (γ − 1)ω, as needed to relatable to the de Broglie phase as is non-radial motion. obtain a modulation spatial period having the de Broglie This is contradicted by observation, where, for example, wavelength in the low-velocity limit. For small rela- the de Broglie phase is directly relatable to radial mo- tive velocity between the zitter particles and such that tion in the two-slit diffraction of electrons. However up 2 βr >> γ − 1 ≈ β /2, radial motion leads to much more to this point only retarded electromagnetic interactions rapid oscillation of the inverse square law magnetic force have been considered. between zitter particles than is the case for non-radial 8

VI. RADIAL RELATIVE MOTION IN THE where ta is the time advancement from the test zit- TIME-SYMMETRIC ELECTRODYNAMICS ter particle at time t to the source zitter particle at the PICTURE advanced time. Substituting for ta = t + radv/c,

Although the sign reversals of the inverse-square-law magnetic force between zitter particles, that result due βt · βs,adv = to small changes in interparticle separation compared 2 cos (γ − 1)ωt − γωυ · r/c − ω(radv/c) − δφ0 . (31) to atomic scale, seem to rule out an exact correspon- dence with the Bohmian quantum force for s states, up Since the source zitter particle is stationary, r  (t)= to this point only time-retarded electromagnetic interac- adv r(t)= r0−vrt. The modulating factor similar to βt·βs,ret tions have been considered. Since there is no physical of Eq. (23), for the case of a uniformly-translating test a priori basis for ruling out time-advanced interactions , zitter particle moving in the advanced field of the source it is worthwhile investigating what possible contribution zitter particle is thus they might make. For this purpose the time-symmetric inverse-square magnetic interaction is considered for the case of constant-velocity relative motion between the test β · β = and source zitter particles. t s,adv υ r 2 The modulating factor βt · βs for the case of a cos (γ − 1)ωt − γω · /c + ωβrt − φa , (32) uniformly-translating test zitter particle moving in the   advanced field of the source zitter particle is where φa ≡ ωr0/c + δφ0 is a constant. Combining this result with Eq. (23) for motion in the retarded field, the modulating factor for the test zitter particle moving in 2 βt · βs,adv = cos γω(t − υ · r/c ) − ωta − δφ0 , (30) the sum of the retarded and advanced fields is  

2 βt · βs,ret + βt · βs,adv = cos ω (γ − 1)t − γυ · r/c + φd · cos(ωβrt + φs) (33)   

where φs = (φr + φa)/2 and φd = (φr − φa)/2 are vanishing of the total time-symmetric Coulomb field of constants. a stationary charge, such is not the case for the electric Thus, summing the retarded and advanced magnetic acceleration fields of relativistically circulating charges. radial inverse square zitter force components, the total In the case of a vanishing time-symmetric electrostatic can be written as a product of two sinusoidal factors. force, the magnetic force considered here would need to This differs from the situation of the retarded-only force play the role of Coulomb attraction. as illustrated by Eq. (29), where such a factorization is not possible when there is radial relative motion. Fur- thermore, for general uniform relative motion, the first VII. DISCUSSION sinusoidal factor of the right hand side of Eq. (33) has an angular frequency of ω(γ − 1), which as shown above That the quantum force of Bohmian mechanics might leads to a modulation with a spatial period of the de be related to the magnetic force expected between Broglie wavelength in the limit of small relative velocity. charges undergoing relativistic circulatory motions on Therefore, for general motion, unlike in the retarded-only Compton wavelength scales seems consistent with recent case, the time-symmetric Coulomb-like magnetic force work showing quantum potential can always be between zitter particles is sinusoidally modulated with regarded as of concealed motion [13], and spatial period equal to the de Broglie wavelength times that the Bohmian quantum potential is an internal en- −1 (1 + γ )/2. ergy [14]. Also, the quantum potential has been identified Also, in the time-symmetric picture, there is an addi- directly with the zitterbewegung [15, 16]. If interacting tional electromagnetic force available that resembles the charges are moving relativistically, then magnetic forces Coulomb force between stationary charges. It is the elec- can be expected to be comparable in strength to Coulomb tric force due to the advanced electric field. Although interactions and should not be overlooked. time-symmetric electrodynamics usually assumes a time- The restriction to equal-mass particles herein has lim- symmetric field that is the mean of the retarded and ad- ited the applicabilty of the present analysis to hydro- vanced field, taking a difference of the advanced from genic atoms consisting of equal-mass particles, such as the retarded field would agree well with Bohm’s theory, the positronium atom. Consideration of the hydrogen at least in the s states. Although this would cause the atom itself will likely require regarding the proton prop- 9 erly as a composite of other spinning particles. The spin- ics picture where the total field acting is the mean of spin coupling of the electron to the proton is very weak the retarded and advanced fields [28–30]. In the case in the hydrogen atom because the proton intrinsic mag- of a stationary charge, the time-symmetric field as the netic moment is three orders of magnitude weaker than mean of retarded and advanced fields is identical to the the electron’s, resulting in only the hyperfine splitting. retarded-only or undelayed Coulomb field. Similarly, the Nonetheless, if all charged particles must relativistically time-symmetric field of a stationary current loop is iden- circulate on Compton wavelength scales, then magnetic tical to the retarded-only field. But, the time-symmetric couplings that are stronger than usually accounted for field of the point charge circulating around a fixed point is can be expected. non-constant, and distinct from the retarded-only field. Another problem of the present approach which must The time-symmetric picture is thus physically distinct be overcome is that the radial inverse-square magnetic from the retarded-only case, and so it may be possible to force tends to average to zero over a zitter period if the determine whether time-advanced interactions are nec- spins are opposite. Going forward it will be attempted essary to physical description of reality. Application of to obtain a deeper understanding by considering a more classical electrodynamics to explaining quantum behav- complete version of total inverse-square force, rather than ior while taking time-symmetric account of delay may be just the explicitly radial term. There are other inverse- found in [31–33]. The de Broglie wavelength in the low- square terms in Eq. (12) which, although not explicitly velocity limit is derived to order of magnitude in [32], radial, have radial contributions, and so need to be taken without direct insertion of the Planck constant as in the into account in analysis of the radial motion. Non-radial present approach. In [32], furthermore, as in the present forces need to be investigated as well. approach, non-instantaneous electromagnetic forces act Elementary particle internal periodicities and their time-symmetrically. modulations have recently been related to quantum be- havior [17–19]. Electromagnetic forces deriving from ac- celeration fields of highly relativistic charges can plausi- VIII. ERRATA bly mediate such periodicity modulations. Velocity fields may still play a role as well. For sub-luminal relativis- Previous versions of this article (e.g. [34]) accounted tic motion, electric and magnetic velocity fields remain properly for time dilation [35] on the test particle zitter finite but may be large in the plane of the zitterwegung frequency as ω′ = ω/γ, where ω is the zitter frequency in motion. For luminal motion, where the electromagnetic the test particle rest frame and γ = (1 − (v/c)2)−1/2 is field becomes singular, a field-free approach based on a the Lorentz factor. Although this is a correct description Fokker principal obtains a well-defined dynamical of the time dilation effect on frequency, it does not cor- picture [20]. rectly account for the phase of the zitter motion of the Since the early days of atomic physics, it has been rec- non-stationary test particle as observed in the laboratory ognized that the classical electrodynamics representation frame. The latter is properly [12] determined by Lorentz of the hydrogen atom as involving Kepler-like orbital mo- transformation as ωτ = γω(t − υ · x/c2), where τ is the tion with a correction for radiative decay is an approx- time coordinate in the test zitter particle rest frame, and imation. The exact electromagnetic two-body problem x is the test zitter particle position. The present version was considered intractable [21], however, and propaga- (v13) of this article is corrected accordingly beginning at tion delay and magnetic effects due to the electron in- Eq. (20) above. The corrections do not negate the orig- trinsic spin and magnetic moment, when they became inal conclusion of this article, that the time-symmetric known, were considered to be negligible or small correc- electrodynamic interaction between Dirac particles leads tions. Nonetheless, it has been proposed that accounting to a modulation of a Coulomb-like magnetic force that properly for propagation delay, which requires reformu- has similarity to the de Broglie matter wave. The correc- lating the electromagnetic two-body problem in terms tion makes the stated relation closer, in that it naturally of functional rather than ordinary differential equations, involves the relativistic momentum rather than its low- might account for quantum behavior as a phenomenon velocity approximation (as can be seen in Eq. (28) above, emergent from classical electrodynamics [22]. This math- and by comparing with the corresponding equation, Eq. ematically rigorous approach leads to “stiff” equations (26) of [34]. The correction to zitter phase as described of motion that can exhibit dynamical behavior on multi- will have wider benefits when applied to the sequel [36] to ple scales, where zitterbewegung-like relativistic motions the present article, which when corrected will be placed may be present simultaneously with non-relativistic av- as indicated at [37]. erage motions, and so may provide a classical electrody- namics explanation for the zitterwebegung motion that is a feature of the Dirac electron theory [23, 24] IX. CONCLUDING REMARKS It has been suggested that physicality of time-advanced interactions might account for the apparent non-locality A preliminary evaluation of the magnetic interaction of quantum behavior [25–27]. The magnetic force found between relativistically-circulating point charge particles here is also present in the time-symmetric electrodynam- has been performed. It was found that the magnetic in- 10 teraction between such particles has a radially-directed both its non-relativistic and relativistic [39] forms, the component with a strength equal to that of the Coulomb magnetic force must be derived from a vector potential. interaction, except that it is modulated by the rela- The quantum potential is [40] tive orientation of particle circulatory motions, and their propagation-delayed phase difference. Taking account of ~2 ∇2ρ 1 (∇ρ)2 relativistic time dilation due to motion of one of the parti- Q = − − , (A1) 4m ρ 2 ρ2 cles, for non-radial relative motion, it was found that the   spatial period of the modulation caused by the difference of zitterbewegung frequencies of the two particles can be where, if the Schr¨odinger wavefunction ψ ≡ − ~ 2 equated with the de Broglie wavelength times (1+γ 1)/2. R exp(iS/ ), then ρ ≡ R . For radial relative motion the de Broglie wavelength re- For the hydrogen ground state, the normalized lated modulation was found only if the electromagnetic Schr¨odinger wavefunction is [41] field was time symmetric, i.e., is the sum of retarded and advanced fields. This type of interaction might find 2 − application in efforts to explain quantum behavior using ψ = e r/a0 , (A2) a 3/2 only classical physics.  0  or, equivalently, Appendix A: Evaluation of the Quantum Force in the Hydrogen Atom Ground State 4 −2r/a0 ρ = 3 e , (A3) a0 The quantum potential and resulting quantum force is   calculated from the Schr¨odinger wavefunction. In the ap- and the quantum potential evaluates as plication of the Schr¨odinger equation to the hydrogen-like atom, the classical potential term is simply the Coulomb ~2 2 1 potential of a point charge, while the quantum potential Q = − , (A4) 2ma r a depends on the solution. For the ground state solution 0  0  of the time-independent Schr¨odinger equation, and for all s-states, the electron momentum vanishes [38]. In these and the quantum force is cases, the quantum potential is everywhere opposite the Coulomb potential, and so the quantum force exactly ~2 rˆ cancels the Coulomb force. However, the quantum force FQ = −∇Q = . (A5) ma r2 is not usually evaluated explicitly in the literature, since 0   it is not needed to obtain the physical results of inter- 2 2 With a0 = ~ /me , est (e.g., that the electron momentum vanishes), so for completeness it is evaluated in the following for the hy- drogen atom ground state. Also, it serves to illustrate rˆ F = e2 . (A6) that although the magnetic force component described Q r2 here bears a similarity to the quantum force, it cannot   be equated to it directly. More generally, the quantum The quantum force thus exactly cancels the Coulomb force cannot in principle be equatable to the magnetic attraction, in the Schr¨odinger hydrogen atom ground force, because while the quantum potential is scalar in state.

[1] D. Bohm, A Suggested Interpretation of the Quantum [7] M. Rivas, Kinematical Theory of Spinning Particles, Theory in Terms of “Hidden” Variables. I, Phys. Rev. Kluwer Academic, Dordrecht (2001) 85(2), 166-179 (1952) [8] M. Rivas, An interaction Lagrangian for two spin 1/2 [2] D. Hestenes, The Zitterbewegung Interpretation of elementary Dirac particles, J.Phys.A 40:2541-2552 (2007) , Found. Phys. 20(10), 1213-1232 [9] J. D. Jackson, Classical Electrodynamics, Wiley, New (1990) York (1998) [3] P.A.M. Dirac, The Quantum Theory of the Electron, [10] D. Hestenes, Zitterbewegung in Quantum Mechanics, Proc. Roy. Soc. Lon. A 117, 610 (1928) Found. Physics. 40(1), 1-54 (2010) [4] E. Schr¨odinger, Sitzungber. Preuss. Akad. Wiss. Phys.- [11] Rivas, op. cit., Section 6.1.1, 254-264 Math. Kl. 24, 418 (1930) [12] W. Baylis, De Broglie waves as an effect of clock desyn- [5] K. Huang, On the Zitterbewegung of the Dirac Electron, chronization, Can. Jour. Phys. DOI: 10.1139/P07-121 Am. J. Phys. 20 (479) (1952) (Dec 2007) [6] D. Hestenes, op. cit. [13] P. Holland, Quantum potential energy as concealed mo- 11

tion, Found. Phys. 45(2) 134-141, (2015) [28] P. A. M. Dirac, Classical theory of radiating electrons, [14] G. Dennis, M. de Gosson, B. Hiley, Bohm’s Quan- Proc Roy Soc A, 148 (1938) tum Potential as an Internal Energy, (preprint) [29] J. A. Wheeler and R.P. Feynman, Interaction with the arxiv.org/abs/arXiv:1412.5133 (2014) Absorber as the Mechanism of Radiation, Rev. Mod. [15] G. Salesi, Spin and the Madelung Fluid, Mod. Phys. Lett. Phys. 17, 157-181 (1945) A 11, 1815 (1996) [30] A. Schild, Electromagnetic Two-Body Problem, Phys. [16] G. Salesi, E. Recami A Velocity Field and Operator for Rev. 131, 2762 (1963) Spinning Particles in (Nonrelativistic) Quantum Mechan- [31] Jayme De Luca, Variational principle for the Wheeler- ics, Found. Phys. 28: 763-773 (2010) Feynman electrodynamics, Jour. Math. Phys. 50, 062701 [17] Donatello Dolce, Compact Time and Determinism for (2009) Bosons, Found. Phys. 41:178, 178-203 (2011) [32] Jayme De Luca, Double-Slit and Electromagnetic Mod- [18] Donatello Dolce, Gauge Interaction as Periodicity Mod- els to Complete Quantum Mechanics, Jour. Comp. and ulation, Annals of Physics, 327(6), 1562-1592 (2012) Theor. Nanoscience, 8(6), 1040-1051 (2010) [19] D. Dolce, A. Perali, The role of quantum recurrence in [33] Jayme De Luca, Variational electrodynamics of Atoms, superconductivity, carbon nanotubes, and related gauge Progress In Electromagnetics Research B, vol. 53, 147- symmetry breaking, Found. Phys. 44(9), 905-922 (2014) 186 (2013) [20] Jayme De Luca, Minimizers with discontinuous velocities [34] D. Lush, Does Bohm’s Quantum Force Have for the electromagnetic variational method, Phys. Rev. E, a Classical Origin?, Found Phys (2016) 46: 82, 026212 (2010) 1006. doi:10.1007/s10701-016-9990-1 (preprint) [21] Rodney D. Driver, A “Backwards” Two-Body Problem of http://arxiv.org/abs/1409.8271v10 (2016) Classical Relativistic Electrodynamics, Phys. Rev. 178(5) [35] Jackson, op. cit., Equation (11.26), page 528 (1969) [36] D. Lush, Spin-Spin Interaction as Quantum Phase [22] C. K. Raju, The Electrodynamic 2-Body Problem and Measurement and the Mechanism of Atom Formation the Origin of Quantum Mechanics, Found. Phys. 34(6), (preprint) http://arxiv.org/abs/1609.04446v7 (2017) 937-962 (2004) [37] D. Lush, Spin-Spin Interaction as Quantum Phase [23] Jayme De Luca, Stiff three-frequency orbit of the hydro- Measurement and the Mechanism of Atom Formation gen atom, Phys. Rev. E, 026221 (2006) (preprint) http://arxiv.org/abs/1609.04446 (2019) [24] C. K. Raju, Suvrat Raju, Radiative Damping and Func- [38] B. J. Hiley, Is the electron stationary in the ground tional Differential Equations, Mod. Phys. Lett. A, 26(35), state of the Dirac hydrogen atom in Bohms Theory?, 2627-2638 (2011) (preprint) arxiv.org/abs/ arXiv:1412.5887v1 (2014) [25] Huw Price, Time’s Arrow and Archimedes Point: New [39] B. J. Hiley, R. E. Callaghan, Clifford Algebras and Directions for the Physics of Time, Oxford (1996) the Dirac-Bohm Quantum Hamilton-Jacobi Equation, [26] Huw Price, Ken Wharton, A Live Alternative to Quan- Found. Phys. 42 192-208 (2012) tum Spooks, (preprint) arxiv.org/abs/arXiv:1510.06712 [40] D. Bohm, op. cit. (2015) [41] David J. Griffiths, Introduction to Elementary Particles, [27] Huw Price, Ken Wharton, Disentangling the Quantum 162, Wiley-VCH, Weinheim (2008) World, (preprint) arxiv.org/abs/arXiv:1508.01140 (2015)