arxiv:quant-ph/ v1 10 Oct 2002

Transcription

1 A Quatum Radom Walk Search Algorithm Neil Shevi 1, Julia Kempe 1,,3, ad K. Birgitta Whaley 1 Departmets of Chemistry 1 ad Computer Sciece, Uiversity of Califoria, Berkeley, CA CNRS-LRI, UMR 863, Uiversité de Paris-Sud, Orsay, Frace Dated: February 1, 008 arxiv:quat-ph/010064v1 10 Oct 00 Quatum radom walks o graphs have bee show to display may iterestig properties, icludig expoetially fast hittig times whe compared with their classical couterparts. However, it is still uclear how to use these ovel properties to gai a algorithmic speed-up over classical algorithms. I this paper, we preset a quatum search algorithm based o the quatum radom walk architecture that provides such a speed-up. It will be show that this algorithm performs a oracle search o a database of N items with O N calls to the oracle, yieldig a speed-up similar to other quatum search algorithms. It appears that the quatum radom walk formulatio has cosiderable flexibility, presetig iterestig opportuities for developmet of other, possibly ovel quatum algorithms. I. INTRODUCTION Recet studies of quatum radom walks have suggested that they may display differet behavior tha their classical couterparts [1,, 3, 4, 5]. Oe of the promisig features of these quatum radom walks is that they provide a ituitive framework o which to build ovel quatum algorithms. Sice may classical algorithms ca be formulated i terms of radom walks, it is hoped that some of these may be traslated ito quatum algorithms which ru faster tha their classical couterparts. However, previous to a very recet paper by Childs et. al. [6], there had bee o quatum algorithms based o the radom walk model. I this paper we show that a quatum search algorithm ca be derived from a certai kid of quatum radom walk. Optimal quatum search algorithms are already well kow[7, 8, 9]. The search algorithm from a quatum radom walk we preset here shows some differeces from the established search algorithms ad may possess useful properties with respect to robustess to oise ad ease of physical implemetatio. It also provides a ew directio for desig of quatum algorithms from radom walks, which may evetually lead to etirely ew algorithms. Curret research uses two distict models for quatum radom walks, based o either discrete time steps or o cotiuous time evolutio. Discrete time quatum radom walks were itroduced as a possible ew tool for quatum algorithms geeralizig discrete classical Markov chais []. The discrete time walk ca be thought of as a successio of uitary operatios, each of which has a o-zero trasitio amplitude oly betwee eighborig odes of the graph. The relatio of these to classical Markov chais provides cosiderable motivatio for exploratio of discrete radom walks. Withi the field of classical algorithms, the applicatio of classical Markov chais i classical algorithms has bee quite revolutioary, providig ew approximatio ad optimizatio algorithms. By aalogy, it might reasoably be hoped that similar algorithmic advaces could be obtaied for quatum algorithms from developmet of the quatum radom walks. The secod quatum radom walk model is the cotiuous-time quatum radom walk, itroduced i [4, 5, 6]. I the cotiuous-time walk the adacecy matrix of the graph is used to costruct a Hamiltoia which gives rise to a cotiuous time evolutio. This model differs from the discrete time walk i that eve for small times there is a expoetially small probability of trasitio to o-adacet odes. I this paper, we will cosider the discrete-time model oly. The paper is orgaized as follows. Sec. II provides a brief itroductio to discrete-time quatum radom walks. Sec. III describes the radom walk search algorithm ad provides a proof of its correctess. Sec. IV summarizes the similarities ad differeces betwee the radom walk search algorithm ad Grover s search algorithm. Coclusios are preseted i Sec. V. Notatio: Followig stadard computer sciece otatio we will use the followig to characterize the growth of certai fuctios: We will say f = Og if there are positive costats c ad k such that 0 f cg for k. Similarly f = Ωg if 0 cg f for costats c, k 0 ad k. II. BACKGROUND The discrete time radom walk ca be described by the repeated applicatio of a uitary evolutio operator U. This operator acts o a Hilbert space H C H S where H C is the Hilbert space associated with a quatum coi ad H S is the Hilbert space associated with the odes of the graph. The operator U ca be writte as [] U = S C 1

