# Dispersion Elation

This blog has been getting a bit too pop-science for my tastes recently. Card games? Word clouds? Urgh. Let’s do some proper physics. I hope you’re paying attention at the back.

I’d like to talk about dispersion relations. As an illustrative example take the humble simple harmonic oscillator, where the position of something ( $a$) is governed by the ODE $\frac{d^2a}{dt^2} + \omega_a^2a = 0$

where $\omega_a$ is the natural frequency of the oscillator. If left by itself, this ODE states that $a$ will simply oscillate harmonically (hence the name) at frequency $\omega_a$. If we want to know the behaviour of $a$ at different frequencies, we could apply them all at once and observe the resulting motion, which means solving $\frac{d^2a}{dt^2} + \omega_a^2a = \delta(t)$

The $\delta$-function represents the action of an infinitesimally short impulse to $a$ at $t = 0$. The simplest approach is in Fourier space, so after a Fourier transform the ODE becomes $(-\omega^2 + \omega_a^2)\tilde{a}(\omega) = 1$

where $d/dt \rightarrow i\omega$ and the Fourier transform of a delta function is constant in frequency space. Solving for the amplitude of $\tilde{a}$, $\tilde{a}(\omega) = \frac{1}{\omega_a^2 - \omega^2}$

and inverse Fourier-transforming back to obtain the temporal motion, $a(t) = \int_{\infty}^{\infty}\frac{e^{-i\omega t}}{\omega_a^2 - \omega^2}\,d\omega$

The subsequent motion is therefore composed of a sum of oscillations at different frequencies $\omega$, but the amplitude of each is weighted by how close $\omega$ is to the natural frequency. The farther from $\omega_a$, the lower the amplitude of the oscillation. This is a more precise statement of the fact that the oscillator ‘wants’ to move at its natural frequency.

This statement can be generalised. Consider a system governed by some linear differential operator $\hat{L}(\partial_t, \partial_t^2,...)a(t) = 0$

We again take the Fourier transform of both sides, rearrange, and transform back: $a(t) = \int_{\infty}^{\infty}\frac{e^{-i\omega t}}{D(\omega)}\,d\omega$

Here $D(\omega)$ is the dispersion relation corresponding to the operator $\hat{L}$, and in general the motion of $a$ is dominated by frequencies which cause $D$ to be as small as possible. In the SHM case $D = 0$ implies $\omega = \omega_a$.

Let’s add a little complexity to the simple case. Consider two oscillators $a$ and $b$ which are coupled together, i.e. the motion of one drives the motion of the other. The strength of this coupling is given by the constants $c_1$ and $c_2$: $\frac{d^2a}{dt^2} + \omega_a^2a = c_1b$ $\frac{d^2b}{dt^2} + \omega_b^2b = c_2a$

Now with two degrees of freedom the motion is represented by a vector, and the linear operator by a matrix: $\left(\begin{array}{cc} \frac{d^2}{dt^2} + \omega_a^2& -c_1 \\ -c_2 & \frac{d^2}{dt^2} + \omega_b^2 \end{array}\right)\left(\begin{array}{c}a \\ b \end{array}\right) = 0$

For the motion not to be trivial, the determinant of the matrix should be zero, which leads to the coupled dispersion relation $D(\omega, \omega') = (\omega_a^2 - \omega^2)(\omega_b^2 - \omega'^2) - c_1c_2$

where the frequencies of the two oscillators are $\omega$ and $\omega'$. If the couplings $c_1, c_2$ go to zero, then $D(\omega, \omega')$ factorises to the two independent dispersion relations for the two oscillators, $D(\omega, \omega') = D_a(\omega)D_b(\omega')$.

Let’s see what $D$ looks like in this case:

Everything looks normal. $D$ is minimised when both oscillators move at their own resonant frequencies. Let’s add some coupling:

This is interesting. At the resonant frequencies of the oscillators, it’s no longer possible for $D$ to be small. The problem is that we’re missing a chunk of the space that $D$ lives in. Let’s allow one of the frequencies to be imaginary, e.g. set $\omega \rightarrow \omega + i\gamma$.

The dispersion relation is now 3-dimensional, so we can’t visualise it in the same way. Here I’ve picked the isosurface of $D$ evaluated at $D = 10^{-4}$, and rendered it below using Blender:

Now we can see what has happened – our coupled system was forced out of the real plane near resonance, and the dispersion relation has had to develop branches above and below the real plane. Farther from resonance the branches fall back to the real plane, where they merge back into the real dispersion relation.

