INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL

J. Phys. A: Math. Gen. 34 (2001) 34213428 www.iop.org/Journals/ja PII: S0305-4470(01)14888-2

Relativistic integro-differential form of the
LorentzDirac equation in 3D without runaways
Michael Ibison and Harold E Puthoff
Institute for Advanced Studies at Austin, 4030 Braker Lane West, Suite 300 Austin, TX 78759,
USA
E-mail: [email protected] and [email protected]
Received 13 June 2000, in nal form 19 December 2000
Abstract
It is well known that the third-order LorentzDirac equation admits runaway
solutions wherein the energy of the particle grows without limit, even when
there is no external force. These solutions can be denied simply on physical
grounds, and on the basis of careful analysis of the correspondence between
classical and quantum theory. Nonetheless, one would prefer an equation that
did not admit unphysical behaviour at the outset. Such an equationan integro-
differential version of the LorentzDirac equationis currently available either
in only one dimension or in three dimensions (3D) in the non-relativistic limit.
It is shown herein how the LorentzDirac equation may be integrated
without approximation, and is thereby converted to a second-order integro-
differential equation in 3D satisfying the above requirement, i.e. as a result, no
additional constraints on the solutions are required because runaway solutions
are intrinsically absent. The derivation is placed within the historical context
established by standard works on classical electrodynamics by Rohrlich, and
by Jackson.
PACS numbers: 0320, 0350D, 0365B, 1460C
1. Introduction
The LorentzDirac equation (LDE) describes the motion of a classical charged particle subject
to both an external force and self-interaction due to radiation. An undesirable characteristic is
the prediction of an exponential (runaway) acceleration in the absence of an applied force. The
source of the trouble may be traced to the third-order derivative with respect to time. Since
one would prefer a second-order equation anyhow, a natural approach is to convert the original
LDE into a second-order equation by integrating over time. At the same time, one might
take the opportunity to eliminate the runaway solution by a suitable choice for the constant of
integration.
This is the method cited by Jackson [1], as it applies to a non-relativistic (and thereby
linearized) version of the LDE. It is successful in that runaway solutions are absent. The
0305-4470/01/163421+08$30.00 2001 IOP Publishing Ltd Printed in the UK 3421
3422 M Ibison and H E Puthoff
same approach was employed by Rohrlich [2] to the relativistic LDE, but without success; his
resulting equation still permits runaway solutions. The attempt failed because he was unable to
combine the acceleration and radiation parts (times an integrating factor) as a total differential
in proper time.
Jackson and Rohrlich are referred to herein because they are standard texts on classical
theory. However, for an earlier review of the subject that is both lucid and thorough, the reader
is referred to Erber [3]. The rst appearance of the non-relativistic integro-differential form
of the LDE is due to Haag [4] (and subsequently, in English, Plass [5]).
It has been shown that the non-relativistic integro-differential form of the LDE is the
nite-point limit of a nite-size (non-relativistic) model of the electron [68]. Since the latter
is free of runaway solutions, this may be regarded as evidence in favour of the validity of
the integro-differential form, over the original LDE. Also (very importantly) Sharp [9] has
shown that the non-relativistic integro-differential LDE corresponds to the quantum theory
of a non-relativistic charge coupled to the quantized electromagnetic eld (neither of which,
therefore, display runaway solutions). Both of these results point to the need for a relativistic
generalization of the existing non-relativistic integro-differential version of the LDE.
Barut [10] has proposed a method to eliminate both the runaway and pre-acceleration
behaviour of the LDE by making the Abraham vector disappear when the external eld
disappears. However, as pointed out by Blanco [11], such an equation is essentially quite
different fromthe original LDE. Jimenez and Hirsch [12] suggest that the non-relativistic LDE
be supplemented by an external, stochastic, electromagnetic eld, in the spirit of stochastic
electrodynamics (see, for instance, [13])). This, they argue, has the effect of eliminating the
undesirable runaway behaviour without modication of the LDE (to an integro-differential
form). Their programme, though promising, potentially suffers from an externally induced
runaway problem unless the stochastic eld is band-limited (which would be similar to
supposing a nite-sized charge).
2. Runaway solutions of the LorentzDirac equation
The LorentzDirac equation in proper time is [1, 2]
m
0
a m
0

