arxiv: v2 [quant-ph] 12 Mar 2018

Transcription

1 A new time-series model based on quantum walk Norio Konno Department of Applied Mathematics, Faculty of Engineering, Yokohama National University Hodogaya, Yokohama , Japan arxiv: v [quant-ph] Mar 08 Abstract. The quantum walk (QW) was introduced as a quantum counterpart of the classical random walk. A number of non-classical properties of the QW have been shown, e.g., ballistic spreading, anti-bellshaped limit density, localization. Since around 000, extensive research has been conducted in both theoretical aspects as well as the practical application of QWs. However, the application of a QW to the time-series analysis is not known. On the other hand, it is well known that the ARMA or GARCH models have been widely used in economics and finance. These models are studied under some suitable stationarity conditions. In this paper, we propose a new time-series model based on the QW, which does not assume such a stationarity. Therefore, our method would be applicable to the non-stationary time series. Introduction The quantum walk (QW) was introduced as a quantum counterpart of the classical random walk and it has extensively studied since around 000 ([7,, 9, ]). A number of non-classical properties of the QW have been shown, for example, ballistic spreading, anti-bellshaped limit density, and localization. Currently, active research is being conducted in both theoretical aspects as well as the practical application of QWs. The various proposals regarding QW application methods include the strongly correlated electron system [], topological insulators [9, ], and radioactive waste reduction [8, 7, 8]. Several books and reviews related to QWs have been published, for instance, Venegas-Andraca [6, 7], Konno [], Cantero et al. [5], Manouchehri and Wang [6], Portugal [3]. However, the application of a QW to the time-series analysis is not known. On the other hand, it is well known that the ARMA or GARCH models have been widely used in economics and finance, see Box and Jenkins [4], Engle [6], Bollerslev [3]. One usually considers the ARMA or GARCH models under some suitable stationarity conditions. In this paper, we present a new time-series model based on the QW, which does not assume such a stationarity. Therefore, our method would be applicable to the non-stationary time-series. As such a candidate, time series like a cryptocurrency such as Bitcoin can be considered. Concerning Bitcoin, see Nakamoto s very famous paper as founder [0]. Of course it would also be applicable to time series other than economy and finance. We should remark that our method is not related to an Itô s formula for the discrete-time QW proposed by the author []. This paper is the first step for a time-series analysis via the QW. The rest of the manuscript is organized as follows. Section is devoted to the definition of the QW considered here. In Section 3, we propose a new time-series model using the QW. Section 4 deals with a two-state one-dimensional QW case. Finally, Section 5 concludes the paper. Abbr. title: A new time-series model based on quantum walk AMS 000 subject classifications: 60F05, 60G50, 8B4, 8Q99 PACS: Lx, Fb, 0.50.Cw Keywords: Quantum walk, time-series model, stationarity

2 Definition of QW Let us first give a definition of the m-state QW on a d-dimensional integer lattice, Z d, where Z is the set of integers. The Hilbert space of the quantum walker is given by the tensor product where H p is the position space defined by H = H p H c, H p = Span{ x ; x Z d }. Moreover, H c is the coin space with an m-state given by H c = Span{ C, C,..., C m }. Here { x ; x Z d } and { C, C,..., C m } are complete orthonormal bases of H p and H c, respectively. For example, in the case of m = and d = (two-state one-dimensional QW), the quantum walker on a line, Z, has two possibilities (or chiralities), L = C (left) and R = C (right), in each step - it can move to the left or right. To each of these options, we assign a vector of the standard basis of the coin space H c = C = Span{ L, R }, such as L = [ [ 0, R =. 0] ] Here C is the set of complex numbers. In our setting, L corresponds to the fact that the value in time series decreases by. Similarly, R corresponds to the fact that the value in time series increases by. A single step of the QW is determined by the unitary operator U (s) (superscript s means the initial of system) given by U (s) = S (I p U), where S is the shift operator, I p denotes the identity on the position space H p and U is a unitary operator called quantum coin on the coin space H c. If a quantum walker starts from the origin with a coin state C, then the initial state of the QW is given by 0 C H p H c. For instance, in the case of m = and d = (two-state one-dimensional QW), S has the following form S = x Z( x x L L + x+ x R R ). If we choose the Hadamard matrix as the quantum coin, U = [ ], then the QW is the well-known Hadamard walk which has been intensively investigated in the study of the QW. Another equivalent definition of the QW based on the path counting will be presented in Section 4. Define parameter sets Θ R M and Θ R M, where M,M Z >, where R is the set of real numbers and Z > = {,,...}. For each parameter θ k = (θ () k,θ() k,...,θ(m k) k ) Θ k (k =,), we consider the corresponding quantum coin U = U(θ ) = U(θ (),θ(),...,θ(m) ) and the corresponding initial state ϕ = ϕ(θ ) = ϕ(θ (),θ(),...,θ(m) ) at the origin. In Section 4, we will deal with a time-series analysis based on a set of two-state one-dimensional QWs with the following setting. Put Θ = [0,], Θ = [0,π/](M = M = ) and θ = θ (),θ = θ (). Moreover, we define a one-parameter set of quantum coins: U = U = θ () θ () (θ () (θ () ) θ () ), ()

