Normal Modes: Why Coupled Things Beat, and When They Become Chaos

Take two identical pendulums and tie a weak spring between them. Pull one aside, let go, and watch. The first pendulum swings, slows, and falls still — while the second, untouched, somehow winds itself up to full amplitude. Then the second donates everything back, and the first swells again. Energy sloshes between them on and on, a slow heartbeat riding the fast swing, and the motion of either pendulum on its own looks like nothing you could write a formula for. Two coupled masses, and already the picture seems hopeless.

That hopelessness is a lie about the coordinates, not a fact about the system. There exists a change of variables — a rotation of your axes, nothing more — in which the same two pendulums become two completely independent oscillators that never exchange a single joule. In those coordinates each one swings at one fixed frequency forever, clean as a metronome, and the sloshing simply does not happen. The beat you saw was never in the physics. It was an artifact of watching the system in the wrong basis. Find the right basis and the coupling evaporates.

These privileged coordinates are the normal modes, and the whole apparatus of small oscillations — molecular spectra, bridge resonances, the phonons in a crystal — is the claim that any system of coupled oscillators becomes a set of independent ones in the right coordinates. Then we break it. Add one more pendulum on a single arm and the trick collapses, because normal modes only live near equilibrium. Beyond that small neighbourhood lies the double pendulum, and the double pendulum is chaos. One post spans the most-solvable and least-solvable systems in mechanics, separated by a single act of linearisation.

The trick: diagonalise the coupling

Write the physics down and the structure jumps out. Two masses $m$ on springs of stiffness $k$ to the walls, coupled by a spring of stiffness $k_c$. Let $x_1$ and $x_2$ be the displacements from rest. Newton's second law on each mass gives two equations, and the only thing tying them together is the coupling term $k_c(x_1 - x_2)$:

$$\begin{aligned} m\,\ddot{x}_1 &= -k\,x_1 - k_c\,(x_1 - x_2) \\ m\,\ddot{x}_2 &= -k\,x_2 - k_c\,(x_2 - x_1) \end{aligned}$$

Copy LaTeXwith $$ delimitersCopy mathexpression only

The mess lives entirely in the cross-terms. Equation 1 says the acceleration of mass 1 depends on the position of mass 2, and vice versa — the mathematical fingerprint of "they talk to each other". You cannot solve either equation without the other.

The way out is to stop tracking $x_1$ and $x_2$ and track two clever combinations instead. Define the sum coordinate $q_1 = x_1 + x_2$ and the difference coordinate $q_2 = x_1 - x_2$. Add the two lines of equation 1 and the coupling cancels; subtract them and it doubles, but still touches only $q_2$:

$$m\,\ddot{q}_1 = -k\,q_1, \qquad m\,\ddot{q}_2 = -(k + 2k_c)\,q_2$$

Look at what happened. Each new coordinate obeys a single equation mentioning only itself. No cross-terms, no sloshing. $q_1$ and $q_2$ are two independent simple harmonic oscillators that have never heard of each other. The coupling didn't disappear from the world — it disappeared from the coordinates. Each $q$ is a normal mode, and its frequency falls straight out:

$$\omega_1 = \sqrt{\frac{k}{m}} \quad \text{(in phase)}, \qquad \omega_2 = \sqrt{\frac{k + 2k_c}{m}} \quad \text{(out of phase)}$$

Copy LaTeXwith $$ delimitersCopy mathexpression only

The in-phase mode $\omega_1$ is the masses moving together: the coupling spring rides along at constant length, never stretches, stores no energy, and so leaves the frequency untouched — exactly an uncoupled pendulum. The out-of-phase mode $\omega_2$ is the masses moving against each other: now the coupling spring fights the motion, and the extra stiffness $2k_c$ pushes the frequency up. Higher pitch, exactly what you saw switching to Normal mode 2 in the lab.

The coupled system is a single linear operator $\mathbf{K}$ acting on the displacement vector. Diagonalising means finding the basis in which $\mathbf{K}$ has no off-diagonal entries — and in that basis the equations of motion uncouple, by construction. The eigenvalue problem is:

$$\mathbf{K}\,\mathbf{v}_i = m\,\omega_i^2\,\mathbf{v}_i, \qquad \mathbf{K} = \begin{pmatrix} k + k_c & -k_c \\ -k_c & k + k_c \end{pmatrix}$$

