advanced32 min studyLesson 260 in Physics

Use radiation reaction as a controlled approximation

Diagnose runaway solutions and derive a weak reduced-order self-force while keeping near-field energy in the balance.

On this page

What you’ll understand

  • Identify the additional dynamical freedom in the Abraham–Lorentz equation.
  • Derive a first-order reduction for a slowly changing external force.
  • Explain the roles of Schott energy, nonrunaway conditions, and approximation limits.

Radiated energy must come from work or from energy already stored in matter and fields. Turning this balance into a force on an ideal point charge is subtle because its near electromagnetic field is singular. The nonrelativistic Abraham–Lorentz equation displays the difficulty through an extra time derivative and runaway solutions. Reduction of order gives a useful weak-correction model when its small time scale is genuinely separated from the external dynamics; it is not an exact cure valid at arbitrary frequency or field strength.

A third derivative introduces an extra mode

Write the equation in terms of acceleration and its time derivative. The characteristic self-force time is proportional to charge squared divided by mass and speed cubed. Its SI value for an electron is about six times ten to the minus twenty-four seconds. The coefficient is tiny for many macroscopic trajectories, but a small coefficient multiplying the highest derivative is a singular perturbation: the differential equation gains an independent mode rather than merely a small adjustment to every solution.

With the external force set to zero, the acceleration equation admits exponential growth on this short time scale. Ordinary position and velocity initial data do not fix the extra acceleration degree of freedom. Rejecting growing solutions by a future boundary condition can select a bounded branch, but that branch generally responds slightly before a force changes. This preacceleration exposes the idealized point-particle equation's limitation. Choosing an arbitrary numerical initial acceleration is not a physical resolution of the runaway problem.

ma=Fext+mτ0a˙,τ0=q26πϵ0mc3,Fext=0 ⟹ a=a0et/τ0m\mathbf a=\mathbf F_{\mathrm{ext}}+m\tau_0\dot{\mathbf a},\qquad \tau_0=\frac{q^2}{6\pi\epsilon_0mc^3},\qquad \mathbf F_{\mathrm{ext}}=0\ \Longrightarrow\ \mathbf a=\mathbf a_0e^{t/\tau_0}
m is renormalized particle mass in kg, τ₀ is in seconds, and dots use nonrelativistic laboratory time. This is the ideal point-charge Abraham–Lorentz model outside quantum and relativistic regimes.↔ Scroll wide equations horizontally.

Reduce order within a small-parameter expansion

Suppose the external force is a smooth prescribed function of time varying on a scale much longer than the self-force time. To zeroth order acceleration is force divided by mass. Substitute that leading expression only into the already small self-force term and discard corrections of second order in the short time. The resulting acceleration depends on force and its derivative but has no independent exponentially growing mode. This is a perturbative approximation to the physical slow branch.

For forces depending on position and velocity, the derivative is a total derivative evaluated with leading equations of motion; it is not simply an explicit partial time derivative. The relativistic analogue is associated with the Landau–Lifshitz equation. Reduction of order does not authorize abrupt force steps, motion varying on the self-force time, or quantum recoil. Higher terms may form an asymptotic rather than convergent expansion. A bounded reduced solution is useful only where its neglected terms remain small.

a≃F(t)m+τ0mdF(t)dt,τ0T≪1;F=−mω02x ⟹ x¨+τ0ω02x˙+ω02x≃0\mathbf a\simeq\frac{\mathbf F(t)}m+\frac{\tau_0}m\frac{d\mathbf F(t)}{dt},\qquad \frac{\tau_0}{T}\ll1;\qquad F=-m\omega_0^2x\ \Longrightarrow\ \ddot x+\tau_0\omega_0^2\dot x+\omega_0^2x\simeq0
The first formula assumes force prescribed in time and drops terms of relative order (τ₀/T)². For the harmonic force the leading total derivative gives weak damping, valid when τ₀ω₀≪1; ω₀ is in rad/s.↔ Scroll wide equations horizontally.

Instantaneous self-work is not minus radiated power

Dot the Abraham–Lorentz self-force with velocity and use the product rule. One term is minus Larmor power, while another is a total derivative involving acceleration dotted into velocity. This latter term describes exchange with the bound near field, commonly organized as Schott energy. Omitting it incorrectly equates instantaneous mechanical energy loss with emitted power for arbitrary nonperiodic motion. Over a complete steady cycle its net change vanishes, making cycle-averaged work equal to minus average radiation.