3 where θ () Θ = [0,]. We should note that the QW determined by the coin parameter θ () = / becomes the Hadamard walk. Moreover, we define a one-parameter set of initial states: ( ϕ = ϕ θ () [ ) = cosθ () i sinθ () ]. () Here θ () Θ = [0,π/]. From Theorem 4 in Konno [0], we see that a necessary and sufficient condition that the probability distribution of a QW with parameter (θ (),θ() ) for any time n Z is symmetric with respect to the origin is θ () = π/4, where Z = {0,,,...}. If we treat a time-series {x 0,x,x,...} in which some values do not change, e.g., x n = x n+, we take a time-series analysis using a set of three-state one-dimensional QWs. Then we consider the following setting whose quantum coin was introduced by Stefanak et al. [5] as such a candidate. The position space is and the coin space is H p = Span{ x ; x Z} H c = C 3 = Span{ L, S, R }, such as L = 0, S = 0, R = The quantum walker on a line, Z, can move to the left L or right R or stay at its current position S in each step. Then S corresponds to the fact that the value in time series does not change. Put Θ = [0,], Θ = [0,π] (M =, M = ) and θ = θ (),θ = (θ (),θ() ). We introduce a one-parameter set of quantum coins: (θ () ) θ () (θ () ) (θ () ) U = U θ () = θ () (θ () ) (θ () ) θ () (θ () ), (3) (θ () ) θ () (θ () ) (θ () ) where θ () Θ = [0,]. Remark that the QW determined by the coin parameter θ () = / 3 becomes the three-state Grover walk on Z. Moreover, we define a two-parameter set of initial states: ϕ = ϕ θ (),θ() = 3 exp(iθ () ), (4) exp(iθ () where (θ (),θ() ) Θ = [0,π]. We will move to another example. For example, in order to estimate the price of a cryptocurrency, we would need to consider not only price but also volume of the cryptocurrency in which both values have a correlation. So we treat such a data by using a set of four-state two-dimensional QWs with the following setting. The position space is and the coin space is H p = Span{ x ; x Z } ) H c = C 4 = Span{ L, R, D, U }, 3

