# Simulation of the Radiation Reaction Orbits of a Classical Relativistic Charged Particle with

Authors: Aviad Roitgrund; Lawrence Horwitz
Journal: Discrete Dynamics in Nature and Society
 (2010)

---

## Abstract

We review the formulation of the problem of electromagnetic self-interaction of a relativistic charged particle in the framework of the manifestly covariant classical mechanics of Stueckeleberg, Horwitz, and Piron. The gauge fields of this theory, in general, cause the mass of the particle to change. We study the four dynamical off-mass-shell orbit equations which result from the expansion of Green's function in the Lorentz force equation for the self-interaction. It appears that there is an attractor in this system which stabilizes the motion of the
relativistic charged electron. The attractor may acquire fractal characteristics in the presence of an external field and thus become a strange
attractor.


---

## Body

## 1. Introduction

This work is concerned with studying the problem of the classical relativistic charged particle. The existence of radiation due to the accelerated motion of the particle raises the question of how this radiation, acting back on the particle, affects its motion. Rohrlich has described the historical development of this problem, where the first steps were taken by Abraham in 1905, culminating in the work of Dirac, who derived the equation for the ideal point electron in the form(1.1)md2xμds2=ecFνμdxνds+Γμ,
where m is the electron mass, including electromagnetic correction, s is the proper time along the trajectory xμ(s) in spacetime, Fνμ is the covariant form of the electromagnetic force tensor, e is the electron charge, and(1.2)Γμ=23e2c3(d3xμds3-d2xνds2d2xνds2dxμds).Here, the indicesμ, ν, running over 0,1,2,3 label the spacetime variables that represent the action of the Lorentz group; the index raising and lowering Lorentz invariant tensor ημν is of the form diag(-1,+1,+1,+1). The expression for Γμ was originally found by Abraham in 1905 [1], shortly after the discovery of special relativity and is known as the Abraham four-vector of radiation reaction. Dirac's derivation [2] was based on a direct application of Green's functions for the Maxwell fields, obtaining the form (1.1), the Abraham-Lorentz-Dirac equation. In this calculation, Dirac used the difference between retarded and advanced Green's functions so as to eliminate the singularity carried by each. Sokolov and Ternov [3], for example, give a derivation using the retarded Green's function alone and show how the singular term can be absorbed into the mass m.The formula (1.1) contains the so-called singular perturbation problem. There is a small coefficient multiplying a derivative of higher order than that of the unperturbed problem; since the highest-order derivative is the most important in the equation, one sees that dividing by e2, any small deviation in the lower order terms results in a large effect on the orbit, and there are unstable solutions, often called “runaway solutions”. There is extensive literature (see [4–6] and the references therein) on the methods of treating this instability.The existence of these “runaway solutions”, for which the electron undergoes an exponential acceleration with no external force beyond a short initial perturbation of the free motion, is a difficulty for the point electron picture in the framework of the Maxwell theory with the covariant Lorentz force. Rohrlich [7–9] has discussed the idea that the point electron idealization may not be really physical, based on arguments from classical and quantum theory and emphasized that for the corresponding classical problem, a finite size can eliminate this instability. Of course, this argument is valid, but it leaves open the question of the consistency of the Maxwell-Lorentz theory which admits the concept of point charges.Stueckelberg, in 1941 [10–14], proposed a manifestly covariant form of classical and quantum mechanics in which space and time become dynamical observables. They are therefore represented in quantum theory by operators on a Hilbert space on square integrable functions in space and time. The dynamical development of the state is controlled by an invariant parameter τ, which one might call the world time, coinciding with the time on the (on mass shell) freely falling clocks of general relativity.Stueckelberg postulated the existence of an invariant “Hamiltonian”K, for the classical theory, which would generate Hamilton equations for the canonical variables xμ and pμ of the form (we take units for which c=1 unless otherwise specified)(1.3)ẋμ=∂K∂pμ,ṗμ=-∂K∂xμ,
where the dot indicates differentiation with respect to τ. Taking, for example, the model(1.4)K0=pμpμ2M,
we see that the Hamilton equations imply that(1.5)ẋμ=pμM.
It then follows that(1.6)dx⃗dt=p⃗E,
where p0≡E, and we set the velocity of light c=1; this is the correct definition for the velocity of a free relativistic particle. Moreover, it follows that(1.7)ẋμẋμ=dxμdxμdτ2=pμpμM2.
With our choice of metric, dxμdxμ=-ds2 and pμpμ=-m2, where m is the classical experimentally measured mass of the particle (at a given instant of τ). We see from this that(1.8)ds2dτ2=m2M2,
and hence the proper time is not identical to the evolution parameter τ. In the case where m2=M2, it follows that ds=dτ, and we say that the particle is “on shell”.For example, in the case of an external potentialV(x), where we write x≡xμ, the Hamiltonian becomes(1.9)K=pμpμ2M+V(x),
so that, since K is a constant of the motion, m2 varies from point to point with the variations of V(x). It is important to recognize from this discussion that the observable particle mass depends on the state of the system (in the quantum theory, the expectation value of the operator pμpμ provides the expected value of the mass squared).Since, in nature, particles appear with fairly sharp mass values (not necessarily with zero spread), we may assume the existence of some mechanism which will drive the particle's mass back to its original mass-shell value (after the source responsible for the mass change ceases to act) so that the particle's mass shell is defined. We will not take such a mechanism into account explicitly here in developing the dynamical equations. We will assume that if this mechanism is working, it is a relatively smooth function (e.g., a minimum in free energy which is broad enough for our off-shell driving force to work fairly freely) (A relativistic Lee model has been worked out which describes a physical mass shell as a resonance, and therefore a stability point on the spectrum [15], but at this point it is not clear to us how this mechanism works in general.).The Stueckelberg formulation implies the existence of a fifth “electromagnetic” potential, through the requirement of gauge invariance, and there is a generalized Lorentz force which contains a term that drives the particle off shell, whereas the terms corresponding to the electric and magnetic parts of the usual Maxwell fields do not (for the nonrelativistic case, the electric field may change the energy of a charged particle, but not the magnetic field; the electromagnetic field tensor in our case is analogous to the magnetic field, and the new field strengths, derived from theτ dependence of the fields and the additional gauge field, are analogous to the electric field, as we will see).In Section2, we give the structure of the field equations and show that the standard Maxwell theory is properly contained in this more general framework.In Section3, we apply Green's functions to the current source provided by the relativistic particle. We obtain the equations of motion for a relativistic particle which is, in general, off-shell. As in Dirac's result, these equations are of third order in the evolution parameter and therefore are highly unstable. Our results exhibit what appears to be an attractor in this system.We introduce the quantity(1.10)ɛ=1+ẋμẋμ,
which measures the deviation of the motion of the charged particle from its “mass shell”, where we used the metric (-,+,+,+) for the four-vector product.In Section4, we derive an autonomous mass deviation equation (with coefficients independent of τ or xμ and its derivatives) for the development of ɛ in the specific case of a particle under the influence of self-interaction in the absence of external field.When the formal derivation of the autonomous equation is complete, a numerical analysis follows in Section5 and demonstrates the compatibility between the (four component) dynamical equation of motion and the autonomous equation.In Sections6 and 7, we analyze the motion of the (charged) particle under the influence of self-interaction in the absence of external fields. We show that when initial conditions are within a basin of attraction, the trajectory of a particle obeying the equations of motion developed in Sections 3 and 4 forms an attractor.We further show that the attractor has a degree of stochasticity which is dictated by the cutoff frequency one chooses. A closer to zero (and thus more physical) cutoff frequency causes fluctuations to appear in the attractor, while a higher cutoff frequency results in a non-fluctuating attractor. There is an opposite relation between stochasticity and stability. Thus, a stochastic attractor has shorter life span. However, even an attractor with a relatively high cutoff frequency (and hence a longer life span) eventually ceases to exist as the orbit runs away to infinity.Our analysis shows that the effective mass of the electron grows as it propagates towards the vicinity ofɛ=1 (the light cone). This growth can be seen either through the growing amplitude of the effective mass in the case of a relatively low cutoff frequency or through a nonfluctuating growth of the effective mass (absolute) value in the case of a relatively high cutoff frequency.We demonstrate that the growing effective mass stabilizes the effect of the generalized Lorentz force. The demonstration consists of a correlation between the growing effective mass, the growing trajectory stability (measured through calculating the Lyapunov exponent), and the decay of acceleration and velocity.The chaotic nature of the electron's trajectory is also checked, using two methods: the (mentioned above) Lyapunov exponent calculation and a time series analysis. As will be shown, the results do not indicate a chaotic nature of the orbit, even in fluctuating cases.However, the nature of the attractor changes in the presence of an external field, as shown in Section8. In that section we show some results with remarkable stochasticity. This is reflected in the phase space graphs of the dynamical equation of motion, which show fractal behavior. Thus, the simple attractor of the pure self-interaction becomes a strange attractor (due to the fractal nature) in the presence of an external field. Experiments with high-frequency detectors may pick up such stochastic-type signals.


## 2. Equations of Motion

