Two coupled oscillators

Push one swing in a playground and only that swing moves. Now imagine two swings joined by a stretchy rope — push one, and the rope drags the other into motion too, and soon the push you gave to swing 1 has mysteriously ended up making swing 2 move while swing 1 sits nearly still. Neither swing is doing plain simple harmonic motion any more; they are trading energy back and forth. This is the world of coupled oscillators: two identical masses on springs, joined by a third spring, is the simplest system that shows it, and it is the gateway to normal modes, the idea that any messy coupled motion is secretly just two clean, independent oscillations added together. This article rebuilds the two-mass-two-spring model from Newton's second law, shows the algebraic trick — adding and subtracting the two equations of motion — that decouples them, and follows a single worked configuration all the way from spring constants to normal-mode frequencies to the actual x1(t),x2(t)x_1(t), x_2(t) trajectories, including the striking beat pattern that appears when only one mass is disturbed initially.

The core insight

Two masses linked by a coupling spring obey coupled equations of motion that decouple exactly into two independent SHM equations for the sum and difference coordinates q1=x1+x2, q2=x1−x2q_1=x_1+x_2,\ q_2=x_1-x_2, whose frequencies ω1=k/m\omega_1=\sqrt{k/m} and ω2=(k+2kc)/m\omega_2=\sqrt{(k+2k_c)/m} are the system's normal modes.

The core idea

Setting up the system. Take two identical masses mm, each held to a wall by its own spring of stiffness kk, and joined to each other by a coupling spring of stiffness kck_c. Let x1,x2x_1, x_2 be their displacements from equilibrium. The coupling spring's stretch is x1−x2x_1 - x_2; it pulls mass 1 back and pushes mass 2 forward by the same amount, so Newton's second law on each mass gives

mx¨1+(k+kc)x1−kcx2=0,mx¨2+(k+kc)x2−kcx1=0m\ddot x_1 + (k+k_c)x_1 - k_c x_2 = 0, \qquad m\ddot x_2 + (k+k_c)x_2 - k_c x_1 = 0

where: mm is each mass, kk the stiffness of the two identical end springs, kck_c the stiffness of the coupling spring, and x1,x2x_1, x_2 the displacements of mass 1 and mass 2 from equilibrium. Each equation contains both coordinates — that is exactly what "coupled" means, and it is why neither mass on its own executes SHM.

Decoupling by symmetry. The system has a left–right symmetry (swap mass 1 and mass 2 and the equations swap too), and that symmetry is the clue: add the two equations, then subtract them.

m(x¨1+x¨2)=−k(x1+x2),m(x¨1−x¨2)=−(k+2kc)(x1−x2)m(\ddot x_1+\ddot x_2) = -k(x_1+x_2), \qquad m(\ddot x_1-\ddot x_2) = -(k+2k_c)(x_1-x_2)

Naming q1=x1+x2q_1 = x_1+x_2 and q2=x1−x2q_2 = x_1-x_2 turns these into two independent SHM equations, one per coordinate:

q¨1+ω12q1=0,  ω1=k/mq¨2+ω22q2=0,  ω2=(k+2kc)/m\ddot q_1+\omega_1^2 q_1=0,\ \ \omega_1=\sqrt{k/m}\qquad\qquad \ddot q_2+\omega_2^2 q_2=0,\ \ \omega_2=\sqrt{(k+2k_c)/m}

where: q1,q2q_1, q_2 are the normal coordinates, and ω1,ω2\omega_1, \omega_2 are the normal-mode angular frequencies. Physically, ω1=k/m\omega_1=\sqrt{k/m} is the frequency of the mode where both masses move together, in phase (the coupling spring never stretches, so kck_c drops out); ω2=(k+2kc)/m\omega_2=\sqrt{(k+2k_c)/m} is the mode where the masses move in antiphase, mirror-imaging each other (the coupling spring stretches maximally, stiffening the restoring force, hence the higher frequency — always ω2>ω1\omega_2>\omega_1 for kc>0k_c>0).

Recovering the real motion. q1,q2q_1,q_2 are pure cosines, q1(t)=C1cos⁡(ω1t+ϕ1)q_1(t)=C_1\cos(\omega_1 t+\phi_1) and q2(t)=C2cos⁡(ω2t+ϕ2)q_2(t)=C_2\cos(\omega_2 t+\phi_2), and inverting the definitions gives the physical displacements as an equal mixture of both modes:

x1(t)=12[q1(t)+q2(t)],x2(t)=12[q1(t)−q2(t)]x_1(t)=\tfrac12\big[q_1(t)+q_2(t)\big], \qquad x_2(t)=\tfrac12\big[q_1(t)-q_2(t)\big]

The four constants C1,ϕ1,C2,ϕ2C_1,\phi_1,C_2,\phi_2 come from the four initial conditions x1(0),x2(0),x˙1(0),x˙2(0)x_1(0),x_2(0),\dot x_1(0),\dot x_2(0). Unless the system is launched into a single pure mode, each mass's motion is a superposition of two frequencies, not SHM.