0
_
da
d
+ a
2
u
_
= f () (1)
where the force f can depend on explicitly, and implicitly via the position and its derivatives.
a {a

} {a
0
, a} is the proper acceleration, u {u

} {u
0
, u} = {, u} is the proper
velocity, and a
2
= (a
0
)
2
a a, c = 1 and
0
= e
2
/6
0
m
0
is two-thirds the time it takes
for light to travel across the classical electron radius. The notorious runaway solution is most
easily demonstrated in one dimension, wherein the LDE is easily linearized [2]. With the
substitution dx/d = sinh (w()), one obtains from equation (1):
w
0
w = f/m
0
(2)
where f is the ordinary Newton force in the x direction. It is clear that even when there is no
external force, f = 0, w may increase without limit, since w exp (/
0
) is a solution. This
causes dx/d and to increase without limit, giving rise to the interpretation that the particle
has accelerated to the speed of light and has acquired an innite kinetic energy. The reason
for the presence of such solutions may be traced to the intrinsically non-conservative nature of
the equation of motion. It was conceived to account for losses due to radiation, but turns out
to admit gains, presumably by the same mechanism.
Relativistic integro-differential form of the LorentzDirac equation in 3D without runaways 3423
3. The non-relativistic integro-differential equation
When the velocities are small compared with c, 1, d dt , and equation (1) becomes
m
0
a m
0

0
da
dt
= f. (3)
(This non-relativistic formof the LDE is also called the AbrahamLorentz equation.) It suffers
fromthe same runaway solution as equation (2)the relativistic one-dimensional result written
in hyperbolic coordinates. The traditional remedy [1] is to replace equation (3) with the integro-
differential equation
m
0
a =
_

0
ds e
s
f (t + s
0
) (4)
It is readily veried upon substitution that the x that solves this equation is a subset of
those that solve equation (3). It is also clear that, provided f vanishes in the remote future,
the acceleration also vanishes in the remote future. Not only does this prescription eliminate
the runaway solution, but it also restores the boundary condition requirements to those of
a second-order differential equation, e.g. the position and velocity are given at some time.
This time need not be when the force is zero (i.e. the remote past or the remote future); it
may be any time. Though the runaway behaviour is tamed, it is at the expense of an acausal
connection between the applied force and the resulting acceleration. Specically, it is seen
from equation (4) that the acceleration depends on future forces (exhibits pre-acceleration).
However, the temporal range,
0
, of that dependence is such that pre-acceleration is too small
to be observed on classical time scales.
4. Rohrlichs relativistic integro-differential equation
It is carefully argued by Rohrlich [2] that runaway solutions must be denied by imposing a
suitable constraint, i.e. a boundary condition on the acceleration. In this paper, we will be
content with the condition
lim
+
d
2
x

d
2
= 0 lim
t +
d
2
x

dt
2
= 0 lim
t +
d
2
x
dt
2
= 0 (5)
since we require an acceptable prediction of future behaviour based on some initial condition,
given at some nominal but nite time. With the aim of integrating the constraint into the
equation of motion, Rohrlich investigates a formal integration of equation (1),
a

= A

e
/
0
+
_

e
(

)/
0
_
1
m
0

0
f

) + a
2
(

)u

)
_
(6)
where A

is a 4-vector constant of integration. He sets A

= 0, and considers the new

equation as a possible replacement for equation (1). However, as he points out, setting A

to
zero guarantees only that lim

e
/
0
a

= 0 which, clearly, is weaker than the requirement

that the acceleration vanish, equation (5). Therefore, we conclude that equation (6) with
A

= 0 is unsatisfactory, since a supplemental constraint must still be imposed to lter out