The Stueckelberg-Schrödinger equation which governs the evolution of a quantum state over the manifold of spacetime was postulated by Stueckelberg [10–14] to be, for the free particle,(2.1)i∂ψτ∂τ=pμpμ2Mψτ,
where ipμ is represented by ∂/∂xμ≡∂μ (we take ℏ=1).Taking into account fullU(1) gauge invariance, corresponding to the requirement that the theory maintain its form under the replacement of ψ by eie0Λψ, the Stueckelberg-Schrödinger equation (including a compensation field for the τ-derivative of Λ) [16–18] is (2.2)(i∂∂τ+e0a5(x,τ))ψτ(x)=(pμ-e0aμ(x,τ))(pμ-e0aμ(x,τ))2Mψτ(x),
where the gauge fields may depend on τ and e0 is a coupling constant which we will see having the dimension l-1. The corresponding classical Hamiltonian then has the form (2.3)K=(pμ-e0aμ(x,τ))(pμ-e0aμ(x,τ))2M-e0a5(x,τ)
in place of (2.1). Stuckelberg did not take into account this full gauge invariance requirement, working in the analog of what is known in the nonrelativistic case as the Hamilton gauge (where the gauge function Λ is restricted to be independent of time) and therefore had some difficulty in accounting for pair creation and annihilation in an electromagnetic field. The equations of motion for the field variables are given (for both the classical and quantum theories) by(2.4)λ∂αfβα(x,τ)=e0jβ(x,τ),
where α,β=0,1,2,3,5, the last corresponding to the τ index and λ, of dimension l-1, is a factor on the terms fαβfαβ in the Lagrangian associated with (2.2) (with, in addition, degrees of freedom of the fields) required by dimensionality. The field strengths are (2.5)fαβ=∂αaβ-∂βaα,
and the current satisfies the conservation law [16–18](2.6)∂αjα(x,τ)=0.
Integrating over τ on (-∞,∞) and assuming that j5(x,τ) vanishes at |τ|→∞, one finds that(2.7)∂μJμ(x)=0,
where [19](2.8)Jμ(x)=∫-∞∞dτjμ(x,τ).
We identify this Jμ(x) with the Maxwell conserved current. In [20], for example, this expression occurs with (2.9)jμ(x,τ)=ẋ(τ)δ4(x-x(τ)),
and τ is identified with the proper time of particle (an identification which can be made for the motion of a free particle). The conservation of the integrated current then follows from the fact that(2.10)∂μjμ=ẋμ(τ)∂μδ4(x-x(τ))=-ddτδ4(x-x(τ))
is a total derivative. We assume that the world line runs to infinity (at least in the time dimension) and therefore its integral vanishes at the end points [10–14, 20], in accordance with the discussion above.As for the Maxwell case, one can write the current formally in five-dimensional form(2.11)jα=ẋαδ4(x(τ)-x).
For α=5, the factor ẋ5 is unity, and this component therefore represents the event density in spacetime.Integrating theμ components of (2.4) over τ (assuming fμ5(x,τ)→0 pointwise for τ→±∞), we obtain the Maxwell equation with the (dimensionless) Maxwell charge e=e0/λ and the Maxwell fields given by(2.12)Aμ=∫-∞∞aμ(x,τ)dτ.A Hamiltonian in the form (2.3) without τ dependence of the fields, and without the a5 terms, as written by Stuckelberg [10–14], can be recovered in the limit of the zero mode of the fields (with a5=0) in a physical state for which this limit is a good approximation. This state occurs when the Fourier transform of the fields, defined by(2.13)aμ(x,τ)=∫-∞∞dsâμ(x,s)e-isτ,
has support only in a small neighborhood Δs of s=0. The vector potential then takes on the form aμ(x,τ)≈Δsâμ(x,0)=(Δs/2π)Aμ(x), and we identify e=(Δs/2πλ)e0. The zero mode therefore emerges when the inverse correlation length of the field Δs is sufficiently small. We remark that in this limit, the fifth equation obtained from (2.4) decouples. The generalized Lorentz force obtained from this Hamiltonian, using the Hamilton equations, coincides with the usual Lorentz force, and, as we have seen, the generalized Maxwell equations reduce to the usual Maxwell equations. The theory therefore contains the usual Maxwell Lorentz theory in the zero mode; for this reason we have called this generalized theory the “pre-Maxwell” theory.If such a pre-Maxwell theory really underlies the standard Maxwell theory, then there should be some physical mechanism which restricts most observations in the laboratory to be close to the zero mode. For example, in a metal there is a frequency, the plasma frequency, above which there is no transmission of electromagnetic waves. In this case, if the physical universe is imbedded in a medium which does not allow high “frequencies” to pass, the pre-Maxwell theory reduces to the Maxwell theory. Some study has been carried out, for a quite different purpose (of achieving a form of analog gravity), on the properties of the generalized fields in a medium with general dielectric tensor [21]. We will see in the present work that the high level of nonlinearity of this theory in interaction with matter may itself generate an effective reduction to Maxwell-Lorentz theory, with the high frequency chaotic behavior providing the regularization achieved by models of the type discussed by Rohrlich [7–9].We remark that integration overτ does not bring the generalized Lorentz force into the form of the standard Lorentz force, since it is nonlinear, and a convolution remains. If the resulting convolution is trivial, that is, in the zero mode, the two theories then coincide. Hence, we expect to see dynamical effects in the generalized theory which are not present in the standard Maxwell-Lorentz theory.Writing the Hamilton equations(2.14)ẋμ=dxμdτ=∂K∂pμ,ṗμ=dpμdτ=-∂K∂xμ,
for the Hamiltonian (2.3), we find that the generalized Lorentz force is(2.15)Mẍμ=e0fνμẋν+f5μ.
Multiplying this equation by ẋμ, one obtains(2.16)Mẋμẍμ=e0ẋμf5μ.
This equation therefore does not necessarily lead to the simple relation between ds and dτ discussed above in connection with (1.9). The f5μ term has the effect of moving the particle off-shell (as, in the nonrelativistic case, the energy is altered by the electric field). This can be seen by noting that
(2.17)ẋμẋμ=-m2M2,
and the derivative is
(2.18)ẋμẍμ=-1M2ddτm2=1Me0ẋμf5μ,
or in another words,
(2.19)ddτm2=-Me0ẋμf5μ.Let us now define(2.20)ɛ=1+ẋμẋμ=1-ds2dτ2,
where ds2=dt2-dx⃗2 is the square of the proper time. Since ẋμ=(pμ-e0aμ)/M and (pμ-e0aμ)(pμ-e0aμ)=-m2 is interpreted as the gauge-invariant particle mass [16–18], then(2.21)ɛ=1-m2M2
measures the deviation from “mass shell” (on mass shell, ds2=dτ2).


## 3. Derivation of the Differential Equations for the Spacetime Orbit with Off-Shell Corrections