4 such as L = 0 0, R = 0, D = 0, U = The quantum walker on a plane, Z, can move to the left L or right R or down D or up U in each step. Put Θ = [0,], Θ = [0,π] 3 (M =, M = 3) and θ = θ (),θ = (θ (),θ(),θ(3) ). Moreover, we define a one-parameter set of quantum coins introduced and intensively studied by Watabe et al. [8]: θ () θ () θ () θ () θ () θ () U = U θ () θ () θ () θ () θ () θ () θ () =, (5) θ () θ () θ () θ () θ () θ () θ () θ () θ () θ () θ () θ () where θ () Θ = [0,]. Note that the QW determined by the quantum coin parameter θ () = / becomes the four-state Grover walk on Z. Furthermore, we define a three-parameter set of initial states: ϕ = ϕ θ (),θ(),θ(3) = exp(iθ () ) exp(iθ () ). (6) exp(iθ (3) ) Here (θ (),θ(),θ(3) ) Θ = [0,π] 3. In this section, we present three examples defined by Eqs. () and (), Eqs. (3) and (4), and Eqs. (5) and (6), respectively. Of course, there are other options as well. 3 Time-Series Model via QW In this section, we propose a new time-series model based on the QW. We define the time series as a vector D n = {x 0,x, x,...,x n }, where each element x t = (x () t,x () t,...,x (d) t ) R d (t = 0,,,...,n) is an R- valued d-dimensional vector. Each one of the d values corresponds to the input variable measured in the time series. For example, if d = case, then each element x n = (x () n,x () n ) R denotes that x () n is the price and x () n of the volume of a cryptocurrency at time n, respectively. If D n is given, we want to estimate a value x n+ at the next time n + in a framework of our QW whose quantum coin is determined by U = U(θ ) = U(θ (),θ(),...,θ(m) ) and initial state is given by ϕ = ϕ(θ ) = ϕ(θ (),θ(),...,θ(m) ) at the origin. To estimate x n+, we introduce the following timedependent evaluation function V n = V n (θ,θ ) as V n = V n (θ,θ ) = n t t=0 x= t x x t µ t (x), where x = x () + x () + + x (d) for x = (x (),x (),...,x (d) ) R d and µ t (x) is the probability measure for the QW at position x and at time t. Here we assume that x 0 is the zero vector 0 without loss of generality. This assumption corresponds to the fact that our QW starts at the origin. The following is our new algorithm to estimate x n+ from a given data D n = {x 0 = 0,x,...,x n }. 4

5 Step. Find a (θ,n,θ,n) such that (θ,n,θ,n) attains the minimum of V n (θ,θ ). Step. For the (θ,n,θ,n ) given in Step, we compute E(X n+), where E(X n+ ) is the expectation of the position of quantum walker X t at time t = n+. If (θ,n,θ,n) is uniquely determined, then E(X n+ ) is considered as our estimated vector x n+. So we put x n+ = E(X n+). If we have some options of (θ,n,θ,n ), then we put x n+ as the average of E(X n+) over these choices. If V n = V n (θ,θ ) is a constant for any θ,θ, then we put x n+ = x n. Step 3. By repeating above procedures, Steps and, we obtain a sequence of estimated vectors {x,x,...,x n,...}. Therefore, for each given data D n = {x 0 = 0,x, x,...,x n }, we can get an estimated vector x n+ by using the evaluation function V n based on the QW, sequentially. The details are explained for a two-state one-dimensional case in the next section. 4 Two-State One-Dimensional Case In the first half of this section, we briefly give a definition of the two-state QW on Z via a path counting method, which is equivalent to the definition mentioned in Section. The QW is a quantum version of the classical random walk with an additional degree of freedom called chirality. The chirality takes values left and right, and it means the direction of the motion of the walker. At each time step, if the walker has the left chirality, it moves one step to the left, and if it has the right chirality, it moves one step to the right. In this paper, we put [ L = 0 ] [ 0, R = where L and R refer to the left and right chirality state, respectively. For the general setting, the time evolution of the walk is determined by a unitary matrix called quantum coin, U, where U = [ a b c d with a,b,c,d C. The matrix U rotates the chirality before the displacement, which defines the dynamics of the walk. To describe the evolution of our model, we divide the quantum coin U into two matrices: [ ] [ ] a b 0 0 P =, Q = 0 0 c d with U = P + Q. The important point is that P (resp. Q) represents that the walker moves to the left (resp. right) at position x at each time step. In the present paper, we take ϕ = T [α,β] with α, β C and α + β = as the initial qubit state, where T is the transpose operator. Let Ξ n (l,m) denote the sum of all paths starting from the origin in the trajectory consisting of l steps left and m steps right at time n with n = l+m. For example, Ξ (,) = PQ+QP, Ξ 4 (,) = P Q +Q P +PQPQ+QPQP +PQ P +QP Q. The probability that our quantum walker is in position x( Z) at time n ( Z ) starting from the origin with ϕ = T [α,β] with α, β C and α + β = is defined by ], ], P(X n = x) = Ξ n (l,m)ϕ, 5

