arxiv: v1 [quant-ph] 6 Jun 2008

Transcription

1 Bloch vectors for quits Reinhol A. Bertlmann an Philipp Krammer Faculty of Physics, University of Vienna, Boltzmanngasse 5, A-1090 Vienna, Austria arxiv: v1 [quant-ph] 6 Jun 008 We present three ifferent matrix bases that can be use to ecompose ensity matrices of imensional quantum systems, so-calle quits: the generalize Gell- Mann matrix basis, the polarization operator basis, an the Weyl operator basis. Such a ecomposition can be ientifie with a vector the Bloch vector, i.e. a generalization of the well known qubit case an is a convenient expression for comparison with measurable quantities an for explicit calculations avoiing the hanling of large matrices. We present a new metho to ecompose ensity matrices via so calle stanar matrices, consier the important case of an isotropic two quit state an ecompose it accoring to each basis. In case of qutrits we show a representation of an entanglement witness in terms of expectation values of spin 1 measurements, which is appropriate for an experimental realization. I. INTRODUCTION The state of a imensional quantum system a quit is usually escribe by a ensity matrix. For high imensions, where the matrices become large (for composite systems of n particles the matrices are of even much larger imension n n, a simple way to express ensity matrices is of great interest. Since the space of matrices is a vector space, there exist bases of matrices which can be use to ecompose any matrix. For qubits such a basis contains the three Pauli matrices, accoringly, a ensity matrix can be expresse by a 3 imensional vector, the Bloch vector, an any such vector has to lie within the so-calle Bloch ball [1, ]. Unique for qubits is the fact that any point on the sphere, Bloch sphere, an insie the ball correspons to a physical state, i.e. a ensity matrix. The pure states lie on the sphere an the mixe ones insie. In higher imensions there exist ifferent matrix bases that can be use to express quits as ( 1 imensional vectors as well. Different to the qubit case, however, is that the map inuce is not bijective: not every point on the Bloch sphere in imensions 1 correspons to a physical state. Nevertheless the vectors are often also calle Bloch vectors (see in this context, e.g., Refs. [3, 4, 5, 6, 7]. In this paper we want to present an compare three ifferent matrix bases for a Bloch vector ecomposition of quits. In Sec. II we propose the properties of any matrix basis for using it as a practical ecomposition of ensity matrices an recall the general notation of Bloch vectors. In Secs. III V we offer three ifferent matrix bases: the generalize Gell-Mann matrix basis, the polarization operator basis, an the Weyl operator basis. For all these bases we give examples in the imensions of our interest an present the ifferent Electronic aress: reinhol.bertlmann@univie.ac.at Electronic aress: philipp.krammer@univie.ac.at

2 Bloch vector ecompositions of an arbitrary ensity matrix in the stanar matrix notation. Next in Sec. VI, by constructing tensor proucts of states we stuy the isotropic two quit state an present the results for the three matrix ecompositions, i.e. for the three ifferent Bloch vectors. In Sec. VII we focus on the isotropic two quit state an calculate the Hilbert Schmit measure of entanglement (see, e.g., Refs. [8, 9, 10, 11]. Its connection to the optimal entanglement witness is shown, which is etermine in terms of the three matrix bases. An example for the experimental realization of an entanglement witness is given in Sec. VII B. The mathematical an physical avantages/isavantages by using the three ifferent matrix bases are iscusse in Sec. VIII, where also the final conclusions are rawn. II. PRELIMINARIES A quit state is represente by a ensity operator in the Hilbert Schmit space acting on the imensional Hilbert space H that can be written as a matrix the ensity matrix in the stanar basis { k }, with k = 1,,... or k = 0, 1,, Properties of a practical matrix basis. For practical reasons the general properties of a matrix basis which is use for the Bloch vector ecomposition of quits are the following: i The basis inclues the ientity matrix ½ an 1 matrices {A i } of imension which are traceless, i.e. TrA i = 0. ii The matrices of any basis {A i } are orthogonal, i.e. Tr A i A j = N δ ij with N Ê. (1 Bloch vector expansion of a ensity matrix. Since any matrix in the Hilbert-Schmit space of imension can be ecompose with a matrix basis {A i }, we can of course ecompose a quit ensity matrix as well an get the Bloch vector expansion of the ensity matrix, ρ = 1 ½ b Γ, ( where b Γ is a linear combination of all matrices {A i } an the vector b Ê 1 with b i = Γ i = TrρΓ i is calle Bloch vector. The term 1 ½ is fixe because of conition Trρ = 1. Remark. Note that a given ensity matrix ρ can always be ecompose into a Bloch vector, but not any vector σ that is of the form ( is automatically a ensity matrix, even if it satisfies the conitions Trσ = 1 an Trσ 1 since generally it oes not imply σ 0. Each ifferent matrix basis inuces a ifferent Bloch vector lying within a Bloch hypersphere where, however, not every point of the hypersphere correspons to a physical state (with ρ 0; these points are exclue (holes. The geometric character of the Bloch space in higher imensions turns out to be quite complicate an is still of great interest (see Refs. [3, 4, 5, 6, 7]. All ifferent Bloch hyperballs are isomorphic since they correspon to the same ensity matrix ρ. The interesting question is which Bloch hyperball which matrix basis is optimal for a specific purpose, like the calculation of the entanglement egree or the etermination of the geometry of the Hilbert space or the comparison with measurable quantities.

3 3 III. THE GENERALIZED GELL-MANN MATRIX BASIS A. Definition an example The generalize Gell-Mann matrices (GGM are higher imensional extensions of the Pauli matrices (for qubits an the Gell-Mann matrices (for qutrits, they are the stanar SU(N generators (in our case N =. They are efine as three ifferent types of matrices an for simplicity we use here the operator notation; then the ensity matrices follow by simply writing the operators in the stanar basis (see, e.g. Refs. [3, 1]: i ( 1 symmetric GGM ii ( 1 antisymmetric GGM Λ jk s = j k k j, 1 j < k, (3 Λ jk a = i j k i k j, 1 j < k, (4 iii ( 1 iagonal GGM ( l Λ l = j j l l 1 l 1, 1 l 1. (5 l(l 1 j=1 In total we have 1 GGM; it follows from the efinitions that all GGM are Hermitian an traceless. They are orthogonal an form a basis, the generalize Gell-Mann matrix basis (GGB. A proof for the orthogonality of GGB we present in the Appenix A1. Examples. Let us recall the case of imension 3, the 8 Gell-Mann matrices (for a representation see, e.g., Refs. [11, 13] i 3 symmetric Gell-Mann matrices λ 1 s = , λ 13 s = , λ 3 s = , (6 ii 3 antisymmetric Gell-Mann matrices λ 1 a = 0 i 0 i 0 0, λ 13 a = 0 0 i i 0 0, λ 3 a = i 0 i 0, (7 ii iagonal Gell-Mann matrices λ 1 = , λ = (8 0 0

4 4 To see how they generalize for higher imensions we show the case we nee for quits of imension = 4 : i 6 symmetric GGM Λ 1 s = Λ 3 s = ii 6 antisymmetric GGM Λ 1 a = Λ 3 a = 0 i 0 0 i i 0 0 i 0 0,,,, Λ13 s = Λ4 s = Λ13 a = Λ4 a = i 0 i i 0 i 0 0,,,, Λ14 s = Λ34 s = Λ14 a = Λ34 a = , i i i 0 0 i 0, (9,, (10 iii 3 iagonal GGM Λ 1 = , Λ = , Λ3 = (11 Using the GGB we obtain, in general, the following Bloch vector expansion of a ensity matrix: ρ = 1 ½ b Λ, (1 with the Bloch vector b = ( {b jk s }, {bjk a }, {bl }, where the components are orere an for the inices we have the restrictions 1 j < k an 1 l 1. The components are given by b jk s = TrΛ jk s ρ, b jk a = TrΛ jk a ρ an b l = TrΛ l ρ. All Bloch vectors lie within a hypersphere of raius b ( 1/. For example, for qutrits the Bloch vector components are b = ( b 1 s, b 13 s, b 3 s, b 1 a, b 13 a, b 3 a, b 1, b corresponing to the Gell-Mann matrices (6, (7, (8 an b 1/3. As alreay mentione the allowe range of b is restricte. It has an interesting geometric structure which has been calculate analytically for the case of qutrits by stuying imensional planes in the 8 imensional Bloch space [3] or numerically by consiering 3 imensional cross sections [7]. In any case, pure states lie on the surface an the mixe ones insie.

5 5 B. Stanar matrix basis expansion by GGB The stanar matrices are simply the matrices that have only one entry 1 an the other entries 0 an form an orthonormal basis of the Hilbert Schmit space. We write these matrices shortly as operators j k, with j, k = 1,...,. (13 Any matrix can easily be ecompose into a vector via a certain linear combination of the matrices (13. Knowing the expansion of matrices (13 into GGB we can therefore fin the ecomposition of any matrix in terms of the GGB. We fin the following expansion of stanar matrices (13 into GGB : j k = ( 1 Λ jk ( s 1 Λ kj s iλ jk a iλ kj a j 1 Λ j 1 j 1 j n=0 for j < k for j > k (14 1 Λ jn 1 (jn(jn1 ½ for j = k. Proof. The first two cases can be easily verifie. To show the last case we first set up a recurrence relation for l l, which we obtain by eliminating the term l 1 j=1 j j in the two expressions (5 for Λl an Λ l 1 l 1 l 1 l l = Λ l 1 Λ l l 1 l 1, (15 l l an we consier the case l 1 = 1 1 = ( 1 Λ ( 1 Λ 1. (16 From Λ 1 given by Eq. (5 Λ 1 = ( 1 ( 1 j j ( 1, (17 j=1 we get the Bloch vector ecomposition of ( = 1 ( 1 Λ 1 ½, (18 where we have applie 1 j=1 j j = ½. Inserting now ecomposition (18 into relation (16 we gain the Bloch vector expansion for 1 1 an recurrence relation (15 provies an so forth. Thus finally we fin n n = n 1 ( n Λ n 1 n 1 k=0 the relation we ha to prove, where n = j. 1 ( n k 1( n k Λ nk 1 ½, (19

6 6 IV. THE POLARIZATION OPERATOR BASIS A. Definition an examples The polarization operators in the Hilbert-Schmit space of imension are efine as the following matrices [4, 14] : T LM = L 1 s 1 k,l=1 C sm k sm l, LM k l. (0 The use inices have the properties s = 1, L = 0, 1,..., s, M = L, L 1,..., L 1, L, m 1 = s, m = s 1,...,m = s. (1 The coefficients C sm k sm l, LM are ientifie with the usual Clebsch Goran coefficients Cjm j 1 m 1, j m of the angular momentum theory an are isplaye explicitly in tables, e.g., in Ref. [14]. For L = M = 0 the polarization operator is proportional to the ientity matrix [4, 14], T 00 = 1 ½. ( It is shown in Ref. [4] that all polarization operators (except T 00 are traceless, in general not Hermitian, an that orthogonality relation (1 is satisfie Tr T L 1 M 1 T L M = δ L1 L δ M1 M. (3 Therefore the polarization operators (0 form an orthonormal matrix basis the polarization operator basis (POB of the Hilbert Schmit space of imension. Examples. The simplest example is of imension, the qubit. For a qubit the POB is given by the following matrices (s = 1/; L = 0, 1; M = 1, 0, 1 ( ( T 00 = , T =, 0 0 ( ( T 10 = , T =. (4 1 0 For the next higher imension = 3 (s = 1, the case of qutrits, we get 9 polarization

7 7 operators T LM with L = 0, 1, an M = L,..., L an we have T 11 = , T 10 = , T 1 1 = , T = , T 1 = , T 0 = , T 1 = , T = ( Then the ecomposition of any ensity matrix into a Bloch vector by using the POB has, in general, the following form: ρ = 1 s ½ L L=1 M= L b LM T LM = 1 ½ b T, (6 with the Bloch vector b = (b 1 1, b 10, b 11, b, b 1, b 0,..., b LM, where the components are orere an given by b LM = Tr T LM ρ. In general the components b LM are complex since the polarization operators T LM are not Hermitian. All Bloch vectors lie within a hypersphere of raius b ( 1/. In imensions the Bloch vector b = (b 1 1, b 10, b 11 is limite by b 1 an forms a spheroi [4], the pure states occupy the surface an the mixe ones lie in the volume. This ecomposition is fully equivalent to the stanar escription of Bloch vectors with Pauli matrices. In higher imensions, however, the structure of the allowe range of b (ue to the positivity requirement ρ 0 is quite complicate, as can be seen alreay for = 3 (for etails see Ref. [4]. Nevertheless, pure states are on the surface, mixe ones lie within the volume an the maximal mixe one correspons to b = 0, thus b is a kin of measure for the mixeness of a quantum state. B. Stanar matrix basis expansion by POB The stanar matrices (13 can be expane by the POB as [14] i j = L M L 1 s 1 Csm i sm j, LM T LM. (7 Note that M is actually fixe by the conition m j M = m i.

8 8 Proof. Inserting efinition (0 on the right han sie (RHS of equation (7 we fin RHS = ( L 1 s 1 Csm i sm j, LM Csm k sm l, LM k l = k,l L = δ jl δ ik k l k,l = i j, (8 where we use the sum rule for Clebsch Goran coefficients [14] c,γ c 1 b 1 Cbβ aα, cγ Cbβ aα, cγ = δ αα δ ββ. (9 V. WEYL OPERATOR BASIS A. Definition an example Finally we want to iscuss a basis of the Hilbert Schmit space of imension that consists of the following operators: U nm = 1 k=0 e πi kn k (k m mo n, m = 0, 1,..., 1, (30 where we use the stanar basis of the Hilbert space. The operators in notation (30 have been introuce in the context of quantum teleportation of quit states [15] an are often calle Weyl operators in the literature (see e.g. Refs. [16, 17, 18]. The operators (30 are unitary an form an orthonormal basis of the Hilbert Schmit space Tr U nmu lj = δ nl δ mj (31 (a proof is presente in Appenix A3 the Weyl operator basis (WOB. They can be use to create a basis of maximally entangle quit states [16, 19, 0]. Clearly the operator U 00 represents the ientity U 00 = ½. Example. Let us show the example of imension 3, the qutrit case. There the Weyl operators (30 have the following matrix form U 01 = , U 0 = , ( U 10 = 0 e πi/3 0, U 11 = 0 0 e πi/3, U 1 = e πi/3 0 0, 0 0 e πi/3 e πi/ e πi/ U 0 = 0 e πi/3 0, U 1 = 0 0 e πi/3, U = e πi/ e πi/3 e πi/ e πi/3 0

9 9 Using the WOB we can ecompose quite generally any ensity matrix into a Bloch vector ρ = 1 1 ½ n,m=0 b nm U nm = 1 ½ b U, (33 with n, m = 0, 1,..., 1 (b 00 = 0. The components of the Bloch vector b = ( {b nm } are orere an given by b nm = Tr U nm ρ. In general the components b nm are complex since the Weyl operators are not Hermitian an the complex conjugates fulfil the relation b n m = e πi nm b n m, which follows easily from efinition (30 together with the hermiticity of ρ. All Bloch vectors lie within a hypersphere of raius b 1/. For example, for qutrits the Bloch vector is expresse by b = (b 01, b 0, b 10, b 11, b 1, b 0, b 1, b an b /3. In 3 an higher imensions the allowe range of the Bloch vector is quite restricte within the hypersphere an the etaile structure is not known yet. Note that in imensions the WOB as well as the GGB coincies with the Pauli matrix basis an the POB represents a rotate Pauli basis (where σ ± = 1 (σ 1 ± iσ, in particular {U 00, U 01, U 10, U 11 } = {½, σ 1, σ 3, iσ }, (34 { ½, λ 1 s, λ 1 a, λ1} = {½, σ 1, σ, σ 3 }, (35 { } 1 {T 00, T 11, T 10, T 1 1 } = ½, σ 1, σ 3, σ. (36 B. Stanar matrix basis expansion by WOB The stanar matrices (13 can be expresse by the WOB in the following way j k = 1 1 l=0 e πi lj U l (k jmo. (37 Proof. We insert the efinition of the Weyl operators (30 on the right han sie (RHS of Eq. (37, use Eq. (A4 an get RHS = 1 1 l,r=0 e πi = j k 1 l(r j r (r k j mo 1 1 r j, r=0 l=0 e πi l(r j r (r k j mo = j k. (38 VI. ISOTROPIC TWO QUDIT STATE Now we consier bipartite systems in a imensional Hilbert space H A H B. The observables acting in the subsystems H A an H B are usually calle Alice an Bob in quantum communication.

10 10 can be ecom- Quite generally, a ensity matrix of a two quit state acting on HA H B pose in the following way (neglecting the reference to A an B ρ = 1 ½ ½ n i Γ i ½ m i ½ Γ i c ij Γ i Γ j, n i, m i, c ij, (39 where {Γ i } represents some basis in the subspace H. The term c ij Γ i Γ j always can be iagonalize by two inepenent orthogonal transformations on Γ i an Γ j [1]. Altogether there are ( 1 terms. However, for isotropic two quit states the case we consier in our paper the secon an thir term in expression (39 vanish an the fourth term reuces to c ii Γ i Γ i, which implies the vanishing of ( 1 ( 1 = ( 1 terms. Consequently, for an isotropic two quit ensity matrix there remain 1 inepenent terms, which provies the imension of the corresponing Bloch vector. Thus the isotropic two quit Bloch vector is of the same imension lives in the same subspace as the one quit vector, which is a comfortable simplification. Explicitly, the isotropic two quit state ρ α ( is efine as follows [, 3, 4] : ρ ( α = α φ φ 1 α ½, α Ê, 1 1 α 1, (40 where the range of α is etermine by the positivity of the state. The state φ, a Bell state, is maximally entangle an given by φ 1 = j j j, (41 where { j } enotes the stanar basis of the imensional Hilbert space. A. Expansion into GGB Let us first calculate the Bloch vector notation for the Bell state φ φ in the GGB. It is convenient to split the state into two parts where A an B are efine by φ φ = 1 j k j k j,k=1 = A B, (4 A := 1 j k j k 1 k j k j, (43 j<k j<k B := 1 j j j j, (44 an to calculate the two terms separately. j

11 11 For term A we use the stanar matrix expansion (14 for the case j k an get [ A = 1 ( ( ( ( ] Λ jk 4 s iλ jk a Λ jk s iλ jk a Λ jk s iλ jk a Λ jk s iλ jk a = 1 j<k i<j ( Λ jk s Λ jk s j<k Λ jk a Λjk a. (45 For term B we nee the case j = k in expansion (14 an obtain after some calculations (the etails are presente in Appenix A B = 1 1 m=1 Λ m Λ m 1 ½ ½. (46 Thus all together we fin the following GGB Bloch vector notations, for the Bell state (4 φ φ = 1 ½ ½ 1 Λ, (47 an for the isotropic two quit state (40 ρ ( α = 1 ½ ½ α Λ, (48 where we efine Λ := i<j Λ jk s Λjk s i<j Λ jk a Λjk a 1 m=1 Λ m Λ m. (49 B. Expansion into POB Now we calculate the Bell state φ φ in the POB. Using expansion (7 an the sum rule for the Clebsch Goran coefficients [14] C cγ aα,bβ Ccγ aα,b β = c 1 b 1 δ bb δ ββ, (50 we obtain φ φ = 1 α,γ i j i j i,j=1 = 1 ( (L 1(L 1 s 1 L,L i,j = 1 (L 1(L 1 δ L,L T LM T L L 1 M L,L = 1 T LM T LM L C sm i sm j,lm Csm i sm j,l M T LM T L M = 1 ½ ½ 1 T, (51

12 1 where we extracte the unity (recall Eq. ( an efine T := T LM T LM. (5 L,M 0,0 Result (51 provies the POB Bloch vector notation of the isotropic two quit state (40 ρ ( α = 1 ½ ½ α T. (53 C. Expansion into WOB Finally we present the Bell state in the WOB (the etails for our approach using the stanar matrix expression (37 can be foun in the Appenix A4, see also Ref. [16] φ φ = 1 ½ ½ 1 U, (54 with U := 1 l,m=0 U lm U lm, (l, m (0, 0, (55 where negative values of the inex l have to be consiere as mo, an from formula (54 we fin the WOB Bloch vector notation of the isotropic two quit state ρ ( α = 1 ½ ½ α U. (56 VII. APPLICATIONS OF THE MATRIX BASES A. Entangle isotropic two quit states In Ref. [11] the connection between the Hilbert Schmit (HS measure of entanglement [8, 9, 10] an the optimal entanglement witness is investigate. Explicit calculations for both quantities are presente in case of isotropic qutrit states. For higher imensions, the isotropic two quit states, the above quantities are etermine as well but in terms of a rather general matrix basis ecomposition. With the results of the present paper we can calculate all quantities explicitly. Let us recall the basic notations we nee. The HS measure is efine as the minimal HS istance of an entangle state ρ ent to the set of separable states S D(ρ ent := min ρ S ρ ρ ent = ρ 0 ρ ent, (57 where ρ 0 enotes the nearest separable state, the minimum of the HS istance. An entanglement witness A A (A = A A A B, the HS space of operators acting on the Hilbert space of states is a Hermitian operator that etects the entanglement of a state ρ ent via inequalities [10, 5, 6, 7].

13 13 FIG. 1: Illustration of the Bertlmann Narnhofer Thirring Theorem (6 Definition 1. An entanglement witness A is a Hermitian operator with the following properties: The expectation value of A is negative for an entangle state, whereas it is non negative for any separable state. ρ ent, A = Tr ρ ent A < 0, ρ, A = Tr ρa 0 ρ S. (58 The fact, however that there exists an operator satisfying inequalities (58 for any entangle state, i.e. that the efinition is meaningful, has to be prove; it follows from the Hahn Banach Theorem of functional analysis (for a simple geometric approach, see Ref. [11]. An entanglement witness is optimal, enote by A opt, if apart from Eq. (58 there exists a separable state ρ 0 S such that ρ 0, A opt = 0. (59 The operator A opt efines a tangent plane to the set of separable states S an all states ρ p with ρ p, A opt = 0 lie within that plane; see Fig. 1. Let us call the lower one of the inequalities (58 an entanglement witness inequality, short EWI. It etects entanglement whereas a Bell inequality etermines non locality. Rewriting Eq. (58 as ρ, A ρ ent, A 0 ρ S, (60 the maximal violation of the EWI is efine by B(ρ ent = max A, A a½ 1 ( min ρ S ρ, A ρ ent, A, (61 where the maximum is taken over all possible entanglement witnesses A, suitably normalize. Then an interesting connection between the HS measure an the concept of entanglement witnesses is given by the Bertlmann Narnhofer Thirring Theorem, illustrate in Fig. 1 [10].

14 14 Theorem 1. i The maximal violation of the EWI is equal to the minimal istance of ρ ent to the set S B(ρ ent = D(ρ ent. (6 ii The maximal violation of the EWI is attaine for an optimal entanglement witness A opt = ρ 0 ρ ent ρ 0, ρ 0 ρ ent ½ ρ 0 ρ ent. (63 Thus the calculation of the optimal entanglement witness A opt to a given entangle state ρ ent reuces to the etermination of the nearest separable state ρ 0. In special cases ρ 0 is etectable but in general its etection is quite a ifficult task. We are able to fin the nearest separable state by working with Lemma 1, a metho we call guess metho [11]. Lemma 1. A state ρ is equal to the nearest separable state ρ 0 if an only if the operator C = ρ ρ ent ρ, ρ ρ ent ½ ρ ρ ent (64 is an entanglement witness. Lemma 1 probes if a guess ρ is inee correct for the nearest separable state. Then operator C represents the optimal entanglement witness A opt (63. Now let us apply the matrix bases we iscusse in the previous sections an calculate the quantities introuce above. As an entangle state we consier the isotropic two quit (, ent state ρ α, that is the state ρ α ( 1 (40 for < α 1. 1 Starting with the GGB we can express that state in our Bloch vector notation by formula (48. By using Lemma 1 we fin that the nearest separable state is reache at α = 1 1 It provies the HS measure D(ρ ( α,ent = ρ ( 0 = ρ ( α= 1 1 ρ ( 0 ρ ( = 1 ½ ½ 1 ( 1 Λ. (65 α, ent an the optimal entanglement witness (63 A opt (ρ ( α, ent = 1 = 1 ( α 1, ( ½ ½ 1 1 Λ, (67 where we use the HS norm Λ = 1. Clearly, the maximal violation B of the EWI equals the HS measure D B(ρ ( α, ent = ρ ( α, ent, A opt ( 1 = α 1 = D(ρ ( α, 1 ent. (68

15 15 For expressing above quantities by the matrix bases POB an WOB it suffices to calculate the proportionality factors between Λ, T an U. By comparison of the three forms for the isotropic quit state (48, (53 an (56 we fin It provies the following expressions, for the POB an for the WOB Λ = T an T = 1 U. (69 ρ ( 0 = ρ ( = 1 α= 1 1 ½ ½ 1 ( 1 T, (70 A opt (ρ ( α, ent = ½ ½ 1 1 T, (71 ρ ( 0 = ρ ( α= 1 1 A opt (ρ ( α, ent = 1 = 1 ½ ½ 1 ( 1 U, (7 1 1 ½ ½ 1 1 U. (73 Of course, the HS measure D(ρ ( α, ent remains the same expression (66 inepenent of the chosen matrix basis, which can easily be verifie using T = 1 an U = 1. B. Entanglement witness representation for experiments Entanglement witnesses are Hermitian operators an therefore observables that shoul be measurable in a given experimental set up an thus provie an experimental verification of entanglement. The quantity to be measure is the expectation value A = TrAρ (74 of an entanglement witness A for some state ρ. If A < 0 then the state ρ is entangle. But which measurements have to be performe? Obviously it is appropriate to express the entanglement witness in terms of generalize Gell-Mann matrices (3 (5, since they are Hermitian. For = 3 qutrits the Gell-Mann matrices (6 (8 can be expresse in terms of eight physical operators, the observables S x, S y, S z, S x, S y, {S x, S y }, {S y, S z }, {S z, S x } of a spin 1 system, where S = (S x, S y, S z is the spin operator an {S i, S j } = S i S j S j S i (with i, j = x, y, z enotes the corresponing anticommutator. The ecomposition of the Gell-Mann matrices into spin 1 operators is as follows (for a similar expansion, see Ref. [7]: λ 1 s = 1 h ( hs x {S z, S x }, λ 13 s = 1 h ( S x S y, λ 3 s = 1 h ( hs x {S z, S x }, λ 1 a = 1 h ( hs y {S y, S z }, λ 13 a = 1 h {S x, S y }, λ 3 a = 1 h ( hs y {S y, S z }, λ 1 = ½ 1 ( hsz h 3Sx 3S y, λ = 1 ( ½ 3 ( hsz 3 h Sx S y. (75

16 16 All operators can be represente by the following matrices: S x = h , S y = h 0 i 0 i 0 i, S z = h , i Sx = h , Sy = h , {S x, S y } = h 0 0 i 0 0 0, {S y, S z } = h 0 i 0 i 0 i, i i 0 {S z, S x } = h ( Thus we can express any observable on a n qutrit Hilbert space a composite system of n particles with 3 egrees of freeom in terms of above spin operators (76. As an example we want to stuy the entanglement witness for the isotropic two qutrit state, i.e. state (40 for = 3. In this case we obtain for the optimal entanglement witness A iso = 1 3 (½ ½ 34 Λ, (77 (i.e. Eq. (67 for = 3 where the operator Λ is efine in Eq. (49. Expressing the Gell-Mann matrices in Λ (49 by the spin operator ecomposition (75 we fin for the expectation value of the entanglement witness A iso A iso = ½ ½ 4 where Λ = 1 h ( S x S x S y S y S z S z h ( ½ S x ½ S y S x ½ S y ½ Λ, (78 ½ ½ 4 h ( S 4 x Sx Sy Sy h ( S 4 x Sy Sy Sx 1 h ( {S 4 z, S x } {S z, S x } {S y, S z } {S y, S z } {S x, S y } {S x, S y }. (79 Decomposition (79 has to be etermine experimentally by measuring the several expectation values with the set ups on both Alice s an Bob s sie. The avantage of the entanglement witness proceure is that for an experimental outcome A iso < 0 the consiere quantum state is efinitely entangle, whereas in case of Bell inequalities a violation etects nonlocal states. That means by the entanglement witness proceure we are able to etect more entangle states than with Bell inequalities. The amount of measurement steps necessary to etermine an entanglement witness is about the same as in the Bell inequality proceure (see, e.g., Refs. [8, 9, 30, 31].

17 17 VIII. CONCLUSION In this paper we present three ifferent matrix bases which are quite useful to ecompose ensity matrices for higher imensional quits. These are the generalize Gell-Mann matrix basis, the polarization operator basis an the Weyl operator basis. Each ecomposition we ientify with a vector, the so-calle Bloch vector. Consiering just one particle states we observe the following features: The generalize Gell-Mann matrix basis is easy to construct, the matrices correspon to the stanar SU(N generators (N =, but in general (in imensions it is rather unpractical to work with the iagonal matrices (5 ue to their more complicate efinition. On the other han, the Bloch vector itself has real components, which is avantageous, they can be expresse as expectation values of measurable quantities. For example, in 3 imensions the Gell-Mann matrices are Hermitian an the Bloch vector components can be expresse by expectation values of spin 1 operators. The polarization operator basis is also easy to set up, all you nee to know are the Clebsch Goran coefficients which you fin tabulate in the literature. However, the Bloch vector contains complex components. For the Weyl operator basis the corresponing operators are again simple to construct, they are non Hermitian but unitary. The Bloch vector itself has a very simple structure, however, with complex components. Let us note that in imensions all bases are equivalent since they correspon to Pauli matrices or linear combinations thereof. In case of two quits we have stuie the isotropic states explicitly an fin the following: In the generalize Gell-Mann matrix basis the Bloch vector (48 with expression (49 is more complicate to construct, in particular the iagonal part B (46 (see Appenix A. In the polarization operator basis the Bloch vector (53 with expression (5 can be easily set up by the knowlege of the Clebsch Goron coefficient sum rule (50 an in the Weyl operator basis the Bloch vector (56 with efinition (55 is actually most easily to construct. The Hilbert Schmit measure of entanglement can be calculate explicitly for all isotropic two quit states an we want to emphasize its interesting connection to the maximal violation of the entanglement witness inequality, Theorem 1. For the experimental realization of an entanglement witness the generalize Gell-Mann matrix basis is the appropriate one since the generalize Gell-Mann matrices are Hermitian. For a ifferent task, however, the etermination of the geometry of entanglement the Weyl operator basis turns out to be optimal. In our example of the entangle isotropic two qutrit state the entanglement witness can be expresse by experimental quantities, the expectation values of spin 1 measurements. In this way one can experimentally fin out whether a state is entangle or not, i.e., we can obtain rather precise information on the quality of entanglement. Quite generally, the Bloch vector ecomposition into one of the three matrix bases is of particular avantage in the construction of entanglement witnesses. It turns out that if the coefficients of the ecomposition satisfy a certain conition the consiere operator represents an entanglement witness, i.e. satisfies inequalities (58 (for etails see Ref. [3]. Acknowlegments We woul like to thank Beatrix Hiesmayr an Heie Narnhofer for helpful iscussions. This research has been financially supporte by FWF project CoQuS No W110-N16 of the Austrian Science Founation an by F140-N Research Grant of the University of Vienna.

18 18 APPENDIX A 1. Proof of Orthogonality of GGB We want to proof conition (1 for the GGB which consists of the 1 GGM (3, (4, (5 an the unity ½. Since all GGM are Hermitian (thus TrA i A j = TrA i A j = TrA j A i it suffices to proof the following conitions: Proof of conition (A1. Inserting efinition (3 we have Tr Λ jk s Λmn s = Tr Λ jk s Λ mn s = δ jm δ kn (A1 Tr Λ jk a Λ mn a = δ jm δ kn (A Tr Λ l Λ m = δ lm (A3 Tr Λ jk a Λmn s = 0 (A4 Tr Λ jk s Λm = 0 (A5 Tr Λ jk a Λm = 0. (A6 l ( j k k j ( m n n m l l=1 = l ( l j k m n l l j k n m l l k j m n l l k j n m l = δ jn δ km δ jm δ kn δ kn δ jm δ km δ jn = δ jm δ kn, (A7 where we use in the last step that δ jn δ km = 0 since we have j < k an m < n. Proof of conition (A. This case is equivalent to the one before apart from change signs that o not matter Tr Λ jk a Λmn a = δ jn δ km δ jm δ kn δ kn δ jm δ km δ jn = δ jm δ kn. (A8 Proof of conition (A3. Using efinition (5 an enoting C l = l(l 1, (A9 where l m without loss of generality, we get ( l m Tr Λ l Λ m = C l C m p k k n n p lm p l 1 l 1 m 1 m 1 p m p=1 k=1 n=1 l p k k m 1 m 1 p l k=1 = C l C m (l lm δ lm m m p l 1 l 1 n n p n=1 l δ k(m1 l k=1 n=1 m δ n(l1. (A10

19 19 Using the fact that δ k(m1 = 0 for m k an { m l δ n(l1 0 if l = m = l if l < m we obtain n=1 Tr Λ l Λ m = (C l l(l 1 δ lm = δ lm. Proof of conition (A4. Analogously to the proofs (A7 an (A8 we fin (A11 (A1 Tr Λ jk a Λ mn s = i ( δ jn δ km δ jm δ kn δ jm δ kn δ jn δ km = 0. (A13 Proof of conition (A5. Inserting efinitions (3 an (5 gives Tr Λ jk s Λ m ( = C m m p k j m 1 m 1 p m p j k m 1 m 1 p p=1 m p j k n n p n=1 = m δ j(m1 δ k(m1 = 0, m p k j n n p n=1 m δ kl δ jl l=1 (A14 since per efinition we have j < k. Proof of conition (A6. This proof is equivalent to the previous one since constant factors in front of the terms o not matter.. Calculation of term B in GGB To obtain the Bloch vector notation of term B (44 we insert the stanar matrix expansion (14 for the case j = k. We split the tensor proucts in the following way B = 1 (B 1 B B 3 B 4 1 ½ ½, (A15

20 0 where the terms B 1,...,B 4 are introuce by (note that Λ 0 = 0 B 1 = j 1 j 1 j Λj 1 Λ j 1 1 (j n(j n 1 Λjn Λ jn B = B 3 B 4 = 1 = 1 j=1 ( j=1 m=1 j 1 l=0 j 1 n=0 j 1 n l, n,l=0 ( j=1 ( j=1 n(=l=0 j 1 4j(j l(j l 1 Λj 1 Λ jl j 1 4j(j n(j n 1 Λjn Λ j 1 1 (j n(j n 1(j l(j l 1 Λjn Λ jl j 1 j Λj 1 ½ j 1 j ½ Λj 1 j 1 n=0 j 1 l=0 Only the first term B 1 (A16 gives a contribution 1 ( m B 1 = (m 1 m m(m 1 1 Λ jn ½ (j n(j n 1 (A16 (A17 (A18 1 ½ Λ jl. (A19 (j l(j l 1 Λ m Λ m = 1 1 m=1 Λ m Λ m, (A0 whereas the remaining terms vanish: ( 1 m B = 4(m 1p(p 1 m Λ m Λ p m<p, m,p=1 4m(m 1p(p 1 ( 1 p 4(p 1m(m 1 p Λ m Λ p m>p, m,p=1 4p(p 1m(m 1 ( m m = m<p m(m 1p(p 1 p p m>p Λ m Λ p m(m 1p(p 1 = 0, (A1 an in quite the same manner B 3 B 4 = 1 = 1 1 m=1 p=1 m m m(m 1 Λ m ½ = 0, 1 p p ½ Λ p = 0. p(p 1 Thus we fin the following Bloch vector of B (44 B = 1 1 m=1 Λ m Λ m 1 ½ ½. (A (A3

21 1 3. Proof of Orthonormality of WOB For proofs relevant in the WOB we often nee the equivalence 1 e πi nx = n=0 { if x = 0 0 if x 0, x. (A4 So we use Eq. (A4 to proof the orthonormality (31 of the Weyl operators (30 TrU nm U lj = = = 1 1 p=0 k, k=0 1 1 p=0 k, k=0 1 k=0 e πi e πi e πi k(l n δ mj ( kl kn p (k m mo k k ( k j mo p ( kl kn p (k m mo ( k j mo p δ k k = δ nl δ mj. (A5 4. Expansion into WOB Formula (54 for the Bell state in terms of WOB we erive in the following way. We express the stanar matrices by the WOB (37, rewrite the inices an separate the nonvanishing terms φ φ = 1 j k j k j,k=1 = j,k=0 l,l =0 = m,k=0 l,l =0 e πi e πi j(ll U l(k jmo U l (k jmo (k m(ll U lm U l m ( = 1 U 0m U 0m m m ( 1 3 m l,l ; l,l 0,0; ll l,l ; ll = k U lm U l m e πi (k m(ll U lm U l m. (A6 The last term in Eq. (A6 vanishes ue to relation (A4. Ientifying U 00 = ½ an using the notation with negative values of the inex l, which have to be consiere as mo, we gain the formula φ φ = 1 ½ ½ 1 1 l,m=0 U lm U lm, (l, m (0, 0. (A7

22 [1] F. Bloch, Phys. Rev. 70, 460 (1946. [] M. Nielsen an I. Chuang, Quantum Computation an Quantum Information (Cambrige University Press, Cambrige, Englan, 000. [3] G. Kimura, Phys. Lett. A 314, 339 (003. [4] S. Kryszewski an M. Zachcia l, quant-ph/ [5] L. Jakobczyk an M. Siennicki, Phys. Lett. A 86, 383 (001. [6] G. Kimura an A. Kossakowski, Open Sys. Information Dyn. 1, 07 (005. [7] I. P. Menaš, J. Phys. A: Math. Gen. 39, (006. [8] C. Witte an M. Trucks, Phys. Lett. A 57, 14 (1999. [9] M. Ozawa, Phys. Lett. A 68, 158 (000. [10] R. A. Bertlmann, H. Narnhofer, an W. Thirring, Phys. Rev. A 66, (00. [11] R. A. Bertlmann, K. Durstberger, B. C. Hiesmayr, an P. Krammer, Phys. Rev. A 7, (005. [1] C. Rakotonirina, math.gm/ [13] C. M. Caves an G. J. Milburn, Optics Communications 179, 439 (000. [14] D. A. Varshalovich, A. N. Moskalev, an V. K. Khersonskii, Quantum theory of angular momentum (Worl Scientific Publishing Co. Pte. Lt., [15] C. H. Bennett, G. Brassar, C. Crépeau, R. Jozsa, A. Peres, an W. K. Wootters, Phys. Rev. Lett. 70, 1895 (1993. [16] H. Narnhofer, J. Phys. A: Math. Gen. 39, 7051 (006. [17] B. Baumgartner, B. C. Hiesmayr, an H. Narnhofer, Phys. Rev. A 74, 0337 (006. [18] B. Baumgartner, B. C. Hiesmayr, an H. Narnhofer, J. Phys. A: Math. Theor. 40, 7919 (007. [19] R. F. Werner, J. Phys. A: Math. Gen. 34, 7081 (001. [0] K. G. H. Vollbrecht an R. F. Werner, J. Math. Phys. 41, 677 (000. [1] E. M. Henley an W. Thirring, Elementary Quantum Fiel Theory (McGraw Hill, New York, 196. [] M. Horoecki an P. Horoecki, Phys. Rev. A 59, 406 (1999. [3] E. M. Rains, Phys. Rev. A 60, 179 (1999. [4] M. Horoecki, P. Horoecki, an R. Horoecki, in Quantum Information, eite by G. A. et al. (Springer Verlag Berlin, 001, vol. 173 of Springer Tracts in Moern Physics, p [5] M. Horoecki, P. Horoecki, an R. Horoecki, Phys. Lett. A 3, 1 (1996. [6] B. M. Terhal, Phys. Lett. A 71, 319 (000. [7] B. M. Terhal, Theoretical Computer Science 87, 313 (00. [8] D. Collins, N. Gisin, N. Linen, S. Massar, an S. Popescu, Phys. Rev. Lett 88, (00. [9] D. Kaszlikowski, L. C. Kwek, J.-L. Chen, M. Żukowski, an C. H. Oh, Phys. Rev. A 65, (00. [30] S. Gröblacher, T. Jennewein, A. Vaziri, G. Weihs, an A. Zeilinger, New J. Phys. 8, 75 (006. [31] A. Vaziri, G. Weihs, an A. Zeilinger, Phys. Rev. Lett. 89, (00. [3] R. A. Bertlmann an P. Krammer, Phys. Rev. A 77, (008.