arxiv: v1 [quant-ph] 6 Oct 2012

Transcription

1 A State Dstllaton Protocol to Implement Arbtrary Sngle-qubt Rotatons arxv:.98v [quant-ph] 6 Oct Gullaume Duclos-Canc and Krysta M. Svore Départment de Physque Unversté de Sherbrooke Sherbrooke Québec JK R (Canada) Quantum Archtectures and Computaton Group Mcrosoft Research Redmond WA 985 (USA) (Dated: October 9 ) An mportant task requred to buld a scalable fault-tolerant quantum computer s to effcently represent an arbtrary sngle-qubt rotaton by fault-tolerant quantum operatons. Tradtonally the method for decomposng a sngle-qubt untary nto a dscrete set of gates s Solovay-Ktaev decomposton whch n practce produces a sequence of depth O(log c (/ɛ)) where c 3.97 s the state-of-the-art. The proven lower bound s c = however an effcent algorthm that saturates ths bound s unknown. In ths paper we present an alternatve to Solovay-Ktaev decomposton employng state dstllaton technques whch reduces c to between. and.7 dependng on the settng. For a gven sngle-qubt rotaton our protocol sgnfcantly lowers the length of the approxmatng sequence and the number of requred resource states (ancllary qubts). In addton our protocol s robust to nose n the resource states. PACS numbers: 3.67.Lx 3.67.Pp 3.65.Fd Keywords: state dstllaton Solovay-Ktaev decomposton Gven recent progress n quantum algorthms quantum error correcton and quantum hardware a scalable quantum computer s becomng closer and closer to realty. For many proposed quantum computer archtectures e.g. topologcal systems based on the bradng of non-abelan anyons [ 5] or the surface code model based on code deformaton [6 7] so-called Clfford operatons and stablzer state preparatons or measurements can be mplemented effcently and accurately. However these operatons alone are not suffcent for quantum unversalty snce they can be smulated classcally [8 ]. One technque to acheve quantum unversalty s to use magc state dstllaton [ 3] to augment the set wth a sngle non-clfford operaton e.g. the sngle-qubt π/8 gate T. Ths augmented set can be used to approxmate any sngle-qubt untary usng Solovay-Ktaev decomposton [4]. The Solovay-Ktaev theorem [5] states that for any ɛ and sngle-qubt gate U U can be approxmated to precson ɛ usng Θ(log c (/ɛ)) gates drawn from a unversal dscrete gate set where c s a small constant. Stateof-the-art mplementatons of Solovay-Ktaev decomposton result n c 3.97 [4 6] resultng n an average decomposton sequence wth hundreds to thousands of T gates [6]. Each T gate requres a number of copes of a quantum magc state H (Hadamard +-egenstate) where the number depends on the specfc state dstllaton protocol and purty of the state [ 3]. Therefore t s especally mportant to mnmze the number of T gates when decomposng a rotaton snce each T gate requres addtonal ancllary qubts for mplementaton. In ths paper we present an alternatve protocol to Solovay-Ktaev decomposton that allows the mplemen- Gullaume.Duclos-Canc@USherbrooke.ca ksvore@mcrosoft.com taton of sngle-qubt rotatons through the use of dstlled magc H states. We show that the resources requred by our protocol are substantally fewer than the resources requred by state-of-the-art mplementatons of Solovay-Ktaev decomposton [6] n both the number of gates and the number of quantum magc states necessary to apply an arbtrary sngle-qubt untary. I. MAGIC STATE DISTILLATION AND IMPLEMENTATION OF ROTATIONS Solovay-Ktaev decomposton [4 6] enables the approxmaton of any gate usng an approxmately unversal set of elementary gates e.g. {H T Λ(X)} where Λ(X) denotes the controlled-not gate. In partcular one can approxmate any sngle-qubt untary operaton usng the set {H T }. Magc state dstllaton s then used to produce the magc H states necessary to mplement the T gate. We call a state ψ magc f gven n nosy copes of ψ and the ablty to perform perfect Clfford operatons we can obtan or dstll a purer copy of ψ from a Clfford crcut appled to the n nosy copes of ψ. We can then obtan even purer states whch are arbtrarly close to the perfect state by applyng the protocol recursvely [ 3]. These dstlled states can be used to mplement non- Clfford operatons e.g. the T gate as descrbed below. We brefly revew how to perform an arbtrary rotaton about the Z-axs usng a resource state and how to apply the T gate. We assume that Clfford operatons are appled perfectly snce they can be mplemented faulttolerantly and that the resource states are very close to pure. One can use the protocols of [ ] to perform the ntal dstllaton n order to obtan arbtrarly pure resource states for nput. We focus on the + egenstate

