Average Hamiltonian Theory

Average Hamiltonian theory replaces a rapid, usually periodic control cycle by the effective Hamiltonian that generates its stroboscopic evolution. It began as the language of coherent averaging in nuclear magnetic resonance (NMR), where pulse sequences can suppress dipolar broadening while retaining selected spin interactions. The same construction now organizes composite pulses, Hamiltonian engineering, dynamical decoupling, and digital or analog quantum simulation.

“Average” here does not mean an expectation value in a state. It means an operator obtained from the logarithm of a one-cycle propagator. The ordinary time average is only its leading term:

$$\begin{aligned} U(T_c) &= \exp\left( -\frac{i}{\hbar}\overline H T_c \right), \\ \overline H &= \overline H^{(0)} +\overline H^{(1)} +\cdots. \end{aligned}$$

Required background. The Magnus expansion supplies the exponential series; rotating frames explains the exact change of representation; stroboscopic dynamics identifies what a cycle generator predicts. Dynamical decoupling relates pulse timing to noise spectra and filter functions.

Periodic Control Problem

Separate the Hamiltonian into an internal or drift part and an applied control:

$$H(t)=H_{\mathrm{int}}(t)+H_c(t).$$

Let the control propagator satisfy

$$i\hbar\dot U_c(t) = H_c(t)U_c(t), \qquad U_c(0)=I.$$

Factor the full propagator as

$$U(t)=U_c(t)\widetilde U(t).$$

Substitution into the Schrödinger equation gives the exact toggling-frame evolution

$$i\hbar\frac{d\widetilde U}{dt} = \widetilde H(t)\widetilde U(t),$$

with

$$\widetilde H(t) = U_c^\dagger(t) H_{\mathrm{int}}(t) U_c(t).$$

The strong control has disappeared from the generator because it is represented exactly by $U_c(t)$. Its effect is to rotate, switch, or otherwise modulate the internal Hamiltonian.

Cyclic control

Choose a cycle duration $T_c$ for which

$$U_c(T_c)=e^{i\phi_c}I.$$

The irrelevant scalar phase may be set to one. If the control waveform and internal Hamiltonian repeat, then $\widetilde H(t+T_c)=\widetilde H(t)$. Define the exact cycle Hamiltonian by

$$\begin{aligned} \widetilde U(T_c) &= \mathcal T \exp\left[ -\frac{i}{\hbar} \int_0^{T_c}\widetilde H(t)\,dt \right], \\ &= \exp\left( -\frac{i}{\hbar}\overline H T_c \right). \end{aligned}$$

At integer numbers of cycles,

$$\widetilde U(nT_c) = \left[\widetilde U(T_c)\right]^n = \exp\left( -\frac{i}{\hbar}\overline H\,nT_c \right).$$

The exact equality requires a consistent branch of the logarithm. Different branches shift quasienergies by integer multiples of $2\pi\hbar/T_c$. If $U_c(T_c)$ does not close, its net rotation must be included in the one-cycle propagator rather than silently discarded.

Average Hamiltonian theory is therefore naturally stroboscopic. Motion within a cycle is micromotion and generally cannot be reconstructed from $\overline H$ alone.

Piecewise Pulse Sequences

Ideal hard-pulse models divide the cycle into free intervals separated by instantaneous control unitaries. Let $P_j$ be pulse $j$, and define the cumulative control after that pulse by

$$G_j=P_jP_{j-1}\cdots P_1, \qquad G_0=I.$$

During a free interval of duration $\tau_j$ after $j$ pulses, a time-independent internal Hamiltonian becomes

$$\widetilde H_j = G_j^\dagger H_{\mathrm{int}}G_j.$$

Cycle closure requires $G_m=e^{i\phi_c}I$ after the final pulse, and

$$T_c=\sum_{j=0}^{m-1}\tau_j$$

in the zero-width pulse idealization.

Control pulses change the cumulative frame $G_j$, so the same internal Hamiltonian appears as $\widetilde H_j=G_j^\dagger H_{\mathrm{int}}G_j$ in interval $j$. The weighted sum gives $\overline H^{(0)}$; ordered commutators of different intervals give higher terms.

