arxiv:math/ v1 [math.fa] 4 Jan 2007

Similar documents
arxiv:math/ v2 [math.fa] 29 Mar 2007

Notes on matrix arithmetic geometric mean inequalities

Norm inequalities related to the matrix geometric mean

On a Block Matrix Inequality quantifying the Monogamy of the Negativity of Entanglement

arxiv: v1 [math.fa] 19 Jul 2009

Compound matrices and some classical inequalities

A norm inequality for pairs of commuting positive semidefinite matrices

Clarkson Inequalities With Several Operators

Abstract. In this article, several matrix norm inequalities are proved by making use of the Hiroshima 2003 result on majorization relations.

The matrix arithmetic-geometric mean inequality revisited

Multivariate trace inequalities. David Sutter, Mario Berta, Marco Tomamichel

arxiv: v1 [math.fa] 1 Apr 2007

Some inequalities for unitarily invariant norms of matrices

Singular Value and Norm Inequalities Associated with 2 x 2 Positive Semidefinite Block Matrices

Mathematical Methods for Quantum Information Theory. Part I: Matrix Analysis. Koenraad Audenaert (RHUL, UK)

Elementary linear algebra

Matrix Inequalities by Means of Block Matrices 1

MATRIX INEQUALITIES IN STATISTICAL MECHANICS

Trace Inequalities for a Block Hadamard Product

Representations of Lorentz Group

arxiv: v1 [math.fa] 1 Oct 2015

Exercise Sheet 1.

arxiv: v4 [math.sp] 19 Jun 2015

On a decomposition lemma for positive semi-definite block-matrices

arxiv: v1 [math.ra] 8 Apr 2016

Singular Value Inequalities for Real and Imaginary Parts of Matrices

Plan for the rest of the semester. ψ a

CLASSIFICATION OF COMPLETELY POSITIVE MAPS 1. INTRODUCTION

Some Inequalities for Commutators of Bounded Linear Operators in Hilbert Spaces

ON THE HÖLDER CONTINUITY OF MATRIX FUNCTIONS FOR NORMAL MATRICES

j=1 [We will show that the triangle inequality holds for each p-norm in Chapter 3 Section 6.] The 1-norm is A F = tr(a H A).

Eigenvalue inequalities for convex and log-convex functions

EE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 2

2. reverse inequalities via specht ratio To study the Golden-Thompson inequality, Ando-Hiai in [1] developed the following log-majorizationes:

arxiv: v1 [math.fa] 6 Nov 2015

Distance Preserving Maps on Matrices. Chi-Kwong Li Department of Mathematics College of William and Mary Williamsburg, Virginia

arxiv:math/ v2 [math.oa] 8 Oct 2006

Section 3.9. Matrix Norm

Matrix Energy. 1 Graph Energy. Christi DiStefano Gary Davis CSUMS University of Massachusetts at Dartmouth. December 16,

Lecture 18. Ramanujan Graphs continued

Duke University, Department of Electrical and Computer Engineering Optimization for Scientists and Engineers c Alex Bronstein, 2014

Compression, Matrix Range and Completely Positive Map

On the Generalized Reid Inequality and the Numerical Radii

Some applications of matrix inequalities in Rényi entropy

SOME INEQUALITIES FOR COMMUTATORS OF BOUNDED LINEAR OPERATORS IN HILBERT SPACES. S. S. Dragomir

SAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear Algebra

SAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear Algebra

Lecture 5. Ch. 5, Norms for vectors and matrices. Norms for vectors and matrices Why?

CHARACTERIZATIONS. is pd/psd. Possible for all pd/psd matrices! Generating a pd/psd matrix: Choose any B Mn, then

Matrix Theory, Math6304 Lecture Notes from March 22, 2016 taken by Kazem Safari

Entanglement Measures and Monotones

MATH 431: FIRST MIDTERM. Thursday, October 3, 2013.

The Improved Arithmetic-Geometric Mean Inequalities for Matrix Norms

Matrix inequalities from a two variables functional

Matrix Mathematics. Theory, Facts, and Formulas with Application to Linear Systems Theory. Dennis S. Bernstein

Journal of Inequalities in Pure and Applied Mathematics

Distribution for the Standard Eigenvalues of Quaternion Matrices

Optimisation Problems for the Determinant of a Sum of 3 3 Matrices

On Symmetric Norm Inequalities And Hermitian Block-Matrices

Introduction to Numerical Linear Algebra II

Banach Journal of Mathematical Analysis ISSN: (electronic)