Copy LaTeXwith $$ delimitersCopy mathexpression only

Every coupled-oscillator problem you will ever meet — a chain of 3 masses, a 2D lattice of 10⁶ of them, carbon dioxide with its 3 atoms and 4 vibrational modes — is the same sentence: build the stiffness matrix, find its eigenvectors, and the impossible-looking motion resolves into a sum of independent metronomes. Normal modes feel like magic only because the linear algebra has already done the work, handing you the answer in a basis where nothing is coupled.

The coupling never disappears from the world. It disappears from the coordinates — and a normal mode is just the basis in which the mess was hiding.

Beats are two modes interfering

So where did the sloshing come from, if the modes never slosh? It comes from watching a superposition of two modes through the original coordinates $x_1, x_2$. When you displace only the left mass, you have not excited a single mode — you have lit up both at once, in equal measure, because the lopsided start $(x_1, x_2) = (A, 0)$ is the average of mode 1 $(A, A)/2$ and mode 2 $(A, -A)/2$. Each mode then runs at its own frequency, $\omega_1$ and $\omega_2$, and the left mass is their sum:

$$x_1(t) = \frac{A}{2}\cos(\omega_1 t) + \frac{A}{2}\cos(\omega_2 t)$$

Two cosines at nearby frequencies add to a product — the standard trigonometric identity — of a fast carrier at the average frequency and a slow envelope at half the difference:

$$x_1(t) = A\,\cos\!\left(\frac{\omega_2 - \omega_1}{2}\,t\right)\cos\!\left(\frac{\omega_2 + \omega_1}{2}\,t\right)$$

Copy LaTeXwith $$ delimitersCopy mathexpression only

The fast factor is the swing you see; the slow factor is the beat envelope that pumps the amplitude up and down. The beat frequency is $(\omega_2 - \omega_1)/2$ — it depends only on how far apart the two mode frequencies are, which is set by the coupling strength $k_c$. Weak coupling means $\omega_1$ and $\omega_2$ are nearly equal, the envelope is glacially slow, and the energy takes a long time to crawl across. Strong coupling spreads the modes apart and the beat speeds up.

Go back to the lab for a moment and turn the Coupling k_c slider up. The beat gets faster, exactly as equation 4 predicts: a larger $k_c$ pushes $\omega_2$ higher (equation 2), widens the gap $\omega_2 - \omega_1$, and shortens the envelope. Nothing about the modes themselves changed — they are still two clean metronomes — but the rate at which their interference pattern cycles tracks the gap you just widened.

A mode is a closed ellipse in phase space

There is a second way to see that a normal mode is simple, and it lives in phase space — the plane whose axes are position and velocity. A single oscillator at one frequency, plotted as the point $(x, \dot{x})$ over time, traces a closed ellipse: position and velocity trade off, energy conserved, the same loop retraced forever. One frequency, one clean closed curve. That is the geometric signature of a normal mode.

The ellipse is fragile. It stays a perfect closed curve only because the restoring force is linear — proportional to displacement, $F = -kx$. That linearity is itself an approximation: a real pendulum's restoring force is $-mg\sin\theta$, and only for small angles does $\sin\theta \approx \theta$ make it linear. Normal modes are eigenvectors of the stiffness matrix you get by linearising the forces about equilibrium. Stay in the small-amplitude neighbourhood and the ellipses hold, the modes stay independent, the system is exactly solvable. Push the amplitude up and the approximation frays: the ellipse warps, the modes leak into each other, and the clean decomposition fails.

Where it breaks: the double pendulum

Here is the cliff edge. Take the coupled system and make the coupling strong and the amplitudes large by hanging one pendulum directly off the end of another — a double pendulum. The two bobs are now coupled through the rigid arm and through gravity, and there is no weak-spring approximation to hide behind. The equations of motion are still exact Newtonian mechanics, no randomness anywhere. But they are violently nonlinear, full of $\sin(\theta_1 - \theta_2)$ and $\cos(\theta_1 - \theta_2)$ terms that no rotation of axes will ever diagonalise.

For small swings, the double pendulum does have normal modes — two of them, a slow in-phase mode and a fast out-of-phase mode, found by exactly the eigenvalue recipe of equation 3 applied to its linearised stiffness matrix. Release it gently and you would see clean beats, the same heartbeat as the spring-coupled masses. The chaos is not in the equations; it is in the amplitude. Crank the energy up and the system lives far from equilibrium, where the $\sin\theta$ nonlinearity dominates: the eigenvectors stop being constant, modes bleed into one another, and trajectories that start together are stretched apart exponentially.