The toggling frame converts pulse design into Hamiltonian design. Instead of tracking a long product of pulse and free-evolution operators directly, one asks which rotated Hamiltonians appear, for how long, and in what order.

Finite pulses are not instantaneous boundaries. During a pulse of duration $\tau_p$, the internal Hamiltonian and control act simultaneously, so that interval must be included in $\widetilde H(t)$. The hard-pulse approximation requires at least

$$\frac{\tau_p\lVert H_{\mathrm{int}}\rVert}{\hbar}\ll 1$$

on the relevant subspace, together with adequate control-amplitude, detuning, leakage, and bandwidth checks.

Notation and experimental conventions for Ramsey, echo, CPMG, XY-family, and composite-pulse protocols are catalogued in Pulse Sequences.

Zeroth-Order Average Hamiltonian

The leading term is the cycle average of the toggling-frame Hamiltonian:

$$\overline H^{(0)} = \frac{1}{T_c} \int_0^{T_c}\widetilde H(t)\,dt.$$

For ideal piecewise-constant intervals,

$$\overline H^{(0)} = \sum_{j=0}^{m-1} \frac{\tau_j}{T_c}\widetilde H_j.$$

This is an operator-valued weighted average. It can:

• cancel an unwanted term by rotating it through signs or orientations whose weighted sum vanishes;
• retain a desired term that is invariant under every control operation;
• rescale an interaction by changing interval weights;
• transform an anisotropic coupling into one with a chosen symmetry;
• select one component of a larger operator decomposition.

The control cannot distinguish “desired” and “undesired” by intention. Any term transforms according to its operator structure. A sequence that removes dephasing may also remove a wanted signal or gate Hamiltonian unless that term lies in the retained symmetry sector.

Echo as exact commuting cancellation

For static qubit detuning,

$$H_{\mathrm{int}} = \frac{\hbar\delta}{2}Z,$$

an ideal $X$ pulse at half a cycle gives the following toggling Hamiltonian (apply a second $X$ at the endpoint to close the control cycle)

$$\widetilde H(t) = \begin{cases} +\hbar\delta Z/2, &0\lt t\lt T_c/2, \\ -\hbar\delta Z/2, &T_c/2\lt t\lt T_c. \end{cases}$$

The zeroth-order average vanishes. Because the two interval Hamiltonians are proportional to $Z$, all their commutators vanish as well. In this ideal static model, the toggling propagator is exactly the identity. Without the closing pulse the laboratory propagator still contains the net $X$ rotation; cancellation of detuning does not erase that intended action. Bath evolution, time-dependent noise, transverse couplings, and pulse imperfections remove that special exactness.

Higher-Order Corrections

Noncommuting interval Hamiltonians remember their order. Applying the Magnus series to one toggling-frame cycle gives

$$\overline H = \overline H^{(0)} +\overline H^{(1)} +\overline H^{(2)} +\cdots.$$

Average-Hamiltonian indexing starts the time average at order zero. The first correction is

$$\overline H^{(1)} ={} -\frac{i}{2\hbar T_c} \int_0^{T_c}dt_1 \int_0^{t_1}dt_2 \times \left[ \widetilde H(t_1), \widetilde H(t_2) \right].$$

The second correction is

$$\overline H^{(2)} ={} -\frac{1}{6\hbar^2T_c} \int_0^{T_c}dt_1 \int_0^{t_1}dt_2 \int_0^{t_2}dt_3 \times \Big( [\widetilde H_1,[\widetilde H_2,\widetilde H_3]] + [\widetilde H_3,[\widetilde H_2,\widetilde H_1]] \Big),$$

where $\widetilde H_j=\widetilde H(t_j)$ only within this displayed formula.

For piecewise-constant free intervals, the first correction simplifies to

$$\overline H^{(1)} = -\frac{i}{2\hbar T_c} \sum_{j>k} \tau_j\tau_k [\widetilde H_j,\widetilde H_k].$$