6 where n = l+m and x = l+m. Let µ n (x) = P(X n = x). So we have µ 0 (0) =, µ ( ) = Pϕ, µ () = Qϕ, µ ( ) = P ϕ, µ (0) = (PQ+QP)ϕ, µ () = Q ϕ. We define the probability amplitude of the QW in position x at time n by [ ] Ψ L Ψ n (x) = n (x) Ψ R n (x). Then we see that P(X n = x) = Ψ n (x) = Ψ L n(x) + Ψ R n(x). In this section, we focus on the following setting which is essentially equivalent to one given in Eqs. () and (): [ ] [ ] cosθ sinθ cosξ U = U(θ) =, ϕ = ϕ(ξ) = (0 θ,ξ π/). (7) sinθ cosθ i sinξ Note that θ = θ () = θ and θ = θ () = ξ. Let µ n = µ n (θ,ξ) be a probability measure for this QW. So we have a sequence of probability measures {µ 0,µ,...,µ n,...}. Assume that D n = {x 0,x,...,x n } is a set of R-valued time-series data until time n. If D n is given, we want to estimate a next time x n+ by using a framework of our QW with a pair of parameters (θ,ξ). To do so, we introduced the following time-dependent evaluation function V n = V n (θ,ξ) as V n = V n (θ,ξ) = n t t=0 x= t (x x t ) µ t (x). We should note that we can not define the joint distribution of the QW such as P(X t = x, X t+ = y). Here we assume that x 0 = 0 which corresponds to the fact that our QW starts at the origin. The following was our new algorithm to estimate x n+ from a given data {x 0 = 0,x,...,x n }. Step. Find a (θ n,ξ n ) such that (θ n,ξ n ) attains the minimum of V n(θ,ξ). Step. For the (θn,ξ n ) given in Step, we compute E(X n+). If (θ,n,θ,n ) is uniquely determined, then E(X n+ ) is considered as our estimated vector x n+. So we put x n+ = E(X n+). If we have some options of (θ,n,θ,n), then we put x n+ as the average of E(X n+ ) over these choices. If V n = V n (θ,θ ) is a constant for any θ,θ, then we put x n+ = x n. Step 3. By repeating above procedures, Steps and, we obtain a sequence of estimated values {x,x,...,x n,...}. From now on, we consider Steps and for n = 0 and n = in the setting given by Eq. (7). First we treat n = 0 case. We want to obtain x from {x 0 = 0}. This case is trivial, i.e., x = x 0 = 0. In fact, we see that V 0 = V 0 (θ,ξ) = (0 x 0 ) µ 0 (0) = 0. Thus V 0 = V 0 (θ,ξ) = 0 for any θ,ξ. Therefore we can not determine (θ 0,ξ 0). Then we put x = x 0 = 0. Next we deal with n = case. We want to obtain x from {x 0 = 0,x }. 6

7 Step. We begin with V = V (θ,ξ) = = V 0 + x= since V 0 = 0 and µ (0) = 0. Thus we have t t=0 x= t (x x ) µ (x) (x x t ) µ t (x) = ( x ) µ ( )+(0 x ) µ (0)+( x ) µ () = ( x +x + ) µ ( )+ ( x x + ) µ (), V = ( x +x + ) µ ( )+ ( x x + ) µ (). (8) From now on we compute µ ( ) and µ () as follows. [ ][ ] [ ] cosθ sinθ cosξ cosθcosξ +i sinθsinξ Pϕ = =, 0 0 i sinξ 0 [ ][ ] [ ] 0 0 cosξ 0 Qϕ = =. sinθ cosθ i sinξ sinθcosξ i cosθsinξ By using these, we get Note that From Eqs. (8), (9), (0), and (), we have So we get Moreover, µ ( ) = Pϕ = cos θ cos ξ +sin θ sin ξ, (9) µ () = Qϕ = sin θ cos ξ +cos θ sin ξ. (0) µ ( )+µ () =. () V = ( x +x + ) µ ( )+ ( x x + ) µ () = {µ ( )+µ ()} x +{µ ( ) µ ()} x +{µ ( )+µ ()} = x +cos(θ) cos(ξ)x +. V = V (θ,ξ) = x +cos(θ) cos(ξ)x +. () V θ = 4sin(θ) cos(ξ)x, (3) V ξ = 4cos(θ) sin(ξ)x. (4) Here we consider three cases (i) x > 0, (ii) x = 0, (iii) x < 0 as follows. Case (i) x > 0. By Eqs. (), (3), and (4), we consider four cases in the following way. If (θ,ξ ) = (0,0), thenv = V (0,0) = x +x +, ( If (θ,ξ ) = 0, π ) (, thenv = V 0, π ) = x x +, ( π ) ( π ) If (θ,ξ ) =,0, thenv = V,0 = x x +, ( π If (θ,ξ ) =, π ( π, thenv = V ), π ) = x +x +. 7

