Fidelity and trace norm

Transcription

1 Fidelity and trace norm In this lecture we will see two simple examples of semidefinite programs for quantities that are interesting in the theory of quantum information: the first is the trace norm of a given operator and the second is the fidelity between two positive semidefinite operators. When considering both of these semidefinite programming formulations, we will make use of a lemma that will come up several more times throughout the course; it gives a precise characterization of the relationship between the blocks of any 2-by-2 positive semidefinite block operator. 4. Trace norm Suppose X and Y are complex Euclidean spaces and A LY, X is an arbitrary operator. By the singular value decomposition theorem, one may always write A = r s k x k y k 4. k= for r = ranka and for some choice of orthonormal sets {x,..., x r } X and {y,..., y r } Y and positive real numbers s,..., s r which we always assume are ordered from largest to smallest, so s s 2 s r > 0. The values s,..., s r are the singular values of A, and up to their ordering are uniquely determined by A. The trace norm of A, which is denoted A, is defined as the sum of the singular values of A: A = s + + s r. 4.2 There are other, equivalent expressions of the trace norm, such as A = Tr AA = Tr A A 4.3 and A = max { B, A : B LY, X, B }. 4.4

2 CS 867/IC 890 Semidefinite Programming in uantum Information The second expression may alternatively be written as A = max { Re B, A : B LY, X, B } ; 4.5 it always holds that Re B, A B, A, and for an arbitrary choice of B one can choose a complex number γ with γ =, and therefore γb = B, such that Re γb, A = B, A. We ll write down an SDP for the trace norm of a given operator A momentarily, but first let us state and prove the lemma about block operators that was suggested above. Lemma 4.. Let X and Y be complex Euclidean spaces, let P PosX and PosY be positive semidefinite operators, and let X LY, X be an operator. It holds that P X X PosX Y 4.6 if and only if X = PK for K LY, X satisfying K. Proof. Suppose first that X = PK for K LY, X being an operator for which K. It follows that KK X, and therefore 0 PK K P = For the reverse implication, assume PKK P X P X X X. 4.7 P X X PosX Y, 4.8 and define K = P + X Here P + and + refer to the Moore Penrose pseudo-inverses of P and, respectively. For a case such as this, where P and are positive semidefinite, they may be obtained by taking a spectral decomposition of either P or and replacing each nonzero eigenvalue by its reciprocal. It will be proved that X = PK and K. Observe first that, for every Hermitian operator H HermX, the block operator H 0 P X H 0 HPH HX 0 X = 0 X H 2 4.0

3 is positive semidefinite. In particular, for H = Π kerp being the projection onto the kernel of P, one has that the operator 0 Π kerp X X 4. Π kerp is positive semidefinite, which implies that Π kerp X = 0, and therefore Π imp X = X. Through a similar argument, one finds that XΠ im = X. It therefore follows that PK = ΠimPXΠ im = X. 4.2 Next, note that x Px x Xy x 0 P X x 0 = y X x y y 0 y X y for every choice of vectors x X and y Y. Setting x = P + u and y = + v 4.4 for arbitrarily chosen unit vectors u X and v Y, one finds that u Kv u Π v K imp u u Kv u v K u v Π im v and therefore u Kv. As this inequality holds for all unit vectors u and v, it follows that K, as required. Now we can devise an SDP whose optimal value is A for a given operator A LY, X. The primal problem, in a slightly simplified form, will be the following: maximize: Primal problem 2 K, A + 2 K, A X K 0, K Y K LY, X. The optimal value of this primal form is evidently A, because the objective function can alternatively be expressed as Re K, A, and by the lemma above K is free to range over all elements of LY, X with spectral norm at most. 3

4 CS 867/IC 890 Semidefinite Programming in uantum Information The primal problem above can be matched to the rigid formal definition of semidefinite programs we have adopted by first defining Φ TX Y as X X 0 Φ = 4.6 Y 0 Y for all X LX and Y LY and the dots representing unnamed elements of LY, X or LX, Y that get zeroed out by the map. The primal problem above is then equivalent to the primal problem of the semidefinite program Φ, 2 0 A A 0, X 0 0 Y. 4.7 The dual problem is particularly easy to calculate for this semidefinite program because Φ is its own adjoint. After just a bit of simplification we obtain the following problem: minimize: Dual problem TrX + TrY X 0 0 Y 2 0 A A 0, X HermX, Y HermY. We can further simplify this problem by first noting that we must have X PosX and Y PosY for every feasible X and Y; and by setting P = 2X and = 2Y we get the following problem statement: minimize: Dual problem 2 TrP + 2 Tr P A A 0, P PosX, PosY. If you think about a singular value decomposition r A = s k x k y k 4.8 k= of A, then the dual problem makes perfect sense: we can take P = r s k x k x r k and = s k y k y k 4.9 k= k= to obtain a feasible solution whose objective value is A, and by weak duality this is the best we can do. 4

