Some Bipartite States Do Not Arise from Channels arxiv:quant-ph/0303141v3 16 Apr 003 Mary Beth Ruskai Department of Mathematics, Tufts University Medford, Massachusetts 0155 USA marybeth.ruskai@tufts.edu November 6, 018 Dedicated to Charles H. Bennett on his 60th Birthday Abstract It is well-known that the action of a quantum channel on a state can be represented, using an auxiliary space, as the partial trace of an associated bipartite state. Recently, it was observed that for the bipartite state associated with the optimal average input of the channel, the entanglement of formation is simply the entropy of the reduced density matrix minus the Holevo capacity. It is natural to ask if every bipartite state can be associated with some channel in this way. We show that the answer is negative. PACS number 03.67; MSC classification 8P68. 1 Background Recently, Matsumoto, Shimono and Winter (MSW) [4] pointed out an important connection between the channel capacity and of entanglement of formation, which allows one to draw some conclusion about the additivity of the latter from that of the former. There has also been speculation that the additivity of channel capacity and of entanglement of formation, are equivalent. A connection between additivity of entanglement of formation and multiplicativity of the p-norm measure of purity has also been given by Audenaert and Braunstein [1]. Work partially supported by the U.S. National Science Foundation under Grant number DMS-03-148. 1
In view of this it is natural to ask if every bipartite state ρ AB can be associated with a channel in the sense of MSW. In particular is its reduced density matrix, ρ A, the optimal average output of a channel whose capacity is related to the entanglement of formation of ρ AB as described in [4]? We will give a more precise formulation of this statement and show that the answer is negative. Recall that a state is represented by a density matrix ρ, i.e., a positive semidefinite operator with Trρ = 1. The von Neumann entropy of ρ is S(ρ) = Trρlogρ. By an ensemble E = {π i,ρ i } we means a set of density matrices, ρ i and associated probabilities π i. By a channel Φ we mean a completely positive, trace-preserving (CPT) map. The Holevo capacity of the map Φ is C Holv (Φ) = sup E { S[Φ(ρ)] i } π i S[Φ(ρ i )] where the supremum is taken over all ensembles E = {π i,ρ i } and ρ = i π iρ i. The optimal average input is the state ρ opt = i π iρ i associated with the ensemble which attains the supremum in (1). The optimal average output is then Φ(ρ opt ). The key to the MSW construction is the following result of Stinespring [7], which was subsequently used by Lindblad [3] in his work on relative entropy. Theorem 1 Given a CPT map Φ on B(H), one can find an auxiliary space H B, a density matrix ν B and a unitary map U AB on H A H B such that [ Φ(ρ) = Tr B U AB ρ ν ] BU AB () where Tr B denotes the partial trace and we have identified H with H A. Although we are interested in the case in which ρ is a density matrix, the representation () is valid for all operators in B(H). When Φ(ρ) = k A k ρa k, it is easytoconstructarepresentationoftheform(),asdiscussedinsectioniii.dof[6]. In fact, U AB ρ ν BU AB is the block matrix which has the form jk A j ρa k j k. Given a bipartite state γ AB on H A H B, its entanglement of formation satisfies { EoF(γ AB ) = inf π k S ( ) Tr B γ k : {πk γ k } ensemble with γ AB = π k γ k }. (3) k k Although it is customary to define the EoF using ensembles for which all γ k are pure states, there is no loss of generality in allowing arbitrary states. Theorem (Matsumoto, Shimono and Winter) Let Φ be a CPT map with a representation H B,ν B,U AB as in Theorem 1. The Holevo capacity of Φ satisfies C Holv (Φ) = sup {S ( ) } Tr B γ AB EoF(γAB ) : γ AB = U AB ρ ν BU AB (4) where ρ is a density matrix on H = H A. Moreover, the state γ AB which attains this supremum satisfies Tr B γ AB = Φ(ρ opt ) where ρ opt is the optimal input of Φ. (1)
It follows immediately that for the state γ AB It is thus natural to ask the following EoF( γ AB ) = S [ Φ(ρ opt ) ] C Holv (Φ). (5) Question 3 Given a bipartite state γ AB on H A H B, is there a CPT map Φ on B(H A ) such that Tr B γ AB is the optimal average output state of Φ and EoF(γ AB ) = S ( Tr B γ AB ) CHolv (Φ)? (6) We will show that the answer to this question is negative. Our counter-example is based on a result in [5] and [9] which comes from the fact that when ρ = k π kρ k S(ρ) k π k S(ρ k ) = k π k H(ρ k,ρ) (7) where the relative entropy is defined as H(ω,ρ) = Trω [ logω logρ ]. (8) Theorem 4 Let Φ be a CPT map with optimal average input ρ opt. Then C Holv (Φ) = suph[φ(ω),φ(ρ opt )]. (9) ω It follows as an immediate corollary, that if E = {π k,ρ k } is any optimal ensemble for the channel Φ, then H[Φ(ρ k ),Φ(ρ opt )] = C Holv (Φ) is independent of k, i.e., the optimal average output is equi-distant in the sense of relative entropy from all outputs Φ(ρ k ) in the ensemble. Counter-example An affirmative answer to Question 3 above would imply that if the ensemble {π k, ψ k ψ k } is optimal for EoF(γ AB ), then {π i,tr B ψ k ψ k } would be an optimal output ensemble for the corresponding map Φ. It would then follow from the equi-distance corollary to Theorem 8 that H(ω k,γ A ) = S(ω k ) Tr A ω k logγ A = C. (10) where ω k = Tr B ψ k ψ k and C is a constant which is independent of k. (In fact, C is the Holevo capacity of Φ if such a channel exists.) We will present an example of a bipartite qubit state which does not satisfy (10). 3
It may be that counter-examples to (10) can already be found in the literature on entanglement. However, we will take advantage of a result of Wootters [10] to construct a rather simple qubit counter-example which does not require extensive numerical computation. In [10], Wootters obtained an explicit formula for the EoF of a bipartite qubit state using a quantity called the concurrence. Moreover, he showed that the EoF could be achieved using an ensemble of at most four pure states ψ k ψ k, forwhichs(tr B ψ k ψ k ) = EoF(γ AB ). Thus, forsuchanensemble, (10) holds if and only if Tr A [(Tr B ψ k ψ k )logγ A ] (11) is independent of k. We can assume without loss of generality that γ A is diagonal with eigenvalues 1 (1±x). Then Trωlogγ A depends only on the diagonal elements of ω which can be written as 1 (1±d). In fact, Trωlogγ A = 1+d log 1+x + 1 d log 1 x = log 1 [ log(1 x )+dlog 1+x 1 x (1) ] which depends linearly on d for fixed x 0. Thus, when γ A 1 I, (11) is independent of k if and only if all ω k = Tr B ψ k ψ k have the same diagonal elements. However, we also know that S(ω k ) = EoF(γ AB ) is independent of k, which implies that all ω k have the same eigenvalues. ( Thus, ) all the reduced density matrices 1+d e Tr B ψ k ψ k must have the form 1 iθ k t e iθ for some fixed t. It is not hard k t 1 d to find an example for which this does not hold. Let β 0 = 1 ( 00 + 11 ) and β 3 = (σ z I) β 0 = 1 ( 00 11 ) be the indicated maximally entangled Bell states, and γ AB = 5 β 8 0 β 0 + 1 β 16 3 β 3 + 1 01 01 + 1 10 10. (13) 4 16 The reduced density matrix is γ A = Tr B γ = 1(I + 3 σ 16 z) 1 I. One can show that the optimal EoF decomposition of (13) has the form γ AB = 4 k=1 π k ψ k ψ k with ψ 1 = a 0 β 0 +a 3 β 3 sothatω 1 = Tr B ( ψ 1 ψ 1 )isdiagonalinthebasis 0, 1 (i.e, t = 0 in the matrix above). However, the remaining ψ k are superpositions which contain Bell states β k and the product states 01, 10 in a form which necessarily yields a reduced density matrix ω k which is not diagonal. By the discussion above, this implies a negative answer to Question 3. To actually find the entanglement of formation and optimal decomposition of γ AB, let γ = 5 β 8 0 β 0 + 1 β 16 3 β 3 + 1 10 10 + 1 01 01 (14) 4 16 4
be the density matrix with all spins flipped. The concurrence µ can be expressed in terms of the eigenvalues of ( γ γ γ) 1/. Following Wootters [10], one finds µ = 5 8 1 16 1 8 1 8 = 5 16. Let h(x) = 1+x log 1+x 1 x log 1 x. (15) Then, proceeding as described in [10], one finds EoF(γ AB ) = h( ( 31 ) ( 1 µ ) = h 0.1689 < 0.9745 h The optimal ensemble has weights 16 π 1 = 0.157, π = π 3 = π 4 = 0.84 associated with the projections for the following pure states ψ 1 = 0.8101 β 0 +0.5863 β 3 ψ = 0.7870 β 0 +0.1087 β 3 +0.543 01 +0.716 10 3 16 ) = S(γ A ). ψ 3 = 0.7870 β 0 0.1087 β 3 +0.543e iπ/3 01 +0.716e iπ/3 10 ψ 4 = 0.7870 β 0 +0.1087 β 3 0.543e iπ/3 01 0.716e iπ/3 10 One can easily verify that the diagonal elements of ω 1 = Tr B ψ 1 ψ 1 are not equal to those of ω k = Tr B ψ k ψ k for k =,3,4. One can also compute the reduced density matrices ω k and see that ω 1 is diagonal in the basis 0, 1, but the others are not. Alternatively, one can observe that if a reduced density matrix is diagonal in the basis 0, 1, then any purification must have the form with φ 1, φ orthogonal. By rewriting Ψ = a 0 φ 1 +b 1 φ. (16) ψ = 0.4796 00 0.6333 11 +.543 01 +.716 10 = 0 ( 0.4796 0.543 1 ) + 1 (.716 0 0.6333 1 ), one sees that ψ does not have the form (16) with orthogonal ( φ k. The actual ) reduced density matrices have the form ω 1 = [ 1 I + 31 σ ] 0.9750 0 16 z =, ( ) 0 0.050 0.551 e iθ k0.4743 ω k = e iθ with θ k 0.4743 0.4749 = π,θ 3 = π,θ 3 4 = π. 3 5
3 Remarks on representations In the canonical method of constructing a representation of the form (), the reference state ν B is pure and d B the dimension of H B is equal to the number of Kraus operators. Then the lifted state γ AB = U AB ρ ν B U AB for which Φ(ρ) = Tr Bγ AB has rank at most d = d A, the dimension of the original Hilbert space. This is clearly a very restricted class of bipartite states. Moreover, generically, d B > d since many maps require the maximum number d of Kraus operators. Thus the canonical representation yields only bipartite states which are far more singular than typical states. Some maps Φ(ρ) can be represented in the form () using a mixed reference state ν B. Indeed, given a mixed bipartite state ν B and unitary operator U AB, () can be used to define a channel Φ. Unfortunately, relatively little is known about such mixed state representations. It has been suggested that one might be able to represent a channel Φ in the form () using a space of dimension d B = d if mixed states are used. However, this is known not to be true in general [8, 11]. Now suppose that a qubit channel Φ can be represented (possibly using a mixed referencestateν B )intheform()usinganauxiliaryqubitspaceh B sothatd B =. Then theargument given before(11)canbeusedtoshow thatfortheoptimal input distribution {π k ρ k } TrΦ(ρ k )logφ(ρ opt ) (17) must be independent of k. One can easily check that the 3-state examples given in [] do not have this property. Most non-unital qubit maps for which the translation of the image of the Bloch sphere lacks symmetry will also violate (17). This implies that such channels require an auxiliary space with d B >. One expects d B = 4, consistent with the fact that such maps also require 4 Kraus operators. This gives another, somewhat indirect, proof that some CPT maps require an auxiliary space with d B > d. It would be of some interest to characterize the maps which admit mixed state representations with d B = d, as well as the bipartite states corresponding to the optimal average output. Acknowledgment: It is a pleasure to thank Professor Christopher King for stimulating and helpful discussions. 6
References [1] K. M.R. Audenaert and S. L. Braunstein On strong superadditivity of the entanglement of formation quant-ph/0303045. [] C. King, M. Nathanson and M.B. Ruskai Qubit channels can require more than two inputs to achieve capacity Phys. Rev. Lett. 88, 057901 (00). [3] G. Lindblad, Completely positive maps and entropy inequalities Commun. Math. Phys. 40, 147-151 (1975). [4] K. Matsumoto, T. Shimono and A. Winter, Remarks on additivity of the Holevo channel capacity and of the entanglement of formation quantph/00614 [5] M. Ohya, D. Petz and N. Watanabe, On capacities of quantum channels Prob. Math. Stats. 17, 170 196 (1997). [6] M.B. Ruskai, Inequalities for quantum entropy: a review with conditions for equality J. Math. Phys. 43, 4358 4375 (00). [7] W.F. Stinespring, Positive functions on C -algebras Proc. Amer. Math. Soc. 6, 11 16 (1955). [8] B.M. Terhal, I.L. Chuang, D.P. DiVincenzo, M. Grassl, J.A. Smolin Simulating quantum operations with mixed environments Phys.Rev. A60, 881 884 (1999). [9] B. Schumacher and M. D. Westmoreland, Optimal signal ensembles Phys. Rev. A 63, 0308 (001). [10] W. Wootters, Entanglement of formation of an arbitrary state of two qubits Phys. Rev. Lett. 80, 45 48 (1998). [11] C. Zalka and E. Rieffel Quantum operations that cannot be implemented using a small mixed environment J. Math. Phys. 43, 4376 4381 (00). 7