Ps7 2010d Solution
Ps7 2010d Solution
MIT 2.111/8.411/6.898/18.435
Quantum Information Science I
October 28, 2010
P1: (Continuous time Grover with noise) Consider a continuous time Grover algorithm on n qubits
employing two Hamiltonians, an oracle Hamiltonian, HO = EO |xihx|, and a driving Hamiltonian
HD = ED |ψi hψ|, where |xi is the target state and |ψi is state with all the qubits being |0i + |1i. Here,
we explore the impact of a simple noise model on this algorithm.
(a) Ideally, HD has the same strength as the oracle HO . However, we the strength of HD may
fluctuate.
√ Assuming ED = EO (1 + δ), calculate the probability of being in the target state at
time π 2n /(2EO ).
Answer:
H = HO + HD = EO |xihx| + ED |ψihψ|
⊗n
1 1 X
|ψi = √ (|0i + |1i) =√ |yi
2 2n y
exp(−iHt) = exp(−ia0 EO tI −ia1 EO tZ −ia2 EO tX) = exp(−ia0 EO t)exp(−it0 (a01 EO Z +a02 EO X))
p p p
where a01 = a1 / a21 + a22 , a02 = a2 / a21 + a22 , t0 = a21 + a22 t.
Hence,
√ 2
a01 + a02 2n − 1 2n − 1 2
0 0 2 2
cos (EO t ) = −πδ/4 , sin (EO t ) = 1 − π δ /32 , √ =1− δ
2n 4
1 π2
p(δ) = 1 − δ 2 (2n − 1)( + n ) = 1 − λδ 2
4 2 16
π2
where λ = (2n − 1)( 14 + 2n 16 ).
(b) Suppose we run the algorithm on an imperfect quantum computer in which each qubit experiences
1
a δ that has a random value given by the probability distribution √2πσ exp(−δ 2 /(2σ 2 )).
Calculate how many times one needs to run the algorithm to know what the target state is with
probability 2/3, as a function of σ.
Answer:
The probability of getting the right answer each time is
Z∞
1
p(δ) √ exp(−δ 2 /(2σ 2 )) dx = 1 − λσ 2
2πσ
−∞
Therefore, after running the algorithm N times, the probability of NOT getting the target state
goes down exponentially as (λσ 2 )N . In order to succeed with probability 2/3, N ≥ log(λσ2 ) 31 .
P2: (Compositions of Hamiltonian operations) Consider a physical system with four energy levels
which are addressable, |0i, |1i, |2i, and |3i. You are provided with controls which turn one of two
Hamiltonians Hb , which couples {|0i, |2i}, and Ha , which couples {|2i, |3i} and {|0i, |1i}. Specifically,
0 α 0 0 0 0 β 0
α∗ 0 0 0 0 0 0 0
Ha (α) = Hb (β) = ∗ . (1)
0 0 0 α β 0 0 0
0 0 α∗ 0 0 0 0 0
Note that α and β are both complex numbers. These may be visualized as transitions between some
subset of energy levels, eg in the hyperfine levels of an atomic system:
|0ñ
|3ñ
Ha Hb
|1ñ
|2ñ
Your mission is to perform qubit rotations in the {|0i, |1i} qubit subspace, while leaving {|2i, |3i} alone.
3
(a) Suppose H1 and H2 are Hamiltonians such that tr|H1 | ≤ and tr|H2 | ≤ . Prove that
e−iH1 e−iH2 eiH1 eiH2 = e−iHc + O(3 ), where Hc = i[H1 , H2 ] = i(H1 H2 − H2 H1 ).
Answer:
Up to second order in
1 0 0 0
0 0 0 0
Hc = i[Hb |β=1 , Hb |β=i ]/2 =
0 0 −1 0
0 0 0 0
0 1 0 0
1 0 0 0
Hd = i[Ha |α=i , Hc ] =
0 0 0 −1
0 0 −1 0
Define He = i[Ha , Hc ]
0 −iα 0 0
iα∗ 0 0 0
He = i[Ha , Hc ] =
0 0 0 iα
0 0 −iα∗ 0
Note that here we can take the sum of Ha |α=γ and He |α=iγ because they commute. eitH01 =
eitHa |α=γ /2 eitHe |α=iγ
h
/2
. i
θ
(c) Let Rx (θ) = exp −i 2 (|0ih1| + |1ih0|) be a rotation about the x̂ axis of the {|0i, |1i} qubit sub-
space (it acts as identity on the {|2i, |3i} subspace). Give a sequence of individual Hamiltonian
evolutions, eg U = eit1 Ha |α=1 eit2 Hb |β=i · · · , turning on Ha and Hb sequentially (only one Hamil-
tonian on at a time), with specified values of α, β, and pulse durations, such that U = Rx (θ)
exactly.
Answer:
0 0 1 0 0 0 −i 0
0 0 0 0 0 0 0 0
Hb |β=1 = = X02 , Hb |β=−i = = Y02
1 0 0 0 i 0 0 0
0 0 0 0 0 0 0 0
4
Composing X02 and Y02 we can get any operation on |0i and |2i. For example
eiθ/2 0 0 0
0 1 0 0
eiY02 π/4 eiX02 θ/2 e−iY02 π/4 = eiZ02 θ/2 =
0 0 e−iθ/2 0 = e
iθ/4
Uz (−θ/2)01 ⊕ e−iθ/4 Uz (θ/2)23
0 0 0 1
−iHa |α=−i π/4 1 1 −1 1 1 −1
Ha |α=−i = Y01 ⊕ Y23 , e =√ ⊕√
2 1 1 01 2 1 1 23
eiHa |α=−i π/4 eiZ02 θ/2 e−iHa |α=−i π/4 = eiθ/4 Ux (θ/2)01 ⊕ e−iθ/4 Ux (−θ/2)23 = U1
But
Therefore, the total pulse sequence for applying an X operation on |0i and |1i is
Rx (θ) = eiHa |α=−i π/4 eiHb |β=−i π/4 eiHb |β=1 θ/2 e−iHb |β=−i π/4 e−iHa |α=−i π/4 e−iHa |α=1 θ/4
h i
(d) Do the same for Rz (θ) = exp −i θ2 (|0ih0| − |1ih1|) , such that you now have “pulse sequences”
for performing arbitrary operations on the {|0i, |1i} qubit.
Answer:
Rotation around Z axis can be obtained from rotaion around X axis as
Rz (θ) = eiHa |α=−i π/4 Rx (θ)eiHa |α=−i π/4
= eiθ/4 (eiY01 π/4 Ux (θ)01 e−iY01 π/4 ) ⊕ e−iθ/4 (eiY23 π/4 I23 e−iY23 π/4 )
= eiθ/4 Uz (θ)01 ⊕ e−iθ/4 I23