the unphysical behaviour.
5. An integrating factor for the LorentzDirac equation
A fully relativistic integro-differential form of the LorentzDirac equation that does not admit
runaway solutions (and therefore does not require supplemental constraints) is possible if a
3424 M Ibison and H E Puthoff
suitable integrating factor for the original LDE can be found. If it exists, an integrating factor
S
_
S

()
_
satisfying
d
d
(Sa) =
1

0
m
0
Sf (7)
will permitvia the integration of equation (7)the imposition of boundary conditions
equation (5) on the acceleration. For this integrating factor to exist, by carrying out the
differentiation in equation (7) and comparing with equation (1) left multiplied by S, it must be
true that
_
dS
d
+
S

0
_
a = a
2
Su (8)
where none of the elements of S can depend on the acceleration a. A substitution into
equation (8) of
S = Re
/
0
(9)
where R
_
R

()
_
, removes the exponential decay factor to give the requirement that R
satises
dR
d
a = a
2
Ru. (10)
There are only three independent equations in equation (1) because the product of both sides
with the 4-velocity is identically zero. Consequently, for any b

(), R

= b

()u

sets each
side of equation (7) to zero, and so cannot be a candidate for the integrating factor. It follows
that R cannot have a unique solution, since any candidate solution R

= C

(say) will
generate a family of solutions just by addition of this null solution: R

= C

+ b

()u

.
Of course, whatever form is chosen, that choice cannot impact the equation of motion for each
component of x

.
With the sign convention {q

} {q
0
, q}, a particularly simple solution of equation (1)
for the integrating factor is
R =
_
R

_
=
_
_
_
_
u
0
u
1
u
2
u
3
u
1
u
0
0 0
u
2
0 u
0
0
u
3
0 0 u
0
_
_
_
_
=
_
c

+ u
0

_
(11)
where {c

} {1, 0, 0, 0} is a unit timelike vector. With this denition, one easily sees that
equation (10) is satised, and, in particular, that the two terms are
dR

d
a

= a
2
R

= a
2
c

. (12)
Recalling equation (9), it follows that the LorentzDirac equation, equation (1), may be written
as
d
d
_
e
/
0
Ra
_
=
e
/
0

0
m
0
Rf (13)
where R is given by equation (11), and the inverse of R, denoted here by

R, is
R
1

_

R

_
=
1
u
0
_
_
_
_
u
2
0
u
0
u
1
u
0
u
2
u
0
u
3
u
0
u
1
u
2
1
+ 1 u
1
u
2
u
1
u
3
u
0
u
2
u
1
u
2
u
2
2
+ 1 u
2
u
3
u
0
u
3
u
1
u
3
u
2
u
3
u
2
3
+ 1
_
_
_
_
=
_
1
u
0
_

_
+ 2u

_
. (14)
Relativistic integro-differential form of the LorentzDirac equation in 3D without runaways 3425
R does not behave like a tensor under boosts, and is therefore not a Lorentz tensor. However, it
does behave like a tensor under spatial rotations and space and time translations, and is therefore
a Euclidean tensor. Nonetheless, the Lorentz invariance of the LorentzDirac equation is
preserved. This can be seen more readily if equation (13) is written as
m
0
a m
0

0
_
da
d
+

R
dR
d
a
_
= f (15)
whereupon it is apparent that the requirement is not that R should be a Lorentz tensor, but that

R
dR
d
a should be a true 4-vector. The latter is guaranteed by design. Specically it is equal to
a
2
u, in conformity with equation (1), as may be conrmed using equations (11) and (14).
6. Integration and imposition of the boundary condition
Formally, the rst integral of equation (13) is
e
/
0
R()a() e

c
/
0
R(
c
)a(
c
) =
1

0
m
0
_

c
d