We now review the derivation of the radiation reaction formula in the Stueckelberg formalism (see also [19]). Calculating the self-interaction contribution, one must include the effects of the force acting upon the particle due to its own field (fself) in addition to the fields generated by other electromagnetic sources (fext). Therefore, the generalized Lorentz force, using (2.14) takes the form:(3.1)Mẍμ=e0ẋνfextνμ+e0ẋνfselfνμ+e0fext5μ+e0fself5μ,
where the derivatives (dot) are with respect to the universal time τ and the fields are evaluated on the event's trajectory. Multiplying (3.1) by ẋμ, we get the projected equation (2.16) in the form(3.2)M2ɛ̇=e0ẋμfext5μ+e0ẋμfself5μ.The field generated by the current is given by the pre-Maxwell equations (2.4), and choosing for it the generalized Lorentz gauge ∂αaα=0, we get(3.3)λ∂α∂αaβ(x,τ)=λ(σ∂τ2-∂t2+∇2)aβ=-e0jβ(x,τ),
where σ=±1 corresponds to the possible choices of metric for the O(4,1) or O(3,2) symmetry of the homogeneous field equations.Green's functions for (3.3) can be constructed from the inverse Fourier transform(3.4)G(x,τ)=1(2π)5∫d4kdkei(kμxμ+σkτ)kμkμ+σk2.
Integrating this expression over all τ gives Green's function for the standard Maxwell field (treatment of the zeros in the denominator must be reconsidered after the τ integration). Assuming that the radiation reaction acts causally in τ, we will use here the τ-retarded Green's function. In this calculation of the Lorentz force, Dirac used the difference between advanced and retarded Green's functions in order to cancel the singularities that they contain. One can, alternatively, use the retarded Green's function and “renormalize” the mass in order to eliminate the singularity [3]. In this analysis, we follow the latter procedure.Theτ-retarded Green's function [19] is given by multiplying the principal part of the integral (3.4) by θ(τ). Carrying out the integrations (on a complex contour in k we consider the case σ=+1 in what follows), one finds (this Green's function differs from the t-retarded Green's function, constructed on a complex contour in k0)(3.5)G(x,τ)=2θ(τ)(2π)3{tan-1(-x2-τ2/τ)(-x2-τ2)3/2-τx2(x2+τ2),x2+τ2<0,121(x2+τ2)3/2ln⁡|τ-τ2+x2τ+τ2+x2|-τx2(x2+τ2),x2+τ2>0,
where we have written x2≡xμxμ.With the help of this Green's function, the solutions of (3.3) for the self-fields (substituting the current from (2.11)) are(3.6)aselfμ(x,τ)=e0λ∫d4x′dτ′G(x-x′,τ-τ′)ẋμ(τ′)δ4(x′-x(τ′))=e0λ∫dτ′ẋμ(τ′)G(x-x(τ′),τ-τ′),aself5(x,τ)=e0λ∫d4x′dτ′G(x-x(τ′),τ-τ′).The Green's function is written as a scalar, acting in the same way on all five components of the sourcejα; to assure that the resulting field is in Lorentz gauge, however, it should be written as a five by five matrix, with the factor δβα-kαkβ/k2(k5=k) included in the integrand. Since here we compute only the gauge-invariant field strengths, this extra term will not influence any of the results. It then follows that the generalized Lorentz force for the self-action (the force of the fields generated by the world line on a point xμ(τ) of the trajectory), along with the effect of external fields, is (3.7)Mẍμ=e02λ∫dτ′(ẋν(τ)ẋν(τ′)∂μ-ẋν(τ)ẋμ(τ′)∂ν)G(x-x(τ′))|x=x(τ)+e02λ∫dτ′(∂μ-ẋμ(τ′)∂τ)G(x-x(τ′))|x=x(τ)+e0(fextνμẋν+fext5μ).
We define u≡(xμ(τ)-xμ(τ′))(xμ(τ)-xμ(τ′)), so that(3.8)∂μ=2(xμ(τ)-xμ(τ′))∂∂u.
Equation (3.7) then becomes(3.9)Mẍμ=2e02λ∫dτ′[ẋν(τ)ẋν(τ′)(xμ(τ)-xμ(τ′))-ẋν(τ)ẋμ(τ′)(xν(τ)-xν(τ′))]∂∂uG(x-x(τ′))|x=x(τ)+e02λ∫dτ′[2(xμ(τ)-xμ(τ′))∂∂u-ẋμ(τ′)∂τ]G(x-x(τ′),τ-τ′)|x=x(τ)+e0(fextνμẋν+fext5μ).In the self-interaction problem whereτ→τ′,xμ(τ′)-xμ(τ)→0, Green's function is very divergent. Therefore one can expand all expressions in τ′′=τ-τ′ assuming that the dominant contribution is from the neighborhood of small τ′′. The divergent terms are later absorbed into the mass and charge definitions leading to renormalization (effective mass and charge). Expanding the integrands in Taylor series around the most singular point τ=τ′ and keeping the lowest-order terms, the variable u reduces to(3.10)u≅ẋμẋμτ′′2-ẋμẍμτ′′3+13ẋμx⃛μτ′′4+14ẍμẍμτ′′4.We now recall the definition of the off-shell deviationɛ given in (2.20), along with its derivatives:(3.11)ẋμẋμ=-1+ɛ,ẋμẍμ=12ɛ̇,ẋμx⃛μ+ẍμẍμ=12ɛ̈.
Next, we define(3.12)ω≡112ẋμx⃛μ+18ɛ̈,Δ≡-12ɛ̇τ′′+ωτ′′2.
Using these definitions along with those of (3.10) and (3.11), we find(3.13)u+τ2τ′′2≅ɛ+Δ.
We then expand Green's function to leading orders:(3.14)∂G∂u≅Θ(τ′′)f1(ɛ+Δ)(2π)3τ′′5=Θ(τ′′)(2π)3[f1(ɛ)τ′′5-ɛ̇f1′(ɛ)2τ′′4+(ωf1′(ɛ)+18ɛ̇2f1′′(ɛ))1τ′′3],∂G∂τ′′≅2Θ(τ′′)f2(ɛ+Δ)(2π)3τ′′4+2δ(τ′′)f3(ɛ+Δ)(2π)3τ′′3=2Θ(τ′′)(2π)3[f2(ɛ)τ′′4-ɛ̇f2′(ɛ)2τ′′3+(ωf2′(ɛ)+18ɛ̇2f2′′(ɛ))1τ′′2]+2δ(τ′′)(2π)3[f3(ɛ)τ′′3-ɛ̇f3′(ɛ)2τ′′2+(ωf3′(ɛ)+18ɛ̇2f3′′(ɛ))1τ′′],
where f′≡df/dɛ. For ɛ<0(3.15)f1(ɛ)=3tan-1(-ɛ)(-ɛ)5/2-3ɛ2(1-ɛ)+2ɛ(1-ɛ)2,f2(ɛ)=3tan-1(-ɛ)(-ɛ)5/2-1ɛ2-2-ɛɛ2(1-ɛ),f3(ɛ)=tan-1(-ɛ)(-ɛ)3/2+1ɛ(1-ɛ),
and for ɛ>0(3.16)f1(ɛ)=(3/2)ln⁡|(1+ɛ)/(1-ɛ)|ɛ5/2-3ɛ2(1-ɛ)+2ɛ(1-ɛ)2,f2(ɛ)=(3/2)ln⁡|(1+ɛ)/(1-ɛ)|ɛ5/2-1ɛ2-2-ɛɛ2(1-ɛ),f3(ɛ)=-(1/2)ln⁡|(1+ɛ)/(1-ɛ)|ɛ3/2+1ɛ(1-ɛ).
For either sign of ɛ, when ɛ≈0,(3.17)f1(ɛ)≈85+247ɛ+163ɛ2+O(ɛ3),f2(ɛ)≈-25-47ɛ-23ɛ2+O(ɛ3),f3(ɛ)≈23+45ɛ+67ɛ2+O(ɛ3).
One sees that the derivatives in (3.14) have no singularity in ɛ at ɛ=0.From (3.6), we have(3.18)fself5μ(x(τ),τ)=e∫[2(xμ(τ)-xμ(τ′))∂∂u-ẋμ(τ′)∂τ]×G(x-x(τ′),τ-τ′)|x=x(τ)dτ′.
We then expand xμ(τ)-xμ(τ′) and ẋμ(τ)-ẋμ(τ′) in power series in τ′′ and write the integrals formally with infinite limits.Substituting (3.18) into (3.2), we obtain (note that xμ and its derivatives are evaluated at the point τ and are not subject to the τ′′ integration), after integrating by parts using δ(τ′′)=(∂/∂τ′′)θ(τ′′),(3.19)M2ɛ̇=2e02λ(2π)3∫-∞∞dτ′′[g1τ′′4(ɛ-1)-g2τ′′3ɛ̇2+g3τ′′2ẋνx⃛ν-h1τ′′3(ɛ-1)ɛ̇+h22τ′′2ɛ̇2+h3τ′′2(ɛ-1)ω+h48τ′′2ɛ̇2(ɛ-1)]Θ(τ′′)+e0ẋμfext5μ,
where we have defined(3.20)g1=f1-f2-3f3,g2=12f1-f2-2f3,g3=16f1-12f2-12f3,h1=12f1′-12f2′-f3′,h2=14f1′-12f2′-12f3′,h3=(f1′-f2′-f3′),h4=f1′′-f2′′-f3′′.The integrals are divergent at the lower boundτ′′=0 imposed by the θ-function; we therefore take these integrals to a cut-off μ>0. Equation (3.19) then becomes(3.21)M2ɛ̇=2e02λ(2π)3[g13μ3(ɛ-1)-g24μ2ɛ̇+g3μẋνx⃛ν-h12μ2(ɛ-1)ɛ̇+h22μɛ̇2+h3μ(ɛ-1)ω+h48μɛ̇2(ɛ-1)]+e0ẋμfext5μ.Following a similar procedure, we obtain from(3.7)(3.22)Mẍμ=2e02λ(2π)3[-12((1-ɛ)ẍμ+ɛ̇2ẋμ)(f12μ2-ɛ̇f1′2μ)+f13μ(ẋνx⃛νẋμ+(1-ɛ)x⃛μ)+g13μ3ẋμ-g22μ2ẍμ+g3μx⃛μ-h1ɛ̇2μ2ẋμ+h2ɛ̇μẍμ+h3ωμẋμ+h4ɛ̇28μẋμ]+e0(fextνμẋν+fext5μ).Substituting (3.21) for the coefficients of the ẋμ terms in the second line of (3.22), we find(3.23)M(ɛ,ɛ̇)ẍμ=-12M(ɛ,ɛ̇)(1-ɛ)ɛ̇ẋμ+2e02λ(2π)3μF(ɛ)[x⃛μ+1(1-ɛ)ẋνx⃛νẋμ]+e0ẋμẋνfext5ν1-ɛ+e0fextνμẋν+e0fext5μ,
where(3.24)F(ɛ)=f13(1-ɛ)+g3.Here, the coefficients ofẍμ have been grouped formally into a renormalized (off-shell) mass term, defined (as done in the standard radiation reaction problem) as(3.25)M(ɛ,ɛ̇)=M+e′22μ[f1(1-ɛ)2+g2]-e′2[14f1′(1-ɛ)+h2]ɛ̇.
We have here a Lorentz effective charge as(3.26)e′2=2e02λ(2π)3μ,
corresponding to a renormalization depending as well on the cut-off. We will call this quantity e in the following formulas, since, in this context, it can be identified with the measured electric charge through the radiation reaction formulas.We remark that one can change variables, with the help of (2.20) (here, for simplicity, assuming ɛ<1), to obtain a differential equation in which all derivatives with respect to τ are replaced by derivatives with respect to the proper time s. The coefficient of the second derivative of xμ with respect to s, as well as “effective mass” for the proper-time equation, is then given by(3.27)MS(ɛ,ɛ̇)=23F(ɛ)(1-ɛ)[M+e22c3μ[f1(1-ɛ)2+g2]-e2c3[14f1(1-ɛ)+h2+32F(ɛ)1-ɛ]ɛ̇].Note that the renormalized mass depends onɛ(τ); for this quantity to act as a mass, ɛ must be slowly varying on some interval on the orbit of the evolution compared to all other motions. The computer analysis we give below indeed shows that there are large intervals of almost constant ɛ. In case, as at some points, ɛ may be rapidly varying, one may consider the definition (3.25) as formal; clearly, however, if M(ɛ,ɛ̇) is large, ẍμ will be suppressed (e.g., for ɛ close to unity, where MS(ɛ,ɛ̇) goes as ɛ̇/(1-ɛ)2 and M(ɛ,ɛ̇) as ɛ̇/(1-ɛ)3; note that F(ɛ) goes as (1-ɛ)-3).We now obtain from  (3.23)(3.28)M(ɛ,ɛ̇)ẍμ=-12M(ɛ,ɛ̇)(1-ɛ)ɛ̇ẋμ+F(ɛ)e2[x⃛μ+1(1-ɛ)ẋνx⃛νẋμ]+e0fextνμẋν+e0(ẋμẋν1-ɛ+δνμ)fext5ν.We remark that when one multiplies this equation byẋμ, it becomes an identity (all of the terms except for fextνμẋν may be grouped to be proportional to (ẋμẋν/(1-ɛ)+δνμ)); one must use (3.21) to compute the off-shell mass shift ɛ corresponding to the longitudinal degree of freedom in the direction of the four-velocity of the particle. Equation (3.28) determines the motion orthogonal to the four-velocity. Equations (3.21) and (3.28) are the fundamental dynamical equations governing the off-shell orbit.We remark that as in [19] it can be shown that (3.28) reduces to the ordinary (Abraham-Lorentz-Dirac) radiation reaction formula for small, slowly changing ɛ and that no instability, no radiation, and no acceleration of the electron occur when it is on shell. There is therefore no “runaway solution” for the exact mass-shell limit of this theory; the unstable Dirac result is approximate for ɛ close to, but not precisely, zero.


## 4. Theɛ Evolution

