Peeter Joot's (OLD) Blog.

Math, physics, perl, and programming obscurity.

PHY450H1S. Relativistic Electrodynamics Lecture 26 (Taught by Prof. Erich Poppitz). Radiation reaction force for a dipole system.

Posted by peeterjoot on April 11, 2011

[Click here for a PDF of this post with nicer formatting (especially if my latex to wordpress script has left FORMULA DOES NOT PARSE errors.)]


Covering chapter 8 section 65 material from the text [1].

Covering pp. 181-195: (182-189) [Tuesday, Mar. 29]; the EM potentials to order (v/c)^2 (190-193); the “Darwin Lagrangian. and Hamiltonian for a system of nonrelativistic charged particles to order (v/c)^2 and its many uses in physics (194-195) [Wednesday, Mar. 30]

Next week (last topic): attempt to go to the next order (v/c)^3 – radiation damping, the limitations of classical electrodynamics, and the relevant time/length/energy scales.


A system of N charged particles m_a, q_a ; a \in [1, N] closed system and nonrelativistic, v_a/c \ll 1. In this case we can incorporate EM effects in a Largrangian ONLY involving particles (EM field not a dynamical DOF). In general case, this works to O((v/c)^2), because at O((v/c)) system radiation effects occur.

In a specific case, when

\begin{aligned}\frac{m_1}{q_1} = \frac{m_2}{q_2} = \frac{m_3}{q_3} = \cdots\end{aligned} \hspace{\stretch{1}}(2.1)

we can do that (meaning use a Lagrangian with particles only) to O((v/c)^4) because of specific symmetries in such a system.

The Lagrangian for our particle after the gauge transformation is

\begin{aligned}\mathcal{L}_a = \frac{1}{{2}} m_a \mathbf{v}_a^2 + \frac{m_a}{8} \frac{\mathbf{v}_a^4}{c^2} -\sum_{b \ne a} \frac{q_a q_b}{{\left\lvert{\mathbf{x}_a(t) - \mathbf{x}_b(t)}\right\rvert}}+\sum_b q_a q_b \frac{\mathbf{v}_a \cdot \mathbf{v}_b + (\mathbf{n} \cdot \mathbf{v}_a) (\mathbf{n} \cdot \mathbf{v}_b)}{2 c^2 {\left\lvert{\mathbf{x} - \mathbf{x}_b}\right\rvert} }.\end{aligned} \hspace{\stretch{1}}(2.2)

Next time we’ll probably get to the Lagrangian for the entire system. It was hinted that this is called the Darwin Lagrangian (after Charles Darwin’s grandson).

We find for whole system

\begin{aligned}\mathcal{L} =  \sum_a \mathcal{L}_a + \frac{1}{{2}} \sum_a \mathcal{L}_a (interaction)\end{aligned} \hspace{\stretch{1}}(2.3)

\begin{aligned}\mathcal{L} = \frac{1}{{2}} \sum_a m_a \mathbf{v}_a^2 + \sum_a \frac{m_a}{8} \frac{\mathbf{v}_a^4}{c^2} -\sum_{ a < b} \frac{q_a q_b}{{\left\lvert{\mathbf{x}_a(t) - \mathbf{x}_b(t)}\right\rvert}}+\sum_b q_a q_b \frac{\mathbf{v}_a \cdot \mathbf{v}_b + (\mathbf{n} \cdot \mathbf{v}_a) (\mathbf{n} \cdot \mathbf{v}_b)}{2 c^2 {\left\lvert{\mathbf{x} - \mathbf{x}_b}\right\rvert} }.\end{aligned} \hspace{\stretch{1}}(2.4)

This is the Darwin Lagrangian (also Charles). The Darwin Hamiltonian is then

\begin{aligned}H = \sum_a \frac{p_a}{2 m_a} \mathbf{v}_a^2 + \sum_a \frac{p_a^4}{8 m_a^3 c^2} +\sum_{ a < b} \frac{q_a q_b}{{\left\lvert{\mathbf{x}_a(t) - \mathbf{x}_b(t)}\right\rvert}}- \sum_{ a < b } \frac{q_a q_b}{ 2 c^2 m_a m_b } \frac{\mathbf{p}_a \cdot \mathbf{p}_b + (\mathbf{n}_{a b} \cdot \mathbf{p}_a) (\mathbf{n}_{ a b } \cdot \mathbf{p}_b)}{{\left\lvert{\mathbf{x} - \mathbf{x}_b}\right\rvert} }.\end{aligned} \hspace{\stretch{1}}(2.5)

Incorporating radiation effects as a friction term.

To O((v/c)^3) obvious problem due to radiation (system not closed). We’ll incorporate radiation via a function term in the EOM

Again consider the dipole system

\begin{aligned}m \dot{d}{z} &= -k z \\ \omega^2 &= \frac{k}{m}\end{aligned} \hspace{\stretch{1}}(3.6)


\begin{aligned}m \dot{d}{z} = -\omega^2 m z\end{aligned} \hspace{\stretch{1}}(3.8)


\begin{aligned}\frac{d{{}}}{dt}\left( \frac{m}{2} \dot{z}^2 + \frac{ m \omega^2 }{2} z^2 \right) = 0 \end{aligned} \hspace{\stretch{1}}(3.9)

(because there’s no radiation).

The energy radiated per unit time averaged per period is

\begin{aligned}P = \frac{ 2 e^2 }{ 3 c^3} \left\langle{{ \dot{d}{z}^2 }}\right\rangle\end{aligned} \hspace{\stretch{1}}(3.10)