Let’s see what happens as we increase and decrease the coupling strengths: Higher coupling strengths cause to ‘open up’. My laptop wasn’t too happy to be forced to render this animation.

Stronger coupling pushes $D$ farther out of the real plane, forcing the frequency of the oscillator to have larger imaginary component.

This is all very pretty, but let’s not forget the purpose of this post: hardcore physics. Back to the maths, let’s suppose that $\gamma$ isn’t too large, and so the two oscillators are near their resonant frequencies.

Explicitly let $a$ have frequency $\omega = \omega_a + i\gamma$ and let $b$ have frequency $\omega' = \omega_b + i\gamma$. Substituting into the dispersion relation, $\gamma^2(\gamma - 2i\omega_a)(\gamma - 2i\omega_b) = c_1c_2$

Assuming $\gamma \ll \omega_a, \omega_b$, then after multiplying out and discarding the higher-order terms, $\gamma \approx \sqrt{\frac{-c_1c_2}{4\omega_a\omega_b}}$

This is a very useful expression. If the couplings are the same sign, then $\gamma$ is imaginary and the effect of the coupling is to alter the frequencies of the oscillators. If the couplings are opposite signs then $\gamma$ is real, and so represents an exponentially growing solution.

To check this, I numerically evaluated the ODEs above and calculated the expected growth curve, plotted below in yellow.

You can see that the agreement is excellent for this low growth rate, and at later times the oscillating solutions almost exactly follow the curve.

A second interesting thing is that the oscillators adjust their frequencies such that they are approximately $\pi/2$ out of phase. To see why this might occur, let’s examine growth of the total oscillator energy $E$, i.e. $\frac{d}{dt}\left(\dot{a}^2 + \omega_a^2a^2 + \dot{b}^2 + \omega_b^2b^2\right)$

From the ODEs, we can do a bit of algebra to show that $\frac{dE}{dt} = \frac{d}{dt}\left(c_1\dot{a}b + c_2a\dot{b}\right)$

If we assume that the oscillators are at their own natural frequencies, and that they are $\phi$ out of phase, then after a little more algebra we can take the time-average and show $\langle \frac{dE}{dt} \rangle = \frac{1}{2}(c_2\omega_b - c_1\omega_a)\sin\phi$

Therefore, for maximum growth in energy $\sin\phi = 1$ and $\phi = \pi/2$. Also, for $c_1, c_2$ both the same sign and the natural frequencies equal, there is no net growth in oscillator energy.

Finally there are some other secondary effects at play here. If the natural frequencies are different and the coupling not very strong, the growth doesn’t occur forever:

You can see that what’s happened is the oscillators have gotten out of a growing phase with respect to one another. They’d like to be at $\phi = \pi/2$, but the fact that one would like to oscillate a bit more quickly means that they can’t sustain that relative phase difference forever.

However, bump up the coupling a bit and the oscillators will happily stick at $\phi = \pi/2$ even when their frequencies are different. The growth is so rapid that there’s no time for the oscillator to revert to it’s natural frequency, and so the oscillators are said to be ‘phase-locked’. Very roughly, this occurs when $\gamma > |\omega_a - \omega_b|$.

You might well wonder what the point of all of this is. My personal interest arises from plasma physics, where more or less everything reduces to an oscillation of some sort. One oscillator is the electron fluid, which oscillates at the plasma frequency $\omega_p$. A second common oscillator is some form of electromagnetic wave at frequency $\omega_{EM}$. Typically these frequencies are very different, but can be coupled through the action of a third high-amplitude electromagnetic wave $\omega_0$, commonly provided by a laser pulse. In this case the coupling strengths $c_1 , c_2$ are proportional to the laser intensity and the density of the plasma. When the laser intensity is very large, all sorts of interesting nonlinear oscillator phenomena come into play. One example of which is Stimulated Raman Scattering (SRS) where a plasma resonantly scatters a laser pulse if $\omega_{EM} = \omega_0 - \omega_p$.

The investigation of SRS in certain kinds of laser-plasma interactions formed a core of my PhD thesis so I thought it best to write a little about it, even tangentially, while I can claim a reason to.

## 3 thoughts on “Dispersion Elation”

1. bear24rw says:

What colormap and graphing package was used to make those heatmaps?

Like

2. Floyd Kelly says:

Just read all of what you published – holy crap!

Like