Skip to main content

From Impulse to Phase Shift — the Derivation

β: This English translation is in beta — the Traditional-Chinese original is the authoritative version.

Prerequisites: oscillator_phase (limit-cycle and phase/amplitude geometry), lti_vs_ltv (why an oscillator is LTV with respect to noise) | Up next: isf_definition (the rigorous ISF definition and the multi-node picture) → convolution_derivation (superposing arbitrary noise by convolution)

This page answers a very concrete question: when a current impulse is injected into some node of an oscillator, by how much does the oscillator's phase shift? The answer is the operational definition at the heart of ISF theory:

Δϕ=Γ(ω0τ)qmaxΔq\Delta\phi=\frac{\Gamma(\omega_0\tau)}{q_{max}}\,\Delta q

We will not simply memorize this line — we will derive it from the capacitor relation q=Cvq=Cv, one step at a time, stating at every step which physics is invoked, which approximation is made, what the units are, and how to run the dimension check.

Physical intuition (conclusion first): the noise current first becomes a packet of charge Δq\Delta q on the node; that charge nudges the oscillator's instantaneous state slightly off its steady orbit; how much of that nudge turns into phase and how much into amplitude depends on at which phase of the waveform you deliver the kick — and that conversion ratio is precisely the ISF Γ\Gamma. Kick where the waveform slope is large and the phase gets pushed a lot; kick at the peak and you mostly perturb the amplitude.

Step 1: the current impulse becomes charge

Current integrated over time is charge. A very narrow, very short current pulse i(t)i(t) injected into the node deposits a total charge

Δq=i(t)dt.\Delta q=\int i(t)\,dt .
  • Physics used: charge conservation / the definition of current, i=dq/dti=dq/dt.
  • Unit check: [A][s]=[C][\text{A}]\cdot[\text{s}]=[\text{C}] ✓ (amperes times seconds are coulombs).
  • Why this is legitimate: when the pulse width is far smaller than the period TT, the injection can be treated as dumping Δq\Delta q onto the node "instantaneously".

Step 2: charge becomes a voltage step

The node carries a total capacitance CnodeC_{node}. Charge appearing instantaneously makes the node voltage jump by a step ([P1] Eq.(9), p.181):

ΔV=ΔqCnode.\Delta V=\frac{\Delta q}{C_{node}} .
  • Physics used: the capacitor relation q=Cvq=Cv; differentiate for a small change, Δq=CΔV\Delta q=C\,\Delta V.
  • Unit check: [C]/[F]=[C]/[C/V]=[V][\text{C}]/[\text{F}]=[\text{C}]/[\text{C/V}]=[\text{V}] ✓.
  • Key observation: in a parallel LC, a current impulse changes only the capacitor voltage (instantaneously); it cannot change the inductor current (an inductor current cannot jump). In state space, the perturbation is therefore a horizontal displacement along the voltage axis.

Step 3: the voltage step pushes the state off the limit cycle — but only the tangential component becomes phase

Draw the oscillator's state in a 2-D plane (say, horizontal axis = capacitor voltage vv, vertical axis = ww, proportional to the inductor current). In steady state, the state point circulates along the limit cycle. The voltage jump from Step 2 shoves the state point sideways:

  • The radial component, pushing away from the cycle → changes the amplitude AA. But the oscillator has an amplitude-restoring mechanism, so this part is slowly pulled back and leaves no permanent residue.
  • The tangential component, along the cycle → changes the phase ϕ\phi. Phase has no restoring force, so this part stays forever.

The same ΔV\Delta V, kicked at a different phase τ\tau, splits differently between tangential and radial. Collect the "ΔV\Delta V-to-Δϕ\Delta\phi conversion ratio" into a periodic function that depends only on the injection phase — that is the ISF:

Δϕ=(phase-conversion ratio)depends only on ω0τ×ΔV.\Delta\phi=\underbrace{\big(\text{phase-conversion ratio}\big)}_{\text{depends only on }\omega_0\tau}\times\Delta V .