2 Z(θ) X(Z) m ψ Z(X)( m θ) ψ FIG.. Crcut randomly mplementng a rotaton of angle ±θ around the Z(X)-axs. of the Hadamard operaton H H = cos π 8 + sn π 8 and omt the cost of the ntal dstllaton of the H state. We concentrate on sngle-qubt states found n ether the XZ- or XY -plane of the Bloch sphere. Note that one can easly move a state from one plane to the other by the applcaton of the Clfford HSHX operaton. Suppose we start wth states Z(θ) and ψ: Z(θ) = + e θ ψ = a + b. The crcut to mplement a rotaton around the Z-axs presented n Fg. leads to the two-qubt state Z(θ) ψ = a + b + ae θ + be θ Λ(X) a + b + ae θ + be θ. Upon measurement of the frst qubt n the computatonal bass we obtan m= a + be θ m= ae θ + b = a + be θ each wth probablty /. Thus the angle of rotaton s chosen at random to be θ or θ up to global phase. An analogous crcut performs a rotaton about the X-axs. Smlar crcuts can be found n [7]. As an mportant example of ths procedure consder the XY -plane verson of the H state: Z(π/4) = HSHX H = + e π/4. Usng the crcut n Fgure we can mplement a Z- rotaton of angle ±π/4 producng at random ether the T gate or ts adjont T. In ths partcular case we can determnstcally apply the desred gate T or T by applyng the phase gate S snce ST = T. In general however ths determnstc correcton wll not be possble. II. NEW STATES FROM H STATES In ths secton we show that we can use a very smple two-qubt Clfford crcut to obtan other non-stablzer states usng only H states as an ntal resource and then show that these states enable the approxmaton of H X ( ) H H + ( H ) FIG.. Two-qubt crcut used to obtan new H states from ntal resource states H. Upon measurng the outcome the output state s H +. Upon measurng the outcome the output state s H. any sngle-qubt rotaton. We assume that we are provded wth perfect copes of H. We would lke to mnmze the number of H states requred to mplement an arbtrary sngle-qubt rotaton snce these dstlled states can be costly to produce. Consder the crcut of Fg.. One can easly verfy that t measures the party of the two nput qubts and decodes the resultng state nto the second qubt. We begn by consderng the two nputs to be H states. We defne θ = π 8 and H = H = cos θ + sn θ. The crcut begns as: H H = cos θ + sn θ + cos θ sn θ ( + ) Λ(X) cos θ + sn θ + cos θ sn θ ( + ). Upon measurement of the frst qubt we have m= cos θ + sn θ cos 4 θ + sn 4 θ m= ( + ). We defne θ such that cos θ + sn θ = cos θ + sn θ cos 4 θ + sn 4 θ from whch we deduce cot θ = cot θ. We defne H = cos θ + sn θ a non-stablzer state obtaned from H states Clfford operatons and measurements. If the outcome of the measurement s then we obtan a stablzer state and dscard the output (see Fg. ). The two measurement outcomes occur wth respectve probabltes p = cos 4 θ + sn 4 θ = 3 4 p = p = 4. We now recurse on ths protocol usng the nonstablzer states produced by the prevous round of the protocol as part of the nput to the crcut of Fg.. We defne H = cos θ + sn θ cot θ = cot + θ.

