Accurate and decoherence-protected adiabatic quantum computation Daniel Lidar, USC PRL 100, 160506 (2008) JMP 50, 102106 (2009) Workshop on Quantum Algorithms, Computational Models, and Foundations of Quantum Mechanics Vancouver, July 23, 2010 $: Motivation Protect Adiabatic Quantum Computation (AQC) against decoherence and control errors Find rigorous bounds on …delity of AQC in presence of coupling to environment Main result of this talk: Theorem: AQC can be performed with a …delity approaching size, using only 1 as a power law in the system 2-local Hamiltonians, in the presence of 1-local noise, assuming access to (dynamical decoupling) pulses whose width and intervals shrink as a power law in the system size. The power law exponent is linear in the dynamical critical exponent of the closed system. Main reference: DAL, “Towards Fault Tolerant Adiabatic Quantum Computation”, Phys. Rev. Lett. 100, 160506 (2008) Motivation Protect Adiabatic Quantum Computation (AQC) against decoherence and control errors Find rigorous bounds on …delity of AQC in presence of coupling to environment Main result of this talk: Theorem: AQC can be performed with a …delity approaching size, using only 1 as a power law in the system 2-local Hamiltonians, in the presence of 1-local noise, assuming access to (dynamical decoupling) pulses whose width and intervals shrink as a power law in the system size. The power law exponent is linear in the dynamical critical exponent of the closed system. Main reference: DAL, “Towards Fault Tolerant Adiabatic Quantum Computation”, Phys. Rev. Lett. 100, 160506 (2008) In other words: You can get arbitrarily accurate open-system AQC for the following price: Design your adiabatic evolution Hamiltonian (2-local) so that it is analytic and has N zero derivatives at the initial and …nal times. Apply dynamical decoupling pulses generated by global magnetic …elds in the x and z directions. The shorter the pulses and pulse intervals, the higher the …delity. This result doesn’t depend on bath temperature. To ensure that adiabatic evolution and dynamical decoupling are compatible, encode into a simple (distance-2) quantum error detection code. Scheme for dynamical-decoupling-protected universal adiabatic QC 1 2 j+1 n σ x ⊗σ x σ z ⊗σ z Xj Zj switch interactions adiabatically, locally 1 i+1 j+1 σ x ⊗σ x XiX j σ z ⊗σ z ZiZ j n Scheme for dynamical-decoupling-protected universal adiabatic QC 1 2 j+1 n σ x ⊗σ x σ z ⊗σ z Xj Zj switch fields rapidly but globally Bz Bz Bz Bz 1 i+1 j+1 n σ x ⊗σ x Bx Bx XiX j σ z ⊗σ z ZiZ j Bx Proof Strategy (and talk outline) Let S distance between desired ground state and actual system state, at …nal time T Lemma 1: the distance inequality : S “decoherence distance” due to open system + non-unitary evolution} | {z dD “adiabatic distance” due to closed system deviations from adiabaticity} | {z ad and show that both these distances can be made arbitrarily small. How? Lemma 2 — adiabatic distance ad: For slow enough evolution, can be made arbitrarily small using analytic interpolation and almost-constant boundary conditions (or using only C 2 interpolation provided evolving more slowly, for the same error) Lemma 3 (with the help of Lemmas 4 & 5) — decoherence distance dD: Can be made arbitrarily small using dynamical decoupling (with …nitewidth pulses and …nite pulse intervals) Tools: Distance Measure and Operator Norm Distance between states: trace distance 1 D[ 1; 2] k 1 2 k1 2 X kAk1 TrjAj = sing.val.(A); jAj p AyA When applied to pure states i = j iih ij I’ll write D[ 1; 2]. Operator norm kAk sup kj ik=1 q h jAyAj i = max sing.val.(A) Distances X 0 : H SB = 0 X ad : perfectly adiabatic φad (T ) φad (T ) ⊗ ρ B (T ) 0 = ρad (T ) dS ρ (T ) , d = system+bath Distances X 0 : H SB = 0 X ad : perfectly adiabatic dS φad (T ) φad (T ) ⊗ ρ B (T ) 0 = ρad (T ) ρ (T ) dad dD ρ 0 (T ) , d = system+bath Distances 0 φad (T ) φad (T ) = ρS,ad (T ) δS ρS (T ) TrB TrB X 0 : H SB = 0 X ad : perfectly adiabatic dS φad (T ) φad (T ) ⊗ ρ B (T ) 0 = ρad (T ) ρ (T ) dad dD ρ 0 (T ) , d = system+bath ,δ = system only Distances δS 0 φad (T ) φad (T ) = ρS,ad (T ) ρS (T ) TrB TrB dS φad (T ) φad (T ) ⊗ ρ B (T ) 0 = ρad (T ) TrB X 0 : H SB = 0 X ad : perfectly adiabatic ρ (T ) dad dD ρ 0 (T ) φad (T ) δ ad TrB ψ (T ) , d = system+bath ,δ = system only Distances δS 0 φad (T ) φad (T ) = ρS,ad (T ) ρS (T ) TrB TrB dS φad (T ) φad (T ) ⊗ ρ B (T ) 0 = ρad (T ) TrB X 0 : H SB = 0 X ad : perfectly adiabatic ρ (T ) dad dD ρ 0 (T ) φad (T ) δ ad , d = system+bath ,δ = system only TrB ψ (T ) δS ≤ δad + dD Distances δS 0 φad (T ) φad (T ) = ρS,ad (T ) ρS (T ) TrB TrB dS φad (T ) φad (T ) ⊗ ρ B (T ) 0 (T ) = ρad TrB X 0 : H SB = 0 X ad : perfectly adiabatic ρ (T ) dad dD ,δ = system only ρ 0 (T ) φad (T ) δ ad TrB , d = system+bath open system: reduce via dynamical decoupling ψ (T ) closed system: reduce via improved adiabatic theorem δS ≤ δad + dD Part 1: Closed System AQC Dimensionless time: = t=T 2 [0; 1]; T = …nal time True …nal state j (1)i is the solution of the (rescaled) Schrödinger equation: dj i = iT Hadj i d Goal of AQC: simultaneously minimize T (n) and the error ad Questions... D[ (1); ad(1)]: What determines T (n)? How to make ad small? Textbook criterion: to have jh (1)j ad(1)ij2 need T 1 2 1 maxs jhexcited( )jH_ ad( )j ad( )ij : 2 min gap ( ) Gap dependence on n determines T (n). For e¢ cient AQC algorithms: (n) 0n z; z = dynamical critical exponent Unfortunately this criterion is not quite right... Adiabatic Distance for Closed Systems Depending on the di¤erentiability of Had one can prove di¤erent versions of the adiabatic theorem (Lemma 2). Jansen, Ruskai, & Seiler [J. Math. Phys. 48, 102111 (2007)]. Assume: – the ground state manifold of Had is gapped. – Had( ) is twice di¤erentiable on [0; 1] – H_ ad(0) = H_ ad(1) = 0 Then (norm is operator norm): T r kH_ adk2 3 =) 2 < r ad DAL, A. Rezakhani, A, Hamma, JMP 50, 102106 (2009): Assume Had (τ ): F has a non-degenerate and gapped ground state F is analytic in a strip of width γ around [0,1] F its first N derivatives vanish at τ = 0, 1 Let r > 1. Then T = r N kH_ adk2 3 =) ad < (N + 1) +1r N =) Closed-system AQC is resilient against control errors which preserve gap but cause Had(s) to deviate from its intended path, as long as …nal Hamiltonian is correct. This is a form of inherent fault tolerance to control errors not shared by the circuit model! How to accomplish this? Design Hamiltonian according to criteria above. Part 2: Open System AQC Joint System-Bath Evolution H (t) = HS (t) HS (t) = HSB = | IB + IS {z H0 H (t) | ad {z } implements AQC X S B + HB + HSB ; } HC(t) | {z } implements DD [We’ll see later how to enforce [Had(t); HC(t0)] = 0 8t; t0:] Joint system-bath propagator and joint state: U (t) = T exp[ i Z t 0 H (t0)dt0 ]; (t) = U (t) (0)U (t)y Decoupled joint system-bath state: 0(t) = U 0(t) (0)U 0(t)y Decoupling distance: dD D[ (T ); 0(T )] S dD + ad: Already saw ad can be made arbitrarily small using analytic interpolation. Goal: Simultaneously minimize dD using dynamical decoupling. This is an optimization problem: generically decoherence worsens with increasing T , while closed-system adiabaticity improves. Dynamical Decoupling (DD) = a sequence of pulses applied to the system, sometimes forming a group G , designed to reduce the e¤ective system-bath coupling. Implemented via HC(t). The sequence ZOO, in increasing order of performance quality: PDD = a periodic repetition of a basic sequence RDD = a random pulse sequence CDD = a concatenated sequence (recursively structured) UDD = a sequence optimized to cancel pure qubit dephasing with the smallest possible number of pulses QDD = a sequence optimized to cancel general qubit decoherence with the smallest possible number of pulses Dynamical Decoupling Theory “Symmetrizing group” of pulses { gi } and their inverses are applied in series: ( g N† fg N )" ( g 2†fg 2 )( g1†fg1 ) ≈ exp(−iτ ∑ i gi † H SB gi ) f ≡ exp(−iH SBτ ) Choose the pulses so that: (1) H SB 6 H eff ≡ ∑ i gi † H SB gi = 0 first order Magnus expansion Dynamical Decoupling Condition For a qubit the Pauli group G={X,Y,Z,I } (π pulses around all three axes) removes an arbitrary HSB : (XfX)(YfY)(ZfZ)(IfI) = XfZfXfZf Periodic DD: periodic repetition of the universal DD pulse sequence DD Parameters: K = no. of pulses w = pulse width = pulse interval c = K (w + ) = cycle time L = number of PDD cycles T = L c = total time Transform to interaction picture de…ned by Had + HB , i.e., ~ (t) = U y (t) U ad y UB (t)U (t) The “e¤ective error Hamiltonian”: ~ (t) U e itHe (t) He (t) can be calculated using Dyson or (better) Magnus expansion. Note that in decoupled limit: t=j ~ (t) HSB!0 ! UC(t) = U c I , i.e., He (j c) = 0: Deviation from ideality at …nal time is thus quanti…ed by “error phase”: T kHe (T )k How is the “error phase” related to the “decoupling distance” dD? “error phase”: T kHe (T )k Then distance between desired and actual (Lemma 4): D[ (T ); 0(T )] = dD min[1; (e if 1)=2] 1: Minimize error phase ) mininize decoupling distance. [DAL, P. Zanardi & K. Khodjasteh, Phys. Rev. A 78, 012308 (2008)] Error Phase Bound for Periodic DD [K. Khodjasteh & DAL, Phys. Rev. A 78, 012355 (2008)] Let J and assume J c < kHSB k kHad + HB k (absolute convergence condition of Magnus expansion). Then (Lemma 5): (T ) JT w | + {z w} error due to …nite pulse width (JT )2 exp(2 + + JT min[1; ( L 2 {z | c) c error due to terms not removed by …rst order DD procedure Use this to bound decoupling distance. 1 1)]: } Joint AQC-DD Optimization for Local Hamiltonians For local Hamiltonians Had and HB : = kHad + HB k O(n2) Recall closed system adiabaticity theorem, ensuring error small if N is large: T = r r N N kH_ adk2 3 1 3z+4 n 3 =) ad < (N + 1) +1r N KL( + w) 0 where we used gap condition for e¢ cient AQC algorithms: (n) 0n z; z = dynamical critical exponent Given and …xed parameters of the problem are J , 0, and z . Need to ensure that each of the terms upper bounding function of n. This is the case that if pulse interval n (3z+3+ 1)= 0; with 1 > 0 and 2 > 0. (T ) vanishes as a and pulse width w scale as w n (6z+5+ 1+ 2)=J For then, using Lemmas 4 & 5 we have proven Lemma 3: dD . n n!1 + (J= 0)2n (1+ 1) + (J= 0)n 1 ! 0 =) Using PDD with properly chosen parameters we can obtain arbitrarily accurate AQC. 2 Shortcoming: pulse intervals and widths must shrink with n as a power law... Could perhaps be remedied by employing concatenated DD. Need to ensure that each of the terms upper bounding function of n. This is the case that if pulse interval n (3z+3+ 1)= 0; with 1 > 0 and 2 > 0. (T ) vanishes as a and pulse width w scale as w n (6z+5+ 1+ 2)=J For then, using Lemmas 4 & 5 we have proven Lemma 3: dD . n n!1 + (J= 0)2n (1+ 1) + (J= 0)n 1 ! 0 =) Using PDD with properly chosen parameters we can obtain arbitrarily accurate AQC. 2 Shortcoming: pulse intervals and widths must shrink with n as a power law... Could perhaps be remedied by employing concatenated DD. In conclusion we’ve (almost) proven the Theorem: AQC can be performed with a …delity approaching 1 as a power law in the system size, [using only 2-local Hamiltonians, in the presence of 1-local noise], assuming access to (dynamical decoupling) pulses whose width and intervals shrink as a power law in the system size. The latter power law exponent is linear in the dynamical critical exponent of the closed system. Part 3: Seamlessly Combining AQC & DD Need to show how to achieve “non-interference” condition [Had(t); HC(t0)] = 0 8t; t0: Can be done using encoding. Example — Stabilizer-normalizer construction: DD pulses are stabilizer elements. AQC implemented via normalizer elements. 2-Local Universal AQC Resistant Against 1-local Noise A 2-local Hamiltonian that is universal for AQC (J.D. Biamonte and P.J. Love, arXiv:0704.1287): univ (t) Had = X X hi (t) i + i; 2fx;zg i;j; 2fx;zg 1-local noise (linear decoherence model): lin = HSB X n X =x;y;z j=1 j Bj Jij (t) i j : lin = Decoupling group that decouples HSB P G = fI; X; Y; Zg; X = Pn =x;y;z j=1 j n O x, j Bj : etc. j=1 This requires only global pulses. G is the stabilizer of an [[n; n weight binary string): 2; 2]] stabilizer code C (n even, x =even p C = fj xi = (jxi + jnot xi) = 2g E.g. n = 4: p p j00iL = (j0000i + j1111i) = 2; j10iL = (j0011i + j1100i) = 2 p p j01iL = (j0101i + j1010i) = 2; j11iL = (j1001i + j0110i) = 2 p Encoded single-qubit generators for C = fj xi = (jxi + jnot xi) = 2g: Xj = x1 xj+1 Zj = zj+1 zn Encoded two-qubit generators: XiXj = xi+1 xj+1 ZiZj = zi+1 zj+1 Thus universal AQC can be combined with DD using only 2-local z z interactions over C . i j x x i j and Universality for the [[n,n-2,2]] code 1 2 j+1 n σ x ⊗σ x σ z ⊗σ z Xj Zj Bz Bz Bz Bz 1 i+1 j+1 n σ x ⊗σ x Bx Bx XiX j σ z ⊗σ z ZiZ j Bx Physical examples where X , Z (as pulses for DD) and xi xj; zi zj (as Hamiltonians for AQC) are available and controllable: Capacitively coupled ‡ux qubits (D.V. Averin and C. Bruder, Phys. Rev. Lett. 91, 057003 (2003)) Spin models implemented with polar molecules (A. Micheli, G. Brennen, and P. Zoller, Nature Phys. 2, 341 (2006)) Conclusions 1 as a power law in the system size n, using only 2-local Hamiltonians, in the presence of 1-local noise, assuming Theorem: AQC can be performed with a …delity approaching access to (dynamical decoupling) pulses whose width and intervals shrink as a power law in the system size. The latter power law exponent is linear in the dynamical critical exponent of the closed system. Open questions: - Can a similar result be shown with n-independent pulse width and interval? - The distance bound is rather crude because of use of the triangle inequality; can it be tightened by directly treating adiabatic evolution in the open system? - A fault-tolerance threshold for AQC? Phys. Rev. Lett. 100, 160506 (2008) Conclusions 1 as a power law in the system size n, using only 2-local Hamiltonians, in the presence of 1-local noise, assuming Theorem: AQC can be performed with a …delity approaching access to (dynamical decoupling) pulses whose width and intervals shrink as a power law in the system size. The latter power law exponent is linear in the dynamical critical exponent of the closed system. Open questions: - Can a similar result be shown with n-independent pulse width and interval? - The distance bound is rather crude because of use of the triangle inequality; can it be tightened by directly treating adiabatic evolution in the open system? - A fault-tolerance threshold for AQC? Phys. Rev. Lett. 100, 160506 (2008) Odds and Ends Energy gap stabilization Jordan, Farhi & Shor (PRA 74, 052322 (2006)) introduced a scheme that can protect AQC against 1-local noise using a energy gap against local excitations. Such gaps can be engineered into the Hamiltonians using error detecting codes However, the resulting universal Hamiltonians are at least 4-local Their scheme is fully compatible with the DD scheme presented here Control Error ● Suppose final Hamiltonian is slightly off: ● Resulting ground state is: Gap against local operators ensures denominator gap protection ● Why not Use Standard Quantum Error Correction? Not known how to embed fault-tolerant QEC into adiabatic evolution: Requires feedback, which may break adiabaticity (A. Allahverdyan & G. Mahler, arXiv:0804.1643) Embedded FT-QEC circuit will face di¤erent error model than what it was designed for We used dynamical decoupling (open-loop, feedback-free) to deal with coupling to environment, and showed it can deal with the relevant error model. And, we used an analytic interpolation to boost …delity of closed-system adiabatic evolution. Is AQC “Better” than Other QC Paradigms? No: Computationally equivalent to circuit model (Aharonov et al., :quant-ph/0405098; A. Mizel, DAL, M. Mitchell, PRL 99, 070502 (2007)) Doesn’t have a fault tolerance theory to back it up Maybe: Perhaps more easily implemented in certain systems, in particular solid state (quantum dots, superconducting) Inherently protected against leakage if kB T is smaller than gap Inherently protected against path deformations: only end points matter Open System Evolution Consider the uncoupled setting HSB = 0, denoted by the superscript 0. ideal adiabatic state: 0S;ad(t) = j ad(t)ih ad(t)j 0 (t) = j S 0 (t) C 0 (t) B actual state under Had(t) : state under HC(t) : state under HB (t) : decoupled joint state: ideal adiabatic joint state: 0 (t) = 0S (t) 0 (t) ad (t)ih (t)j 0 (t) 0 (t) B C 0 (t) 0 (t) 0 (t) B S;ad C Let d ( ) denote distances in the joint (system) Hilbert space. target distance: S decoupling distance: dD adiabatic distance: dad D[ S (T ); 0S;ad(T )] D[ (T ); 0(T )] D[ 0(T ); 0ad(T )] = ad Proof of the Distance Inequality Partial trace can only decrease distance: D[ S (T ); 0S;ad(T )] {z } | D[ (T ); 0ad(T )] S Triangle inequality: D[ (T ); 0ad(T )] D[ (T ); 0(T )] + D[ 0(T ); 0ad(T )] {z dD | } =) desired distance inequality (Lemma 1): S dD + ad: | {z ad } Dynamical Decoupling Inspired by Hahn spin echo (1950) π/2-Y Pulses Signal (|0〉 + |1〉) π-Z τ τ (|0〉 + |1〉) FID – T2* T Echo Dynamical Decoupling (DD) = a sequence of pulses applied to the system, sometimes forming a group G , designed to reduce the e¤ective system-bath coupling. Implemented via HC(t). The sequence ZOO, in increasing order of performance quality: PDD = a periodic repetition of a basic sequence RDD = a random pulse sequence CDD = a concatenated sequence (recursively structured) UDD = a sequence optimized to cancel pure qubit dephasing with the smallest possible number of pulses QDD = a sequence optimized to cancel general qubit decoherence with the smallest possible number of pulses Subsystem Code Construction for [Had(t); HC(t0)] = 0 The decoupling group G induces a decomposition of the system Hilbert space HS via its group algebra CG and its commutant CG 0 : HS = CG = M J M CnJ InJ CdJ ; MdJ ; J CG 0 = M MnJ IdJ : J nJ = multiplicity of irrep J ; dJ = dimension Adiabatic state is encoded into a left factor CJ log nJ qubits. AQC is enacted via CG 0. CnJ : an nJ -dimensional code CJ storing DD pulses act on the right factors, enacted via the elements of DD pulses project each S in system-bath Hamiltonian to Non-interference condition is satis…ed because L CG . J [CG; CG 0] = 0. J; InJ IdJ .
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- Workshop on Quantum Algorithms, Computational Models and Foundations of Quantum Mechanics /
- Accurate and decoherence-protected adiabatic quantum...
Open Collections
Workshop on Quantum Algorithms, Computational Models and Foundations of Quantum Mechanics
Accurate and decoherence-protected adiabatic quantum computation Lidar, Daniel 2010
Page Metadata
Item Metadata
Title | Accurate and decoherence-protected adiabatic quantum computation |
Creator |
Lidar, Daniel |
Date Issued | 2010-07-23 |
Description | In the closed system setting I will show how to obtain extremely accurate adiabatic QC by proper choice of the interpolation between the initial and final Hamiltonians. Namely, given an analytic interpolation whose first N initial and final time derivatives vanish, the error can be made to be smaller than 1/N^N, with an evolution time which scales as N and the square of the norm of the time-derivative of the Hamiltonian, divided by the cube of the gap (joint work with Ali Rezakhani and Alioscia Hamma). In the open system setting I will describe a method for protecting adiabatic QC by use of a hybrid encoding-dynamical decoupling scheme. This strategy can be used to protect spin-based universal adiabatic QC against arbitrary 1-local noise using only global magnetic fields. By combining error bounds for the closed and open system settings, I will show that in principle the method is scalable to arbitrarily large computations. |
Type |
Moving Image Other |
Language | eng |
Collection |
Workshop on Quantum Algorithms, Computational Models and Foundations of Quantum Mechanics |
Date Available | 2010-11-22 |
Provider | Vancouver : University of British Columbia Library |
DOI | 10.14288/1.0103158 |
Affiliation |
Non UBC |
Peer Review Status | Reviewed |
Scholarly Level | Faculty |
URI | http://hdl.handle.net/2429/30072 |
Aggregated Source Repository | DSpace |
Download
- Media
- [if-you-see-this-DO-NOT-CLICK]
- [if-you-see-this-DO-NOT-CLICK]
- Lidar-AQCDD.pdf [ 2.12MB ]
- [if-you-see-this-DO-NOT-CLICK]
- WS Jul 23 Lidar.mp4 [ 150.76MB ]
- Metadata
- JSON: 1.0103158.json
- JSON-LD: 1.0103158+ld.json
- RDF/XML (Pretty): 1.0103158.xml
- RDF/JSON: 1.0103158+rdf.json
- Turtle: 1.0103158+rdf-turtle.txt
- N-Triples: 1.0103158+rdf-ntriples.txt
- Original Record: 1.0103158 +original-record.json
- Full Text
- 1.0103158.txt
- Citation
- 1.0103158.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Country | Views | Downloads |
---|---|---|
United States | 18 | 2 |
France | 4 | 1 |
China | 4 | 29 |
Germany | 2 | 0 |
Canada | 2 | 0 |
City | Views | Downloads |
---|---|---|
Lewes | 8 | 0 |
Unknown | 6 | 1 |
Ashburn | 3 | 0 |
Shenzhen | 2 | 13 |
San Francisco | 2 | 0 |
Plano | 2 | 0 |
Beijing | 2 | 16 |
Wilmington | 1 | 0 |
Ottawa | 1 | 0 |
Montreal | 1 | 0 |
Sunnyvale | 1 | 0 |
Seattle | 1 | 0 |
{[{ mDataHeader[type] }]} | {[{ month[type] }]} | {[{ tData[type] }]} |
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.29986.1-0103158/manifest