Figure
Two masses on springs attached to walls and to each other by a coupling spring showing the in phase and antiphase normal modes
Two normal modes of the coupled system mode 1: in phase, ω₁ = √(k/m) mode 2: antiphase, ω₂ = √((k+2k_c)/m) mode 1 → → m m mode 2 → ← m m both masses move the same way together, coupling spring never stretches so ω₁ does not depend on k_c masses move in mirror image directions, coupling spring stretches most, raising ω₂ above ω₁
Derivation ADecoupling the two equations of motion into normal modes

Starting from Newton's second law applied separately to each of the two coupled masses, show that the two simultaneous, mutually dependent differential equations for x1(t)x_1(t) and x2(t)x_2(t) can be converted, by a simple linear change of coordinates, into two completely independent simple-harmonic-oscillator equations, and read off the resulting normal-mode frequencies ω1\omega_1 and ω2\omega_2 in terms of mm, kk, and kck_c.

Step 1. Write Newton's second law for each mass. The coupling spring stretches by x1−x2x_1-x_2, exerting force −kc(x1−x2)-k_c(x_1-x_2) on mass 1 and the equal-and-opposite +kc(x1−x2)+k_c(x_1-x_2) on mass 2, alongside each mass's own end spring:

mx¨1=−kx1−kc(x1−x2),mx¨2=−kx2+kc(x1−x2)m\ddot x_1 = -kx_1-k_c(x_1-x_2), \qquad m\ddot x_2=-kx_2+k_c(x_1-x_2)

Step 2. Rearrange both into standard form, collecting the coefficients of x1x_1 and x2x_2 on each line:

mx¨1+(k+kc)x1−kcx2=0,mx¨2+(k+kc)x2−kcx1=0m\ddot x_1+(k+k_c)x_1-k_cx_2=0, \qquad m\ddot x_2+(k+k_c)x_2-k_cx_1=0

Both equations mix x1x_1 and x2x_2, so this pair cannot be solved one line at a time — they are genuinely coupled.

Step 3. Add the two equations. The cross terms −kcx2-k_cx_2 and −kcx1-k_cx_1 combine with the like terms to give

m(x¨1+x¨2)+k(x1+x2)=0m(\ddot x_1+\ddot x_2)+k(x_1+x_2)=0

which involves only the combination x1+x2x_1+x_2.

Step 4. Subtract the second equation from the first. The −kx1-kx_1 and −kx2-kx_2 terms combine into −k(x1−x2)-k(x_1-x_2), and the coupling terms −kcx2-k_cx_2 and +kcx1+k_cx_1 both add constructively:

m(x¨1−x¨2)+(k+2kc)(x1−x2)=0m(\ddot x_1-\ddot x_2)+(k+2k_c)(x_1-x_2)=0

which involves only the combination x1−x2x_1-x_2.

Step 5. Define q1=x1+x2q_1=x_1+x_2 and q2=x1−x2q_2=x_1-x_2. Steps 3 and 4 become mq¨1+kq1=0m\ddot q_1+kq_1=0 and mq¨2+(k+2kc)q2=0m\ddot q_2+(k+2k_c)q_2=0 — two completely separate SHM equations, each in only one coordinate. Reading off ω2=(coefficient)/m\omega^2 = (\text{coefficient})/m gives the normal-mode frequencies ω1=k/m\omega_1=\sqrt{k/m} and ω2=(k+2kc)/m\omega_2=\sqrt{(k+2k_c)/m}, and since q1,q2q_1,q_2 are independent, each obeys its own separate initial-condition-fixed cosine solution.

Worked example

Two identical masses m=0.50 kgm=0.50\ \text{kg} are connected by end springs k=20 N/mk=20\ \text{N/m} and a coupling spring kc=5.0 N/mk_c=5.0\ \text{N/m}. The system is released from rest with mass 1 displaced by x1(0)=A0=0.02 mx_1(0)=A_0=0.02\ \text{m} and mass 2 undisturbed, x2(0)=0x_2(0)=0. Find the normal-mode frequencies and the resulting x1(t),x2(t)x_1(t), x_2(t).

Step 1. Normal-mode frequencies follow directly from the formulas derived above:

ω1=k/m=20/0.5=40≈6.32 rad/s,ω2=(k+2kc)/m=30/0.5=60≈7.75 rad/s\omega_1=\sqrt{k/m}=\sqrt{20/0.5}=\sqrt{40}\approx 6.32\ \text{rad/s}, \qquad \omega_2=\sqrt{(k+2k_c)/m}=\sqrt{30/0.5}=\sqrt{60}\approx 7.75\ \text{rad/s}

Step 2. Convert the initial conditions to normal coordinates: q1(0)=x1(0)+x2(0)=A0q_1(0)=x_1(0)+x_2(0)=A_0 and q2(0)=x1(0)−x2(0)=A0q_2(0)=x_1(0)-x_2(0)=A_0. Released from rest means x˙1(0)=x˙2(0)=0\dot x_1(0)=\dot x_2(0)=0, so q˙1(0)=q˙2(0)=0\dot q_1(0)=\dot q_2(0)=0 too — each normal coordinate starts at maximum displacement with zero velocity, so both are pure cosines with zero phase:

q1(t)=A0cos⁡ω1t,q2(t)=A0cos⁡ω2tq_1(t)=A_0\cos\omega_1 t, \qquad q_2(t)=A_0\cos\omega_2 t

Step 3. Invert to the physical displacements using x1=12(q1+q2)x_1=\tfrac12(q_1+q_2), x2=12(q1−q2)x_2=\tfrac12(q_1-q_2):

x1(t)=A02[cos⁡(6.32t)+cos⁡(7.75t)],x2(t)=A02[cos⁡(6.32t)−cos⁡(7.75t)]x_1(t)=\frac{A_0}{2}\big[\cos(6.32t)+\cos(7.75t)\big], \qquad x_2(t)=\frac{A_0}{2}\big[\cos(6.32t)-\cos(7.75t)\big]

with A0=0.02 mA_0=0.02\ \text{m} and tt in seconds.

Step 4. Using the sum-to-product identity, x1(t)=A0cos⁡ ⁣(ω2−ω12t)cos⁡ ⁣(ω1+ω22t)x_1(t)=A_0\cos\!\left(\tfrac{\omega_2-\omega_1}{2}t\right)\cos\!\left(\tfrac{\omega_1+\omega_2}{2}t\right) — a fast carrier oscillation at the average frequency ≈7.03 rad/s\approx 7.03\ \text{rad/s}, modulated by a slow envelope at half the frequency difference, ≈0.71 rad/s\approx 0.71\ \text{rad/s}. This is the beat phenomenon: energy sloshes from mass 1 to mass 2 and back, with mass 1's amplitude dying as mass 2's grows, at beat frequency ∣ω2−ω1∣≈1.43 rad/s|\omega_2-\omega_1|\approx 1.43\ \text{rad/s}.

Common traps

  • Assuming each mass individually undergoes SHM. It does not, in general — only the normal coordinates q1,q2q_1,q_2 are simple harmonic; x1(t)x_1(t) and x2(t)x_2(t) are sums of two different frequencies unless a single pure mode is excited.
  • Getting the coupling-spring sign wrong: it must pull mass 1 and push mass 2 in opposite senses, −kc(x1−x2)-k_c(x_1-x_2) and +kc(x1−x2)+k_c(x_1-x_2) — using the same sign in both equations of motion breaks Newton's third law and gives unphysical results.
  • Confusing normal coordinates with physical displacements. q1,q2q_1,q_2 are algebraic constructs chosen to decouple the maths; only x1,x2=12(q1±q2)x_1,x_2=\tfrac12(q_1\pm q_2) are what an observer actually measures moving.
  • Forgetting that ω2>ω1\omega_2>\omega_1 is a consequence of this particular spring layout (k+2kc>kk+2k_c>k), not a universal law — a different way of attaching the coupling spring changes the formula for ω2\omega_2.

Ask the tutor: if you want to see how the setup changes for a coupling spring attached differently, or want the initial-condition algebra worked through for a different starting displacement, ask the AI tutor to walk through it step by step.

Why do we bother defining normal coordinates instead of just solving the two coupled equations directly?

The two original equations for x1,x2x_1,x_2 each contain both variables, so you cannot solve one without the other — there is no clean "integrate this alone" step. The combinations q1=x1+x2q_1=x_1+x_2 and q2=x1−x2q_2=x_1-x_2 are chosen precisely because adding and subtracting the equations makes the cross terms cancel or double up cleanly, leaving two completely separate, ordinary SHM equations you already know how to solve. It is a coordinate trick, not new physics.

Why is the antiphase frequency ω2 always higher than the in phase frequency ω1?

In the in-phase mode both masses move together by the same amount at every instant, so the coupling spring's length never changes and contributes zero extra restoring force — only the end springs act, giving ω1=k/m\omega_1=\sqrt{k/m}. In the antiphase mode the masses move oppositely, so the coupling spring stretches by twice as much as either mass's displacement, adding extra stiffness 2kc2k_c to the effective restoring force and pushing the frequency up to ω2=(k+2kc)/m\omega_2=\sqrt{(k+2k_c)/m}, which is always greater than ω1\omega_1 whenever kc>0k_c>0.

Does the beat pattern mean energy is being lost from the system?

No — total mechanical energy is conserved throughout; nothing is dissipated because there is no damping in this idealised model. What happens is that energy is continuously exchanged between the two masses: as mass 1's oscillation amplitude shrinks, mass 2's grows correspondingly, and the sum of their energies stays constant. The beat envelope you see in x1(t)x_1(t) and x2(t)x_2(t) is this energy transfer made visible, at the slow beat frequency ∣ω2−ω1∣|\omega_2-\omega_1|.

Before this

Where this leads