We now derive an equation for the evolution of the off-shell deviation,ɛ, when the external field is removed. We then use this equation to prove that a fixed mass shell is consistent only if the particle is not accelerating, and therefore no runaway solution occurs. Using the definitions(4.1)F1(ɛ)=g1(ɛ-1)3μ2,F2(ɛ)=g2+2(ɛ-1)h14μ,F3(ɛ)=g3+112(ɛ-1)h3,F4(ɛ)=12h2+18(ɛ-1)h4,F5(ɛ)=18(ɛ-1)h3
in (3.21), in the absence of external fields, we write(4.2)ẋμx⃛μ=1F3(ɛ){M2e2ɛ̇-F1(ɛ)+F2(ɛ)ɛ̇-F4(ɛ)ɛ̇2-F5(ɛ)ɛ̈}.
Differentiating with respect to τ, we find that(4.3)ẋμx....μ+ẍμx⃛μ=1F3{F2′ɛ̇2+ɛ̈(M2e2+F2)-F1′ɛ̇-F4′ɛ̇3-(2F4+F5′)ɛ̇ɛ̈-F5ɛ⃛}-F3′F32{F2ɛ̇+M2e2ɛ̇-F1-F4ɛ̇2-F5ɛ̈}ɛ̇≡H-F5F3ɛ⃛.
Together with (4.4), which is the τ derivative of the last equation in (3.11), one finds (using (4.3)) that(4.4)ẋμx....μ+3ẍμx⃛μ=12ɛ⃛,
which is the τ derivative of the last equation in (3.11), one finds using(4.3)(4.5)ẍμx⃛μ=(14+F52F3)-12H(ɛ,ɛ̇,ɛ̈).
Multiplying (3.28) by ẍμ (with no external fields) and using (4.2) and (4.5), we obtain(4.6)(1+2F5F3)ɛ⃛-A(ɛ)ɛ̈+B(ɛ)ɛ̇2+C(ɛ)ɛ̇-D(ɛ)+E(ɛ)ɛ̇3+I(ɛ)ɛ̇ɛ̈=0,
where(4.7)A=2F3(M2e2+F2)+2M(ɛ,ɛ̇)e2F(ɛ)+4M(ɛ,ɛ̇)F5F3e2F(ɛ),B=2F3′F32(F2+M2e2)-2F2′F3+21-ɛ1F3(M2e2+F2)-M(ɛ,ɛ̇)e2F(ɛ)(1-ɛ)-4F4M(ɛ,ɛ̇)F3e2F(ɛ),C=4M(ɛ,ɛ̇)e2F(ɛ)F3(M2e2+F2)-2F32F1F3′-2F1(1-ɛ)F3+2F3F1′,D=4M(ɛ,ɛ̇)F1e2F(ɛ)F3,E=-2F4(1-ɛ)F3-2(F4F3′F32-F4′F3),I=2(F5′F3+2F4F3-F5F3′F32)-2F5(1-ɛ)F3.
We will call this the autonomous equation.


## 5. The Mass Correction Equation and the Orbit Equation: A Compatibility Check

One needs to verify that the off-shell mass correction equation (4.6) gives the same results as the four-component orbit equation (3.22).In the absence of external fields, (3.22) can also be written as(5.1)aμ(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3)=f13μ(ẋνx⃛νẋμ+(1-ɛ)x⃛μ)+g3μx⃛μ+h3ẋνx⃛ν3μẋμ,
where(5.2)aμ(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3)≡M⋅λ(2π)32e02ẍμ+12((1-ɛ)ẍμɛ̇2ẋμ)(f12μ2-ɛ̇f1′2μ)-g13μ3ẋμ+g22μ2ẍμ+h1ɛ2μ2ẋμ-h2ɛμẍμ-h3ẍνẍν4μẋμ-h4ɛ̇28μẋμ.
In order to obtain x⃛μ as a function of the components of the first and second derivatives, one may treat this system as a nonhomogeneous, linear system:(5.3)Ax⃛=a,
where A is a 4 × 4 matrix.The solutions for this system, which were calculated using the computer packageMathematica, yield initial values for x⃛μ for a set of initial numeric values assigned in ẋμ and ẍμ.Differentiating (5.1) with respect to τ, we find(5.4)ȧμ(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3,x⃛0,x⃛1,x⃛2,x⃛3)=f13μ(ẋνx....νẋμ+(1-ɛ)x....μ)+g3μx....μ+h3ẋνx....ν3μẋμ.The initial valuesx⃛0,x⃛1,x⃛2,x⃛3 of (5.4) are the solutions to (5.1) with initial values for ẋ0,…,ẋ3 and ẍ0,…,ẍ3 (of course (5.1) and (5.4) have the same initial values of ẋ0,…,ẋ3 and ẍ0,…,ẍ3). Equations (5.4) and (5.1) have the same form, and therefore the solutions for x⃛μ and x....μ also have the same form, with the only difference being in replacing aμ (5.2) with ȧμ(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3,x⃛0,x⃛1,x⃛2,x⃛3).Now, we have to show that (4.4) is satisfied. One has(5.5)ẋμx....μ+3ẍμx⃛μ=12ɛ⃛,
where x⃛μ(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3) are the solutions of (5.1), x....μ(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3,x⃛0,x⃛1,x⃛2,x⃛3) are the solutions of (5.4), and ɛ⃛(ɛ,ɛ̇,ɛ̈) is the solution of the autonomous equation with initial conditions which are consistent with the ones in (5.1) and (5.4) that satisfy (3.11).In Tables1 and 2, we see two sets of parameters and initial conditions which were chosen for the components of ẋμ and ẍμ. For each set, the orbit equation yields its results for x⃛μ and xμ···· and the autonomous equation yields its results for ɛ⃛. In both sets the results satisfy (5.5).Table 1
First set. The parameters:μ=3*10-3 GeV-1, M⁄α=4.84473*105  GeV2.
Initial conditions for the orbitThe third and forth derivatives (calculated with the initialequationconditions for the orbit equation)ẋ0=1.6ẍ0=15x⃛0=-394197.84731x....0=5.85464*109ẋ1=0.1ẍ1=10x⃛1=20879.93019x....1=5.62709*108ẋ2=0.3ẍ2=-13x⃛2=-153331.92947x....2=7.54373*108ẋ3=0.9ẍ3=-23x⃛3=379634.21816x....3=6.46893*109The orbit equation resultɛ⃛=2·(ẋμx⃛μ+3ẍμx⃛μ)=-1.726245753395*109Initial conditions for the autonomous equationɛ=1-ẋ02+ẋ12+ẋ22+ẋ32=-0.65ɛ̇=2·(-ẋ0·ẍ0+ẋ1·ẍ1+ẋ2·ẍ2+ẋ3·ẍ3)=-95.2ɛ̈=2·(-ẋ0·x⃛0+ẋ1·x⃛1+ẋ2·x⃛2+ẋ3·x⃛3-ẍ0·ẍ0+ẍ1·ẍ1+ẍ2·ẍ2+ẍ3·ẍ3)=1.8581*106The autonomous equation resultɛ⃛=-1.726245753395*109Table 2
Second set. The parameters:μ=3*105  GeV-1 Hz, M/α=4.84473*10-10  GeV2.
Initial conditions for the orbitThe third and forth derivatives (calculated with the initialequationconditions for the orbit equation)ẋ0=1.001ẍ0=1.6476x⃛0=22.5575x....0=124.9841ẋ1=0.1ẍ1=4x⃛1=11.9631x....1=58.5639ẋ2=0.2ẍ2=1x⃛2=6.2052x....2=33.0308ẋ3=0.3ẍ3=2x⃛3=10.5736x....3=55.5532The orbit equation resultɛ⃛=2·(ẋμx....μ+3ẍμx⃛μ)=36.272327445316Initial conditions for the autonomous equationɛ=1-ẋ02+ẋ12+ẋ22+ẋ32=0.13800ɛ̇=2·(-ẋ0·ẍ0+ẋ1·ẍ1+ẋ2·ẍ2+ẋ3·ẍ3)=-0.89850ɛ̈=2·(-ẋ0·x⃛0+ẋ1·x⃛1+ẋ2·x⃛2+ẋ3·x⃛3-ẍ0·ẍ0+ẍ1·ẍ1+ẍ2·ẍ2+ẍ3·ẍ3)=2.62951The autonomous equation resultɛ⃛=36.272327445316As a further test we compared trajectories of the orbit equations and the autonomous equation using the Runge-Kutta method.


### 5.1. The Numerical Solution Techniques

We investigated the dynamical behavior of the system using two methods: the first is the completely “computerized” method and the second is the “manual” fourth-order Runge-Kutta algorithm.The first method is based uponMathematica's numerical equation-solving set of commands. This method was only used with the autonomous equation (as the orbit equations are too complicated and thus require a different approach). The second “manual” method was used to numerically solve both the orbit equation and the autonomous equation (also using Mathematica's computation power), as follows.The four dynamical equations (or orbit equations), (3.22), can be written as a set of eight first-order equations of the form (5.6)dẋ0dτ=ẍ0,dẋ1dτ=ẍ1,dẋ2dτ=ẍ2,dẋ3dτ=ẍ3,dẍ0dτ=f0(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3),dẍ1dτ=f1(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3),dẍ2dτ=f2(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3),dẍ3dτ=f3(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3).We assigned the expressions obtained earlier for the four third-order derivativesx⃛0-x⃛3 into the functions f0-f3, respectively, and used the Runge-Kutta method to calculate the trajectory for the system of equations. The algorithm was programmed in Mathematica.As for the autonomous equation, (4.6), it can be written as a set of three first-order equations of the form(5.7)dɛdτ=ɛ̇,dɛ̇dτ=ɛ̈,dɛ̈dτ=Q(ɛ,ɛ̇,ɛ̈),
with
(5.8)Q=ɛ⃛=11+2(F5/F3)(A(ɛ,ɛ̇)ɛ̈-B(ɛ,ɛ̇)ɛ̇2-C(ɛ,ɛ̇)ɛ̇+D(ɛ,ɛ̇)-E(ɛ)ɛ̇3-I(ɛ)ɛ̇ɛ̈).Again, the initial valuesɛ,ɛ̇,ɛ̈ for the autonomous quation were chosen to match those of the orbit equation, according to (3.11). Thus, they obey the three conditions(5.9)ɛ=ẋμẋμ+1,ɛ̇=2ẋμẍμ,ɛ̈=2(ẋμx⃛μ+ẍμẍμ).With these initial values inserted, we calculated the trajectory for the autonomous equation using the Runge-Kutta method.As explained above, we used two sets of initial values, which appear in Tables1 and 2, and evolve to the trajectories shown in Figures 1 and 2, respectively. Figures 1(a) and 2(a) show the trajectories of the autonomous equation for each set, and Figures 1(b) and 2(b) show those of the orbit equation. The identical results demonstrate the compatibility between the orbit equation and the autonomous equation.(a) First set, the autonomous equation trajectory. (b) First set, the orbit equation trajectory.
(a)(b)(a) Second set, the autonomous equation trajectory. (b) Second set, the orbit equation trajectory.
(a)(b)Mathematica's numerical equation-solving commands produced identical results to the Runge-Kutta method for the autonomous equation in both sets of the initial conditions.


