Normal Modes

Strike a bell and it does not ring one note — it rings a chord. Pluck a guitar string anywhere, tap a bridge deck, kick the leg of a table, send an infrared photon through the atmosphere toward a CO2_2 molecule: in every case the thing that vibrates back at you is not a single wiggle but a superposition of pure tones, each with its own frequency and its own frozen pattern of motion. Those pure tones are the normal modes, and they are the deepest organizing idea in all of vibration physics. A normal mode is a collective motion in which every part of the system oscillates at one common frequency, with fixed amplitude ratios — the pattern never changes shape, it only breathes in scale. A system with NN degrees of freedom has exactly NN of these special motions, and here is the miracle: every possible motion, however messy, is a sum of them. The modes are independent — each conserves its own energy, forever, never leaking to a neighbour — so an impossibly coupled NN-body problem collapses into a list of NN separate simple harmonic oscillators. Complexity in, arithmetic out. This article is about the concept: what a mode is, why fixed shape forces one shared frequency, how a struck system divides its energy among its modes, and what the modes of real systems look like — chains of masses, standing waves on a string, the three vibration patterns of CO2_2 that make it a greenhouse gas, the swaying modes of skyscrapers. The concrete two-mass workhorse is solved in Coupled Oscillators; the systematic matrix machinery lives in The Normal-Mode Eigenvalue Problem. Here we learn to think in modes — the skill that later becomes phonons, photons, and the mode expansions of quantum field theory.

Key takeaways
  • A normal mode is a motion in which every coordinate of a coupled system oscillates sinusoidally at one shared frequency with fixed amplitude ratios: the shape is frozen, only its overall scale oscillates.
  • NN degrees of freedom → exactly NN normal modes, each with its own frequency ωn\omega_n and shape a(n)\mathbf{a}^{(n)}.
  • The general motion is a superposition of all NN modes, q(t)=∑na(n)(Ancos⁡ωnt+Bnsin⁡ωnt)\mathbf{q}(t) = \sum_n \mathbf{a}^{(n)}(A_n\cos\omega_n t + B_n\sin\omega_n t), with the 2N2N constants fixed by initial positions and velocities.
  • In a linear system each mode carries its own energy and never trades it with another mode — the modes are the coordinates in which the coupled system becomes NN independent simple harmonic oscillators.
  • Mode shapes of a uniform chain of beads are discrete sine patterns, xj∝sin⁡(nπj/(N+1))x_j \propto \sin(n\pi j/(N+1)), with frequencies ωn=2k/m sin⁡(nπ/2(N+1))\omega_n = 2\sqrt{k/m}\,\sin\big(n\pi/2(N+1)\big); as N→∞N \to \infty these become the standing waves of a string.
  • Real-world galleries: CO2_2's symmetric stretch, antisymmetric stretch, and (doubly degenerate) bend; the sway modes of buildings; the vibrational modes every satellite is tested against.
  • Key idea: a normal mode is a fixed pattern that oscillates at one common frequency; the NN modes of an NN-coordinate system are independent, energy-conserving building blocks whose superposition is every motion the system can perform.
The core insight

A normal mode is a collective motion with one shared frequency and a frozen shape; a linear system with NN degrees of freedom has exactly NN of them, they oscillate independently — each keeping its own energy — and every possible motion is a superposition of modes with amplitudes and phases set by the initial conditions.

Where this lives today

Modal analysis is a daily industrial ritual. Every satellite and launch-vehicle payload is bolted to a shaker table and swept through its normal modes before flight, because a mode of the structure landing on a mode of the rocket is how payloads break; finite-element packages compute galleries of thousands of modes for aircraft wings, EV battery packs, and turbine blades as a matter of course. At the heaviest end, Taipei 101 carries a 660-tonne steel pendulum near its top — a tuned mass damper deliberately built as a second oscillator matched to the tower's lowest sway mode, soaking up its energy exactly by the coupled-resonant-transfer physics of this cluster of articles. The concept's reach now extends to the most exotic objects measured. When two black holes merge, the remnant rings down through its quasinormal modes — damped oscillations of spacetime itself with frequencies set by the final mass and spin; LIGO–Virgo–KAGRA fit these ringdown tones to test general relativity, a program literally called black-hole spectroscopy. And in quantum hardware, the normal modes of nanomechanical resonators and microwave cavities are cooled to their ground states and manipulated one phonon or photon at a time: each mode, being an independent harmonic oscillator, quantizes into its own ladder of energy levels. That last sentence is the door this article opens — a field, classical or quantum, is nothing but an infinite gallery of normal modes, and its quanta are the modes' excitations.

The physical picture

Watch a system move in a normal mode and the striking thing is its discipline. Every mass passes through equilibrium at the same instant; every mass reaches its extreme at the same instant; the ratios between the displacements never change. Freeze the motion at any moment and you always see the same photograph, merely rescaled. That photograph is the mode shape; the rate at which it breathes is the mode frequency. Nothing else about the motion is negotiable — start a system exactly in a mode shape, at rest, and it will oscillate in that shape at that frequency forever. Here is the full portrait gallery for three equal beads coupled by four equal springs between two walls — three degrees of freedom, hence exactly three modes:

