arxiv:quant-ph/ v4 2 Nov 2005

Transcription

1 Quantum Information and Computation, Vol. 0, No. 0 (003) c Rinton Press A new algorithm for producing quantum circuits using KAK decompositions arxiv:quant-ph/ v4 Nov 005 Yumi Nakajima a, Yasuhito Kawano, and Hiroshi Sekigawa NTT Communication Science Laboratories, Nippon Telegraph and Telephone Corporation 3-1 Morinosato-Wakamiya, Atsugi, Kanagawa , JAPAN Received (received date) Revised (revised date) We provide a new algorithm that translates a unitary matri into a quantum circuit according to the G = KAK theorem in Lie group theory. With our algorithm, any matri decomposition corresponding to type-aiii KAK decompositions can be derived according to the given Cartan involution. Our algorithm contains, as its special cases, Cosine- Sine decomposition (CSD) and Khaneja-Glaser decomposition (KGD) in the sense that it derives the same quantum circuits as the ones obtained by them if we select suitable Cartan involutions and square root matrices. The selections of Cartan involutions for computing CSD and KGD will be shown eplicitly. As an eample, we show eplicitly that our method can automatically reproduce the well-known efficient quantum circuit for the n-qubit quantum Fourier transform. Keywords: G = KAK theorem, KGD, CSD, decomposition of the QFT Communicated by: to be filled by the Editorial 1 Introduction Decomposing a unitary matri into an efficient sequence of elementary gates is a fundamental problem in designing quantum circuits. There are two types of decomposition: One is eact decomposition where an arbitrary unitary matri is decomposed precisely into a sequence of elementary gates, such as arbitrary single-qubit rotations and the CNOT. The other involves approimate strategies by which an arbitrary unitary matri is decomposed approimately into a sequence of a fied set of elementary gates, as shown in Solovay-Kitaev theorem (cf. [1], Appendi 3). In this paper, we treat eact decomposition. Cosine-Sine decomposition (CSD), which is a well-known algorithm in numerical linear algebra, was the first algorithm utilized for this purpose []. CSD applies the well-known algorithm for computing generalized singular value decomposition (GSVD). In CSD, we first divide an input matri g into four square matrices and then apply SVD to each matri. Then, we have g = ( ) g11 g 1 g 1 g = ( )( )( ) u1 0 c s v1 0 0 u s c 0 v = UΣV, (1) where c = diag[cos(ζ 1 ),cos(ζ ),,cos(ζ n 1)], s = diag[sin(ζ 1 ),sin(ζ ),,sin(ζ n 1)], then g 11 = u 1 cv 1, g 1 = u 1 sv, g 1 = u sv 1, and g = u cv are SVD. Recursively pera nakajima@theory.brl.ntt.co.jp 1

2 A new algorithm for producing quantum circuits using KAK decompositions forming decomposition (1) leads to a quantum circuit. The CSD-based algorithms are easy to implement on a computer because algorithms for calculating GSVD are well-known, and software libraries including GSVD are available. Some CSD-based algorithms [3 6] have been investigated with the aim of improving Barenco s result that an arbitrary n n unitary matri is composed of O(n 4 n ) elementary gates [7]. And improvement to O(4 n ) elementary gates has been reported by Möttönen et al. [4] and by Shende et al [6]. On the other hand, Khaneja and Glaser provided another kind of decomposition [8], which waslaternamedkgd.kgdlies withinthe frameworkoftheg = KAK theorem(cf. [9], Theorem8.6)in Lie grouptheory. This theoremshowsthat an element g SU( n ) is decomposed into matri products k 1 ak for some k 1,k ep(k) and a ep(h). Here, su( n ) = k m is a Cartan decomposition in Lie algebra su( n ), k and m = k are orthogonal vector spaces contained in su( n ), and h is a maimal Abelian subalgebra (a Cartan subalgebra) contained in m (cf. [15], VI.). Matrices k 1, a, and k are not uniquely determined from g. They depend on the selections of the bases of k, m, and h; besides, they are not determined even if bases are selected. Khaneja and Glaser provided a particular selection of bases of k, m, and h in Ref. [8] so that the selection matches an NMR system, and they provedthat a time-optimal control on a two-qubit NMR quantum computer can be obtained from the decomposition [10]. Thus, KGD can be regarded as the G = KAK theorem on the particular bases. It should be noted that KGD does not give a unique translation of the input matri into a quantum circuit. Bullock [11] showed that CSD can also be regarded in the framework of the G = KAK theorem; i.e., CSD uses the type-aiii KAK decomposition with the global Cartan decomposition Θ defined as Θ(X) = σ 1z Xσ 1z for X SU( n ), where σ jz denotes that the operation defined as the Pauli matri σ z acts on the j-th qubit. He also introduced a method that translates matrices U, Σ, and V in (1) into k 1 ep(k), a ep(h), and k ep(k), respectively, where k and h are the ones defined in KGD. Here, KGD correspondsto G = KAK decomposition with the selection of Θ defined as Θ(X) = σ nz Xσ nz. We can thus produce a KGD-based quantum circuit by combining Bullock s translation and the CSD-based algorithms. We introduce a new algorithm that translates a n n unitary matri into a quantum circuit according to the G = KAK theorem. The algorithm can derive any matri decomposition corresponding to type-aiii KAK decompositions for the given global Cartan involution Θ. The algorithm contains, as its special cases, both CSD and KGD in the sense that it derives the same quantum circuits as the ones calculated by them if we select suitable Cartan involutions and square root matrices. Here, we select Θ(X) as σ 1z Xσ 1z for CSD and as σ nz Xσ nz for KGD, where X SU( n ). The strategy utilized in our algorithm is related to those used in Refs. [1, 13]. However, those strategies provided methods for computing type- AII KAK decomposition; no translation between type-aii decompositions and type-aiii decompositions was provided. Furthermore, the method utilized in Ref. [1] is different from ours in the square root matri calculations, i.e., methods for calculating m from m (where g = km is a global Cartan decomposition of the input matri g). In the method proposed in Ref. [1], first, a square root matri is calculated in Lie algebralevel. And then it is translated into an element in Lie group level via eponential mapping. In contrast, with our method, a square root matri is calculated directly at the Lie group level. Although our algorithm contains CSD and derives any matri decomposition correspond-