arxiv:quant-ph/ v1 24 Aug 2006

Spectral inequalities and equalities involving products of matrices

arxiv:math/ v1 [math.fa] 4 Nov 2000

Mathematical Methods wk 2: Linear Operators

HERMITE-HADAMARD TYPE INEQUALITIES FOR OPERATOR GEOMETRICALLY CONVEX FUNCTIONS

RANKS OF QUANTUM STATES WITH PRESCRIBED REDUCED STATES

Math 489AB Exercises for Chapter 2 Fall Section 2.3

Lecture Notes 1: Vector spaces

Lecture 2: Linear operators

Lecture Note 5: Semidefinite Programming for Stability Analysis

FUNCTIONAL ANALYSIS LECTURE NOTES: ADJOINTS IN HILBERT SPACES

Journal of Inequalities in Pure and Applied Mathematics

arxiv: v1 [math.fa] 19 Aug 2017

Definition (T -invariant subspace) Example. Example

LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM

Interpolating the arithmetic geometric mean inequality and its operator version

A generalization of Araki s log-majorization

Lecture 1 and 2: Random Spanning Trees

Linear and non-linear programming

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2017 LECTURE 5

1 Positive and completely positive linear maps

Embedding constants of trilinear Schatten von Neumann classes

A connection between number theory and linear algebra

Throughout these notes we assume V, W are finite dimensional inner product spaces over C.

Mathematik-Bericht 2010/5. An eigenvalue estimate and its application to non-selfadjoint Jacobi and Schrödinger operators

Singular Value Decomposition and Polar Form

2. Linear algebra. matrices and vectors. linear equations. range and nullspace of matrices. function of vectors, gradient and Hessian

Extensions of interpolation between the arithmetic-geometric mean inequality for matrices

GENERALIZED NORMS INEQUALITIES FOR ABSOLUTE VALUE OPERATORS

Math Linear Algebra II. 1. Inner Products and Norms

October 25, 2013 INNER PRODUCT SPACES

POSITIVE MAP AS DIFFERENCE OF TWO COMPLETELY POSITIVE OR SUPER-POSITIVE MAPS

Complex symmetric operators

ALUTHGE ITERATION IN SEMISIMPLE LIE GROUP. 1. Introduction Given 0 < λ < 1, the λ-aluthge transform of X C n n [4]:

0.1 Rational Canonical Forms

MAT 445/ INTRODUCTION TO REPRESENTATION THEORY

The Singular Value Decomposition and Least Squares Problems

Lecture notes: Applied linear algebra Part 1. Version 2

Positive semidefinite matrix approximation with a trace constraint

Transcription:

Tr[ABA] p = Tr[B 1/2 A 2 B 1/2 ] p. On the Araki-Lieb-Thirring inequality arxiv:math/0701129v1 [math.fa] 4 Jan 2007 Koenraad M.R. Audenaert Institute for Mathematical Sciences, Imperial College London 53 Prince s Gate, London SW7 2PG, United Kingdom March 4, 2008, 3:59 Abstract We prove an inequality that complements the famous Araki-Lieb-Thirring (ALT) inequality for positive matrices A and B, by giving a lower bound on the quantity Tr[A r B r A r ] q in terms of Tr[ABA] rq for 0 r 1 and q 0, whereas the ALT inequality gives an upper bound. Secondly, we also prove a generalisation of the ALT inequality to general matrices. 1 Introduction A famous inequality with a lot of applications in mathematics and mathematical physics is the Araki-Lieb- Thirring inequality [1, 3]: Theorem 1 (Araki-Lieb-Thirring) For A, B 0, q 0, and for 0 r 1, the following inequality holds: Tr[A r B r A r ] q Tr[ABA] rq, (1) while for r 1, the inequality is reversed. In this paper we do two things. In Section 2 we obtain complementary inequalities. That is, for 0 r 1 we obtain upper bounds on Tr[ABA] rq (in terms of the quantity Tr[A r B r A r ] q ), and lower bounds for r 1. Such bounds may be useful, for example, to obtain estimates on the error incurred by going from Tr[ABA] rq to Tr[A r B r A r ] q. Second, in Section 3, we find a generalisation of the ALT inequality to general matrices. 2 A complementary inequality In this Section we want to obtain upper bounds on Tr[ABA] rq in terms of the quantity Tr[A r B r A r ] q (for 0 r 1). A remark one can make right away is that both quantities Tr[ABA] rq and Tr[A r B r A r ] q have the same degrees of homogeneity in A and B. Any upper bound on Tr[ABA] rq that only depends on Tr[A r B r A r ] q should therefore be linear in the latter quantity. Unfortunately, numerical calculations show that the required proportionality factor should be infinitely large to accomodate all possible A and B. This means that extra ingredients are needed to obtain a reasonable upper bound. We have found such an upper bound by including certain norms of A and B as additional ingredients. In fact, we have found a whole family of such bounds. The simplest bound in this family, but also the weakest, is given by Proposition 1 Let A, B 0. For q 0 and r 0, we have the upper bound while for q 0 the inequality is reversed. Tr[ABA] rq A 2rq Tr B rq, (2) Here,. denotes the operator norm, which for positive semidefinite (PSD) matrices is nothing but the largest eigenvalue. Proof. Put p = rq 0. Note first that

