Two-Level Dynamics & Rabi Oscillations

4 hours ~9 min read

Two-Level Dynamics & Rabi Oscillations

This is the last lesson of Term 1, and it is where the physics becomes the machine. We take the propagator U(t)=eiHtU(t)=e^{-iHt} from the previous lesson and apply it to the one Hamiltonian every qubit platform implements: a driven two-level system. Out falls the Rabi oscillation — population sloshing coherently between 0|0\rangle and 1|1\rangle — and, with it, the punchline of the entire term: a quantum gate is just time evolution run for a chosen duration. Driving H=Ω2XH=\tfrac{\Omega}{2}X for the right time turns 0|0\rangle into 1|1\rangle; that pulse is the XX gate. Everything in Term 2 (gates and circuits) and Term 4.4 (how superconducting qubits actually do this) rests here.

Learning Objectives

After this lesson you will be able to:

  1. Write the two-level drive Hamiltonian and identify the on-resonance limit H=Ω2XH = \tfrac{\Omega}{2}X.
  2. Derive the Rabi formula P1(t)=sin2(Ωt/2)P_1(t) = \sin^2(\Omega t/2) for a qubit starting in 0|0\rangle, two independent ways (closed-form exponential and energy eigenbasis).
  3. Generalize to nonzero detuning: P1(t)=Ω2Ωeff2sin2(Ωefft/2)P_1(t) = \frac{\Omega^2}{\Omega_{\text{eff}}^2}\sin^2(\Omega_{\text{eff}}t/2) with Ωeff=Ω2+Δ2\Omega_{\text{eff}} = \sqrt{\Omega^2 + \Delta^2}.
  4. Interpret Rabi dynamics as Bloch-vector precession about the drive axis.
  5. Show that evolving under Ω2X\tfrac{\Omega}{2}X for time tt realizes Rx(Ωt)R_x(\Omega t), so a π\pi-pulse is an XX gate — connecting analog dynamics to the digital gate model.

Intuition

Picture a qubit as a tiny magnet (a Bloch vector) sitting in a static field along z^\hat z that splits the energies of 0|0\rangle and 1|1\rangle by ω0\omega_0. Left alone, the vector just precesses about z^\hat z and the populations never change0|0\rangle and 1|1\rangle are energy eigenstates (this was Example 1 of the previous lesson). To move population between the levels you have to push the qubit off the zz-axis: you apply an oscillating drive (a microwave or laser tone) that, viewed in the right rotating frame, becomes a constant field along an axis in the equatorial plane. The Bloch vector then precesses about that tilted axis — and as it swings, its z^\hat z-projection runs from 0|0\rangle down to 1|1\rangle and back. That coherent round-trip of population is a Rabi oscillation.

The drive strength Ω\Omega (the Rabi frequency) sets how fast the vector swings. If the drive is perfectly on resonance with the qubit's energy gap, the rotation axis lies flat on the equator and the population swings all the way to 1|1\rangle. If the drive is detuned, the axis tilts up out of the plane, the precession cone is narrower, and the qubit never fully reaches 1|1\rangle. So the two knobs are: how hard you drive (amplitude Ω\Omega) and how far off you are in frequency (detuning Δ\Delta). Stop the drive after exactly half a swing and you have flipped the qubit — that is an XX gate, built out of physics.


Theory

Symbol note. HH is the Hamiltonian here, never the Hadamard. =1\hbar = 1, rotation gates use the half-angle convention Rσ(θ)=eiθσ/2R_\sigma(\theta)=e^{-i\theta\sigma/2} (Appendix C, E).

The driven two-level Hamiltonian

A qubit with bare energy splitting ω0\omega_0 has free Hamiltonian H0=ω02ZH_0 = \tfrac{\omega_0}{2}Z (energies ±ω02\pm\tfrac{\omega_0}{2} for 0,1|0\rangle,|1\rangle). Couple it to a classical field oscillating at drive frequency ωd\omega_d with amplitude (Rabi frequency) Ω\Omega; the coupling is a transverse term:

H(t)=ω02Z+Ωcos(ωdt)X. H(t) = \frac{\omega_0}{2}Z + \Omega\cos(\omega_d t)\,X .

This is time-dependent, so the clean eiHte^{-iHt} machinery does not apply directly. The standard move is to pass to the frame rotating at the drive frequency ωd\omega_d — i.e. substitute ψ=eiωd2Ztψ~|\psi\rangle = e^{-i\frac{\omega_d}{2}Zt}|\tilde\psi\rangle — and drop the fast-oscillating counter-rotating term (the rotating-wave approximation, RWA). The transformed Hamiltonian becomes time-independent (derivation deferred to Term 4.4, where it is the working model of a real device):

HRWA=Δ2Z+Ω2X,Δ:=ω0ωd \boxed{\,H_{\text{RWA}} = \frac{\Delta}{2}\,Z + \frac{\Omega}{2}\,X,\qquad \Delta := \omega_0 - \omega_d\,}

with detuning Δ\Delta (how far the drive is off resonance) and Rabi frequency Ω\Omega (how hard we drive). This single, constant, Hermitian 2×22\times2 operator generates all the dynamics below. The whole point of the rotating frame is that it converts an oscillating drive into the kind of constant-HH problem we solved in 1.6.1.

On resonance: deriving P1(t)=sin2(Ωt/2)P_1(t)=\sin^2(\Omega t/2)

Set the drive exactly on resonance, Δ=0\Delta = 0. Then

H=Ω2X. H = \frac{\Omega}{2}X .

Start in ψ(0)=0|\psi(0)\rangle = |0\rangle and apply the propagator. Because X2=IX^2 = I, the exponential closes into trig (Appendix E, §3):

U(t)=eiΩ2Xt=cos ⁣Ωt2Iisin ⁣Ωt2X=(cosΩt2isinΩt2isinΩt2cosΩt2). U(t) = e^{-i\frac{\Omega}{2}Xt} = \cos\!\frac{\Omega t}{2}\,I - i\sin\!\frac{\Omega t}{2}\,X = \begin{pmatrix}\cos\frac{\Omega t}{2} & -i\sin\frac{\Omega t}{2}\\[2pt] -i\sin\frac{\Omega t}{2} & \cos\frac{\Omega t}{2}\end{pmatrix}.

Apply it to 0|0\rangle:

ψ(t)=U(t)0=cos ⁣Ωt20isin ⁣Ωt21. |\psi(t)\rangle = U(t)|0\rangle = \cos\!\frac{\Omega t}{2}\,|0\rangle - i\sin\!\frac{\Omega t}{2}\,|1\rangle .

The Born rule (Appendix C) gives the excited-state population:

P1(t)=1ψ(t)2=sin2 ⁣Ωt2=1cosΩt2 \boxed{\,P_1(t) = |\langle 1|\psi(t)\rangle|^2 = \sin^2\!\frac{\Omega t}{2} = \frac{1-\cos\Omega t}{2}\,}

and P0(t)=cos2(Ωt/2)P_0(t) = \cos^2(\Omega t/2), so P0+P1=1P_0 + P_1 = 1 at all times (unitarity). The population oscillates fully between 00 and 11 at angular frequency Ω\Omega — these are Rabi oscillations. The period is T=2π/ΩT = 2\pi/\Omega; the qubit is fully inverted (P1=1P_1=1) at t=π/Ωt = \pi/\Omega.

Second derivation (energy eigenbasis), as a cross-check. The eigenstates of H=Ω2XH=\tfrac{\Omega}{2}X are ±|\pm\rangle with energies E±=±Ω2E_\pm = \pm\tfrac{\Omega}{2}. Expand 0=12(++)|0\rangle = \tfrac1{\sqrt2}(|+\rangle+|-\rangle) and evolve each by its phase (the recipe from 1.6.1):

ψ(t)=12(eiΩt/2++e+iΩt/2). |\psi(t)\rangle = \tfrac1{\sqrt2}\big(e^{-i\Omega t/2}|+\rangle + e^{+i\Omega t/2}|-\rangle\big).

Now 1±=±12\langle 1|\pm\rangle = \pm\tfrac1{\sqrt2}, so

1ψ(t)=12(12eiΩt/212e+iΩt/2)=12(eiΩt/2e+iΩt/2)=isin ⁣Ωt2, \langle 1|\psi(t)\rangle = \tfrac1{\sqrt2}\Big(\tfrac{1}{\sqrt2}e^{-i\Omega t/2} - \tfrac{1}{\sqrt2}e^{+i\Omega t/2}\Big) = \tfrac12\big(e^{-i\Omega t/2}-e^{+i\Omega t/2}\big) = -i\sin\!\frac{\Omega t}{2},

giving P1(t)=sin2(Ωt/2)P_1(t) = \sin^2(\Omega t/2) — identical. ∎ The oscillation frequency Ω=E+E\Omega = E_+ - E_- is exactly the Bohr frequency (energy difference) of the drive Hamiltonian, as the general theory predicted.

Off resonance: detuning and the effective Rabi frequency

Restore the detuning, H=Δ2Z+Ω2X=12nσH = \tfrac{\Delta}{2}Z + \tfrac{\Omega}{2}X = \tfrac12\,\vec n\cdot\vec\sigma with n=(Ω,0,Δ)\vec n = (\Omega,\,0,\,\Delta). Write Ωeff=n=Ω2+Δ2\Omega_{\text{eff}} = |\vec n| = \sqrt{\Omega^2+\Delta^2} and the unit axis n^=n/Ωeff\hat n = \vec n/\Omega_{\text{eff}}. Then H=Ωeff2(n^σ)H = \tfrac{\Omega_{\text{eff}}}{2}\,(\hat n\cdot\vec\sigma), and since (n^σ)2=I(\hat n\cdot\vec\sigma)^2 = I the half-angle exponential applies again (Appendix E, §3):

U(t)=eiΩeff2(n^σ)t=cos ⁣Ωefft2Iisin ⁣Ωefft2(n^σ). U(t) = e^{-i\frac{\Omega_{\text{eff}}}{2}(\hat n\cdot\vec\sigma)\,t} = \cos\!\frac{\Omega_{\text{eff}}t}{2}\,I - i\sin\!\frac{\Omega_{\text{eff}}t}{2}\,(\hat n\cdot\vec\sigma).

The amplitude to land in 1|1\rangle from 0|0\rangle comes from the off-diagonal piece of n^σ\hat n\cdot\vec\sigma, which carries the factor Ω/Ωeff\Omega/\Omega_{\text{eff}}. Squaring gives the generalized Rabi formula:

P1(t)=Ω2Ω2+Δ2sin2 ⁣(Ωefft2),Ωeff=Ω2+Δ2 \boxed{\,P_1(t) = \frac{\Omega^2}{\Omega^2+\Delta^2}\,\sin^2\!\Big(\tfrac{\Omega_{\text{eff}}\,t}{2}\Big), \qquad \Omega_{\text{eff}} = \sqrt{\Omega^2+\Delta^2}\,}

Read off the physics: detuning makes the oscillation faster (frequency ΩeffΩ\Omega_{\text{eff}}\ge\Omega) but shallower — the maximum reachable population is Ω2/(Ω2+Δ2)<1\Omega^2/(\Omega^2+\Delta^2) < 1. You can no longer fully flip the qubit off resonance. Resonance (Δ=0\Delta=0) is special precisely because it is the only setting that achieves complete population transfer. (Sweeping Δ\Delta and watching this contrast collapse is how experimentalists find a qubit's frequency — a Rabi/chevron measurement, Term 4.4.)

Bloch-sphere picture: precession about the drive axis

The Liouville/Ehrenfest machinery from 1.6.1 makes the geometry exact. For a single qubit, H=12nσH = \tfrac12\vec n\cdot\vec\sigma and ρ=12(I+rσ)\rho=\tfrac12(I+\vec r\cdot\vec\sigma) (Appendix E, §1). Feeding ρ˙=i[H,ρ]\dot\rho = -i[H,\rho] through the Pauli commutators [σj,σk]=2iεjklσl[\sigma_j,\sigma_k]=2i\varepsilon_{jkl}\sigma_l yields the Bloch equation:

drdt=n×r \boxed{\,\frac{d\vec r}{dt} = \vec n\times\vec r\,}

This is exactly the equation of a vector precessing about the axis n\vec n at angular speed n=Ωeff|\vec n| = \Omega_{\text{eff}} — the classical-magnet picture is literally correct. On resonance n=(Ω,0,0)\vec n = (\Omega,0,0) lies on the equator along x^\hat x: starting from the north pole r(0)=(0,0,1)\vec r(0)=(0,0,1) (=0=|0\rangle), the vector precesses in the yyzz plane down through the south pole (1|1\rangle) and back. Its height is rz(t)=cosΩt=12P1(t)r_z(t) = \cos\Omega t = 1 - 2P_1(t), reproducing P1=sin2(Ωt/2)P_1=\sin^2(\Omega t/2). With detuning the axis tilts up to n=(Ω,0,Δ)\vec n=(\Omega,0,\Delta); the precession cone no longer reaches the south pole, so rzr_z never hits 1-1 — the geometric reason the population transfer is incomplete.

The key computing connection: a gate is timed evolution

Compare the on-resonance propagator with the rotation gate from Appendix E:

U(t)=eiΩ2Xt=ei(Ωt)X/2=Rx(Ωt). U(t) = e^{-i\frac{\Omega}{2}Xt} = e^{-i\,(\Omega t)\,X/2} = R_x(\Omega t).

So letting the qubit evolve under the drive Ω2X\tfrac{\Omega}{2}X for a duration tt implements the gate Rx(θ)R_x(\theta) with rotation angle θ=Ωt\theta = \Omega t. The "gate angle" is the Rabi frequency times the pulse length — you dial it with time. Special pulses:

| Pulse | Duration tt | Angle θ=Ωt\theta=\Omega t | Gate Rx(θ)R_x(\theta) | Effect on 0|0\rangle | |---|---|---|---|---| | π/2\pi/2-pulse | π/(2Ω)\pi/(2\Omega) | π/2\pi/2 | Rx(π/2)R_x(\pi/2) | 12(0i1)\to \tfrac1{\sqrt2}(|0\rangle - i|1\rangle) (equator) | | π\pi-pulse | π/Ω\pi/\Omega | π\pi | Rx(π)=iXR_x(\pi) = -iX | i11\to -i|1\rangle \equiv |1\rangle | | 2π2\pi-pulse | 2π/Ω2\pi/\Omega | 2π2\pi | Rx(2π)=IR_x(2\pi) = -I | 00\to -|0\rangle \equiv |0\rangle |

A π\pi-pulse is an XX gate (NOT gate): Rx(π)=cosπ2Iisinπ2X=iXR_x(\pi) = \cos\tfrac\pi2 I - i\sin\tfrac\pi2 X = -iX, which equals XX up to the global phase i-i — and global phase is physically invisible (Appendix C). It sends 01|0\rangle\to|1\rangle, the very first peak of the Rabi curve P1(π/Ω)=sin2(π/2)=1P_1(\pi/\Omega)=\sin^2(\pi/2)=1. This is not an analogy — it is how XX is physically realized on real hardware, by sending a calibrated π\pi-pulse. The "digital" gate model of Term 2.1 is an abstraction layered on top of exactly this analog control; how the pulse is engineered on a superconducting qubit (and what limits its fidelity) is the subject of Term 4.4. Choosing a different drive axis (a YY-quadrature drive gives Ω2YRy\tfrac{\Omega}{2}Y\Rightarrow R_y; the Δ\Delta term gives RzR_z) lets a single physical mechanism realize any single-qubit rotation — which, by the Euler/ZYZ decomposition (Appendix E, §3), is all of SU(2)SU(2).


Worked Examples

Example 1 — A π\pi-pulse is an XX gate

Drive on resonance, H=Ω2XH=\tfrac{\Omega}{2}X, for the time tπ=π/Ωt_\pi = \pi/\Omega (so θ=Ωtπ=π\theta = \Omega t_\pi = \pi):

U(tπ)=cos ⁣π2Iisin ⁣π2X=iX=(0ii0). U(t_\pi) = \cos\!\tfrac\pi2\,I - i\sin\!\tfrac\pi2\,X = -iX = \begin{pmatrix}0 & -i\\ -i & 0\end{pmatrix}.

Acting on 0|0\rangle: U(tπ)0=i1U(t_\pi)|0\rangle = -i|1\rangle. Up to the irrelevant global phase i-i, this is 1|1\rangle — a complete bit flip, with P1(tπ)=sin2(π/2)=1P_1(t_\pi) = \sin^2(\pi/2) = 1. Acting on 1|1\rangle: U(tπ)1=i00U(t_\pi)|1\rangle = -i|0\rangle \equiv |0\rangle. So the π\pi-pulse swaps 01|0\rangle\leftrightarrow|1\rangle exactly as the XX (NOT) gate does. Since X=iRx(π)X = i\,R_x(\pi), the standard XX gate and the physical π\pi-pulse differ only by global phase — operationally the same gate.

Example 2 — Detuned drive: incomplete transfer

Drive with Ω=1\Omega = 1 and detuning Δ=1\Delta = 1 (so the drive is off resonance by exactly the Rabi frequency). Then Ωeff=12+12=2\Omega_{\text{eff}} = \sqrt{1^2+1^2} = \sqrt2 and the maximum reachable excited population is

P1max=Ω2Ω2+Δ2=11+1=12. P_1^{\max} = \frac{\Omega^2}{\Omega^2+\Delta^2} = \frac{1}{1+1} = \frac12 .

The full curve is P1(t)=12sin2 ⁣(22t)P_1(t) = \tfrac12\sin^2\!\big(\tfrac{\sqrt2}{2}t\big). Compared to the resonant case it oscillates faster (frequency 2\sqrt2 instead of 11) but only ever reaches halfway up — a π\pi-pulse duration no longer produces a clean bit flip. This is why qubit drives must be carefully calibrated to resonance: a control error in Δ\Delta both speeds up and caps the rotation, corrupting the intended gate. (Detuning errors are a standard coherent-error model analyzed in Term 4.)


Hands-on (Python)

We simulate the dynamics with scipy.linalg.expm, plot the Rabi oscillation, verify the π\pi-pulse realizes XX, and overlay a detuned curve to see incomplete transfer.

import numpy as np
from scipy.linalg import expm
import matplotlib.pyplot as plt

# Pauli matrices (Appendix E)
I = np.eye(2, dtype=complex)
X = np.array([[0, 1], [1, 0]], dtype=complex)
Z = np.array([[1, 0], [0, -1]], dtype=complex)

ket0 = np.array([1, 0], dtype=complex)
ket1 = np.array([0, 1], dtype=complex)

def H_drive(Omega, Delta):
    """Rotating-frame two-level Hamiltonian H = (Δ/2)Z + (Ω/2)X."""
    return 0.5 * Delta * Z + 0.5 * Omega * X

def P1(t, Omega, Delta):
    """Excited-state population starting from |0>, via U(t)=exp(-iHt)."""
    U = expm(-1j * H_drive(Omega, Delta) * t)
    psi = U @ ket0
    return np.abs(ket1 @ psi) ** 2          # Born rule |<1|ψ(t)>|²

# --- Resonant Rabi oscillations: P1(t) = sin²(Ωt/2) ---
Omega = 2.0
ts = np.linspace(0, 3 * np.pi / Omega, 500)
P1_num   = np.array([P1(t, Omega, Delta=0.0) for t in ts])
P1_exact = np.sin(Omega * ts / 2) ** 2

print("resonant: max |numeric - sin²(Ωt/2)| =", np.max(np.abs(P1_num - P1_exact)))  # ~1e-15
# --- The π-pulse maps |0> -> |1> and equals X up to global phase ---
t_pi = np.pi / Omega                         # θ = Ω t_pi = π
U_pi = expm(-1j * H_drive(Omega, 0.0) * t_pi)
psi_pi = U_pi @ ket0
print("ψ(t_π) =", np.round(psi_pi, 3))       # [0.-0.j  0.-1.j] = -i|1>
print("P1(t_π) =", round(P1(t_pi, Omega, 0.0), 12))   # 1.0  (full inversion)

# U_pi equals X up to a global phase: U_pi = -i X  ⇒  i*U_pi == X
print("U_π == X up to global phase? ", np.allclose(1j * U_pi, X))   # True
# --- Detuning: faster but shallower oscillations (incomplete transfer) ---
fig, ax = plt.subplots(figsize=(7.5, 3.4))
for Delta, style in [(0.0, "-"), (Omega, "--"), (2 * Omega, ":")]:
    Oeff = np.sqrt(Omega**2 + Delta**2)
    cap = Omega**2 / (Omega**2 + Delta**2)   # P1_max = Ω²/Ω_eff²
    curve = np.array([P1(t, Omega, Delta) for t in ts])
    ax.plot(ts, curve, style, lw=2,
            label=fr"$\Delta={Delta/Omega:.0f}\,\Omega$  ($P_1^{{\max}}={cap:.2f}$)")

ax.axhline(1.0, color="gray", ls=":", lw=0.8)
ax.set_xlabel("time  t"); ax.set_ylabel(r"$P_1(t)$")
ax.set_title(r"Rabi oscillations: resonant vs. detuned ($\Omega=2$)")
ax.legend(loc="upper right"); fig.tight_layout()
plt.show()
# On resonance the population reaches 1 (a clean X gate at t=π/Ω);
# detuned drives oscillate faster (Ω_eff) but cap below 1.
# --- The gate view: evolving (Ω/2)X for time t IS the gate R_x(Ωt) ---
def Rx(theta):                               # Appendix E half-angle convention
    return np.cos(theta/2) * I - 1j * np.sin(theta/2) * X

t = 0.37
U_t = expm(-1j * H_drive(Omega, 0.0) * t)
print("U(t) == R_x(Ω t)? ", np.allclose(U_t, Rx(Omega * t)))   # True

Braket note (optional). Because eiΩ2Xt=Rx(Ωt)e^{-i\frac{\Omega}{2}Xt} = R_x(\Omega t), the identical dynamics is a single gate in the SDK — from braket.circuits import Circuit; Circuit().rx(0, Omega * t) on a LocalSimulator — and the π\pi-pulse is just Circuit().x(0). That is the whole bridge from this analog picture to the digital circuit model. We build Braket circuits properly in Term 2.1; the pulse-level control beneath the gate is Term 4.4.


Exercises

E1 (easy). A qubit has Rabi frequency Ω=2π×50\Omega = 2\pi\times 50 MHz on resonance. How long is a π\pi-pulse (an XX gate)? A π/2\pi/2-pulse?

Solution

tπ=π/Ω=π/(2π×50MHz)=1/(100MHz)=10 nst_\pi = \pi/\Omega = \pi/(2\pi\times 50\,\text{MHz}) = 1/(100\,\text{MHz}) = 10\ \text{ns}. The π/2\pi/2-pulse is half that, tπ/2=π/(2Ω)=5 nst_{\pi/2} = \pi/(2\Omega) = 5\ \text{ns}. (These ns-scale gate times are realistic for superconducting qubits, Term 4.4.)

E2 (easy). On resonance, at what times is the qubit in an equal superposition, P0=P1=12P_0 = P_1 = \tfrac12?

Solution

$P_1 = \sin^2(\Omega t/2) = \tfrac12 \Rightarrow \sin(\Omega t/2) = \tfrac1{\sqrt2} \Rightarrow \Omega t/2 = \pi/4 + k\pi/2,i.e., i.e. t = \frac{\pi}{2\Omega}(2k+1)forinteger for integer k\ge0$. The first is the π/2\pi/2-pulse t=π/(2Ω)t=\pi/(2\Omega), which puts the Bloch vector on the equator.

E3 (medium). Derive the generalized Rabi formula's maximum, P1max=Ω2/(Ω2+Δ2)P_1^{\max} = \Omega^2/(\Omega^2+\Delta^2), and explain geometrically why detuning prevents full inversion.

Solution

P1(t)=Ω2Ω2+Δ2sin2(Ωefft/2)P_1(t) = \frac{\Omega^2}{\Omega^2+\Delta^2}\sin^2(\Omega_{\text{eff}}t/2) is maximized when sin2()=1\sin^2(\cdot)=1, giving P1max=Ω2/(Ω2+Δ2)P_1^{\max} = \Omega^2/(\Omega^2+\Delta^2). Geometrically the Bloch vector precesses about n=(Ω,0,Δ)\vec n=(\Omega,0,\Delta), which makes angle α\alpha with z^\hat z where cosα=Δ/Ωeff\cos\alpha = \Delta/\Omega_{\text{eff}}. Starting at the north pole, the vector sweeps a cone of half-angle α\alpha about n\vec n; its lowest point reaches rz=cos2α=12sin2αr_z = \cos 2\alpha = 1 - 2\sin^2\alpha, and 1rzmin=2P1max1-r_z^{\min} = 2P_1^{\max} recovers the formula. Only Δ=0\Delta=0 (α=π/2\alpha=\pi/2, axis on the equator) lets the cone reach the south pole. ∎

E4 (medium). Show Rx(π)=iXR_x(\pi) = -iX and explain in one sentence why the π\pi-pulse counts as the XX gate despite the prefactor.

Solution

Rx(π)=cosπ2Iisinπ2X=0Ii1X=iXR_x(\pi) = \cos\tfrac\pi2\,I - i\sin\tfrac\pi2\,X = 0\cdot I - i\cdot 1\cdot X = -iX. The i-i is a global phase, which has no observable consequence (Appendix C) — iX-iX and XX produce identical measurement statistics on every state — so the π\pi-pulse is the XX gate physically. ∎

E5 (hard). From H=Ω2XH = \tfrac{\Omega}{2}X and ρ=12(I+rσ)\rho=\tfrac12(I+\vec r\cdot\vec\sigma), derive the Bloch equation r˙=(Ωx^)×r\dot{\vec r} = (\Omega\hat x)\times\vec r and solve it for r(0)=(0,0,1)\vec r(0)=(0,0,1).

Solution

ρ˙=i[H,ρ]=iΩ4[X,rσ]\dot\rho = -i[H,\rho] = -\tfrac{i\Omega}{4}[X, \vec r\cdot\vec\sigma]. Using [X,Y]=2iZ[X,Y]=2iZ, [X,Z]=2iY[X,Z]=-2iY, [X,X]=0[X,X]=0: [X,rσ]=ry[X,Y]+rz[X,Z]=2i(ryZrzY)[X,\vec r\cdot\vec\sigma] = r_y[X,Y]+r_z[X,Z] = 2i(r_y Z - r_z Y). So ρ˙=Ω2(ryZrzY) ⁣ ⁣122=12(r˙σ)\dot\rho = \tfrac{\Omega}{2}(r_y Z - r_z Y)\!\cdot\!\tfrac12\cdot 2 = \tfrac12(\dot{\vec r}\cdot\vec\sigma) forces r˙x=0, r˙y=Ωrz, r˙z=+Ωry\dot r_x = 0,\ \dot r_y = -\Omega r_z,\ \dot r_z = +\Omega r_y, i.e. r˙=Ωx^×r\dot{\vec r} = \Omega\hat x\times\vec r. With r(0)=(0,0,1)\vec r(0)=(0,0,1): rx(t)=0r_x(t)=0, rz(t)=cosΩtr_z(t)=\cos\Omega t, ry(t)=sinΩtr_y(t)=-\sin\Omega t — precession about x^\hat x, and P1=(1rz)/2=sin2(Ωt/2)P_1 = (1-r_z)/2 = \sin^2(\Omega t/2) as required. ∎

E6 (hard). A 2π2\pi-pulse on resonance returns every population to its start (P1(2π/Ω)=0P_1(2\pi/\Omega)=0), yet U(2π/Ω)=IIU(2\pi/\Omega) = -I \neq I. Construct a setting in which this 1-1 is physically observable, and explain why it does not contradict "global phase is unobservable."

Solution

U(2π/Ω)=Rx(2π)=cosπIisinπX=IU(2\pi/\Omega) = R_x(2\pi) = \cos\pi\,I - i\sin\pi\,X = -I. On a single, isolated qubit the 1-1 is a global phase and unobservable. But make it relative: take two qubits in 12(0+1)A0B\tfrac1{\sqrt2}(|0\rangle+|1\rangle)_A|0\rangle_B and apply the 2π2\pi drive to qubit BB conditioned on A=1A=|1\rangle (a controlled rotation). Then the 1A|1\rangle_A branch picks up 1-1 while the 0A|0\rangle_A branch does not, yielding 12(01)A0B=A0B\tfrac1{\sqrt2}(|0\rangle - |1\rangle)_A|0\rangle_B = |-\rangle_A|0\rangle_B — a measurable phase flip on AA, detectable by measuring AA in the XX-basis. The rule "global phase is unobservable" concerns an overall phase on the whole state; once the 1-1 multiplies only part of a superposition it becomes a relative phase, which is physical. (This 4π4\pi-periodicity is the spin-12\tfrac12 signature of SU(2)SU(2), and the mechanism behind controlled-phase gates, Term 2.1.) ∎


Checkpoint

  1. Write the rotating-frame two-level Hamiltonian and name each parameter.
  2. Derive P1(t)P_1(t) for H=Ω2XH=\tfrac{\Omega}{2}X starting from 0|0\rangle. What is the period?
  3. Give the detuned Rabi formula and the effective frequency; what is the maximum P1P_1?
  4. Geometrically, what is a Rabi oscillation on the Bloch sphere?
  5. Which pulse realizes an XX gate, and why is that the central message of this term?
Answers
  1. H=Δ2Z+Ω2XH = \tfrac{\Delta}{2}Z + \tfrac{\Omega}{2}X: Δ=ω0ωd\Delta=\omega_0-\omega_d is the detuning, Ω\Omega is the Rabi frequency (drive amplitude).
  2. U(t)=eiΩ2Xt=cosΩt2IisinΩt2XU(t)=e^{-i\frac{\Omega}{2}Xt} = \cos\frac{\Omega t}{2}I - i\sin\frac{\Omega t}{2}X, so ψ(t)=cosΩt20isinΩt21|\psi(t)\rangle = \cos\frac{\Omega t}{2}|0\rangle - i\sin\frac{\Omega t}{2}|1\rangle and P1(t)=sin2(Ωt/2)P_1(t)=\sin^2(\Omega t/2); period T=2π/ΩT=2\pi/\Omega.
  3. P1(t)=Ω2Ω2+Δ2sin2(Ωefft/2)P_1(t) = \frac{\Omega^2}{\Omega^2+\Delta^2}\sin^2(\Omega_{\text{eff}}t/2) with Ωeff=Ω2+Δ2\Omega_{\text{eff}}=\sqrt{\Omega^2+\Delta^2}; maximum Ω2/(Ω2+Δ2)1\Omega^2/(\Omega^2+\Delta^2) \le 1, equal to 11 only on resonance.
  4. Precession of the Bloch vector about the drive axis n=(Ω,0,Δ)\vec n=(\Omega,0,\Delta) at rate Ωeff\Omega_{\text{eff}} (Bloch equation r˙=n×r\dot{\vec r}=\vec n\times\vec r); on resonance the axis is x^\hat x and the vector swings from the north pole to the south pole and back.
  5. The on-resonance π\pi-pulse (t=π/Ωt=\pi/\Omega): Rx(π)=iXR_x(\pi)=-iX, an XX gate up to global phase. This shows a quantum gate is just time evolution run for a chosen duration — the physics (Term 1) is the gate model (Term 2).

Further Reading

  • [Sak] Sakurai & Napolitano, §5.5 — time-dependent two-level problems, the rotating-wave approximation, and the Rabi formula derived in full (with \hbar explicit).
  • [NC] Nielsen & Chuang, §4.2 and §7.7 — single-qubit rotations Rx,Ry,RzR_x,R_y,R_z, and how physical Hamiltonians realize gates on real hardware.
  • [Gri] Griffiths & Schroeter, §11.1–11.2 — two-level systems driven by a time-dependent perturbation; Rabi flopping.
  • [Pre] Preskill, Ph219, Ch. 3 — qubit dynamics, the Bloch equation, and control.

← Prev: Time Evolution · Up: Term 1 · Next: Single-Qubit Gates (Term 2)

Ready to measure your state?

6 exercises · 10 checkpoint questions

Start the quiz