### 5.2. About the Parameters Used

The calculation carried out uses the parametersM/α and μ, where α=μe2 (the orbit equation contains the factor α on the right-hand side, so that M/α is a convenient parameter for the numerical analysis). We chose to work with natural units: ℏ=c=1, and so energy has units of GeV, and mass (m and M) is also in units of GeV. Length (ℏ/mc) is in units of GeV-1, and time (ℏ/mc2) is in units of  GeV-1 (both proper and universal). The velocity is therefore without units, acceleration has units of GeV, and so on. The cutoff frequency, μ, has units of time, GeV-1, and the charge, e=2e02/λ(2π)3μ, has no units. The value e2 can be identified with the measured electric charge squared. It depends on three bare parameters: e0, μ, and λ. In our examples we chose these bare parameters so that e would be equal to the electron's charge. The initial conditions were chosen to be close to the mass shell (ɛ=1-m2/M2≈0).


## 5.1. The Numerical Solution Techniques

We investigated the dynamical behavior of the system using two methods: the first is the completely “computerized” method and the second is the “manual” fourth-order Runge-Kutta algorithm.The first method is based uponMathematica's numerical equation-solving set of commands. This method was only used with the autonomous equation (as the orbit equations are too complicated and thus require a different approach). The second “manual” method was used to numerically solve both the orbit equation and the autonomous equation (also using Mathematica's computation power), as follows.The four dynamical equations (or orbit equations), (3.22), can be written as a set of eight first-order equations of the form (5.6)dẋ0dτ=ẍ0,dẋ1dτ=ẍ1,dẋ2dτ=ẍ2,dẋ3dτ=ẍ3,dẍ0dτ=f0(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3),dẍ1dτ=f1(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3),dẍ2dτ=f2(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3),dẍ3dτ=f3(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3).We assigned the expressions obtained earlier for the four third-order derivativesx⃛0-x⃛3 into the functions f0-f3, respectively, and used the Runge-Kutta method to calculate the trajectory for the system of equations. The algorithm was programmed in Mathematica.As for the autonomous equation, (4.6), it can be written as a set of three first-order equations of the form(5.7)dɛdτ=ɛ̇,dɛ̇dτ=ɛ̈,dɛ̈dτ=Q(ɛ,ɛ̇,ɛ̈),
with
(5.8)Q=ɛ⃛=11+2(F5/F3)(A(ɛ,ɛ̇)ɛ̈-B(ɛ,ɛ̇)ɛ̇2-C(ɛ,ɛ̇)ɛ̇+D(ɛ,ɛ̇)-E(ɛ)ɛ̇3-I(ɛ)ɛ̇ɛ̈).Again, the initial valuesɛ,ɛ̇,ɛ̈ for the autonomous quation were chosen to match those of the orbit equation, according to (3.11). Thus, they obey the three conditions(5.9)ɛ=ẋμẋμ+1,ɛ̇=2ẋμẍμ,ɛ̈=2(ẋμx⃛μ+ẍμẍμ).With these initial values inserted, we calculated the trajectory for the autonomous equation using the Runge-Kutta method.As explained above, we used two sets of initial values, which appear in Tables1 and 2, and evolve to the trajectories shown in Figures 1 and 2, respectively. Figures 1(a) and 2(a) show the trajectories of the autonomous equation for each set, and Figures 1(b) and 2(b) show those of the orbit equation. The identical results demonstrate the compatibility between the orbit equation and the autonomous equation.(a) First set, the autonomous equation trajectory. (b) First set, the orbit equation trajectory.
(a)(b)(a) Second set, the autonomous equation trajectory. (b) Second set, the orbit equation trajectory.
(a)(b)Mathematica's numerical equation-solving commands produced identical results to the Runge-Kutta method for the autonomous equation in both sets of the initial conditions.


## 5.2. About the Parameters Used

The calculation carried out uses the parametersM/α and μ, where α=μe2 (the orbit equation contains the factor α on the right-hand side, so that M/α is a convenient parameter for the numerical analysis). We chose to work with natural units: ℏ=c=1, and so energy has units of GeV, and mass (m and M) is also in units of GeV. Length (ℏ/mc) is in units of GeV-1, and time (ℏ/mc2) is in units of  GeV-1 (both proper and universal). The velocity is therefore without units, acceleration has units of GeV, and so on. The cutoff frequency, μ, has units of time, GeV-1, and the charge, e=2e02/λ(2π)3μ, has no units. The value e2 can be identified with the measured electric charge squared. It depends on three bare parameters: e0, μ, and λ. In our examples we chose these bare parameters so that e would be equal to the electron's charge. The initial conditions were chosen to be close to the mass shell (ɛ=1-m2/M2≈0).


## 6. The Autonomous Equation and the (Universal Time) Effective MassM(ɛ,ɛ̇)

In both the orbit equation and the autonomous equation, initial conditions within a basin of attraction will evolve to an orbit of an attractor. In the phase space ofɛ, ɛ̇, ɛ̈, and ɛ⃛, the orbit of the attractor propagates towards the vicinity of the coordinate (ɛ,ɛ̇,ɛ̈,ɛ⃛)=(1,0,0,0), which will simply be named unity.The nature of the orbit depends on the cutoff frequencyμ. Higher cutoff frequency results in a nonfluctuating orbit, while a lower, more physical frequency (closer to zero) results in a fluctuating one. The fluctuations appear as loop cycles in the phase space graphs.


### 6.1. Analysis of a Cycling Attractor (with a Lower, More “Physical” Cutoff Frequency)

In this subsection, we bring an example of a fluctuant attractor's orbit and present its properties. We start with Figures3(a) and 3(b), which are the orbit's phase space graphs (ɛ, ɛ̇) and (ɛ̇, ɛ̈), respectively.(a) and (b) The orbit of an attractor in two phase space graphs. Initial conditions:ɛ=0.01, ɛ̇=55, ɛ̈=106, M/α=556.85  GeV2, and μ=10-5  GeV-1. (c) The Lyapunov exponent λ. (d) The correlation between two close orbits: the attractor and an orbit which starts close to it (distance 0.00001). One can clearly see the correlation between the orbits, which is not a chaotic pattern. (e) and (f) ɛ and the effective mass as a function of the number of time iterations. As ɛ becomes close to unity the mass fluctuation amplitude grows significantly.
(a)(b)(c)(d)(e)(f)In Figure3(c) we show the global Lyapunov exponent of the attractor. We have computed the global Lyapunov exponent by studying the average separation of orbits associated with nearby initial conditions. Segments of relative stability and instability can be observed.Despite the fluctuations which create temporary instabilities, in general the Lyapunov exponent decays. A time series of the correlation between the attractor and an orbit with an initial coordinate at a distance of 0.00001 (Figure3(d)) shows the same phenomenon: the distance between the orbits generally decreases, despite the temporary fluctuations. Thus, no evidence of chaotic behavior is found.Figures3(e) and 3(f) show the attractor's ɛ and the effective mass as a function of the universal time iterations. The effective mass fluctuation amplitude grows significantly near unity.We will analyze three consecutive cycles (fluctuations) in our chosen attractor and observe the phase space trajectory, its stability, and the mass behavior.


#### 6.1.1. First Cycle

In Figures4(a) and 4(b), we see ɛ,ɛ̇ and ɛ̇,ɛ̈ graphs of this cycle, respectively. Figures 4(c) and 4(d) show ɛ and the effective mass as functions of the number of iterations. Figure 4(e) shows the Lyapunov exponent. The orbit in Figure 4(a) enters a region of folding, which develops to a loop. In Figure 4(b), one can see an apparent attractor inducing motion beginning at the lower convex portion of the orbit which then turns back in a characteristic way to reach the last point visible on this graph.(a) and (b). Theɛ,ɛ̇ and the ɛ̇,ɛ̈ graphs, respectively. One can see the loop cycle in both phase space graphs. (c) and (d) ɛ and the (universal time) effective mass M(ɛ,ɛ̇) as functions of the number of iterations. One can see how ɛ temporarily decreases during the cycle. The effective mass significantly grows (in absolute value) as the cycle terminates. (e) The Lyapunov exponent fluctuates during the cycle (relative instability) but generally decreases.
(a)(b)(c)(d)(e)The instability in the Lyapunov exponent (Figure4(e)) is accompanied by a temporary decrease in the ɛ graph (Figure 4(c)). This instability is the climax of the loop cycle in the ɛ,ɛ̇ and the ɛ̇,ɛ̈ graphs. As the loop cycle terminates, ɛ climbs to a higher value than the one it had when the cycle began, and the exponent continues to decay.In Figure4(d) we see that the effective mass fluctuates with a growing amplitude. As the cycle reaches towards its end, the effective mass continues to a higher amplitude and has its highest absolute values.


#### 6.1.2. Second Cycle

The attractor's orbit continues to a second cycle, as can be seen in theɛ,ɛ̇ and ɛ̇,ɛ̈ graphs (Figures 5(a) and 5(b), resp.). This second cycle is again accompanied by a fluctuation (relative instability) in the Lyapunov exponent, as can be seen in Figure 5(e). The mass, as before, fluctuates with growing amplitude (we can observe the growing values on the y-axis in Figure 5(d) in comparison to Figure 4(d)). As the mass reaches a relatively low absolute value, the Lyapunov exponent rises to a relatively high one (the “jump” area of the exponent graph).(a) and (b) Theɛ,ɛ̇ and the ɛ̇,ɛ̈ phase space graphs, respectively. (c) and (d) ɛ and the (universal time) effective mass M(ɛ,ɛ̇) as functions of the number of iterations. As in the previous cycle, the mass fluctuation amplitude grows. (e) Again, the Lyapunov exponent fluctuates during the cycle and returns to decay as it ends.
(a)(b)(c)(d)(e)


#### 6.1.3. Third Cycle

