Dynamical Decoupling

A qubit left alone does not stay still. Its phase drifts because the surrounding material — nuclear spins, trapped charges, flux noise, stray fields — pushes on it with a signal that wanders slowly in time. Dynamical decoupling is the art of turning that wandering against itself: apply a fast train of strong control pulses that periodically flip the qubit, so that whatever the environment did during one interval is undone during the next. The accumulated coupling integrates to (almost) zero, and the coherence that would otherwise leak away survives. The oldest instance is the nuclear-magnetic-resonance spin echo of Erwin Hahn (1950): one π pulse halfway through a free evolution refocuses any static frequency offset, resurrecting a signal that inhomogeneous broadening had scrambled. Carr, Purcell, Meiboom and Gill turned the single echo into a robust multi-pulse train; four decades later Viola, Knill and Lloyd recast the whole idea in the language of open quantum systems, and Uhrig found the optimal pulse timings for a sharply cut-off bath. The physics rests on one honest assumption — a separation of timescales. The bath must vary slowly compared with the spacing between pulses, so a fast cycle can catch it before it changes. When that holds, dynamical decoupling is spectacularly effective and, unlike quantum error correction, needs no measurement, no ancillas, and no feedback — only well-timed unitaries. When it fails — against fast, memoryless, white noise — no pulse train can help, because there is no temporal correlation left to exploit. This article makes both the promise and the limit precise.

Key takeaways
  • Dynamical decoupling (DD) suppresses decoherence with a train of fast, strong control pulses that flip the qubit so its coupling to a slowly fluctuating environment averages toward zero over each cycle.
  • The founding case is the Hahn spin echo: one π pulse at the midpoint cancels any static dephasing exactly, recovering coherence lost to T2∗T_2^* inhomogeneous broadening.
  • In the toggling frame set by the pulses, the system–bath coupling becomes time-dependent; a good sequence makes its cycle average — the zeroth-order Magnus term Hˉ(0)\bar H^{(0)} — vanish, pushing loss to order Tc2T_c^2 and beyond.
  • Equivalently, DD is a high-pass filter: the pulse train imprints a filter function F(ωT)F(\omega T) with a deep notch at zero frequency, rejecting the slow noise that dominates solid-state qubits.
  • DD is open-loop and coherent — no syndrome measurement, no ancillas — which makes it complementary to, not a substitute for, active error correction.
  • It cannot beat white (Markovian) noise and cannot lift T2T_2 past the relaxation ceiling 2T12T_1; Uhrig DD reaches NN-th order suppression of a hard-cutoff bath with only NN pulses.
  • Key idea: in the toggling frame defined by the control pulses the system–bath coupling becomes a time-dependent operator whose cycle-averaged value can be engineered to vanish; driving that average Hamiltonian to zero, and its Magnus corrections to high order, freezes decoherence over each cycle.
The core insight

In the toggling frame defined by the control pulses, the system–bath coupling becomes a time-dependent operator whose cycle-averaged value can be engineered to vanish; pushing this average Hamiltonian to zero (and its Magnus corrections to high order) freezes decoherence over each cycle.

Where this lives today

Dynamical decoupling has moved from an NMR curiosity to standard infrastructure on nearly every quantum processor. On superconducting transmons and semiconductor spin qubits, CPMG and XY-8 trains routinely convert an inhomogeneous T2∗T_2^* of a few microseconds into a decoupled T2T_2 approaching the 2T12T_1 ceiling, because the dominant 1/f1/f charge and flux noise lives exactly in the low-frequency band DD notches out. The same filter that protects also measures: sweeping the inter-pulse spacing turns a DD sequence into a tunable spectrometer, and "DD noise spectroscopy" is now the workhorse method for reconstructing S(ω)S(\omega) and diagnosing the defects — two-level systems, nuclear-spin baths — that limit a device. In nitrogen-vacancy centers in diamond it does double duty: extending the sensor spin's coherence into the millisecond range and, by tuning the pulse spacing into resonance with a target nuclear precession, enabling single-nuclear-spin NMR and nanoscale magnetometry. On gate-model processors DD is compiled into idle windows to protect spectator qubits, and it composes with quantum error correction rather than competing with it — suppressing the correlated, non-Markovian errors codes handle poorly and helping push physical fidelities below the fault-tolerance threshold. Open problems remain live: optimizing sequences against realistic pulse errors, decoupling many-qubit crosstalk without disturbing entangling gates, and "dynamically corrected gates" that decouple during, not just between, operations.

The physical picture

Picture a crowd of runners set loose together on a circular track, each carrying a slightly different, unknown pace because the wind pushes each one a little differently. Left alone they fan out — the pack blurs into a smear, and the coherent "signal" of a tight group is lost. Now, at the halfway whistle, spin every runner around to face backward. Each retraces their lead at the same pace they built it, and at the final whistle the entire pack crosses the start line together again. The individual paces never had to be known or equal; they only had to hold steady across the flip. That is the spin echo, and the first figure is exactly this refocusing drawn as accumulated phase against time: divergent lines during dephasing, a π pulse that mirrors them, convergence to zero at the echo.