Normal modes of 2 masses

2

exactly as many modes as degrees of freedom — always

In-phase frequency

√(k/m)

coupling spring never stretches, so it does not appear

Beat frequency

(ω₂−ω₁)/2

set entirely by the coupling strength k_c

The same idea in other rooms

Normal modes are not a pendulum curiosity. They are the universal language for any system of coupled linear oscillators, and the same eigenvalue problem (equation 3) appears in fields that never speak to each other.

In molecular spectroscopy, a molecule of $N$ atoms is $3N$ coupled oscillators — every atom tugged by chemical bonds to its neighbours. Strip out the six rigid-body motions (three translations, three rotations) and diagonalise the bond stiffness matrix, and you get $3N-6$ vibrational normal modes, each at its own frequency. A linear molecule loses one rotation and so keeps $3N-5$ of them: carbon dioxide, three atoms in a line, has four. Those frequencies are exactly what an infrared (IR) spectrometer reads: shine IR light through the sample, and the molecule absorbs at the photon energies that match its mode frequencies. The "fingerprint region" of an IR spectrum is a direct readout of the eigenvalues of equation 3 for that molecule. Carbon dioxide's bending mode near 15 µm — one of those four eigenmodes, and the band that makes it a greenhouse gas — sits almost exactly where the warm Earth radiates hardest.

In electronics, replace masses-on-springs with capacitors-and-inductors. Two LC resonant circuits coupled by a shared inductor obey the same matrix equation as equation 1, with charge playing the role of displacement. They beat, they have in-phase and out-of-phase modes, and the band-pass filter at the front of every radio is two coupled LC tanks tuned so their normal-mode frequencies straddle the channel you want. The engineer tuning a coupled-resonator filter is solving an eigenvalue problem with a screwdriver.

And in solid-state physics, push $N$ to Avogadro's number. A crystal is a 3D lattice of atoms on springs — $10^{23}$ coupled oscillators — and its normal modes are the phonons, the quantised lattice vibrations that carry sound and heat through every solid object you have ever touched. The same diagonalisation, just very large.

Which is what makes the Fermi-Pasta-Ulam-Tsingou problem of 1955 such a beautiful shock. Enrico Fermi, John Pasta, Stanislaw Ulam and Mary Tsingou ran one of the first scientific computations ever, on the MANIAC computer at Los Alamos: a chain of a few dozen masses on springs (the most-cited run used 32) with a weak nonlinear term added. The expectation, from statistical mechanics, was unanimous. Start all the energy in the lowest normal mode, let the nonlinearity couple the modes, and energy should spread until every mode holds its fair share — the chain should thermalise, reaching equipartition the way a hot object reaches uniform temperature. Instead the energy refused. It sloshed among a few low modes, then — astonishingly — came almost all the way back to the first mode, recurring nearly to its starting state. The normal modes were far more robust than anyone had any right to expect.

The coupled pendulums, the IR spectrum, the radio filter, the crystal, and the FPUT chain are one idea wearing five costumes: a system of coupled oscillators is a matrix, its modes are eigenvectors, and the motion is a superposition of independent metronomes — right up until the amplitude grows large enough that the matrix stops being constant, the eigenvectors start to drift, and the linear picture surrenders the system to chaos.

Reading further

• Goldstein, Poole & Safko, Classical Mechanics, 3rd ed., chapter 6 — the canonical derivation of small oscillations as the eigenvalue problem of the stiffness and mass matrices; equation 3 in its full Lagrangian generality.
• A. P. French, Vibrations and Waves — the clearest physical-intuition treatment of coupled oscillators, beats, and normal modes; builds the two-pendulum case slowly and visually before any matrices appear.
• Fermi, Pasta, Ulam & Tsingou, Studies of Nonlinear Problems (Los Alamos report LA-1940, 1955) — the original computation; the chain that refused to thermalise and quietly opened the study of solitons and integrable systems.
• Strogatz, Nonlinear Dynamics and Chaos, chapters 6–9 — the bridge from linear normal modes to the nonlinear phase portraits, limit cycles, and chaos that take over once the small-amplitude approximation fails.