The consecutive cycle behaves in the same way as the two previous ones, as can be seen in Figures6(a)–6(e). The occurrence of the present loop cycle is seen when ɛ has a down spike (around the iteration value 2500 in Figure 6(c)). In this area the Lyapunov exponent has an up spike (Figure 6(e)). The effective mass, in contrast to the Lyapunov exponent, decreases during the cycle but rises to a higher amplitude as it ends (Figure 6(d)).(a) and (b) Theɛ,ɛ̇ and the ɛ̇,ɛ̈ phase space graphs, respectively. (c) and (d) ɛ and the (universal time) effective mass M(ɛ,ɛ̇) as functions of the number of iterations. As in the previous cycles, the mass fluctuation amplitude significantly grows towards the cycle’s end. One should also observe the relatively low value of the effective mass around the iteration value of 2500 (the climax of the cycle). (e) Again, the Lyapunov exponent fluctuates during the cycle and returns to decay afterwards.
(a)(b)(c)(d)(e)The above results suggest that as long as the electron “survives” the fluctuations (or cycles) and its orbit doesn't run away to infinity, the effective mass will reach higher absolute values in the end of each cycle. The effective mass, in general, fluctuates with increasing amplitude. As can be seen from the Lyapunov exponent graphs, the attractor is relatively unstable during the cycles but regains stability between them.A  correlation  is  clearly  observed  between  a  relatively  high  absolute  value  of  the  effective  mass  and  a  relatively  low  value  of  the  Lyapunov  exponent  (which  indicates  a  relative  stability).We also saw correlation between our attractor and nearby orbits, as seen in the time series associated with them, even in the vicinity of unity (the light cone). This correlation shows the absence of chaotic behavior.


### 6.2. A Smooth (Less Physical) Attractor

In this subsection we give an example of a smooth attractor, as can be seen in Figures7(a)–7(e).(a) and (b) An attractor without fluctuations (cycles). (a)ɛ,ɛ̇ graph and (b) ɛ̇,ɛ̈ graph. (c) Lyapunov exponent and (d) ɛ itself, as functions of the (universal) time iterations. (e) The effective mass steadily grows (in absolute value) as ɛ moves towards unity.
(a)(b)(c)(d)(e)The only change in the initial conditions and parameters between this smooth orbit and the former, fluctuant, one is the cutoff frequency, which has a higher value now: 10-1 GeV-1 (instead of 10-5 GeV-1 in the fluctuant case). This orbit has no loop cycles in its phase space graphs. It has, however, a higher life span than the previous, fluctuant attractor, and thus it stays a longer time in the vicinity of unity. Nevertheless, also in this case the orbit will, at some point, run away to infinity.Although a higher cutoff frequency is less “physical”, the advantage of showing an attractor with relatively high cutoff frequency is that the fluctuations are removed, and the “pure” correlation between the growing effective mass and the rising stability (in the form of declining and even negative Lyapunov exponent) can be seen. Indeed, the attractor's Lyapunov exponent in Figure7(c) turns to be unambiguously negative and smooth, which may indicate the existence of a fixed point.The attractor'sɛ and effective mass as functions of the time iterations appear in Figures 7(d) and 7(e), respectively. There are no fluctuations. We can therefore conclude that the pattern is even clearer in the present example: as the attractor is steadily approaching the light cone, the effective mass is steadily growing (in absolute value), without fluctuations, and the Lyapunov exponent is negative and stable.


## 6.1. Analysis of a Cycling Attractor (with a Lower, More “Physical” Cutoff Frequency)

In this subsection, we bring an example of a fluctuant attractor's orbit and present its properties. We start with Figures3(a) and 3(b), which are the orbit's phase space graphs (ɛ, ɛ̇) and (ɛ̇, ɛ̈), respectively.(a) and (b) The orbit of an attractor in two phase space graphs. Initial conditions:ɛ=0.01, ɛ̇=55, ɛ̈=106, M/α=556.85  GeV2, and μ=10-5  GeV-1. (c) The Lyapunov exponent λ. (d) The correlation between two close orbits: the attractor and an orbit which starts close to it (distance 0.00001). One can clearly see the correlation between the orbits, which is not a chaotic pattern. (e) and (f) ɛ and the effective mass as a function of the number of time iterations. As ɛ becomes close to unity the mass fluctuation amplitude grows significantly.
(a)(b)(c)(d)(e)(f)In Figure3(c) we show the global Lyapunov exponent of the attractor. We have computed the global Lyapunov exponent by studying the average separation of orbits associated with nearby initial conditions. Segments of relative stability and instability can be observed.Despite the fluctuations which create temporary instabilities, in general the Lyapunov exponent decays. A time series of the correlation between the attractor and an orbit with an initial coordinate at a distance of 0.00001 (Figure3(d)) shows the same phenomenon: the distance between the orbits generally decreases, despite the temporary fluctuations. Thus, no evidence of chaotic behavior is found.Figures3(e) and 3(f) show the attractor's ɛ and the effective mass as a function of the universal time iterations. The effective mass fluctuation amplitude grows significantly near unity.We will analyze three consecutive cycles (fluctuations) in our chosen attractor and observe the phase space trajectory, its stability, and the mass behavior.


### 6.1.1. First Cycle

In Figures4(a) and 4(b), we see ɛ,ɛ̇ and ɛ̇,ɛ̈ graphs of this cycle, respectively. Figures 4(c) and 4(d) show ɛ and the effective mass as functions of the number of iterations. Figure 4(e) shows the Lyapunov exponent. The orbit in Figure 4(a) enters a region of folding, which develops to a loop. In Figure 4(b), one can see an apparent attractor inducing motion beginning at the lower convex portion of the orbit which then turns back in a characteristic way to reach the last point visible on this graph.(a) and (b). Theɛ,ɛ̇ and the ɛ̇,ɛ̈ graphs, respectively. One can see the loop cycle in both phase space graphs. (c) and (d) ɛ and the (universal time) effective mass M(ɛ,ɛ̇) as functions of the number of iterations. One can see how ɛ temporarily decreases during the cycle. The effective mass significantly grows (in absolute value) as the cycle terminates. (e) The Lyapunov exponent fluctuates during the cycle (relative instability) but generally decreases.
(a)(b)(c)(d)(e)The instability in the Lyapunov exponent (Figure4(e)) is accompanied by a temporary decrease in the ɛ graph (Figure 4(c)). This instability is the climax of the loop cycle in the ɛ,ɛ̇ and the ɛ̇,ɛ̈ graphs. As the loop cycle terminates, ɛ climbs to a higher value than the one it had when the cycle began, and the exponent continues to decay.In Figure4(d) we see that the effective mass fluctuates with a growing amplitude. As the cycle reaches towards its end, the effective mass continues to a higher amplitude and has its highest absolute values.


### 6.1.2. Second Cycle

The attractor's orbit continues to a second cycle, as can be seen in theɛ,ɛ̇ and ɛ̇,ɛ̈ graphs (Figures 5(a) and 5(b), resp.). This second cycle is again accompanied by a fluctuation (relative instability) in the Lyapunov exponent, as can be seen in Figure 5(e). The mass, as before, fluctuates with growing amplitude (we can observe the growing values on the y-axis in Figure 5(d) in comparison to Figure 4(d)). As the mass reaches a relatively low absolute value, the Lyapunov exponent rises to a relatively high one (the “jump” area of the exponent graph).(a) and (b) Theɛ,ɛ̇ and the ɛ̇,ɛ̈ phase space graphs, respectively. (c) and (d) ɛ and the (universal time) effective mass M(ɛ,ɛ̇) as functions of the number of iterations. As in the previous cycle, the mass fluctuation amplitude grows. (e) Again, the Lyapunov exponent fluctuates during the cycle and returns to decay as it ends.
(a)(b)(c)(d)(e)


### 6.1.3. Third Cycle

The consecutive cycle behaves in the same way as the two previous ones, as can be seen in Figures6(a)–6(e). The occurrence of the present loop cycle is seen when ɛ has a down spike (around the iteration value 2500 in Figure 6(c)). In this area the Lyapunov exponent has an up spike (Figure 6(e)). The effective mass, in contrast to the Lyapunov exponent, decreases during the cycle but rises to a higher amplitude as it ends (Figure 6(d)).(a) and (b) Theɛ,ɛ̇ and the ɛ̇,ɛ̈ phase space graphs, respectively. (c) and (d) ɛ and the (universal time) effective mass M(ɛ,ɛ̇) as functions of the number of iterations. As in the previous cycles, the mass fluctuation amplitude significantly grows towards the cycle’s end. One should also observe the relatively low value of the effective mass around the iteration value of 2500 (the climax of the cycle). (e) Again, the Lyapunov exponent fluctuates during the cycle and returns to decay afterwards.
(a)(b)(c)(d)(e)The above results suggest that as long as the electron “survives” the fluctuations (or cycles) and its orbit doesn't run away to infinity, the effective mass will reach higher absolute values in the end of each cycle. The effective mass, in general, fluctuates with increasing amplitude. As can be seen from the Lyapunov exponent graphs, the attractor is relatively unstable during the cycles but regains stability between them.A  correlation  is  clearly  observed  between  a  relatively  high  absolute  value  of  the  effective  mass  and  a  relatively  low  value  of  the  Lyapunov  exponent  (which  indicates  a  relative  stability).We also saw correlation between our attractor and nearby orbits, as seen in the time series associated with them, even in the vicinity of unity (the light cone). This correlation shows the absence of chaotic behavior.


## 6.1.1. First Cycle

In Figures4(a) and 4(b), we see ɛ,ɛ̇ and ɛ̇,ɛ̈ graphs of this cycle, respectively. Figures 4(c) and 4(d) show ɛ and the effective mass as functions of the number of iterations. Figure 4(e) shows the Lyapunov exponent. The orbit in Figure 4(a) enters a region of folding, which develops to a loop. In Figure 4(b), one can see an apparent attractor inducing motion beginning at the lower convex portion of the orbit which then turns back in a characteristic way to reach the last point visible on this graph.(a) and (b). Theɛ,ɛ̇ and the ɛ̇,ɛ̈ graphs, respectively. One can see the loop cycle in both phase space graphs. (c) and (d) ɛ and the (universal time) effective mass M(ɛ,ɛ̇) as functions of the number of iterations. One can see how ɛ temporarily decreases during the cycle. The effective mass significantly grows (in absolute value) as the cycle terminates. (e) The Lyapunov exponent fluctuates during the cycle (relative instability) but generally decreases.
(a)(b)(c)(d)(e)The instability in the Lyapunov exponent (Figure4(e)) is accompanied by a temporary decrease in the ɛ graph (Figure 4(c)). This instability is the climax of the loop cycle in the ɛ,ɛ̇ and the ɛ̇,ɛ̈ graphs. As the loop cycle terminates, ɛ climbs to a higher value than the one it had when the cycle began, and the exponent continues to decay.In Figure4(d) we see that the effective mass fluctuates with a growing amplitude. As the cycle reaches towards its end, the effective mass continues to a higher amplitude and has its highest absolute values.


