arxiv: v3 [math.ra] 22 Aug 2014

Similar documents
b jσ(j), Keywords: Decomposable numerical range, principal character AMS Subject Classification: 15A60

Spectral inequalities and equalities involving products of matrices

Maximizing the numerical radii of matrices by permuting their entries

for all A M m and B M n Key words. Complex matrix, linear preserver, spectral norm, Ky Fan k-norm, Schatten p-norm, tensor product.

CHI-KWONG LI AND YIU-TUNG POON. Dedicated to Professor John Conway on the occasion of his retirement.

LINEAR MAPS PRESERVING NUMERICAL RADIUS OF TENSOR PRODUCTS OF MATRICES

Some inequalities for sum and product of positive semide nite matrices

RANKS OF QUANTUM STATES WITH PRESCRIBED REDUCED STATES

Matrix Inequalities by Means of Block Matrices 1

Linear Operators Preserving the Numerical Range (Radius) on Triangular Matrices

Central Groupoids, Central Digraphs, and Zero-One Matrices A Satisfying A 2 = J

FACTORING A QUADRATIC OPERATOR AS A PRODUCT OF TWO POSITIVE CONTRACTIONS

arxiv: v4 [quant-ph] 8 Mar 2013

Generalized Interlacing Inequalities

Uniqueness of the Solutions of Some Completion Problems

Recent study in preserver problems. Chi-Kwong Li Department of Mathematics College of William and Mary Williamsburg, VA 23185

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

Convexity of the Joint Numerical Range

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

Math 408 Advanced Linear Algebra

A Note on Eigenvalues of Perturbed Hermitian Matrices

CANONICAL FORMS, HIGHER RANK NUMERICAL RANGES, TOTALLY ISOTROPIC SUBSPACES, AND MATRIX EQUATIONS

Interlacing Inequalities for Totally Nonnegative Matrices

A new criterion and a special class of k-positive maps

arxiv: v2 [quant-ph] 5 Dec 2013

PRODUCT OF OPERATORS AND NUMERICAL RANGE

an Off-Diagonal Block of

arxiv: v1 [math.co] 10 Aug 2016

Linear and Multilinear Algebra. Linear maps preserving rank of tensor products of matrices

A survey on Linear Preservers of Numerical Ranges and Radii. Chi-Kwong Li

PROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS

Isometries Between Matrix Algebras

Fuzhen Zhang s Publication and Presentation List (updated Dec. 2007)

Decomposable Numerical Ranges on Orthonormal Tensors

DAVIS WIELANDT SHELLS OF NORMAL OPERATORS

Positive definite preserving linear transformations on symmetric matrix spaces

On the Schur Complement of Diagonally Dominant Matrices

Schur Multiplicative Maps on Matrices

SOME INEQUALITIES FOR THE KHATRI-RAO PRODUCT OF MATRICES

Yimin Wei a,b,,1, Xiezhang Li c,2, Fanbin Bu d, Fuzhen Zhang e. Abstract

AMS526: Numerical Analysis I (Numerical Linear Algebra)

Jordan Canonical Form of A Partitioned Complex Matrix and Its Application to Real Quaternion Matrices

The automorphism group of separable states in quantum information theory. University of Illinois at Chicago, Chicago, Illinois ,

GENERALIZED NUMERICAL RANGES AND QUANTUM ERROR CORRECTION

Massachusetts Institute of Technology Department of Economics Statistics. Lecture Notes on Matrix Algebra

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

An Update on a Few Permanent Conjectures

1 Linear Algebra Problems

arxiv: v1 [math.co] 9 Aug 2016

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

Matrices A brief introduction

Fuzhen Zhang s Publication List

Some inequalities for unitarily invariant norms of matrices

1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i )

Clarkson Inequalities With Several Operators

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

arxiv: v2 [math.fa] 26 Aug 2014

Schur complements and matrix inequalities in the Löwner ordering

Singular Value Inequalities for Real and Imaginary Parts of Matrices

ELLIPTICAL RANGE THEOREMS FOR GENERALIZED NUMERICAL RANGES OF QUADRATIC OPERATORS

Lecture 15 Review of Matrix Theory III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Fuzhen Zhang s Publication and Presentation List

We first repeat some well known facts about condition numbers for normwise and componentwise perturbations. Consider the matrix

Matrices A brief introduction

Preservers of eigenvalue inclusion sets

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