Figure
The three normal modes of a three-bead chain — each pattern breathes in scale at its own frequency
three beads, four springs — exactly three modes mode 1 — ω₁ = 0.77√(k/m) — everything together, gentlest mode 2 — ω₂ = 1.41√(k/m) — middle bead never moves rest mode 3 — ω₃ = 1.85√(k/m) — neighbours opposed, stiffest and fastest arrows show the frozen displacement pattern — amplitudes in the exact ratio sin(nπj/4), a discrete sine curve

Read the gallery from top to bottom and a universal trend emerges. The lowest mode is the most collective: all beads swing together, the springs between them barely stretch, and the restoring force per unit displacement is smallest — hence the gentlest frequency. Each successive mode adds one more sign change across the chain, one more "node," stretching the connecting springs harder and driving the frequency up. The highest mode is the most jagged — every neighbour opposed — and the stiffest. Counting nodes orders the spectrum, and this staircase of increasingly wrinkled patterns is exactly what will become the harmonics of a guitar string. Now the second half of the concept, the half that earns modes their keep. Start the system in anything other than a mode shape — strike one bead, say — and the motion looks hopelessly complicated: no bead traces a sinusoid, the pattern writhes, energy wanders from bead to bead. But underneath, the modes are simply running independently, each at its own frequency, and the mess is their sum:

Figure
A messy motion decomposed — what one mass does is exactly the sum of two clean mode sinusoids
what mass 1 actually does — complicated, never repeating exactly mode 1 contribution — pure tone at ω₁ mode 2 contribution — pure tone at ω₂ = + the messy trace is nothing but two clean mode sinusoids added — the modes are the atoms of the motion

The modes are the atoms of motion: indivisible, independent, and universal. Learning a system means learning its mode gallery — after that, every question about its dynamics is a question about a handful of independent simple harmonic oscillators.

What makes a motion a mode

Let us make the defining property precise. A system of NN coupled coordinates moves in a normal mode when every coordinate executes simple harmonic motion at one common frequency, with a fixed set of amplitude ratios and a common phase:

qi(t)=ai cos⁡(ωt−φ),i=1,…,N(1)q_i(t) = a_i\,\cos(\omega t - \varphi), \qquad i = 1, \dots, N \qquad (1)

