arXiv:1405.5749v2 [math-ph] 2 Jun 2014 +1 etePuisi arcs einclude We matrices. spin Pauli the be hntesi arcsaegvnby given are matrices spin the Then h al pnmtie r ntr n emta ihegnaus+ a +1 eigenvalues with hermitian and unitary are matrices spin Pauli The hi rnce rdc.Let spin Pauli product. of Kronecker anticommutators their and commutators the investigate We Introduction 1 provided. are u operators mutually Bose operators, and Hamilton operators spin Fermi matrices, such pro eigenvalue of to functions Applications tial studied. are product Kronecker the Abstract. h matrices The σ omttr ftePuisi arcsvns,i.e. vanish, matrices spin Pauli the of commutators 3 σ / pnOeaos al ru,Cmuaos Anti-Commutato Commutators, Group, Pauli Operators, Spin 1 and 2 = iσ nvriyo oansug ukadPr 06 ot Africa, South 2006, Park Auckland Johannesburg, of University 2 − n h omttr r ie by given are commutators the and al pnmtie,Puigop omttr,anti-commutator commutators, group, Pauli matrices, spin Pauli 1 / iσ σ .FrtePuisi arcsw n that find we matrices spin Pauli the For 2. nvriyo ot fia rtra ot Africa, South Pretoria, Africa, South of University 1 1 [ σ , † = [ 1 σ iσ nentoa colfrSinicComputing, Scientific for School International σ , il-asSteeb Willi-Hans 1  rnce rdc n Applications and product Kronecker 2 σ , ∗ , 2 0 1 1 0 eateto ahmtclSciences, Mathematical of Department 2 = ] 2 iσ ] + 3 e-mail: e-mail: 0 = omabsso h ipeLealgebra Lie simple the of basis a form  iσ σ , 3 2 , , [email protected] 2 [ [email protected] σ [ σ = 2 2 s σ , σ , 1  3 = † 3 0 σ 2 = ] i ] n oikHardy Yorick and 0 + 2 1 = σ 0 = − 0 1 iσ i , I 2  s 2 1 where , 2 , , σ , = [ [ σ σ 1 2 σ 3 3 3 σ , σ , 2 = , I 2 1 1 s  ] 2 = ] 3 ste2 the is + σ = 0 0 1 0 = 1 σ ∗ 2 1 iσ 2 − σ 2 = 3 2 1 . su × iheigenvalues with lm,exponen- blems,  iσ bae bases, nbiased ntmatrix. unit 2 2.Teanti- The (2). 3 arcsand matrices , σ 2 σ 3 nd rs, = and s iσ − 1. 1 , where 0 is the 2 2 zero matrix. Here we study the commutators and anticommu- 2 × tators of the Kronecker product of the Pauli spin matrices [1]. Thus we study ⊗ the commutator and anticommutator of the 2n 2n unitary matrices of the form × n j0 ( i) σj − t t=1 O where j 0, 1, 2, 3 and jt 0, 1, 2, 3 . These matrices are elements of the 0 ∈ { } ∈ { } Pauli group [2], [3]. Furthermore the square of the matrices

n

σjt t=1 O is the 2n 2n unit matrix. Then × n n 1 1 Π = I n + σj , Π = I n σj 1 2 2 t 2 2 2 − t t=1 ! t=1 ! O O are 2n 2n projection matrices. Whether the commutator or anticommutator of × such matrices vanishes is helpful for the eigenvalue problem of such matrices and for the calculation of the exponential function of such matrices. Other applica- tions discussed concern spin-Hamilton operators and the projection to sub-Hilbert space and mutually unbiased bases. Finally application with Fermi operators are described.

2 Kronecker Product of Pauli Spin Matrices

Let us first give some examples where the commutator or the anticommutator vanishes. Consider first the 4 4 matrices × σ := σ σ , σ := σ σ , σ := σ σ . 12 1 ⊗ 2 23 2 ⊗ 3 31 3 ⊗ 1 Then we find that the commutators vanish, i.e.

[σ12, σ23]=04, [σ23, σ31]=04, [σ31, σ12]=04 where 0 is the 4 4 zero matrix. For the anticommutators we find the hermitian 4 × invertible matrices which can expressed as Kronecker products of the Pauli spin matrices