From the basic inequality A 2 A 2 11 = A 2 11 follows B 1/2 A 2 B 1/2 A 2 B. If p is between 0 and 1, we may take the p-th power of both sides (x x p is then operator monotone). Taking the trace of both sides then yields (2). If p is larger than 1, we may take the Schatten p-norm of both sides, by Weyl-monotonicity of unitarily invariant (UI) norms. Taking the p-th power of both sides again yields (2). For reasons that will immediately become clear, we call this inequality the water -inequality, to express the fact that it is rather weak, and not very spiritual. In contrast, we call the Araki-Lieb-Thirring inequality the wine -inequality, because it is too strong for our purposes: it gives a lower bound, rather than an upper bound. We can obtain better upper bounds by cutting the wine with the water. Fixing A and B, some t obviously must exist, with 0 t 1, such that the following holds: Tr[ABA] rq (water) t (wine) 1 t = ( A 2rq TrB rq) t (Tr[A r B r A r ] q ) 1 t. Of course, this would be a rather pointless (and disappointing) exercise if there were some A and B for which the smallest valid value of t would be 1, because then inequality (2) would be the only upper bound valid for all A and B. Fortunately, numerical experiments revealed the fact (which we will prove below) that here any value of t between 1 r and 1 yields an upper bound, for any A and B. This yields the promised family of inequalities, of which the sharpest and most relevant one is the one with t = 1 r. Theorem 2 For A, B 0, q 0 and 0 r 1, For r 1, the inequality is reversed. Tr[ABA] rq ( A 2rq TrB rq (Tr[A r B r A r ] q ) r. (3) This inequality is sharp, just like the original ALT inequality, as can be seen by taking scalar A and B. An equivalent formulation of inequality (3) is Yet another formulation is obtained if one notes the equality by which we get (ABA) r q ( A 2r B r q A r B r A r r q. (4) (ABA) r q = ABA r rq = B 1/2 A 2 B 1/2 r rq = (B 1/2 A 2 B 1/2 ) r q, (ABA) r q = (B 1/2 A 2 B 1/2 ) r q ( B 1/2 2r A 2r q B r/2 A 2r B r/2 r q = ( B r A 2r q A r B r A r r q. (5) It is this last formulation that we will consider in the following proof. Proof. Let first 0 r 1. Then, for X 0, X 1 r = X 1 r. Since B 0 we can write B = B r B 1 r B 1 r B r. Thus ABA B 1 r AB r A = B 1 r A 1 r (A r B r A r )A 1 r. (6) Consider first the easiest case q =. Taking the operator norm of both sides of (6) then gives ABA B 1 r A 1 r (A r B r A r )A 1 r B 1 r A 1 r 2 A r B r A r, where the last inequality follows from submultiplicativity of the operator norm. One obtains (5) for q = by taking the r-th power of both sides, and noting (again) X r = X r.

To prove (5) for general q, let us first take the r-th power of both sides of (6), which preserves the ordering because of operator monotonicity of x x r for 0 r 1: (ABA) r B r 1 r ( A 1 r (A r B r A r )A 1 r) r. Thus, on taking the q-norm of both sides (or q-quasinorm if 0 < q < 1), (ABA) r q B r 1 r ( A 1 r (A r B r A r )A 1 r) r q. (7) We can now apply a generalisation of Hölder s inequality, by which for all positive real numbers s, t, u such that 1/s + 1/t = 1/u we have ([2], Eq. (IV.43)) XY u 1/u X s 1/s Y t 1/t, for all X, Y and for all UI norms.. In fact, this inequality extends to UI quasinorms like the Schatten q-quasinorms for 0 < q < 1. For X, Y 0, two successive applications of this inequality yield (XY X) u 1/u X 2s 1/s Y t 1/t. We apply the latter inequality to the second factor of the right-hand side (RHS) of (7), with the substitutions X = A 1 r, Y = A r B r A r, u = r, s = r/(1 r) (the positivity of which requires r to lie between 0 and 1), t = 1, and. =. q, giving ( A 1 r (A r B r A r )A 1 r) r 1/r q A 2r (1 r)/r q A r B r A r q. Taking the r-th power of both sides and substituting in (7) yields (5). The case r 1 follows very easily from the case 0 r 1 by making in (3) the substitutions A = A r, B = B r, r = 1/r, q = qr. Taking the r -th power of both sides, rearranging factors, and subsequently dropping primes yields (3) for r 1. As a special case of Theorem 2 we consider the comparison between Tr[AB] and AB 1, for A, B 0. The lower bound AB 1 Tr[AB] is well-known and easy to prove. Indeed, since A 1/2 BA 1/2 is normal, Tr[AB] = Tr[A 1/2 BA 1/2 ] = A 1/2 BA 1/2 1 A 1/2 A 1/2 B 1 = AB 1. To obtain an upper bound, (3) with r = 1/2 and q = 1 yields AB 1 = Tr(AB 2 A) 1/2 ( A Tr B Tr(A 1/2 BA 1/2 ) ) 1/2 = ( A TrB Tr(AB)) 1/2. (8) 3 Generalisations of the Araki-Lieb-Thirring inequality In this Section we want to find a generalisation of the ALT inequality to general matrices A and B. As a first step, we generalise it to the case of general A while keeping B 0. Proposition 2 For any matrix A and B 0, for q 1, and for any UI norm, (ABA ) q A q B q A q. (9)

Using the facts X Y p p = X p p+ Y p p and X q p = X q p, yields the statement of the Theorem. Proof. Let the polar decomposition of A be A = U A, where A = (A A) 1/2 is the modulus of A and U is unitary. Then, for any UI norm, (ABA ) q = (U A B A U ) q = ( A B A ) q A q B q A q, where in the last step we used the ALT inequality proper. The second step is to generalise this statement to the case where B is Hermitian. First we need a Lemma. Lemma 1 For X, Y 0, and any UI norm, X Y X + Y. ( ) ( ) X 0 X + Y X Y Proof. The matrix is unitarily equivalent with /2. Since X, Y 0, the 0 Y X Y X + Y latter matrix is also PSD, whence there exists a contraction K such that X Y = (X + Y ) 1/2 K(X + Y ) 1/2. Moreover, as X Y and X + Y are Hermitian, so is K. Therefore X Y = (X + Y ) 1/2 K(X + Y ) 1/2 (X + Y )K, where we have used the fact that AB BA whenever AB is normal. Every contraction K can be written as a convex combination of unitaries, so by convexity of norms, and the fact that we re considering UI norms, we have (X + Y )K X + Y. Proposition 3 For any matrix A and Hermitian B, for q 1, and for any UI norm, ABA q A q B q A q. (10) Proof. Let the Jordan decomposition of B be B = B + B, where B + and B (the positive and negative part, respectively) are both PSD. Then the two terms in the right-hand side of ABA = AB + A AB A are also PSD. Thus ABA q (AB + A + AB A ) q. This follows from the Lemma applied to the norm q 1/q. But as AB + A + AB A = A B A, we get ABA q (A B A ) q A q B q A q, where in the last step we used Proposition 2. If we now specialise to Schatten p-norms, we can drop the conditions on B: Theorem 3 For general matrices A and B, and for p, q 1, ABA q p A q B q + B q p A q. (11) 2 ( ) ( A 0 0 B Proof. We use Proposition 3 with A = 0 A and B = B 0 ( ) ( ) 0 X X X = 0, 0 0 X ). Noting that this yields, after dropping primes, ( ) ( AB A q 0 0 ABA q A q B q A q 0 p 0 A q B q A q ) ). p

4 Acknowledgements This work was supported by the Institute for Mathematical Sciences, Imperial College London, and is part of the QIP-IRC (www.qipirc.org) supported by EPSRC (GR/S82176/0). References [1] H. Araki, On an inequality of Lieb and Thirring, Lett. Math. Phys. 19, 167 170 (1990). [2] R. Bhatia, Matrix Analysis, Springer, Heidelberg (1997). [3] E. Lieb and W. Thirring, in E. Lieb, B. Simon and A. Wightman (eds), Studies in Mathematical Physics, pp. 301 302, Princeton Press (1976).