arxiv:cond-mat/ v1 20 Sep 2004

Size: px

Start display at page:

Download "arxiv:cond-mat/ v1 20 Sep 2004"

Charla Potter
5 years ago
Views:

1 One-dimensional continuous-time quantum walks DANIEL BEN-AVRAHAM ERIK M. BOLLT CHRISTINO TAMON September, 4 arxiv:cond-mat/4954 v Sep 4 Abstract We survey the equations of continuous-time quantum walks on simple one-dimensional lattices, which include the finite and infinite lines and the finite cycle, and compare them with the classical continuous-time Markov chains. The focus of our expository article is on analyzing these processes using the Laplace transform on the stochastic recurrences. The resulting time evolution equations, classical versus quantum, are strikingly similar in form, although dissimilar in behavior. We also provide comparisons with analyses performed using spectral methods. Introduction The theory of Markov chains on countable structures is an important area in mathematics and physics. A quantum analogue of continuous-time Markov chains on the infinite line is well-studied in physics (for example, it can be found in [], Chapters 3 and 6. More recently, it was studied by Aharonov et al. [] and by Farhi and Gutmann []. The latter work placed the problem in the context of quantum algorithms for search problems on graphs. Here the symmetric stochastic matrix of the graph is viewed as a Hamiltonian of the quantum process. Using Schrödinger s equation with this Hamiltonian, we obtain a quantum walk on the underlying graph, instead of a classical random walk. Recent works on continuous-time quantum walks on finite graphs include the analyses of mixing and hitting times on the n-cube [6, 4], of mixing times on circulant graphs and Cayley graphs of the symmetric group [6, 3], and of hitting times on path-like graphs [8, 9]. Most of these are structural results based on spectral analysis of the underlying graphs, such as the n-cube, circulant and Cayley graphs, and (weighted paths. For example, Moore and Russell [6] proved that the mixing time of a quantum walk on the n-cube is asymptotically faster than a classical random walk; Kempe [4] proved that the hitting time for vertices on opposite ends of the n-cube is exponentially faster than in a classical random walk. Ahmadi et al. [6] and Gerhardt and Watrous [3] proved that circulants and the Cayley graph of the symmetric group lack the uniform mixing property found in classical random walks. A recent work of Childs et al. [9] gave intriguing evidence that continuous-time quantum walk is a powerful method for designing new quantum algorithms. They analyzed diffusion processes on onedimensional structures (finite path and infinite line using spectral methods. Another work by Childs and Dept. Physics, Clarkson University. qd@clarkson.edu. Supported by NSF grant NSF PHY-494. Dept. Mathematics and Computer Science, and Dept. Physics, Clarkson University. bolltem@clarkson.edu. Supported by NSF grant NSF DMS Dept. Mathematics and Computer Science, and Center for Quantum Device Technology, Clarkson University. tino@clarkson.edu. Supported by NSF grant DMR-46.

2 Goldstone [] explored the application of continuous-time quantum walks to perform Grover search on spatial lattices. There is an alternate theory of discrete quantum walks on graphs, which we will not discuss here. This alternate model was studied in Aharonov et al. [4] and Ambainis et al. [5], but had appeared earlier in work by Meyer [5]. The work by Ambainis et al. [5] had also focused on one-dimensional lattices. Recently, Ambainis [] developed an optimal (discrete quantum walk algorithm for the fundamental problem of Element Distinctness. This offers another idea for developing quantum algorithms. We survey and (rederive equations for the continuous-time classical and quantum walks on onedimensional lattices using the Laplace transform that works directly with the recurrences. The Laplace transform is a well-known tool in stochastic processes (see [7] and it might offer a useful alternative to the Fourier transform in certain settings.. Stochastic walks on graphs Let G (V,E be a simple (no self-loops, countable, undirected graph with adjacency matrix A. Let D be a diagonal matrix whose j-th entry is the degree of the j-th vertex of G. The Laplacian of G is defined as H A D. Suppose that P(t is a time-dependent probability distribution of a stochastic (particle process on G. The classical evolution of the continuous-time walk is given by the Kolmogorov equation P (t HP(t. ( The solution this equation, modulo some conditions, is P(t e th P(, which can be solved by diagonalizing the symmetric matrix H. This spectral approach requires full knowledge of the spectrum of H. A quantum analogue of this classical walk uses the Schrödinger equation in place of the Kolmogorov equation. Let ψ : V (G C be the time-independent amplitude of the quantum process on G. Then, the wave evolution is iħ d ψ(t Hψ(t. ( dt Assuming ħ for simplicity, the solution of this equation is ψ(t e ith ψ(, which, again, is solvable via spectral techniques. The classical behavior of this quantum process is given by the probability distribution P(t whose j-th entry is P j (t ψ j (t, where ψ j (t j ψ(t. The average probability of vertex T j is defined as P(j lim T T P j(tdt (see [4] Figure : Examples of some one-dimensional lattices. From left to right: Z, P 4, Z 4. The table in Figure shows the known equations for continuous-time stochastic walks on the infinite (integer line Z {...,,,,,,...}, the finite cycle Z N {,...,N } on N vertices, and the finite path P N {,...,N} on N + vertices, in terms of the two kinds of Bessel functions I( and J(. We assume here that the particle is initially at. The plots in Figures 3 and 4 show the dissimilar behavior of the classical versus quantum walks.

3 Graph Classical walk Quantum walk P j (t probability on vertex j at time t ψ j (t amplitude on vertex j at time t Z e t I j (t (i j J j (t Z N P N α ±j(modn α ±j(mod N e t I α (t e t I α (t α ±j(modn α ±j(mod N (i α J α (t (i α J α (t Figure : The equations of the continuous-time classical versus quantum walks on the infinite line, finite cycle, and the finite line, assuming the particle starts at position.. Laplace transform The Laplace transform of a time-dependent function P(t, denoted ˆP(s L{P(t}, is defined as L{P(t} e st P(t dt. (3 The only basic properties of the Laplace transform which we will need are (see [3]: Linearity: L{aP(t + bq(t} a ˆP(s + b ˆQ(s Derivative: L{P (t} s ˆP(s P( Shifting: L{e at P(t} ˆP(s a The relevant Inverse Laplace transform involving the Bessel functions are (for ν > : ˆP(s (s s a ν s a P(t a ν I ν (at (Eqn in [3] (4 ˆP(s ( s + a s ν s + a P(t a ν J ν (at (Eqn in [3] (5 The infinite line Classical process. The Kolmogorov equation for the infinite line is P j(t P j(t P j (t + P j+(t, (6 with initial value P j ( δ,j. The Laplace transform of (6 is ˆP j+ (s (s + ˆP j (s + ˆP j (s P j (. (7 3

4 The solution of q (s + q + is q ± (s + ± (s +. A natural guess of the solution is { Aq ˆP j j (s + if j < Aq j (8 if j > When j, we get A ( + s q. Thus, for j Z, q j ((s + (s + ˆP j (s j. (9 ( + s q (s + Using the Inverse Laplace transform (4, after shifting S s +, we get P j (t e t I j (t. ( This is a probability function, since e t/(z+/z k zk I k (t, if z (see Eqn in [3]. Quantum process. The Schrödinger equation for the infinite line is The Laplace transform of ( is iψ j (t ψ j(t + ψ j+(t. ( ˆψ j+ (s i(s ˆψ j (s ψ j ( + ˆψ j (s ( The solutions of q isq + are q ± i(s ± s +, where q + q. A guess for the solution is { Aq ˆψ j j (s + if j Aq j if j > (3 When j, we get A (s + iq. Thus, q ˆψ j s j (s (s + iq (i j ( + s j. (4 s + The Inverse Laplace transform (5 yields, for j Z, ψ j (t (i j J j (t, (5 This is a probability function, since J (z + k J k (z (see Eqn in [3]. Spectral analysis. Let H be a Hamiltonian defined as j H k if j k ±, and otherwise. For each p [π,π], define p so that j p e ipj. (6 π The eigenvalue equation H p λ p p has the solution λ p cos(p. Thus, the amplitude of position j when the particle starts at position is j e iht π π Childs et al. [9] gave a more generalized analysis along these lines. π e ipj e itcos(p dp (i j J j (t (see Eqn. 9.. in [3] (7 4

5 3 The finite cycle Classical process. If A is the adjacency matrix of the finite cycle, let H A I be its Laplacian matrix. The Kolmogorov equation for the finite cycle is Applying the Laplace transform to (8, we get P j(t P j(t P j (t + P j+(t. (8 (s + ˆP j (s P j ( ˆP j (s + ˆP j+ (s. (9 For convenience, define the extra condition ˆP (s ˆP N (s +, so that ˆP j+ (s (s + ˆP j (s + ˆP j (s holds for j Z N. The cycle condition is ˆP N (s ˆP (s. We guess the solution to be ˆP j (s Aq j + + Bqj, ( where q ± is the solution to x (s + x +, i.e., q ± (s + ± (s +, with q + q. Using the cycle condition, we get Aq N + + Bq N A + B A(q N + B( q N B Aq N +. ( Using the extra condition and (, we get A ((q + q (q N +. Thus, for j Z N, ˆP j (s Aq j + + Bqj A(qj + + qnj + (q + q (q j + qnj ( q N (q + q [ ((s + (s + kn+j (s + ( q kn+j The Inverse Laplace transform (4, after shifting, yields, for j Z N, P j (t e t [ I kn+j (t + I (k+nj (t ] + q (k+nj + ((s + (s + (k+nj (s + α ±j(modn ] e t I α (t. (. Quantum process. Since the finite cycle is a regular graph, instead of the Laplacian, we use the adjacency matrix directly. In a continuous-time quantum walk, this simply introduces an irrelevant phase factor in the final expression. The Schrödinger equation, in this case, is The Laplace transform of (3 is iψ j (t ψ j(t + ψ j+(t. (3 ˆψ j+ (s i(s ˆψ j (s ψ j ( + ˆψ j (s (4 The cycle boundary condition is ˆψ N (s ˆψ (s. For convenience, define ˆψ (s ˆψ N (s i. (5 5

6 The solutions of q isq + are q ± i(s ± s +, with q + q. A solution guess, for j Z N, is ˆψ j (s Aq j + + Bqj. (6 The cycle boundary condition yields B Aq N +. By (5, we get A i((q + q (q N +. Thus, for j Z N, ˆψ j (s Aq j + + Bqj A(qj + + qnj + i (q + q (q j + qnj ( q N [ ((i( s + s kn+j s + The Inverse Laplace transform (5 gives, for j Z N, ψ j (t i (q + q ( q kn+j [ ] (i kn+j J kn+j (t + (i (k+nj J (k+nj (t + q (k+nj + ((i( s + s (k+nj s + α ±j(modn ]. (i α J α (t. (7 Spectral analysis. The normalized adjacency matrix H of Z N is the circulant matrix /... / / /... H /.... ( /... / It is well-known that all N N circulant matrices are unitarily diagonalized by the Fourier matrix F N V (ω N, where ω N e πi/n and V (ω N is the Vandermonde matrix defined as V (ω N [j,k] ω jk N, for j,k {,,...,N }. The eigenvalues of H are λ j (ωj N + ωj(n N cos(πj/n, for j,,...,n. Thus, the wave amplitude at vertex j at time t is ψ j (t N N From earlier analysis, we get the following Bessel equation N e itcos(πk/n e πijk/n N e itcos(πk/n ω jk N. (9 α ±j(mod N (i α J α (t. (3 It is an open question if there exists a time t R + such that for all j Z N we have ψ j (t /N, i.e., uniformity is achieved at some time t. For N,3,4, it is known that instantaneous exact uniform mixing is achieved (see [6, 6]. 6

7 4 The finite path Classical process. Let A be the normalized adjacency matrix of the finite path, where A is a stochastic matrix with the probability transitions are proportional to the degrees of the vertices. Let H A I be its Laplacian. Then, the Kolmogorov equation, in this case, is P j(t P j(t P j (t + P j+(t, (3 for < j < N, with initial condition P j ( δ j, and boundary conditions The Laplace transform of (3 is P (t P (t P (t, P N(t P N (t P N (t. (3 ˆP j+ (s (s + ˆP j (s + ˆP j (s, < j < N, (33 and two boundary equations ( + s ˆP (s ˆP (s, and ( + s ˆP N (s ˆP N (s. A guess of the solution is ˆP j (s Aq j + + Bqj, j N, (34 where q ± (s+ ± (s +. The boundary equations give B A /(q + q and A Bq N. Combining these last two equations, we get Thus, for j,,...,n, A (q + q q N ( q N. (35 ˆP j (s Aq j + + Bqj A(qj + + qnj + (q + q (q Nk+j + q N(k+j. The Inverse Laplace transform (4, after shifting, yields, for j,,...,n, P j (t e t [ I Nk+j (t + I N(k+j (t ] (q j + qnj (q + q ( q N α ±j(mod N e t I α (t. (36 Quantum process. The Schrödinger equation for the finite path is iψ j(t ψ j(t + ψ j+(t, (37 for < j < N, with initial condition ψ j ( δ,j and boundary conditions iψ (t ψ (t, iψ N(t ψ N (t. (38 The Laplace transform of (37 is ˆψ j+ (s is ˆψ j (s + ˆψ j (s, < j < N, (39 7

8 and two boundary equations is ˆψ (s i ˆψ (s, and is ˆψ N (s ˆψ N (s. The solutions of q isq + are q ± i(s ± s +. A guess of the solution is ˆψ j (s Aq j + + Bqj, j N. (4 From the boundary equations, we get B A i/(q + q and B Aq N +. Thus, For j,,...,n, A i (q + q q N ( q N. (4 ˆψ j (s Aq j + + Bqj A(qj + + qnj + i (q + q (q j + qnj ( q N [ ((i( s + s Nk+j s + i (q + q The Inverse Laplace transform (5 yields, for j,,...,n, ψ j (t ( q Nk+j [ ] (i Nk+j J Nk+j (t + (i N(k+j J N(k+j (t + q (k+nj + ((i( s + s (k+nj s + α ±j(mod N ]. (i α J α (t. (4 Spectral analysis. The spectrum of a path on n vertices is given by Spitzer [7]. For j {,,..., N}, the eigenvalue λ j and its eigenvector v j are given by ( ( (j + π (j + π λ j cos, v j (l N + N + sin (l +. (43 N + The probability of measuring vertex at time t is given by 4 P (t (N + 4 (N + j j,k sin ( (j + π N + sin 4 ( (j + π N + sin ( (k + π N + e it(λ j λ k. (44 Since all eigenvalues are distinct, the average probability of measuring the starting vertex is 4 ( ( (j + π (k + π P( (N + sin sin T lim e it(λ j λ k dt N + N + T T j,k. Equating this with (4, we obtain a Bessel-like equation: T lim T T (i a J a (t dt a (mod N 4 (N + N ( (k + π sin 4 N + (45 8

9 5 Conclusions This expository survey reviews equations for the continuous-time quantum walks on one-dimensional lattices. The focus was on analysis based on the Laplace transform which works directly with the stochastic recurrences. It would be interesting to extend this analysis to higher-dimensional or to regular graph-theoretic settings. Another interesting direction is to consider lattices with defects and weighted graphs [9]. Acknowledgments We thank the anonymous referees for pointing out several omissions and for helpful comments that improve the presentation of this manuscript. References [] A. Ambainis, Quantum walk algorithm for element distinctness, quant-ph/3. [] Y. Aharonov, L. Davidovich, N. Zagury (993, Quantum Random Walks, Phys. Rev. Lett., 48(: [3] M. Abramowitz, I. Stegun (974, Handbook of Mathematical Functions With Formulas, Graphs, and Mathematical Tables, Dover. [4] D. Aharonov, A. Ambainis, J. Kempe, and U. Vazirani (, Quantum Walks on Graphs, in Proc. 33rd ACM Ann. Symp. Theory Computing (ACM Press, [5] A. Ambainis, E. Bach, A. Nayak, A. Viswanath, and J. Watrous (, One-dimensional Quantum Walks, in Proc. 33rd ACM Ann. Symp. Theory Computing (ACM Press, [6] A. Ahmadi, R. Belk, C. Tamon, C. Wendler (3, On Mixing in Continuous-time Quantum Walks on Some Circulant Graphs, Quantum Information and Computation, 3(6:6-68. [7] K. L. Chung (967, Markov Chains with Stationary Transition Probabilities, nd. ed., Springer. [8] A. Childs, E. Farhi, and S. Gutmann (, An example of the difference between quantum and classical random walks, Quantum Information Processing, 35. [9] A. Childs, E. Deotto, R. Cleve, E. Farhi, S. Gutmann, D. Spielman (3, Exponential algorithmic speedup by quantum walk, Proc. 35th ACM Ann. Symp. Theory of Computing (ACM Press, [] A. Childs, J. Goldstone (3, Spatial search by quantum walk, quant-ph/3654. [] E. Farhi and S. Gutmann (998, Quantum computation and decision trees, Phys. Rev. A, 58: [] R. Feynman, R. Leighton, M. Sands (965, The Feynman Lectures on Physics, vol. III, Addison- Wesley. [3] H. Gerhardt and J. Watrous (3, Continuous-time quantum walks on the symmetric group, in Proc. 7th Int. Workshop Rand. Approx. in Computer Science, LNCS 764, Springer,

10 [4] J. Kempe (3, Quantum Random Walks Hit Exponentially Faster, in Proc. 7th Int. Workshop Rand. Approx. in Computer Science, LNCS 764, Springer, [5] D. Meyer (996, From quantum cellular automata to quantum lattice gases, J. Stat. Phys., 85: [6] C. Moore and A. Russell (, Quantum Walks on the Hypercube, in Proc. 6th Int. Workshop Rand. Approx. in Computer Science, LNCS 483, Springer, [7] F. Spitzer (976, Principles of Random Walk, nd ed., Springer Figure 3: Stochastic walks on the infinite line Z: (a plot of P (t in the continuous-time random walks for t [,5]. (b plot of ψ (t in a continuous-time quantum walk for t [,5]. Both processes exhibit exponential decay, but with the quantum walk showing an oscillatory behavior.

11 Figure 4: Stochastic walks on the finite cycle Z 5, each approximated using terms: (a plot of P (t in the continuous-time random walks for t [,5]. (b plot of ψ (t in the continuous-time quantum walk for t [, 5]. The classical walk settles quickly to /5, while the quantum walk exhibit a short-term chaotic behavior and a long-term oscillatory behavior below..

A note on graphs resistant to quantum uniform mixing

A note on graphs resistant to quantum uniform mixing William Adamczak Kevin Andrew Peter Hernberg Christino Tamon Abstract Continuous-time quantum walks on graphs is a generalization of continuous-time