0 1 0 0 − 1 0 00 [σ12, σ23]+ =2  −  = 2σ3 σ1 = 2σ31 0 0 01 − ⊗ −  0 0 10      0 0 0 1 0 0 1 0 [σ23, σ31]+ =2i  −  = 2σ1 σ2 = 2σ12 0 1 0 0 − ⊗ −  10 0 0   −   0 01 0  0 00 1 [σ31, σ12]+ =2i  −  = 2σ2 σ3 = 2σ23. 100 0 − ⊗ − −  0 10 0      Consider now the three hermitian and unitary 8 8 matrices × σ := σ σ σ , σ := σ σ σ , σ := σ σ σ . 123 1 ⊗ 2 ⊗ 3 312 3 ⊗ 1 ⊗ 2 231 2 ⊗ 3 ⊗ 1 The commutators are non-zero and we find the skew-hermitian invertible 8 8 × matrices

02 02 σ1 02 02 02 02 σ1 [σ123, σ312]=2  −  =2iσ2 σ3 σ1 =2iσ231 σ1 02 02 02 ⊗ ⊗ −  0 σ 0 0   2 1 2 2    02 02 02 σ3 02 02 σ3 02 [σ312, σ231]=2  −  =2iσ1 σ2 σ3 =2iσ123 02 σ3 02 02 ⊗ ⊗  σ 0 0 0   − 3 2 2 2    02 σ2 02 02 σ2 02 02 02 [σ231, σ123]=2i   =2iσ3 σ1 σ2 =2iσ312. 02 02 02 σ2 ⊗ ⊗ −  0 0 σ 0   2 2 − 2 2    However the anti-commutators vanish, i.e.

[σ123, σ312]+ =08, [σ312, σ231]+ =08, [σ231, σ123]+ =08.

As mentioned above for the Pauli spin matrices σ1, σ2, σ3 we find that the anti- commutators vanish and the commutators are given by [σ1, σ2]= iσ3,[σ2, σ3]= iσ1,

[σ3, σ1]= iσ2. Consider now the three unitary and hermitian matrices

σ := σ σ , σ := σ σ , σ := σ σ . 11 1 ⊗ 1 22 2 ⊗ 2 33 3 ⊗ 3 Then the commutators vanish, i.e.

[σ11, σ22]=04, [σ22, σ33]=04, [σ33, σ11]=04 and the anticommutators can be written as Kronecker products, i.e.

[σ , σ ] = 2σ σ = 2σ 11 22 + − 3 ⊗ 3 − 33 [σ , σ ] = 2σ σ = 2σ 22 33 + − 1 ⊗ 1 − 11 [σ , σ ] = 2σ σ = 2σ . 33 11 + − 2 ⊗ 2 − 22 Now consider the 8 8 hermitian and unitary matrices × σ := σ σ σ , σ := σ σ σ , σ := σ σ σ . 111 1 ⊗ 1 ⊗ 1 222 2 ⊗ 2 ⊗ 2 333 3 ⊗ 3 ⊗ 3 Here the anticommutators vanish, i.e.

[σ111, σ222]+ =08, [σ222, σ333]+ =08, [σ333, σ111]+ =08 and the commutators can be written as Kronecker products, i.e.

[σ , σ ]= 2iσ σ σ = 2iσ 111 222 − 3 ⊗ 3 ⊗ 3 − 333 [σ , σ ]= 2iσ σ σ = 2iσ 222 333 − 1 ⊗ 1 ⊗ 1 − 111 [σ , σ ]= 2iσ σ σ = 2iσ . 333 111 − 2 ⊗ 2 ⊗ 2 − 222 Consider now the general case of the three unitary and hermitian matrices

σ11...1, σ22...2, σ33...3 with n Kronecker products. If n is odd the three matrices form a basis of a simple Lie algebra. For n odd the anti-commutators vanish. If n is even the commutators

vanish and the anti-commutators can be expressed as Kronecker products of σ1,

σ2 and σ3. Another useful case is that

[σ σ , σ σ ] = [σ σ , σ σ ]=0 . 1 ⊗ 2 3 ⊗ 3 2 ⊗ 1 3 ⊗ 3 4 This implies that

[σ σ σ σ σ σ σ σ , σ σ σ ]=0 n 0 ⊗ 0 ⊗ 1 ⊗ 0 ⊗··· 0 ⊗ 2 ⊗ 0 ⊗···⊗ 0 3 ⊗ 3 ⊗···⊗ 3 2 with σ and σ at the j’th and k’th position (j = k) with j, k =1, 2,...,n. 1 2 6 From the commutators given above we can also infer that the Hamilton operator

