# Advanced Classical Mechanics/The Eikonal Approximation and Classical Particle Motion

This wavepacket somewhat obeys the eikonal approximation because the potential energy, U(x), is nearly uniform over the length of the wavepacket. Moreover, the wavepacket is sufficiently large compared to the wavelength that dispersion is only moderate.

## History

This article appears in the "zeroth" (unrefereed) issue of WikiJournal of Science

In 1834, William Rowan Hamilton published a paper[1] in which he displayed an understanding between the mathematics particle motion and the propagation of linear waves in the eikonal approximation. Over a century later, Herbert Goldstein would remark that Hamilton would have postulated Schrodinger's equation had there been experimental evidence that particles were waves.[2] For most modern physicists, the connection is understood through the Ehrenfest theorem of quantum mechanics.

Hamilton's discovery was that wavepackets of all linear wave equations obey what we now call Hamiltonian equations of motion, provided the wavelength is sufficiently short compared with the distance over which the medium varies.

It should be noted that not only did Hamilton's original calculation differ significantly than as presented below, the methods used in the simple one-dimensional wave do not allow one to recover Hamilton's theory. With more than one spacial variable, it is necessary to perform a different calculation along the lines of Wikipedia:Eikonal_approximation

## The terms in a simple linear wave equation

Let us write a linear partial differential equation that is designed to illustrate the eikonal approximation to wavepacket motion:

${\displaystyle {\frac {\partial ^{2}\,\psi }{\partial t^{2}}}=a(x)\cdot {\frac {\partial ^{2}\psi }{\partial x^{2}}}+b(x)\cdot {\frac {\partial ^{2}\psi }{\partial x\partial t}}+c(x)\cdot \psi }$
The wake creation is nonlinear, but perhaps the linear approximation holds behind the ducks.

This PDE is linear in wave amplitude, ${\displaystyle \psi =\psi (x,t)}$, and the three coefficients (a, b, c) depend on position.[3]