There is no contribution from two times in the same interval because $[\widetilde H_j,\widetilde H_j]=0$. Reversing the interval order changes the sign of each pairwise commutator while leaving $\overline H^{(0)}$ unchanged. Two pulse cycles can therefore have the same time fractions and different effective Hamiltonians.

Expansion parameter

If $\lVert\widetilde H(t)\rVert$ is characterized by an energy scale $h$, then the natural cycle parameter is

$$\lambda=\frac{hT_c}{\hbar}.$$

Dimensionally,

$$\overline H^{(n)} \sim h\lambda^n.$$

Fast cycling means $\lambda\ll1$ on the retained subspace. A standard sufficient Magnus convergence condition is

$$\frac{1}{\hbar} \int_0^{T_c} \lVert\widetilde H(t)\rVert\,dt \lt \pi.$$

It is sufficient rather than necessary, and it can be unhelpful for unbounded or extensive many-body Hamiltonians. In those settings, locality, energy density, resonances, and prethermal time windows can matter more than the global norm. A truncated average Hamiltonian should be validated against the exact one-cycle propagator or a converged numerical calculation whenever feasible.

Accumulation over many cycles

An exact $\overline H$ reproduces all integer-cycle evolution. A truncation leaves a one-cycle error that can accumulate with the number of repetitions. For unitary cycle maps $U,V$, telescoping gives the rigorous operator-norm bound $\|U^n-V^n\|\le\min(2,n\|U-V\|)$. Thus a per-cycle error $\epsilon_c$ guarantees at most $n\epsilon_c$ after $n$ repetitions; cancellations may improve the actual result. Small $T_c$ alone does not certify arbitrarily long evolution.

Sequence Symmetry

Symmetry can remove whole classes of commutator corrections without evaluating them term by term.

If the complete toggling-frame cycle is symmetric about its midpoint,

$$\widetilde H(t) = \widetilde H(T_c-t),$$

then the odd average-Hamiltonian orders vanish:

$$\overline H^{(1)} = \overline H^{(3)} = \overline H^{(5)} = \cdots = 0.$$

In Magnus notation these correspond to the even Magnus exponents $\Omega_2,\Omega_4,\ldots$. The distinction is purely indexing: $\overline H^{(0)}$ comes from $\Omega_1$.

A practical way to create this symmetry is a palindromic cycle: traverse a set of toggling-frame intervals and then retrace them in reverse order with matching durations. The statement applies to the full Hamiltonian during the full cycle. Treating pulses as finite can break a symmetry that appears exact in a delta-pulse diagram unless the pulse shapes are also arranged symmetrically.

Time antisymmetry and phase-cycling symmetries can impose other selection rules, but they must be derived for the actual toggling-frame Hamiltonian rather than inferred from a visually balanced pulse diagram.

Group Averaging and the Commutant

Suppose the control visits every element of a finite unitary group $\mathcal G$ for equal durations. Zeroth-order averaging defines

$$\Pi_{\mathcal G}(H) = \frac{1}{\lvert\mathcal G\rvert} \sum_{g\in\mathcal G} g^\dagger Hg.$$

For any $h\in\mathcal G$,

$$h^\dagger \Pi_{\mathcal G}(H) h = \Pi_{\mathcal G}(H).$$

Thus the group average lies in the commutant of the control representation: it retains exactly the operator components invariant under conjugation by the group. With the Hilbert–Schmidt inner product in finite dimensions, $\Pi_{\mathcal G}$ is the projector onto that invariant operator subspace.

Single-qubit Pauli average

Write a general qubit–bath interaction as

$$H ={} I\otimes B_0 +X\otimes B_x + Y\otimes B_y +Z\otimes B_z.$$

Average over

$$\mathcal G=\{I,X,Y,Z\},$$

with global phases ignored. Conjugation by the Pauli operators changes the signs of the nonidentity Pauli components. Their four contributions cancel, leaving

$$\Pi_{\mathcal G}(H) = I\otimes B_0.$$

This is the algebraic core of universal first-order qubit decoupling. It does not by itself guarantee practical coherence protection: higher commutators, finite pulses, bath evolution, control errors, and noise above the cycle rate remain.

Selective Averaging in NMR

For two spins $I$ and $S$, a secular dipolar coupling aligned with the laboratory $z$ axis has the schematic form

$$H_{D,z} = d\left( 2I_zS_z -I_xS_x -I_yS_y \right).$$

Ideal rotations can permute the distinguished axis, producing

$$\begin{aligned} H_{D,x} &= d\left( 2I_xS_x-I_yS_y-I_zS_z \right), \\ H_{D,y} &= d\left( 2I_yS_y-I_zS_z-I_xS_x \right). \end{aligned}$$

The three orientations obey

$$H_{D,x}+H_{D,y}+H_{D,z}=0.$$

Equal time in these orientations therefore removes this ideal dipolar interaction from $\overline H^{(0)}$. This identity is an algebraic core of multiple-pulse dipolar averaging in solid-state NMR. It is not a complete specification of WAHUHA or any laboratory sequence: real designs must also track chemical shifts, scalar couplings, pulse phases, receiver timing, finite widths, and higher-order terms.

The broader NMR goal is selective averaging. One engineers a control cycle whose invariant operator sector contains useful chemical-shift or coupling information while strong broadening terms average away. Average Hamiltonian theory supplies both the leading selection rule and the corrections that limit spectral resolution.

Choosing the Cycle Origin

The exact cycle Hamiltonian is also a Floquet Hamiltonian. Changing the cycle origin $t_0$ changes

$$U(t_0+T_c,t_0)$$

by a micromotion conjugation. Exact quasienergies are unchanged, but the matrix representation of the effective Hamiltonian and any finite-order truncation can depend on that choice. A reported average Hamiltonian should therefore state the frame, cycle origin, and logarithm branch.

Dynamical Decoupling Preview

Dynamical decoupling uses control cycles to suppress system–environment coupling. At zeroth order, group or sign averaging can remove selected error operators. Higher average-Hamiltonian terms diagnose coherent residual couplings generated by noncommuting toggling-frame pieces.

That operator calculation is only one view of decoupling. For classical or quantum noise with temporal structure, pulse timing determines a frequency-domain filter. A sequence with a vanishing static average can still transmit noise near its filter peaks, and no finite-bandwidth control removes arbitrary Markovian relaxation. The canonical treatment of spin echo, CPMG, Uhrig timing, filter functions, and sensing tradeoffs is Dynamical Decoupling.

For applications inside quantum circuits, dynamical decoupling in quantum information explains how pulse placement interacts with gates and idle intervals.

Use average Hamiltonian theory when the cycle is rapid enough that coherent operator corrections are the natural diagnostic. Use filter-function or open-system methods when the noise spectrum, bath correlation time, dissipation, or measured coherence envelope is central. Many experiments require both.

Exercises

1. Symmetrize a two-step cycle

Consider the palindromic cycle

$$H_A\ \text{for }\tau/2, \quad H_B\ \text{for }\tau, \quad H_A\ \text{for }\tau/2.$$

Show directly from the piecewise formula that $\overline H^{(1)}=0$.

Solution

Label the three intervals $0,1,2$ with

$$\begin{aligned} H_0&=H_A, & H_1&=H_B, \\ H_2&=H_A. \end{aligned}$$

and durations $\tau_0=\tau_2=\tau/2$, $\tau_1=\tau$. The pair $(1,0)$ contributes

$$\frac{\tau^2}{2}[H_B,H_A].$$

The pair $(2,1)$ contributes

$$\frac{\tau^2}{2}[H_A,H_B] = -\frac{\tau^2}{2}[H_B,H_A].$$

The pair $(2,0)$ vanishes because $[H_A,H_A]=0$. The sum is zero, so $\overline H^{(1)}=0$. This is the first nontrivial instance of the general time-symmetry rule.

2. Pauli-group decoupling

Verify explicitly that

$$\frac14 \sum_{g\in\{I,X,Y,Z\}} g^\dagger Xg=0,$$

and similarly for $Y$ and $Z$. What happens to a desired qubit Hamiltonian proportional to $Z$ under the same average?

Solution

Conjugating $X$ gives

$$\begin{array}{c|cccc} g & I & X & Y & Z\\ \hline g^\dagger Xg & X & X & -X & -X \end{array}$$

so the average vanishes. Cyclic permutations give the same cancellation for $Y$ and $Z$. Therefore every traceless single-qubit operator is removed, while the identity survives.

A desired Hamiltonian proportional to $Z$ is also removed. The group average recognizes operator symmetry, not experimental intention. To preserve a nontrivial gate while decoupling errors, one must modify the control construction, encode the system, or use dynamically corrected gates.

3. The special case of two spin-one-half particles

For one isolated spin-one-half pair, put $I_\alpha=\sigma_\alpha/2$ and $S_\alpha=\sigma_\alpha/2$ on their respective factors. Do the three ideal dipolar orientations above commute? What changes for a network with overlapping bonds?

Solution

For distinct axes, $[\sigma_x\otimes\sigma_x,\sigma_y\otimes\sigma_y]=0$, and cyclic permutations give the other commutators. Therefore the three $H_{D,\alpha}$ commute for this isolated pair. Equal-duration evolution in all three orientations cancels exactly in the ideal toggling frame, since their sum is zero.

For overlapping bonds, however,

$$[\sigma_1^x\sigma_2^x,\sigma_2^y\sigma_3^y] =2i\sigma_1^x\sigma_2^z\sigma_3^y\ne0.$$

A spin network can therefore retain higher-order terms even though each bond’s zeroth-order average vanishes. Finite pulses and additional interactions also spoil the isolated-pair exactness.

4. The drift during a finite π pulse

A resonant rectangular pulse has $H_c=\hbar\Omega X/2$, duration $\tau_p=\pi/\Omega$, and a static drift $H_{\mathrm{int}}=\hbar\delta Z/2$. Find the drift’s toggling-frame time average over the pulse.

Solution

The control rotation gives $U_c^\dagger ZU_c=Z\cos\Omega t+Y\sin\Omega t$. Therefore

$$\frac1{\tau_p}\int_0^{\tau_p}\widetilde H(t)\,dt =\frac{\hbar\delta}{\pi}Y.$$

The cosine averages to zero but the sine does not. The leading error exponent has magnitude $|\delta|/\Omega$, even when the pulse area is exactly $\pi$. Treating the pulse as instantaneous would miss this residual rotation.

References

• S. Blanes, F. Casas, J. A. Oteo, and J. Ros, “The Magnus expansion and some of its applications,” Physics Reports 470, 151–238 (2009), doi:10.1016/j.physrep.2008.11.001.
• U. Haeberlen and J. S. Waugh, “Coherent averaging effects in magnetic resonance,” Physical Review 175, 453–467 (1968), doi:10.1103/PhysRev.175.453.
• U. Haeberlen, High Resolution NMR in Solids: Selective Averaging, Academic Press (1976).
• W. Magnus, “On the exponential solution of differential equations for a linear operator,” Communications on Pure and Applied Mathematics 7, 649–673 (1954), doi:10.1002/cpa.3160070404.
• M. M. Maricq, “Application of average Hamiltonian theory to the NMR of solids,” Physical Review B 25, 6622–6632 (1982), doi:10.1103/PhysRevB.25.6622.
• L. Viola, E. Knill, and S. Lloyd, “Dynamical decoupling of open quantum systems,” Physical Review Letters 82, 2417–2421 (1999), doi:10.1103/PhysRevLett.82.2417.
• L. Viola, S. Lloyd, and E. Knill, “Universal control of decoupled quantum systems,” Physical Review Letters 83, 4888–4891 (1999), doi:10.1103/PhysRevLett.83.4888.
• J. S. Waugh, L. M. Huber, and U. Haeberlen, “Approach to high-resolution NMR in solids,” Physical Review Letters 20, 180–182 (1968), doi:10.1103/PhysRevLett.20.180.