Figure
Hahn spin echo as phase refocusing: spins with different detunings δ fan out during the first interval, the π_x pulse at time τ reflects every accumulated phase through zero, and by t = 2τ they reconverge to a single echo — the net static phase is φ_net = (δ+b)τ − (δ+b)τ = 0 for any offset
The spin echo: a π pulse reflects every phase through zero, so any static offset refocuses at 2τ accumulated phase φ t = 0 t = τ t = 2τ π_x pulse dephasing rephasing echo: φ_net = 0 φ → −φ

Dynamical decoupling generalizes the single echo in two complementary ways, and both are worth carrying in your head. The first is geometric: instead of one flip, chop the evolution into many short intervals and rotate the qubit between them so that the bath operator it couples to is presented in a set of orientations that sum to zero — averaging a vector over the vertices of a symmetric figure. Because the bath barely moves during one fast cycle, this cancellation survives and the residual error scales with a high power of the small cycle time. The second view is spectral, and it is the one an experimentalist actually uses. Free evolution lets the qubit accumulate phase from every noise frequency, weighted heavily toward DC; a pulse train modulates that response into a filter function with a deep notch at zero frequency and a passband near the pulse-repetition rate. Slow noise falls in the notch and is rejected. The second figure is the punchline: how sharply that notch cuts, and where the passband sits, for real sequences.

Figure
Filter-function view on log–log axes: DD is a high-pass filter with F(0)=0. A spin-echo or CPMG filter rises as (ωT)² out of the DC notch, while Uhrig DD with N=5 pulses rises as (ωT)^(2N+2) = (ωT)^12, a far deeper notch. The passband sits at ω ≈ Nπ/T, e.g. N=8 pulses over T=80 μs give f ≈ N/2T = 50 kHz
The punchline: DD is a high-pass filter — deeper notch, higher order, more coherence noise frequency ω T (log scale) filter |F(ωT)|² (log) 1 10⁻⁴ passband ω* ≈ Nπ/T CPMG / echo: |F|² ∝ (ωT)² UDD N=5: |F|² ∝ (ωT)¹² noise S(ω) ∝ 1/ω DC notch: F(0) = 0 N=8, T=80 μs ⇒ f* ≈ N/2T = 50 kHz

The core equations

Everything begins with the joint description of qubit plus environment. The total Hamiltonian splits into free parts and a coupling, and the coupling is written as a sum of products — a system operator times a bath operator — because that factored form is what the toggling frame will act on.

H=HS+HB+HSB,HSB=∑αSα⊗Bα(1)H = H_S + H_B + H_{SB}, \qquad H_{SB} = \sum_\alpha S_\alpha \otimes B_\alpha \qquad (1)

where: HS, HBH_S,\,H_B — free Hamiltonians of the qubit and the bath; HSBH_{SB} — their coupling; SαS_\alpha — Hermitian qubit operators (e.g. Pauli σz\sigma_z for pure dephasing); BαB_\alpha — bath operators, assumed to fluctuate slowly. The tensor product ⊗\otimes separates the two Hilbert spaces the pulses will treat differently: the pulses act on SαS_\alpha, never on BαB_\alpha. The cleanest and most important special case is pure dephasing, where the coupling commutes with the qubit's energy operator, so populations are untouched and only the relative phase between ∣0⟩|0\rangle and ∣1⟩|1\rangle decays. This is the noise DD is built to defeat.

HSBdeph=σz⊗B,HSdeph=12 δ σz(2)H_{SB}^{\text{deph}} = \sigma_z \otimes B, \qquad H_S^{\text{deph}} = \tfrac{1}{2}\,\delta\,\sigma_z \qquad (2)

where: σz\sigma_z — Pauli-zz operator, which commutes with HSBH_{SB} so no energy is exchanged; BB — the bath (or classical noise) operator multiplying σz\sigma_z; δ\delta — a static frequency detuning, the inhomogeneous offset responsible for the fast T2∗T_2^* decay a single echo removes. To see what the pulses do, move into the frame that rotates with the control alone — the toggling frame. Conjugating the coupling by the control propagator turns a static interaction into a time-dependent one whose time dependence is entirely under the experimenter's control.

H~SB(t)=Uc†(t) HSB Uc(t),Uc(t)=Texp⁡ ⁣(−i ⁣∫0tHc(s) ds)(3)\tilde H_{SB}(t) = U_c^\dagger(t)\, H_{SB}\, U_c(t), \qquad U_c(t)=\mathcal{T}\exp\!\left(-i\!\int_0^t H_c(s)\,\mathrm{d}s\right) \qquad (3)

