arxiv:quant-ph/ v1 21 Jan 1998

Transcription

1 Noliear quatum mechaics implies polyomial-time solutio for NP-complete ad #P problems arxiv:quat-ph/98004v Ja 998 Daiel S. Abrams Departmet of Physics, MIT -8b Cambridge, MA 039 Seth Lloyd d Arbeloff Laboratory for Iformatio Scieces ad Techology Departmet of Mechaical Egieerig, MIT 3-60 Cambridge, MA 039 (slloyd@mit.edu) If quatum states exhibit small oliearities durig time evolutio, the quatum computers ca be used to solve NP-complete problems i polyomial time. We provide algorithms that solve NP-complete ad #P oracle problems by exploitig oliear quatum logic gates. It is argued that virtually ay determiistic oliear quatum theory will iclude such gates, ad the method is explicitly demostrated usig the Weiberg model of oliear quatum mechaics. This work supported i part by grat # N from the Office of Naval Research, ad by ARO ad DARPA uder grat # DAAH to QUIC, the Quatum Iformatio ad Computatio iitiative, ad by a DARPA grat to NMRQC, the Nuclear Magetic Resoace Quatum Computig iitiative.

2 Computers are physical devices: like all other physical systems, their behavior is determied by physical laws. This seemigly obvious statemet has importat implicatios, because it idicates that as our uderstadig of physical pheomea expads, the theoretical limits to the power of computig machies may grow accordigly. Recetly, it has bee show that quatum computers ca i theory exploit quatum pheomea to perform tasks that classical computers apparetly caot, such as factorig large umbers i polyomial time[], searchig databases of size M i time M [], or simulatig the detailed behavior of other quatum systems i less tha expoetial time ad space [3][4][5]. The realizatio that quatum mechaics could be used to build a fudametally more powerful type of computig machie has led to a huge amout of recet activity i the field of quatum computatio; for a review, see Ekert[7] or DiVicezo[6]. It has bee suggested [8][9][0][] [] that uder some circumstaces the superpositio priciple of quatum mechaics might be violated - that is, that the time evolutio of quatum systems might be (slightly) oliear. Such oliearity is purely hypothetical: all kow experimets cofirm the liearity of quatum mechaics to a high degree of accuracy[3][4][5][6]. Further, oliear quatum theories have ofte bee cotroversial ad frequetly have had theoretical difficulties[7] [8][9]. Nevertheless, the implicatios of oliear quatum mechaics o the theory of computatio are profoud. I particular, we show that it is geerally possible to exploit oliear time evolutio so that the classes of problems NP ad #P (icludig oracle problems) may be solved i polyomial time. A experimetal questio - that is, the exact liearity of quatum mechaics - could thereby determie the aswer to what may have previously appeared to be a purely mathematical oe. This paper therefore establishes a ew lik betwee physical law ad the theoretical power of computig machies. The class NP is the set of problems for which, oce give a potetial aswer, oe ca determie i polyomial time if the potetial aswer is i fact a solutio. These iclude all problems i the class P (those that ca be solved i polyomial time) as well as the NP-complete problems, e.g., the travelig salesma, satisfiability, ad sub-graph isomorphism, for which o kow polyomial time algorithms exist. We phrase our algorithm i terms of a oracle (or black box ), which calculates a fuctio that maps bits ito a sigle bit (i.e., it takes a iput betwee 0 ad ad returs either 0 or ). We eed to determie if there exists a iput value x for which f(x) = ; with a polyomial time algorithm to solve this problem, it is easy to solve all problems i the class NP. Ideed, this oracle problem is i fact a harder problem tha those i NP, because it clearly requires expoetial time o a classical Turig machie. A simple algorithm that solves the NP oracle problem ca be thought of as a extesio of Grover s data-base search algorithm[] to a oliear regime. Suppose that it is possible to perform a oliear operatio o a sigle qubit that has a positive Lyapuov expoet over some (ot expoetially small) regio: i.e., somewhere o the uit sphere there exists a lie of fiite extet alog which applicatio of the operatio causes earby poits to move apart expoetially