where: qiq_i — the ii-th coordinate (a bead's displacement, a bond stretch, an angle); aia_i — the fixed amplitude of coordinate ii, whose ratios define the mode shape; ω\omega — the single shared mode frequency; φ\varphi — one common phase for the whole pattern. Two things are worth staring at. Every coordinate shares the same cos⁡(ωt−φ)\cos(\omega t - \varphi): the pattern (a1,…,aN)(a_1, \dots, a_N) is multiplied by one global oscillating factor, which is exactly the statement that the shape only breathes in scale. And the shape is not arbitrary — only NN special shapes are self-consistent, one per mode. Why does a frozen shape force a single shared frequency? Because of how restoring forces work in a linear system, the two properties are actually the same property.

Derivation AOne frozen shape forces one shared frequency

We prove the defining equivalence of normal modes: if a linear system moves so that its shape never changes — all coordinates locked in fixed proportion — then every coordinate must oscillate sinusoidally at one common frequency, and the shape must be one the dynamics itself selects. The strategy is separation of variables in miniature: write the assumed motion as a fixed pattern times a single unknown function of time, feed it to the equations of motion, and watch the equations force that function to be a sinusoid. The step to watch is the "constant of separation" argument — a function of time cannot equal a different function of time unless both are constant — which is the same logic that later splits partial differential equations into space and time parts.

Step 1. Assume the shape is frozen: qi(t)=aif(t)q_i(t) = a_i f(t) for all ii, with fixed numbers aia_i and a single scale function f(t)f(t) shared by every coordinate. Why this is the right formalization: "the pattern only breathes in scale" means precisely that all coordinates track one common time dependence, each scaled by its own fixed amplitude.

Step 2. Feed this into the linearized equations of motion. Near stable equilibrium every coordinate obeys an equation of the form miq¨i=−∑jKij qjm_i \ddot{q}_i = -\sum_j K_{ij}\,q_j, where the constants KijK_{ij} encode the springs and couplings. Substituting the ansatz and dividing by miaif(t)m_i a_i f(t):

f¨(t)f(t)=− ∑jKij ajmi aifor every i\frac{\ddot{f}(t)}{f(t)} = -\,\frac{\sum_j K_{ij}\,a_j}{m_i\,a_i} \quad \text{for every } i

Why divide this way: the left side depends only on time, the right side only on the constants of the problem and the chosen shape — the equality must hold at every instant, for every coordinate.

Step 3. Apply the separation argument. The left side is one function of time; the right side is a constant (a different constant for each ii, potentially). For the equality to hold at all times, the left side must equal a single constant, call it −ω2-\omega^2; and then every right-hand side must equal that same constant. Why the constant is negative for stable systems: near a stable equilibrium the restoring forces push back toward the pattern's zero, making f¨\ddot{f} and ff have opposite signs.

Step 4. Read off both conclusions at once. From the time part, f¨=−ω2f\ddot{f} = -\omega^2 f, so f(t)=cos⁡(ωt−φ)f(t) = \cos(\omega t - \varphi): every coordinate oscillates sinusoidally at the same frequency ω\omega — equation (1). From the shape part, ∑jKijaj=ω2miai\sum_j K_{ij} a_j = \omega^2 m_i a_i for every ii: the shape cannot be arbitrary, but must reproduce itself under the action of the spring network, each coordinate's restoring force per unit inertia coming out in the same proportion ω2\omega^2. Why only special shapes survive: this self-consistency condition is an eigenvalue problem in disguise — for an NN-coordinate system it has exactly NN independent solutions, which is where the count "N degrees of freedom, N modes" comes from. The systematic solution method, det⁡(K−ω2M)=0\det(K - \omega^2 M) = 0, is developed in The Normal-Mode Eigenvalue Problem. *For one coordinate the condition collapses to ka=ω2mak a = \omega^2 m a, i.e. ω=k/m\omega = \sqrt{k/m} — the ordinary SHO; and any shape multiplied by a non-sinusoidal f(t)f(t) fails Step 3, confirming that frozen shape and single shared frequency are inseparable.*

The simplest system with a genuine mode gallery is the two-mass, three-spring chain, worked out completely in Coupled Oscillators by the sum-and-difference trick. Quoting its result: the two shapes are the symmetric pattern (1,1)(1,1), in which the coupling spring never stretches, and the antisymmetric pattern (1,−1)(1,-1), in which it is doubly engaged, with frequencies

ω1=km    for (1,1),ω2=k+2κm    for (1,−1)(2)\omega_1 = \sqrt{\frac{k}{m}} \;\; \text{for } (1,1), \qquad \omega_2 = \sqrt{\frac{k + 2\kappa}{m}} \;\; \text{for } (1,-1) \qquad (2)

where: kk — each wall spring; κ\kappa — the coupling spring; mm — each mass. When all three springs are equal (κ=k\kappa = k) the frequencies are k/m\sqrt{k/m} and 3k/m\sqrt{3k/m}. Note the pattern that generalizes: the collective, node-free shape is slow; the wrinkled shape that fights the coupling is fast — the same staircase as the three-bead gallery above.

Superposition: every motion is a chord

Modes would be a curiosity if systems only ever moved in one of them. Their real power is completeness. Because the equations of motion are linear, any sum of solutions is a solution — and it is a theorem of the eigenvalue machinery that the NN mode shapes span every possible configuration. So the most general motion of the system is a chord: all NN modes sounding at once, each with its own amplitude and phase,

q(t)  =  ∑n=1Na(n) (Ancos⁡ωnt+Bnsin⁡ωnt)(3)\mathbf{q}(t) \;=\; \sum_{n=1}^{N} \mathbf{a}^{(n)}\,\big(A_n\cos\omega_n t + B_n\sin\omega_n t\big) \qquad (3)

where: a(n)\mathbf{a}^{(n)} — the nn-th mode shape; ωn\omega_n — its frequency; An,BnA_n, B_n — 2N2N constants in total, fixed by the NN initial positions and NN initial velocities. This is the classical mechanic's Fourier theorem: initial conditions in, mode amplitudes out, and from then on each term evolves in serene isolation. The construction of the special coordinates that make this decomposition automatic — and that diagonalize the energy — is the business of the Normal Coordinates article; the concept is what matters here: solving a linear vibrating system means finding its modes and dealing the initial conditions out to them. One consequence deserves its own paragraph, because it is the deepest fact in this article: the modes do not talk to each other. Each mode's amplitude, once set by the initial conditions, is set forever; each mode's energy is separately conserved. The complicated wandering of energy you see — from bead to bead, from one end of a struck structure to the other — is entirely an artifact of watching the wrong coordinates. In mode coordinates nothing wanders at all.

Go deeper · Derivation B · 2 minEach mode keeps its own energy
Derivation BEach mode keeps its own energy

We show, concretely on the two-mass system, that the total energy splits exactly into one energy per mode, with no cross terms, and that each mode's share is separately constant in time. The strategy is direct computation: write the kinetic and potential energies in terms of the mode amplitudes (the sum and difference coordinates of the 2-DOF system) and watch the cross terms cancel. The cancellation is not an accident of this example — it happens because the mode shapes are orthogonal in the appropriate sense, a property proved in general in The Normal-Mode Eigenvalue Problem — but seeing it happen by hand is what makes it believable. Watch the factor 14\tfrac14: it comes from the change of variables x1=(s+d)/2x_1 = (s+d)/2, and forgetting it breaks the energy bookkeeping.

Step 1. Set up the two-mass, three-spring system with coordinates x1,x2x_1, x_2 and mode combinations s=x1+x2s = x_1 + x_2, d=x1−x2d = x_1 - x_2, so x1=12(s+d)x_1 = \tfrac12(s+d), x2=12(s−d)x_2 = \tfrac12(s-d). The energy is

E=12mx˙12+12mx˙22+12kx12+12kx22+12κ (x1−x2)2E = \tfrac12 m\dot{x}_1^2 + \tfrac12 m\dot{x}_2^2 + \tfrac12 k x_1^2 + \tfrac12 k x_2^2 + \tfrac12 \kappa\,(x_1 - x_2)^2

Why the last term: only the coupling spring cares about the relative displacement; the wall springs each store energy in their own mass's displacement.

Step 2. Convert the kinetic energy. Squaring and adding, x˙12+x˙22=12(s˙2+d˙2)\dot{x}_1^2 + \dot{x}_2^2 = \tfrac12(\dot{s}^2 + \dot{d}^2) — the cross term s˙d˙\dot{s}\dot{d} appears with coefficient +12+\tfrac12 from x˙12\dot{x}_1^2 and −12-\tfrac12 from x˙22\dot{x}_2^2 and cancels. So T=14ms˙2+14md˙2T = \tfrac14 m\dot{s}^2 + \tfrac14 m\dot{d}^2. Why the cross terms cancel: the shapes (1,1)(1,1) and (1,−1)(1,-1) are orthogonal — the cancellation is the dot product (1,1)⋅(1,−1)=0(1,1)\cdot(1,-1) = 0 doing its work.

Step 3. Convert the potential energy the same way: x12+x22=12(s2+d2)x_1^2 + x_2^2 = \tfrac12(s^2 + d^2) and (x1−x2)2=d2(x_1 - x_2)^2 = d^2, giving

V=14k s2+14(k+2κ) d2V = \tfrac14 k\,s^2 + \tfrac14 (k + 2\kappa)\,d^2

Why the coupling lands only in dd: the symmetric mode never stretches the middle spring, so κ\kappa cannot store energy in ss — the same physics that kept ω1\omega_1 at the bare frequency.

Step 4. Collect the pieces. The energy is a clean sum of two single-oscillator energies with the correct mode frequencies from (2):

E=14m(s˙2+ω12s2)⏟E1  +  14m(d˙2+ω22d2)⏟E2(4)E = \underbrace{\tfrac14 m\big(\dot{s}^2 + \omega_1^2 s^2\big)}_{E_1} \;+\; \underbrace{\tfrac14 m\big(\dot{d}^2 + \omega_2^2 d^2\big)}_{E_2} \qquad (4)

where: E1,E2E_1, E_2 — the mode energies. Since s(t)s(t) and d(t)d(t) each obey their own SHO equation, each bracket is the conserved energy of an independent oscillator: E1E_1 and E2E_2 are separately constant, for all time, whatever the initial conditions. The numeric example below checks E1+E2E_1 + E_2 against the directly computed initial spring energy and gets exact agreement; and in the <a class="wiki-link" href="/wiki/aq-mech-coupled-oscillators" data-track="wiki_article.body.wiki_link">Coupled Oscillators</a> energy slosh, E1E_1 and E2E_2 stay constant even while the two masses trade all their energy — the trade is interference between modes, not flow between them.

Chains of beads. Stretch the two-mass system to NN equal beads between two walls and the mode shapes organize themselves into a stunning pattern: the amplitude of bead jj in mode nn is a sampled sine curve, aj(n)∝sin⁡(nπj/(N+1))a_j^{(n)} \propto \sin(n\pi j/(N+1)) — the discrete sines visible in the three-bead gallery. The frequencies follow a closed formula worth deriving, because it contains the guitar string.

Go deeper · Derivation C · 2 minThe modes of a chain of N beads
Derivation CThe modes of a chain of N beads

We derive the complete mode spectrum of NN equal masses mm joined by equal springs kk between two walls: frequencies ωn=2k/m sin⁡(nπ/2(N+1))\omega_n = 2\sqrt{k/m}\,\sin\big(n\pi/2(N+1)\big) with sine-curve shapes. The strategy is an educated guess with a self-consistency check — try a sampled sine wave as the shape and verify the equations of motion accept it. This works because a sine pattern has a magic property under second differences: shifting a sine left and right and averaging returns the same sine, rescaled. The step to watch is the boundary condition — it is the walls that quantize the allowed wavelengths, exactly as they will for the continuous string.

Step 1. Write the equation of motion of interior bead jj, pulled by its two neighbouring springs: mx¨j=k(xj+1−xj)−k(xj−xj−1)=k(xj+1−2xj+xj−1)m\ddot{x}_j = k(x_{j+1} - x_j) - k(x_j - x_{j-1}) = k(x_{j+1} - 2x_j + x_{j-1}), with the wall conditions expressed as fictitious stationary beads x0=xN+1=0x_0 = x_{N+1} = 0. Why the combination xj+1−2xj+xj−1x_{j+1} - 2x_j + x_{j-1}: each spring pulls toward its far end, so the net force measures how much bead jj sits away from the average of its neighbours — a discrete second derivative.

Step 2. Guess the mode: xj(t)=sin⁡(jθ)cos⁡ωtx_j(t) = \sin(j\theta)\cos\omega t for some angle θ\theta per step. Use the identity sin⁡((j+1)θ)+sin⁡((j−1)θ)=2sin⁡(jθ)cos⁡θ\sin\big((j{+}1)\theta\big) + \sin\big((j{-}1)\theta\big) = 2\sin(j\theta)\cos\theta:

−mω2sin⁡(jθ)=k (2cos⁡θ−2)sin⁡(jθ)-m\omega^2 \sin(j\theta) = k\,\big(2\cos\theta - 2\big)\sin(j\theta)

Why the guess survives: the sine pattern reproduces itself under the neighbour sum — every bead's equation reduces to the same scalar condition, which is precisely the self-consistency a mode requires.

Step 3. Solve the scalar condition using 1−cos⁡θ=2sin⁡2(θ/2)1 - \cos\theta = 2\sin^2(\theta/2):

ω2=2km (1−cos⁡θ)=4km sin⁡2 ⁣θ2⟹ω=2km ∣sin⁡θ2∣\omega^2 = \frac{2k}{m}\,(1 - \cos\theta) = \frac{4k}{m}\,\sin^2\!\frac{\theta}{2} \quad\Longrightarrow\quad \omega = 2\sqrt{\frac{k}{m}}\,\Big|\sin\frac{\theta}{2}\Big|

Step 4. Impose the walls. x0=0x_0 = 0 holds automatically; xN+1=0x_{N+1} = 0 requires sin⁡((N+1)θ)=0\sin\big((N{+}1)\theta\big) = 0, i.e. θ=nπ/(N+1)\theta = n\pi/(N+1). The allowed modes are

ωn=2km sin⁡ ⁣(nπ2(N+1)),n=1,2,…,N(5)\omega_n = 2\sqrt{\frac{k}{m}}\,\sin\!\Big(\frac{n\pi}{2(N+1)}\Big), \qquad n = 1, 2, \dots, N \qquad (5)

where: nn — the mode number, equal to the number of half-wavelengths of the sine pattern across the chain. Exactly NN distinct modes emerge: n=N+1n = N+1 makes every sin⁡(jθ)\sin(j\theta) vanish (no motion), and higher nn merely repeats earlier shapes with a sign — the count comes out right without being imposed. Note the ceiling: no mode of the chain, however wrinkled, exceeds ωmax⁡=2k/m\omega_{\max} = 2\sqrt{k/m}, the frequency of neighbours in full opposition.

Step 5. Take the continuum limit. For n≪Nn \ll N, expand the sine: ωn≈k/m  nπ/(N+1)\omega_n \approx \sqrt{k/m}\;n\pi/(N+1). With bead spacing aa, chain length L=(N+1)aL = (N+1)a, and wave speed v=ak/mv = a\sqrt{k/m}, this is ωn≈nπv/L\omega_n \approx n\pi v/L — in ordinary frequency,

fn=nv2L(6)f_n = \frac{n v}{2L} \qquad (6)

where: vv — the speed of waves on the effective string; LL — its length. These are exactly the harmonics of a string fixed at both ends: the fundamental and its integer multiples. The bead chain's low modes are the string's standing waves, sampled at NN points. Two quick checks confirm the formula: N=1N = 1 gives ω=2k/msin⁡(π/4)=2k/m\omega = 2\sqrt{k/m}\sin(\pi/4) = \sqrt{2k/m}, correct for one mass between two springs; N=2N = 2 gives ω/k/m=1\omega/\sqrt{k/m} = 1 and 3\sqrt{3}, matching the equal-spring case of equation (2).

Figure
The continuum limit — beads on a chain already sit on the standing waves of a string
n=1 n=2 n=3 twelve beads on a string — the mode shapes are already discrete sine curves as N grows the beads fill in the curve: the chain modes become the standing waves of a continuous string

Molecules. A molecule of nn atoms has 3n3n degrees of freedom; removing overall translations and rotations leaves 3n−63n - 6 vibrational modes (3n−53n - 5 if the molecule is linear). CO2_2, linear and triatomic, has exactly four: the symmetric stretch (both oxygens breathe outward, carbon still) at about 1330 cm−11330\ \mathrm{cm^{-1}}, the antisymmetric stretch (oxygens shift one way, carbon the other) at 2349 cm−12349\ \mathrm{cm^{-1}}, and a bend at 667 cm−1667\ \mathrm{cm^{-1}} that counts twice, because the molecule can flex in two perpendicular planes — a degenerate pair. These are not abstractions: infrared light at 667667 and 2349 cm−12349\ \mathrm{cm^{-1}} is absorbed by exactly these modes, which is the molecular mechanism of the greenhouse effect. The symmetric stretch, leaving the molecule's dipole unchanged, is invisible to infrared and shows up in Raman scattering instead — mode shapes decide selection rules.

Figure
The four vibrational modes of carbon dioxide — mode shapes with their measured wavenumbers
CO₂ — linear molecule, 3n − 5 = 4 vibrational modes symmetric stretch — ~1330 cm⁻¹ (Raman only) carbon at rest antisymmetric stretch — 2349 cm⁻¹ (strong infrared) carbon recoils bend — 667 cm⁻¹ (infrared; doubly degenerate: two perpendicular planes) flexes the line the bend and antisymmetric stretch absorb Earth's infrared — the mode gallery of one molecule sets the greenhouse effect

Structures. A skyscraper, a bridge, a turbine blade, a satellite — each has its mode gallery, computed by finite-element software as a giant eigenvalue problem and measured by shaking the real object. A tall building's lowest modes are gentle side-to-side sways with periods of several seconds; higher modes twist and flex it. Engineers keep the mode frequencies away from earthquake spectra and wind-shedding frequencies — precisely because, as Resonance explains, driving any single mode near its frequency amplifies it by its quality factor.

Limiting cases

CaseResultMeaning
N=1N = 1one mode, ω=k/m\omega = \sqrt{k/m}a normal mode is simple harmonic motion; the whole theory reduces to the SHO
N=2N = 2, identical massesshapes (1,1)(1,1) and (1,−1)(1,-1), eq. (2)symmetry alone fixes the shapes — the fully worked case of Coupled Oscillators
N→∞N \to \infty chainωn→nπv/L\omega_n \to n\pi v/L: string harmonicsdiscrete modes become standing waves; mechanics hands off to wave physics
unanchored systemone mode with ω=0\omega = 0zero restoring force along free translation; the "mode" is uniform drift, not oscillation
symmetric system (e.g. CO2_2 bend)two or more modes share one frequencydegeneracy: any combination within the degenerate pair is also a mode
large amplitudeanharmonic terms couple the modesmode independence is a linear privilege; energy begins to leak between modes

Worked examples

1 · Dealing a strike out to the modes (confidence builder)

Two 1.0 kg1.0\ \mathrm{kg} masses sit between walls with springs k=25 N/mk = 25\ \mathrm{N/m} at the walls and a coupling spring κ=12 N/m\kappa = 12\ \mathrm{N/m} between them. Mass 1 is pulled aside 8.0 cm8.0\ \mathrm{cm} and both are released from rest. Find the mode frequencies, decompose the motion into modes, and compute each mode's energy. Does the motion ever repeat?

Step 1. Mode frequencies from (2): ω1=25/1.0=5.0 rad/s\omega_1 = \sqrt{25/1.0} = 5.0\ \mathrm{rad/s} for the shape (1,1)(1,1); ω2=(25+24)/1.0=49=7.0 rad/s\omega_2 = \sqrt{(25 + 24)/1.0} = \sqrt{49} = 7.0\ \mathrm{rad/s} for (1,−1)(1,-1).

Step 2. Decompose the initial state. With s=x1+x2s = x_1 + x_2 and d=x1−x2d = x_1 - x_2: s(0)=d(0)=8.0 cms(0) = d(0) = 8.0\ \mathrm{cm}, all velocities zero. Both modes are excited with equal amplitude — a strike is always a chord.

Step 3. Write the motion from (3): s(t)=8.0cos⁡5ts(t) = 8.0\cos 5t and d(t)=8.0cos⁡7td(t) = 8.0\cos 7t (in cm), so x1(t)=4.0cos⁡5t+4.0cos⁡7tx_1(t) = 4.0\cos 5t + 4.0\cos 7t and x2(t)=4.0cos⁡5t−4.0cos⁡7tx_2(t) = 4.0\cos 5t - 4.0\cos 7t. Neither mass moves sinusoidally, yet the system is just two pure tones playing at once.

Step 4. Mode energies from (4): E1=14m ω12s02=0.25×1.0×25×(0.080)2=0.040 JE_1 = \tfrac14 m\,\omega_1^2 s_0^2 = 0.25 \times 1.0 \times 25 \times (0.080)^2 = \boxed{0.040\ \mathrm{J}} and E2=0.25×1.0×49×(0.080)2=0.078 JE_2 = 0.25 \times 1.0 \times 49 \times (0.080)^2 = \boxed{0.078\ \mathrm{J}}. Cross-check against the directly computed initial energy: the stretched wall spring stores 12(25)(0.080)2=0.080 J\tfrac12(25)(0.080)^2 = 0.080\ \mathrm{J} and the coupling spring 12(12)(0.080)2=0.038 J\tfrac12(12)(0.080)^2 = 0.038\ \mathrm{J}, total 0.118 J=E1+E20.118\ \mathrm{J} = E_1 + E_2. ✓ Those two numbers, 0.0400.040 and 0.078 J0.078\ \mathrm{J}, now remain constant forever — even while the masses trade energy back and forth every π\pi seconds.

Step 5. Periodicity check: the motion repeats when both cos⁡5t\cos 5t and cos⁡7t\cos 7t complete whole cycles, i.e. after t=2π st = 2\pi\ \mathrm{s} (5 and 7 cycles respectively) — periodic, because 7/57/5 is rational. Detune the springs so the ratio is irrational and the exact initial state never recurs.

2 · Weighing a chemical bond with a mode frequency (exam level)

Model CO2_2 as three masses on a line — mO=16.0 um_{\mathrm{O}} = 16.0\,u, mC=12.0 um_{\mathrm{C}} = 12.0\,u — joined by two identical bond "springs" kk. The antisymmetric stretch is measured at ν~3=2349 cm−1\tilde\nu_3 = 2349\ \mathrm{cm^{-1}}. Find the bond stiffness kk, predict the symmetric-stretch wavenumber, and compare with the measured ≈1333 cm−1\approx 1333\ \mathrm{cm^{-1}}.

Step 1. Convert the wavenumber to an angular frequency: ω3=2πc ν~3=2π(2.998×1010 cm/s)(2349 cm−1)=4.42×1014 rad/s\omega_3 = 2\pi c\,\tilde\nu_3 = 2\pi (2.998\times10^{10}\ \mathrm{cm/s})(2349\ \mathrm{cm^{-1}}) = 4.42\times10^{14}\ \mathrm{rad/s}.

Step 2. The mode frequencies of the linear three-mass model (standard results of the eigenvalue machinery) are ω1=k/mO\omega_1 = \sqrt{k/m_{\mathrm{O}}} for the symmetric stretch — carbon sits still, so each oxygen oscillates on its own bond spring — and ω3=k (1/mO+2/mC)\omega_3 = \sqrt{k\,(1/m_{\mathrm{O}} + 2/m_{\mathrm{C}})} for the antisymmetric stretch, where the light carbon's recoil against both springs raises the frequency.

Step 3. Solve for the bond stiffness. With mO=2.657×10−26 kgm_{\mathrm{O}} = 2.657\times10^{-26}\ \mathrm{kg} and mC=1.994×10−26 kgm_{\mathrm{C}} = 1.994\times10^{-26}\ \mathrm{kg}: 1/mO+2/mC=3.76×1025+1.00×1026=1.38×1026 kg−11/m_{\mathrm{O}} + 2/m_{\mathrm{C}} = 3.76\times10^{25} + 1.00\times10^{26} = 1.38\times10^{26}\ \mathrm{kg^{-1}}, so k=ω32/(1.38×1026)=(1.96×1029)/(1.38×1026)=1.4×103 N/mk = \omega_3^2 / (1.38\times10^{26}) = (1.96\times10^{29})/(1.38\times10^{26}) = \boxed{1.4\times10^{3}\ \mathrm{N/m}}. A chemical double bond is in fact softer than a car's suspension spring by more than an order of magnitude; what makes it ring at 1014 rad/s10^{14}\ \mathrm{rad/s} is not stiffness but the minute mass it has to move — measured with light.

Step 4. Predict the symmetric stretch: ω3/ω1=1+2mO/mC=1+32/12=1.91\omega_3/\omega_1 = \sqrt{1 + 2m_{\mathrm{O}}/m_{\mathrm{C}}} = \sqrt{1 + 32/12} = 1.91, so ν~1=2349/1.91=1227 cm−1\tilde\nu_1 = 2349/1.91 = \boxed{1227\ \mathrm{cm^{-1}}}.

Step 5. Compare honestly: the measured value is ≈1333 cm−1\approx 1333\ \mathrm{cm^{-1}}, so the two-spring model comes in about 8% low. The discrepancy is physics, not arithmetic — real CO2_2 has a bond–bond interaction (stretching one bond stiffens the other) that a nearest-neighbour spring model omits, and the bending mode at 667 cm−1667\ \mathrm{cm^{-1}} needs a separate bending stiffness entirely. The model captures the mode structure exactly and the numbers to ten percent — and its failure pinpoints the next term in the molecular force field.

Origins · 1 minDaniel Bernoulli and how this was discovered
DB
Daniel Bernoulli
1700–1782 · Basel and St Petersburg

Daniel Bernoulli was the first person to see that a complex vibration is a sum of simple ones. In 1733 he analyzed the hanging chain — a rope swinging from a fixed point — and found its family of oscillation shapes, each with its own frequency: the first mode gallery ever computed, complete with the function later named after Bessel. Twenty years later, in the great vibrating-string controversy with d'Alembert and Euler, Bernoulli staked out the claim that every motion of a string is a superposition of its sinusoidal modes — pure physical insight, since the mathematics to prove it (Fourier's) lay half a century in the future; his rivals' objections were sharpened into exactly the questions Fourier analysis eventually answered. Between these bookends he wrote Hydrodynamica (1738), founding fluid dynamics and the kinetic picture of gases. Modal superposition — the working language of everything from molecular spectroscopy to structural engineering — began as Bernoulli's minority opinion.