n n J J Hˆ = 12 (σ σ + σ σ )+ 3 (σ σ ) 4 1,j 1,j+1 2,j 2,j+1 4 3,j 3,j+1 j=1 j=1 X X commutes with σ σ σ for both both open end boundary conditions and 3 ⊗ 3 ⊗···⊗ 3 periodic boundary conditions, where

σα,j = σ σ σα σ σ 0 ⊗··· 0 ⊗ ⊗ 0 ⊗···⊗ 0

with σα (α =1, 2, 3) at the j-th position. 3 Pauli Group, Commutator and Anticommuta- tor

The n-qubit Pauli group is defined by

⊗n n := I , σ , σ , σ 1, i P { 2 1 2 3 } ⊗ { ± ± } where σ , σ , σ are the 2 2 Pauli matrices and σ I is the 2 2 identity 1 2 3 × 0 ≡ 2 × matrix. The dimension of the Hilbert space under consideration is dim = 2n. H Thus each element of the Pauli group n is (up to an overall phase 1, i) a Kro- P ± ± necker product of Pauli matrices and 2 2 identity matrices acting on n qubits. × The order of the Pauli group is 22n+2. Thus for n = 1 we have the order 16.

Let j ,...,jn,k ,...,kn 0, 1, 2, 3 and 0 0 ∈{ } n n j0 k0 A =( i) σj , B =( i) σk . − t − t t=1 t=1 O O Using the fact that 3

σjσk = δj,kI2 + i ǫj,k,lσl l X=1 we find, using ǫj ,k ,l = ǫk ,j ,l , t t t − t t t

n 3 j0+k0 AB =( i) δj ,k I + i ǫj ,k ,l σl , − t t 2 t t t t t=1 l ! O Xt=1

n 3 j0+k0 BA =( i) δj ,k I i ǫj ,k ,l σl . − t t 2 − t t t t t=1 l ! O Xt=1 Now, noting that

3 3 1−δj ,k δj ,k I i ǫj ,k ,l σl =( 1) t t δj ,k I + i ǫj ,k ,l σl , t t 2 − t t t t − t t 2 t t t t l l ! Xt=1 Xt=1 we find the following expressions for the commutator and anticommutator of A and B

n n 3 j0+k0 1−δj ,k [A, B]=( i) 1 ( 1) t t δj ,k I + i ǫj ,k ,l σl , − − − t t 2 t t t t t=1 ! t=1 l ! Y O Xt=1 n n 3 j0+k0 1−δj ,k [A, B] =( i) 1+ ( 1) t t δj ,k I + i ǫj ,k ,l σl . + − − t t 2 t t t t t=1 ! t=1 l ! Y O Xt=1 Thus we obtain

n n 1−δ 1−δ [A, B]= 1 ( 1) jt,kt AB, [A, B] = 1+ ( 1) jt,kt AB. − − + − t=1 ! t=1 ! Y Y Since AB = 0, the condition [A, B] = 0 yields 6 n δ n [A, B] = 0 iff ( 1) jt,kt =( 1) . − − t=1 Y

Consequently [A, B] = 0 if and only if the number of coincidences jt = kt is even if n is even, and odd when n is odd. Similarly [A, B]+ = 0 if and only if the number of coincidences jt = kt is even if n is odd, and odd when n is even.

4 Applications

First we look at the eigenvalue problem. Let A, B be two nonzero n n matrices. × Let Av = λv (v = 0) be the eigenvalue equation. If [A, B]=0n, then A(Bv) = 6 λ(Bv). Consequently if Bv = 0, then Bv is an eigenvector of the matrix A. Now 6 let us assume that the anti-commutator vanishes, i.e [A, B]+ =0n. Then we obtain A(Bv)= λ(Bv). Thus if Bv = 0, then Bv is an eigenvector of A corresponding − 6 to the eigenvalue λ. An application is given in section 2 with − A = σ σ σ , B = σ σ σ 1 ⊗ 1 ⊗ 1 3 ⊗ 3 ⊗ 3 with [A, B] = 0 . Thus since +1 is a eigenvalue of A we find that 1 is an + 8 − eigenvalue of A with the eigenvector (σ σ σ )v where Av = v. 3 ⊗ 3 ⊗ 3 One of the main calculations in quantum theory in the Hilbert space Cn is to find eABe−A, where A, B are n n matrices. This is utilized in the solution of the × Heisenberg equation of motion. Now it is well-known that

