Linear Algebra for Normal Modes (IPhO Math Tools)

IPhO math tools lesson on normal modes: turning coupled oscillations into an eigenvalue problem and using symmetry and 2 by 2 shortcuts under exam time.

  • International Physics Olympiad preparation
On this page

Normal modes are where linear algebra pays back instantly: a messy set of coupled equations becomes “find eigenvalues and eigenvectors”. In olympiad scripts, you often get full method credit by writing the matrix form and the determinant condition cleanly.

1. Definitions (Must Know)

A. Vectors and matrices (notation)

  • A column vector of coordinates:
    x = (x₁; x₂; ⋮; x_N)
  • A matrix mapping vectors to vectors: Ax.
  • Transpose: A^T.
  • Identity matrix: I.

B. Eigenvalues and eigenvectors

An eigenvector v ≠ 0 of A satisfies

Av = λ v.

Nontrivial solutions exist only if

det(A-λI) = 0.

C. Symmetric matrices

If A = A^T (real symmetric), then:

  • eigenvalues are real,
  • eigenvectors can be chosen orthogonal.

This is why normal-mode problems (built from energies) are usually clean.

D. Normal modes (physics meaning)

A normal mode is a motion where all coordinates oscillate sinusoidally at one frequency with fixed amplitude ratios:

xᵢ(t) = aᵢ cos(ω t + φ).

The vector a is the mode shape.

2. Key Ideas (What Earns Marks)

  • Write coupled oscillations as a matrix equation:
    Mx double dot + Kx = 0.
  • Try a harmonic form x = ae^(iω t) to get
    (K-ω²M)a = 0.
  • Nontrivial modes require:
    det(K-ω²M) = 0.
  • Use symmetry to guess eigenvectors first, then find eigenvalues quickly.
  • Check limiting cases (coupling spring constant to zero, one mass very large, and so on).

3. Detailed Explanations

A. Where M and K come from

For small oscillations in generalized coordinates x, the energy is typically quadratic:

T = 1/2x dot ^TMx dot, V = (1/2)x^TKx.

Lagrange’s equations give:

Mx double dot + Kx = 0.

Here:

  • M is the (symmetric) mass matrix.
  • K is the (symmetric) stiffness matrix.

B. Reducing to an eigenvalue problem

Assume a mode solution x(t) = ae^(iω t).

Then x double dot = -ω²x, so:

(K-ω²M)a = 0.

This is a generalized eigenvalue problem:

  • eigenvalues are ω²,
  • eigenvectors are mode shapes a.

C. The 2 by 2 shortcut (trace and determinant)

For a 2 by 2 matrix A, eigenvalues satisfy:

λ²-(tr A)λ + det(A) = 0.

So:

λ = (tr A± square root of ((tr A)²-4 det(A)))/2.

If M = mI, then ω² are just eigenvalues of K/m.

D. Orthogonality (why modes decouple)

For symmetric M and K, different modes satisfy an M-weighted orthogonality:

aᵢ^TMaⱼ = 0 (i ≠ j).

This is why you can expand a general motion as a sum of modes and solve each one independently.

E. Limiting checks that catch algebra mistakes

  • If a coupling constant goes to zero, modes should reduce to independent oscillators.
  • If all masses scale by a factor, frequencies scale as ω ∝ 1/square root of m.
  • If all spring constants scale by a factor, frequencies scale as ω ∝ square root of k.

4. Common Mistakes

  • Writing the coupling forces with the wrong sign (a coupling spring exerts opposite forces on the two masses).
  • Forgetting a “zero mode” when the system has translational symmetry (free system has a rigid-translation mode with ω = 0).
  • Solving for ω and forgetting the square root (solving for ω² is often the clean step).
  • Not checking the k_c → 0 or m₂ → ∞ limits.
  • Mixing mode shapes with initial conditions (mode shapes are geometry; constants come later).

5. Exam Tips

  • Draw the mode shapes first if symmetry is strong. For two identical masses, “in-phase” and “out-of-phase” are almost always eigenvectors.
  • Write the determinant condition explicitly: det(K-ω²M) = 0.
  • If you only need frequencies, you may not need normalized eigenvectors.
  • Use energy language if it is quicker: V quadratic directly gives K.