/
0
R(

)f (

a() = e
(
c
)/
0
R
1
()R(
c
)a(
c
)
1

0
m
0
_

c
d

e
(

)/
0
R
1
()R(

)f (

)
(16)
where
c
is the time at which the proper acceleration is presumed known. We are now in a
position to impose the requirement that the acceleration in the remote future,
c
= +when
the force has long since vanishedis zero. With a(
c
) = 0, equation (16) becomes
a() =
1

0
m
0
_

e
(

)/
0
R
1
()R(

)f (

). (17)
Upon the change of variable s = (

)/
0
, this is
m
0
a() =
_

0
ds e
s
R
1
()R ( + s
0
) f ( + s
0
) (18)
which may be recognized as a relativistic version of the non-relativistic form, equation (4). It
is easily seen that, having isolated the second derivative on the left-hand side, the acceleration
is guaranteed to vanish in the remote future if the force also vanishes then. Therefore, the
solution is evidently free of runaways. Furthermore, it is evident that solutions of this equation
are a subset of the solutions of the original LorentzDirac equation, equation (1). Therefore, it
can be concluded that the integro-differential equation (18) is the physically correct equation
of motion for a classical charged particle; it retains the properties of the original LorentzDirac
equation without the unphysical behaviour.
Since it is not immediately evident from equation (18), here we conrm that, as required,
the acceleration is orthogonal to the velocity. Taking the 4-vector product of equation (18)
with the velocity gives
u

()a

() =
_

0
ds e
s
u

()

R

()R

( + s
0
) f

( + s
0
). (19)
Using equation (14) one nds that
u

= u

_
1
u
0
_

_
+ 2u

_
= c

. (20)
3426 M Ibison and H E Puthoff
Inserting this into equation (19) and then using equation (11) gives
u

() a

() =
_

0
ds e
s
R

0
( + s
0
) f

( + s
0
)
=
_

0
ds e
s
u

( + s
0
) f

( + s
0
) = 0 (21)
where the last step follows because the 4-force is required to be orthogonal to the velocity.
7. Proper-time vector form
The 3-vector formof equation (18) is obtained as follows. Given f

= {u f, u
0
f}, where f is
the ordinary Newton force vector (i.e. borrowed from dp/dt = f), and u
0
= =

1 + u u,
then, using equation (11), one obtains
R

=
_
c

+ u
0

_
f

= u

f
0
+ u
0
f

=
_
0,
_
uu
T
u
2
0
_
f
_
= {0, u (u f) f}. (22)
Denoting the 3-space part by w u (u f) f, equation (18) can be written as
m
0
a() = sub
33
_
R
1
()
_
_

0
ds e
s
w( + s
0
) (23)
where a is the proper acceleration, and where the sub operation extracts the 3 3 (spatial)
sub-matrix. Using equation (14) the latter is easily seen to be
sub
33
_
R
1
_
=
1
u
0
_
1 + uu
T
_
(24)
whereupon equation (23) gives the integro-differential version of the LDEin proper-time vector
form:
m
0
a =
1
_
1 + uu
T
_
_

0
ds e
s
_

2
uu
T
_
f
=
1
_
1 + uu
T
_
_

0
ds e
s
(f u (u f)) (25)
where the functions in the integrand are to be evaluated at + s
0
. In particular, if f is the
Lorentz force, f = e(E + u B/ ), then the proper acceleration is
m
0
a = e
1
_
1 + uu
T
_
_

0
ds e
s
(E u (u E) + u B). (26)
To write the proper acceleration in terms of vector cross-products, it is useful to dene an
intermediate quantity

f
_

0
ds e
s
(f u (u f)) (27)
where once again the functions in the integrand are to be evaluated at + s
0
. With this
substitution, an alternative form for equation (26) is therefore
m
0
a =

f + u
_
u

f
_
/. (28)
Relativistic integro-differential form of the LorentzDirac equation in 3D without runaways 3427
8. Proper-time series expansion in
0
A series expansion of the integrand in ascending powers of
0
can be expected to converge
rapidly if the projection of the force,
_

2
uu
T
_
f, is slowly varying on the time scale of the
classical time
0
. From equation (18), one has
m
0
a = R
1

n=0
_

0
d
d
_
n
(Rf ) (29)
where all functions are now evaluated at time . In vector form this is
m
0
a =
1
_
1 + uu
T
_

n=0
_

0
d
d
_
n
_

2
uu
T
_
f. (30)
9. Ordinary-time vector form
The integro-differential form of the LDE can be cast as a 3-vector equation in ordinary time
as follows. From equation (17), one has
R

() =
1

0
m
0
_

e
(

)/
0
R

)f

) (31)
the left-hand side of which is
R