History
1733
— Daniel Bernoulli publishes the modes of the hanging chain: the first computed gallery of normal modes, with amplitude patterns given by what would later be called Bessel functions.
1747
— d'Alembert derives the wave equation for the vibrating string, framing the continuum problem the mode picture must explain.
1753
— Bernoulli asserts that every motion of a string is a superposition of its simple modes, igniting a decades-long dispute with Euler and d'Alembert about what functions a sine series can represent.
1822
— Fourier's Théorie analytique de la chaleur supplies the missing mathematics, turning Bernoulli's physical claim into the theorem underlying all mode decompositions.
1912
— Debye counts the normal modes of a crystal to predict its heat capacity — the step that, with quantization, turns classical modes into phonons.

Going further

References. Morin, Introduction to Classical Mechanics, Ch. 4 (Oscillations) — coupled oscillators and normal modes · Goldstein, Poole & Safko, Classical Mechanics, 3rd ed., Ch. 6 (Oscillations) — eigenvalue formulation · Landau & Lifshitz, Mechanics, 3rd ed., §23–24 (Oscillations of systems with several degrees of freedom) · Taylor, Classical Mechanics, Ch. 11 (Coupled Oscillators and Normal Modes) The two-mass mode frequencies ω=k/m\omega=\sqrt{k/m} and 3k/m\sqrt{3k/m} and the mass-chain dispersion ωn=2k/msin⁡[nπ/2(N+1)]\omega_n=2\sqrt{k/m}\sin[n\pi/2(N+1)] are standard results; molecular mode counts 3n−63n-6/3n−53n-5 are textbook. Worked-example numbers (mode energies 0.0400.040 and 0.0780.078 J summing to the direct 0.1180.118 J; CO2_2 bond constant k≈1.4×103k \approx 1.4\times10^3 N/m from ν~3=2349 cm−1\tilde\nu_3 = 2349\ \mathrm{cm^{-1}}, predicted ν~1=1227\tilde\nu_1 = 1227 vs measured ≈1333 cm−1\approx 1333\ \mathrm{cm^{-1}}) are arithmetic consequences of the boxed formulas; CO2_2 wavenumbers are standard spectroscopic values.