8 Note that x +x + > x x +, since x > 0. Therefore we see that ( (θ,ξ ) = 0, π ) ( π ),,0 and V = V (θ,ξ ) = x x +. Case (ii) x = 0. In this case, we see that V = V (θ,ξ) = for any θ,ξ. Case (iii) x < 0. As in the case of (i), we see that ( π (θ,ξ ) = (0,0),, π ) and V = V (θ,ξ ) = x +x +. Step. We compute E(X ) for the QW with (θ,ξ) = (θ,ξ ) as follows. We begin with Thus we have E(X ) = x= xµ (x) = ( )µ ( )+µ () = ( ) P ϕ + Q ϕ = cos θ cos(θ) cos(ξ). E(X ) = cos θ cos(θ) cos(ξ). (5) By using Eq. (5), we consider our estimated value x in the following way. Case (i) x > 0. Furthermore, we consider two cases: (a) if (θ,ξ ) = (0,π/), then Eq. (5) gives (b) if (θ,ξ ) = (π/,0), then Eq. (5) gives E(X ) =. (6) E(X ) = 0. (7) In this case, we can not determine x. So we assume that two outcomes, (θ,ξ ) = (0,π/) and (θ,ξ ) = (π/,0), are equally likely. That is, each event is selected with probability /. Thus Eqs. (6) and (7) imply Finally we obtain x =. x = +0 =. Case (ii) x = 0. We can not determine x. Thus we put x = x (= 0). Case (iii) x < 0. As in the case of (i), we consider two cases: 8

9 (a) if (θ,ξ ) = (0,0), then Eq. (5) gives (b) if (θ,ξ ) = (π/,π/), then Eq. (5) gives E(X ) =. (8) E(X ) = 0. (9) In this case also, we can not determine x. Thus we assume that two outcomes, (θ,ξ ) = (0,0) and (θ,ξ ) = (π/,π/), are equally likely. From Eqs. (8) and (9), we see x = ( ) +0 =. Therefore we get x =. To summarize n = case, we have the following Table. x 0 x x x 0 = 0 x > 0 x = x 0 = 0 x = 0 x = x = 0 x 0 = 0 x < 0 x = Table Step 3. By Steps and, for a given data {x 0 = 0,x }, we obtained a sequence of estimated values {x = x 0 = 0,x }. To summarize this, we present the following Table. x x 0 (x > 0) 0 0 (x = 0) 0 (x < 0) Table If we consider general n( ) case, then the following expression of E(X n ) given by Proposition in Konno [0] is useful to compute E(X n ) for n 3 and θ [0,π/): E(X n ) = (cosθ) (n ) [ncos(θ) + [ n ] k= k k γ= δ= ( sin θ cos θ (n k) {ncos(θ)+γ +δ} γδ ) γ+δ k k n k n k γ δ γ δ ] cos(ξ), where [x] is the greatest integer that is less than or equal to x R. For example, we have E(X 3 ) = {( 3cos 4 θ +sin 4 θ ) cos(θ)+sin θ sin (θ) } cos(ξ). We should remark that θ = π/ is a trivial case. Furthermore, in orderto determine x n forlargen, it wouldbe better to useanumericalmethod compared with the analytical method discussed as in the second half of this section [4]. 9