2 where S is a permutatio matrix which performs a cotrolled shift based o the state of the quatum coi, ad C is a uitary matrix which flips the quatum coi. This operatio ca be visualized by aalogy to a classical radom walk. For each iteratio of a discrete time classical radom walk o a graph, a coi is flipped. The walker the moves to a adacet ode specified by the outcome of the coi flip. A equivalet process occurs i the quatum radom walk, with the modificatio that the coi is a quatum coi ad ca therefore exist i a superpositio of states. This modificatio ca lead to dramatic differeces i behavior betwee the classical ad quatum radom walks. However, it should be oted that if the state of the coi is measured after each flip, the the quatum radom walk reverts to a classical radom walk ad similarly if the state of the odes is measured after every step. A importat feature of the discrete time quatum radom walk that has sigificace for its use i developmet of quatum algorithms, is that by virtue of its defiitio this walk will be efficietly implemetable o a quatum computer wheever its classical couterpart is. By efficiet we mea that the walk ca be simulated by a circuit with a umber of gates that is polyomial i the umber of bits qubits. This is due to the very similar structure of both these walks. To illustrate this, assume we have a efficiet way to implemet the classical radom walk o the uderlyig graph, i.e. to perform the coi-flip ad subsequet shift. The shift is coditioal o the outcome of the coi-flip which determies the directio of the ext step, i.e. we have a classical efficiet circuit which performs a cotrolled shift o the basis states. It is straightforward [10] to traslate this circuit ito a quatum circuit which performs the uitary cotrolled shift of Eq. 1. Similarly if there is a efficiet procedure to flip the classical coi of the radom walk, there will be a efficiet way to implemet a quatum coi. Hece implemetatio of the discrete time radom walk is automatically efficiet if the uderlyig classical walk is efficietly implemetable. Note that if o measuremet is made, the quatum walk is cotrolled by a uitary operator rather tha a stochastic oe. This implies that there is o limitig statioary distributio [, 11]. Nevertheless, several recet works have show that cosistet otios of mixig time ca be formulated, ad have show polyomial speed up i these quatum mixig times relative to the classical aalog [, 11]. Aother quatity for which quatum walks have show speed-up relative to their classical aalogs is the hittig time [1, 13]. Uder certai coditios this speed-up ca be expoetial compared to the classical aalogue. We refer the reader to the recet papers [], [3], [11] ad [1] for some results obtaied from discrete time quatum radom walks. Our radom walk search algorithm will be based o a radom walk o the -cube, i.e. the hypercube of dimesio [11, 1]. The hypercube is a graph with N = odes, each of which ca be labelled by a -bit biary strig. Two odes o the hypercube described by bitstrigs x ad y are are coected by a edge if x y = 1, where x is the Hammig weight of x. I other words, if x ad y differ by oly a sigle bit flip, the the two correspodig odes o the graph are coected. Thus, each of the odes o the -cube has degree i.e. it is coected to other odes, so the Hilbert space of the algorithm is H = H H. Each state i H ca be described by a bit strig x, which specifies the positio o the hypercube, ad a directio d, which specifies the state of the coi. The shift operator, S, maps a state d, x oto the state d, x e d, where e d is the d th basis vector o the hypercube. S ca be writte explicitly as, S = d, x e d d, x d=0 x To completely specify the uitary evolutio operator U, the coi operator, C, must also be chose. Normally, the coi operator is chose such that the same coi is applied to each ode o the graph. This is the case i previous studies of discrete quatum walks o the lie [, 3, 14] ad o the hypercube [11, 1]. I other words, the coi operator C ca be writte as, C = C 0 I 3 where C 0 is a uitary operator actig o H C. If C is separable accordig to Eq. 3 the the eigestates of U are simply the tesor product of the eigestates of a modified coi C k ad of the Fourier modes of the hypercube labelled by -bit strigs k [11]. Oe frequetly chose separable coi is Grover s diffusio operator o the coi space, give by C 0 = G = I + s C s C, 4 where s C is the equal superpositio over all directios, i.e. s C = 1 d=1 d [11]. This coi operator is ivariat to all permutatios of the directios, so it preserves the permutatio symmetry of the hypercube. The use of the Grover diffusio operator as a coi for the hypercube was proposed i [11], where it was poited out that this operator is the permutatio ivariat operator farthest away from the idetity operator [11]. So, heuristically, it should provide the most efficiet mixig over states, from ay give iitial state. The o-trivial eigevalues ad

