Chapter 10

Case study: phase-locked loops

Open any book on frequency standards and you meet the same figure: boxes labeled oscillator, multiplier, detector, arrows chasing each other in a ring. This chapter is about reading those diagrams — what each box is for, and how the closed ring reshapes every noise you met in chapters 1–9. The diagrams here are alive: click the blocks.

Why lock anything?

Because no single oscillator has every virtue. A quartz crystal is magnificently quiet close to the carrier but comes at one fixed, low frequency — and multiplying it up to microwave frequencies multiplies its phase noise too, ruining its floor far from the carrier. A microwave VCO is tunable and clean far out, but its phase wanders badly close in. A phase-locked loop welds the two into a single output that follows the reference below a chosen loop bandwidth and the VCO above it: each oscillator is trusted exactly where it is good. That one sentence is the whole chapter; the demo in the noise-shaping section below lets you find the best weld point yourself.

Anatomy of a phase-locked loop

Four blocks and a ring of arrows. The VCO's output is compared against the reference in a phase detector; the phase error is filtered and fed back to the VCO's tuning voltage so that the error shrinks. Optionally a $\div N$ divider sits in the feedback path so a low reference frequency can discipline a high output frequency. Click each block (or tab to it) to see its job and its signal:

Interactive block diagram: the generic PLL

Click a block to highlight it and read what it is for; the mini-plot shows that block's characteristic signal.

output error control V Reference 100 MHz crystal Phase detector V = Kd·Δφ Loop filter F(f) VCO K₀ Hz/V ÷N divider

Noise shaping: the loop bandwidth decides whom to trust

Now the mathematics, compactly. Around the ring, phase error is converted to volts ($K_d$), filtered ($F$), and turned back into phase by the VCO — which, being an integrator from volts to phase, contributes a factor $1/(if)$. The open-loop gain is therefore

$$ G(f) \;=\; \frac{K}{i f}\, F(f), \qquad K = K_d K_0 \ (\times\, 1/N \text{ with a divider}), $$