1 1 eABe−A = B + [A, B]+ [A, [A, B]] + [A, [A, [A, B]]] + 2! 3! ···

If the commutator [A, B] vanishes, we find eABe−A = B. Suppose B is normal with spectral decomposition m

B = λjΠj j=1 X where λj are the m distinct eigenvalues of B and Πj are the projections onto the

corresponding eigenspaces Vj. There exists a unitary n n matrix U such that × m ∗ Λ := UBU = λjIdim(Vj ). j=1 M From eABe−A = B we obtain Λ = UeABe−AU ∗ = (UeAU ∗)Λ(UeAU ∗)−1. Since

the eigenvalues λj are distinct we have

m A ∗ e = U Pj U j=1 ! M

where each Pj is an invertible dim(Vj) dim(Vj) matrix. When × n n j0 k0 A =( i) σj , B =( i) σk − t − t t=1 t=1 O O 2n n−1 (with A, B = I ) we find m = 2 and dim(V1) = dim(V2)=2 so that A = 6 ∗ ∗ V [I n−1 ( I n−1 )]V and B = U[I n−1 ( I n−1 )]U for some unitary U and V . 2 ⊕ − 2 2 ⊕ − 2 Thus eA = V [eI (1/e)I]V ∗ = U ∗[P P ]U. ⊕ 1 ⊕ 2 It follows that eABe−A = B if (UV )[eI (1/e)I](UV )∗ is a direct sum of two ⊕ 2n−1 2n−1 matrices. × There is also a lesser known expansion using the anti-commutator ([4], [5])

1 1 eABeA = B + [A, B] + [A, [A, B] ] + [A, [A, [A, B] ] ] + + 2! + + 3! + + + ···

It follows that

1 1 eABe−A =(B + [A, B] + [A, [A, B] ] + [A, [A, [A, B] ] ] + )e−2A + 2! + + 3! + + + ···

1 1 eABe−A = e2A(B [A, B] + [A, [A, B] ] [A, [A, [A, B] ] ] + ) − + 2! + + − 3! + + + ··· If the anti-commutator of A and B vanishes we obtain

eABe−A = Be−2A and eABe−A = e2AB.

Let A = σ = σ σ σ , B = σ = σ σ σ . 111 1 ⊗ 1 ⊗ 1 222 2 ⊗ 2 ⊗ 2 2 Then [A, B]+ =08 and A = I8 so that

eABe−A = e2AB = (cosh(2)I + sinh(2)σ σ σ )(σ σ σ ). 8 1 ⊗ 1 ⊗ 1 2 ⊗ 2 ⊗ 2 Another application is for spin-Hamilton operators and projection matrices. Let 2 A be an hermitian d d matrix with A = Id. Then × 1 1 Π = (Id + A), Π = (Id A) + 2 − 2 −

are projection matrices which can be used to decompose the Hilbert space Cd into invariant sub Hilbert spaces. Consider for example the spin-Hamilton operators

2

Hˆ = (σj σj I + I σj σj) ⊗ ⊗ 2 2 ⊗ ⊗ j=1 X and 2

Kˆ = (σj σj I + I σj σj + σj I σj). ⊗ ⊗ 2 2 ⊗ ⊗ ⊗ 2 ⊗ j=1 X Then both Hˆ and Kˆ commute with the operator σ σ σ which is an element 3 ⊗ 3 ⊗ 3 of the Pauli group with (σ σ σ )2 = I . Thus we have projection matrices 3 ⊗ 3 ⊗ 3 8 1 1 Π = (I + σ σ σ ), Π = (I σ σ σ ) + 2 8 3 ⊗ 3 ⊗ 3 − 2 8 − 3 ⊗ 3 ⊗ 3 which decomposes the Hilbert space C8 into two four-dimensional sub Hilbert spaces. Then the eigenvalue problem can be solved in these sub Hilbert spaces.

Two orthogonal bases in the Hilbert space Cd

= e ,..., ed , = f ,..., fd A { 1 } B { 1 } are called unbiased if for every 1 j, k d ([6], [7], [8]) ≤ ≤ 1 ej, fk = . |h i| √d