The wave is linear because the terms only vary as the first power of amplitude, ${\displaystyle \psi =\psi (x,t)}$. In virtually all mechanical systems, there is some non-linearity in the wave equation, for example with term like ${\displaystyle d(x)\psi ^{2}}$. But at sufficiently low amplitude, ${\displaystyle d(x)\psi ^{2}< and the nonlinearity may be ignored. Not only are nonlinear equations more difficult to solve, their solutions lack certain properties that yield deep insights into the system.[4]

In the image to the right, linearity permits us to simplify the analysis, making the linear approximation of low amplitude waves. To the extent that this approximation holds, we may apply the superposition principle and calculate the wake of a single duck, ignoring the other other ducks' wakes as well as the gentle wave through which they are swimming.

### Homogeneous case

In general, the coefficients could depend on both position and time (x and t), but for simplicity we assume only a spatial dependence. These coefficients represent the medium that supports the wave, and the dependence on x indicates that this inhomogeneous in position (meaning that the medium changes as one moves along x). Let us begin with the completely homogeneous case, using the o-subscript to denote constants. This "trial solution",

${\displaystyle \psi (x,t)=A_{0}\exp \left(ik_{0}x-i\omega _{0}t\right)}$

substituted into,

${\displaystyle {\frac {\partial ^{2}\,\psi }{\partial t^{2}}}=a_{0}\cdot {\frac {\partial ^{2}\psi }{\partial x^{2}}}+b_{0}\cdot {\frac {\partial ^{2}\psi }{\partial x\partial t}}+c_{0}\cdot \psi }$

yields this polynomial for angular frequency, ${\displaystyle \omega _{0}}$, and wavenumber, ${\displaystyle k_{0}}$:

${\displaystyle \omega _{0}^{2}+k_{0}\omega _{0}b_{0}+\left(c_{0}-k_{0}^{2}a_{0}\right)=0}$

Our wave equation was chosen so that the polynomial is an easily solved quadratic equation that supports two waves or "branches":

${\displaystyle \omega _{0}=-{\frac {1}{2}}k_{0}b_{0}\pm {\sqrt {{\frac {1}{4}}k_{0}^{2}b_{0}^{2}+k_{0}^{2}a_{0}^{2}-c_{0}^{2}}}}$

### A dispersion relation for the weakly inhomogeneous case

If the system is weakly inhomogeneous, the coefficients (a, b, c) are nearly uniform over many wavelengths of a wavepacket. Hence, over the duration of the wavepacket this dispersion relation should be approximately true:

${\displaystyle \omega (k,x)=-{\frac {1}{2}}k(x)b(x)\pm }$${\displaystyle {\sqrt {{\frac {1}{4}}k(x)^{2}b(x)^{2}+k(x)^{2}a(x)^{2}-c(x)^{2}}}}$

Let ${\displaystyle L}$ represent a length over which these coefficients remain nearly unchanged, ${\displaystyle \lambda =2\pi /k}$ be the wavelength, and ${\displaystyle \Delta x}$ be the size of the wavepacket. The eikonal approximation requires that the wavepacket is large compared to one wavelength, yet small compared with the distance one must travel before the coefficients (a, b, c) change significantly:

${\displaystyle \Delta L>>\Delta x>>\lambda .}$

This inequality is the eikonal approximation. Also, we have assumed that the medium does not change with respect to time, which implies that the frequency, ω, does not change as it moves through the medium.[5][6] Thus, k must evolve in such a way that ω remains constant. From a well known relationship for differentials of functions of two or more variables, we have:

${\displaystyle \Delta \omega (x,k)=0={\frac {\partial \omega }{\partial x}}\Delta x+{\frac {\partial \omega }{\partial k}}\Delta k}$

Using the concept of group velocity, we have

${\displaystyle {\frac {\Delta x}{\Delta t}}={\frac {dx}{dt}}={\frac {\partial \omega }{\partial k}}}$,  and  ${\displaystyle {\frac {\Delta k}{\Delta t}}={\frac {dk}{dt}}=-{\frac {\partial \omega }{\partial x}}}$

Here we have replaced Δk, Δx, and Δt, by dk, dx, and dt, respectively.

## Application to the Schrödinger equation

According to De Broglie's matter wave hypothesis, energy is, E = ħω = V + p2/2m, where V=V(x) is potential energy, and momentum is, p = mv = ħk. Expressing angular frequency in terms of this expression for energy, we have:

${\displaystyle \hbar \omega (k,x)={\frac {\hbar ^{2}k^{2}}{2m}}+V(x)}$

It is straightforward to recover Newtonian mechanics from this. For example, the Newtonian force is, F = –∂V/∂x, and

${\displaystyle {\frac {\partial \omega }{\partial k}}={\frac {\hbar k}{m}}={\frac {p}{m}}=v={\frac {dx}{dt}}}$   and
${\displaystyle -{\frac {\partial \omega }{\partial x}}=-{\frac {1}{\hbar }}{\frac {\partial V}{\partial x}}={\frac {F}{\hbar }}={\frac {mdv/dt}{\hbar }}={\frac {1}{\hbar }}{\frac {dp}{dt}}={\frac {dk}{dt}}}$ (as expected)

## Conclusion

Had Hamilton been asked in 1834 to find a wave that obeyed Newton's second law of motion for a conservative force, he would have written the Schrödinger equation and said that ħ needs to be a small number. If you had asked "how small", he would have said "very small" and wondered why you wanted such a wave.

## Statement of the theorem in modern form

Suppose we have a linear partial differential equation involving time, ${\displaystyle t}$ , and the spatial coordinates ${\displaystyle x_{i}}$. If the coefficients do not depend on ${\displaystyle t}$ or ${\displaystyle x_{i}}$, then this partial differential equation is a wave equation if it supports infinite plane wave solutions of the form,

${\displaystyle \psi ({\vec {r}},t)=Ae^{i{\vec {k}}\cdot {\vec {r}}-i\omega t}}$

where the dispersion relation is,

${\displaystyle \omega =\omega ({\vec {k}},{\vec {r}},t)}$

This latter form of the dispersion relation has no meaning for waves in a medium that is truly homogeneous in space in time. But now suppose the medium varies so gently with respect to ${\displaystyle {\vec {r}}}$ and ${\displaystyle t}$ that localized plane waves can exist in the form of wave packets. These wave packets must be many wavelength in size, yet small enough that the medium remains nearly uniform over the extent of the wavepacket. The medium must also vary sufficiently slowly with respect to time. This regime of gently variation in properties of the medium that supports these wave packets is associated with geometrical optics and the eikonal equation. In this approximation, wavelength, angular frequency, and wavenumber of a wave packet all obey the following approximate equations of motion:[7]

${\displaystyle {\dot {q_{i}}}={\frac {\partial \omega }{\partial k_{i}}}}$

${\displaystyle {\dot {k_{i}}}=-{\frac {\partial \omega }{\partial q_{i}}}}$

${\displaystyle {\dot {\omega }}=-{\frac {\partial \omega }{\partial t}}}$

#### Applications

• Plasma physics. These equations also describe the propagation of wave energy and have application in plasma physics applications where microwaves might be used to heat the core of the plasma.[8]
• Derivation of Schrödinger's Equation. If a wavepacket were to obey Newton's laws of physics, with the force being the gradient of a scalar potential, ${\displaystyle \ {\vec {F}}=-\nabla V}$, then the dispersion equation must be of the form,

${\displaystyle \hbar \,\cdot \,\omega ({\vec {k}},{\vec {r}},t)={\frac {\hbar ^{2}k^{2}}{2m}}+V({\vec {r}},t)\,,}$

where ${\displaystyle \hbar }$ is any constant small enough so that the wavepackets act as particles. The Schrödinger equation follows directly.

## References

1. Hamilton, ON A GENERAL METHOD IN DYNAMICS (Philosophical Transactions of the Royal Society, part II for 1834, pp. 247–308.) http://www.maths.tcd.ie/pub/HistMath/People/Hamilton/Dynamics/GenMeth.pdf
2. H. Goldstein, Classical Mechanics, Addison Wesley )Cambridge, Mass.) Chapter 9.
3. The fact that (a,b,c) don't depend on time suggests that the medium supporting the wave is static, for example maintained at constant temperature, pressure, etc. The case of time-dependent coefficients is not difficult to analyze, if the system is linear.
4. For example, for some wave equations the integral over all space of |ψ|^2 does not vary with time. In mechanical waves this might correspond to energy conservation. In quantum theory, the amplitude can be "normalized" so that the integral sums to one. When properly designed, the sum of all probabilities equals one, and this insight leads to the postulate that |ψ|^2 is probability density in Schrödinger's equation.
5. Note, for example that eigensolutions exist of the form ψ=u(x)exp(-iωt)
6. Weinberg S, 1962 Eikonal method in magnetoyhdrodynamics, Phys. Rev. 126 1899-909
7. M. A. Lieberman and A. J. Lichtenberg, Principles of Plasma Discharges and Materials Processing, 2nd ed., Wiley Interscience, Hoboken, NJ: Wiley, 2005. Page 506