3.1 A discrete chain of masses

Imagine a horizontal row of NN identical small masses, each of mass mm, connected to neighbours by identical springs of stiffness κ\kappa. Equilibrium spacing aa. Let yn(t)y_n(t) denote the transverse displacement of the nn-th mass from equilibrium.

The force on mass nn from the spring on its left is κ(yn1yn)\kappa\, (y_{n-1} - y_n) (pulling toward yn1y_{n-1}); from the spring on its right, κ(yn+1yn)\kappa\, (y_{n+1} - y_n). Newton’s second law gives

my¨n  =  κ(yn+12yn+yn1).m\, \ddot y_n \;=\; \kappa\, (y_{n+1} - 2 y_n + y_{n-1}).

This is a coupled system of NN second-order ODEs. Without the coupling (the term in parentheses), each mass would just oscillate independently at ω0=κ/m\omega_0 = \sqrt{\kappa/m}. With the coupling, the masses talk to each other: a kick on one propagates to its neighbours.

A travelling pulse

Watch what happens if you set the leftmost mass moving and leave the rest at rest. The first mass pulls on the second; the second responds, and starts pulling on the third; the third responds in turn. After some time, the disturbance has reached the far end of the chain. Energy that was originally on the left has been transported, mass by mass, to the right.

The speed of that transport depends on κ\kappa, mm, and aa. By dimensional analysis,

c    aκ/m.c \;\sim\; a \sqrt{\kappa / m}.

(Why? The units of κ/m\kappa/m are 1/time21/\text{time}^2; aa has units of length; so aκ/ma\sqrt{\kappa/m} has units of length/time. There is no other dimensionally consistent combination.) We will derive the exact relation — same form, prefactor 11 — in the next lesson.

Normal modes

Because the equations are linear and the chain is uniform, the system has NN normal modes, each oscillating sinusoidally in time at its own frequency. For a chain with fixed ends, the kk-th mode has spatial shape

yn(k)    sin ⁣(knπN+1),y_n^{(k)} \;\propto\; \sin\!\left( \frac{k\, n\, \pi}{N+1} \right),

with k=1,2,,Nk = 1, 2, \ldots, N. The mode frequencies are

ωk  =  2κ/msin ⁣(kπ2(N+1)).\omega_k \;=\; 2\sqrt{\kappa/m}\, \sin\!\left( \frac{k\pi}{2(N+1)} \right).

For large NN and small kk, sin(kπ/2(N+1))kπ/(2(N+1))\sin(k\pi / 2(N+1)) \approx k\pi/(2(N+1)), so

ωk    kπN+1κ/m  =  ckwavenumber,\omega_k \;\approx\; \frac{k\pi}{N+1}\, \sqrt{\kappa/m} \;=\; c\, k_\text{wavenumber},

with kwavenumber=kπ/Lk_\text{wavenumber} = k\pi / L for total length L=(N+1)aL = (N+1) a. The low modes of the discrete chain look exactly like the low modes of a continuous string. We will see in the next lesson that this is no accident.

The continuum limit

Now imagine taking NN \to \infty and a0a \to 0 while keeping the total length L=NaL = N a fixed and the linear mass density μ=m/a\mu = m/a fixed. The chain becomes a continuous string. The discrete second difference yn+12yn+yn1y_{n+1} - 2 y_n + y_{n-1} approaches a22y/x2a^2 \partial^2 y / \partial x^2. The discrete equation of motion becomes the wave equation — the topic of the next lesson.