For the Pauli spin matrices σ3, σ1, σ2 the normalized eigenvectors

1 0 = , B3 0 1      1 1 1 1 1 = , B √2 1 √2 1     −  1 1 1 1 2 = , B √2 i √2 i     −  each form an orthonormal basis in C2. Furthermore this is a set of mutually unbiased bases. Consider now σ σ , σ σ , σ σ . Then 3 ⊗ 3 1 ⊗ 1 2 ⊗ 2 1 1 1 0 0 1 0 0 = , , , B3 0 ⊗ 0 0 ⊗ 1 1 ⊗ 0 1 ⊗ 1                  1 1 1 1 1 1 1 1 1 1 1 1 = , , , B1 2 1 ⊗ 1 2 1 ⊗ 1 2 1 ⊗ 1 2 1 ⊗ 1         −   −     −   −  1 1 1 1 1 1 1 1 1 1 1 1 = , , , B2 2 i ⊗ i 2 i ⊗ i 2 i ⊗ i 2 i ⊗ i         −   −     −   −  provide mutually unbiased bases in C4.

† † Another application is with Fermi operators. Let c1, c2, c1, c2 be Fermi creation and annihilation operators with the anticommutation relations

† [cj,ck]+ = δjkI, j,k =1, 2

† † and [cj,ck]+ = [cj,ck]+ = 0 for j, k = 1, 2. Using the Pauli spin matrices we can form the operators

c X =( c† c† ) σ 1 = c†c + c†c 1 1 2 1 c 1 2 2 1  2  † † c1 † † X2 =( c1 c2 ) σ2 = ic1c2 + ic2c1 c2 −   † † c1 † † X3 =( c1 c2 ) σ3 = c1c1 c2c2. c2 −   Then we find the well-known result for the commutators

[X1,X2]=2iX3, [X2,X3]=2iX1, [X3,X1]=2iX2.

For the anti-commutators we find

[X1,X2]+ =0, [X2,X3]+ =0, [X3,X1]+ =0. This can be extended to higher dimensions. Considering Bose creation and anni- † † hilation operators b1, b2, b1, b2 and

b Y =( b† b† ) σ 1 = b† b + b† b 1 1 2 1 b 1 2 2 1  2  † † b1 † † Y2 =( b1 b2 ) σ2 = ib1b2 + ib2b1 b2 −   † † b1 † † Y3 =( b1 b2 ) σ3 = b1b1 b2b2 b2 −  

† with the commutation relations [bj, bk]= δjkI (j, k =1, 2) provides the well-known result [Y1,Y2]=2iY3, [Y2,Y3]=2iY1, [Y3,Y1]=2iY2 for the commutators.

5 Conclusion

We studied the Pauli group and found the conditions when the commutator or anticommutator of two elements of the Pauli group vanish. Six applications for these conditions have been discussed. Acknowledgment The authors are supported by the National Research Foundation (NRF), South Africa. This work is based upon research supported by the National Research Foundation. Any opinion, findings and conclusions or recommendations expressed in this material are those of the author(s) and therefore the NRF do not accept any liability in regard thereto.

References

[1] W.-H. Steeb and Y. Hardy, Matrix Calculus and Kronecker Product: A Practi- cal Approach to Linear and Multilinear Algebra, 2nd edition, World Scientific, Singapore (2011)

[2] M. A. Nielsen and I. L. Chuang, Quantum Computing and Quantum Infor- mation, Cambridge University Press, Cambridge (2000)

[3] W.-H. Steeb and Y. Hardy, Quantum Mechanics using Computer Algebra, 2nd edition, World Scientific, Singapore (2010) [4] I. Mendaˇs and P. Milutinovi´c, “Anticommutator analogue of the Baker- Hausdorff lemma”, J. Phys. A: Math. Gen. 22, L687-L689 (1989)

[5] W.-H. Steeb, I. Tanski and Y. Hardy, Groups, Lie Groups and Lie Algebras with Applications, World Scientific, Singapore (2012)

[6] M. Matolcsi, “A Fourier analytic approach to the problem of mutually unbi- ased bases”, arXiv:1009.2407v1

[7] M. Revzen, “Geometrical view of the Mean King Problem”, arXiv:1205.5406v1

[8] A. Kalev, “A geometrical relation between symmetric operators and mutually unbiased operators”, arXiv:1305.6044v1