6. Worked Examples

A. Two identical masses with a coupling spring

Two masses m are attached to walls by springs k and coupled by a spring k_c:

  • left wall to mass 1: k
  • between masses: k_c
  • mass 2 to right wall: k

Find the normal-mode frequencies.

Click here to show/hide solution

Let displacements be x₁,x₂ from equilibrium.

Forces:

mx double dot ₁ = -kx₁-k_c(x₁-x₂), mx double dot ₂ = -kx₂-k_c(x₂-x₁).

Write x = (x₁; x₂) so

mx double dot + Kx = 0, K = (k + k_c, -k_c; -k_c, k + k_c) .

Since M = mI, we solve Ka = ω² m a.

Use symmetry guesses:

  1. In-phase mode a ∝ (1; 1):
K (1; 1) = (k; k) = k (1; 1) .

So ω₁² = k/m.

  1. Out-of-phase mode a ∝ (1; -1):
K (1; -1) = (k + 2k_c) (1; -1) .

So ω₂² = (k + 2k_c)/m.

Therefore:

ω₁ = square root of (k/m), ω₂ = square root of ((k + 2k_c)/m) .

B. Two masses connected by one spring (free system)

Two masses m₁ and m₂ are connected by a spring of constant k in space (no walls). Find the mode frequencies.

Click here to show/hide solution

Let displacements be x₁,x₂ along the spring axis.

The spring force depends on extension (x₁-x₂):

m₁x double dot ₁ = -k(x₁-x₂), m₂x double dot ₂ = -k(x₂-x₁).

There is a translation symmetry: if x₁ = x₂ always, the spring is never stretched and the system can move as a rigid body. That gives a zero-frequency mode:

ω₁ = 0.

For the nontrivial mode, use the relative coordinate u = x₁-x₂. Subtract the equations:

u double dot = x double dot ₁-x double dot ₂ = -k(1/m₁ + 1/m₂)u.

So the oscillation frequency is:

ω₂ = square root of (k(1/m₁ + 1/m₂)) .

7. Mind Stretchers

A. Three-mass chain between fixed walls

Three identical masses m are connected in a line between two fixed walls by four identical springs k:

wall, spring, mass, spring, mass, spring, mass, spring, wall.

Find the three normal-mode frequencies and a convenient set of mode shapes.

Click here to show/hide hint

For an N-mass fixed-end chain, trial eigenvectors are discrete sines:

aⱼ⁽ⁿ⁾ = sin(jnπ/(N + 1)), j = 1,2,…,N.

Then ωₙ² = (2k/m)(1- cos(nπ/(N + 1))).

Click here to show/hide answer

Here N = 3, so

ωₙ² = (2k/m)(1- cos(nπ/4)), n = 1,2,3.

Compute:

  • n = 1: cos(π/4) = square root of 2/2
    ω₁² = (2k/m)(1-(square root of 2)/2)
  • n = 2: cos(π/2) = 0
    ω₂² = 2k/m
  • n = 3: cos(3π/4) = - square root of 2/2
    ω₃² = (2k/m)(1 + (square root of 2)/2)

One convenient set of (unnormalized) mode shapes from aⱼ⁽ⁿ⁾ = sin(jnπ/4):

  • n = 1: (sin(π/4), sin(π/2), sin(3π/4)) ∝ (square root of 2/2,1, square root of 2/2)
  • n = 2: (sin(π/2), sin(π), sin(3π/2)) ∝ (1,0,-1)
  • n = 3: (sin(3π/4), sin(3π/2), sin(9π/4)) ∝ (square root of 2/2,-1, square root of 2/2)

8. Practice

Practice

To get fast at modes:

  • write Mx double dot + Kx = 0 from energy or forces
  • try symmetry mode shapes before you touch a determinant
  • always check the coupling-goes-to-zero limit
Syllabus and review details

No official syllabus alignment is listed for this lesson.