A reduced harmonic model therefore resembles an ordinary weakly damped oscillator. The available oscillator lab has positive stiffness, unit initial displacement, and zero initial velocity. Its damping is phenomenological and adjustable; it does not derive a charge self-field or display runaway modes. Use it to compare the qualitative slow-branch decay with the undamped limit. Changing its damping control recomputes the whole motion at retained time, so independent comparisons require a fresh reset and both parameters specified.

Fself⋅v=mτ0ddt(a⋅v)−mτ0a2,ddt(K+ES)=Fext⋅v−P,ES=−mτ0a⋅v\mathbf F_{\mathrm{self}}\cdot\mathbf v=m\tau_0\frac{d}{dt}(\mathbf a\cdot\mathbf v)-m\tau_0a^2,\qquad \frac{d}{dt}(K+E_S)=\mathbf F_{\mathrm{ext}}\cdot\mathbf v-P,\quad E_S=-m\tau_0\mathbf a\cdot\mathbf v
K=mv²/2 is kinetic energy, P=mτ₀a² is Larmor power, and ES is bound-field Schott energy relative to the chosen bookkeeping reference. The identity is for the nonrelativistic Abraham–Lorentz model.↔ Scroll wide equations horizontally.

Reduce a smooth sinusoidal force

A prescribed one-dimensional force has F(t)/m=10 cos(ωt) m/s² and τ₀ω=0.02. Find reduced-order acceleration at ωt=0 and π/2, and estimate the relative order of omitted terms.

  1. Differentiate the leading force divided by mass: its derivative is −10ω sin(ωt). Multiplying by τ₀ gives the correction −0.20 sin(ωt) m/s².

  2. The reduced result is a≈10 cos(ωt)−0.20 sin(ωt). At phase zero it is 10 m/s²; at phase π/2 the leading force is zero but the derivative correction gives −0.20 m/s².

  3. The expansion parameter is 0.02, so nominal next terms are of relative order 0.02²=0.0004 compared with the characteristic 10 m/s² amplitude. Relative error must not be measured against the instantaneous leading force at its zero.

The result

The acceleration is 10 m/s² and approximately −0.20 m/s² at the two phases. The nominal neglected scale is about 0.004 m/s², subject to smooth forcing and the assumed classical regime.

Explore the model

Interactive lab

Interactive lab · oscillator

Compare a weakly damped slow branch

For each fresh run Reset and set BOTH damping ratio ζ=0, natural frequency ω₀=1; then Reset and set ζ=0.05, ω₀=1. Predict an undamped oscillation versus a slowly decaying envelope before Play. Both runs start x=1, v=0. The unit-mass oscillator illustrates damping, not a self-force calculation or the Abraham–Lorentz runaway branch.

Compare a weakly damped slow branchUnit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. Current readings: Displacement x(t): 1; Velocity v(t): 0; Regime: underdamped. All parameters are adjustable with the labeled controls below.02.557.510−1.15−0.5700.571.15Time · sDisplacementunderdamped
Displacement x(t)Current time
Displacement x(t)
1
Velocity v(t)
0
Regime
underdamped
t = 0 s

Unit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA reduced weak harmonic self-force behaves like damping; the underlying approximation still needs a short self-force time.

Check your understanding

Why is it inconsistent to solve the reduced-order approximation as though it retained the independent runaway acceleration mode?

Try it yourself

For a weak harmonic reduced model take ω₀=2 rad/s and τ₀=0.01 s. Find the velocity damping coefficient, the amplitude-envelope decay rate, and the small parameter.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Compare x¨+τ₀ω₀²x˙+ω₀²x=0 with x¨+2Γx˙+ω₀²x=0.

Compare with the solution

The velocity coefficient is τ₀ω₀²=0.04 s⁻¹. Thus Γ=0.02 s⁻¹ and the weakly damped amplitude envelope is exp(−0.02t), with a 50 s e-folding time. The expansion parameter τ₀ω₀=0.02 is small. These illustrative parameters do not represent an electron self-force time; a physical particle calculation must use its own charge and mass.

Take it with you

  • A small highest-derivative coefficient can introduce a rapidly growing extra mode.
  • Reduction of order is a weak smooth-dynamics approximation with a stated small parameter.
  • Bound-field energy exchange separates instantaneous self-work from radiated power.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.