ON WEIGHTED PARTIAL ORDERINGS ON THE SET OF RECTANGULAR COMPLEX MATRICES

MATH2210 Notebook 2 Spring 2018

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

MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2

Banach Journal of Mathematical Analysis ISSN: (electronic)

Algebra C Numerical Linear Algebra Sample Exam Problems

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

Norm Inequalities of Positive Semi-Definite Matrices

University, Fort Lauderdale, FL, USA. Published online: 30 Jun To link to this article:

ON WEIGHTED PARTIAL ORDERINGS ON THE SET OF RECTANGULAR COMPLEX MATRICES

Quantum Computing Lecture 2. Review of Linear Algebra

LinGloss. A glossary of linear algebra

Ranks and determinants of the sum of matrices from unitary orbits

Numerical Ranges and Dilations. Dedicated to Professor Yik-Hoi Au-Yeung on the occasion of his retirement

Trace Inequalities for a Block Hadamard Product

THREE LECTURES ON QUASIDETERMINANTS

Lecture notes on Quantum Computing. Chapter 1 Mathematical Background

Symmetric Norm Inequalities And Positive Semi-Definite Block-Matrices

Multiple eigenvalues

ACI-matrices all of whose completions have the same rank

Foundations of Matrix Analysis

Lecture 2 INF-MAT : , LU, symmetric LU, Positve (semi)definite, Cholesky, Semi-Cholesky

Introduction to Numerical Linear Algebra II

The skew-symmetric orthogonal solutions of the matrix equation AX = B

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2017 LECTURE 5

Operators with numerical range in a closed halfplane

The matrix arithmetic-geometric mean inequality revisited

The Ultimate Estimate of the Upper Norm Bound for the Summation of Operators

arxiv: v1 [math.na] 1 Sep 2018

Some Convexity Features Associated with Unitary Orbits

1 Determinants. 1.1 Determinant

Multiplicative Perturbation Bounds of the Group Inverse and Oblique Projection

Matrix Theory, Math6304 Lecture Notes from October 25, 2012

Mathematics. EC / EE / IN / ME / CE. for

Transcription:

arxiv:1407.0331v3 [math.ra] 22 Aug 2014 Positivity of Partitioned Hermitian Matrices with Unitarily Invariant Norms Abstract Chi-Kwong Li a, Fuzhen Zhang b a Department of Mathematics, College of William and Mary Williamsburg, VA 23187, USA b Division of Math, Science, and Technology, Nova Southeastern University 3301 College Ave, Fort Lauderdale, FL 33314, USA In memory of Robert C. Thompson We give a short proof of a recent result of Drury on the positivity of a 3 3 matrix of the form ( R i R j tr 1 i,j 3 for any rectangular complex (or real matrices R 1,R 2,R 3 so that the multiplication R i R j is compatible for all i,j, where tr denotes the trace norm. We then give a complete analysis of the problem when the trace norm is replaced by other unitarily invariant norms. AMS Classification: Primary: 15A60, 15B47; Secondary: 47A30, 47B65. Keywords: polar decomposition, positive semi-definite matrix, trace norm, unitarily invariant norm 1 Introduction Let A = (A ij 1 i,j m be a partitioned positive semi-definite matrix. We are interested in the question when ( A ij is also positive semi-definite, 1 i,j m where is a unitarily invariant norm; see [5, 7] and the references therein. Let R 1,R 2,R 3 be rectangular complex matrices so that the multiplication Rj R k is compatible for all (j,k pairs. Drury [1] recently showed that the 3 3 matrix M = (m jk is positive semi-definite if m jk = tr( R j R k = R j R k tr, where for a (rectangular matrix X, X denotes its conjugate transpose, X = (X X 1/2 denotes its polar factor, and X tr = tr X denotes its trace norm, E-mails: ckli@math.wm.edu E-mails: zhang@nova.edu 1

i.e. the sum of all singular values s j (X of X. We will denote by M n the set of n n complex matrices. Drury s theorem is a generalization of a result of Marcus and Watkins (see [5, p.238] asserting that if A = (a ij M p is positive semi-definite, 1 p 3, then so is ( a ij (by taking absolute value of each entry. It is known that for p 4, ( a ij need not be positive semi-definite in general. In this note, we present a short proof of Drury s result. We then give a complete analysis of the problem when the trace norm is replaced by other unitarily invariant norms. 2 Results and proofs Theorem 1 Suppose that A = (A ij 1 i,j 3 is apartitioned positive semi-definite matrix, where each A ii is a square matrix, i = 1,2,3. Then (tr A ij M 3 is also positive semi-definite. Proof. We may assume that all A ij are of size n n by adding zero rows (at the bottom and zero columns (on the right. Let A 12 have polar decomposition P 12 U. We may replace A by (I n U I n A(I n U I n and assume that A 11,A 22,A 33,A 12 = A 21 are positive semi-definite. Next, in the modified matrix, assume that A 23 has polar decomposition P 23 V. We may replace A by (I n I n VA(I n I n V and assume that A 11,A 22,A 33,A 12 = A 21,A 23 = A 32 are all positive semi-definite. Now, in the resulting matrix, suppose that A 13 has polar decomposition P 13 W. Let W = X DX, D = diag(µ 1,...,µ n for some unitary X and complex units µ 1,...,µ n. Then for X = X X X, we have P 11 P12 P13 XA X D = (Ãij = P 12 P22 P23, D P13 P23 P33 where P ij = XP ij X. If we remove all the off-diagonal entries of P ij to get a diagonal matrix Q ij, then tr Q ij = tr P ij = tr Ãij = tr A ij, 1 i,j 3. We claim that the matrix (Q ij 1 i,j 3 is positive semi-definite. It will then follow that (tr A ij = (tr Q ij is positive semi-definite; see, e.g. [4, 7]. To prove our claim, note that if we take the (r,r entries of à ij to form a matrix Q r M 3, then Q r is a principal submatrix of (Ãij and is positive semidefinite for each r = 1,...,n. By the result of 3 3 matrices, we can change the (1,3,(3,1 entries to their absolute values to get a positive semi-definite Q r M 3. Because (Q ij 1 i,j 3 is permutationally similar to Q 1 Q n and thus is positive semi-definite, we get the desired conclusion. 2

The conclusion of Theorem 1 may not hold if we replace the trace norm by other norms on matrices (that can be defined on A ij for all 1 i,j n. For example, suppose p is the Schatten p-norm defined by X p = (tr X p 1/p for p [1, ], where X is the operator norm of X and X 1 is the trace ( norm. Take the 4 4 positive semi-definite matrix A = I2 I 2 I 2 into (A ij 1 i,j 3 with A 11 = A 22 = (1 and A 33 = I 2. Then 1 0 1 ( A ij p = 0 1 1 1 1 2 1/p I 2 and partition it is not positive semi-definite if p > 1. So, it is interesting to determine the types of norms on matrices such that Theorem 1 is valid. In the following, we give a complete answer of the problem for unitarily invariant norms, i.e. norms on matrices such that UAV = A for any matrix A and any unitary matrices U and V of appropriate sizes. One may see [2] and its references for some general background of unitarily invariant norms. Suppose that A = (A ij 1 i,j m is positive semi-definite. Through block permutation, we may assume that the diagonal blocks A 11,...,A mm have sizes in ascending order. Let m m be the size of A mm. Suppose is a unitarily invariant norm on M m. Extend the definition of to other blocks by setting A ij = Ãij, where Ãij M m is obtained by adding zero rows and columns to A ij. We can then consider ( A ij for any A = (A ij. Ifm = 2, then ( A ij is positivesemi-definite; see[3]. Ifm 4, then wecan choose a positive semi-definite matrix B = (b ij M m such that the matrix B 0 obtained by taking the absolute values of the entries of B is not positive semidefinite; see the example of Thompson in [5]. Let A = (A ij such that A ij has (1,1 entry equal to b ij and all other entries equal to zero. Then ( A ij = γb 0 is not positive semi-definite, where γ = E 11 and {E 11,E 12,...,E nn } is the standard basis for M n. For m = 3, we have the following. Theorem 2 Consider the set M(n 1,n 2,n of block matrices in the form A = (A ij 1 i,j 3, where A 11 M n1, A 22 M n2, and A 33 M n with n 1 n 2 n. Let be a unitarily invariant norm on M n and k = min{n 1 +n 2,n}. The following conditions are equivalent. (a The matrix ( A ij is positive semi-definite whenever (A ij M(n 1,n 2,n is positive semi-definite. (b E 11 + +E kk = E 11 + + E kk = k E 11. Proof. We may normalize so that E 11 = 1. By the result in [6] (see also [2], there is a compact set S of real vectors v = (v 1,...,v n with v 1 3

v n 0 such that for every matrix B M n, B = max{ B v : v S} with B v = n v j s j (B, where s 1 (B s n (B are the singular values of B. Because E 11 = 1, we see that 1 = max{v 1 : (v 1,...,v n S}. Suppose (b holds. Then S contains a vector ˆv whose first k components are equal to 1 so that E 11 + +E kk = E 11 + +E kk ˆv = k. Hence, for any matrix B M n with rank not larger than k, we have s j (B = B ˆv B s j (BE jj = s j (B. That is, if rank(b k, then B tr = B. For B with rank larger than k, s j (B B ˆv B. Let A = (A ij 1 i,j 3 M(n 1,n 2,n be positive semi-definite. If k = n, then A ij has rank at most k and A ij = A ij tr for all 1 i,j 3. By Theorem 1, ( A ij is positive semi-definite. Suppose n 1 +n 2 k < n. Then rank( A13 A 23 ( k, and there is a unitary V M n such that the last n k columns of A13 A 23 V are zero. We may replace A by à = U AU without changing ( A ij, where ( U = I n1+n 2 V, and assume that the last n k columns of A13 A 23 V are zero. Suppose  = (Âij is obtained by deleting the last n k rows and columns of A. Then rank(âij k for every Âij. Because Â33 M k is a principal submatrix of A 33, from the above discussion on B, we have For (i,j (3,3, we have tr Â33 s j (A 33 A 33. tr Âij = Âij = A ij. Thus, ( A ij = (tr Âij +diag(0,0,ε, 4

where ε = A 33 tr Â33 0. By Theorem 1, (tr Âij is positive semidefinite. It follows that ( A ij is positive semi-definite. Suppose (b does not hold. Let s be the largest integer such that E 11 + + E ss = s E 11. Then s < k. Choose ñ 1 n 1,ñ 2 n 2 such that ñ 1 + ñ 2 = s + 1 = ñ n. Thus, ñ 1,ñ 2 s. By the choice of s, Iñ1 = ñ 1, Iñ2 = ñ 2, and s Iñ < s + 1. Let Iñ = s + δ for some δ [0,1. Construct the block matrix in which à = (Ãij 1 i,j 3 Mñ1+ñ 2+ñ à 11 = Iñ1 à 22 = Iñ2, à 33 = Iñ, à t 21 = Ã12 = Oñ1 ñ 2, à t 31 = Ã13 = ( Iñ1 Oñ1 ñ 2, à t 32 = Ã23 = ( Oñ2 ñ 1 Iñ2. Then à is positive semi-definite and Observe that 0 ñ 1  = ( A ij = ñ1 0 ñ 2 ñ 2. ñ 2 s+δ ñ 1 det(â = ñ 1ñ 2 (ñ 1 +ñ 2 1+δ ñ 1 ñ 2 2 ñ2 1ñ2 = ñ 1 ñ 2 (δ 1 < 0. Hence,  is not positive semi-definite. If needed, we can add zero rows and columns to à to get a matrix (A ij with A 11 M n1,a 22 M n2,a 33 M n. The matrix (A ij is positive semi-definite while  remains the same. Acknowledgment Li was an affiliate member of the Institute for Quantum Computing, University of Waterloo, an honorary professor of the Shanghai University, and University of Hong Kong. The research of Li was supported by the USA NSF and HK RCG. This project was done while Li and Zhang were participating the 2014 Summer International Program at the Shanghai University. Li and Zhang would like to thank the support of the Shanghai University. 5

References [1] S.W. Drury, Positive semidefiniteness of a matrix related to partitioning, Linear Algebra Appl. 446 (2014 369-376. [2] C.K. Li, Some aspects of the theory of norms, Linear Algebra Appl. 212 (1994 71-100. [3] C.K. Li and R. Mathias, Generalizations of Ky Fan s dominance theorem, SIAM J. Matrix Anal. Appl. 19 (1998 99-106. [4] M. Marcus and S.M. Katz, Matrices of Schur functions, Duke Math. J. 36 (1969 343-352. [5] M. Marcus and W. Watkins, Partitioned Hermitian matrices, Duke Math. J. 38 (1971 237-249. [6] R. Mathias, A Representation of Unitary Similarity Invariant Norms, Technical Report 515, Dept. of Mathematical Sciences, Johns Hopkins Univ., Blatimore, 1989. [7] F. Zhang, Positivity of matrices with generalized matrix functions, Acta Mathematica Sinica, Vol. 28, No. 9, pp. 1779-1786, 2012. 6