where: Uc(t)U_c(t) — unitary generated by the control Hamiltonian Hc(t)H_c(t) alone; H~SB(t)\tilde H_{SB}(t) — the coupling seen in the frame co-rotating with the pulses; T\mathcal{T} — the time-ordering operator, needed because HcH_c at different times need not commute. Over one full pulse cycle the toggling-frame evolution can be summarized by a single effective generator, the average Hamiltonian, whose leading term is just the time average of the toggled coupling. Making this term vanish is the entire design goal of first-order DD.

Hˉ(0)=1Tc∫0TcH~SB(t) dt(4)\bar H^{(0)} = \frac{1}{T_c}\int_0^{T_c} \tilde H_{SB}(t)\,\mathrm{d}t \qquad (4)

where: TcT_c — the duration of one decoupling cycle; Hˉ(0)\bar H^{(0)} — the zeroth-order term of the Magnus expansion of the cycle propagator. When Hˉ(0)=0\bar H^{(0)}=0 the qubit feels no environment to leading order, and residual decoherence is demoted to O(Tc2)\mathcal{O}(T_c^2) through the higher Magnus terms. Why any traceless coupling can be killed is a one-line fact of group theory: averaging a traceless operator over conjugation by the single-qubit Pauli group projects it onto the identity, whose traceless part is zero. This is the algebraic engine behind universal decoupling and the XY-4 cycle.

14∑g∈{I,X,Y,Z}g† Sα g=12 Tr(Sα) 1=0(traceless Sα)(5)\frac{1}{4}\sum_{g\in\{I,X,Y,Z\}} g^\dagger\, S_\alpha\, g = \tfrac{1}{2}\,\mathrm{Tr}(S_\alpha)\,\mathbb{1} = 0 \quad(\text{traceless }S_\alpha) \qquad (5)

where: gg — elements of the decoupling group (here the four single-qubit Paulis); SαS_\alpha — any traceless qubit coupling operator; 1\mathbb{1} — the identity. Because XZX=−ZXZX=-Z and YZY=−ZYZY=-Z while IZI=ZZZ=ZIZI=ZZZ=Z, the four conjugates of σz\sigma_z are {Z,−Z,−Z,Z}\{Z,-Z,-Z,Z\} and sum to zero. The spectral face of the same physics is the filter function, the tool used to predict and measure coherence. For Gaussian dephasing the coherence decays as e−W(T)e^{-W(T)}, with WW an overlap integral of the noise power spectrum against the sequence's filter.

W(T)=1π∫0∞dω  S(ω) ∣F(ωT)∣2ω2,∣ρ01(T)∣=∣ρ01(0)∣ e−W(T)(6)W(T) = \frac{1}{\pi}\int_0^\infty \mathrm{d}\omega\; S(\omega)\, \frac{|F(\omega T)|^2}{\omega^2}, \qquad |\rho_{01}(T)| = |\rho_{01}(0)|\,e^{-W(T)} \qquad (6)

where: W(T)W(T) — the decay exponent; S(ω)S(\omega) — the one-sided power spectral density of the dephasing field; F(ωT)F(\omega T) — the filter function fixed by the pulse timings; ρ01\rho_{01} — the off-diagonal coherence. The factor 1/ω21/\omega^2 makes free evolution most sensitive to low frequencies; the filter's job is to suppress FF there. Uhrig's insight was that the pulse times themselves are free parameters to be optimized. For a bath with a sharp high-frequency cutoff, a specific non-uniform placement makes the filter vanish to the highest possible order at low frequency.

tj=T sin⁡2 ⁣(jπ2N+2),j=1,…,N(7)t_j = T\,\sin^2\!\left(\frac{j\pi}{2N+2}\right), \quad j=1,\dots,N \qquad (7)

where: tjt_j — the time of the jj-th π pulse within total evolution time TT; NN — the number of pulses. The timings are symmetric about T/2T/2 and reproduce the single echo at N=1N=1 (then t1=T/2t_1=T/2). The payoff of that choice is a filter that is flat to order 2N+22N+2 at the origin — each pulse buys two extra powers of suppression, the sharpest low-frequency rejection any NN-pulse sequence can achieve for a hard cutoff.

∣FN(ωT)∣2∼(ωT)2N+2(ωT→0)(8)|F_N(\omega T)|^2 \sim (\omega T)^{2N+2} \qquad (\omega T \to 0) \qquad (8)

where: FNF_N — the UDD filter with NN pulses; the exponent 2N+22N+2 is the order to which slow noise is annihilated. Uniform CPMG rises only as (ωT)2(\omega T)^2 regardless of NN, so for a truly hard-cutoff bath UDD wins decisively. Finally, the ceiling. DD attacks dephasing, but it cannot touch energy relaxation, which changes populations and sets a hard bound on how long any coherence can survive.

1T2=12T1+1Tϕ    ⟹    T2≤2 T1(9)\frac{1}{T_2} = \frac{1}{2T_1} + \frac{1}{T_\phi} \;\;\Longrightarrow\;\; T_2 \le 2\,T_1 \qquad (9)

