Zitterbewegung Revisited: Insights into Spacetime

Zitterbewegung Revisited: Insights into Spacetime Damian Pikor, Paweł Kurzawski Abstract Standard physics treats Zitterbewegung - trembling motion predicted by the Dirac theory - as a mathematical artifact. This is because it suggests an internal structure for the electron, whereas experiments consistently show the electron to be a point particle. In this paper, we propose a revision of this view. We postulate that Zitterbewegung is a real phenomenon, not as a motion of the particle’s charge, but as a persistent vacuum resonance; an imprint of its creation. We argue this dynamic constitutes a stable, internal phase geometry associated with the particle. The central point of the model is the strict separation of the point-like particle from this associated geometry. We formalize this by introducing two distinct form factors. The charge form factor ( FEM ) describes the particle itself and remains operationally point-like ( FEM ≡ 1), ensuring full consistency with scattering data. Simultaneously, we define a new, internal observable – the phase form factor ( Fϕ ) – which describes the coherence of this associated phase geometry at the Compton scale ( ¯ λC ). This model is calibrated against the anomalous magnetic moment of the muon (∆ aµ ). This calibration then predicts the observed "silence" of the electron anomaly (∆ ae ) and a significant, testable anomaly for the tau lepton (∆ aτ ). This approach resolves the fundamental paradox, linking the g− 2anomaly to the Dirac equation’s kinematic substructure. The model predicts that this phase structure ( Fϕ ), while invisible to scattering, is accessible via coherence-sensitive protocols, such as interferometry. Key words: Zitterbewegung, Lepton Geometry, Dirac Equation, Compton Scale, Anomalous Magnetic Moment, g-2, SMEFT. 1 1. Introduction The relativistic Dirac equation is a cornerstone of lepton dynamics, successfully describing spin and predicting antimatter [1-3]. However, it also yields one of physics’ most enduring puzzles: Zitterbewegung (ZBW), an ultra-rapid trembling of the free particle’s position operator. This motion is characterized by the ZBW frequency ωZ = 2 mc2/ℏ and an amplitude on the order of the reduced Compton wavelength, ¯ λC = ℏ/ ( mc ). This prediction is in direct conflict with observation. High-energy scattering experiments and precision QED tests consistently constrain the lepton’s charge radius to be much smaller than its Compton wavelength, ≪¯ λC [4, 5]. Because ZBW implies an internal structure at the ¯ λC scale, it is traditionally dismissed as an unphysical mathematical artifact [1, 3, 6]. In this work, we revise this view. We postulate that the ZBW paradox arises from a misinterpretation of the phenomenon. Our central hypothesis is that ZBW is a real, physical process, but it is not the motion of the particle’s charge. Instead, we posit that ZBW is a persistent, dynamic resonance of the vacuum itself, imprinted in the particle’s vicinity during its creation (e.g., pair production). The particle remains point-like but is inextricably "dressed" by this associated, coherent phase-geometry. The key to our model is a strict sector separation hypothesis. We formalize this by defining two distinct form factors for the lepton: 1. The Charge Form Factor ( FEM ( q2 )): This is the standard observable probed in scattering experiments. In our model, the charge density is a point-like delta function, ρEM (x) = δ(3) (x). Consequently, FEM ( q2 ) ≡ 1. This "geometric silence" in the charge sector ensures full consistency with scattering data and precision QED tests [4, 5]. 2. The Phase Form Factor ( Fϕ ( q2 )): This is a new, internal observable we introduce. It does not describe charge, but rather the spatial coherence of the associated ZBW vacuum resonance. Unlike FEM , this phase-form factor possesses a non-trivial structure, Fϕ(q2)= 1, with a characteristic scale defined by ¯ λC. This model resolves the paradox: FEM ≡ 1explains why scattering experiments see a point particle, while Fϕ = 1 allows ZBW to be a real phenomenon with its characteristic Compton scale. Most importantly, this framework is not just a conceptual model; it has concrete, falsifiable consequences. We propose that this underlying phase-geometry is the source of the anomalous magnetic moment (∆ aℓ ) of the leptons. We use the persistent, high-precision discrepancy in the muon anomalous magnetic moment, ∆ aµ , as a calibration input for our model. By anchoring our geometry to the observed ∆ aµ , the model becomes predictive. It correctly "post-dicts" the observed null result for the electron anomaly (∆ae≈0) and predicts a specific, large, and potentially measurable anomaly for the tau lepton, ∆aτ. This article is organized as follows: Section 2 details the formal, covariant theory of the phase-geometric world-surface and its observables. Section 3 presents the results of our numerical analysis, showing how the model is calibrated on ∆ aµ and its predictions for ∆ ae and ∆ aτ across the UV and Planck-scale parameter space. Section 4 synthesizes these findings, discussing the model’s implications, its connection to the r.m.s. phase radius r(ϕ) rms , and its falsifiability via interferometry. Detailed derivations and numerical protocols are provided in the Appendices. 2 2. Theoretical Framework 2.1. Assumptions and Sector Separation We introduce a strict separation of observables into a charge sector and a phase-geometric sector. In the rest frame, the charge density is point-like, ρEM (x) = δ(3) (x), yielding FEM ( q ) ≡ 1[7, 8]. This ensures operational point-likeness in scattering. The new sector possesses a radial coherence density ρϕ ( r )with a natural scale ¯ λC and an independent phase form factor Fϕ(q2). We define the phase form factor in the rest frame via the spherical transform Fϕ(q) = 4πZ∞ 0 dr r2ρϕ(r)sin(qr) qr and treat it relativistically as a function solely of the invariant q2 = qµqµ , which ensures manifest covariance. The r.m.s. phase radius is determined by the standard slope-radius relation ⟨r2⟩ϕ = − 6 d Fϕ/ d q2|q2=0 , with r(ϕ) rms = p⟨r2⟩ϕ . This formally maps the form-factor apparatus onto coherence observables. 2.2. Covariant Resonance World-Surface (Σ) The coherence carrier is defined as a smooth, two-dimensional resonance world-surface Σ ⊂R1,3 with an embedding Xµ ( σa ), induced metric γab = ηµν∂aXµ∂bXν , two normals nµ I , a second fundamental form KI ab , curvatures H and KG , and the Laplace-Beltrami operator ∆Σ. All observables are built exclusively from Poincaré scalars. The covariant phase density in spacetime is defined distributionally as ϱϕ(x;u) = NZΣ d2σ√−γ δ(4)(x−X(σ)). This is projected onto the hyperplane orthogonal to the four-velocity uµvia ρϕ(x) = Zdτ ϱϕ(xµ=τuµ+xµ ⊥;u), which guarantees a frame-independent definition of the radial density ρϕ ( r ). A sketch of the proof of covariance is provided in Appendix A. 2.3. Action Functional on Σ The minimal, manifestly covariant action functional combines surface tension, curvature elasticity, and a phase field ϑ on Σ, with a "soft" resonance condition at the Compton scale ℓ∼¯ λC . The full action S [ X, ϑ ]and the derivation of the Helfrich/Willmore-type shape equation [9, 10], its boundary conditions, and stability analysis (positive second variation) are detailed in Appendix B. 2.4. Projection and Definition of ρϕ(r) The projection from ϱϕ ( x ; u )onto the rest frame using the projector Pµν = ηµν−uµuν leads to an isotropic radial density ρϕ ( r )satisfying the normalization 4 πR∞ 0 d r r2ρϕ ( r ) = 1. 3 This defines the moments ⟨rn⟩ϕ and the spherical transform Fϕ . The procedure is reparametrization-invariant, and UV/IR regulators on Σensure a unique projection and the existence of moments, as is standard in analyses of surfaces with bending energy. This provides a direct interface to interferometric measurements via the Fϕ(q2)slope. 2.5. Geometric Ansatz: "Gabriel’s Horn" For an axially-symmetric surface of revolution with profile r ( x ) = a/x , the mean curvature H asymptotically realizes the resonance condition 2 ℓH ≃ 1and selects the scale r∼ℓ . A small-slope approximation H≃1 2 ( −r′′ +1 /r )with a self-similar ansatz r ( x ) = ax−α yields an energy minimum at α≃ 1, justifying the inverse profile. After non-dimensionalization ( ˆr = r/ℓ, ˆx = x/ℓ ), this anchors the geometry to ¯ λC . The profile requires UV/IR regulators [ rmin, rmax ], and the pushforward of the surface measure leads to a normalizable density. The full derivation of the normalization constant C−1is provided in Appendix C. 2.6. Spectral Ansatz: Bessel + Localization Complementarily, we employ a spectral ansatz ρϕ ( r ) = A r J2 1 ( kr ) W ( r ; rmin, rmax )with a C∞ window function W , where k = c′/ℓ sets the Compton scale. The normalization constant is A−1 = 4 πR d r r3J2 1 ( kr ) W ( r ). The J1 ( z ) ∼z/ 2behavior at r→ 0ensures UV regularity, while the asymptotic behavior necessitates the IR-damping window W . The calculation of Fϕ ( q2 )is realized as a spherical Bessel transform j0 ( qr ) = sin ( qr ) / ( qr ) using methods detailed in Appendix D. 2.7. Low-q2Expansion and "Shape-Blindness" The expansion of j0(z) = 1 −z2/6+... yields the low-q2expansion for the form factor: Fϕ(q2)=1−q2 6⟨r2⟩ϕ+(q2)2 120 ⟨r4⟩ϕ+O((q2)3). This defines the slope-radius relation dFϕ/dq2|0 = −⟨r2⟩ϕ/ 6. The "shape-blindness" property implies that at O ( q2 ), Fϕ depends only on the second moment ⟨r2⟩ϕ . Any two normalized ansätze (e.g., Geometric and Spectral) with the same ⟨r2⟩ϕ will have an identical slope and r.m.s. radius, with differences only appearing at O (( q2 ) 2 )via the ⟨r4⟩ϕ moment, IR anchoring with ⟨r2⟩ϕ = 1 . 5(Compton units) and the demonstration of shape-blindness (agreement in the slope O ( q2 )between ansätze after matching ⟨r2⟩ϕ ) are illustrated through data and plots in Appendix G. 2.8. Parameter Estimation and Anchoring The ansatz parameters (e.g., a, k ) are determined by a unified optimization problem. The procedure minimizes the geometric action Sgeo (see Appendix B) subject to the constraints of normalization (4 πR d r r2ρϕ ( r )=1) and anchoring to the Compton scale ℓ = ¯ λC . A critical consistency test ("slope-match") ensures that both ansätze yield the same ⟨r2⟩ϕ and thus the same low-q2physics. 4 2.9. Ansätze Comparison: IR Equivalence and UV Divergence As established in Sec. 2.7, the Geometric and Spectral ansätze are equivalent in the lowq2 (IR) regime, provided their second moments ⟨r2⟩ϕ are matched. Their divergence occurs at highq2 (UV), which is dictated by the larger tails of their respective density profiles and the choice of regularization window W ( r ). This ensures that while the lowq2 observables ( r(ϕ) rms ,∆ aℓ ) are robust, the highq2 behavior remains distinct. A numerical verification of this IR equivalence is provided in Appendix C. 3. Numerical Analysis and Phenomenological Projections 3.1. UV Grid and Planck-Line Scans We performed a numerical analysis based on the framework detailed in Section 2 and Appendix D. We scanned the parameter space in two primary ways: 1. UV Grid Scan (54 nodes): We performed a Cartesian scan over the regulator space (Λ , p ). This scan confirms a robust separation of roles. The muon anchor ∆ aµ = 2 . 5 × 10 −9 is held constant by construction, and the electron anomaly remains negligible, with a mean ⟨ ∆ ae⟩ ≈ 1 . 084 × 10 −17 (std. dev. 1 . 93 × 10 −19 ). In contrast, the tau anomaly exhibits a wide, stable plateau, with a mean ⟨ ∆ aτ⟩ ≈ 2 . 463 × 10 −7 (std. dev. 4.40 ×10−9), as shown in smeft_uv_summary.csv. 2. Planck-Line Trajectory (Λ(ℓ) = c/ℓ): We probed the scale-dependence by sampling a trajectory along the "Planck-Line." This scan ( smeft_planck_summary.csv ) confirms that the light-lepton anchors (∆ ae, ∆ aµ ) are invariant. The tau channel, however, exhibits the characteristic geometry-controlled crossover, interpolating from the classical plateau (max. ∆ aτ≈ 2 . 867 × 10 −4 ) to a deep compression window (min. ∆aτ≈5.842 ×10−9). These scans numerically validate the "IR-anchored, UV-sensitive" nature of the model. The heavy-lepton channel ( τ ) acts as a short-distance probe, while the light-lepton channels (e, µ) are protected by the low-q2"shape-blindness" of the geometry. 3.2. SMEFT Projection and Consistency We project these dipole predictions into the Warsaw basis of the Standard Model Effective Field Theory (SMEFT) [11]. We map the anomalies ∆ aℓ to the electromagnetic dipole operator coefficient Cℓeγ and its SU (2) L×U (1) Y decomposition ( CℓeB, CℓeW ). The full conventions are detailed in Appendix E. This mapping confirms the results of our scans at the operator level. • UV Grid: Cµeγ is constant with zero variance (mean 9 . 610 × 10 −11 , std. dev. 3 . 88 × 10 −26 ). Ceeγ is ultra-small (mean 1 . 084 × 10 −17 , std. dev. 1 . 93 × 10 −19 ). Cτeγ shows a stable plateau (mean 2.463 ×10−7, std. dev. 4.40 ×10−9). • Planck-Line: Cµeγ and Ceeγ remain fixed, while Cτeγ spans orders of magnitude, from a minimum of 5.842 ×10−9to a maximum of 2.867 ×10−4. 5 Critically, we verify the lowq2 consistency at every node by calculating the reconstruction residuals ( εℓ ). As shown in Appendix D (and smeft_planck_summary.csv ), these residuals are zero (e.g., dat_resid mean 1 . 07 × 10 −18 , std. dev. 3 . 70 × 10 −18 ) for all leptons across all scans, certifying that the low-q2slope-radius relation is exactly preserved. 3.3. Neutrino and Z′Portal Consistency We performed two unit tests to ensure our framework is consistent with other phenomenological constraints (details in Appendix E). 1. Neutrino Unit Test: We defined an effective neutrino magnetic moment proxy, ξν , and determined the maximum allowed value consistent with experimental limits. Across the 54-node UV grid ( smeft_nu_summary.csv ), this value is uniformly small and stable: ξmax ν has a mean of 6 . 940 × 10 −11 with a standard deviation of only 6 . 347 × 10 −13 (min 6 . 531 × 10 −11 , max 6 . 959 × 10 −11 ). This confirms that our model induces no tension with neutrino constraints, demonstrating robust sector separation. 2. Leptophilic Z′ Portal: We confirmed that a modular Z′ portal can be introduced to account for the g− 2anomaly. A scan over ( mZ′, g′ )shows viable parameter space (e.g., mZ′≃ 0 . 261 GeV , g′≃ 1 . 78 × 10 −3 ) that matches the ∆ aµ target while passing a proxy "trident" screen [12, 13]. This demonstrates that the Z′ can be treated as an orthogonal effector, leaving the core ZBW-driven ∆ ae/ ∆ aτ predictions intact. 4. Discussion and Conclusion We have proposed a revision of the Zitterbewegung paradox, postulating that ZBW is a real vacuum resonance that endows leptons with a "phase-geometry" at the Compton scale, ¯ λC . This geometry is distinct from the point-like charge structure ( FEM ≡ 1) and is described by a new observable, the phase form factor Fϕ(q2). We demonstrated that this phase-geometry, when calibrated to the observed muon g− 2, provides a coherent framework for lepton anomalies. The model’s "shape-blindness" at lowq2 naturally explains the IR-anchoring of the light leptons, resulting in a negligible electron anomaly (mean ∆ ae≈ 1 . 084 × 10 −17 ) while remaining consistent with the muon anchor (∆aµ= 2.5×10−9). The model’s primary predictive power lies in the UV-sensitive tau lepton. The framework predicts a significant tau anomaly (∆ aτ ) that varies from a "classical plateau" of ∆ aτ≈ 2 . 463 × 10 −7 (on the UV grid) down to a "compressed" value of ∆ aτ≈ 5 . 842 × 10 −9 (on the Planck-line). This provides two distinct, falsifiable experimental regimes. Furthermore, the lowq2 slope of the phase form factor defines a measurable r.m.s. phase radius, r(ϕ) rms . Our g− 2calibration is consistent with the initial prediction from the simplest ZBW models, yielding r(ϕ) rms ≈p3/2¯ λC . This prediction is not accessible via scattering but is testable in coherence-sensitive experiments, such as electron interferometry. A concrete experimental protocol for this measurement is proposed in Appendix F, At small q , we recover rrms ( ϕ ) = p3/2¯ λC from the slope of the visibility; details of V ( q )vs q2 6 fits, comparisons with the curve 1 −q2 6⟨r2⟩ϕ , and residual tables for the UV-grid and Planck-line are found in Appendix G. This work provides a framework that resolves a foundational paradox of the Dirac equation and links it directly to modern precision phenomenology. It recasts the g− 2anomaly not as a sign of new particles, but as the first measurement of the particle’s inherent phase-geometric structure. The model is numerically robust (as shown by zero-residuals in Appendix D), consistent with SMEFT mapping and other constraints, and offers a clear, falsifiable experimental program for both precision ∆ aτ measurements [14, 15] and interferometric r(ϕ) rms extraction. 7 A. Proof Sketch: Covariance of the Construction The proof that all observables, particularly Fϕ ( q2 ), are Lorentz invariants follows from the fact that they are constructed entirely from Poincaré scalars. Given a Lorentz transformation Λ ∈SO+ (1 , 3), the objects transform as: x′µ = Λ µνxν , u′µ = Λ µνuν , and X′µ= ΛµνXν. 1. Invariance of Measures: The metric ηµν is invariant by definition (Λ Tη Λ = η ). This guarantees that all contractions (like uµuµ = 1) and the induced metric γ′ ab = ηµν∂aX′µ∂bX′ν = γab are invariant. Consequently, the surface element √−γ′ d 2σ = √−γd2σis a scalar. 2. Invariance of Distributions: The 4D Dirac delta transforms as a scalar density: δ(4)(x′−X′)=δ(4)(Λ(x−X))=δ(4)(x−X)due to |det Λ|= 1. 3. Invariance of Projection: The projector P′µν = ηµν−u′µu′ ν transforms covariantly. The squared radial distance r′2=x′µ ⊥x′ ⊥µis invariant. Therefore, the covariant density ϱ′ ϕ ( x′ ; u′ ) = ϱϕ ( x ; u )retains its numerical value. Since the projection onto the rest frame is defined entirely by these invariant objects, the resulting radial density ρϕ ( r )and its moments ⟨rn⟩ϕ are Lorentz scalars. The form factor Fϕ , being a function of ρϕ ( r )and q , must therefore be a function of the only available Poincaré invariant, q2=qµqµ, in a fully relativistic treatment. B. Action Functional and Variational Principle B.1. Action Functional As stated in Section 2, the minimal, covariant action functional on Σis given by: S[X, ϑ] = ZΣ d2σ√−γT+κ 2(2H−C0)2+¯κ 2KG+h 2γab∂aϑ ∂bϑ+µ(2ℓH −1)2 where T is tension, κ is bending modulus, ¯κ is the Gaussian curvature modulus, C0 is spontaneous curvature, h is phase field rigidity, and µ is the resonance weight for the scale ℓ∼¯ λC. B.2. Variational Equations and Boundary Conditions Varying S with respect to ϑ yields the Laplace-Beltrami equation on Σ: ∇a ( √−γγab∂bϑ ) = 0. Varying S with respect to the embedding Xµ yields the generalized Helfrich/Willmore shape equation, which includes pressure from the ϑfield and the resonance term: 2TH−2κ(∆ΣH+ 2H(H2−KG))+hPφ+ 4µℓ(2ℓH −1)H= 0 (+ boundary terms) where Pφ is the effective pressure from the ϑ field. The stability of the solution (e.g., the "Gabriel’s Horn" ansatz) is ensured by analyzing the second variation, δ2S . The resonance term µ (2 ℓH − 1) 2 generates a positive "mass" term ∝µℓ2 , which ensures δ2S > 0and stabilizes the solution against normal fluctuations for a realistic parameter window (T, κ, µ > 0). 8 C. Ansatz Derivations and Low-q2Equivalence C.1. Gabriel’s Horn Normalization For the geometric ansatz r ( x ) = a/x , the unnormalized density (derived from the pushforward of the surface measure dΣ=2πr√1+r′2dx) is: ˜ρϕ(r)∝dΣ |dr|=2πr√1+r′2 |r′|=2πrp1+(r4/a2) r2/a = 2πap1+r4/a2 r We define the normalized density ρϕ ( r ) = C·˜ρϕ ( r )on the compact support [ rmin, rmax ]. The normalization constant Cis fixed by 4πRrmax rmin dr r2ρϕ(r)=1, which gives: 1 = 8π2aC Zrmax rmin dr rp1+r4/a2 Using the substitution u = r2 , this integral becomes 4 π2CRumax umin d u√a2+u2 . Solving this yields the exact normalization constant C−1: C−1= 2π2hu√u2+a2+a2ln(u+√u2+a2)iu=r2 max u=r2 min This can be written as C−1 = 2 π2 [ F ( rmax ; a ) −F ( rmin ; a )], where F ( r ; a ) = r2√r4+a2 + a2ln(r2+√r4+a2). C.2. Numerical Test of Shape-Blindness To confirm the lowq2 equivalence of the Geometric and Spectral ansätze (the "shapeblindness" principle), we numerically compute both form factors with their second moments ⟨r2⟩ϕ matched. As shown in Figure 1, while the form factors diverge at high q2 , they are identical at q2→0, possessing the exact same value and slope. D. Numerical and Reproducibility Specifications D.1. Numerical Algorithms and Regulators The Fϕ ( q2 )integrals are oscillatory and were computed using adaptive quadrature for the Geometric ansatz. For the Spectral ( J2 1 ) ansatz, specialized fast Hankel transform algorithms (QDHT/FFTLog) were used for efficiency. The numerical scans were performed with the following specifications: •UV Grid: 54-node Cartesian scan over (Λ, p). •Planck-Line: Trajectory scan Λ(ℓ) = c/ℓ. • Regulators: All ansätze were regularized on a compact support [ rmin, rmax ]with a C∞(smooth) window function W(r). • Consistency: Convergence was verified by ensuring Fϕ (0) ≈ 1and by checking the lowq2 slope-match. The residuals (e.g., in smeft_planck_summary.csv ) are ≈ 0(mean and std. dev.), confirming the "shape-blindness" principle to machine precision. 9 [13] C. Fox et al., "Dressing Lµ−Lτin Color," Phys. Rev. D 89, 095033 (2014). [14] A. Crivellin et al., "LHC tau-pair production constraints on aτ and dτ ," SciPost Phys. 16, 048 (2023). [15] H. Gay, N. Schulman et al., "Strategy to measure tau g-2 via photon fusion in LHC proton collisions," Phys. Rev. D 110, 092016 (2024). 16