We’ll modify the EOM

\begin{aligned}m \dot{d}{z} = -\omega^2 m z + f_{\text{radiation}}\end{aligned} \hspace{\stretch{1}}(3.11)

Employing an integration factor \dot{z} we have

\begin{aligned}m \dot{d}{z} \dot{z} = -\omega^2 m z \dot{z} + f_{\text{radiation}} \dot{z}\end{aligned} \hspace{\stretch{1}}(3.12)


\begin{aligned}\frac{d{{}}}{dt} \left( m \dot{z}^2 + \omega^2 m z^2 \right) = f_{\text{radiation}} \dot{z}\end{aligned} \hspace{\stretch{1}}(3.13)

Observe that the last expression, force times velocity, has the form of power

\begin{aligned}m \frac{d^2 z}{dt^2} \frac{dz}{dt} = \frac{d{{}}}{dt} \left( \frac{m}{2} \left( \frac{dz}{dt} \right)^2 \right)\end{aligned} \hspace{\stretch{1}}(3.14)

So we can make an identification with the time rate of energy lost by the system due to radiation

\begin{aligned}\frac{d{{}}}{dt} \left( m \dot{z}^2 + \omega^2 m z^2 \right) \equiv \frac{d{{\mathcal{E}}}}{dt}.\end{aligned} \hspace{\stretch{1}}(3.15)

Average over period both sides

\begin{aligned}\left\langle{{ \frac{d{{\mathcal{E}}}}{dt} }}\right\rangle = \left\langle{{ f_{\text{radiation}} \dot{z} }}\right\rangle=- \frac{2 e^2 }{3 c^3} \left\langle{{\dot{d}{z}^2}}\right\rangle\end{aligned} \hspace{\stretch{1}}(3.16)

We demand this last equality, by requiring the energy change rate to equal that of the dipole power (but negative since it is a loss) that we previously calculated.


\begin{aligned}f_{\text{radiation}} = \frac{2 e^2 }{3 c^3} \dddot{z}\end{aligned} \hspace{\stretch{1}}(3.17)


We need to show

\begin{aligned}\left\langle{{ f_{\text{radiation}} }}\right\rangle= - \frac{2 e^2 }{3 c^3} \left\langle{{\dot{d}{z}^2}}\right\rangle\end{aligned} \hspace{\stretch{1}}(3.18)

We have

\begin{aligned}\frac{2 e^2}{3 c^3} \left\langle{{ \dddot{z} \dot{z} }}\right\rangle &= \frac{2 e^2}{3 c^3} \frac{1}{{T}} \int_0^T dt \dddot{z} \dot{z} \\ &= \frac{2 e^2}{3 c^3} \frac{1}{{T}} \int_0^T dt {\frac{d{{}}}{dt} (\dot{d}{z} \dot{z}) }-\frac{2 e^2}{3 c^3} \frac{1}{{T}} \int_0^T dt (\dot{d}{z})^2\end{aligned}

We first used (\dot{d}{z} \dot{z})' = \dddot{z} \dot{z} + (\dot{d}{z})^2. The first integral above is zero since the derivative of \dot{d}{z} \dot{z} = (-\omega^2 z_0 \sin\omega t)(\omega z_0 \cos\omega t) = -\omega^3 z_0^2 \sin(2 \omega t)/2 is also periodic, and vanishes when integrated over the interval.

\begin{aligned}\frac{2 e^2}{3 c^3} \left\langle{{ \dddot{z} \dot{z} }}\right\rangle =-\frac{2 e^2}{3 c^3} \left\langle{{ (\dot{d}{z})^2 }}\right\rangle\end{aligned} \hspace{\stretch{1}}(3.19)

We can therefore write

\begin{aligned}m \dot{d}{z} = -m \omega^2 z + \frac{2 e^2}{3 c^3} \dddot{z}\end{aligned} \hspace{\stretch{1}}(3.20)

Our “frictional” correction is the radiation reaction force proportional to the third derivative of the position.

Rearranging slightly, this is

\begin{aligned}\dot{d}{z} = - \omega^2 z + \frac{2}{3 c} \left( \frac{e^2}{m c^2} \right) \dddot{z} = - \omega^2 z + \frac{2}{3 c} \frac{r_e}{c} \dddot{z},\end{aligned} \hspace{\stretch{1}}(3.21)

where r_e \sim 10^{-13} \text{cm} is the “classical radius” of the electron. In our frictional term we have r_e/c, the time for light to cross the classical radius of the electron.

There are lots of problems with this. One of the easiest is with \omega = 0. Then we have

\begin{aligned}\dot{d}{z} = \frac{2}{3} \frac{r_e}{c} \dddot{z}\end{aligned} \hspace{\stretch{1}}(3.22)

with solution

\begin{aligned}z \sim e^{\alpha t},\end{aligned} \hspace{\stretch{1}}(3.23)


\begin{aligned}\alpha \sim \frac{c}{r_e} \sim \frac{1}{{\tau_e}}.\end{aligned} \hspace{\stretch{1}}(3.24)

This is a self accelerating system! Note that we can also get into this trouble with \omega \ne 0, but those examples are harder to find (see: [2]).

FIXME: borrow this text again to give that section a read.

The sensible point of view is that this third term (f_{\text{rad}}) should be taken seriously only if it is small compared to the first two terms.


[1] L.D. Landau and E.M. Lifshitz. The classical theory of fields. Butterworth-Heinemann, 1980.

[2] D.J. Griffith. Introduction to Electrodynamics. Prentice-Hall, 1981.


Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

%d bloggers like this: