Lean, informalized

Nielsen & Chuang · Chapter 6, Quantum search algorithms · p. 261

Exercise 6.11 - Continuous time quantum search with several solutions, found by a single measurement

A Hamiltonian that finds one of several solutions.

Hover over a blue or green term in the text, or a node in the graph, to find it in the other. Select it to read its declaration both in Lean and in words.

The Exercise

A search problem runs over NN items, labelled 0,…,N−10, \ldots, N-1, and a set TT of M≥1M \geq 1 of them are solutions. In the continuous time version the items are the computational basis states of a system with NN levels. The system starts in the uniform superposition of all of them and evolves under a Hamiltonian for a time tt; then it is measured in the computational basis. The task is to find a Hamiltonian, which may depend on TT, and a time after which the measurement returns a solution.

A Qudit and Its Basis

A quantum system SS is a physical system whose state space is a finite dimensional complex Hilbert space.✦ A pure state |ψ⟩|\psi\rangle of SS is a unit vector in the state space of SS, ‖ψ‖=1\lVert \psi \rVert = 1.✦

The qudit with NN levels is the quantum system whose state space is ℂN\mathbb{C}^N with its standard inner product.✦ The computational basis {|i⟩}\{|i\rangle\} of the qudit with NN levels, for i=0,…,N−1i = 0, \ldots, N-1, is the orthonormal basis of standard unit vectors, |i⟩|i\rangle having a one in position ii and zeros elsewhere, and the computational basis state |i⟩|i\rangle is the pure state of the qudit with that vector.

For N≥1N \geq 1 the uniform superposition |s⟩|s\rangle is the pure state of the qudit in which every level enters with amplitude 1/N1/\sqrt{N}.✦ For a nonempty set TT of MM levels the superposition over the solutions |β⟩|\beta\rangle is the pure state in which each level of TT enters with amplitude 1/M1/\sqrt{M} and no other level enters: |s⟩=1N∑k=0N−1|k⟩,|β⟩=1M∑x∈T|x⟩.|s\rangle = \frac{1}{\sqrt{N}} \sum_{k=0}^{N-1} |k\rangle, \qquad |\beta\rangle = \frac{1}{\sqrt{M}} \sum_{x \in T} |x\rangle.(1)✦

The Search Hamiltonian

An observable AA on SS is a self adjoint operator on the state space of SS, A†=AA^\dagger = A.✦ The projector Hamiltonian of a pure state |ψ⟩|\psi\rangle of SS is the observable |ψ⟩⟨ψ||\psi\rangle\langle\psi|, the projector onto |ψ⟩|\psi\rangle.✦ The search Hamiltonian of pure states |φ⟩|\varphi\rangle and |ψ⟩|\psi\rangle of SS is the observable |φ⟩⟨φ|+|ψ⟩⟨ψ||\varphi\rangle\langle\varphi| + |\psi\rangle\langle\psi|, the sum of their projector Hamiltonians. For the solutions TT the Hamiltonian is the search Hamiltonian of |β⟩|\beta\rangle and |s⟩|s\rangle: H=|β⟩⟨β|+|s⟩⟨s|.H = |\beta\rangle\langle\beta| + |s\rangle\langle s|.(2)✦

Evolution and Measurement

A unitary UU of SS is a closed system evolution of SS: an operator on the state space of SS with U†U=UU†=1U^\dagger U = U U^\dagger = 1.✦ It carries a pure state |ψ⟩|\psi\rangle of SS to the pure state U|ψ⟩U|\psi\rangle.✦ A Hamiltonian HH on SS is an observable on SS, the energy of the system. The propagator of HH from time t1t_1 to time t2t_2, with Planck’s constant ℏ\hbar, is the unitary U(t1,t2)=exp⁡[−iH(t2−t1)/ℏ],U(t_1, t_2) = \exp\bigl[-iH(t_2 - t_1)/\hbar\bigr],(3)✦ which takes the state at t1t_1 to the state at t2t_2, |ψ(t2)⟩=U(t1,t2)|ψ(t1)⟩|\psi(t_2)\rangle = U(t_1, t_2)|\psi(t_1)\rangle.

A state ρ\rho of SS is a positive operator of unit trace on the state space of SS, and a pure state |ψ⟩|\psi\rangle determines the state |ψ⟩⟨ψ||\psi\rangle\langle\psi|, the projector onto its vector.

A measurement on SS with a finite set of outcomes is a family of measurement operators MiM_i on the state space of SS, one for each outcome ii, satisfying the completeness relation ∑iMi†Mi=1\sum_i M_i^\dagger M_i = 1.✦ The positive operator valued measure (POVM) of this measurement is the family of effects Ei=Mi†MiE_i = M_i^\dagger M_i, one for each outcome ii. The probability of outcome ii when the measurement is performed on a state ρ\rho is the Born rule value p(i)=Tr⁡(MiρMi†)=Tr⁡(Eiρ).p(i) = \operatorname{Tr}\bigl(M_i \rho M_i^\dagger\bigr) = \operatorname{Tr}(E_i \rho).(4)✦

The projective measurement in an orthonormal basis {|bi⟩}\{|b_i\rangle\} of SS is the measurement whose operator for outcome ii is the projector |bi⟩⟨bi||b_i\rangle\langle b_i|.✦ The computational basis measurement of the qudit is the projective measurement in its computational basis, with outcomes 0,…,N−10, \ldots, N-1 and measurement operators |i⟩⟨i||i\rangle\langle i|.✦

The Claim

With ℏ=1\hbar = 1 the propagator of HH from 00 to tt is e−iHte^{-iHt}. The Hamiltonian maps the plane spanned by |s⟩|s\rangle and |β⟩|\beta\rangle into itself, and with the overlap ⟨β|s⟩=M/N\langle\beta|s\rangle = \sqrt{M/N} the evolution of |s⟩|s\rangle is e−iHt|s⟩=e−it(cos⁡(M/Nt)|s⟩−isin⁡(M/Nt)|β⟩).e^{-iHt}|s\rangle = e^{-it}\Bigl(\cos\bigl(\sqrt{M/N}\,t\bigr)|s\rangle - i \sin\bigl(\sqrt{M/N}\,t\bigr)|\beta\rangle\Bigr).(5) At t=π2N/Mt = \frac{\pi}{2}\sqrt{N/M} the cosine vanishes, and the qudit is in |β⟩|\beta\rangle up to a phase.

Theorem 1. Let N≥1N \geq 1, let TT be a nonempty set of MM levels of the qudit with NN levels, and let t=π2N/Mt = \frac{\pi}{2}\sqrt{N/M}. Let |ψ⟩=U(0,t)|s⟩|\psi\rangle = U(0, t)|s\rangle be the uniform superposition evolved by the propagator of the Hamiltonian (2) with ℏ=1\hbar = 1, and let p(x)p(x) be the probability (4) of outcome xx when the computational basis measurement is performed on the state |ψ⟩⟨ψ||\psi\rangle\langle\psi|. Then the outcome lies in TT with probability one: ∑x∈Tp(x)=1.\sum_{x \in T} p(x) = 1.(6)✦

What Is Checked

The claim is stated over the qudit with NN levels, for every N≥1N \geq 1 and every nonempty set TT of its levels, and over two constructions on it: the propagator of the search Hamiltonian HH applied to |s⟩|s\rangle, and the computational basis measurement with its Born rule probabilities.

The library proves that |s⟩|s\rangle and |β⟩|\beta\rangle are unit vectors and that HH, a sum of two projectors, is self adjoint. For the measurement, the projectors |i⟩⟨i||i\rangle\langle i| satisfy the completeness relation, and the effects EiE_i are positive with sum 11.