The figure below draws this decomposition on the state plane. The same packet of charge produces a horizontal ΔV=Δq/C\Delta V=\Delta q/C (the current kicks only the capacitor-voltage axis; the inductor current cannot jump), which the dashed lines split into a tangential component along the cycle (green — permanent phase) and a radial component off the cycle (red — decays). Injected near the zero crossing, the horizontal kick lands almost entirely on the tangent (the phase gets pushed a lot); injected near the peak, it lands almost entirely on the radius (only the amplitude changes). The lower-left corner of each panel prints the value of Γ(θ)=sinθ\Gamma(\theta)=-\sin\thetathe tangential fraction is exactly Γ(θ)\Gamma(\theta), which makes the ideal-LC Γ=sinθ\Gamma=-\sin\theta geometrically inevitable (the algebraic derivation of Γ=sinθ\Gamma=-\sin\theta is in isf_definition).

Impulse decomposition on the ideal-LC state plane: the same horizontal ΔV=Δq/C is split by the dashed lines into a tangential component (→Δφ, permanent phase) and a radial component (→ΔA, decays); near the zero crossing the kick is almost entirely tangential, near the peak almost entirely radial, and the tangential fraction is exactly Γ(θ)=−sin θ. Pedagogical geometry for an ideal lossless LC, not transistor-level.

  • Math used: project the perturbation vector onto the tangent direction (the phase direction) of the limit cycle.
  • Why only the tangential part is kept: see phase_vs_amplitude_noise and oscillator_phase — amplitude is pulled back by the restoring mechanism; phase is not.

Step 4: normalize by qmaxq_{max} to get the dimensionless ISF

Chain Steps 1–3 together: ΔϕΔV=Δq/Cnode\Delta\phi\propto\Delta V=\Delta q/C_{node}. Hajimiri–Lee chose to normalize by the node's maximum charge swing qmax=CnodeVmaxq_{max}=C_{node}V_{max}, writing that "phase-conversion ratio" as the dimensionless function Γ(ω0τ)\Gamma(\omega_0\tau) ([P1] Eq.(10)–(11), p.182):

 Δϕ=Γ(ω0τ)qmaxΔq \boxed{\ \Delta\phi=\frac{\Gamma(\omega_0\tau)}{q_{max}}\,\Delta q\ }
  • Why the qmaxq_{max} normalization: it turns Γ\Gamma into a dimensionless, amplitude-independent "shape" that describes only where the waveform is sensitive. The actual magnitude of the phase shift is set by Δq/qmax\Delta q/q_{max} (the injected charge relative to the signal charge).
  • Why Γ\Gamma is dimensionless: Δϕ\Delta\phi is in rad (dimensionless) and Δq/qmax\Delta q/q_{max} is dimensionless, so Γ\Gamma must be dimensionless. Dimension check: [rad]=Γ[C]/[C][\text{rad}]=\Gamma\cdot[\text{C}]/[\text{C}] Γ\Rightarrow\Gamma dimensionless ✓.
  • Small-signal assumption: Step 3 treats the projection as linear, which requires Δqqmax\Delta q\ll q_{max} (the kick must not knock the oscillator over). [P1] Fig. 6 confirms, on an actual Colpitts and a 5-stage ring, that Δϕ\Delta\phi is proportional to Δq\Delta q for small charge.

The corresponding impulse response (setting up the next chapter)

Because the phase step is retained forever, writing it as an impulse response brings along a unit step u(tτ)u(t-\tau) ([P1] Eq.(10)):

hϕ(t,τ)=Γ(ω0τ)qmaxu(tτ).h_\phi(t,\tau)=\frac{\Gamma(\omega_0\tau)}{q_{max}}\,u(t-\tau).

Note that it depends on the absolute injection instant τ\tau (through Γ(ω0τ)\Gamma(\omega_0\tau)), not merely on tτt-\tau — the very signature of an LTV (linear time-variant) system. The next chapter, convolution_derivation, uses it to superpose an arbitrary noise current. The complete ISF definition and the multi-node discussion are in isf_definition.

Numerical example (building a feel for the numbers)

Example A: qmax=1q_{max}=1 pC, Δq=1\Delta q=1 fC, Γ=0.5\Gamma=0.5, f0=5f_0=5 GHz.