## 6.1.2. Second Cycle

The attractor's orbit continues to a second cycle, as can be seen in theɛ,ɛ̇ and ɛ̇,ɛ̈ graphs (Figures 5(a) and 5(b), resp.). This second cycle is again accompanied by a fluctuation (relative instability) in the Lyapunov exponent, as can be seen in Figure 5(e). The mass, as before, fluctuates with growing amplitude (we can observe the growing values on the y-axis in Figure 5(d) in comparison to Figure 4(d)). As the mass reaches a relatively low absolute value, the Lyapunov exponent rises to a relatively high one (the “jump” area of the exponent graph).(a) and (b) Theɛ,ɛ̇ and the ɛ̇,ɛ̈ phase space graphs, respectively. (c) and (d) ɛ and the (universal time) effective mass M(ɛ,ɛ̇) as functions of the number of iterations. As in the previous cycle, the mass fluctuation amplitude grows. (e) Again, the Lyapunov exponent fluctuates during the cycle and returns to decay as it ends.
(a)(b)(c)(d)(e)


## 6.1.3. Third Cycle

The consecutive cycle behaves in the same way as the two previous ones, as can be seen in Figures6(a)–6(e). The occurrence of the present loop cycle is seen when ɛ has a down spike (around the iteration value 2500 in Figure 6(c)). In this area the Lyapunov exponent has an up spike (Figure 6(e)). The effective mass, in contrast to the Lyapunov exponent, decreases during the cycle but rises to a higher amplitude as it ends (Figure 6(d)).(a) and (b) Theɛ,ɛ̇ and the ɛ̇,ɛ̈ phase space graphs, respectively. (c) and (d) ɛ and the (universal time) effective mass M(ɛ,ɛ̇) as functions of the number of iterations. As in the previous cycles, the mass fluctuation amplitude significantly grows towards the cycle’s end. One should also observe the relatively low value of the effective mass around the iteration value of 2500 (the climax of the cycle). (e) Again, the Lyapunov exponent fluctuates during the cycle and returns to decay afterwards.
(a)(b)(c)(d)(e)The above results suggest that as long as the electron “survives” the fluctuations (or cycles) and its orbit doesn't run away to infinity, the effective mass will reach higher absolute values in the end of each cycle. The effective mass, in general, fluctuates with increasing amplitude. As can be seen from the Lyapunov exponent graphs, the attractor is relatively unstable during the cycles but regains stability between them.A  correlation  is  clearly  observed  between  a  relatively  high  absolute  value  of  the  effective  mass  and  a  relatively  low  value  of  the  Lyapunov  exponent  (which  indicates  a  relative  stability).We also saw correlation between our attractor and nearby orbits, as seen in the time series associated with them, even in the vicinity of unity (the light cone). This correlation shows the absence of chaotic behavior.


## 6.2. A Smooth (Less Physical) Attractor

In this subsection we give an example of a smooth attractor, as can be seen in Figures7(a)–7(e).(a) and (b) An attractor without fluctuations (cycles). (a)ɛ,ɛ̇ graph and (b) ɛ̇,ɛ̈ graph. (c) Lyapunov exponent and (d) ɛ itself, as functions of the (universal) time iterations. (e) The effective mass steadily grows (in absolute value) as ɛ moves towards unity.
(a)(b)(c)(d)(e)The only change in the initial conditions and parameters between this smooth orbit and the former, fluctuant, one is the cutoff frequency, which has a higher value now: 10-1 GeV-1 (instead of 10-5 GeV-1 in the fluctuant case). This orbit has no loop cycles in its phase space graphs. It has, however, a higher life span than the previous, fluctuant attractor, and thus it stays a longer time in the vicinity of unity. Nevertheless, also in this case the orbit will, at some point, run away to infinity.Although a higher cutoff frequency is less “physical”, the advantage of showing an attractor with relatively high cutoff frequency is that the fluctuations are removed, and the “pure” correlation between the growing effective mass and the rising stability (in the form of declining and even negative Lyapunov exponent) can be seen. Indeed, the attractor's Lyapunov exponent in Figure7(c) turns to be unambiguously negative and smooth, which may indicate the existence of a fixed point.The attractor'sɛ and effective mass as functions of the time iterations appear in Figures 7(d) and 7(e), respectively. There are no fluctuations. We can therefore conclude that the pattern is even clearer in the present example: as the attractor is steadily approaching the light cone, the effective mass is steadily growing (in absolute value), without fluctuations, and the Lyapunov exponent is negative and stable.


## 7. The Dynamical Equations

In this section, we will demonstrate some characteristics the orbit equation exhibits, which were previously found to be properties of the autonomous equation for the off-mass shell.As a first example, we choose a fluctuant, two-loop-cycle orbit (Figures8(a)–8(c), 9(a)–9(c), 10(a)–10(e)) with given initial conditions, and analyze it. In Figures 8(a)–8(c) we see (ẋ1, ẍ1) graphs of the orbit. An attractor can be clearly seen. The orbit reaches velocities which are very large (These velocities, in which ɛ is close to unity, correspond to very short segments of the world line due to the factor 1-ɛ connecting dτ and ds.)See also Figure9. Figures 8(a) and 9(a) a two-loop-cycle orbit, as shown in ẋ1,ẍ1 graph (Figure 8(a)) and ɛ,ɛ̇ graph (Figure 9(a)). Initial conditions: ẋ0=0.8, ẋ1=0.1, ẋ2=0.1, ẋ3=0.1, ẍ0=-10000, ẍ1=4, ẍ2=5, ẍ3=7,M/α=4.8·107  GeV2, and μ=10-5  GeV-1. Figures 8(b) and 9(b). The first loop cycle. ɛ has relatively low value (Figure 9(b)) when the value of ẋ1 is relatively high (Figure 8(b)). Figures 8(c) and 9(c). The second loop cycle, smaller than the first in the ɛ,ɛ̇ graph and larger than the first in the ẋ1,ẍ1 graph.
(a)(b)(c)See the caption of Figure8 which applies to both Figures 8 and 9.
(a)(b)(c)(a)ɛ during the first loop cycle. (b) The ẋ1 velocity during the first loop cycle. (c) The ẍ1 acceleration during the first loop cycle. (d) The effective mass M(ɛ,ɛ̇) during the first loop cycle. (e) The Lyapunov exponent during the first loop cycle.
(a)(b)(c)(d)(e)Figures9(a)–9(c) show the (ɛ,ɛ̇) graphs of the orbit. While the second loop cycle in the (ɛ,ɛ̇) graph becomes smaller than the first and closer to unity, the second loop cycle in the velocity-acceleration graph becomes larger than the first one.The fact that the second loop cycle in the velocity-acceleration graph is larger than the first means that during the attraction towards the light cone the velocity reaches more extreme values. One can look at this with connection to the behavior shown in the autonomous equation case, in which each loop cycle causes relative instability.Figures10(a)–10(e) show ɛ, the ẋ1 component of the velocity, the ẍ1 component of the acceleration, the effective mass M(ɛ,ɛ̇), and the Lyapunov exponent, respectively, as functions of the time iterations during the first loop cycle (Figures 8(b) and 9(b)). When the fluctuation occurs (around the iteration value of 2200), the Lyapunov exponent rises (relative instability). This is accompanied by a sudden rise in the velocity and the acceleration. However, the acceleration suddenly and very sharply jumps down and becomes negative. The powerful force which is responsible for the acceleration (brutal) change is the generalized Lorentz force. It is exactly at this point that the effective mass starts to grow (in absolute value), in order to reach a higher amplitude in its sinusoidal orbit. The effective mass is now growing, and at the same time the acceleration is decaying. The acceleration becomes less negative but slows the electron down, and the fluctuation terminates. Along with the slowing down, stability is restored (i.e., the Lyapunov exponent decays again).To repeat, the effective mass markedly grows (in absolute value) immediately after the sharp change of the acceleration, which is caused by the generalized Lorentz force. Actually, the effective mass increases its fluctuation amplitude. As the mass grows the acceleration decays, and stability is restored.The  growth  of  the  effective  mass  should  therefore  be  regarded  as  a  “stabilizing  reaction”  of  the  electron  against  the  generalized  Lorentz  force  laid  upon  it. The effective mass generated in the vicinity of unity (the light cone) is very large (The effective mass value at ɛ≈0.999 is about 300 times greater than its value at ɛ≈0.84.)As a second example, of which the motivation for presentation was discussed in the former section, we present a smooth, nonfluctuant orbit. This orbit has the same initial conditions and parameters as the former, fluctuant one. However, its cutoff frequency is higher: 10-1 GeV-1. The higher cutoff frequency removes the fluctuations and shows a clear connection between the effective mass growth and the decay of the velocity and the acceleration. Figures 11(a)–11(e) present ɛ, the acceleration, the velocity, Lyapunov exponent, and the effective mass of this orbit, respectively. The effective mass (absolute value) growth (Figure 11(e)) is accompanied by decaying acceleration and velocity (Figures 11(b) and 11(c), resp.) and growing stability (Figure 11(d)).(a) and (b)ɛ and the acceleration component ẍ1 as functions of the number of time iterations. (c) and (d) The velocity component ẋ1 and the Lyapunov exponent as functions of the number of time iterations. (e) The effective (universal time) mass.
(a)(b)(c)(d)(e)The growth of effective mass and stability on the one hand and the decay of acceleration and velocity on the other hand go together. This, again, should be regarded as a stabilizing effect, which occurs also near the light cone.


## 8. An External Force Added