and standard feedback algebra (the same $1/(1+G)$ bookkeeping as chapter 9's servo section) gives the closed-loop output phase as a weighted blend of the two oscillators' phases:

$$ \varphi_{\mathrm{out}} \;=\; H\,\varphi_{\mathrm{ref}} \;+\; (1 - H)\,\varphi_{\mathrm{vco}}, \qquad H = \frac{G}{1+G}. $$

To dock the notation explicitly onto chapter 9: its disturbance suppression $1/(1+G)$ is this chapter's $1-H$, and its sensor-noise injection $G/(1+G)$ is this chapter's $H$ — same algebra, with the "sensor" here being the reference-plus-phase-detector chain. Because of the VCO's built-in $1/f$, $|G|$ is large at low $f$ and small at high $f$: $H$ is a low-pass (the loop follows the reference, but only slowly) and $1-H$ is a high-pass (the loop cannot fix fast VCO wiggles). Note the flip side: inside the bandwidth the loop writes the reference's phase noise — and the divider's, and the phase detector's — straight onto the VCO. For the workhorse integrator-plus-lead ("type 2") loop filter the closed-loop response takes the textbook second-order form (with $s = i\,2\pi f$ and $x = f/f_n$):

$$ H(s) = \frac{2\zeta\omega_n s + \omega_n^2}{s^2 + 2\zeta\omega_n s + \omega_n^2}, $$

$$ |H|^2 = \frac{1 + 4\zeta^2 x^2}{(1-x^2)^2 + 4\zeta^2 x^2}, \qquad |1-H|^2 = \frac{x^4}{(1-x^2)^2 + 4\zeta^2 x^2}, $$

with natural frequency $f_n \approx$ the loop bandwidth $f_L$ and damping $\zeta$. The composite spectrum is then computed exactly like a chapter-9 noise budget — in linear power units, never in dB:

// closed-loop composite (linear units, then 10·log10 for dBc/Hz)
const x = f / fn, z2 = 4 * zeta * zeta * x * x;
const den = (1 - x * x) ** 2 + z2;
const H2 = (1 + z2) / den;                       // |H|²  : low-pass
const E2 = x ** 4 / den;                         // |1−H|²: high-pass
const Ltot = N * N * Lref(f) * H2 + Lvco(f) * E2;
The headline demo: welding two oscillators with one knob

A 100 MHz crystal, multiplied by $N$ to the 10 GHz carrier, locks a 10 GHz VCO. Drag the loop bandwidth $f_L$ and watch the closed-loop total trade the (multiplied) reference for the VCO. The RMS phase error is minimized when $f_L$ sits near the crossing of the two free-running curves; the readout computes it as $\varphi_{\mathrm{rms}} = \sqrt{\smash{\int} S_\varphi\,df} = \sqrt{\smash{2\int} \mathcal{L}(f)\,df}$, using chapter 6's $\mathcal{L} = S_\varphi/2$ in reverse. Lower the damping $\zeta$ and a servo bump grows at $f_L$ — chapter 9's peaking, in a PLL.

Three discoveries the sliders hand you. Too narrow: the VCO's close-in wander leaks through $|1-H|^2$ and the phase error explodes — at $f_L = 100\,$Hz it exceeds a radian, and a real loop would be slipping cycles. Too wide: you copy the multiplied reference's noise floor, $20\log_{10}N$ above the bare crystal's, out to offsets where the VCO was already quieter. Low $\zeta$: the loop resonates at $f_n$ and $|H|$ peaks (about $+5.8\,$dB at $\zeta = 0.3$), stacking a bump onto the output — when you see a bump on a beat note, you are usually looking at somebody's loop bandwidth.

The 20 log N tax

Why does multiplying a beautiful crystal wreck its floor? Phase is counted in cycles, and a frequency multiplier by definition produces $N$ output cycles per input cycle — so every phase excursion, deterministic or random, is multiplied by $N$ too: $\varphi \to N\varphi$, hence $S_\varphi \to N^2 S_\varphi$, hence

$$ \mathcal{L}(f) \;\longrightarrow\; \mathcal{L}(f) + 20\log_{10} N \quad \text{(everywhere, at every offset } f\text{)}. $$

This is the rule chapter 6 leaned on when it declared frequency multiplication "noise-neutral" in Leeson's model ($\mathcal{L} \propto \nu_0^2$ exactly cancels the $+20\log_{10}N$) — here is where the tax itself comes from. Worked number: a 5 MHz OCXO with $\mathcal{L} = -120$ dBc/Hz at 1 kHz, multiplied to 10 GHz ($N = 2000$, $20\log_{10}2000 \approx 66$ dB), arrives at $-54$ dBc/Hz. This single line is why microwave synthesis is hard, and it is exactly what the $N$ selector in the demo above does: the whole reference curve rides up by $20\log_{10}N$, dragging the optimal loop bandwidth down with it. The converse is a gift: frequency dividers reduce phase noise by $20\log_{10}N$ — down to the divider's own additive noise floor, which is why a good divider is one of the quietest boxes in a synthesizer.

Lab note: real loops have dirt
The blocks themselves add noise, and inside the bandwidth the loop faithfully writes it onto the output. A good analogue (mixer) phase detector has an additive floor near $-150$ to $-160$ dBc/Hz; digital phase-frequency detectors are noisier and add flicker. Expect reference spurs — bright coherent sidebands (chapter 2) at the phase-comparison frequency and its harmonics, from detector ripple leaking into the VCO tune line. And every spur entering the loop is shaped by the same $H$ or $1-H$ as the noise: spurs on the reference are copied inside the bandwidth, spurs on the VCO are suppressed there. An in-loop measurement of any of this is subject to chapter 8's caveat: the error signal flatters the loop; witness the output against an independent reference.

An atomic clock is a lock loop

Now the classic diagram this chapter promised to teach you to read: the atomic frequency standard. A quartz oscillator is multiplied up to the caesium hyperfine transition near 9.19 GHz, the atoms are interrogated, and a control voltage steers the quartz so the interrogation stays centered on the atomic line. The output of the clock is the quartz — the atoms never emit your signal; they only referee it. Click the blocks:

Interactive block diagram: the atomic frequency standard

The same five roles as the PLL above, wearing atomic clothes. Click each block; the modulator + synchronous detector pair is the trick worth understanding, and it gets its own demo below.

ω ≟ ω₀ ω₀/2π = 9 192 631 770 Hz error signal V clock output sin Ωt, Ω = ω/M Quartz VCXO Ω = Ω(V) ×M multiplier Modulator Atomic resonator Q = ω₀/Δω  ·  S/N Synchronous detector Integrator loop filter dither osc fₘ
Is this a phase-locked loop?
Structurally, yes: a steered oscillator, an error signal, feedback, and exactly the same $H$ / $1-H$ noise shaping — inside the loop bandwidth the output inherits the discriminator (atoms + detection), outside it the quartz flywheel. Strictly, though, it is a frequency-lock loop: the dithered atomic line senses the frequency offset $\nu - \nu_0$ from the transition, not an accumulated phase difference against a running reference oscillator. The loop algebra of this chapter carries over with $\varphi$ replaced by $\nu$.

The S-curve: how a symmetric line steers a loop

The atomic resonance is a symmetric curve: its height tells you how far $|\,\nu - \nu_0|$ you are from center, but not on which side — useless for feedback, which needs a signed error. The fix is chapter 8's lock-in trick: dither the interrogation frequency across the line at $f_m$ and demodulate the atomic response at $f_m$. Sitting left of center, a positive dither swing increases the signal; sitting right of center, it decreases it — the component at $f_m$ flips sign as you cross the line. The demodulated output versus detuning is an antisymmetric S-curve, linear through zero: a proper error signal, conjured out of a symmetric bump.

Dither, demodulate, lock: the S-curve at work

Drag the detuning $\delta$ and watch the modulated response (top right) flip phase as you cross the line; the S-curve (bottom left) is its honestly computed first harmonic. Then tick "close the loop" and press run lock: $\delta(t)$ relaxes to the line center against noise. Untick it while running to watch the free-running quartz wander off. Start the lock from $|\delta| > 1.5$ linewidths and notice how weakly the loop pulls — the S-curve dies outside the line: a capture range.

Play with the dither depth. The center slope — the discriminator sensitivity, in chapter 8's language literally volts per hertz — is greatest when the dither swings by about half a linewidth (the numerical optimum here is an amplitude $a \approx 0.35\,\Delta\nu$, i.e. swinging roughly between the half-maximum points; anything from $0.2$ to $0.6\,\Delta\nu$ is within about 20% of it). Too small a dither and the modulated signal vanishes into detector noise; too large and the slope collapses — you are averaging over the whole line, effectively broadening it.

And here the whole course lands in one formula. Within the loop bandwidth, the clock output's fractional frequency stability is set by the atoms and the detection SNR, not by the quartz:

$$ \sigma_y(\tau) \;\approx\; \frac{\Delta\nu}{\nu_0}\, \frac{1}{\mathrm{SNR}}\, \tau^{-1/2}. $$

The discriminator splits the line to $\delta\nu \approx \Delta\nu/\mathrm{SNR}$ in one second (chapter 8: noise referred to the input via the S-curve slope), and successive interrogations average down like white frequency noise — the $\tau^{-1/2}$ regime you met on every Allan-deviation plot in chapter 7 is this loop at work. High $Q = \omega_0/\Delta\omega$ and high SNR: that is the entire recipe for an atomic clock, and the reason the definition of the second lives in a box like this.

The same diagram everywhere

Once you can read one of these rings, you can read them all. Five roles, four labs:

roleRF synthesizeratomic clock PDH laser lockLIGO arm lock
flywheel oscillatorVCOquartz VCXO laserlaser + test masses
frequency translation÷N divider ×M multiplierEOM RF sidebands RF sidebands on the light
discriminatormixer phase detector dither + synchronous detector reflection photodiode + mixer photodiode + demodulation
referencecrystal oscillator atomic transitionoptical cavity arm-cavity resonance
integratorloop filterintegrator servo amplifiersuspension actuators + servo

The moral: the loop algebra — $H$, $1-H$, the servo bump, the $20\log_{10}N$ bookkeeping — is identical in every column. All the physics hides in the discriminator: its sensitivity (slope, V/Hz) and its noise floor set what the lock can possibly deliver, exactly as the atomic line's $Q$ and SNR set a clock's $\sigma_y(\tau)$.

Exercises

Exercise 10.1 — quartz to caesium

A 5 MHz quartz oscillator is multiplied to the caesium transition at 9 192 631 770 Hz. What is $M$, and by how many dB does the phase noise rise? If the quartz shows $-130$ dBc/Hz at 1 kHz offset, what does the multiplied signal show there? A decent 9 GHz dielectric-resonator oscillator (DRO) manages $-95$ dBc/Hz at 1 kHz — qualitatively, where should the loop bandwidth go if you lock the DRO to the multiplied quartz?

Solution

$M = 9.192631770\times10^9 / 5\times10^6 \approx 1839$, so $20\log_{10}M \approx 65$ dB. The multiplied quartz sits at $-130 + 65 = -65$ dBc/Hz at 1 kHz — a full 30 dB worse than the DRO there. But close to the carrier (below $\sim$10–100 Hz) the multiplied quartz still wins, because the DRO's $1/f^3$ wander keeps climbing. So the loop bandwidth belongs near the crossing of the two curves: low, of order tens to a few hundred Hz — trust the quartz only very close in, and let the DRO speak everywhere else.

Exercise 10.2 — the VCO is an integrator

A VCO's frequency obeys $\nu = \nu_0 + K_0 V(t)$. Show that its voltage-to-phase transfer function is $2\pi K_0 / s$, and hence that even a flat loop filter $F(f) = F_0$ gives a first-order low-pass $H$. What is the loop bandwidth?

Solution

Chapter 1: frequency is the slope of phase, so $\dot\varphi = 2\pi\,\delta\nu = 2\pi K_0 V(t)$, and integrating (in Laplace language, dividing by $s$) gives $\varphi(s) = 2\pi K_0 V(s)/s$ — an integrator, with no capacitor in sight. The open-loop gain is then $G(s) = K_d F_0 \cdot 2\pi K_0/s \equiv K/s$, and $H = G/(1+G) = 1/(1 + s/K)$: a one-pole low-pass with 3 dB bandwidth $f_L = K/2\pi = K_d K_0 F_0$ in Hz. The $1/f$ that makes the loop a low-pass is supplied by the physics of the VCO itself.

Exercise 10.3 — splitting the line

For a Lorentzian resonance of FWHM $\Delta\nu$ and peak signal $S_0$, the discriminator slope is at best of order $S_0/\Delta\nu$ (check this against the S-curve demo). With $\mathrm{SNR} = 1000$ in 1 s of averaging, what frequency uncertainty do you reach in 1 s? Convert to $\sigma_y(1\,\mathrm{s})$ for a caesium beam tube with $\Delta\nu \approx 100$ Hz.

Solution

An error-signal noise of $S_0/\mathrm{SNR}$ divided by a slope $\sim S_0/\Delta\nu$ gives $\delta\nu \approx \Delta\nu/\mathrm{SNR} = 100\,\mathrm{Hz}/1000 = 0.1$ Hz — you split the line a thousand ways. Fractionally, $\sigma_y(1\,\mathrm{s}) \approx 0.1 / 9.19\times10^9 \approx 1\times10^{-11}$, in the right ballpark for a classic beam standard. Longer averaging improves it as $\tau^{-1/2}$ (chapter 7), until systematic shifts of the line itself take over.

Exercise 10.4 — accurate but not stable

Why does an integrator in the loop filter make an atomic clock accurate (zero mean frequency error, inherited from the atoms) while doing nothing to remove its short-term instability?

Solution

An integrator has unbounded gain at DC: any nonzero average error signal keeps charging it, sliding the quartz until the average is exactly zero. So the mean output frequency is pinned to the atomic line regardless of the quartz's calibration or drift — that is accuracy. But a loop only reshapes noise (chapter 9): within the bandwidth the output carries the discriminator's noise — atomic-detection SNR, white frequency noise, the $\tau^{-1/2}$ of exercise 10.3 — and above the bandwidth the quartz free-runs. Zero mean error says nothing about the variance around that mean; the integral of $S(f)\,|W(f)|^2$ is moved around in frequency, not destroyed.

Exercise 10.5 — the bump on the beat note

Your cavity-locked laser's beat note against an independent reference shows a noise bump at 1 MHz that was absent unlocked. Which loop parameter do you reach for, and what do you not do?

Solution

The bump marks the loop's unity-gain frequency, where $|1+G|$ dips below 1 because the phase margin is small — chapter 9's servo bump, here in a laser lock. Recover phase near 1 MHz: shorten delay, remove a low-pass or resonance that eats phase there, or add a lead (differentiator) stage; alternatively lower the gain so unity gain retreats to a frequency with more margin. What you do not do is crank up the gain — that pushes $G$ closer to $-1$ at the bump and makes it taller.

Where you started, where you are

Ten chapters ago this course opened with one equation and three knobs: $V(t) = A\cos(2\pi\nu_0 t + \varphi)$. Everything since has been that equation under stress. The knobs wiggled deliberately gave sidebands (chapter 2); wiggled randomly they needed ensembles (chapter 3) and spectral densities (chapter 4); the spectra turned out to be power laws with personalities (chapter 5), which on the phase knob became $\mathcal{L}(f)$, linewidths and jitter (chapter 6), compressed into $\sigma_y(\tau)$ for clocks (chapter 7), priced in per-root-hertz for sensors (chapter 8), stacked into budgets and filter functions for LIGO (chapter 9), and finally reshaped at will by a feedback ring (this chapter). If one sentence survives: noise is a spectral density; every measurement is a filter function; every quoted number is their integral. The next time an Allan plot, a dBc/Hz spec or a lock-loop block diagram crosses your desk, it should read like prose. Go twiddle the real knobs.