QGSS'20 Lectures
QGSS'20 Lectures
QGSS'20 Lectures
2020
geoeffneten
Introduction to Quantum
-5--2--2 Computing -0 -
Overview :
-
,
states
.
,
are O or 1
superposition
" " "
•
in ,
a can in i. e. in 0 1
, ,
→
super positions
all ow to
perform calcalations on states at the some time
many
up
to all states
only get superposition)
"
(
→ we can one answer answer in the
it is not THAT
easy to design quantum algorithms ,
but we can use interferenc effects
( other out while the
" "
Diracnotation
used to describe quantum states Led b E? ( 2- dimensionat vector with enties)
. :
a.
E
complex
-
Ket :
las -
I :) irre;"!!!!!
-
bra : < bl =
Ibst :
( Y) ! (bi bi)
-
bra.net : < BIÄÖÄD ten bit azbz "
"
=L albi c- E ↳
complex number)
" "
n
) an X
)
Gbi anbei
-
Ket bra - :
Iaxb , =
( azbi aab;
↳ 2×2 matrix) -
•
we define the states los : =
( I) and 11 ) : =
(9) ,
which are orthogonal : < 0113=1.0+01 :O
(to > TM ) =
Ii)
Measurement
FR.
choose we orthogonal bases to describe & measure quantum states ( →
projects measurement)
during a meas .
on to the basis { los ,
In } ,
the into other
to > or Its →
as those are the eigenstates of Tz ,
we all this a Z -
measurement
e.
there are
infinity many different bases ,
but other commons ones are
{ Its : :
¥ ( los ths)
,
Is :-.
fz (los -
Hs)
} and { Iti ) Jan : :
Host ih ) , 1- is :<
f. (los -
it ) } ,
corresponding to the
eigen
states of a and sry .
respectirely .
: a in meas . on
'
KXIYY
Plx ) =
.
? Plxi) 1 :
example s: 14>
ffn Host VI 11 ) is meas in the basis { 10 ? 11 > }
•
-
: .
:
'
} →
pa) :
}
= =
-
14 > =
¥ ( to > HD -
?
?
→ Ptt) = 1414>1 1¥ KOHN ) -
¥ -
Ho > -
Hs) /
=
O →
expected , as AM > = o
↳ PH =
# Y? -
orthogonal
Bloch sphere :
-
e -
where y E [ 0,2T] desaibes the relative phase and OE [9T] determines the
probability to
all normal ized states be illustrated the surface of here with radius
sp
a
can
pure on
(%) *
, o>
: : : :: : f)
"
-
-
÷
.
.
-0
Ez :(§ )
t
↳ : :
T ö
'
y
:
→
-
X
- It is :
O :
Ez ,
y
:
Ez →
F (Y ) : n>
I i)
- :
O :
Ez .
y
:
Ez → E-
(%)
Beck On the Bloch sphere angles twice big in Hilbert
space
: are as as
. .
e.
g.
los & les are orthogonal ,
but on the Bloch sphere their angle
is 180° .
For a
general state 14 ) : Co Ost . .
.
→ -0 is the
corresponds to a
projectin on Z axis
-
and analogons ly for X& Y
!
I. Quantum Circuits
FELLNER
circuit model :
sequence carry ,
ifeng.rithn.FR -
Singlegebügeltes
•
Classical example : NOT I -
Do -
J
•
quantum example s : as
quantum theory is unitary ,
quantum gates are
represen ed ! by
unitary matrices : UTU = 11 recall
:( I;)
"
Dirac notation : U
u ; = UOOIOXOI tun 10×11
( )
O 1 +4^011×01 11×11
JX 10×11 11×01 um
*
= +
-
=
1 0 Dirac notation
↳ 5×10 > =
. : = =
,
gate ,
e.
g.
N) -
⑦ -
Iss rotation around x
-
axis
T
by
-
Tz =
(G !) =
10×01-11×11
↳ rz.lt > =
(J ! ) # / ! ) # (1) =/ .
: -
7 ,
TH -
s: (10×01-11×11) # .
(los IN )
-
axisbyu
-
sry
:
| ! Ö) : i -
Ox Oz -
bit & phase flip
% , sry & Tz are the so -
called Pauli matrices and a?
-11=1 ! !) (does
nothing)
together with
identity 11 they forma basis of 2×2 matrices
(→ any 1-
qabit rotation can be written as a linear combination of them)
-
Hadamard gate :
one of the most
important gates for quantum circuits
H
:# ( ! ! ) Ii =
( IOXOI +10×111-11×01 -
11×11)
↳ HIO> =
In / ! ! ) / ! ) -
:
# (I) : It >
,
HI 1) =
# ( loxott 10×111-11×01 -
11×111.117 :# 1107 1) =/ - -
)
S H -
is
applied to charge from Z to Y basis
Multipamrtitenqaantamstatesm
•
we use tensor products to describe
state
multiple states : las ① Ibs =
(G) (!!) ⑦ :
| !! !)
•
example :
system A is in 11 ? and
system B is in state 1%
| !)
↳ remark : states of this form are called uncorrdated ,
but there are also bipartite states that cannot
be written as 14h .
⑦ Iy >ß .
Those states are correlated and sometimes even
II | !)
a so
-
called Bell state ,
used for teleportation cryptographg
.
,
Bell tests etc ,
.
IEEE gates
classical
example XOR
J III × ⑦ irreversible
Agieren the output we cannot
• -
→
y
:
re the
cover
input )
BUT :
as
quantum theory is
unitary ,
we
only unsicher
unitary and therefore reversible
gates
quantum example
(0 0)=100×001
• : 1 O O O
( µ,
^
* 101×01 It I 10×1^1 tl ^^ ° !
0 0 0 1
0 0 1 0
f) | ! )
↳
CNOT 100 > : CNOT .
= : too > + CNOT
1103×4 =
111 >xy
xy y ,
infantil
circuit
Y) Iaa
:
0 1 01
10 1 1 ± reversible XOR
11 1 0
on
systems A. B cannot be written as 14h
.
to } ,
it is
entangled
Bell states
e-
There are
for so -
an orthonormal basis :
IT
"
-
¥ (1003-111) ,
14
"
>
:'-( 101 > -
Ho )
write 14J > Morini ) 14
"
>
in
general we can =
gefiel
Bellstatesmpiy.FI#JlyijsIj7s- AB
initial state
ltaoxdrdlijh.rs 14J >
Iij >
n.rs
HA (NOT"
→
1015 ( 101 > t 111>3/5 →
( 101 ) +1103)4147
1103 400 ) -
110DM ( too > -1×1%2=14^0 >
→
opposite direction : Bell measurement
n
-2mm
µ : Alice 10ns : =
dass TßH ? to Bob .
¥ ( 100>aßt 111>AB ) .
initieoftheysteni.IO/7s
# ⑦ 14 >AB
"
=
(al 000>SAß t a 1011 > saß tß 1100> SAß ß 1111 > SAB )
+
=
ITIL ( 1007, t 111)
"
) ( a 10 >ßt ßH >B) ( > satt10 >
101 t
sa ) ⑦( al 1) ßIO ) )
+
( 100 } , -
Hossa ) ④ ( HSB
a -
ßIO )
=
f- [ 14^7*0×10 > B t
140^4*4107,3 )
14 Isa (Tz 107,3 ) ④ Tz ]
"
+ Ox t 14^3 SA
⑤ to > B )
•
Protocol :
①
②
-
FIETE)
""
1 . Alice performsameas .
.
2 She sends her Classical
- . . . .
. . .
14003 _
.
. _ . . . . . .
.
. . . . . . . .
.
. . . . . . . . . . . . _ . . . .
.
.
5×0 - -
3. Bob
apples riss to his
'
⑦
14^07 Gz Icf>ß 10 Tz u
14 " S 5×410/3,3 11 Tz ?
"
to >
pole
"
• "
: -
re
,
-
the color indicates the relative phase compared to state to>
¥ to > NEIN
ÖÖÖ: ④
000
states
"
•
for n
qubits ,
there are 2 basis ,
e.
g. for n =3 we have 001 .
010
.
100
011 , 101, 110
those basis
we
plot states as
eqaally distributed point on a
sphere ,
^ ^^
"
Ü
"
"
the
"
pole aligned
on
pole ,
1 on
parallel st the each latitade anstand and increasing from North South
"
of 1 's to
'
on s
,
. number on is
Eagle for y :3 :
→
size & color 0
f the blobs as before
000>
vii. no )
es '
Fill ""
:
""> +
*.
i.
gßE-_NT]
Eigentum!gerithmslf
21.07.2020
Oranien
I) Deutsch Jozsa -
I. Deutsch Jozsa
algorithmentheorie -
nennen =n .
Orontes .
goal determine
:
some property of the oracle
asing the minimal number of que
ries
"
•
on a classical computer ,
such an Oracle is
give by a fct f :{0,13
.
→
{0,15
bit Oracle
IIII!! }
be
unitary which performs
0g : can seea as a
J
,
I the
µ map Oglxsly > =
Ixslytoflx) >
→
for f : { 0,13
"
→ { 0,13 we can contract Ug :
,
Fängt "
"
ancilla-
qubit
In
{ Jz
WHO ) Hs) Ixsly > if fix) :O
# ( Ix )
-
=
IM UN los)
Ixslg ) ifflx) 1
-
= =
-
,
f)
"
=
1×31 )
y
Ix >
Hadamard
Tur m
on in
qubits
_
.
-
recall that It to> =
Its =
# Hatte>) ,
HH ) : ts
¥40 > -
)
k
f- Ä To > + HI
"
#
"
113)
for x E { 0,1} : Ixs -
THI -
Igs
=
( to > t
f- e) 4) =
¥ -
""
=
¥ Ego "
EN l k>
.
inner product
"
for '
HÄM
HD
{
{ 0,13 Ixo>
!!)
µ
x c- : -
"
↳
"
"→
:
's - " " =
# E.„
lxn.is -
TH -
↳ ever
ylyis is eitler It > or ↳
Igs must be a
superposition of all possible
1×3=101 > bit
"
e.
g .
: 2
strings
J lys
'
% :
It> ⑤ I -
> =
Elias -
101 > tho> -
Hes)
Fur
Deutsch Jozsa algorithmen
-
g- -
function f :{0,1J
"
We
give a { 0,1J realized by oracle of which
→
are →
,
an ,
# O '
is
'
Classical :
we need to ask the oracle at least twice but if we
get twice the ,
some
Output ,
we need to ask again .
. .
.
at most It 1=2^71
"
→
queries ,
n : # input bits ,
N 2
> : # realizable bit strings
"
demonstrative different to Throw the fair
example : 2 ways a coin → is coin ?
quantumsolution.im needs
only one
query !
"
÷:
"
"
° 1
2
14 )
"
•
14 > = 100 07 = 10 ) oo
;
- o
Mi
. .
.
X. Yo
IY > =3 ¥!
"
H '"
>
? ""
e =
.
",
n
"
.
Elf ! ! }! ÄK
"
¥
""
# E.
" "
•
Mir H -
-
Mas =
H -
H Ixs = >
"„
"
"" "
=
! [ II.go.y.nl 1) -
""
-
= :
Ck
OIÄ .EE/klf=koo....l2
"
Plug oo . . . oo . .
-014371 ! !! nen
1 h:
{
= ,
00
. . .
0
'
.gg !!!
orthogonal)
=/ In EN) ! {
1
O
.
-
.
if f balanced
{
"
=
+2
"
, if f) EO
2
if flxl Et
-
°
,
if flxl balanced PTT
II-rbafgorith.mn
g- U n - Umm
rother : find
elements
" " ×
.
t.gl#=1
→
classical
alg . needs on
average Nz = OCN) time
folx) :{1Qifdse
"
goal :
,
gieren an
Ug with
f :
{ 0,13 → {0,13 , ,
F , ,
Ug
105%105*1×3
Ich → Iws
Ug
: -
:
.
GETEILT :
:
↳ Ug : 1- 2 IWXWI ↳
Ugo 210×010 =
"
-
11
DH : It
:*
107 -
=
."
:
:
:-.
Ist
repeat r times
Clair y
= w ( with high prob ) .
:# ¥,
"
Ugo .
H = .
H -
Hot! H "
= 21 SXSI - 11
"
Get I be the plane spanned by Iss and Iws and Kt lat ) be the state
orthogonal to ↳ in E :
Int) : =
¥77,1 ×>
VII Iwt
ftp./ws--:cosEIwtstsinEzIwsdefineO-s.t.sinzE:Izn-
Is > = > t
protocol :
0=2
arcsinfznn
- -
Ä
tun [ " s>
2.) Apply Ug -11 -
21WXwl →
reflection at lat>
±:):::c: : :.ie?::::::.
" "
after r
applications of 2.) & 3.) the state , is
oH:→ Ups>
rotatedbyr.O-nea.ge
E EVT ( ru)
↳
chooser.si .
rot ! E r
-
-
¥ E -
ü
=
I -
probability
plus ) > 1- sin ! =
1- ¥ µ
:
'Ä µ. ) n
Amplitude amplification
nun zum -
c-
1.) Is > : = H
⑦ "
10)
"
2.) Ug Is > = (1- 21wxwl) Is ) →
flip amplitude of ↳
amplitude amplitude
÷ -
I-E.in
EF I :p i
-
: :
-
-
.
¥: *
I-E.it
:
;] ¥ ;
-
:
-
-
÷ ü:{i :
t.vn
± -
? Iis
yidds :
|
as
Ip (
di
3- E) •
|
KIND a) 145.2 E ! lj ) 7. <ht ? a- Ii)
? li>
¥:-, ¥
-
a:
- - . -
""
"
- - - -
g.
-
III. ¥ =
,
.
=3.se is
=
ca>
SHIN
-
Iws Er ?! Iwie →
tut > IFF ¥ '
.. .. . . "
' "
Iss =
↳
sin ¥ VII =
angle becomes
lager !
r =
4.at#rmjy-E--O(VN-T)
→ we can see this
speed up also when
looking at the amplitude :
1.) 2.)
⇐
TÄTE Ä "
" " ":
§ §
lo ) h ) Iwie Iwz) IN-19
"
IÄÄ
¥
in
µ
-
;
goes faster
"
←
:÷i
ShorbAlgorithmiandQPEd
Lecture3
Problem
given a function that is periodic find its period period
for la k kf
strict definition f fly sety iff
se
y
np
21T
easy example
hard example 5
µwh g
classically 0 exp 43
c n 13µg a
P has n bits
quantum Shor's algorithm 0 n'Clog n loglog n
little faster than 0hr3
Reason why this works 1 quantum Fouriertransform
2 modular exponentiation
Implication difficulty of factoring on a classical computer is basis
for security more tomorrow
f m
lot ID
ta
lo 107
105 IT
117 117
IT QFT be
p t
Fourierbasis computational basis
N I
zitisey
analogous to inverse discrete
Z ENT y Fourier transform
y o
eg 1 qubit case N 2
To
ta My ta IZ Y
lo HD
it in El n
lo C D 117
107 117
eatin8 231y y
IIe 11 of
gy o eg ly I't
4 22 l t 2 lo t 20 o
4
23 2417
12403 12
0323
N I Y Cycyrys Yun
Tse 210ft 2n y
tip foe
12 Yu
1
y
oe.TT
imi hiyn
shy
e4 IYu
Yuya ynk
fro e
i i
in
eating
go g
TIE
k 4 term to my
62 term to t Iti F ly
notice i D foe Nz Az An
ye 0 MHz kn l
2
OOH Uz Rn 2
Hq
Final
form't lost e EH lo e h
be 157
QFT la lost e EH lo e Em
eatiEhs
do
Notice 2mi E 2Ti
E tf
21T 21
z
since e2TiZ
1 forany integer Z
Tilt
e2Ti I 2
2
eati eti I
in eminently
Ez eatin'Enly
0
9 n k
Let's write y Yk2 y Yi Yn
zn y 2 yw
u
then y Yh 2k and we have
gn
n
ta 1µg
eat Emily ya yn
Yeo
n 21T t 12h
If II ly.y
e yn
remember
o j
to
127 1 to I Tn lost b
if
lo EI In
Is e E E e
e too to Hints
phase is qubit dependent
e2TiH2t l lo oo needto add up more components with
1more 1
G UROTulsej e lap
WROTE to 2k
211424
f
applies phase e
THI DURAI ITURot URoTn
i
tii
to e'Time seipp I
steps
i
step 3 lo t e e me Mµy
n
im
c
step n lost e e ly
1
exp
Ii
Fat KINI
t
My
recall se 2 se 2 2kt t I nu
t
age Zen
term above is 21T
exp if
after step ni lost e 711 being sen
reverted
order fromderivedform but ok otherwise
Part II Quantum Phase Estimation
U ly ei04147
Can we extract 04 given the ability to prepare 147
107 i Dt t.TT
F
c
c q l
I
THX Duet 1
c
o I
2 3
step o to 147
step 1 107147 117147
t
step 2 lo 147 t 117 e04147
04 Go ID 143
step 3 f 107 1 117 47 t e
I to Hei 04 t H l e
4
147
measure qubit O
if 04 0 prob o cos z4
prob I sin Et
To measure precisely using this bit of
information need to do
lots of measurements
155
04 10 Pablo Prob D 09999,7 615
eg
04 100 n 0.9924 O 007596
TF
i
he I ITEEFI 7u
2
step 0 1050 147
0
step I lo t 1177 147
GET
Before proceeding note 212 147 UE U 147
wite 04147
42
e 147
I E
step z lo tei I co Eat Ip
2g
D t e y
Us app
n
127 lo I In lost b
iz
lo e In
he
iF
T
TTT
ITEEFI Fug Tue
24
Measurement should yield
10
Lecture 4 shores Algorithm II Factoring to period finding
Abe Asfaw
Quickreviey
Yesterday covered QFT and QPE
G QFT basis change from computational basisto Fourier basis
quick demo
circuit implementation
1
wrote
Expezia
o
target
fine control
CROI 11 x e l1xj7
Eder ofqubits Hln 7
lo7tg
Form of QFT n qubits NE 2 Ise 124 sen
QFT Ise L
in
lo te Tp
x oxz sen
lo t e 117 se ooses and
o Iii sea Coo osei
212thse
lo t e 117 toe
from H gate
fromoUR0Tn fromUROTz
Circuit implements o o
Recall se Ge uz seu
2 x 2 Zzz1 Inn
Ise lost exp 22in a 117
frompUROTny four H gate
9
lazy lost expf sent
fifteen it 12 sea ID
to
define se 0K22g Rn
2n o 12 2221 Isen
then above is se and we have
j
d
Iaa lost exp 2I 117
i
z In
se
Isen to exp
000 Own
2 Quantum Phase Estimation
QPE protocol
E'T 0147 2h0 where
given UI QPE gives
n is the number of qubits used to estimate O
n
youI lo'T
T.i
1 THE 1
EEEnys.no
U lol
4
Etiology to this point
up
have QFT in input
subtleties i notice the zito register
in the QPE protocol
b 3 quotient 1 355
3
remainder 2 2
5 2 mod 3
x I 2 3 4 5 6 7 8 9
2e I 2 0 I 2 0 I 2 0 mod 3
generally
y mod _3kt
y for some KEZ
2e 3 se
a Ark mod N
sett 5 0 mod 15
4 I att 3 5
p g
4 Quantum circuit for factoring 15 subtleties later for now
I I I follow aloyg
g
1M l
I t
t
DH QTT l I
i
1 I e I i
1 c l i l
w p
1
10504 T f T i
To l iz f3 4 s 5
be Iw be lwtofa.CH
04107
04
fan Cx a Cmod N
step 0 107
iz z
04
step 1 H 410504 103
O 2 Z we have
Ya 11741137g t 1274142 t 1374177g
I 102411 t
Recall N I 2
ay
1
QFT Ise 127
if 2 e ly
g op 2
QFT Ee In
M
ly
foe
2
QFT 132 Eno e 34
y
38
IT y0
e ly
H
QFT'T 7 a
I IIE ly
Y
QFT 111 a Ezo e ly
5
t
QFT 115 a ta of ly
step 5 measure the foe register get 107 on 147 or 187 or 1127
with equalprobability of 14
Analyze measured result i meas results peak near
j Nyforsomejez
10 is trivial If we measure 103 restart X
147 j kn 4 r 4 if j L r 4 is even good
42 2
mod 15 4
see a mod N 13
5
set I j g od C sett N 5
se I 3 j g od se l N 3
18 8 r 2 and j I 4 and
j 16g or r
j 2
d
works like above
131 mod 15 2
set 1 3 3
ged Catt
N
se i I
J
1 partial
god se l N
ark 11 to modN
How to implement
Ufa n
fa n se A mod N
be 12 best 25hm
q2 mod N
22 a2
af Ma2 mod N
se ooo
Jez
Nu
n 2
m
91 at a'and
o o
depth is polynomial
in n
9
can see Shor's algo is QPE in disguise
– We need to somehow make perfect logical qubits out of noisy physical qubits
• Repetition code
• Surface code
• Syndrome measurements
• Decoding
• Logical operations
– For d repetitions
𝑑/2 𝑑/2
𝑑 𝑛 𝑝
𝑃= 𝑝 (1 − 𝑝)𝑑−𝑛 ~
𝑛 (1 − 𝑝)
𝑛=0
– The example we just considered was one of communication, with errors occurring during transmission.
– The ‘noisy operation’ here is just doing nothing (with some errors)
– Here we don’t completely decode and re-encode, but just do it enough to find and fix the errors.
– This method works great for bits, but not for qubits
– We do need to measure, to get information about errors. But we must avoid learning
about the encoded information
– For this, note that we don’t actually need the bit values. We only need to know which
ones have a different value to the rest.
– They are are initialized in state |0ۧ, and are used as the target for two CX gates
– The net effect is to measure the observable 𝑍𝑗 𝑍𝑗+1 : the Z basis parity of the two qubits
cx 00ۧ = 00ۧ
cx 01ۧ = 00ۧ
cx 10ۧ = 11ۧ
cx 11ۧ = 10ۧ
– The results are used to identify where errors likely occurred and how to remove their effects
– Errors create pairs of ‘defects’. Majority voting can be used to find a minimal pairing.
– Bit flips create defects with space-like separation, measurement errors with time-like
– A likely set of errors corresponds to the ‘minimum weight perfect matching’ of the graph
– The code only allows us to detect and correct bit fips, and
only bit flips
– Example 1: z basis states are product states, the others are entangled
|0ۧ → |000ۧ
1 1
|+ۧ = (|0ۧ + 1ۧ → (|000ۧ + |111ۧ)
2 2
– Example 2: Distinguishing encoded z basis requires a single qubit measurement, distinguishing x basis
states requires d
– Example 3: Flipping between z basis states requires d gates, but the x basis only takes one
– It is not good when things are too easy, because they are easy for errors too!
– In the surface code we use a 2D lattice of code qubits, and define observables for plaquettes and vertices
– These are similar to the two qubit measurements in the repetition code
– Valid states are those with trivial outcome for all plaquette syndrome measurements:
Even parity on all plaquettes
– Valid states are those with trivial outcome for all plaquette syndrome measurements:
Even parity on all plaquettes
– Valid states are those with trivial outcome for all plaquette syndrome measurements:
Even parity on all plaquettes
– These can’t be a different encoded state: they are only a few bit flips away from our encoded 0 state
• Pick a vertex
• Any line from top to bottom (passing along edges) has even parity
– Unlike the repetition code, distinguishing encoded 0 and 1 requires some effort (which is good!)
– Uses the edges has a nice advantage: we can think of them as large (unenforced) plaquettes
– In this case they look at the |+ۧ and |−ۧ states, and count the parity of the
number of |−ۧs
– They are encoded using suitable patterns of |+ۧ and |−ۧ states
for the code qubits
– Satisfies 𝐴𝑣 |𝜓ۧ= |𝜓ۧ and 𝐵𝑣 |𝜓ۧ= |𝜓ۧ, so will give the 0 outcome for all stabilizer measurements
– Though such things can be hard to make, we create and protect them with the syndrome measurements
𝐸 =𝑎𝐼+𝑏𝑋+𝑐𝑌+𝑑𝑍
– This creates a superposition of different types of error on the surface code
– Again we start with the simple and unrealistic case: errors only
between measurement
– Correcting according to the wrong class causes an operation on the encoded qubit (without our knowing)
– What is the probability of such an error, P, given the probability on the qubits of the code, p?
– Who will be among the first to run and study this code?
– Maybe you!
Zlatko K. Minev
IBM Quantum
IBM T.J. Watson Research Center, Yorktown Heights, NY
Image copyright:
@zlatko_minev zlatko-minev.com ZKM unless otherwise noted
What is a real qubit?
Why?
Illustration: IBM Qiskit Textbook Zlatko Minev — Qiskit Global Summer School 2020 (2)
On the road ahead
Qubit in the cloud cQED: Transmon qubit Unveiling the
quantum oscillator
Qubit from atom / oscillator Classical circuits & the LC Transmon qubit
Don't need to know much going in, but we will go far Advanced
material
Examples: simplest, most practical examples
Step by step
Ask questions!
Icons here by FreePik. Dangerous bend symbol based on Donald Knuth. Zlatko Minev — Qiskit Global Summer School 2020 (4)
Avoid firehose of information
|1⟩
|0⟩
* Images: Minev, arXiv:1902.10355; Illustration on right: IBM Qiskit Textbook; Zlatko Minev — Qiskit Global Summer School 2020 (7)
Qubit: idea and reality
|1⟩
|0⟩
* Images: Minev, arXiv:1902.10355; Illustration on right: IBM, Qiskit Textbook Zlatko Minev — Qiskit Global Summer School 2020 (8)
Qubit: idea and reality
|1⟩
|0⟩
* Images: Minev, arXiv:1902.10355; Illustration on right: IBM, Qiskit Textbook Zlatko Minev — Qiskit Global Summer School 2020 (9)
Operation at
15 mK (-273.13 °C)
Zlatko Minev — Qiskit Global Summer School 2020 (10)
qubit
qubit
readout
qubit
readout
qubit
readout qubit
readout
…
readout
…
* Laptop image: rawpixel.com; Photos: IBM Zlatko Minev — Qiskit Global Summer School 2020 (13)
Qubit
Concepts
|1⟩
|0⟩
* Bloch sphere is a mere geometrical representation of SO(3), but the density matrix ρ is in SU(2), a double cover of SO(3).
* A density matrix operator lives not in the Hilbert space H but in the Liouville space H ⊗ H.
Images: Minev, arXiv:1902.10355; atom art: Indoleces. Zlatko Minev — Qiskit Global Summer School 2020 (15)
Origins of quantum
e.g., atom cloud
individual
atom
…
|3⟩
Helium
|2⟩
Neon
|1⟩
Sodium
|0⟩
Mercury
Quantization of energy
* The notion of an energy level was proposed by Bohr in 1913.
Image credit: NMSU, N. Vogt
Zlatko Minev — Qiskit Global Summer School 2020 (21)
The light of atoms
Quantization of
energy
Low-loss
No level here
Level diagram allows for
qubit-specific control and
readout
capacitor inductor
Magnetic flux
Radius
Zlatko Minev — Qiskit Global Summer School 2020 (26)
non-linear
inductor
anharmonicity
Magnetic flux
Radius
Zlatko Minev — Qiskit Global Summer School 2020 (27)
Big-picture connections
|1⟩
|1⟩
|0⟩ |0⟩
Combine …
Superconductivity
Nominally zero intrinsic dissipation and heat
Nominal temperature far below energy level splitting
Non-linear, robust Josephson tunnel junction effect
* Low energy dynamics. Microwave image: vectorpocket; Nobel photograph: Jonathunder; Josephson photograph: Nobel Foundation archive.
Wavelength (m)
Wavelength
scale
Frequency (Hz)
Effective
temperature
Introduction
Zlatko Minev — Qiskit Global Summer School 2020 (35)
To do quantum,
Classical oscillator
Qubit inductor
Non-Linear
e.g., Josephson Junction
Voltage
Electric field
Capacitance
Initial conditions
Inductance
Initial conditions
Magnetic flux
For kinetic inductors,
~98% of qubit inductive energy is not
in stored magnetic fields, but in
kinetic inductance
Zlatko Minev — Qiskit Global Summer School 2020 (43)
Power and energy
Universal
Constitutive
elemental
relationships
inductance
value
inductor
capacitor
ground
As we will see later, for the Lagrangian description in flux basis, KVL acts as a set of holonomic
constraints and KCL as the equations of motion
Zlatko Minev — Qiskit Global Summer School 2020 (48)
Oscillator analogy
Velocity
Equilibrium position
Resonance
frequency Spring: Svjo
Zlatko Minev — Qiskit Global Summer School 2020 (49)
Lagrangian and Hamiltonian*
Euler-Lagrange equations, F = ma
Henri Poincaré
with
Dangerous bend
ahead on quantization
These look
a lot alike!
Derivation: Dirac derives the quantum form of the Poisson bracket — the commutator — from
merely assuming that
1. Classical Poisson bracket rules hold (by analogy, the new theory must be consistent with the old!)
2. The dynamical variables do not commute; i.e., xp ≠ px
3. The Poisson bracket has a single result and single unique meaning
Problems:
• Operator ordering ambiguities.
Consider A and B are polynomials in x and p; e.g., x2×p or p×x2
• Curvilinear coordinate systems (potentially transmon if cos(Φ/φ0) considered wrapped)
• Quantum gravity …
Zlatko Minev — Qiskit Global Summer School 2020 (59)
Photo: F. Schmutzer
Classical Quantum
Hamiltonian
Phase space
(Heisenberg
picture)
Energy
Classically forbidden
region
Energy /
Scaled
wavefunction
amplitude
Classically forbidden
region
Energy /
Scaled
wavefunction
amplitude
Classically forbidden
region
be imaginary?
Will discuss
answers on my
blog sometime
soon
Zlatko Minev — Qiskit Global Summer School 2020 (74)
Quantum Harmonic Oscillator Applets
Wigner function
https://demonstrations.wolfram.com/WignerFunctionOfHarmonicOscillator/
Coherent states
https://demonstrations.wolfram.com/CoherentStatesOfTheHarmonicOscillator/
Graph: Minev; Spring-mass: Svjo Zlatko Minev — Qiskit Global Summer School 2020 (77)
The LC quantum harmonic oscillator
…
wavefunction amplitude
Energy / Scaled
Classically
forbidden region
Zlatko Minev — Qiskit Global Summer School 2020 (78)
The road behind and ahead
Qubit in the cloud cQED: Transmon qubit Unveiling the
quantum oscillator
Qubit from atom / oscillator Classical circuits & the LC Transmon qubit
Tightly integrated lab work with Dr. Nick Bronn and Co.!
More depth on qubit control
Run experiments on real devices
Thank you!
Zlatko K. Minev
@zlatko_minev zlatko-minev.com IBM Quantum
Superconducting Qubits II
Transmon & Measurements
Introduction
Zlatko K. Minev
IBM Quantum
IBM T.J. Watson Research Center, Yorktown Heights, NY
Image copyright:
@zlatko_minev zlatko-minev.com ZKM unless otherwise noted
The road ahead
Transmon qubit Qubit and readout (/qubit) Lab
Magnetic flux
Energy
Magnetic flux
…
wavefunction amplitude
Energy / Scaled
Classically
forbidden region
Graph: Minev; Spring: Svjo Zlatko Minev — Qiskit Global Summer School 2020 (9)
Phase-space (Husimi Q function)
Josephson
tunnel junction
Non-linear inductor
generalized
magnetic
flux
EJ
Energy
Classically
forbidden region
-EJ
-" 0 "
≈ 3.3 × 10−16 Wb Reduced magnetic flux Φ/!0
Image: Minev et al., EPR to appear (2020) Zlatko Minev — Qiskit Global Summer School 2020 (15)
Semi-classical intuition: phase space picture
EJ
Energy
Classically
forbidden region
-EJ
-" 0 "
Reduced magnetic flux Φ/!0
Zlatko Minev — Qiskit Global Summer School 2020 (17)
Hand-written notes
Hand-written notes
The Transmon qubit
EJ
RWA or 1st
order PT
Energy
Due to nl & commutator Due to non-linearity (nl) Classically
forbidden region
where -EJ
-" 0 "
“Lamb shift” Anharmonicity Zero-point Reduced magnetic flux Φ/!0
due to ZPF fluctuations Zlatko Minev — Qiskit Global Summer School 2020 (20)
The Kerr Hamiltonian of the transmon
EJ
Energy
Solution known (SHO) New from nonlinearity
Classically
forbidden region
-EJ
-" 0 "
Reduced magnetic flux (Φ/!0)
To first order perturbation theory the eigenstates do not change! Only the energy changes. Dispersive.
Zlatko Minev — Qiskit Global Summer School 2020 (21)
Exploring a real transmon qubit
EJ
Energy
Solution known (SHO) New from nonlinearity
Classically
forbidden region
-EJ
-" 0 "
Experimental parameters Reduced magnetic flux (Φ/!0)
Energy
Dispersive, states didn’t change to 1st order
“Lamb shift” due to ZPF
Classically
forbidden region
-EJ
Parameters used in figure (of a measured qubit) -" 0 "
Reduced magnetic flux (Φ/!0)
EJ
Energy
-EJ
-" 0 "
Reduced magnetic flux (Φ/!0)
Zlatko Minev — Qiskit Global Summer School 2020 (24)
Energy diagram and transition spectrum
Energy levels Transition spectrum
ωq
Energy
0
0
Φ/!0
Energy
Fock number Qubit Pauli Z Qubit Pauli X and Y
operator operator operators
0
0
Φ/!0
Minev et al., EPR to appear (2020); Dissertation arXiv: 1902.10355 Zlatko Minev — Qiskit Global Summer School 2020 (30)
Energy participation-ratio: a bridge
Control and beyond Zlatko Minev — Qiskit Global Summer School 2020 (33)
Single-qubit quantum gates
IBM Quantum
Zlatko Minev — Qiskit Global Summer School 2020 (34)
Qubit control
input-output line
qubit (control, shaped signals, but also environment)
Zlatko Minev — Qiskit Global Summer School 2020 (35)
Qubit control: overview
qubit
qubit
input-output line
(control, shaped signals, but also environment)
Zlatko Minev — Qiskit Global Summer School 2020 (39)
Control, noise, and dissipation go hand-in-hand
Can lead to uncontrolled, random bit and
phase flip.
T1 : energy
T2 : (coherence) transverse
Noisy environment, and
always zero-point quantum fluctuations
* This is a major topic in condensed matter physics;
we will only touch on it. Zlatko Minev — Qiskit Global Summer School 2020 (40)
Coherence in superconducting circuits
Isolation vs. control
Fluxonia
(Chicago, UMD)
(Yale) (Yale)
2-D Transmon
(IBM)
2-D Transmon
(IBM) multi-qubit
systems
Year
2020
Image reproduced Reagor (2015), an update of Devoret and Schoelkopf (2013), and updated
Zlatko Minev — Qiskit Global Summer School 2020 (41)
Is the junction phase/flux
compact or not?
Barrel image:freepik.com Zlatko Minev — Qiskit Global Summer School 2020 (47)
Two basic classes of measurements
Demolition Non-demolition
Caution: Classical
measurements can
also be disturbing!
Minev et al., EPR to appear (2020); Dissertation arXiv: 1902.10355 Zlatko Minev — Qiskit Global Summer School 2020 (59)
Energy-participation eigenmode approach
Linearize
Minev et al., EPR to appear (2020); Dissertation arXiv: 1902.10355 Zlatko Minev — Qiskit Global Summer School 2020 (60)
Linearize
All linear
Easy to solve, classical problem
Eigenmode solutions
Minev et al., EPR to appear (2020); Dissertation arXiv: 1902.10355 Zlatko Minev — Qiskit Global Summer School 2020 (61)
Linearize
Solutions are eigenmodes (normal mode),
just more SHO. Their Eigenfrequencies are
“Bare”
coordinates
Canonical transformation
eigenvectors (found classically)
Use in It follows
See Minev et al., EPR to appear (2020); Dissertation arXiv: 1902.10355 Zlatko Minev — Qiskit Global Summer School 2020 (65)
Breakdown of the Hamiltonian
Transition
spectrum !q !q χqc
Recall that harmonic
oscillator has evenly spaced
levels
ωq -Δq ωc
Lorentzian peak
(response of cavity on irradiation)
distinguishability
Tightly integrated lab work with Dr. Nick Bronn and Co.!
Run experiments on real devices
Thank you!
Zlatko K. Minev
@zlatko_minev zlatko-minev.com IBM Quantum
The important thing is not to stop questioning.
Curiosity has its own reason for existence.
Antonio Mezzacapo
Qiskit Global Summer School 2020
IBM Quantum
Molecular Hamiltonians
Molecular Hamiltonians IBM Quantum
Quantum
Speedup
Hardware Practical
economy applications
The parity set: tells us the set of qubits that encodes the parity of the fermionic modes
with index less than i. The parity of this set of qubits will give us the global phase
The update set: the set of qubits that must be flipped when the fermionic mode i change
occupation
The flip set: the set of qubits that determines whether qubit i has the same or inverted
parity with respect to fermionic mode i
Fermion to qubit mappings IBM Quantum
H2 Hamiltonian (4 fermionic modes): 4–qubit to 1–qubit
Independent Pauli 3 single-qubit Pauli
Symmetries: with commutation rules
Variational circuits
Variational circuits IBM Quantum
Variational circuits IBM Quantum
Variational circuits IBM Quantum
Fermionic Hamiltonian problem
Energy measurement
on quantum chip
Solution
Variational circuits IBM Quantum
Variational circuits IBM Quantum
Variational circuits IBM Quantum
Variational circuits IBM Quantum
0
10
Entanglers: 0
Entanglers: 1
Entanglers: 2
Entanglers: 3
-1 Entanglers: 4
10 Entanglers: 6
Entanglers: 8
Energy Error
10 -2
Chemical accuracy
-3
10
10 -4
10 -6 10 -5 10 -4 10 -3 10 -2
Noise Strength
Variational circuits IBM Quantum
QAOA on fixed
Qubit architectures
UCCSD
Hardware Efficient:
QAOA
Entanglers + Euler Rotations
Optimization
Optimization IBM Quantum
Fermionic Hamiltonian problem
Energy measurement
on quantum chip
Solution
Optimization IBM Quantum
Abhinav Kandala
akandala@us.ibm.com
IBM Quantum
0.2
Energy (Hartree)
-0.2
-0.4
-0.6
-0.8
-1
-1.2
0 1 2 3 4
Interatomic distance (Angstrom)
Energy (Hartree) 0
-0.2
-0.6
-0.8
-1
∆E Rate ∝ Exp(-∆E/kT)
-1.2
0 1 2 3 4
4 © 2017 IBM Corporation
Interatomic distance (Angstrom)
IBM Quantum
N
X XN XM
1 2 ZA X 1
He = ri +
i=1
2 i=1
r iA j>i
r ij
A=1
H| Gi = EG | Gi
10-9 m 10-3 m
Classical Quantum
0 0
Quantum Bits (“Qubits”) can be in
a superposition of 0 and 1
1 1
10110…1001
10110…1001 11010…1010
00100…1101
10110…1111
N bit number: 2N possible states Can be in a superposition of all
Always in 1 of those states 2N states
IBM Quantum
Quantum Circuits
Decompose unitary N
U 2 U(2 )
| i = U |0 . . . 0i
†
U= CNOT0,2 T0 T2 . . . X0 X1
Universal gate set: CNOT (two-qubit gate) and single qubit rotations
David DiVincenzo, Two-Bit Gates are Universal for Quantum Computation, PRA 15, 1015 (1995)
IBM Quantum
Short depth circuits
100 µm
Superconductor (Al)
Insulator (Al2O3)
Superconductor (Al)
Superconductor (Al)
|2⟩
Energy
|2⟩
Energy
|1⟩
ℏ𝜔0
|0⟩
|2⟩
Energy
|1⟩
ℏ𝜔0
|0⟩
Q2
𝜔2
Q1
Device characterization
Frequency Single-qubit U2
Qubit T1 (µs) T2 (µs) (GHz) Readout error error rate CNOT error rate
cx0_1: 1.380e-2,
Q0 64.2 74.54 5.286 2.45E-02 6.49E-04 cx0_2: 1.430e-2
cx1_0: 1.380e-2,
Q1 77.35 83.17 5.238 1.75E-02 5.72E-04 cx1_2: 1.713e-2
cx2_0: 1.430e-2,
cx2_1: 1.713e-2,
cx2_3: 1.470e-2,
Q2 50.9 41.91 5.031 2.75E-02 8.32E-04 cx2_4: 1.708e-2
cx3_2: 1.470e-2,
Q3 60.04 47.72 5.296 3.35E-02 9.79E-04 cx3_4: 1.245e-2
cx4_2: 1.708e-2,
Q4 57.19 64.79 5.084 1.45E-02 4.17E-04 cx4_3: 1.245e-2
Device characterization
Frequency Single-qubit U2
Qubit T1 (µs) T2 (µs) (GHz) Readout error error rate CNOT error rate
cx0_1: 1.380e-2,
Q0 64.2 74.54 5.286 2.45E-02 6.49E-04 cx0_2: 1.430e-2
• Incoherent
• Leakage
• Coherent
• Measurement
Incoherent errors: Energy relaxation
• Initialize qubit in |1> state and probe its
decay to |0>
• P(|1>) = exp{-t/T1}
• Characteristic time: T2
• Leakage out of
computational subspace
Gaussian Pulse:
-1
Signal (a.u.)
-1.2
-1.4
-1.6
-1.8
0 5 10 15 20
Number of Repetitions
Prepare ‘0’
Measured {‘0’ : 99% ‘1’: 1%}
Cavity
Prepare ‘1’
Readout Tone Frequency Measured {‘0’ : 5% ‘1’: 95%}
N
X XN XM
1 2 ZA X 1
He = ri +
i=1
2 i=1
r iA j>i
r ij
A=1
H| Gi = EG | Gi
~
h (✓)|H| ~
(✓)i
EG
~ (✓)i
h (✓)| ~
<latexit sha1_base64="RvrZaf+U6pdsiZ0ltLPkNA3tPyM=">AAACXXicdVFLSwMxGMyu7/qqelDwEiyCXsquFNSbIKLHClaFppRs+m0bzGbX5FuhrP2T3vTiXzHdVlCxA4Fh5ps8JlGmpMUgePf8ufmFxaXllcrq2vrGZnVr+96muRHQEqlKzWPELSipoYUSFTxmBngSKXiIni7H/sMLGCtTfYfDDDoJ72sZS8HRSd0qsthwUTDFdV8Ba1p5xF5AMBwA8uPXm9e/EjPl5GhmZFaAsj4806tucT3qVmtB/bwEnZDTxpSchzSsByVqZIpmt/rGeqnIE9AoFLe2HQYZdgpuUAq3dYXlFjIunngf2o5qnoDtFGU7I3rolB6NU+OWRlqqPxMFT6wdJpGbTDgO7F9vLP7ntXOMzzqF1FmOoMXkoDhXFFM6rpr2pAGBaugIF0a6u1Ix4K5udB9ScSV8v5TOJvcn9TCoh7eN2sXetI5lsk8OyBEJySm5IDekSVpEkA+PeCtexfv0F/w1f2My6nvTzA75BX/3C3JUtt8=</latexit>
~
E(✓)
X
H= h↵ (↵)
↵
Map problem onto Paulis
✓~ Use classical optimizer
Tapering of qubits
Tapering of qubits
Software implemented Z-
gates
(Perfect, zero time)
Variational parameters:
• Pulse envelope
• Pulse length
• Amplitude
• Phase
42 © 2017 IBM Corporation
IBM Quantum
• Sequence of interleaved
arbitrary single qubit rotations
and naturally available
entangling gates
Post
rotations
Post
rotations
Post
rotations
preparation is perfect,
estimating <O> will still
1212 6 12
1010
1010
5
1010
10
have error 10
8 8 4 8
Counts
Counts
Counts
Counts
Counts
Counts
Counts
Counts
6 6 3 6
5 5 5 5 5
4 4 2 4
2 2 1 2
0 0 0 0 000 0 0
-1 -1 -1.2
-1.2 -1.1
-1.1 -1 -1 -1.2
-1.2 -1.1
-1.1 -1 -1 -1.2
-1.2
-1.2 -1.1
-1.1
-1.1 -1-1-1 -1.2 -1.1 -1 -1.2 -1.1 -1
gy Energy
Energy Energy
Energy Energy
Energy
Energy
Can affect accuracy of
Energy Energy
arXiv:1712.09913
Classical optimizer: Simultaneous Perturbation Stochastic
Approximation (SPSA) a b
arXiv:1712.09913
IBM Quantum
• Number of shots
• System initialization
Energy (Hartree)
/2 -7.
-7
-14
0
-7.
-14.2 -7
Energy (Hartree)
q,i,
- /2
j
-7.
-14.4
- c 0.
0
-14.6
E (Hartree)
0 100 200 0 100 200 0 100 200 0 100 200 0 100 200
0.
Iteration Iteration Iteration Iteration Iteration
-14.8 0
0.
-15
-15.2
Final Exp. Result
-15.4
Exact
-15.6
0 50 100 150 200 250
Iteration A. Kandala, A. Mezzacapo, et al,
Nature 549, 242-246 (2017)
Application to quantum chemistry : H2
Energy (Hartree)
(0.11279)ZZ+(0.18093)XX -1
d=0.735 Å
Dissociation d = 4 Å -1.2
H = (-0.70461)II+(0.00012)ZI+(0.00012)IZ+ -1.4
(1.6673e-10)ZZ+(0.33438)XX
-1.6
-1.8
Energy (Hartree)
C
R
2-
(0.11279)ZZ+(0.18093)XX -0.2
4
Q7 Q2
-0.4 H H
Dissociation d = 4 Å
-0.6
H = (-0.70461)II+(0.00012)ZI+(0.00012)IZ+
(1.6673e-10)ZZ+(0.33438)XX
-0.8
-1
a
-1.2
0 1 2 3 4
Interatomic distance (Angstrom)
IBM Quantum
C
5
C
50 20
6-
R
R
20
1-
R
1-
3
3
CR2-1
CR2-1
0 Q5 Q4 Q3 Q5 Q4 Q3 Q5 CR4-5 Q4 Q3
0 -7 0 -13
0
Energy (Hartree)
Energy (Hartree)
Energy (Hartree)
C
C
R
R
2-
2-
2-
-0.2
4
4
Q7 Q2 Q7 Q2 -13.5 Q7 Q2
-7.2
H H Be
-0.4 H H
-14
-7.4 Li H
-0.6
-14.5
-7.6
-0.8
-15
-1 -7.8
a b -15.5 c
-1.2 -8
0 1 2 3 4 1 2 3 4 5 1 2 3 4 5
Interatomic distance (Angstrom) Interatomic distance (Angstrom) Interatomic distance (Angstrom)
H2: 2 qubits LiH: 4 qubits BeH2: 6 qubits
5 pauli terms, 2 sets 100 pauli terms, 25 sets 165 pauli terms, 44 sets
10 0
Entanglers: 0
Entanglers: 1
Entanglers: 2
Entanglers: 3
Entanglers: 4
10 -1 Entanglers: 6
Entanglers: 8
Energy Error
10 -2
Chemical accuracy
10 -3
-4
10
10 -6 10 -5 10 -4 10 -3 10 -2
Noise Strength
A. Kandala, A. Mezzacapo, et al,
61 Nature 549, 242-246 (2017)
© 2017 IBM Corporation
The longitude problem IBM Quantum
H1 H2 H3 H4
1730-1735 1737-1739 1740-1759 1755-1759
Invention of
Bimetallic strip
62 http://collections.rmg.co.uk
© 2017 IBM Corporation
Error mitigation
IBM Quantum
Zero-noise extrapolation
Expectation value of observable of interest:
n
X
⇤ k
EK ( ) = E + ak + Rn+1 ( , L, T )
k=1
Assume experimentalist has exquisite control over incoherent noise (T1, T2)
EK ( ) = E ⇤ + aK + O( 2 )
⇤ 2 2
EK (c ) = E + aK c + O(c )
2 cEK ( ) EK (c ) ⇤ 2
ÊK ( )= = E + O( )
c 1
N
With N measurements, can reduce error in estimate to O( )
How
64
does one measure E(ci ) ?? © 2017 IBM Corporation
IBM Quantum
0.5 0.5
0 0
-0.5 -0.5
-1 -1
1 1
0 0 0
-0.5 -1 -1 -0.5
0 0 -0.5 0
0.5
|0> |0>
0.5 0.5 -1
|0>
-1 1 -1 1 1
Li H Li H
Kandala et al Kandala et al
Nature
67 549, 242-246 (2017) Nature 567, 491–495 (2019) © 2017 IBM Corporation
IBM Quantum
C
5
C
50 20
6-
R
R
20
1-
R
1-
3
3
CR2-1
CR2-1
0 Q5 Q4 Q3 Q5 Q4 Q3 Q5 CR4-5 Q4 Q3
0 -7 0 -13
0
Energy (Hartree)
Energy (Hartree)
Energy (Hartree)
Techniques go
C
C
R
R
2-
2-
2-
-0.2
4
4
Q7 Q2 Q7 Q2 -13.5 Q7 Q2
-7.2
H H Be
beyond just
-0.4 H H
-14
-7.4 Li H
-0.6
chemistry …
-14.5
-7.6
-0.8
-15
-1 -7.8
a b -15.5 c
-1.2 -8
0 1 2 3 4 1 2 3 4 5 1 2 3 4 5
Interatomic distance (Angstrom) Interatomic distance (Angstrom) Interatomic distance (Angstrom)
68 © 2017 IBM Corporation
THANK YOU