10 5 Conclusions and Future Work In this paper, we have proposed a new time-series method based on the QW. Additionally, we have discussed the method for a two-state one-dimensional case. On of the interesting problems would be to compare our method with the ARMA and/or GARCH models for some real data [4], since our model does not impose a stationarity, but so do the ARMA and GARCH models. Our model presented here is based on the QW. Recently the author extended the QW to a new walk called quaternionic quantum walk (QQW) determined by a unitary matrix whose component is quaternion [3]. In general, the behavior of QQW is different from usual QW [5, 4], it is interesting to compare our method with a time-series one based on the QQW model. Moreoveran extension from the QQW time-series model to the Clifford algebra time-series one would be also attractive. Acknowledgment. The author would like to thank Hiwon Yoon, Song-Ju Kim, Shinya Kawata, Masato Takei, Takashi Komatsu for useful discussions. References [] Aharonov, Y., Davidovich, L., Zagury, N.: Quantum random walks. Phys. Rev. A 48, (993) [] Ambainis, A., Bach, E., Nayak, A., Vishwanath, A., Watrous, J.: One-dimensional quantum walks. In: Proceedings of the 33rd Annual ACM Symposium on Theory of Computing, pp.37 49, 00. [3] Bollerslev, T.: Generalized autoregressive conditional heteroskedasticity. Journal of Ecomometrics 3, (986) [4] Box, G. E. P., Jenkins, G. M.: Time Series Analysis: Forecasting and Control. Holden-Day, CA 970. [5] Cantero, M. J., Grünbaum, F. A., Moral, L., Velázquez, L.: The CGMV method for quantum walks. Quantum Inf. Process., 49 9 (0) [6] Engle, R. F.: Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation. Econometrica 50, (98) [7] Gudder, S. P.: Qunatum Probability. Academic Press Inc. CA (988) [8] Ichihara, A., Matsuoka, L., Kurosaki, Y., Yokoyama, K.: An analytic formula for describing the transient rotational dynamics of diatomic molecules in an optical frequency comb. Chin. J. Phys. 5, (03) [9] Kitagawa, T.: Topological phenomena in quantum walks: elementary introduction to the physics of topological phases. Quantum Inf. Process., (0) [0] Konno, N.: Quantum random walks in one dimension. Quantum Inf. Process., (00) [] Konno, N.: Quantum Walks. In: Quantum Potential Theory, Franz, U., and Schürmann, M., Eds., Lecture Notes in Mathematics: Vol. 954, pp , Springer-Verlag, Heidelberg (008) [] Konno, N.: A note on Itô s formula for discrete-time quantum walk. J. Compu. Theo. Nanosci. 0, (03) [3] Konno, N.: Quaternionic quantum walks. Quantum Stud.: Math. Found., (05) [4] Konno, N., Kawata, S., Kim, S-J, Yoon, H.: In preparation. [5] Konno, N., Mitsuhashi, H., Sato, I.: The discrete-time quaternionic quantum walk on a graph. Quantum Inf. Process. 5, (06) [6] Manouchehri, K., Wang, J.: Physical Implementation of Quantum Walks. Springer, Berlin (03) [7] Matsuoka, L., Kasajima, T., Hashimoto, M., Yokoyama, K.: Numerical study on quantum walks implemented on cascade rotational transitions in a diatomic molecule. J. Korean Phys. Soc. 59, (0) [8] Matsuoka, L., Yokoyama, K.: Physical implementation of quantum cellular automaton in a diatomic molecule. J. Compu. Theo. Nanosci. 0, (03) [9] Meyer, D. A.: From quantum cellular automata to quantum lattice gases. J. Statist. Phys. 85, (996) [0] Nakamoto, S.: Bitcoin: A peer-to-peer electronic cash system. (008) [] Obuse, H., Kawakami, N.: Topological phases and delocalization of quantum walks in random environments. Phys. Rev. B 84, 9539 (0) 0

11 [] Oka, T., Konno, N., Arita, R., Aoki, H.: Breakdown of an electric-field driven system: a mapping to a quantum walk. Phys. Rev. Lett. 94, 0060 (005) [3] Portugal, R.: Quantum Walks and Search Algorithms. Springer, Berlin (03) [4] Saito, K.: Probability distributions of quaternionic quantum walks. arxiv: (07) [5] Stefanak, M., Bezdekova, I., Jex, I.: Continuous deformations of the Grover walk preserving localization. Eur. Phys. J. D, 4 (0) [6] Venegas-Andraca, S. E.: Quantum Walks for Computer Scientists. Morgan and Claypool (008) [7] Venegas-Andraca, S. E.: Quantum walks: a comprehensive review. Quantum Inf. Process., (0) [8] Watabe, K., Kobayashi, N., Katori, M., Konno, N.: Limit distributions of two-dimensional quantum walks. Phys. Rev. A (008)