3 Yumi Nakajima, Yasuhito Kawano, and Hiroshi Sekigawa 3 ing to type-aiii KAK decompositions according to the given Cartan involution, the efficiency for calculating a decomposition is not sacrificed. The reason is as follows: Roughly speaking, to decompose g into k1ak, the CSD-based algorithms apply SVD to four n 1 n 1 matrices (g 11, g 1, g 1, and g in Eq.(1)), while our algorithm applies eigenvalue decomposition to n n matri. Therefore, the efficiencies for computing SVD on four n 1 n 1 matrices and for computing eigenvalue decomposition on one n n matri are the same. In addition, our algorithm might have an advantage over CSD when we would like to determine a class of quantum circuits for a given class of matrices. The reason is as follows: In CSD-based algorithms, it is difficult to formulate a class of matrices u 1, u, v 1, and v such that relation (1) holds for a given class of input matrices. Actually, to reproduce the well-known QFT circuit by using CSD [,14], Tucci changes the rows and columns of the QFT matrices beforehand and makes each submatri hold a convenient form, which can be written by the (n 1)-qubit QFT. It would not be possible to describe the general form of the decomposition when the input matri does not have a convenient form like QFT. On the other hand, our algorithm does not require such a preliminary change of rows and columns. All matrices appearing through our algorithm can be described using the input matri g, the given global Cartan involution Θ, and the eigenvalues and eigenvectors of these matri products. This will be shown eplicitly as an eample of the QFT decomposition (Section 4). The paper is organized as follows: In the following section, we cover some preliminaries about notations, the G = KAK theorem, and KGD. Section 3 presents our algorithm for computing the decomposition follows from the G = KAK theorem. Section 4 presents decompositions of the n-qubit QFT using our algorithm and CSD-based algorithms. We show that we can produce the well-known QFT circuit by using these matri decompositions. Preliminaries.1 Notations Let σ, σ y, andσ z denotethe Paulimatricesand I s be a s s identity matri(i = 1 1 ). We use σ jα to denote the Pauli matri acting on the j-th qubit; σ jα = I (j 1) σ α I (n j), (α =,y, or z). Let U CNOT denote the standard CNOT gate, H denote a Hadamard gate, and R (ζ) = ep( iζσ ). All these notations follow those in Ref. [1].. G = KAK theorem The G = KAK theorem for compact groups (cf. [9], Theorem 8.6) provides a framework for decomposing g SU( n ) into the following matri products: g = k 1 ak, k 1,k ep(k),a ep(h) ep(m). () Here, su( n ) = k m is a Cartan decomposition, where k and m = k are orthogonal vector spaces, and h is a Cartan subalgebra, that is, a maimal Abelian subalgebra contained in m. Let θ denote the Cartan involution of its Lie algebra su( n ); i.e., (i) θ = I n (θ I n ) and (ii) θ is an automorphism of the Lie algebra su( n ). And let the global Cartan involution (cf. [15], p. 36) of SU( n ) be Θ. Then k and m have the following property: θ() = { if k if m, Θ(X) = { X if X ep(k) X if X ep(m). (3)

4 4 A new algorithm for producing quantum circuits using KAK decompositions k and ep(k) m Fig. 1. Patterns of an element of k, ep(k) and m for a three-qubit system, where each square represents an element of an 8 8 matri. The white elements are always zero, and the black elements take some value that depends on the input matri. Three types of k-algebra, named AI, AII, and AIII, arise for su( n ). Here, AI, AII, and AIII correspondto k = so( n ), k = sp( n ), and s[u(p) u(q)] (p+q = n ), respectively (cf. [9], p. 518)..3 Khaneja-Glaser decomposition (KGD) Khaneja and Glaser provided a particular selection of bases of k, m, and h (cf. [8], Notation 3 and 5) so that the selection matches an NMR system. Notice that we use h instead of h(n). Here, generators of k, m, and h, are denoted to be tensor products of the Pauli matrices; k = span{a σ z /,B I,iσ nz / A,B su( n 1 )}, (4) m = span{a σ /,B σ y /,iσ n /,iσ ny / A,B su( n 1 )}. (5) Here, generators of k and m have a specific operation on the last qubit; i.e., σ z or I for generators of k and σ or σ y for generators of m. Thus, to determine k and m, the Cartan involution θ and the global Cartan involution Θ can be chosen as follows: θ() = σ nz σ nz, Θ(X) = σ nz Xσ nz. (6) Since θ(σ z ) = σ z, θ(i) = I, θ(σ ) = σ, θ(σ y ) = σ y, we can check that the above k and m satisfy relation (3) when θ and Θ are chosen as (6). For the number of qubits n 3, k and m have specific patterns, as shown in Fig. 1, because all generators defined in (4) and (5) have these patterns. Note that in contrast to an element of ep(k) taking the same pattern as an element of k, an element of ep(m) does not take the same pattern as m. This property enables us to apply the KAK decomposition recursively, as shown in Fig.. 3 Our algorithm 3.1 Basic ideas We provide a new constructive algorithm that computes a decomposition based on the G = KAK theorem. Here, we choose a Cartan subalgebra h, which is different from the h used in KGD. Since Cartan subalgebras are Abelian, they can translate each other by h = AdT ep(k) (h). Here, T is fied for given h and h. It should be noted that G = KAK decomposition for the fied input g SU( n ) is not unique. We provide the following theorem.

5 Yumi Nakajima, Yasuhito Kawano, and Hiroshi Sekigawa 5 Theorem 1 Let g SU( n ) be the input matri. If g has a global Cartan decomposition g = km (k ep(k),m ep(m)), then m is uniquely determined by m = Θ(g )g. Proof. From (3), Θ(g )g = Θ(m k )km = Θ(m )Θ(k )km = mk km = m. Theorem 1 shows that the fied global Cartan involution Θ only determines m. Therefore, arbitrariness remains in the selection of m, and also k. Furthermore, a in () has also arbitrariness because it follows from a decomposition of m, m = k a k, where k ep(k). (See, [15], AII.3). Using Theorem 1, we compute the decomposition in () as follows: 1. Compute m = Θ(g )g.. Decompose m = pbp such that p ep(k) and b ep( h). Such decomposition always eists because m ep(m) (cf. [15], Proposition 7.9). We compute p using eigenvectors of m. We show eamples in Section Find y such that y = b and y ep( h). The y can be computed by replacing the diagonal blocks of b, when we choose a suitable h. We show eamples of selections of h in Section Compute m = pyp. Here, m ep(m) because Θ(m) = Θ(p)Θ(y)Θ(p ) = py p = m. 5. Compute k = gm. Then, k always satisfies k ep(k) because (m ) = g Θ(g) and Θ(gm ) = Θ(g)Θ(m ) = g(m ) m = gm. Steps 4 provide a method for computing the square root of a matri to find m from m. After these procedures, we obtain g = kpyp = kyp. (7) Here, Θ( k) = k, Θ(y) = y, and Θ(p ) = p, so that the decomposition follows the G = KAK theorem. All matrices that appear through the algorithm can be described using g and Θ. This enables us to show that algorithm can automatically reproduce the well-known QFT circuit (See, Section 4.). 3. Methods for performing steps and 3 based on KGD Above, we provided a method for computing G = KAK decomposition. We did not provide a concretemethodforcomputingp, b, andy instepsand3inouralgorithm. Toshoweamples of such concrete methods, we fied Θ and h. Here, we treat eamples that compute KGD. We choose Θ as in (6) and show two methods that compute p, b, y for particular selections of h; i.e., h 1 = span { j j iσ j = 0,, n 1 1}and h = span{ j j i(σ σ +σ y σ y ) j = 0,, n 1}. We should note that the second selection is not always possible because it demands that all eigenvalues appearing in m in step 1 should be duplicated twice Eample 1 Let Θ be as in (6), and h 1 as span { j j iσ j = 0,, n 1 1}. Then, we compute p, b, and y in steps and 3 as follows:

6 6 A new algorithm for producing quantum circuits using KAK decompositions (i) Compute eigenvalue decomposition of m. Let m = pd p be eigenvalue decomposition and µ 1,,µ N be the columns of p, where N = n. Then, d is a diagonal matri in which diagonals have eigenvalues of m and all µ j s are eigenvectors of m and mutually orthogonal. (ii) Normalize all µ j s by ν j 1 = µ j +σ nz µ j µ j +σ nz µ j, ν j = µ j σ nz µ j µ j σ nz µ j. (8) (iii) For all ν j s that are associated with imaginary eigenvalues, (a) let W 1, W, and W 3 be sets of vectors such that W 1 = {ν j σ nz ν j = ν j }, W = {ν j σ nz ν j = ν j }, W 3 = {ν j σ nz ν j ±ν j }. (b) For each w W 3, compute ν + = w+σ nzw w+σ nz w, ν = w σ nzw w σ nz w. (9) Here, denotes the length of a vector. Then, if all elements in W 1 and u + are linearly independent, then W 1 = W 1 {ν + }; if all elements in W and u are linearly independent, then W = W {ν }. (iv) Repeat steps (a) and (b) for all µ j s that are associated with positive real eigenvalues. (v) Repeat steps (a) and (b) for all µ j s that are associated with negative real eigenvalues. (vi) Let p = (υ 1,υ,,υ N ), where υ j 1 W 1 and υ j W, for j = 1,,N/. The computation procedure follows from Appendi A. Since σ nz υ j 1 = υ j 1 and σ nz υ j = υ j, we can easily check that Θ(p) = p. Then, b = p m p satisfy Θ(b) = b and b ep( h), where b = n 1 j=0 j j R (ζ j ) (0 ζ j < π). (vii) Compute y by replacing all R (ζ j ) appears in b with R (ζ j ). Since Θ(R (ζ j )) = R (ζ j), R (ζ j) = R (ζ j ), then y satisfies y = b and y ep( h 1 ). In step (vi), one may notice that, when R (π) appears in b, then we can use R y (π/) instead of R (π/) as a replacement rule. Since k and p are elements of ep(k) that has the specific pattern as shown in Fig. 1, they have the following decomposition: k = g (0) g (1) 1 1 1, p = g (0) 0 0 +g (1) 1 1, (10) where g (j) 1,g(j) SU( n 1 ) for k,p SU( n ) (j = 0 or 1). Here, g (0) 1 and g (0) are composed of nonzero elements (black squares in Fig. 1) of odd rows, and g (1) 1 and g (1) are composed of nonzero elements of even rows.

7 Yumi Nakajima, Yasuhito Kawano, and Hiroshi Sekigawa 7 g = k ~ y p g ( j) R ( ζ j ) g ( j) 1 Fig.. Image of a decomposition when we choose Θ as in (6) and a Cartan subalgebra h 1 as span { j j iσ j = 0,, n 1 1}: Thematricesusesthe same notation asinfig.1to representthe properties. In the quantum circuit, SU(4) (l,j = 1 or ) are applied selectively; that is, g (0) 1 and g (0) are applied when the third qubit is 0, the symbol ofthe control qubit represents the uniformlycontrolled rotations [5,16]. g (j 1) l whereas g (1) 1 and g (1) are applied when it is 1. Fig.showstheimageofadecompositionin()forathree-qubitsystem; thatis, wechoose Θ as in (6) and a Cartan subalgebra as h 1. Applying the decomposition in (7) recursively to elements g (j) 1,g(j) SU( n 1 ) (j=0 or 1), we obtain a sequence of uniformly controlled rotations like in Fig. 13 in Ref. [5], ecept that R is used instead of R y in our case. The full decomposition of these uniformly controlled rotations into elementary gates has been provided by Möttönen et al. [5,16]. Also, if we change the order of qubits and apply the quantum Multipleor decomposition to k and p in (), the produced circuit is the same as that in Fig. in Ref. [6], ecept that rotation R y is used instead of R. Therefore, the number of elementary gates needed to compose g SU( n ) in our method is O(4 n ), which is the same as in Refs. [4,6]. 3.. Eample Here, we show another eample of methods for computing p, b, and y in Section 3.1. Here, we choose h = span { j j i(σ σ + σ y σ y ) j = 0,, n 1 1}. The Cartan involution Θ is the same as in (6). The decomposition of this type is chosen as an eample of a decomposition of the QFT. In this case, we only change steps (vi) and (vii) in Section 3..1 as follows: (vi) Let p = (υ 1,υ,,υ N ), for j = 1,,N/4, ν 4j 3 W 1 and it is associated with positive eigenvalues, ν 4j W and it is associated with negative eigenvalues, ν 4j 1 W 1 and it is associated with negative eigenvalues, ν 4j W and it is associated with positive eigenvalues. Then, p also satisfies Θ(p) = p, and b = p m p is a block-diagonal matri, in which the diagonals are constructed from the 4 4 matri cos(ζ j ) isin(ζ j ) 0 0 isin(ζ j ) cos(ζ j ) 0 = ep(iζ j(σ σ +σ y σ y )) Here, the middle part of the above matri is R (ζ j ).

8 8 A new algorithm for producing quantum circuits using KAK decompositions (vii) Compute y by replacing R (ζ j ) in b with R (ζ j ). 4 Decompositions of the QFT using G = KAK decomposition In this section, we show that we can automatically reproduce the well known QFT circuit using our method. All the matrices that appear through our algorithm can be described using the input matri g and Θ. In contrast, it is difficult to describe all the matrices that appear through the CSD-based algorithm because, as shown in (1), the input matri g has to be divided into four square matrices and SVD has to be applied to each partitioned matri. Furthermore, we have to choose a suitable decomposition for each partitioned matri g 11, g 1, g 1, and g such that Eq. (1) holds, which makes it difficult to formulate U, Σ, and V in (1). Fortunately, the n-qubit QFT is a very special matri that has the following property: If we permute the order of qubits, then each partitioned matri can be described using (n 1)-qubit QFT. Using the feature, we provide a decomposition of the QFT by CSD-based algorithm. This is shown in Section Notation The QFT on n qubits, F n, is a n n matri such that F n = ( ) ( ) 1 πi n ω(j 1)(l 1) n, where ω n = ep jl n. (11) We define Q n as a n n permutation matri: Q n = χ n n 1 χn χn 1, where χk j is the SWAP gate applied to the j-th and the k-th qubits. Let H 1 = H I (n 1), then (11) is written as F n = 1 ( ) Fn 1 Ω n 1 F n 1 Q F n 1 Ω n 1 F n = H 1 D n (I F n 1 )Q n, (1) n 1 where ( ) I (n 1) 0 D n =, Ω 0 Ω n 1 = diag (1,ω n,,ω n 1 1 n ). n 1 This notation follows from Section in Ref. [17]. 4. Decomposition of the QFT by our method Following Section 3.1, we compute a decomposition as follows: 1. Compute m = Θ(F n )F n. Let S be (I F n 1 )Q n. Since Θ(H 1 ) = H 1, Θ(D n ) = D n, and σ 1z S = Sσ nz, we have m = S σ 1z σ nz S. All column vectors of S are then eigenvectors of m because σ 1z σ nz is a diagonal matri in which diagonal elements are eigenvalues of m.. Decompose m = pbp such that Θ(p) = p and Θ(b) = b. We define p = S Q n (This selection was done so that p satisfies Θ(p) = p and follows Section 3..), then b = p m p = Q n σ 1zσ nz Q n = I (n ) diag(1, 1, 1, 1).

9 Yumi Nakajima, Yasuhito Kawano, and Hiroshi Sekigawa 9 3. Choose y such that Θ(y) = y. Following the step (vi) in Section 3.., we have y = I (n ) 0 0 i 0 0 i 0 0 = I (n ) ep(π(σ σ +σ y σ y )/4) This is obtained by replacing each 4 4 diagonal block of v, b = diag(1, 1, 1, 1), with y = ep(π(σ σ +σ y σ y )/4). We can easily check that Θ(y) = y and y = b, because Θ(y) = y and y = b. 4. Compute m = pyp. 5. Compute k = gm = F n m. Then, k = kp = (F n py p )p = F n py, so we have the following decomposition: F n = kyp = (H 1 D n Q n y )y(q n S) = H 1D n S = H 1 D n (I F n 1 )Q n. (13) We apply a similar decomposition to F j, for j = n 1,n,,. Net, we show a decomposition of D n. D n is controlled-ω j (where j = n 1,n,,), so it suffices to consider the decomposition of Ω j. Since Ω j ep(k) (it follows from Θ(Ω j ) = Ω j ), we apply the decomposition in (10) to Ω j. Consider Ω 3 as an eample, then we have Ω 3 = g (0) g (1) 1 1 1, where g (0) 1 = diag (1,ω,ω 4,ω 6 ) and g (1) 1 = diag(ω,ω 3,ω 5,ω 7 ) = ω diag (1,ω,ω 4,ω 6 ). Then, we have Ω 3 = diag (1,ω,ω 4,ω 6 ) diag (1,ω). Similarly, since diag(1,ω,ω 4,ω 6 ) is also an element ofep(k), it is decomposed into diag(1,ω 4 ) diag(1,ω ). Therefore, Ω n 1 is composed of n 1 single-qubit rotations as follows: ( ) ( ) ( ) ( ) Ω n 1 = 0 ω n 1 1 n 0 ωn j 1 0 ωn. (14) 0 ω n The circuit obtained from the above decomposition is shown in Fig. 3. We apply a similar decomposition to F j (for j = n 1,n,,). Finally, we have a full QFT decomposition composed of n Hadamard gates, 1 n(n 1) controlled-rotations, and n SWAP gates. Here, SWAP gates that appeared in a sequence of permutations Q n Q n 1 Q were optimized. It is known that a controlled-rotations can be implemented by three single-qubit rotation and two CNOTs, so that the number of elementary gates that appear in Fig. 3 is O(n ). 4.3 Decomposition of the QFT by CSD If we apply Q n to the input matri F n beforehand, then we have the following decomposition: F n Q n = UΣV ( ) ( ) I n I (n 1) I = 0 I n 1 (n 1) I (n 1) I (n 1) ( ) Fn 1 0,(15) 0 Ω n 1 F n 1 = (Hσ z I (n 1) )(σ z I (n 1) )D n (I F n 1 ) = H 1 S. (16) We can easily check that (15) satisfies the definition of CSD, where all θ j s appear in c and s in (1) are π/4. Although we have to apply Q j beforehand when we apply CSD recursively to F j (j = 1,,n 1), we obtain the well-known QFT circuit. Using the feature, Tucci reproduced the well-known QFT circuit using a CSD-based algorithm [, 14].

10 10 A new algorithm for producing quantum circuits using KAK decompositions 1 R n H F = n F n 1 n 3 R n n R n n 1 R n Q n D n Fig. 3. Decomposition of the QFT by our algorithm. Here, R j n = diag (1,ω j n), i.e., a single-qubit rotation, and Q n is composed of n 1 SWAP gates. 4.4 Decomposition of the QFT by the quantum Shannon decomposition (QSD) The QSD [6] is a method that combines CSD and the quantum multipleor decomposition. First, we compute CSD for the input matri and then apply the quantum multipleor decomposition (cf. [6], Theorem 1), ( ) u1 0 0 u = ( )( )( ) v 0 d 0 w 0 0 v 0 d, (17) 0 w to U and V in (1). Here, u 1 u = vd v and w = dv u. In the QFT, we apply the decomposition in (17) to V because U is an identity matri in (15). Then, we have v = I (n 1) and d = Ω n 1, and w = Ω n 1 F n 1 in (17). Here, Ω n 1 is a n 1 n 1 diagonal matri whose (j,j)-th component is iω j 1 n+1, (j = 1,,n 1 ). Therefore, the decomposition of QFT by the QSD is as follows: F n Q n = ep(iσ y δ )ep( iσ z δ 3 )(I v 4 ), (18) where v 4 = Ω n 1 F n 1 and δ and δ 3 are n 1 n 1 diagonal matrices. Each element of δ is π/ n 1 and each (j,j)-th element of δ 3 is jπ/ n+1. Furthermore, Ω n 1 is composed of n 1 single-qubit rotations as follows: ) ( ) ( ) ( ) Ωn 1 = i( 0 ωn+1 n 1 0 ωn+1 j 1 0 ωn+1. 0 ω n+1 Here,(18) is also equal to the well-known QFT decomposition(1) after optimization as shown in Fig. 4. To obtain the decomposition in (18), note that we have to apply Q j (j = 0,,n) beforehand for F j. 5 Conclusion We introduced a new algorithm for computing any type-aiii KAK decomposition according to the given global Cartan involution Θ. Recursively performing the decomposition leads us to a quantum circuit composed of uniformly controlled rotations. The algorithm can derive any matri decomposition corresponding to the type-aiii KAK decomposition, and it contains

11 Yumi Nakajima, Yasuhito Kawano, and Hiroshi Sekigawa 11 σ z H H F n = 1 D D D = 1 Ω n 1 ep( I v4) ep( iσ ) y δ ep( σ δ3) Fig. 4. Decomposition of the QFT by the QSD. Here, we move Q n to the right-hand-side by inverting it. The dark-gray block in the middle circuit can be simplified. We show the simplified block by the same color in the rightmost circuit. Here, D = Ω n 1. CSD and KGD as its special cases. This is because our algorithm contains arbitrariness in selecting the Cartan subalgebra h and a square root matri m for the given Cartan involution Θ, where m is a matri derived from a global Cartan decomposition g = km. We also showed two methods for computing a square root matri. Although the correctness of our algorithm depends on Lie group theory, the main methods involved are eigenvalue decomposition and a simple replacement rule. Thus, we can compute a decomposition without knowledge of Lie group theory. As an eample, we showed that our method automatically reproduces the well-known QFT circuit for arbitrary n-qubits. When using CSD-based algorithms, we have to apply some permutations beforehand in order to reproduce the circuit. The same technique can not always be used to describe the canonical form of the decomposition for a given matri because matrices do not always have a convenient form like QFT. Our algorithm might be useful in showing the effectiveness of G = KAK matri decompositions for other particular input matrices, because all matrices appearing through our algorithm can be described using a given input matri g and a given Cartan involution Θ. References 1. M. Nielsen, I. Chuang, Quantum Computation and Quantum Information, (Cambridge University Press, Cambridge, England, 000).. R. Tucci, A rudimentary quantum compiler, quant-ph/ K. Svore, Compiling quantum circuits into elementary unitary operations, QIP004, Waterloo, Canada (004). 4. M. Möttönen, J. Vartiainen, V. Bergholm, and M. Salomaa, Quantum circuits for general multiqubit gates, Phys. Rev. Lett. 93, (004) quant-ph/ M. Möttönen, and J. Vartiainen Decompositions of general quantum gates, quant-ph/ V. Shende, S. Bullock, and I. Markov, Synthesis of quantum logic circuits, quant-ph/ , to appear in IEEE Trans. on Computer-Aided Design of Integrated Circuits and Systems. 7. A. Barenco et al., Elementary gates for quantum computation, Phys. Rev. A 5, pp (1995) quant-ph/ N. Khaneja and S. Glaser, Cartan decomposition of SU( n ), constructive controllability of spin systems and universal quantum computing, Chem. Phys. 67, pp (00)

12 1 A new algorithm for producing quantum circuits using KAK decompositions quant-ph/ S. Helgason, Differential geometry, Lie groups, and Symmetric spaces, (Academic Press, 1978). 10. N. Khaneja, R. Blockett, and S. Glaser, Time optimal control in spin systems, Phys. Rev. A 63, (001) quant-ph/ S. Bullock, Note on the Khaneja Glaser decomposition, Quantum Inf. Comput. 4, pp (004) quant-ph/ S. Bullock, G. Brennen, and D. P. O Leary, Time reversal and n-qubit canonical decomposition, J. Math. Phys. 46, (005) quant-ph/ S. Bullock and I. Markov, Arbitrary two-qubit computation in 3 elementary gates, Phys. Rev. A 68, (003) 14. R. Tucci, Quantum fast Fourier transform viewed as a special case of recursive application of Cosine-Sine decomposition, quant-ph/ A. W. Knapp,Lie groups beyond an introduction, nd ed., (Birkhäuser, 005). 16. M. Möttönen, J. Vartiainen, V. Bergholm, and M. Salomaa, Transformation of quantum states using uniformly controlled rotations, Quantum Inf. Comput. 5, pp (005) quant-ph/ G. Golub and C. Van Loan, Matri computations, 3rd ed., (The Johns Hopkins University Press, 1996). Appendi A Properties of eigenvalues of m = Θ(G )G In Section 3..1, we show that b is a block-diagonal matri, each block of which is R (ζ), where h = span{ j j iσ j = 0,, n 1 1} m. This appendi provide two lemmas and proofs for the eigenvalues of m = Θ(G )G, where G SU( n ) to prove that. Let X = m = Θ(G )G. Lemma A.1 (Properties of comple eigenvalues of X) 1. The number of the comple eigenvalues of X, with multiplicity counted, is even. They are of the form α 1, α 1, α, α,..., α t, α t, repeated with multiplicity.. There eists unit length vectors u 1,..., u t that are mutually orthogonal and satisfy σ nz u j 1 = u j 1 and σ nz u j = u j. Furthermore, there eists ζ 1,..., ζ t R satisfying Xu j 1 = cos(ζ j )u j 1 isin(ζ j )u j and Xu j = isin(ζ j )u j 1 +cos(ζ j )u j. Proof. Let α 1,..., α n be the eigenvalues of X, repeated with multiplicity. Then, ecept for the order, the elements of the sequence α 1,..., α n are equal to α 1,..., α n because Θ(X) = X. Therefore, the number of the comple eigenvalues, counted with multiplicity, is even, and we denote it as t. We can write the comple eigenvalues as α 1, α 1,...,α t, α t. That is, the first statement holds. Let ube aneigenvectorofx correspondingtothe eigenvalueα, then σ nz uisaneigenvector of X corresponding to the eigenvalue α because X(σ nz u) = σ nz X u = α j (σ nz u). Now, let W be the eigenspace corresponding to a comple eigenvalue α and β 1,..., β r be an orthonormalbasisofw. Then, letσ nz β 1,..., σ nz β r be anorthonormalbasisoftheeigenspace corresponding to the eigenvalue α. Therefore, u j 1 = β j +σ nz β j and u j = β j σ nz β j for j = 1,,r, spantheeigenspaceofx correspondingtotheeigenvaluesαandα. Furthermore, they are eigenvectors of σ nz because σ nz u j 1 = u j 1 and σ nz u j = u j. On the other

13 Yumi Nakajima, Yasuhito Kawano, and Hiroshi Sekigawa 13 hand, the following relations hold: Xu j 1 = αβ j +ασ nz β j = α+α Xu j = αβ j ασ nz β j = α α u j 1 + α α u j = Re(α)u j 1 +iim(α)u j, u j 1 + α+α u j = iim(α)u j 1 +Re(α)u j. Therefore, put ζ = arg(α). Then cos(ζ) = Re(α) and sin(ζ) = Im(α). That is, Xu j 1 = cos(ζ)u j 1 isin(ζ)u j, Xu j = isin(ζ)u j 1 +cos(ζ)u j. Similar arguments hold for the other comple eigenvalues. Lemma A. (Properties of real eigenvalues of X) Real eigenvalues of X are ±1. 1. The multiplicity of the eigenvalue 1 is even. There eists an orthonormal basis u t+j (j = 1,,µ) of the eigenspace corresponding to the eigenvalue 1 that satisfies the σ nz u t+j 1 = u t+j 1 and σ nz u t+j = u t+j.. The multiplicity of the eigenvalue 1 is even. There eists an orthonormal basis u t+µ+j (j = 1,,ν) of the eigenspace corresponding to the eigenvalue 1 that satisfies σ nz u t+µ+j 1 = u t+µ+j 1 and σ nz u t+µ+j = u t+µ+j. Proof. Lemma A.1 implies that the product of all the comple eigenvalues is 1, and thus the product of all the real eigenvalues is 1. Therefore, the real eigenvalues of X are 1 or 1, thus both of the multiplicities of the eigenvalues 1 and 1 are even. Let W 1 and W be the eigenspaces of X corresponding to the eigenvalues of 1 and 1, respectively. Then σ nz W 1 W 1 and σ nz W W hold because Θ(X) = X. Put W = W 1 W. Then W is the direct sum of the eigenspaces for the comple eigenvalues of X. The trace of σ nz W is 0 because Lemma A.1 implies that the trace of σ nz W is 0, and this implies that the multiplicities of the eigenvalues of 1 and 1 of σ nz W are equal. We write W 1 = W 11 W 1 and W = W 1 W, where W 11 and W 1 are the eigenspaces of σ nz corresponding to the eigenvalue 1, and W 1 and W are the eigenspaces of σ nz corresponding to the eigenvalue 1. Let dimw ij be d ij. Then we have d 11 +d 1 = d 1 +d. We can make similar arguments forx = Θ(G)G. We write W 1 = W 11 W 1 and W = W 1 W, where W 11 and W 1 are the eigenspaces of σ nz corresponding to the eigenvalue 1, and W 1 and W are the eigenspaces of σ nz corresponding to the eigenvalue 1. Let dimw ij be d ij. Then we have d 11 +d 1 = d 1 +d. From Lemma A.3 below, we have d 11 = d 11, d 1 = d 1, d 1 = d, and d = d 1. Then, we have d 11 = d 1 and d 1 = d. Therefore, the statements of the lemma hold. Lemma A.3 Let X = Θ(G )G and X = Θ(G)G. 1. If Xu = u and σ nz u = u, then X (Gu) = Gu and σ nz (Gu) = Gu.. If Xu = u and σ nz u = u, then X (Gu) = Gu and σ nz (Gu) = Gu. 3. If Xu = u and σ nz u = u, then X (Gu) = Gu and σ nz (Gu) = Gu.

14 14 A new algorithm for producing quantum circuits using KAK decompositions 4. If Xu = u and σ nz u = u, then X (Gu) = Gu and σ nz (Gu) = Gu. Proof. First, we prove the statements for the eigenvalues and eigenvectors of X. Since G X G = Θ(X) holds, we have (X (Gu),Gu) = (G X Gu,u) = (Θ(X)u,u) = (X u,u) = (u,xu). Thus, the equation Xu = ǫu, where ǫ = ±1, implies (X (Gu),Gu) = ǫ(u,u). On the other hand, the Cauchy-Schwarz inequality implies (X (Gu),Gu) X (Gu) Gu. Since the right-hand side is equal to Gu = u = (u,u) = (X (Gu),Gu), the equality (X (Gu),Gu) = X (Gu) Gu holds. Therefore, we have X (Gu) = α(gu) for some α C; that is, Gu is an eigenvector of X corresponding to the eigenvalue α. The equality α = ǫ follows ǫ(u,u) = (X (Gu),Gu) = (αgu,gu) = α(gu,gu) = α(u,u). To prove the statements for the eigenvalues and eigenvectors of σ nz, we use the following equations: (σ nz (Gu),Gu) = (G σ nz Gu,u) = (σ nz Xu,u) = (Xu,σ nz u). The equations Xu = ǫu and σ nz u = ǫ u, where ǫ, ǫ = ±1, imply (σ nz (Gu),Gu) = ǫǫ (u,u). From similar arguments for the eigenvalues and eigenvectors of X, we have σ nz (Gu) = ǫǫ (Gu).