3 at a rate e λ θ proportioal to their origial separatio θ. We ca exploit this behavior to solve NP problems if we begi with a ordiary quatum computer (i.e., oe that ca perform the usual quatum logic operatios such as cotrolled rotatio gates) ad use the algorithm described below: Step. Begi by performig a π/ rotatio o each of the first qubits to obtai the state ψ = i, 0 () i=0 Step. Use the oracle (oly oce) to calculate f(i): ψ = i, f(i) () i=0 Step 3. Reverse the π/ rotatio o each of the first qubits. Each state i>the maps ito a superpositio over all possible i>, with amplitude ±. I particular, each state i>cotributes + of its amplitude to the state >, for a total cotributio of amplitude from each i>. At least of these states correspod to a particular value of f(i)=a, ad thus the state ,a>has amplitude at least /. A measuremet o the first qubits will therefore yield the state >with probability at least /4. The system will the be i the state ψ = { + s + s s 0 + s } (3) where s is the umber of solutios i for which f(i)=. The last qubit ow cotais all the iformatio that we eed; however, for small s, a measuremet of the last qubit will almost always retur 0>, yieldig o iformatio. We wish to distiguish betwee the cases s=0 ad s>0. Step 4. Repeatedly apply the oliear operatio to drive the states represetig these two cases apart at a expoetial rate: evetually, at a time determied by a polyomial fuctio of the umber of qubits, the umber of solutios s, ad the rate of spreadig (Lyapuov expoet) λ, the two cases will become macroscopically distiguishable. Step 5. Make a measuremet o the last qubit to determie the solutio. If the agular extet of the oliear regio is small, it may be ecessary to repeat the algorithm several times i order to determie the solutio with high probability. I geeral, the algorithm will require O((π/η) ) trials, where η is the agular extet of the oliear regio. If η is sufficietly large, the oracle may eed to be called oly oce. To solve problems i the class #P, we eed to determie the exact umber of solutios s. This is approximately foud by coutig the umber of times that 3

4 the oliear operator was applied. To determie s exactly, oe proceeds with fier ad fier estimates by rotatig the fial qubit such that the curret best estimate is cetered i the oliear regio; i this way, applyig the oliear operator separates states with s ear this value so that they are distiguishable. With oly a polyomial umber of iteratios, oe determies the value s exactly. The previous algorithm ca therefore by used to solve arbitrary problems (icludig oracle probems) i the class #P as well. The above algorithm has oe disadvatage i that it requires expoetial precisio. We describe below aother algorithm, which (at least i the case of the Weiberg model) appears to be robust agaist small errors. I order to simplify the descriptio, we assume (for ow) that there is at most a sigle value x for which f(x) =. We begi as before with a quatum computer that ca perform the usual quatum logic operatios, ad that ca i additio perform a simple oliear operator whose form will be described below. We coceptually divide the qubits ito iput qubits (represetig the iput, a iteger betwee 0 ad ) ad oe flag qubit. I order to solve the stated oracle problem, we use the followig algorithm: Step. We begi by performig a π/ rotatio o each of the first qubits to obtai the state ψ = i, 0 (4) Step. Use the first N qubits as the iput to the oracle ad store the output i the (N+) st qubit (the flag qubit). (Note that this is the oly time that the algorithm eeds to call the oracle; hece, we require oly evaluatio of the fuctio.) The system is ow i the state i=0 ψ = i, f(i) (5) i=0 Step 3. Cosider the first qubit separately, ad group all the states of the superpositio ito pairs based o the value of qubits... That is, the qubits.. defie subspaces of dimesio 4 = ( dimesios for qubit ) * ( dimesios for the flag qubit). Withi each subspace, the computer will be i oe of the followig states (where we write the value of the first qubit followed by the value of the flag qubit, ad igore the ormalizatio costats): (a) 00 + (b) (6) (c) (At the start of the computatio, most of the superpositio will be i the third state, because the flag qubit is >i at most oly of the compoets.) 4

5 A distictly oliear trasformatio N is the applied to these two qubits (we show below how virtually ay determiistic oliear operator ca be recast ito this form): (a) (b) (7) (c) This trasformatio is like a AND gate - it igores the idex qubit ad places the flag qubit i the state >if ad oly if either of the origial compoets had the state >for the flag qubit. Step 4. The previous step is the repeated usig each of the first qubits as the idex (ad the remaiig - qubits to defie the subspaces). After each iteratio, the umber of compoets i the superpositio that have a >for the flag qubit doubles. After iteratios, the flag qubit is o loger etagled with the first qubits: it is either i the state >for every compoet of the superpositio or the state 0>for every compoet of the superpositio. Step 5. Measure the flag qubit to determie the solutio. Thus, if oe ca perform the two qubit oliear trasformatio N oe ca fid the aswer to a NP-complete problem with certaity i polyomial (i fact liear) time, ad usig oly a sigle evaluatio of the oracle. It may be objected that the oliear operator N appears arbitrary ad uatural: ideed, it was selected exactly so as to be able to solve the stated problem. However, the apparetly arbitrary operatio N ca be built usig ordiary uitary operatios ad much simpler ad more atural sigle qubit oliear operators (that is, to the extet that ay oliear operatio i quatum mechaics ca be cosidered atural ). Oe possible techique for geeratig the trasformatio would be to use the followig steps: first, act o the two qubits with the uitary operator This trasforms the states above as follows (8) There is oe subtlety regardig oliear quatum mechaics which we should address here. Whe the superpositio priciple is abadoed, it is ot immediately clear how etagled qubits will evolve. We follow the Weiberg model, i which the time evolutio for a joit system composed of two subsystems is specified i terms of a preferred basis of vectors for the tesor product Hilbert space. For the purpose of usig the oliear dyamics to perform quatum logic, we specify the joit dyamics i terms of the basis { b>}={ >, 0...0>,...,..>} for each subsystem. The Weiberg prescriptio is as follows: write the joit state for the system Ψ as α b b ψ b ad the act o each ψ b idepedetly with b the oliear trasformatio N. 5