Builds on: Simple Harmonic Motion · Small Oscillations About Equilibrium · Complex Exponential Representation of Oscillations · Amplitude, Phase, and Angular Frequency Leads to: Driven Oscillations · Resonance · Damped Harmonic Oscillator · Energy in Simple Harmonic Motion

Is a normal mode a motion, a shape, or a coordinate? The article seems to use all three.

All three faces belong to one object, used in different sentences. The motion is the full time history, equation (1): everything oscillating together at ωn\omega_n. The shape is the frozen amplitude pattern a(n)\mathbf{a}^{(n)} — what you would photograph. And associated with each mode there is also a coordinate: the single number ξn(t)\xi_n(t) measuring how much of the shape is currently present, which evolves like an independent SHO. Physicists slide between the three because the map is one-to-one: give the shape, and the frequency and coordinate come with it. The coordinate viewpoint — where the system's state is the list (ξ1,…,ξN)(\xi_1, \dots, \xi_N) and the dynamics is NN separate oscillators — is developed properly in the Normal Coordinates article; this article's superposition formula (3) is exactly that list written back in ordinary coordinates.

Why exactly N modes for N degrees of freedom — not more, not fewer?

The self-consistency condition of the first derivation — the shape must reproduce itself under the spring network, ∑jKijaj=ω2miai\sum_j K_{ij}a_j = \omega^2 m_i a_i — is an N×NN\times N eigenvalue problem, and such a problem has exactly NN independent solutions (counting degenerate ones separately). Intuitively: fewer than NN modes could not span all possible initial configurations, so some strikes would have no decomposition; more than NN would be redundant, since NN independent shapes already fill the configuration space. The bead chain shows the count working mechanically: sine patterns with n=1,…,Nn = 1, \dots, N are all distinct, n=N+1n = N+1 moves nothing, and higher nn repeats old shapes. The general proof — via det⁡(K−ω2M)=0\det(K - \omega^2 M) = 0 being degree NN in ω2\omega^2 — lives in The Normal-Mode Eigenvalue Problem.