where: T2T_2 — the transverse (coherence) time DD extends; T1T_1 — the longitudinal (relaxation) time; TϕT_\phi — the pure-dephasing time DD can in principle send to infinity. Even perfect decoupling (Tϕ→∞T_\phi\to\infty) leaves T2=2T1T_2=2T_1; no open-loop single-qubit sequence beats it.

Derivation ASpin echo from average Hamiltonian theory

The goal is to show, with no approximation beyond a static bath, that the single-π-pulse cycle drives the zeroth-order average Hamiltonian of pure dephasing exactly to zero — and that the first Magnus correction vanishes too, so the echo is perfect for a frozen environment and its residual comes only from how much the bath moved during the cycle. We take HSB=σz⊗BH_{SB}=\sigma_z\otimes B with BB constant over one cycle and the control sequence "free-evolve τ\tau, apply πx\pi_x, free-evolve τ\tau, apply πx\pi_x" that returns control to the identity, then compute the toggling-frame coupling piece by piece.

Step 1. In the first interval the control has done nothing, so H~SB(t)=σz⊗B\tilde H_{SB}(t)=\sigma_z\otimes B for 0<t<τ0<t<\tau. After the first πx\pi_x the toggling frame conjugates the coupling by σx\sigma_x, and using σxσzσx=−σz\sigma_x\sigma_z\sigma_x=-\sigma_z we get H~SB(t)=−σz⊗B\tilde H_{SB}(t)=-\sigma_z\otimes B for τ<t<2τ\tau<t<2\tau. Why: the pulse does not change the bath, it relabels the qubit axis, and flipping the sign of σz\sigma_z is precisely how a π rotation about xx reverses phase accumulation.

Step 2. Average over the cycle: Hˉ(0)=12τ∫02τH~SB dt=12τ[(σz⊗B)τ+(−σz⊗B)τ]=12(σz−σz)⊗B=0\bar H^{(0)}=\frac{1}{2\tau}\int_0^{2\tau}\tilde H_{SB}\,\mathrm{d}t = \frac{1}{2\tau}\big[(\sigma_z\otimes B)\tau + (-\sigma_z\otimes B)\tau\big] = \tfrac{1}{2}(\sigma_z-\sigma_z)\otimes B = 0. Why: the two halves carry equal and opposite coupling for exactly equal times, so the accumulated phase kick from the first half is undone by the second — this is the echo, stated as an average Hamiltonian.

Step 3. Check the next order. The first Magnus correction is built from the commutator [H~SB(t1),H~SB(t2)][\tilde H_{SB}(t_1),\tilde H_{SB}(t_2)], but both toggling-frame couplings are proportional to σz⊗B\sigma_z\otimes B and therefore commute, so Hˉ(1)=0\bar H^{(1)}=0 as well. Why: with a single dephasing axis nothing fails to commute, so the only surviving error is the higher-order effect of BB actually changing between the two halves — exactly the slow-noise leakage the filter function will quantify. *As a check: the classical static result is ϕnet=(δ+b)τ−(δ+b)τ=0\phi_{\text{net}}=(\delta+b)\tau-(\delta+b)\tau=0, confirming that any static offset — known or unknown — is cancelled exactly.*

Go deeper · Derivation B · 2 minThe CPMG filter function and the rejection of slow noise
Derivation BThe CPMG filter function and the rejection of slow noise

The goal is to derive, from a classical model of a fluctuating detuning, the frequency-domain formula behind every dephasing DD experiment — the coherence decay exponent as an overlap of the noise spectrum with a sequence-specific filter — and to show that any symmetric pulse train forces that filter to zero at zero frequency, so static and slow noise are rejected exactly. We model the dephasing as H(t)=12 η(t) σzH(t)=\tfrac12\,\eta(t)\,\sigma_z with a zero-mean, stationary, Gaussian field η(t)\eta(t), track the phase accumulated through a train of sign flips, and move to the frequency domain.

Step 1. The accumulated phase is ϕ(T)=∫0Ty(t) η(t) dt\phi(T)=\int_0^T y(t)\,\eta(t)\,\mathrm{d}t, where y(t)=±1y(t)=\pm1 is the sign of σz\sigma_z evolution, starting at +1+1 and flipping at each π-pulse time tjt_j. Why: each π pulse reverses which computational state runs "fast", so the effective coupling to η\eta is a square wave — the pulse sequence is literally encoded in the switching function y(t)y(t).