Phase step:

Δϕ=ΓΔqqmax=0.5×(1×1015C)1×1012C=5×104 rad.\Delta\phi=\frac{\Gamma\,\Delta q}{q_{max}}=\frac{0.5\times(1\times10^{-15}\,\text{C})}{1\times10^{-12}\,\text{C}}=5\times10^{-4}\ \text{rad}.

In degrees: Δϕ=5×104×180π0.0286\Delta\phi=5\times10^{-4}\times\dfrac{180}{\pi}\approx0.0286^\circ.

Converted to a timing error (using Δt=Δϕ/(2πf0)\Delta t=\Delta\phi/(2\pi f_0)):

Δt=5×104 rad2π×5×109 Hz=5×1043.1416×1010 s1.59×1014 s=15.9 fs.\Delta t=\frac{5\times10^{-4}\ \text{rad}}{2\pi\times5\times10^{9}\ \text{Hz}}=\frac{5\times10^{-4}}{3.1416\times10^{10}}\ \text{s}\approx1.59\times10^{-14}\ \text{s}=15.9\ \text{fs}.
  • Dimension check: [rad]/[rad/s]=[s][\text{rad}]/[\text{rad/s}]=[\text{s}] ✓ (note that 2πf02\pi f_0 carries units of rad/s).
  • Feel for the numbers: at 5 GHz (a 200 ps period), a single 1 fC packet of charge (about 6240 electrons) causes only ~16 fs of timing error even at the most sensitive phase. One kick is tiny — but the noise kicks continuously, and the integral accumulates (next chapter).

Quick verification with the built-in functions:

from simulations.common.isf_utils import impulse_to_phase_step
from simulations.common.noise_utils import phase_to_time_error

dphi = impulse_to_phase_step(delta_q=1e-15, gamma_value=0.5, qmax=1e-12)
dt = phase_to_time_error(dphi, f0=5e9)
print(dphi, "rad", dt*1e15, "fs") # -> 0.0005 rad 15.92 fs

(Full scripts: simulations/common/isf_utils.py, simulations/common/noise_utils.py.)

Seeing is believing: measuring Γ\Gamma directly

lab_04 injects a small charge at a sweep of phases in simulation, measures the persistent phase offset, and back-solves for the ISF; the result agrees almost perfectly with the ideal-LC Γ(θ)=sinθ\Gamma(\theta)=-\sin\theta (maximum error about 0.001):

Numerically extracted ISF versus the theoretical -sin(θ)

Where it applies, where it fails

ConditionWhen it holdsWhat happens when it fails
Small signal Δqqmax\Delta q\ll q_{max}Δϕ\Delta\phi is linearly proportional to Δq\Delta qlarge injection → nonlinearity, AM–PM, and the ISF itself is altered
Amplitude perturbations decayonly the phase needs trackingbreaks down with strong AM–PM or without a stable limit cycle
Pulse far narrower than the periodcan be treated as an instantaneous injectionwide pulses require the integral form, Eq.(11)
The correct Γ\Gamma is knownpredictions are accurateΓ\Gamma must be extracted by transient/adjoint simulation (see effective_isf)

Key takeaways

  • Noise current → charge Δq\Delta q → voltage jump ΔV=Δq/C\Delta V=\Delta q/C → projected through the ISF into phase Δϕ\Delta\phi.
  • Δϕ=Γ(ω0τ)Δq/qmax\Delta\phi=\Gamma(\omega_0\tau)\,\Delta q/q_{max}; Γ\Gamma is dimensionless, 2π2\pi-periodic, and depends on the injection phase.
  • qmaxq_{max} normalizes Γ\Gamma into a "shape"; the magnitude of the phase shift is set by Δq/qmax\Delta q/q_{max}.
  • A single 1 fC at 5 GHz, with Γ=0.5\Gamma=0.5 and qmax=1q_{max}=1 pC → 16 fs.
  • Sources: [P1] Eq.(9) p.181, Eqs.(10),(11) p.182; verification figure in lab_04.

Further reading