Energy Transfer Between Coupled Oscillators
In 1665, confined to his room by a minor illness, Christiaan Huygens stared at two of his pendulum clocks hanging from the same wooden beam and noticed something eerie: no matter how he started them, within half an hour their pendulums settled into perfect anti-phase lockstep. He called it an "odd kind of sympathy." The clocks were talking to each other — tiny vibrations carried through the beam were shuttling energy back and forth between them. Huygens had discovered, two centuries early, the physics of energy transfer between coupled oscillators. Take any two identical oscillators — pendulums, masses on springs, LC circuits — and connect them weakly. Start one moving while the other rests. Then wait. The first oscillator slowly dies down while the second, untouched, swings up from nothing; then the flow reverses, and the energy sloshes back. The whole exchange is governed by a single small number: the frequency splitting between the system's two normal modes. Energy migrates completely in a time , and it migrates completely only if the two oscillators are perfectly in tune. That last clause is the deep one. The same mathematics — beating between two nearly degenerate modes — reappears as Rabi oscillations of a quantum two-level atom, the design equation of radio-frequency filters, energy hopping between vibrating bonds in a molecule, and the two-qubit gates of a superconducting quantum computer. Learn the pendulums and you have learned them all.
- Two identical oscillators joined by a weak coupling no longer oscillate independently: the true free motions are the symmetric and antisymmetric normal modes, with slightly different frequencies and .
- Starting all the energy in oscillator 1 excites both modes equally; the small splitting makes them drift out of phase, and the energy sloshes fully into oscillator 2 and back — beating.
- The masses wiggle fast, at the mean carrier frequency ; the amplitude envelope modulates slowly, at ; complete transfer takes .
- Total energy is conserved — only its partition oscillates: , .
- Complete transfer demands identical (resonant) oscillators. Detuning caps the transferred fraction at — the classical Rabi formula.
- The identical mathematics governs quantum Rabi flopping, coupled LC filters, molecular vibrational energy flow, and two-qubit iSWAP gates.
- Key idea: decompose the coupled motion into normal modes — the slow beat at the mode splitting carries the energy back and forth, with complete transfer only at exact resonance.
Decompose the coupled motion into normal modes; the slow beat at the mode-splitting frequency governs the back-and-forth energy flow, with complete transfer only at exact resonance.
Where this lives today
The purest modern incarnation of this article is sitting in a dilution refrigerator. In a superconducting quantum processor, two qubits coupled by a tunable element are exactly the detuned-oscillator problem (8): bring the qubits to resonance and excitation sloshes between them at the coupling rate; wait precisely half a slosh, , and you have executed an iSWAP gate — the energy-transfer beat, frozen at its midpoint, is a native two-qubit operation. Detune the qubits and the exchange shuts off as the Lorentzian (9) collapses — which is precisely how the gate is turned off. Gate calibration sheets at IBM, Google, and Rigetti are, in effect, measurements of and . The same physics runs at every scale. Wireless chargers push power through the resonant coupling of two coils, engineered so the "transfer time" beats the dissipation time — the condition of the limiting-cases table as a design spec. In photonics, two evanescently coupled waveguides swap light along the propagation direction (distance playing the role of time), the basis of directional couplers in every lithium-niobate modulator. And in chemistry, Fermi resonances shuttle vibrational energy between near-degenerate molecular modes, setting the speed of intramolecular energy redistribution that decides which bond of an excited molecule breaks first.
The physical picture
Hang two identical pendulums from a common support and join them with a soft, floppy spring. Pull pendulum 1 aside, hold pendulum 2 dead still, and release. For the first few swings nothing seems to happen — pendulum 1 swings, pendulum 2 barely quivers. But watch for a minute. Swing by swing, pendulum 1 loses amplitude and pendulum 2 gains it, until — startlingly — pendulum 1 hangs motionless and pendulum 2 has stolen the entire motion. Then the theft reverses. The energy is a fluid sloshing between two tanks connected by a thin pipe, and the narrowness of the pipe — the weakness of the coupling — sets the slow sloshing rate.
Why does the energy move at all? Because "pendulum 1 swinging, pendulum 2 still" is not a natural motion of the coupled system. The coupled pair has exactly two natural motions — its normal modes. In the symmetric mode both pendulums swing together, in phase, and the coupling spring between them never stretches; it might as well not exist, so this mode keeps the bare frequency . In the antisymmetric mode they swing in perfect opposition, working the spring hard on every cycle, which stiffens the restoring force and raises the frequency to . A localized excitation — all the motion in pendulum 1 — is secretly a 50/50 mixture of these two modes. And a mixture of two clocks ticking at slightly different rates cannot hold its shape: the modes drift apart in phase at the difference frequency , and where they initially added constructively in pendulum 1, half a beat later they add constructively in pendulum 2 instead. The pattern of energy has migrated even though neither mode has changed its amplitude at all. The result, traced in time, is the signature waveform of this whole subject — a fast carrier oscillation trapped inside a slow envelope:
Keep the two frequencies straight from the start, because conflating them is the classic error. The carrier — the visible wiggle — runs at the mean frequency , fast. The envelope — the slow swelling and dying of each oscillator's amplitude — runs at , slow. Energy, which goes as amplitude squared, completes its journey from oscillator 1 to oscillator 2 in half an envelope period: . Weaken the coupling and shrinks, the pipe narrows, and the slosh takes longer — but it still, eventually, completes. What can stop complete transfer is not weak coupling but detuning: make the two oscillators even slightly unequal and part of the energy stays locked in the first one forever. That resonance condition is the punchline of this article.
Two modes, one beat
This article rides on the normal-mode machinery built in Coupled oscillators — there you will find the general eigenvalue problem, and in Normal modes its full anatomy. Here we need only the simplest case: two equal masses , each on its own spring of stiffness , joined to each other by a weak spring . Newton's second law for each mass reads
where: — displacements from equilibrium; — the coupling force, proportional to the stretch of the connecting spring, pulling the masses toward each other's displacement. Each equation would be ordinary SHM if not for the last term, which feeds motion of one mass into the dynamics of the other. That cross-feed is the pipe through which energy will flow. The system decouples the moment you stop watching the masses and start watching their sum and difference. Define the normal coordinates
where: — the symmetric (center-of-mass-like) coordinate, blind to the coupling spring, which never stretches when the masses move together; — the antisymmetric (stretch) coordinate, which works the coupling spring twice over. Adding the two equations of motion kills the coupling terms; subtracting them doubles the coupling. Each combination obeys its own independent SHM, with frequencies
where: — the symmetric-mode frequency, exactly the bare because the coupling spring is a spectator; — the antisymmetric-mode frequency, raised because opposing motion stretches the coupling spring by and adds of effective stiffness. The mode that ignores the spring is slower; the mode that fights it is faster. Everything about energy transfer hides in how close these two frequencies are. For weak coupling, expand the square root to first order in :
where: — the mode splitting, linear in the coupling stiffness. Halve the coupling and you halve the splitting; the beat slows in exact proportion. This single small number is the metronome of the entire energy exchange — nothing else about the coupling survives into the final answer.
We are deriving the celebrated beating solution — , — for two identical weakly coupled oscillators released with all the motion in oscillator 1. The strategy is the standard normal-mode route: write the general motion as a superposition of the two independent mode oscillations, fix the four constants from the localized initial conditions, then reorganize the resulting sum of two cosines with sum-to-product identities so that the fast and slow time scales separate explicitly. The step to watch is the initial condition: a localized start forces the two modes to enter with equal amplitudes — that equality is precisely what makes later perfect cancellation (and hence complete transfer) possible, and it is exactly what detuning will destroy.
Step 1. Write the general solution as independent SHM in each normal coordinate:
Why: equations (1) added and subtracted give — two uncoupled oscillators, each with its own amplitude and phase to be fixed by initial data.
Step 2. Invert (2) to return to the physical coordinates:
Step 3. Impose the localized start: , , . In mode language with zero initial velocities, so and . Why the equal weights matter: "all the motion in oscillator 1" is not a mode — it is the fifty-fifty superposition of both modes. Each mode alone spreads its motion over both masses equally; only their interference localizes it.
Step 4. Substitute back:
Why the signs differ: mass 1 carries the modes in phase, mass 2 carries them with opposite signs — so wherever the two cosines agree, mass 1 moves and mass 2 stands still, and vice versa.
Step 5. Apply the sum-to-product identities with and :
where: the second factor in each product is the fast carrier at ; the first is the slow envelope at . The envelopes are a cosine and a sine of the same slow argument — when one is full the other is empty, in perfect complement. Verified: at (no coupling) the solution reduces to , — independent oscillators, as it must. Matches Taylor ch. 11 and Morin's coupled-pendulum result.
Equation (5) is worth reading slowly. Oscillator 1 starts at full amplitude () and dies to exactly zero when ; at that same instant oscillator 2's envelope reaches its maximum. The instant of total silence in oscillator 1 is the instant of total uproar in oscillator 2. Note also the carrier phases: oscillates as a sine against 's cosine — the receiving oscillator runs a quarter cycle behind the donor, which is the microscopic mechanism of the transfer: each fast cycle, oscillator 1's spring force does a little net positive work on oscillator 2 through the coupling, precisely because of that 90° lag.
The energy ledger
Go deeper · Derivation B · 1 minEnergy partition and its conservation
We are deriving the energy bookkeeping of the beat — , , with constant total — from the beating solution (5). The route is the two-timescale argument: because , each envelope is effectively frozen over any single fast cycle, so each oscillator behaves at every instant like an ordinary SHM whose amplitude happens to be the current envelope value; the SHM energy formula from Energy in harmonic motion then applies cycle by cycle. Watch the small print: these are cycle-averaged energies, and the tiny energy parked in the coupling spring (of order relative to the total) is neglected — both approximations are controlled by the weak-coupling condition.
Step 1. Read off the slowly varying amplitudes from (5):
Why this is legitimate: over one fast period the envelope argument advances by only — the envelope is constant to first approximation, so each mass executes clean SHM at the current amplitude.
Step 2. Apply the SHM energy formula to each oscillator, with for weak coupling:
Step 3. Insert the envelopes and define , the energy the system started with:
Step 4. Check the sum: , constant. Why the partition oscillates at , not : — squaring the envelope doubles its frequency, so the energy completes a full out-and-back cycle in and a one-way trip in half that. *Verified: at , and — complete transfer, consistent with Landau-Lifshitz §23-24.*
The one-way travel time is the number an experimenter actually measures with a stopwatch:
where: the second form uses the weak-coupling splitting (4). Read it as an impedance statement: stiffer coupling moves energy faster, in inverse proportion. For the classroom pendulum demo, is tens of seconds — dozens of fast swings — which is exactly what makes the effect so theatrical: the audience has time to notice that pendulum 1 has stopped. And note what is not conserved: nothing about oscillator 1 individually. Its energy runs from to zero and back with no half-measures. Only the total — and, separately, the energy in each normal mode — stays constant. Energy conservation lives in mode space; the drama in position space is pure interference.
Detuning, resonance, and the classical Rabi formula
Now break the symmetry. Let the two oscillators have slightly different bare frequencies, and — two pendulums of slightly unequal length, two LC circuits with mismatched capacitors. Everything above came from the equal-weight superposition forced by the localized start; detuning unbalances the weights, and the energy exchange becomes both faster and shallower:
The cleanest way to see it is to strip out the fast carrier entirely and follow only the envelopes. Write each oscillator's motion as , where the complex amplitudes vary slowly compared with the carrier — the same complex-exponential bookkeeping that tames every oscillation problem. Substituting into the coupled equations and dropping second derivatives of the slow envelopes gives
where: — the detuning between the bare frequencies; — the coupling rate, equal to for our spring-coupled pair at resonance. Anyone who has met the Schrödinger equation for a two-level atom should feel déjà vu: equation (8) is that equation, with the mode amplitudes as the "wavefunction" and the detuning and coupling as the level splitting and drive. Classical beats and quantum Rabi flopping are one calculation wearing two costumes.
Go deeper · Derivation C · 2 minIncomplete transfer under detuning
We are deriving the maximum fraction of energy that can ever reach oscillator 2 when the oscillators are detuned by — the classical Rabi formula . The route: solve the linear envelope system (8) with everything starting in oscillator 1, which is cleanest in terms of the system's eigenfrequency; the answer will show that detuning both speeds up the exchange and caps it below completeness. Watch the structure of the result — the cap is a Lorentzian in the detuning, the universal resonance line shape, and it does not care about the sign of : too long a pendulum blocks transfer exactly as effectively as too short.
Step 1. Look for normal solutions of the envelope system. The 2×2 eigenvalue problem gives
Why grows with detuning: the eigenvalues of the matrix with diagonal and off-diagonal sit at plus-or-minus the quadrature sum — detuning pushes the two envelope eigenfrequencies apart, so the exchange oscillation actually gets faster.
Step 2. Solve with the localized initial condition , . Projecting onto the two eigenvectors and recombining:
Why the prefactor appears: with detuning the two eigenvectors are no longer the symmetric and antisymmetric combinations — they lean toward the oscillators they are closer to in frequency. The localized start then projects unequally onto them, and unequal amplitudes can never fully cancel in oscillator 1 or fully add in oscillator 2.
Step 3. Read off the peak of the sine-squared:
where: — the largest fraction of the initial energy that ever visits oscillator 2, reached at . It is a Lorentzian of width in the detuning: match the oscillators to within their coupling rate and transfer is nearly complete; detune by much more and the energy barely leaks.
Step 4. Check the resonant limit: gives and — exactly the energy partition (6), with complete transfer at . Why this consistency check matters: it fixes the identification and confirms that the envelope formalism contains everything the exact solution knew. *Verified: identical in form to the two-level Rabi flopping formula; reduces to (6) at zero detuning; symmetric under as required.*
Which equation answers which question
| you want | use | notes |
|---|---|---|
| the two mode frequencies | , , eq. (3) | symmetric mode ignores the coupling spring |
| the beat/splitting rate | , eq. (4) | linear in coupling; the master small number |
| the full motion from a localized start | eq. (5) | carrier at , envelope at |
| energy in each oscillator vs time | , , eq. (6) | cycle-averaged; total constant |
| time for complete transfer | , eq. (7) | half the envelope period |
| exchange with detuning | and , eqs. (8)-(9) | faster but shallower; Lorentzian cap |
Limiting cases
| Case | Result | Meaning |
|---|---|---|
| weak coupling, identical oscillators | : clean beats, 100% transfer in | the textbook regime — many fast cycles per slosh |
| zero coupling | , | independent oscillators; energy never moves |
| strong coupling | : far apart | no two-timescale separation — just two distinct modes, no "beating" story |
| detuned oscillators | exchange faster () but forever incomplete | |
| light damping | many transfer cycles before decay | beats visible, slowly dying — see Damped harmonic oscillator |
| heavy damping | energy dissipates before one slosh completes | beating suppressed entirely |
Worked examples
1 · The classroom demo, by the numbers (confidence builder)
Two carts of mass each sit on separate springs of stiffness and are joined by a weak spring . Cart 1 is displaced and released with cart 2 at rest. Find the two mode frequencies, the beat splitting, the complete-transfer time, and the number of fast oscillations per transfer.
Step 1. Bare and symmetric-mode frequency: .
Step 2. Antisymmetric mode: .
Step 3. Splitting, exact and approximate: ; the weak-coupling formula gives — within 2%, confirming the linearization.
Step 4. Transfer time: .
Step 5. Fast cycles per transfer: the carrier period is , so the energy needs to change hands — slow enough to watch, fast enough for a lecture.
2 · Detuning kills the demo (exam level)
Two carts on springs, and each, coupled by . (a) Find the transfer time when the carts are identical. (b) A student swaps cart 2 for one of mass — a 5% error. What fraction of the energy now reaches cart 2, and when does the exchange peak?
Step 1. Identical case: , , so and .
Step 2. The detuning: cart 2's new bare frequency is , so . Only the magnitude matters. A 5% mass error has produced a detuning of the same size as the whole coupling effect.
Step 3. Coupling rate: at resonance the exchange formula must reproduce Step 1, which fixes .
Step 4. Transferred fraction from (9): — a fifth of the energy never leaves cart 1, from a 5% mass mismatch.
Step 5. Peak time: , so the maximum occurs at — earlier than the resonant 10.2 s. Detuning speeds the exchange up while making it incomplete, exactly the Lorentzian trade-off of equation (9).
Origins · 1 minEnrico Fermi and how this was discovered
Fermi appears in this story for what may be the most consequential "failed" calculation in mechanics. In 1953, with Pasta, Ulam, and programmer Mary Tsingou at Los Alamos, he set the MANIAC computer loose on a chain of 64 oscillators with weak nonlinear coupling, expecting the energy — started in one mode — to spread and equipartition among all the modes, the founding assumption of statistical mechanics. It refused. The energy sloshed among a handful of low modes and then, astonishingly, reassembled almost perfectly in the initial mode: recurrence, not thermalization. The Fermi-Pasta-Ulam problem, published after Fermi's death in 1955, is the direct many-body descendant of this article's two-oscillator slosh — and its resolution seeded two fields, solitons and chaos theory (via the KAM theorem). Fermi called it "a minor discovery." It was the first numerical experiment in the history of physics.
Going further
Texts. Taylor, Classical Mechanics, ch. 11 (Coupled Oscillators and Normal Modes) · Morin, Introduction to Classical Mechanics, coupled-oscillators chapter · French, Vibrations and Waves, ch. 5 · Goldstein, Poole & Safko, Classical Mechanics, 3rd ed., ch. 6 · Landau & Lifshitz, Mechanics, §23-24.
Provenance. Verifier should rederive the normal-mode frequencies , and the splitting from the coupled equations (Taylor ch. 11, Morin coupled-oscillators chapter); confirm the beating solution and energy partition against Landau-Lifshitz §23-24; and check the detuned transfer fraction against the standard Rabi formula. All expressions are symbolic textbook results; no fitted numbers used. Worked-example numbers on this page ( and s for the carts; and s for the detuned pair) are illustrative arithmetic recomputed from the boxed formulas.
Builds on: Small oscillations about equilibrium · Simple harmonic motion · Amplitude, phase, and angular frequency · Energy in simple harmonic motion · Complex exponential representation of oscillations. Leads to: Resonance · Driven oscillations · Anharmonic pendulum and amplitude dependence · Damped harmonic oscillator.
Why does starting only oscillator 1 excite both normal modes, and why equally?
Because "oscillator 1 displaced, oscillator 2 at rest" is not a shape the system can hold — each normal mode is a pattern spread over both masses. To write the localized start in mode language you must solve and : both modes enter with amplitude , exactly equal. Physically, the symmetric mode wants to move both masses together and the antisymmetric wants to move them oppositely; only the fifty-fifty sum moves mass 1 while cancelling at mass 2. The equality is not a coincidence of this example — any strictly localized excitation of two identical coupled oscillators projects equally onto the two modes, and that equal weighting is precisely what allows the later perfect cancellation that empties oscillator 1 completely.
Is the beat frequency δ or δ/2? Different books seem to disagree
Both, depending on what you track — this is a bookkeeping convention, not a physics disagreement. The amplitude envelope of oscillator 1 is , oscillating at angular frequency . But the energy goes as the envelope squared, , which oscillates at — squaring doubles the frequency because the energy cannot tell a positive envelope from a negative one. What the ear hears in acoustic beats, and what the eye sees in the pendulum demo, is the loudness — the energy — so the perceived beat rate is , and one complete one-way energy transfer takes . Books quoting are describing the envelope; books quoting are describing what you observe.
Where exactly does the energy sit while it is in transit between the oscillators?
Almost all of it is always in one oscillator or the other — the coupling spring itself holds only a fraction of order of the total, negligible by the weak-coupling assumption. Mid-transfer, when , each oscillator simply holds half the energy: . The transfer mechanism is ordinary work: oscillator 2 lags oscillator 1 by a quarter carrier cycle (the against in equation (5)), and with that phase relationship the coupling force on mass 2 points along mass 2's velocity more often than against it — positive net work every cycle. When the energy flows back, the lag reverses. There is no mysterious reservoir; just a weak force with the right phase, applied patiently over many cycles.
How is this "the classical analogue of Rabi oscillations"? What maps to what?
Write the slowly varying envelope equations (8) next to the Schrödinger equation for a two-level atom in a resonant field and they are term-for-term identical: the envelope pair plays the two-level amplitudes, the detuning plays the atom-field frequency mismatch, and the coupling rate plays the Rabi frequency. "Probability of finding the atom excited" becomes "fraction of energy in oscillator 2," and both obey . The analogy is exact because both systems are two coupled modes evolving linearly. What does not map is interpretation: the classical is a definite energy fraction you could measure continuously, while the quantum version is a probability revealed only on measurement. The mathematics is shared; the ontology is not.
Do the beats survive damping, and how many transfers can I actually see in a real demo?
The criterion is the race between the damping rate and the splitting . Each oscillator loses energy as (amplitude as ), so the number of visible transfer cycles is roughly . Lecture pendulums with and a beat period of thirty seconds show two or three clean sloshes before the whole picture fades. If the energy dies before completing even one trip — the beating regime disappears, which is why you cannot demonstrate this with heavily damped oscillators like door closers. Engineers state the same condition positively: for wireless power transfer to work, the coupling rate must beat the loss rate, — the "strong coupling regime" that every cavity-QED and optomechanics paper advertises is exactly this inequality.
Can I speed up the transfer by making the coupling spring stiffer? What breaks if I overdo it?
Yes — , so doubling halves the transfer time. What breaks at large is the story, not the mathematics. The two-timescale picture — fast carrier, slow envelope — requires , i.e. . Push toward and the "beat" becomes as fast as the oscillation itself: the system no longer looks like energy sloshing between two recognizable oscillators, but simply like two very different normal modes doing their own thing. The exact normal-mode solution remains valid at any coupling strength; it is the interpretation — "the energy is in oscillator 1 now" — that loses meaning when the modes are far apart, because on any timescale you can define, both masses are moving comparably.
My two "identical" lab pendulums never hand over all the energy. How close is close enough?
Equation (9) is your tolerance spec. Complete transfer needs the detuning small compared with the coupling: requires , i.e. the bare frequencies matched to well within the mode splitting itself. Since is maybe a percent in a typical demo, you need the pendulum frequencies matched to a fraction of a percent — for , lengths matched to a few parts in a thousand, millimetres on a metre. That is why the demo is fussy: worked example 2 shows a 5% mass mismatch already trapping 20% of the energy. The practical recipe is to tune one pendulum's length while watching the first oscillator's minimum amplitude and minimizing it — you are experimentally zeroing using as the signal.
Does the energy REALLY vanish from oscillator 1 at the beat node, or is a little always left?
For truly identical, undamped oscillators the zero is exact — at the envelope passes through zero and oscillator 1 is instantaneously at rest at its equilibrium point with zero velocity: no kinetic, no potential energy. (Its envelope then goes negative, which just means the returning oscillation is phase-flipped.) In any real experiment three things blur the perfect null: residual detuning caps the transfer via and leaves oscillator 1 wiggling at the node; damping shrinks the whole picture; and the small energy resident in the coupling spring, of order , never belongs to either oscillator. But these are perturbations. The interference argument — two equal-amplitude modes in antiphase at mass 1 — cancels exactly, and the cleaner the experiment, the deader the pendulum at the node.