If every mode is perfectly periodic, why does the article say the general motion never repeats?

A sum of periodic functions is periodic only if their periods fit into a common one — that is, if all frequency ratios are rational. Two modes at 55 and 7 rad/s7\ \mathrm{rad/s} recur together every 2π2\pi seconds (5 cycles meets 7 cycles), which is why the worked example repeats. But generic mode frequencies are incommensurate — ω2/ω1\omega_2/\omega_1 irrational — and then the modes never simultaneously return to their starting phases: the motion comes arbitrarily close to its initial state infinitely often but never exactly repeats. This is not chaos — the motion is perfectly regular and predictable, just quasi-periodic. Real systems are effectively always in this case, which is why a struck bell's waveform, examined closely, never quite recurs even while each partial rings at its own steady frequency.

Where does mode independence actually break down in real systems?

At the same place the whole theory begins: the linear approximation. Modes decouple because the equations of motion are linear, which is the small-oscillation expansion of the true potential; push the amplitudes up and the neglected cubic and quartic terms act as weak couplings between modes, letting energy leak from one to another. The classic demonstration is the Fermi–Pasta–Ulam–Tsingou computation of 1955: energy fed into the lowest mode of a slightly nonlinear chain wandered among a few modes and then — famously, unexpectedly — almost entirely returned, instead of spreading evenly. That surprise seeded both chaos theory and soliton physics. In practice: molecular vibrations at high excitation exchange energy between modes (intramolecular energy redistribution), and structural engineers treat mode coupling as the signature that deformations have left the linear regime. Damping also couples what pure elasticity keeps separate.