Step 2. For Gaussian noise the coherence is ⟨e−iϕ⟩=e−12⟨ϕ2⟩≡e−W(T)\langle e^{-i\phi}\rangle = e^{-\frac12\langle\phi^2\rangle}\equiv e^{-W(T)}, so all we need is the phase variance. Writing ⟨η(t)η(t′)⟩\langle\eta(t)\eta(t')\rangle in terms of its spectrum S(ω)S(\omega) and doing the two time integrals gives ⟨ϕ2⟩=2π∫0∞dω S(ω) ∣Y(ω,T)∣2\langle\phi^2\rangle=\frac{2}{\pi}\int_0^\infty \mathrm{d}\omega\, S(\omega)\,|Y(\omega,T)|^2, where YY is the Fourier transform of y(t)y(t). Why: the decay is set entirely by the second moment of a Gaussian, and correlations turn into a spectral density, converting a messy time-domain average into a clean product in frequency.

Step 3. Collect the transform into the standard filter F(ωT)=1+(−1)N+1eiωT+2∑j=1N(−1)jeiωtjF(\omega T)=1+(-1)^{N+1}e^{i\omega T}+2\sum_{j=1}^{N}(-1)^j e^{i\omega t_j}, so that W(T)=1π∫0∞dω S(ω) ∣F(ωT)∣2/ω2W(T)=\frac{1}{\pi}\int_0^\infty \mathrm{d}\omega\, S(\omega)\,|F(\omega T)|^2/\omega^2. Why: the endpoints and each pulse contribute one term; the alternating signs are the sign flips of y(t)y(t), and the 1/ω21/\omega^2 is the transform of integrating η\eta into a phase.

Step 4. Evaluate at ω→0\omega\to0. A symmetric train has equal total time spent at y=+1y=+1 and y=−1y=-1, so ∫0Ty(t) dt=0\int_0^T y(t)\,\mathrm{d}t=0, which is exactly F(0)=0F(0)=0. Why: F(0)F(0) is the DC response, i.e. the net area under the switching function; balancing the plus and minus segments annihilates a perfectly static offset, recovering the spin echo and explaining why slow noise is suppressed while noise near the passband leaks through. *As a check: setting N=1, t1=T/2N=1,\ t_1=T/2 gives F=1+eiωT−2eiωT/2=−4eiωT/2sin⁡2(ωT/4)F=1+e^{i\omega T}-2e^{i\omega T/2}=-4e^{i\omega T/2}\sin^2(\omega T/4), so ∣F∣2=16sin⁡4(ωT/4)|F|^2=16\sin^4(\omega T/4) with F(0)=0F(0)=0 — the Hahn-echo filter.*

Go deeper · Derivation C · 1 minPauli-group averaging and universal single-qubit decoupling
Derivation CPauli-group averaging and universal single-qubit decoupling

The goal is to prove that a four-pulse cycle can decouple a qubit from an arbitrary single-qubit coupling — not just the dephasing σz\sigma_z but any Hermitian traceless SαS_\alpha — by showing that averaging over conjugation by the Pauli group projects every traceless operator to zero. This "universal" claim is what lifts DD from a dephasing trick to a general Hamiltonian-engineering tool, and it is what the XY-4 sequence realizes in hardware. We use the Pauli multiplication rules and trace-orthogonality.

Step 1. Write any single-qubit coupling operator as S=c01+cxσx+cyσy+czσzS=c_0\mathbb{1}+c_x\sigma_x+c_y\sigma_y+c_z\sigma_z; Hermiticity makes the cc's real and, being a coupling, SS is traceless so c0=0c_0=0. Why: the four Paulis are a basis for 2×22\times2 Hermitian operators, so decoupling all of them decouples everything — no case is left out.

Step 2. Conjugate by each element of {I,X,Y,Z}\{I,X,Y,Z\} and sum. Each Pauli anticommutes with two others and commutes with itself and the identity, so for σz\sigma_z the four conjugates are IZI=ZIZI=Z, XZX=−ZXZX=-Z, YZY=−ZYZY=-Z, ZZZ=ZZZZ=Z, giving {Z,−Z,−Z,Z}\{Z,-Z,-Z,Z\} whose sum is zero; the same balance holds for σx\sigma_x and σy\sigma_y. Why: every traceless Pauli is anticommuted (sign-flipped) by exactly two group elements and left alone by two, so the group average cancels it term by term.

Step 3. Hence 14∑gg†S g=12Tr(S) 1=0\frac14\sum_g g^\dagger S\,g=\frac12\mathrm{Tr}(S)\,\mathbb{1}=0 for traceless SS. A cycle whose toggling frames visit the four Pauli axes therefore sets Hˉ(0)=0\bar H^{(0)}=0 for any single-qubit coupling — this is XY-4. Why: realizing the four group elements as successive toggling frames turns the abstract group average into the physical cycle average of equation (4), so the algebra and the pulse sequence are the same statement. *As a check: the conjugation signs XZX=YZY=−ZXZX=YZY=-Z and IZI=ZZZ=ZIZI=ZZZ=Z sum to zero, and the projection 14∑gg†Sg=12Tr(S)1\frac14\sum_g g^\dagger S g=\frac12\mathrm{Tr}(S)\mathbb{1} is the standard Pauli-twirl identity.*

Limiting cases

The value of dynamical decoupling is decided entirely by how the pulse spacing compares with the bath correlation time τc\tau_c and by the shape of the noise spectrum S(ω)S(\omega). The four regimes below are the ones worth memorizing.

CaseResultMeaning
Slow / quasi-static dephasingτinterval≪τc\tau_{\text{interval}} \ll \tau_cThe ideal regime: a single echo removes static T2∗T_2^* broadening, and a CPMG train extends coherence toward 2T12T_1; suppression order grows with pulse number.
Fast / white (Markovian) noiseS(ω)≈const, τc→0S(\omega)\approx\text{const},\ \tau_c\to0No benefit: the bath decorrelates between any two pulses, the notch captures negligible spectral weight, and decay stays exponential at the Markovian rate.
Hard high-frequency cutoffS(ω)=0S(\omega)=0 for ω>ωc\omega>\omega_cUhrig DD wins: non-uniform timings give NN-th-order suppression, ∣FN∣2∼(ωT)2N+2|F_N|^2\sim(\omega T)^{2N+2}, with only NN pulses.
Soft 1/ωα1/\omega^\alpha spectraS(ω)∝ω−αS(\omega)\propto\omega^{-\alpha}Uniform CPMG usually beats UDD, because the filter must reject a broad band rather than approach a sharp edge; optimal spacing depends on α\alpha.

Worked examples

1 · Uhrig DD pulse timings (fully worked)

You want to protect a qubit over a total idle time T=100 μsT=100\ \mu\text{s} against a bath with a sharp high-frequency cutoff, using an N=5N=5 pulse Uhrig sequence. Where do the five π pulses go, and what suppression order do they buy?

Step 1. Apply equation (7) with N=5N=5, so 2N+2=122N+2=12 and tj=Tsin⁡2 ⁣(jπ/12)t_j=T\sin^2\!\big(j\pi/12\big). The arguments jπ/12j\pi/12 are 15°,30°,45°,60°,75°15°,30°,45°,60°,75° — a clean set to evaluate by hand.

Step 2. Compute each: sin⁡215°=0.0670\sin^2 15°=0.0670, sin⁡230°=0.25\sin^2 30°=0.25, sin⁡245°=0.5\sin^2 45°=0.5, sin⁡260°=0.75\sin^2 60°=0.75, sin⁡275°=0.9330\sin^2 75°=0.9330. Multiplying by T=100 μsT=100\ \mu\text{s} gives pulse times 6.70, 25.0, 50.0, 75.0, 93.3 μs.

Step 3. Sanity-check the structure. The timings are symmetric about T/2T/2: t1+t5=6.70+93.3=100=Tt_1+t_5=6.70+93.3=100=T, t2+t4=25+75=100t_2+t_4=25+75=100, and t3=50t_3=50 sits exactly at the midpoint — the hallmark of a time-symmetric (palindromic) cycle that also kills the first Magnus correction. Compared with a uniform 5-pulse CPMG (which suppresses only as (ωT)2(\omega T)^2), this placement makes the filter vanish as (ωT)12(\omega T)^{12} near DC. Five pulses at 6.70/25.0/50.0/75.0/93.3 μs give 6th-order coherence protection, ∣F∣2∝(ωT)12|F|^2\propto(\omega T)^{12}.

2 · The 2T12T_1 ceiling in a superconducting qubit (fully worked)

A transmon has T1=100 μsT_1=100\ \mu\text{s} and a pure-dephasing time Tϕ=40 μsT_\phi=40\ \mu\text{s} dominated by slow 1/f1/f flux noise. What is T2T_2 before decoupling, and how far can DD push it?

Step 1. Use equation (9): 1/T2=1/(2T1)+1/Tϕ=1/200+1/40=0.005+0.025=0.030 μs−11/T_2 = 1/(2T_1)+1/T_\phi = 1/200 + 1/40 = 0.005+0.025 = 0.030\ \mu\text{s}^{-1}, so T2≈33.3 μsT_2 \approx 33.3\ \mu\text{s} before any decoupling.

Step 2. DD attacks the slow part of TϕT_\phi. In the ideal limit it drives 1/Tϕ→01/T_\phi\to0, leaving 1/T2→1/(2T1)1/T_2\to 1/(2T_1), i.e. T2→2T1=200 μsT_2\to 2T_1 = 200\ \mu\text{s} — a factor of six improvement over the un-decoupled value.

Step 3. Read the ceiling. Even with infinitely many perfect pulses, T2T_2 cannot exceed 200 μs200\ \mu\text{s}, because relaxation changes populations and phase manipulation cannot restore a decayed excitation. To go further you must attack T1T_1 itself — better materials, lower temperature, Purcell filtering — not more pulses. DD lifts T2T_2 from 33.3 μs toward the hard bound 2T1=200 μs2T_1=200\ \mu\text{s}, and no further.

3 · Where the CPMG passband sits (noise spectroscopy)

An N=8N=8 pulse CPMG train runs over T=80 μsT=80\ \mu\text{s}. At which noise frequency is the qubit most sensitive, and how would you use that to measure a spectrum?

Step 1. For a uniform CPMG train the inter-pulse spacing is τ=T/N\tau=T/N, and the filter's first passband peaks near ω∗≈π/τ=Nπ/T\omega^*\approx\pi/\tau=N\pi/T, i.e. f∗≈N/(2T)f^*\approx N/(2T).

Step 2. Plug in: f∗≈8/(2×80 μs)=8/(160 μs)=50 kHzf^*\approx 8/(2\times80\ \mu\text{s}) = 8/(160\ \mu\text{s}) = \mathbf{50\ kHz}. The train is a narrow band-pass probe centered at 50 kHz.

Step 3. Sweep it. Changing τ\tau (via NN or TT) slides f∗f^* across the spectrum; the measured decay WW at each setting samples S(ω)S(\omega) near ω∗\omega^*, so a family of CPMG runs reconstructs the bath spectral density — the standard "DD noise spectroscopy" that tells you which sequence to deploy. The 8-pulse train probes S(ω)S(\omega) at 50 kHz; scanning τ\tau maps the whole spectrum.

Physics of the portrait

The whole subject grows from one accidental discovery in a magnetic-resonance lab: that a signal thought irretrievably lost to dephasing can be brought back with a single well-timed pulse.

Origins · 1 minErwin Hahn and how this was discovered
EH
Erwin Hahn
1921–2016 · nuclear magnetic resonance and the spin echo

Erwin Hahn discovered the spin echo in 1950, as a young researcher at the University of Illinois, and it was very nearly a surprise even to him. Studying free-induction decay in nuclear magnetic resonance, he found that a second radio-frequency pulse applied after the signal had apparently vanished caused it to reappear — an "echo" — as long after the second pulse as the second had followed the first. The lost coherence had not been destroyed; it had merely dephased into an invisible fan of spins precessing at slightly different rates because the static field was not perfectly uniform. The pulse reversed the fan, and the spins re-bunched. Hahn's insight was that the refocusing works for any static inhomogeneity, known or unknown, because the reversal is exact — what distinguishes irreversible decoherence from reversible dephasing is precisely whether the environment changed during the interval. That distinction is the conceptual seed of everything in this article. Carr and Purcell (1954), then Meiboom and Gill (1958), turned the single echo into robust multi-pulse trains; and in 1999 Viola, Knill and Lloyd generalized the idea into the modern theory of dynamical decoupling for open quantum systems, where a measured NMR trick became a design principle for quantum computers.

History
1950
— Hahn discovers the spin echo in nuclear magnetic resonance
1954
— Carr and Purcell build the multi-pulse echo train
1958
— Meiboom and Gill add phase-alternated pulses for robustness (CPMG)
1999
— Viola, Knill and Lloyd formulate dynamical decoupling for open quantum systems
2007
— Uhrig derives the optimal non-uniform pulse timings (UDD)

Going further

Texts. Nielsen & Chuang, Quantum Computation and Quantum Information, §8 (open systems and quantum operations) · Preskill, Lecture Notes on Quantum Computation, Ch. 3 (decoherence) · Wilde, Quantum Information Theory, Ch. 4 (quantum channels) · Watrous, The Theory of Quantum Information (channels and CP maps). Primary sources. E. L. Hahn, "Spin Echoes," Phys. Rev. 80, 580 (1950) · Carr & Purcell, Phys. Rev. 94, 630 (1954); Meiboom & Gill, Rev. Sci. Instrum. 29, 688 (1958) · Viola, Knill & Lloyd, "Dynamical Decoupling of Open Quantum Systems," Phys. Rev. Lett. 82, 2417 (1999) · Uhrig, "Keeping a Quantum Bit Alive by Optimized π-Pulse Sequences," Phys. Rev. Lett. 98, 100504 (2007). The key results in this article are mutually consistent: the echo cancellation ϕnet=0\phi_{\text{net}}=0 and the vanishing Hˉ(0)\bar H^{(0)}; the Pauli conjugation signs XZX=YZY=−ZXZX=YZY=-Z, IZI=ZZZ=ZIZI=ZZZ=Z summing to zero via the Pauli-twirl projection 14∑gg†Sg=12Tr(S)1\frac14\sum_g g^\dagger S g=\frac12\mathrm{Tr}(S)\mathbb{1}; the Hahn-echo filter ∣F(ωT)∣2=16sin⁡4(ωT/4)|F(\omega T)|^2=16\sin^4(\omega T/4) from N=1, t1=T/2N=1,\,t_1=T/2, with F(0)=0F(0)=0; the UDD timings tj=Tsin⁡2[jπ/(2N+2)]t_j=T\sin^2[j\pi/(2N+2)] and the (ωT)2N+2(\omega T)^{2N+2} low-frequency scaling; and the ceiling T2≤2T1T_2\le2T_1 from 1/T2=1/(2T1)+1/Tϕ1/T_2=1/(2T_1)+1/T_\phi. The worked numbers (UDD times 6.70/25/50/75/93.3 μs; T2:33.3→200T_2:33.3\to200 μs; CPMG f∗=50f^*=50 kHz) follow directly from these relations.

Builds on. Open Quantum Systems and Reduced Dynamics · Decoherence · Relaxation and Dephasing: T1 and T2 Times · Markovian vs Non-Markovian Dynamics. Leads to. Markovian vs Non-Markovian Dynamics · Master Equation Derivations: Born–Markov and Secular Approximations · Phase Damping Channel.

Ask the tutor: Ask me to draw the CPMG filter function for a given pulse number and overlay it on a 1/f1/f noise spectrum, and we will read off exactly how much coherence you gain.

How is dynamical decoupling different from quantum error correction?

They solve the same problem by opposite philosophies. Error correction is closed-loop and active: it spreads one logical qubit across many physical ones, measures syndromes with ancillas, and applies feedback to correct detected errors. Dynamical decoupling is open-loop and passive: no encoding, no ancillas, no measurement, no feedback — just a pre-scheduled train of unitary pulses that averages the environmental coupling toward zero. DD exploits the temporal structure of the noise (that it is slow); error correction exploits redundancy and works even against fast, uncorrelated errors. In practice they are complementary: DD lowers the physical error rate — especially the correlated, non-Markovian errors codes handle badly — helping bring fidelities below the fault-tolerance threshold that error correction then works from.

Why does a static frequency offset cancel but slow noise only partly cancel?

A perfectly static offset is the same in the first half of an echo and the second, so the sign flip induced by the π pulse makes the two contributions exactly equal and opposite — they cancel to all orders, which is why F(0)=0F(0)=0 exactly. Slow noise is almost static: it barely changes across one cycle, so most of it cancels, but the small change between the halves survives. Formally, the leftover scales with the time-derivative of the noise, i.e. with spectral weight at nonzero frequency, which the filter function F(ωT)F(\omega T) passes with a small but nonzero amplitude near DC. The faster the noise wanders relative to the pulse spacing, the more leaks through — and in the white-noise limit the cancellation fails entirely.

When should I use UDD instead of CPMG?

Use Uhrig DD when the bath has a genuinely hard high-frequency cutoff — the noise spectrum drops sharply to zero above some ωc\omega_c. There UDD's non-uniform timings make the filter vanish as (ωT)2N+2(\omega T)^{2N+2}, so each added pulse buys two extra orders of low-frequency suppression, and it beats uniform CPMG decisively. But most solid-state noise is soft: 1/f1/f-like spectra that fall off gradually and have significant weight over a broad band. Against those, uniform CPMG usually wins, because the filter must reject a whole band rather than sit below a sharp edge, and the even spacing does that better. The honest workflow is to measure S(ω)S(\omega) first (via DD noise spectroscopy) and then choose the sequence whose filter best notches the spectrum you actually have.

What sets the maximum number of pulses I can usefully apply?

Two competing effects. More pulses push the filter's first passband to higher frequency (ω∗∼Nπ/T\omega^*\sim N\pi/T) and deepen the low-frequency notch, so coherence improves — up to a point. But every real pulse has finite duration and calibration error, and those imperfections inject their own decoherence that accumulates with pulse number. Beyond an optimum NN, added pulse error outweighs added noise suppression and the coherence gets worse. This is exactly why robust sequences exist: XY-4 alternates XX and YY π pulses so flip-angle and detuning errors partly self-cancel; XY-8 adds time-reversal symmetry to kill further error orders; KDD uses phase-shifted composite pulses for broadband robustness. The best sequence balances decoupling order against accumulated pulse error for your hardware.

Can dynamical decoupling protect a running gate, not just an idle qubit?

Yes, and this is an active research direction called dynamically corrected gates. Plain DD protects a qubit during idle windows, where the intended evolution is the identity, so any pulse train that returns the control to identity is fair game. Performing a nontrivial gate while decoupling is harder, because the pulses must average the system–bath coupling to zero and implement the target unitary simultaneously. The trick is to design composite pulse shapes whose control propagator realizes the desired gate while its toggling frame still visits enough orientations to cancel the leading coupling. It works best against the same slow, low-frequency noise that ordinary DD defeats, and it is one of the ways decoupling is being folded directly into compiled circuits rather than inserted only between operations.

Is the filter-function picture equivalent to the average-Hamiltonian picture?

They are two windows onto the same physics. Average Hamiltonian theory is a time-domain operator statement: it organizes one cycle's propagator into an effective generator via the Magnus expansion and asks you to zero its terms order by order — ideal for designing sequences and proving universality. The filter function is a frequency-domain statistical statement: for Gaussian dephasing it gives the exact decay as an overlap of the spectrum with ∣F(ωT)∣2|F(\omega T)|^2 — ideal for predicting decay curves and noise spectroscopy. They agree where both apply: Hˉ(0)=0\bar H^{(0)}=0 for a symmetric train is exactly F(0)=0F(0)=0, and higher Magnus orders correspond to how fast FF rises out of the DC notch. Use whichever suits the question.

Before this

Where this leads