5 4.2 Fidelity The fidelity between two positive semidefinite operators P, PosX is defined as FP, = P Sometimes the fidelity is defined as this quantity squared, but we ll use the definition above, without the square. The alternative definition P FP, = Tr P 4.2 is obtained by expanding the trace norm using the expression 4.3. We will now describe a semidefinite program whose optimal value coincides with FP,. While this task is easily accomplished by applying the semidefinite program for the trace norm to the operator P, we will aim for something simpler that does not require the computation of operator square roots. The primal problem of our semidefinite program will be as follows: maximize: Primal problem 2 TrX + 2 TrX P X X 0, X LX. Notice that the objective function can alternatively be expressed as ReTrX. The lemma proved in the previous section reveals that the optimal value of the primal problem is the fidelity between P and : the operator X ranges precisely over those operators equal to PK for K LX satisfying K, and taking the maximum over all such K yields { } max Re Tr PK : K LX, K { } = max Re P, K : K LX, K { = max Re K, P } 4.22 : K LX, K = P. The primal problem above can be matched to the primal problem of a formally specified semidefinite program in a way that is quite similar to the trace norm 5

6 CS 867/IC 890 Semidefinite Programming in uantum Information semidefinite program from earlier in the lecture. More specifically, we may define a map Φ TX X exactly as in 4.6, taking Y = X so that each block is an element of LX, and then observe that the primal problem above is in agreement with the primal problem for the semidefinite program Φ, 2 0, 0 P The dual problem, after just a bit of simplification similar to the dual problem in the trace norm semidefinite program, is as follows: Dual problem minimize: 2 Y, P + Z, 2 Y 0, Z Y, Z LX. The primal problem is feasible although not strictly feasible in case P or is singular and the dual is strictly feasible, so strong duality holds by Slater s theorem. This dual problem can simplified further using the following proposition. Proposition 4.2. Let Y, Z HermX. It holds that if and only if Y, Z > 0 and Z Y. Y PosX X 4.24 Z Proof. Suppose Y, Z > 0 and Z Y. It holds that Y = Z 0 Y 0 Y 0 Z Y Y and therefore Y PosX X Z Conversely, suppose that Y PosX X Z 6

7 It holds that 0 u v Y Z u = u Yu u v v u + v Zv 4.28 v for all u, v X. If Y were not positive definite, there would exist a unit vector v for which v Yv = 0, and one could then set to obtain u = Z + v Z v Zv u, v + v, u = Z +, 4.30 which is absurd. Therefore Y > 0, and by inverting the expression above, we have Y 0 0 Y Y 0 Z Y = Y Z 0 This implies Z Y and therefore Z > 0 as required. Now, given that is positive semidefinite, it holds that, Z, Y whenever Z Y, so there would be no point in choosing any Z other than Y when aiming to minimize the dual objective function subject to that constraint. The dual problem above can therefore be phrased as follows: Dual problem minimize: 2 Y, P + 2 Y, Y > 0, Y PosX. Because we have strong duality for our semidefinite program, we obtain the following alternative characterization of the fidelity: FP, = inf Y>0 2 Y, P + 2 Y, This characterization is equivalent to Alberti s theorem, which states that FP, 2 = inf Y>0 Y, P Y, The equivalence can be established through the use of the arithmetic-geometric mean inequality. The following two theorems establish well-known properties of the fidelity function, but through proofs based on the semidefinite programming characterization we have obtained. 7

8 CS 867/IC 890 Semidefinite Programming in uantum Information Theorem 4.3. Let P 0, P, 0, PosX be positive semidefinite operators, for X a complex Euclidean space. It holds that FP 0 + P, 0 + FP 0, 0 + FP, Proof. By considering the primal problem associated with our semidefinite program, and noting that Slater s theorem implies that an optimal primal solution is always obtained, we see that we may choose operators X 0, X LX such that the block operators P0 X 0 P X X0 and 0 X 4.35 are both positive semidefinite, and such that TrX 0 = FP 0, 0 and TrX = FP, The sum of two positive semidefinite operators is positive semidefinite, and therefore P0 + P X 0 + X P0 X X 0 + X = 0 P X 0 + X0 + 0 X 4.37 is positive semidefinite. By again considering the primal problem of our semidefinite program, this time for FP 0 + P, 0 +, we find that as required. FP 0 + P, 0 + TrX 0 + X = FP 0, 0 + FP,, 4.38 Note that by combining this theorem with the fact that FλP, λ = λ FP,, we find that the fidelity is jointly concave: FλP 0 + λp, λy 0 + λy λ FP 0, Y 0 + λ FP, Y Theorem 4.4. Let X and Y be complex Euclidean spaces, let P, PosX, and let Φ CX, Y be a channel. It holds that FP, FΦP, Φ Proof. One may choose X LX so that P X X 4.4 8

9 is positive semidefinite and satisfies ReTrX = FP,. By the complete positivity of Φ, the block operator ΦP ΦX ΦP ΦX ΦX = Φ ΦX Φ 4.42 is positive semidefinite as well. By again considering the primal problem of the semidefinite program, this time for FΦP, Φ, along with the fact that Φ is trace-preserving, we find that as required. FΦP, Φ ReTrΦX = ReTrX = FP,,