When an external force of the form(8.1)e0fext5μ=Acos⁡(ωτ)
is added, the orbit equation (3.22) is(8.2)aμ(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3,τ)=f13μ(ẋνx⃛νẋμ+(1-ɛ)x⃛μ)+g3μx⃛μ+h3ẋνx⃛ν3μẋμ.This time we obtain the initial conditionsx⃛μ as a function of the components of the first and second derivatives and also of τ explicitly. The expressions for x⃛μ have the same form as the τ-independent solutions, with the only difference being in replacing aμ from (5.2) with(8.3)aμ(ẋ0,ẋ1,ẋ2,ẋ3,ẍ0,ẍ1,ẍ2,ẍ3,τ)≡(M⋅ẍμ-Acos⁡(ωτ))⋅λ(2π)32e02+12((1-ɛ)ẍμ+ɛ̇2ẋμ)(f12μ2-ɛ̇f1′2μ)-g13μ3ẋμ+g22μ2ẍμ+h1ɛ2μ2ẋμ-h2ɛμẍμ-h3ẍνẍν4μẋμ-h4ɛ̇28μẋμ.Again, the values forẋμ and ẍμ were chosen, and the values for x⃛μ were calculated using the orbit equation with the added perturbation (8.2). We used the initial values for ẋμ,ẍμ,x⃛μ built in this way to program a Runge-Kutta algorithm as before. The program's output is a trajectory of a self-interaction attractor, as a function of the external perturbation. In Figure 12(a), we present an (ɛ,ɛ̇) graph of the attractor's pure self-interaction, with the external perturbation turned off.(a) The orbit without an external field. Initial conditions:ẋ0=3.2, ẋ1=0.4, ẋ2=0.2, ẋ3=0.00049, ẍ0=2.83, ẍ1=0.3, ẍ2=0.1, ẍ3=-0.0011, M/α=105  GeV2, and μ=10-5  GeV-1. (b) and (c) The external perturbation added. The edges of the cycling polygonal orbit (b) form an oval frame, which can be seen as the lines connecting the coordinates are removed (c). Time step: 10-9 GeV-1. (d) The propagation of the perturbate attractor after many time iterations is in the path of the pure self-interaction attractor of (a).
(a)(b)(c)(d)In Figure12(b), the external perturbation, which we chose to be(8.4)e0fext5μ=1.6⋅10-9cos⁡(1014τ),
was turned on. The time step in this figure is of order 10-9 GeV-1. The figure shows a cycling polygonal orbit whose edges form an oval frame. The oval frame can be seen more clearly in Figure 12(c), where we removed the orbit lines which connect the coordinates.As a large number of iterations are done, the polygonal orbit propagates in the path of the unperturbed attractor (Figure12(d)).When we use a different time step which is ten times smaller than the one in Figures12(b)–12(d), we see again a cycling orbit whose edges form an oval (Figure 13(a)), but this time the oval frame is “filled” with a different polygonal pattern and looks like a “stretched” version of the previous one. As for the scale in Figure 13(a), it is about 12 times smaller than the one in Figures 12(b)–12(d). As before, there is a sinusoidal orbit when more iterations are taken (Figure 13(b)) and the propagation is in the path of the pure self-interaction attractor of Figure 12(a).The smaller scale polygonal orbit and its propagation. One can observe the cycling polygonal orbit “filling” the oval frame in a star-like pattern (a). There is also a sinusoidal propagation (b). Time step: 10-10 GeV-1. The scale is about 12 times smaller than the one in Figures 12(b)–12(d).
(a)(b)Two more similar reductions of the time step result in smaller-scale oval patterns, which are similar to the larger ones (Figures14(a), 14(b) and 15(a), 15(b)). Here, again, the sinusoidal nature appears. This self-similarity on smaller time steps (and smaller scales) is a fractal feature and hence may suggest the existence of a strange attractor.The cycling orbit filling the oval frame and the beginning of its propagation. Time step: 10-11 GeV-1. The scale is about 6 times smaller than the one in Figures 13(a) and 13(b).
(a)(b)Time step: 10-12 GeV-1. The scale is about 13 times smaller than the one in Figures 14(a) and 14(b).
(a)(b)


## 9. Summary and Discussion

We have shown that a generalized Lorentz force can be derived directly from the Hamiltonian for the classical relativistic motion of a charged particle in interaction with a generalized electromagnetic field by means of the Hamilton equations.We computed the charged particle self-interaction and obtained differential equations for its motion by using the Green's function for the electromagnetic fields, and by taking the motion of the particle as the current source. We studied the motion both in the absence and in the presence of an external field. We introduced the quantity(9.1)ɛ=1+ẋμẋμ,
which measures the deviation of the charged particle from its “mass shell”.


### 9.1. In the Absence of an External Field

In the absence of an external field, we obtained an autonomous equation for the development ofɛ. This equation describes the mass deviation which is caused by pure self-interaction. Alongside with the autonomous equation, we also derived the four-component equation of motion (which we termed the orbit equation). Their compatibility was demonstrated by solving them numerically using consistent initial values.Both the autonomous equation and the orbit equation show the existence of an attractor. We investigated the accumulation area of the attractor—the vicinity to the phase space coordinate(ɛ,ɛ̇,ɛ̈)=(1,0,0). We demonstrated that a lower and more “physical” cutoff frequency results in fluctuations in the orbit, as opposed to a higher (less physical) cutoff frequency, which results in a smooth orbit. Each fluctuation (or cycle, as it is seen in a phase space graph), causes a temporary instability, and a fluctuating attractor has shorter life span than a nonfluctuating one.A lower cutoff frequency leads toM(ɛ,ɛ̇) fluctuating with increasing amplitude as the orbit propagates towards the light cone. Within each such fluctuation the Lyapunov exponent jumps when the effective mass is relatively close to zero and decays when the absolute value of the effective mass increases. This phenomenon was demonstrated in three consecutive cycles.A higher cutoff frequency removes the cycles/fluctuations and causes the absolute value of the effective mass to grow with increasing rate asɛ propagates towards the light cone. In this scenario, the Lyapunov exponent decays steadily and without jumps (Figure 11(d)) and may become negative (Figure 7(c)).In both cases mentioned above, a very large effective mass is generated asɛ≈1 (Figures 3(f), 7(e), and 11(e)). This feature is possibly associated with effective nonzero size.The numerical solution of the orbit equations shows that the velocities and accelerations attained by the electron become more extreme with each fluctuation asɛ propagates towards unity. The acceleration, which is shown in Figure 10(c), increases at the beginning of a fluctuation and then suddenly and forcefully changes direction and reduces the velocity of the electron (Figure 10(b)). This powerful and sudden change is caused by the generalized Lorentz force. The slowing down of the electron is accompanied by a markedly growing effective mass, which stabilizes the electron against that force  (Figures 10(d) and 10(e)). The velocity decay along with the growing effective mass and stability is also seen in the nonfluctuating case (Figures 11(a)–11(e)).We also checked whether the attractor has chaotic nature, using a Lyapunov exponent calculation and a time series analysis. The results did not indicate chaotic nature.


### 9.2. In the Presence of External Field

In the presence of an external field, the attractor acquires a remarkable stochasticity which may suggest fractal nature. The propagation nature of the attractor changes from a flat line propagation in the pure self-interaction case to a spring like, sinusoidal propagation of a cycling polygon in the external perturbation case (Figures12, 13, and 14). However, the general path of the propagation remains the path of the pure self-interaction propagation (Figure 12(d)). A polygonal pattern also appears as smaller time steps are taken and the scales are smaller. This feature of self-similarity in smaller scales is of fractal nature and may suggest the existence of a strange attractor. Its stochastic type signals may be found in experiments with high-frequency detectors.


## 9.1. In the Absence of an External Field

In the absence of an external field, we obtained an autonomous equation for the development ofɛ. This equation describes the mass deviation which is caused by pure self-interaction. Alongside with the autonomous equation, we also derived the four-component equation of motion (which we termed the orbit equation). Their compatibility was demonstrated by solving them numerically using consistent initial values.Both the autonomous equation and the orbit equation show the existence of an attractor. We investigated the accumulation area of the attractor—the vicinity to the phase space coordinate(ɛ,ɛ̇,ɛ̈)=(1,0,0). We demonstrated that a lower and more “physical” cutoff frequency results in fluctuations in the orbit, as opposed to a higher (less physical) cutoff frequency, which results in a smooth orbit. Each fluctuation (or cycle, as it is seen in a phase space graph), causes a temporary instability, and a fluctuating attractor has shorter life span than a nonfluctuating one.A lower cutoff frequency leads toM(ɛ,ɛ̇) fluctuating with increasing amplitude as the orbit propagates towards the light cone. Within each such fluctuation the Lyapunov exponent jumps when the effective mass is relatively close to zero and decays when the absolute value of the effective mass increases. This phenomenon was demonstrated in three consecutive cycles.A higher cutoff frequency removes the cycles/fluctuations and causes the absolute value of the effective mass to grow with increasing rate asɛ propagates towards the light cone. In this scenario, the Lyapunov exponent decays steadily and without jumps (Figure 11(d)) and may become negative (Figure 7(c)).In both cases mentioned above, a very large effective mass is generated asɛ≈1 (Figures 3(f), 7(e), and 11(e)). This feature is possibly associated with effective nonzero size.The numerical solution of the orbit equations shows that the velocities and accelerations attained by the electron become more extreme with each fluctuation asɛ propagates towards unity. The acceleration, which is shown in Figure 10(c), increases at the beginning of a fluctuation and then suddenly and forcefully changes direction and reduces the velocity of the electron (Figure 10(b)). This powerful and sudden change is caused by the generalized Lorentz force. The slowing down of the electron is accompanied by a markedly growing effective mass, which stabilizes the electron against that force  (Figures 10(d) and 10(e)). The velocity decay along with the growing effective mass and stability is also seen in the nonfluctuating case (Figures 11(a)–11(e)).We also checked whether the attractor has chaotic nature, using a Lyapunov exponent calculation and a time series analysis. The results did not indicate chaotic nature.


## 9.2. In the Presence of External Field

In the presence of an external field, the attractor acquires a remarkable stochasticity which may suggest fractal nature. The propagation nature of the attractor changes from a flat line propagation in the pure self-interaction case to a spring like, sinusoidal propagation of a cycling polygon in the external perturbation case (Figures12, 13, and 14). However, the general path of the propagation remains the path of the pure self-interaction propagation (Figure 12(d)). A polygonal pattern also appears as smaller time steps are taken and the scales are smaller. This feature of self-similarity in smaller scales is of fractal nature and may suggest the existence of a strange attractor. Its stochastic type signals may be found in experiments with high-frequency detectors.



---

*Source: 602784-2010-12-09.xml*