3 3 FIG. 3. Process of obtanng other non-stablzer states from ntal H states. A copy of H and H probablstcally yeld a copy of H or H + usng the crcut of Fg.. Each step along the ladder costs one copy of H except the frst one whch costs two. θ θ If we use as nput a copy of H and a copy of H we have H H = cos θ cos θ + sn θ sn θ + sn θ cos θ + cos θ sn θ Λ(X) cos θ cos θ + sn θ sn θ + sn θ cos θ + cos θ sn θ. TABLE I. Rotaton by angle θ H = cos θ + sn θ state. mplementable usng an Upon measurement of the frst qubt we have where m= (cos θ + sn θ ) m= (cos θ + sn θ ) cot θ = cot θ cot θ = cot + θ = cot θ + cot θ = cot θ tan θ = cot θ = cot θ. Thus f we measure we obtan the state H + and f we measure we obtan the state H. The probablty of measurng s gven by p = cos θ cos θ + sn θ sn θ. Note that 3 4 p < cos π 8 = We can vew ths recursve process as a sem-nfnte random walk wth based non-homogeneous probabltes as Fg. 3 llustrates. Every tme a step s taken along ths ladder of states one copy of H s consumed except at the frst node of the ladder (the producton of state H ) when we requre two copes of H: f the outcome s measured at the frst node we dscard the output and start wth two new copes of H. Table I lsts the rotatons obtaned from the frst few recursons usng the H states and Fgs. 4 5 llustrate. Note that there s a factor of two dfference between the angle θ nvolved n the descrpton of the state and the rotaton appled e.g. the H state s over θ = π 8 and can be used to mplement a π 4 rotaton. Also as < θ < π 4 ( ) the dscontnuty of the cotangent s never a problem. III. NUMERICAL STUDY OF RANDOM Z-ROTATIONS In ths secton we numercally study the cost of mplementng sngle-qubt rotatons. Although the crcut n 4 FIG. 4. Red dots: Rotaton by angle θ mplementable usng an H = cos θ + sn θ state. Fg. randomly apples θ or θ we show that our protocol stll results n an effcent applcaton of the desred Z-rotaton. Assume that we have the deal case where θ = θ +. In ths hypothetcal scenaro f we try to apply some rotaton usng H and t fals then we can correct the gate by applyng a rotaton usng H. If ths gate also fals then we follow wth H and so on. There are two crucal facts to pont out. Frst the probablty of falng n tmes n a row scales as / n.e. t decays exponentally wth n such that the expected number of teratons s well-behaved. Second and very mportantly f a T gate fals (usng a H state) we can always correct determnstcally by applyng the phase gate S (whch s a Z-rotaton by e π/ ). Unfortunately θ θ + ; however ths assumpton s not too far from the truth (see Fg. 5). Snce we only requre that a gate be approxmated to a gven precson we can actually apply any Z-rotaton rapdly and wth good accuracy. Frst apply the H -rotaton that gets you closest to your target angle and then recurse. The next gate to apply wll depend on whether the current gate succeeded or not. Due to the property dscussed n the prevous paragraph we wll rapdly converge towards our target angle. We now present smulaton results of obtanng the Z