3 3 eigevectors of U are give by [11], e ±iω k = 1 k ± i k k 5 v k, v k = 1 k x / { d, x 1/ k if kd = 1 i/ k if k d = 0 x,d 6 Note that the equal superpositio over all states, ψ 0 = s C s S, where s S is the equal superpositio over the odes, is a eigevector of U with eigevalue 1. So repeated applicatio of U leaves the state ψ 0 uchaged. I order to create a search algorithm usig the quatum radom walk architecture, we ow cosider a small perturbatio of the uitary operator U. Specifically, we cosider markig a sigle arbitrary ode by applyig a special coi to that ode. I this respect, the coi operator ow takes o the fuctio of a oracle. Yet, cotrary to the stadard oracle which flips the phase of the marked state i the stadard settig of a search algorithm, the oracle acts istead by applyig a markig coi, C 1, to the marked ode ad a differet coi, C 0, to the umarked odes. Without loss of geerality, we ca assume that the marked ode correspods to the all-zero strig x target = 0. The our coi operator becomes C = C 0 I + C 1 C The markig coi, C 1, ca be ay uitary matrix. For simplicity, we will cosider here the case where C 1 = I. The our perturbed uitary evolutio operator, U, is give by U = S C = S G I G I 0 0 = U S s C s C Aalysis of the effects of this perturbatio leads directly to the defiitio of the radom walk search algorithm, as will be described i the ext sectio. III. RANDOM WALK SEARCH ALGORITHM A. Overview of the Algorithm We defie the search space of the algorithm to be the set of all -bit biary strigs, x = {0, 1}. We cosider the fuctio f x = {0, 1}, such that f x = 1 for exactly oe iput x target. Our goal is to fid x target. Usig the mappig of -bit biary strig to odes o the hypercube, this search problem is the equivalet to searchig for a sigle marked ode amogst the N = odes o the -cube. For purposes of the proof, we have set the marked ode to be x target = 0, but the locatio of the marked ode has o sigificace. The radom walk search algorithm is implemeted as follows: 1. Iitialize the quatum computer to the equal superpositio over all states, ψ 0 = s C s S. This ca be accomplished efficietly o the ode-space by applyig sigle-bit Hadamard operatios to the 0 state. A similar procedure works for the directio space.. Give a coi oracle, C, which applies the coi C 0 = G to the umarked states ad the coi C 1 = I to the marked state, apply the perturbed evolutio operator, U = S C, t f = π times. 3. Measure the state of the computer i the d, x basis. It is our claim that with probability 1 O 1/, the outcome of the measuremet will be the marked state. By repeatig the algorithm a costat umber of times, we ca determie the marked state with a arbitrarily small degree of error. I the remaider of this sectio we provide a proof of this algorithm. The geeral outlie of the proof that we will preset is the followig. We eed to determie the result of the operatio U t o the iitial state ψ 0. To do this, we will first simplify the problem by showig that the perturbed walk o the hypercube ca be collapsed to a walk o the lie Theorem 1. Next, by costructig two approximate

4 eigevectors of U, ψ 0 ad ψ 1, we will show that there are exactly two eigevalues of U that are relevat i.e. the iitial state ψ 0 has high overlap with the space spaed by the correspodig eigevectors, see Theorem ad Theorem 3. We deote these eigevalues by e iω 0 ad e iω 0. We will the show that the correspodig eigevectors ω 0 ad ω 0 ca be well-approximated by liear combiatios of the iitial state ψ 0 ad the secod state ψ 1 Theorem 4. As a result, our radom walk search algorithm ca be approximated by a two-dimesioal rotatio i the ω 0, ω 0 -plae away from the iitial state ψ 0 1/ ω 0 + ω 0 ad towards ψ 1 i/ ω 0 + ω 0, which costitutes a very close approximatio to the target state x target. Fially, we show that each applicatio of the evolutio operator U correspods to a rotatio agle of approximately 1/ Theorem 5. Hece, the N search is completed after approximately π steps, i.e., after O calls to the oracle, where N = is the umber of odes. 4 B. Proof of Correctess I geeral, aalytic determiatio of the eigespectrum of a large matrix is a dautig task, so we will take advatage of the symmetries iheret i U to simplify the problem. Let us first show that the perturbed radom walk o the hypercube ca be collapsed oto a radom walk o the lie. Let P i be the permutatio operator which swaps the bits i ad, i both the ode space ad the coi space. I other words, give a state d, x, uder the permutatio operator, P i, the i th ad th bits of x are swapped ad the directios d = i ad d = are swapped. Clearly the uperturbed evolutio operator U commutes with P i sice every directio i the uperturbed walk is equivalet. Theorem 1 U commutes with P i. Proof. s P i U P i = P i UP i P i S C s C 0 0 P i = U d=0 P i d, e d d, 0 P i 9 = U So, [U, P i ] = 0. Because the iitial state ψ 0 is a eigevector of P i with eigevalue 1 for all i ad, ad [U, P i ] = 0, ay itermediate state ψ t = U t ψ 0 must also be a eigevector of eigevalue 1 with respect to P i. Thus, U t preserves the symmetry of ψ 0 with respect to bit swaps. It is therefore useful to defie basis states, R, 0, L, 1, R, 1,..., R, 1, L, where 1 R, x = x d, x 10 x x =x x d =1 x =x x d =0 1 L, x = x d, x 11 x which are also ivariat to bit swaps P i. These states spa the eigespace of eigevalue 1 of P i. Usig these basis states, we ca proect out all but oe spatial degree of freedom ad effectively reduce the radom walk o the hypercube to a radom walk o the lie. This is illustrated i Fig. 1. The marked ode correspods ow to R, 0. We ca rewrite U, U, ad ψ 0 i this collapsed basis. First ote that the shift operator S i this basis acts as S = R, x L, x L, x + 1 R, x 1 x=0 ad the uperturbed coi acts as C 0 = x=0 cosωx siω x si ω x cosω x x x 13

5 where cosω x = 1 x ad siω x = x x ad where the first part acts o the space spaed by { R, L } ad the secod part acts o the positios { 0,..., } o the lie. Note that the coi of the collapsed walk is ot homogeeous i space ay more. The uitary operator U o the restricted space acts as U = x=0 R, x cosω x+1 L, x siω x+1 R, x x=1 L, x siω. 14 x 1 L, x 1 + cosω x 1 R, x 1 5 Similarly, U = U + U = U L, 1 R, 0 15 Note that the oly differece betwee U ad U is i the sig of the matrix elemet i positio L, 1, R, 0. Fially, ψ 0 = 1 1 R, 0 + L, + x=1 x 1 L, x + x R, x 16 Sice U ad P i are mutually diagoalizable, the eigevectors of U i the reduced space are also bit-flip ivariat. Examiig Eq. 5, it is clear that if we take the equal superpositios of all eigevectors of same eigevalue v k such that k = k, the resultig eigevector will be bit-swap ivariat. Thus we defie, ω k = 1 k k =k v k 17 which are the eigevectors of U with eigevalues e iω k i the collapsed symmetric space. Note that both U ad U are represeted by real matrices; therefore, their eigevalues ad eigevectors will come i complex cougate pairs. Havig determied these geeral properties of the perturbed matrix U, we ow tur to the problem of aalyzig the eigevalue spectrum of U. Let A be the arc o the uit circle cotaiig all complex umbers of uit orm with real part greater tha 1 /3. I other words, { A = z Re z > 1 } 3, z = 1 18 Fig. shows the geometrical represetatio of A together with the eigevalue spectra of the uperturbed ad perturbed matrices for = 8. We will prove that A cotais exactly two eigevalues e iω 0 ad e iω 0 of U. First, we will prove that there are at most two eigevalues with real part greater tha 1 /3. The we will show that there are at least two eigevalues o A. From these facts, it follows that there are exactly two eigevalues of U o A. Theorem There are at most two eigevalues of U with real part greater tha 1 /3. Proof. We will prove by cotradictio. Let us assume that there are three eigevalues, e iω 0,e iω 1, ad e iω, with real part greater tha 1 /3. Let ω 0, ω 1, ad ω be the correspodig eigevectors. The, Re i ω i U ω i = Re i eiω i ω i ω i > 3 / 19 Let us defie Ω to be the subspace spaed by ω 0, ω 1, ω. The we ca write Eq. 19 as the partial trace of U over Ω, Re Tr Ω U > 3 / 0 Let us ow defie ψ = 1 0, R 1, L. We ca expad the ω 0, ω 1, ad ω i terms of ψ 0, ψ ad a residual vector, ω 0 = c 00 ψ 0 + c 01 ψ + c 0 r 0 ω 1 = c 10 ψ 0 + c 11 ψ + c 11 r 1 ω = c 0 ψ 0 + c 1 ψ + c r 1

6 6 where r i is a ormalized vector orthogoal to ψ 0 ad ψ. We ow observe that, due to the basis ivariace of the trace, Eq. 0 holds for ay liear combiatio of ω 0, ω 1, ad ω. Thus, we ca costruct three ew orthoormal vectors, α 0, α 1, ad α by takig liear combiatios of ω 0, ω 1, ad ω such that, α ψ 0 = α ψ = 0 I other words, we ca expad α 0, α 1, ad α as, α 0 = c 00 ψ 0 +c 01 ψ +c 0 r 0 α 1 = c 10 ψ 0 +c 11 ψ +c 1 r 1 α = r 3 Sice α 0, α 1, ad α still form a basis for Ω, from Eq. 0 it follows that, 3 / < Re i α i U α i 4 Sice U is a uitary operator, we kow that Re α i U α i 1 for all α i. Thus, applyig this iequality to the first two terms i the sum, we obtai, Re i α i U α i + Re α U α 5 Sice, U = U + U, we ca write Re α U α = Re α U α + Re α U α 6 Let us first cosider α U α. We ca expad α i terms of the uperturbed eigestates, α = b ω. So, Re α U α = b cosω. However, sice α ψ 0 = 0, there is o cotributio from the eigevalue with value 1. The eigevalue with the ext- largest real part is e iω1 = 1 / + i 1. Thus, Re α U α 1 / 7 Next, we cosider, α U α. Let ψ + = 1 0, R + 1, L. Usig Eq. 15 we ca express U i terms of ψ ad ψ +, U = ψ ψ + ψ ψ + ψ + ψ ψ + ψ + 8 But sice α ψ = 0 see Eq., α U α = ψ + α 0 9 The, Re α U α 1 /. Combiig Eq. 5, Eq. 7, ad Eq. 9, we obtai, Re i α i U α i 3 /. 30 Sice this cotradicts Eq. 4, our assumptio must be false. Theorem 3 There are at least two eigevalues of U o A. Proof. We will costruct two approximate eigevectors of U, ψ 0 ad ψ 1. ψ 0 is give by Eq. 16. Usig Eq. 15, Ad, U ψ 0 = ψ 0 / L, 1 31 ψ 0 U ψ 0 = ψ 0 ψ 0 ψ 0 L, 1 R, 0 ψ 0 = 1 1/ 3

7 So, apart from a small residual, ψ 0 is also almost a eigevector of U with eigevalue 1. Now, we eed to fid a secod approximate eigevector, ψ 1. Let / 1 ψ 1 = 1 x=0 R, x 1 x L, x + 1 /c 33 x 7 where c is a ormalizatio costat, c = / 1 x=0 1 x 34 Usig this defiitio ad Eq. 14, we see that So, U ψ 1 = ψ 1 1 c / R, / 1 + L, / ψ 1 U 1 ψ 1 = 1 c 36 / It is straightforward to show that 1 < c < 1 + / for sufficietly large. Thus, except for a small residual, ψ 1 is almost a eigevector of U with eigevalue 1. Now let us verify that there is at least oe eigevalue of U o A. Let us assume that there are o eigevalues of U o A. The cosω < 1 /3 for all. The usig Eq. 3, 1 1/ = Re ψ 0 U ψ 0 = ψ 0 ω cosω < 1 /3 ψ 0 ω 37 = 1 /3 which is wrog for > 3. Hece our assumptio is false ad there must be at least oe eigevalue of U o A. Now let us assume there is exactly oe eigevalue of U, e iω 0, o A. The, 1 1 = Re ψ 0 U ψ 0 = ψ 0 ω cosω = ψ 0 ω 0 cosω ψ 0 ω cosω ψ 0 ω ψ 0 ω 0 1 /3 38 Rearragig terms, ψ 0 ω If we use ψ 1 as a trial vector ad follow the same argumets, we obtai the iequality But sice ψ 0 ad ψ 1 are orthoormal, this leads to a cotradictio, sice, ψ 1 ω c. 40 / 1 = ω 0 ω 0 ψ 0 ω 0 + ψ 1 ω c / 41

8 which is ot true for large. Hece, there must be at least two eigevalues o the arc A. As oted above, the eigevalues ad eigevectors of U come i complex cougate pairs. I particular the two eigevalues o A must be a complex cougate pair; let e ±iω 0 be the two eigevalues o A. The correspodig eigevectors obey ω 0 = ω 0 if e iω 0 = e iω 0 = 1, the we ca costruct liear combiatios of ω 0 ad ω 0 for which this statemet is true. We will ow show that ±ω 0 ca be well-approximated by liear combiatios of ψ 0 ad ψ 1. Theorem 4 ω 0 = p 0 ψ 0 + p 1 e iη ψ p 0 p 1 r 0 ω 0 = p 0 ψ 0 + p 1 e iη ψ p 0 p 1 r 0 4 where p 0 = ω 0 ψ 0 = ω 0 ψ 0, p 1 = ω 0 ψ 1 = ω 0 ψ 1 ad r 0 is a ormalized vector orthogoal to ψ 0 ad ψ 1. Furthermore, 1/ p 0 1/ 3/ +1 3 ad 1/ p 1 1/ 8c /, with eiη = i +, where = O. / Proof. Sice ψ 0 ad ψ 1 are real vectors, ω 0 ψ 0 = ω 0 ψ 0 1/ ad ω 0 ψ 1 = ω 0 ψ 1 1/. Usig Eq. 3, Rearragig terms, we obtai, 1 1 = Re ψ 0 U ψ 0 = cosω ω ψ 0 = p 0 cosω cosω < p /31 p 0 p 0 1/ Usig Eq. 36 ad the same argumets as above we obtai Up to a global phase ω 0 ca be writte as, ω ψ p 1 1/ 8c. 45 / ω 0 = ω 0 ψ 0 ψ 0 + ω 0 ψ 1 e iη ψ p 0 p 1 r 0 46 which yields Eq. 4 for ω 0 ad ω 0. To estimate e iη ote that sice ω 0 ad ω 0 are eigevectors of a uitary matrix, they must be orthogoal. Cosequetly, Solvig for e iη, we obtai, 0 = ω 0 ω 0 = p 0 + p 1 e iη + 1 p0 p 1 r 0 r 0 47 Re e iη = p 0 1 p 0 p 1 Re r 0 r 0 p 1 48 Assume Re r 0 r 0 0. The, usig 1/1 x 1 + x for small x we get p 0 p 1 Re e iη p c c / p 1 / 3 4 8c / 49

9 which i tur implies that e iη = i + with = O /. A similar reasoig holds if Re r 0 r 0 0. As a last igrediet we eed to boud the agle ω 0. These bouds are provided i the followig fial theorem. Theorem 5 1 c β ω 0 1 c + β, where β = O 3/ 9 Proof. We will approximate e iω 0 = ω 0 U ω 0 by α U α, where α = 1/ ψ 0 + e iη ψ 1. Let us first evaluate e iω 0 α U α. We ca expad U i terms of its eigevectors to obtai, We the ote that from Eq. 4, e iω 0 α U α = eiω 0 ω α e iω 50 ω 0 α = p 0 / + p 1 / 51 So, e iω 0 α U α = eiω e iω 0 p0 0 / + p 1 / e iω p 0 / p 1 / + 1 p 0 / p 1 / c / ω ω 0 ω ω 0 ω α e iω ω α 5 with 1 x 1 x for 0 x 1. Usig the fact that the biomial coefficiets approach the Gaussia distributio for large, such that, x / = x π e / 53 we ca rewrite Eq. 5 takig the leadig order terms i. Recallig that c > 1, we obtai, e iω 0 α U 3/ α = O 54 Eq. 54 is a explicit formula which bouds the distace i the complex plae betwee the eigevalue of iterest, e iω 0, ad the matrix elemet α U α. Fig. 3 shows the geometric represetatio of Eq. 54. Note that, siω 0 Im α U α = Im e iω 0 α U α e iω 0 α U α 55 Next, we evaluate Im α U α usig Eqs. 31 ad 35, Im α U α = Im e iη ψ 0 U ψ 1 e iη ψ 1 U ψ 0 = Im 1 eiη c / ψ 0 R, / 1 + ψ 0 L, / + 1 e iη = Im eiη c = 1 c O The, usig Theorem 4, Eq. 54, ad Eq. 56 we ca write si 1 ω 0 + c + O 3/ = O / / ψ 1 L,

10 Usig si x = x + Ox 3 ad keepig oly leadig order terms solvig for ω 0 gives us, 1 c 3/ O ω 0 1 c 3/ + O. 58 We ca ow quatitatively describe the overall operatio of the algorithm. Startig with iitial state ψ 0, we cosider the state of the computer after t applicatios of U. The usig Theorem 4 we ca expad ψ 0 as ψ 0 = p 0 ω 0 + ω 0 + δ r 59 where δ = 1 p 0 = O / ad r is a residual ormalized vector orthogoal to ω 0 ad ω 0. Now U t ψ 0 = p 0 e iω 0 t ω 0 + e iω 0 t ω 0 + δ r = p 0 cosω 0 t ψ 0 p 0 p 1 si ω 0 t + Re eiω 0 t ψ p 0 p 1 e iω 0 t r 0 + e iω 0 t r δ r = cosω 0 t ψ 0 si ω 0 t ψ 1 + O 3/4 r where r is some residual ormalized vector ot ecessarily orthogoal to ψ 0 ad ψ 1. Startig with ψ 0 ad applyig U for t f = π ω 0 steps, we approximately rotate from ψ 0 to ψ 1. From ψ 1 we ca obtai x target = 0 with high probability p, sice from Eq. 33 ad with 1 + 1/ c 1 + / for large p = 1 c d d, 0 ψ 1 = 1 c R, 0 ψ 1 1/ 1+/ = 1 O 1/ 61 Fially, to obtai t f i terms of, we make use of the bouds o ω 0 provided by Theorem 5: 10 t f = πc 1 ± O 3/ = π 1 + O 1 6 If we set the umber of time steps to be t f = π or the closest iteger the si ω 0t f = si π 1 O 1 = 1 O So the probability to measure x target after t f = π steps is still p success = 1/ O1/. Hece, by repeatig the algorithm a costat umber of times, the probability of error ca be made arbitrarily small. Note the periodic ature of the evolutio uder U Eq. 60; this meas that if we measure at t > t f the probability of success will decrease ad later icrease agai. I summary we arrived at the fial result that the marked state is idetified after O N calls to the oracle. IV. CONNECTION TO GROVER S ALGORITHM The operatio of the radom walk search algorithm is similar i may ways to the operatio of Grover s search algorithm. Both algorithms begi i the equal superpositio state over all bit strigs. Both algorithms make use of the Grover diffusio operator, G, sometimes kow as the Grover iterate. Both algorithms ca be viewed as a rotatio i a two-dimesioal subspace. Both algorithms use a oracle which marks the target state with a phase N of 1. Both algorithms have a ruig time of O. I both algorithms we have to measure at a specific time to obtai maximum probability of success. However, there are several importat differeces betwee the two search algorithms. I this sectio, we call attetio to the ways i which the radom walk search algorithm is distict from Grover s algorithm ad cosider how these differeces affect performace ad implemetatio. It is well-kow that Grover s algorithm ca be mapped exactly oto a rotatio i the two-dimesioal subspace spaed by the equal-superpositio state ψ 0 ad the marked state 0 [7]. Each iteratio i Grover s algorithm correspods to a rotatio i this subspace. I this paper, we have show that the radom walk search algorithm

11 11 ca also be viewed as a rotatio i a two-dimesioal subspace. However, there are two importat distictios. First, the radom walk search algorithm ca oly be approximately mapped oto a two-dimesioal subspace. Ulike Grover s algorithm, this mappig is ot exact. Secod, the two-dimesioal subspace i which the radom walk search algorithm is approximately cotaied is spaed by ψ 0 ad ψ 1, ot by ψ 0 ad 0. Hece, the fial state of the algorithm is ot exactly the pure marked state, 0, as it is i Grover s algorithm. It is a liear combiatio of states which is composed primarily of the marked state, but also possesses small cotributios from its earest eighbors, secod-earest eighbors, etc. Thus, the radom walk search algorithm cotais traces of the uderlyig topology of the hypercube o which it is based. Aother differece betwee the two algorithms is their use of the Grover diffusio operator, G. I Grover s algorithm, this operator is applied to the etire -dimesioal search space correspodig to the ode space i the radom walk search algorithm. O the other had, Grover s diffusio operator, G, i the radom walk algorithm is used as the quatum coi, ad acts oly o the -dimesioal coi space. This fact may be of practical use for certai physical implemetatios sice may physical implemetatios of quatum computers cotai multiple types of qubits, which have differet atural gate sets. We could exploit this variety usig the radom walk search algorithm by choosig the coi space to be represeted by qubits o which it is coveiet to implemet the Grover diffusio operator. Aother similarity betwee the two algorithms is the implemetatio of the oracle. I Grover s algorithm, the oracle marks the target state with a phase of 1. To arrive at this radom walk search algorithm, we chose the markig coi C 1 to be the I coi. This choice was actually motivated because it yielded a result that was ameable to aalysis, ad while the emergece of Grovers algorithm appears atural i hidsight, it was ot obvious at the outset. However, more geerally, it is ot clear whether this choice of marked coi is either optimal or uique. I fact, umerical simulatios have show us that may differet types of markig cois will yield search algorithms. Ufortuately, aalytic treatmet of the quatum radom walk for more complicated cois has prove substatially more difficult tha the istace aalyzed here for C 1 = I. It is a ope questio what costat factor gais might be made by usig differet markig cois to implemet the search. V. CONCLUSIONS I this paper, we have show that the radom walk search algorithm ca search a list of items i time proportioal to N. The lower boud o a quatum search of a N-item list is kow to be Ω [9]. Thus, up to a costat factor, the radom walk search algorithm is optimal. However, although after repetitio of the algorithm a costat umber of times the result is arbitrarily close to the result of Grover s search, the radom walk search algorithm is ot idetically equivalet to Grover s algorithm. I particular, the fial solutio obtaied by the radom walk search still retais some of the uderlyig character of the hypercube o which it was based, with a small admixture of states other tha the solutio at the marked ode. The radom walk search aalyzed here was based o a discrete walk o the hypercube. I geeral a similar methodology ca be applied to ay regular graph, e.g., a two-dimesioal hexagoal lattice with periodic boudary coditios, a three-dimesioal rectagular lattice with periodic boudary coditios, etc. We have umerical evidece idicatig that this methodology will yield quatum search algorithms whe applied to other regular -dimesioal lattices. Future studies will ivestigate the extet of optimality of such search algorithms. The itriguig possibility of fidig ovel algorithms based o the radom walk also remais a ope questio. The results described here idicate that the radom walk search algorithm provides a suggestive framework for ew algorithms. Though the optimality of Grover s algorithm precludes the costructio of a improved oracle-based search algorithm based o a quatum walk, evertheless, may other oracle problems still exist for which a quatum walk may be advatageous. For istace, the lower boud o quatum search holds oly for oracles which provide yes/o iformatio [9]. Our choice of markig coi here has a clear relatio to a idetifiable compoet of Grover s algorithm. I geeral, the markig coi ca be a arbitrary uitary matrix. The markig coi provides a ituitive meas by which to itroduce a large amout of iformatio to a oracle problem. Thus, it is possible that uique cois with iterestig properties may give rise to a etirely ew algorithm. Overall we coclude that the quatum radom walk provides a meas for isight ito existig quatum algorithms ad offers a potetially vast source for developmet of ew algorithms.

12 1 FIG. 1: Collapsig a radom walk o the hypercube to a radom walk o the lie. States o the hypercube are mapped to state o the lie based o their Hammig weight ad the directio i which they poit see text. Ackowledgmets NS thaks the Uiversity of Califoria, Berkeley, for a Berkeley Fellowship. This effort is sposored by the Defese Advaced Research Proects Agecy DARPA ad the Air Force Laboratory, Air Force Material Commad, USAF, uder agreemet umber F We also thak NSF ITR/SY award [1] Y. Aharoov, L. Davidovich, ad N. Zagury, Phys. Rev. A 48, [] D. Aharoov, A. Ambaiis, J. Kempe, ad U. Vazirai, i Proceedigs of ACM Symposium o Theory of Computig STOC ACM, New York, NY, 001, pp , LANL preprit quat-ph/ [3] A. Ambaiis, E. Back, A. Nayak, A. Vishwaath, ad J. Watrous, i Proc. 33rd STOC ACM, New York, NY, 001, pp [4] E. Farhi ad S. Gutma, Phys. Rev. A 58, [5] A. Childs, E. Farhi, ad S. Gutma, quat-ph/ [6] A. Childs, R. Cleve, E. Deotto, E. Farhi, S. Gutma, ad D. Spielma, quat-ph/ [7] L. Grover, i Proc. 8 th Aual ACM Symposium o the Theory of Computatio ACM Press, New York, NY, 1996, pp [8] L. Grover, Phys. Rev. Lett. 79, [9] C. Beett, E. Berstei, G. Brassard, ad U. Vazirai, SIAM J. Comput. 6, [10] M. Nielse ad I. Chuag, Quatum Computatio ad Quatum Iformatio Cambridge Uiversity Press, Cambridge, 000. [11] C. Moore ad A. Russell, Proc. RANDOM, to appear 00, lal-arxive quat-ph/ [1] J. Kempe, quat-ph/ [13] T. Yamasaki, H. Kobayashi, ad H. Imai, quat-ph/ [14] J. Watrous, Joural of Computer ad System Scieces 6,

13 13 FIG. : The results of umerical spectral aalysis of U ad U for = 8. Circles idicate eigevalues of U. Crosses idicate eigevalues of U. FIG. 3: Geometric represetatio of Theorem 5, which proves that the eigevalue, e iω 0, must be located o a disc of radius 8c / cetered at α U α. The positio of the eigevalue is deoted by a cross.