=
_
0, u
du
0
d
u
0
du
d
_
=
_
0, u
2
0
d (u/u
0
)
d
_
=
_
0,
3
d
dt
_
(32)
where = dx/dt is the ordinary velocity, i.e. the left-hand side of equation (31) is already in
the direction of the ordinary acceleration. Furthermore, noting that the product in the integrand
is
R

=
_
0,
2
(( f) f)
_
(33)
then substitution of equations (32) and (33) into equation (17) gives

=
1

0
m
0

3
_

e
(

)/
0

2
(f ( f)) =
1

0
m
0

3
_

t
dt

e
(

)/
0
H(t

) (34)
where the components of H(t

) = (f ( f)) are now redened as functions of ordinary

time. The transformation is complete once the exponential damping factor is explicitly cast as
a function of ordinary time:

=
1

0
m
0

3
_

t
dt

exp
__
t
t

dt

0
(t

)
_
H(t

)
=
1

0
m
0

3
_

0
dt

exp
_

_
t +t

t
dt

0
(t

)
_
H(t + t

). (35)
As for the proper-time form, the variable of integrationcanbe rendereddimensionless, although
here it does not result in a simplication. Letting s = t

/
0
:
d
dt
=
1
m
0

3
_

0
ds exp
_

_
t +
0
s
t
dt

0
(t

)
_
H(t +
0
s)
=
1
m
0

3
_

0
ds exp
_

_
s
0
ds

(t +
0
s

)
_
H(t +
0
s). (36)
If f is the Lorentz force then H = e (E ( E) + B).
3428 M Ibison and H E Puthoff
10. Ordinary-time series expansion in
0
An ordinary-time series expansion of the integrand in ascending powers of
0
can obtained
from equation (36) by integrating by parts. The result is
d
dt
=
1
m
0

n=0
_

0
d
dt
_
n
_

2
(f ( f))
_
(37)
where the functions are of ordinary time, evaluated at time t .
11. Summary
A physically acceptable relativistic equation of motion for a classical charged particle in three
spatial dimensions has been derived that has the properties desired of the original Lorentz
Dirac equation, but without the unphysical behaviour. The exclusion of runaway solutions has
been achieved by nding an integrating factor for the original LorentzDirac equation so that
the acceleration can be written as an integral operator on the force.
References
[1] Jackson J D 1975 Classical Electrodynamics (New York: Wiley) ch 17
[2] Rohrlich F 1965 Classical Charged Particles (Reading, MA: Addison-Wesley)
[3] Erber T 1961 Fortsch. Phys. 9 343
[4] Haag R 1955 Z. Naturf. A 10 752
[5] Plass G N 1961 Rev. Mod. Phys. 33 37
[6] Sorg M 1976 Z. Naturf. A 31 683
[7] Moniz E J and Sharp D H 1977 Phys. Rev. D 15 2850
[8] Levine H, Moniz E J and Sharp D H 1977 Am. J. Phys. 45 75
[9] Sharp D H 1980 Foundations of Radiation Theory and Quantum Electrodynamics ed A O Barut (New York:
Plenum) ch 10
[10] Barut A O 1990 Phys. Lett. A 145 387
[11] Blanco R 1992 Phys. Lett. A 169 115
[12] Jimenez J L and Hirsch J 1986 Nuovo Cimento B 98 87
[13] Boyer T 1980 Foundations of Radiation Theory and Quantum Electrodynamics ed A O Barut (New York:
Plenum) ch 5