6 (a) (b) (c) [ 00 + ] 00 [ ] 0 (9) [ ] [ ] Next, operate o the secod qubit with a simple oe qubit oliear gate that maps both 0>ad >to the state 0>. Thus (a) (b) (c) [ 00 + ] 00 [ ] 00 (0) [ ] A The third fial state is ukow because we have ot bothered to specify how the o-liear gate acts o the state 00>+ 0>- 0>+ >. This omissio thereby allows for flexibility i choosig the gate. Whatever the state A>may be, we ca perform a uitary operatio that will trasform the first qubit ito the pure state 0>while leavig the state 00>i place. The computer is the i oe of the followig states (a) 0 0 (b) 0 0 () (c) 0 (x 0 + y ) A secod o-liear gate + is ow required that will map the state x 0>+ y >to the state >(for the particular values of x ad y which result from the above steps but ot ecessarily for arbitrary x ad y), while leavig the state 0>uchaged. After this gate is applied, the trasformatio resultig from the steps described so far is the: (a) (b) (c) [ 00 + ] 00 [ ] 00 () [ ] 0 The two qubit trasformatio N is the easily obtaied with a NOT gate o the secod qubit ad a π/ rotatio o the first qubit. 6

7 Havig thus show how to geerate N, the questio is ow reduced to that of geeratig the simpler sigle qubit gates ad +. If oe cosiders the state of a qubit as a poit o the uit sphere, the all uitary operatios correspod to rotatios of the sphere; ad while such rotatios ca place two state vectors i ay particular positio o the sphere, they ca ever chage the agle betwee two state vectors. A oliear trasformatio correspods to a stretchig of the sphere, which will i geeral modify this agle. The desired gates ad + are two particular examples of such operatios. Exceptig perhaps certai pathological cases (e.g., discotiuous trasformatios), it is evidet that virtually ay oliear operator, whe used repeatedly i combiatio with ordiary uitary trasformatios (which ca be used to place the two state vectors i a arbitrary positio o the sphere), ca be used to arbitrarily icrease or decrease the agle betwee two states, as eeded to geerate the gates ad +. We describe i detail how these gates ca be obtaied usig the model of oliear quatum mechaics put forth by Weiberg. I Weiberg s model, the Hamiltoia is a real homogeeous o-biliear fuctio h(ψ, ψ ) of degree oe, that is[9] h ψ k = ψ h k ψ k ψ k = h (3) ad state vectors time-evolve accordig to the equatio ψ k t = i h ψ k (4) Followig Weiberg [9], oe ca always perform a caoical homogeeous trasformatio such that a two-state system (i.e., a qubit) ca be described by a Hamiltoia fuctio where h = h(a) (5) = ψ + ψ (6) a = ψ (7) It is easy to verify his solutio to the time depedet oliear Schrodiger equatio (4), which is ψ k (t) = c k e iω k(a)t (8) where 7