What exactly is a zero-frequency mode, and is it really a mode?

Formally it solves the mode condition with ω=0\omega = 0; physically it is a direction of motion with no restoring force. Take the bead chain and remove the walls: shifting every bead equally stretches no spring, so the pattern (1,1,…,1)(1,1,\dots,1) feels zero force. Its "oscillation" equation is ξ¨=0\ddot{\xi} = 0, whose solution is uniform drift, ξ(t)=ξ0+v0t\xi(t) = \xi_0 + v_0 t, not a sinusoid — the free translation of the whole system. Every unanchored object carries such modes (three translations, three rotations for a molecule — exactly the 66 subtracted in the 3n−63n-6 count). They matter practically: in any numerical mode computation, zero (or nearly zero) frequencies flag either genuine rigid-body freedom or a modeling error, and the degeneracy and zero-mode bookkeeping is handled systematically in The Normal-Mode Eigenvalue Problem.

Are mode shapes always sine curves, like in the bead chain?

No — sines are the special privilege of uniform systems. The chain's identical masses and springs make the discrete translation symmetry pick sinusoidal patterns, just as a uniform string gets sin⁡(nπx/L)\sin(n\pi x/L). Change anything — a heavier bead in the middle, a stiffer spring at one end, a tapered beam — and the mode shapes deform into system-specific patterns that only the eigenvalue machinery can find; Bernoulli's hanging chain, with gravity playing the role of position-dependent tension, has Bessel-function shapes. What survives in general is the structure: NN shapes, orthogonal in the mass-weighted sense, ordered by node count, each with its frequency. The sine-shaped case is still the one to memorize, because it is the bridge to Fourier series: for uniform systems, mode decomposition and Fourier decomposition are the same operation.

How do classical normal modes turn into phonons and field quanta?

Through the one fact this article kept emphasizing: each mode is an independent harmonic oscillator. Quantum mechanics knows exactly what to do with a harmonic oscillator of frequency ωn\omega_n — its energy comes in steps of ℏωn\hbar\omega_n. Quantizing a vibrating crystal therefore means quantizing each of its 3N3N modes separately: a state of the crystal is a list of occupation numbers, "how many quanta in each mode," and one quantum of a lattice mode is a phonon. The identical construction applied to the electromagnetic field's modes gives photons. Nothing about the mode concept changes — shapes, frequencies, and independence carry straight over; only the amplitude of each mode becomes a quantum ladder instead of a continuous dial. This is why normal modes are called the classical ancestor of field quanta, and why time invested here pays off all the way into quantum field theory.

Before this

Where this leads