4 (a)ft: ln(c on) = ln(ln(/ɛ)). FIG. 5. Dots: States obtanable by recursvely usng the crcut of Fg. and only H states as ntal nput. Full lne: exponental decay ft θ rotaton Z(φ) where φ s chosen randomly n order to characterze the effcency of the protocol proposed n prevous sectons. The smulaton proceeds as follows:. Set desred accuracy ɛ.. Randomly pck a target rotaton angle < φ < π. 3. Fnd the state H such that θ s close to φ. 4. Smulate an nstance of the ladder to obtan that state and add ts cost to the offlne cost. 5. Apply a rotaton usng the H state and the crcut of Fg. and add one to the onlne cost. 6. Recurse on steps 3 through 5 untl the desred accuracy s reached. We defne the accuracy of the appled rotaton V compared to the target rotaton U as where max ψ D(U ψ ψ U V ψ ψ V ) D(ρ σ) = ) ( (ρ tr σ) (ρ σ) s the trace dstance between states ρ and σ. If U and V are rotatons about the same axs one can show that for small angles of rotaton whch wll always be our case ths reduces to the dfference of rotaton angles ɛ = φ. In [6] the dstance measure used s tr(uv ) D(U V ) =. In the case of rotatons about the same axs t can be reduced to cos( φ) φ/ for small φ. Ths converson between the two measures s mportant snce we later compare performance. To compare the resource cost of our protocol to Solovay-Ktaev decomposton we defne an onlne and offlne cost to apply a untary gate. The onlne cost C on s the expected number of non-clfford gates or H states requred to mplement the untary gate on a (b)ft: ln(c Off ) = ln(ln(/ɛ)). FIG. 6. Target accuraces are chosen such that < ɛ < 4 and the sample sze s.8 4. The clouds of ponts are used to ft the data accordng to a lnear ft. We obtan c.9 and c.7. qubt. The offlne cost C off s the total number of dstlled H states requred to obtan all of the ntermedate H states used to perform the gven untary operaton that s the sum of the H states used for each ladder process. In our resource costs we do not nclude the ntal cost to dstll H states. For Solovay-Ktaev decomposton the offlne cost s always and the onlne cost s the total number of T and T gates n the decomposton. We ran the smulaton for target accuraces rangng between < ɛ < 4 each tme consderng a new random angle to produce a sample of.8 4 nstances of ths protocol. Just lke n the case of Solovay-Ktaev decomposton we suppose that C on ln c ( ɛ ) C off ln c ( ɛ ) where C on and C off are the onlne and offlne costs respectvely such that ln C on c ln ln( ɛ ) ln C off c ln ln( ɛ ). The results are gven n Fg. 6. Fts are also presented from whch we deduce that c.9 and c.7 for our protocol. As dscussed n Secton V both of these scalngs represent a sgnfcant mprovement over the best mplementaton to our knowledge of Solovay-Ktaev decomposton

5 5 [6] whch was tself a sgnfcant mprovement over the prevous mplementaton of [4]. Note that one can mplement any sngle-qubt untary U usng three rotatons around the X- and Z-axes []: U X(α)Z(β)X(γ) for some angles α β γ. We have explctly shown smulaton results for Z-rotatons however X-rotatons can be obtaned at the same cost usng the X-rotaton crcut gven n Fg.. Thus we can use our protocol to produce each of the three rotatons and produce any desred sngle-qubt untary operaton. A. Other states To further reduce the resource costs and ther respectve scalngs we show that we can use dfferent Clfford crcuts to produce new non-clfford states that can be used as ntal resources for the prevously presented protocols. We frst ntroduce three new states and ψ dscuss how to combne them nto the descrbed protocol. Consder the crcut of Fg. 7. It s a Clfford crcut to whch we nput four copes of H. The measurement outcome occurs wth probablty 3( + )/3.3 otherwse the output s dscarded. If the measurement yelds result then the produced state s ψ = cos φ + sn φ ( φ = π + 3 cot ).446. Snce the probablty of success s.3 and that every tral consumes four copes of H the average cost to produce ψ s.5 H states. Ths crcut was desgned to measure the stablzer code presented n Table II. Another nterestng state can be obtaned from the same crcut substtutng one of the nput states by a + state as s llustrated by Fg. 8. The measurement outcome s obtaned wth probablty (6 + )/3.3. The correspondng output state s ψ = cos φ + sn φ ( φ = π ) cot Snce the probablty of success s.3 and every tral consumes three copes of H the average cost to produce ψ s.95 H states. Fg. 9 presents another useful crcut. The measurement outcome s obtaned wth probablty / The correspondng output state s ψ = cos φ + sn φ φ = π ( ) 7 cot H H X H H H H X H H ψ H X Z X FIG. 7. Crcut to produce ψ states. The probablty of success s.3 and every tral consumes four copes of H such that the average cost s.5 to produce a copy of ψ. S ± 3 s + X Z X. s +. X Z X s + X. X Z Z + Z Z Z Z TABLE II. The stablzer code decoded by the crcut of Fg. 7. The probablty of success s.344 and every tral consumes four copes of H such that the average cost to produce ψ s.64 H states. Table III presents the stablzer code n terms of ts generators S that are decoded by the crcut. We wll use these states as nput states to the crcut gven n Fg. where one of these states s used n place of the H nput state. We start wth a copy of ψ and a copy of H. If measurement outcome s obtaned the state s dscarded. Otherwse we obtan ψ = cos φ + sn φ cot φ = cot φ cot θ. Smlarly to the H states we defne ψ j = cos φ j + sn φj cot φ j = cot φj cot θ. If we nput a copy of ψ j and a copy of H we obtan H ψ j Λ(X) cos θ cos φ j + sn θ sn φ j + sn θ cos φ j + cos θ sn φ j. H H X + H H H X H H ψ H X Z X FIG. 8. Crcut to produce ψ states. The probablty of success s.3 and every tral consumes four copes of H such that the average cost s.95 to produce a copy of ψ.

6 6 H H H X X H ψ H X H X FIG. 9. Crcut to produce ψ states. The probablty of success s.344 and every tral consumes four copes of H such that the average cost s.64 to produce a copy of ψ. S ± 3 s + X X X X s + Z. Z. s + Z.. Z Z + Z Z Z Z TABLE III. The stablzer code decoded by the crcut of Fg. 9. such that the output state obtaned s dependng on measurement outcome m= m= ψ j + ψ j. New ladders of states can be obtaned usng the states as nputs n place of the H states. ψ Fg. shows the four ladders. Table IV lsts the rotatons obtaned from the frst few recursons and Fg. llustrates. We see that the set of possble rotatons s more dense. We reproduced the numercal experment of Secton III wth basc offlne costs of.5.95 and.64 for ψ ψ and ψ respectvely. The results are presented n Fg. 3. Snce the set of states s denser we expected mproved scalngs for both the onlne and offlne costs. Ths s ndeed the case we fnd c. and c.75. However the basc offlne costs of our new states ψ are sgnfcantly hgher; for precson 4 even though the onlne cost s smaller usng the new states the offlne cost s stll smaller f we restrct ourselves to the smpler scheme usng only H states. For the protocol usng the new nput states to reduce both the onlne and offlne costs we need to consder precsons smaller then ɛ.8 5 see Fg.. The reason why we do not consder other measurement outcomes n the crcuts consdered n ths secton s that n general potental errors on the H states are amplfed by the crcut. Output states n these cases mght stll prove useful but a careful analyss of the evoluton of errors must be conducted FIG.. Dots: States obtanable by recursvely usng the crcut of Fg. wth ntal resource states H ψ ψ and ψ H 3 Ψ Ψ Ψ Ψ 64 Ψ H FIG.. Dots: Rotatons mplementable usng H ψ ψ ψ states. B. Mnmzng the onlne cost In ths secton we am to mnmze the onlne cost at the prce of potentally ncreasng the offlne cost. The protocol up to ths pont can be summarzed as follows. Suppose one wants to mplement a Z-rotaton of an arbtrary angle φ on a logcal state ψ. One has to mplement a sequence of j rotatons {Z(θ j )} on ψ usng the sequence of states { Hj } such that Z(φ) θ φ φ φ TABLE IV. Rotatons mplementable usng H ψ ψ ψ states.

7 7 ' ' C On C Off C On C Off FIG.. Dashed lne: Offlne (top) and onlne (bottom) costs for the scheme usng only H states. Dotted lnes: Offlne (top) and onlne (bottom) costs for the scheme usng { H ψ ψ ψ } states. The two top curves cross at ɛ.8 5. Ε '' C On '' LnC Off (a)ft: C on =.99. LnLn Ε LnLn Ε (b)ft: ln(c Off ) = ln(ln(/ɛ)). (a)ft: ln(c on) = ln(ln(/ɛ)). FIG. 4. Target accuraces are chosen such that < ɛ < 4 and the sample sze s.8 3. (a) The cloud of ponts s averaged to get the expected number of onlne rotatons. (b) The cloud of pont s used to perform a lnear ft. (b)ft: ln(c Off ) = ln(ln(/ɛ)). FIG. 3. Target accuraces are chosen such that < ɛ < 4 and the sample sze s.8 4. The clouds of ponts are used to ft the data accordng to a lnear ft. We obtan c. and c.75 mprovng on the results of Fg.6. j Z(θ j ). The onlne cost s gven by { H j }. Consder nstead the followng protocol to mplement the same rotaton by angle φ. Prepare offlne the state Z(φ) from a copy of. To acheve ths use the protocol descrbed n the prevous paragraph to rotate to Z(φ). Note that you can mplement t offlne because you are applyng rotatons on an anclla state. Then use Z(φ) onlne to apply the desred rotaton. Wth probablty the rotaton Z(φ) s appled and the onlne cost s. If t fals correct for t by preparng offlne the state Z(φ). Agan wth probablty the overall rotaton Z(φ) s appled and the onlne cost s. If t fals prepare offlne Z(4φ) and so on. The probablty that a number n of teratons s requred before success de- creases exponentally wth n. Ths process s a negatve bnomal of parameter p = and the expected number of onlne rotatons before success goes as p =. We numercally smulated ths process for varous random angles < φ < π and accuraces < ɛ < 4. As Fg. 4 llustrates we note that there s no sgnfcant change n the onlne cost for dfferent values of ɛ. Moreover the expected number of onlne rotatons to apply s roughly two the result one would expect n ths stuaton. The scalng of the offlne cost s the same as that of the prevous scheme as expected. Also even though the scalng s the same the actual values of the offlne cost are bgger. The shft s.3.64 =.59 see Fgs. 3(b) and 4(b) whch corresponds to a factor of e One would expect a slghtly bgger shft of the offlne cost by ln snce we repeat the scheme twce on average. Ths suggests that there mght exst some favorable correlatons between the expected offlne cost of an angle θ and of twce that angle θ. IV. ERRONEOUS STATES In ths secton we determne the effect of errors on our resource H states and ther effect on the produced H states. A pror the errors mght be amplfed by the two-qubt crcut of Fg. however we show ths s not the case. The probablstc and non-homogeneous nature of the

8 FIG. 5. Evoluton of the trace dstance between mperfect ρ a and perfect H states. Crcles: data for p = 4 8. Squares: data for p = 6. Damonds: data for p = 8 6. The full lnes are exponental decay fts: (.8 3 ).3 usng ponts 8 8 (.63 5 ).8 usng ponts 8 and (.6 7 ).4 usng ponts 3 6 for the crcle square and damond data set respectvely. Sample sze s. We conclude that f the ntal resource state has desred accuracy then ths s also true of all derved resource states. presented protocol s not well suted for an analytcal study of the evoluton of errors on H states. Instead we rely on a numercal study for three dfferent types of errors. We use the trace dstance on states ρ and σ D(ρ σ) = tr( (ρ σ) (ρ σ)) to measure the accuracy of the mperfect H states. We assume Clfford operatons are perfect and that errors can only occur on the H states. We consder three types of erroneous states. Frst we assume that the mxed state ρ a s perfectly along the lne jonng the center of the Bloch sphere and the the perfect state.e. ρ a (p) = ( p) H H + p H H where H = sn π 8 cos π 8 s the state orthogonal to H. We denote the mperfect verson of H obtaned from ρ a states as ρ a. If Clfford operatons are perfect we can always brng any mxed state nto ths form usng twrlng []. However for the protocol to be of practcal nterest we requre t to reman stable under the two followng types of errors where we assume that the state s pure but that the rotaton s slghtly off of the desred axs by δ: ρ b (δ) = ( ( π ) ( π ) ) I + sn 4 + δ X + cos 4 + δ Z ρ c (δ) = (I + sn π 4 cos δx + sn π 4 sn δy + cos π ) 4 Z. We numercally generated pseudo-random nstances of the scheme to produce H states for dfferent values of and for dfferent nose strengths. We consdered nstances for each of the three types of errors and nose strengths 4 6 and 8. Fgures 5 and 6 show that the protocol actually reduces the ampltude of possble errors on the resource H states such that f we. (a)fts: (.7 3 ).3 (.3 5 ).9 and (7.5 8 ) (b)fts: (8.8 4 ).3 (7.3 6 ).3 and (5.3 8 ).5 FIG. 6. Dstances between the deal H states and the mperfects states ρ b and ρ c respectvely. Sample sze s. ɛ = n both cases. start wth H meetng our target accuracy all the subsequent derved H states wll also meet t. Ths even suggests that for bgger values of one could use noser H states and stll acheve the desred accuracy. Ths could make a dramatc dfference f t enables one to reduce the number of dstllaton recursons necessary to prepare the H states. We note a very smlar behavor for the three types of errors. The exponental decay of the dstance between erroneous and deal states confrms that the errors are well behaved under the proposed protocol. We note that the bases for the exponental decay of the errors are comparable but smaller than the bass for the exponental decay of the angle mplemented. So for a gven error rate there exts a pont where the angle of rotaton mplemented by H for some s gong to be comparable or smaller to the error on that angle. However ths s not a problem n practce snce for e.g. ɛ 4 we fnd = 5. For ths value of the angle s θ 57. V. COMPARISON TO SOLOVAY-KITAEV DECOMPOSITION In ths secton we compare the performance of the Solovay-Ktaev decomposton of [6] and that of the

9 9 C Z SK C Z Z On C Off C Z SK C 'Z 'Z On C Off (a) ln(csk Z ) = ln(ln(/ɛ)) ln(con Z ) = ln(ln(/ɛ)) ln(coff Z ) = ln(ln(/ɛ)) Ε (a) ln(c SK Z ) = ln(ln(/ɛ)) ln(c On Z ) = ln(ln(/ɛ)) ln(c Off Z ) = ln(ln(/ɛ)) Ε C SK C On C Off ' ' C SK C On C Off (b) ln(c SK ) = ln(ln(/ɛ)) ln(c On ) = ln(ln(/ɛ)) + ln 3 ln(c Off ) = ln(ln(/ɛ)) + ln 3 Ε (b) ln(c SK ) = ln(ln(/ɛ)) ln(c On ) = ln(ln(/ɛ)) + ln 3 ln(c Off ) = ln(ln(/ɛ)) + ln 3 Ε FIG. 7. Full lne: Cost of Solovay Ktaev decomposton (SKD) of (a) random Z-rotatons and (b) random untares as a functon of the precson ɛ. Dotted lne: Offlne costs of (a) random Z-rotatons and (b) random untares. We assume that loosely speakng a random untary s the composton of three random rotatons hence the addtonal factor of three. Dashed lne: Onlne cost of (a) random Z-rotatons and (b) random untares. (a) For practcal values of ɛ the onlne cost s sgnfcantly smaller than SK. The offlne cost s lower than SK when ɛ (b) Agan for practcal values of ɛ the onlne cost s sgnfcantly smaller than SK. The offlne cost s lower than SK when ɛ FIG. 8. Full lne: Cost of Solovay Ktaev decomposton (SKD) of (a) random Z-rotatons and (b) random untares as a functon of the precson ɛ. Dotted lne: Offlne costs usng { H ψ ψ ψ } as ntal resources of (a) random Z-rotatons and (b) random untares. Dashed lne: Onlne cost usng { H ψ ψ ψ } as ntal resources of (a) random Z-rotatons and (b) random untares. (a) For practcal values of ɛ the onlne cost s sgnfcantly smaller than SK. The offlne cost s lower than SK when ɛ (b) Agan for practcal values of ɛ the onlne cost s sgnfcantly smaller than SK. The offlne cost s lower than SK when ɛ.3 6. schemes presented n ths artcle. In order to do ths we frst consder dfferent Z-rotatons of angles π/6 π/8 and π/4 and dfferent accuraces 4 8 and. Table V lsts the expected costs. C on and C off are respectvely the onlne and offlne costs to mplement the gate usng only H states. C on and C off refer to the costs for the verson of the scheme thats uses { H ψ ψ ψ } as ntal resource states. CSK refers to the extrapolated cost to mplement these gates usng the results from [6]. Ths extrapolated cost averages over all untares (c 3.4). Ths s optmstc snce the results of [6] suggest that Z-rotatons are actually harder to mplement (c 4.3). Note that n theory the cost of ths algorthm s O(log c (/ɛ) where c = 3.97 [4]. In all cases the onlne cost s mnmal when our proposed scheme enhanced by { ψ ψ ψ } s used. For rougher precson e.g. 4 the offlne cost mght be such that the total cost s stll mnmal for the Solovay- Ktaev mplementaton of [6]. For fner precson e.g. 8 or our proposed protocol becomes very advantageous snce the cost of the Solovay-Ktaev decomposton becomes prohbtve. We also compare the average behavor of the dfferent schemes as Fgs. 7 and 8 llustrate. We frst start by notng that loosely speakng a random untary s composed of three random rotatons such that the curves presented prevously n Fg. must be shfted by ln 3 n general. Fg. 7 plots the ft for the Solovay-Ktaev decomposton (sold lne) the onlne cost (dashed) and of-

10 θ C ɛ = 4 ɛ = 8 ɛ = π/6 C SK C on C on C off C off π/8 C SK C on C on C off C off π/4 C SK C on C on C off C off TABLE V. C on and C off are respectvely the onlne and offlne costs to mplement the Z-rotaton by angle θ usng only H states to precson ɛ. C on and C off refer to the costs for the verson of the scheme thats uses { H ψ ψ ψ } as ntal resource states. C SK refers to the extrapolated cost to mplement these gates usng the results from [6]. flne cost (dotted). For all practcal accuraces the onlne cost of our proposed scheme s consstently the smallest. However the offlne cost becomes advantageous when ɛ < for Z-rotatons and ɛ <.67 7 for random untares. Fg. 8 plots the same for the scheme wth addtonal ntal resource states. Smlarly the offlne cost becomes advantageous when ɛ < for Z-rotatons and ɛ <.3 6 for random untares. VI. CONCLUSION We have proposed an alternatve protocol to Solovay- Ktaev decomposton that results n sgnfcantly smaller resource costs n both the number of requred resource states and the depth of the crcut. We have shown a sgnfcant mprovement on average n the value of c and n many cases the number of dstlled states and rotatons requred to mplement a sngle-qubt untary gate are reduced. Another advantage of our protocol s that the number of resources requred s a smoother functon of accuracy whereas Solovay-Ktaev decomposton s step-lke n nature because of the recurson process used n practce. However note that our protocols and Solovay-Ktaev decomposton are not exclusve. It mght be that some untares are better mplemented usng Solovay-Ktaev decomposton whle our scheme s better suted for Z-rotatons whch occur among other algorthms n the quantum Fourer transform. As future research there are lkely a varety of other crcuts that enable other ladders of states. One natural extenson would be to use the SH egenstates dstlled usng the protocols of [ 3]. Another extenson would be to perform a systematc study of small Clfford crcuts. Fnally we note that mplementng a rotaton by choosng the state whch results n an angle closest to the target angle s a smple way of achevng our goal but t s surely suboptmal. An mportant research drecton would be to optmze the sequence of angles requred to mplement the desred rotaton. VII. ACKNOWLEDGEMENTS We thank Alex Bocharov for many useful dscussons. [] A. Ktaev Ann. Phys. 33 (3) arxv:quantph/977. [] M. H. Freedman Proc. Natl. Acad. Sc. 95 (998). [3] J. Preskll n Introducton to Quantum Computaton edted by H.-K. Lo S. Popescu and T. Spller (998) arxv:quant-ph/9748. [4] M. H. Freedman M. Larsen and Z. Wang Commun. Math. Phys. 7 () arxv:quant-ph/8. [5] C. Nayak S. Smon A. Stern M. Freedman and S. D. Sarma Rev. Mod. Phys. 8 (8) arxv: [6] H. Bombn and M. A. Martn-Delgado Journal of Physcs A Mathematcal General (9) arxv:74.54 [quant-ph]. [7] A. G. Fowler A. M. Stephens and P. Groszkowsk Phys. Rev. A 8 (9). [8] D. Gottesman The Hesenberg Representaton of Quantum Computers Ph.D. thess Caltech (998) arxv:quant-ph/9876v. [9] S. Aaronson and D. Gottesman Phys. Rev. A (4) arxv:quant-ph/4696. [] M.A.Nelsen and I.L.Chuang Quantum Computaton and Quantum Informaton (Cambrdge Unversty Press Cambrdge UK ). [] S. Bravy and A. Ktaev Phys. Rev. A 7 36 (5). [] A. M. Meer B. Eastn and E. Knll arxv:4.4v (). [3] S. Bravy and J. Haah () [4] C. M. Dawson and M. A. Nelsen (5). [5] A. Ktaev et al. Classcal and Quantum Computaton (Amercan Mathematcal Socety Provdence RI ). [6] A. Bocharov and K. M. Svore arxv:6.33v ().