8 ω (a) = h(a) ah (a) (9) ω (a) = h(a) + ( a)h (a) (0) For oliear h(a), oe sees that the frequecies deped o the magitude of the iitial amplitude i each basis state. Ituitively, oe ca imagie a trasformatio o the uit sphere which, istead of rotatig the sphere at a particular rate, twists the sphere i such a way so that each poit rotates at a rate which depeds upo its agle θ from the axis (clearly, this trasformatio ivolves stretchig of the surface). Oe ca exploit this stretchig of the sphere to build the gate as follows: Step. Perform a rotatio o the first qubit by a agle φ < 45 : 0 cos(φ) 0 si(φ) () si(φ) 0 + cos(φ) () Step. Time-evolve the system accordig to the oliear Hamiltoia h = h(a). Thus 0 cos(φ) 0 si(φ) α cos(φ) 0 β si(φ) (3) si(φ) 0 + cos(φ) γ cos(φ) 0 + δ si(φ) (4) where α, β, γ ad δ are phase factors. Because the iitial amplitudes of the basis states are differet i the two cases, the oliear Hamiltoia will cause the compoets to evolve at differet frequecies. As log as these frequecies are icommesurate, there is a time t at which α=γ=δ= ad β=- (to withi a accuracy ε). (Further, this time t is a polyomial fuctio of the desired accuracy ε.) The et result of these two steps is the Step 3. Reverse the first step. Thus 0 cos(φ) 0 + si(φ) (5) si(φ) 0 + cos(φ) (6) 0 cos(φ) 0 + si(φ) (7) (8) Essetially, we have reduced the agle betwee the two states by a amout φ. By suitable repetitio of this procedure (that is, by choosig φ appropriately for each iteratio), or simply by choosig φ precisely i the first step, the states 0>ad >ca be mapped to withi ε of the state 0>, i a amout of time which is a polyomial fuctio of the desired accuracy. This is the desired behavior for the oliear gate. The procedure ca be modified slightly 8

9 to icrease the agle betwee state vectors ad produce the desired behavior for the gate +. We have thus show explicitly how to solve NP-complete problems usig the Weiberg model, usig a algorithm which did ot require expoetially precise operatios. To solve the problems i the class #P, oe replaces the flag qubit with a strig of log qubits ad modifies the algorithm slightly - so that it adds the umber of solutios i each iteratio rather tha performig what is effectively a oe bit AND. I this case, a measuremet of the fial result reveals the exact umber of solutios. We have demostrated that oliear time evolutio ca i fact be exploited to allow a quatum computer to solve NP-complete ad #P problems i polyomial time. We have show explicitly how to accomplish this expoetial speed-up usig the Weiberg model of oliear quatum mechaics. I cocludig, we would like to ote that we believe that quatum mechaics is i all likelihood exactly liear, ad that the above coclusios might be viewed most profitably as further evidece that this is ideed the case. Nevertheless, the theoretical implicatios ad practical applicatios that would result from a discovery to the cotrary may warrat further ivestigatio ito the matter. D.S.A. ackowledges support from a NDSEG fellowship, ad helpful discussios with J. Jacobso, I. Siger, S. Johso, I. Park, ad T. Wag. Portios of this research were supported by grat # N from the Office of Naval Research, ad by ARO ad DARPA uder grat # DAAH to QUIC, the Quatum Iformatio ad Computatio iitiative, ad by a DARPA grat to NMRQC, the Nuclear Magetic Resoace Quatum Computig iitiative. 9

10 Refereces [] P. Shor, i Proceedigs of the 35th Aual Symposium o Foudatios of Computer Sciece, edited by S. Goldwasser (IEEE Computer Society, Los Alamos, CA, 994), p.4 [] L.K. Grover, Proceedigs of the Twety-Eighth Aual ACM Symposium o the Theory of Computig, p. x+66, -9 [3] R.P. Feyma, It. J. Theor. Phys., 467 (98) [4] S. Lloyd, Sciece 73, 073 (996) [5] D.S. Abrams ad S. Lloyd, Phys. Rev. Lett. 79, 586 (997) [6] D.P. DiVicezo, Sciece (995) [7] A. Ekert ad R. Jozsa, Reviews of Moder Physics (996) [8] S. Weiberg, Phys. Rev. Lett. 6, 485 (989) [9] S. Weiberg, A. of Phys. 94, pg. 336 (989) [0] Phys. Rev. A 56, p (997) [] B.G. Levy, Physics Today, pp. 0 (989) [] O. Bertolami, Physics Letters A, 54, p. 5-9 (99) [3] P.K. Majumder et. al., Phys. Rev. Lett. 65, 93 (990) [4] R.L. Walsworth et. al., Phys. Rev. Lett. 64, 599 (990) [5] T.E. Chupp ad R.J. Hoare, Phys. Rev. Lett. 64, 6 (990) [6] J.J. Bolliger, D.J. Heize, W.M. Itao, S.L. Gilbert, D.J. Wielad, Phys. Rev. Lett. 63, 03 (989) [7] A. Peres, Phys. Rev. Lett. 63, 4 (989) [8] J. Polchiksi, Phys. Rev. Lett. 66, pg. 397 (99) [9] N. Gisi, Phys. Lett. A 3, p. (990) 0