arxiv: v1 [math.ca] 22 Jun 2016

Similar documents
arxiv:math/ v5 [math.ac] 17 Sep 2009

Remarks on BMV conjecture

Herbert Stahl s proof of the BMV conjecture

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

arxiv: v2 [math.fa] 17 May 2016

A determinant characterization of moment sequences with finitely many mass-points

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

arxiv: v1 [math.fa] 1 Oct 2015

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

Banach Journal of Mathematical Analysis ISSN: (electronic)

Fibonacci numbers, Euler s 2-periodic continued fractions and moment sequences

The matrix arithmetic-geometric mean inequality revisited

Spectral inequalities and equalities involving products of matrices

Iteration of the rational function z 1/z and a Hausdorff moment sequence

Lecture notes on Quantum Computing. Chapter 1 Mathematical Background

Recall the convention that, for us, all vectors are column vectors.

A Note on Product Range of 3-by-3 Normal Matrices

MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.

Log-Convexity Properties of Schur Functions and Generalized Hypergeometric Functions of Matrix Argument. Donald St. P. Richards.

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

Math 408 Advanced Linear Algebra

Math 554 Qualifying Exam. You may use any theorems from the textbook. Any other claims must be proved in details.

Scattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions

Elementary linear algebra

Two Results About The Matrix Exponential

arxiv: v3 [math.ra] 22 Aug 2014

arxiv: v1 [math.dg] 29 Dec 2018

On the closures of orbits of fourth order matrix pencils

Some applications of matrix inequalities in Rényi entropy

JUST THE MATHS UNIT NUMBER 1.3. ALGEBRA 3 (Indices and radicals (or surds)) A.J.Hobson

arxiv: v4 [math.ca] 9 May 2012

Lecture notes: Applied linear algebra Part 1. Version 2

arxiv: v1 [math.mg] 15 Jul 2013

LINEAR MAPS PRESERVING NUMERICAL RADIUS OF TENSOR PRODUCTS OF MATRICES

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

Singular Value Inequalities for Real and Imaginary Parts of Matrices

Moore Penrose inverses and commuting elements of C -algebras

arxiv: v1 [math.ca] 7 Jan 2015

Asymptotic behaviour of Hilbert space operators with applications

Journal of Inequalities in Pure and Applied Mathematics

Clarkson Inequalities With Several Operators

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

A Questionable Distance-Regular Graph

arxiv: v3 [math.ra] 10 Jun 2016

Notes on matrix arithmetic geometric mean inequalities

Spectral law of the sum of random matrices

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

Rank minimization via the γ 2 norm

Memoryless output feedback nullification and canonical forms, for time varying systems

Eigenvectors and Reconstruction

Diagonal and Monomial Solutions of the Matrix Equation AXB = C

Journal of Algebra 226, (2000) doi: /jabr , available online at on. Artin Level Modules.

On the adjacency matrix of a block graph

On lower and upper bounds for probabilities of unions and the Borel Cantelli lemma

1 Linear Algebra Problems

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

Quantum Computing Lecture 2. Review of Linear Algebra

Compound matrices and some classical inequalities

The Singular Value Decomposition

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

Eigenvalues and Eigenvectors: An Introduction

HIGHER RANK NUMERICAL RANGES OF NORMAL MATRICES

Rank of Matrices and the Pexider Equation

Orthogonality graph of the algebra of upper triangular matrices *

RANKS OF QUANTUM STATES WITH PRESCRIBED REDUCED STATES

Hadamard matrices and Compact Quantum Groups

Irina Pchelintseva A FIRST-ORDER SPECTRAL PHASE TRANSITION IN A CLASS OF PERIODICALLY MODULATED HERMITIAN JACOBI MATRICES

Block Diagonal Majorization on C 0

Transformations Preserving the Hankel Transform

THE NUMBER OF ORTHOGONAL CONJUGATIONS

Knowledge Discovery and Data Mining 1 (VO) ( )

Matrix Lie groups. and their Lie algebras. Mahmood Alaghmandan. A project in fulfillment of the requirement for the Lie algebra course

ON THE UNIQUENESS PROPERTY FOR PRODUCTS OF SYMMETRIC INVARIANT PROBABILITY MEASURES

Homework 2 Foundations of Computational Math 2 Spring 2019

RANDOM WALKS WHOSE CONCAVE MAJORANTS OFTEN HAVE FEW FACES

Gaps between classes of matrix monotone functions

arxiv: v1 [math.ca] 6 May 2016

arxiv: v1 [physics.gen-ph] 18 Dec 2014

Exponentials of Symmetric Matrices through Tridiagonal Reductions

MAGNETIC BLOCH FUNCTIONS AND VECTOR BUNDLES. TYPICAL DISPERSION LAWS AND THEIR QUANTUM NUMBERS

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

The Singular Value Decomposition (SVD) and Principal Component Analysis (PCA)

Lecture 1 Review: Linear models have the form (in matrix notation) Y = Xβ + ε,

A brief introduction to p-adic numbers

APPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.

TRACIAL POSITIVE LINEAR MAPS OF C*-ALGEBRAS

Short proofs of theorems of Mirsky and Horn on diagonals and eigenvalues of matrices

Quantum Systems Measurement through Product Hamiltonians

AN EXTENDED MATRIX EXPONENTIAL FORMULA. 1. Introduction

Normal Toeplitz Matrices

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

An Approximation Algorithm for Approximation Rank

Norm inequalities related to the matrix geometric mean

Mathematical Methods for Engineers and Scientists 1

PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY MOORE PENROSE INVERSE OF BIDIAGONAL MATRICES. I

Elements of Positive Definite Kernel and Reproducing Kernel Hilbert Space

Angular matrix integrals

On the Hermitian solutions of the

CHAPTER III THE PROOF OF INEQUALITIES

Some bounds for the spectral radius of the Hadamard product of matrices

Transcription:

arxiv:1606.06911v1 [math.ca] 22 Jun 2016 On a special case of the Herbert Stahl theorem Victor Katsnelson Abstract. The BMV conjecture states that for n n Hermitian matrices A and B the function f A,B(t) = tracee ta+b is exponentially convex. Recently the BMV conjecture was proved by Herbert Stahl. The proof of Herbert Stahl is based on ingenious considerations related to Riemann surfaces of algebraic functions. In the present paper we give a purely "matrix" proof of the BMV conjecture for the special case rank A = 1. This proof is based on the Lie product formula for the exponential of the sum of two matrices and does not require complex analysis. Mathematics Subject Classification (2010). Primary 15A15,15A16; Secondary 30F10,44A10. Keywords. BMV conjecture, absolutely monotonic functions, exponentially convex functions, positive definite functions, Lie product formula. 1. Herbert Stahl s Theorem. In the paper [1] a conjecture was formulated which now is commonly known as the BMV conjecture: The BMV Conjecture. Let A and B be Hermitian matrices of size n n. Then the function f A,B (t) = trace{exp[ta+b]} (1.1) of the variable t is representable as a bilateral Laplace transform of a nonnegative measure dσ A,B (λ) compactly supported on the real axis: f A,B (t) = exp(tλ)dσ A,B (λ), t (, ). (1.2) λ (, ) Let us note that the function f A,B (t), considered for t C, is an entire function of exponential type. The indicator diagram of the function f A,B is the closed interval [λ min,λ max ], where λ min and λ max are the least and the greatest eigenvalues of the matrix A respectively. Thus if the function f A,B (t) is representable in the form (1.2) with a non-negative measuredσ A,B (λ), then

2 V. Katsnelson dσ A,B (λ) is actually supported on the interval [λ min,λ max ] and the representation f A,B (t) = exp(tλ) dσ A,B (λ), t C, (1.3) λ [λ min,λ max] holds for every t C. The representability of the function f A,B (t), (1.1), in the form (1.3) with a non-negative dσ A,B is evident if the matrices A and B commute. In this case dσ(λ) is an atomic measure supported on the spectrum of the matrix A. In general case, if the matrices A and B do not commute, the BMV conjecture remained an open question for longer than 35 years. In 2011, Herbert Stahl proved the BMV conjecture. Theorem(H.Stahl) Let A and B be n n hermitian matrices. Then the function f A,B (t) defined by (1.1) is representable as the bilateral Laplace transform (1.3) of a non-negative measure dσ A,B (λ) supported on the closed interval [λ min,λ max ]. The first arxiv version of H.Stahl s Theorem appeared in [2], the latest arxiv version - in [3], the journal publication - in [4]. The proof of Herbert Stahl is based on ingenious considerations related to Riemann surfaces of algebraic functions. In [5],[6] a simplified version of the Herbert Stahl proof is presented. The proof, presented in [5],[6], preserves all the main ideas of Stahl; the simplification consists in technical details. In the present paper we present a proof the BMV conjecture for the special case ranka = 1. Our proof is based on an elementary argument which does not require complex analysis. 2. Exponentially convex functions. Definition 2.1. A function f on R, f : R [0, ), is said to be exponentially convex if 1. For every nonnegative integer N, for every choice of real numbers t 1, t 2,..., t N, and complex numbers ξ 1, ξ 2,...,ξ N, the inequality holds N f(t r +t s )ξ r ξ s 0; (2.1) r,s=1 2. The function f is continuous on R. The class of exponentially convex functions was introduced by S.N.Bernstein, [7], see 15 there. Russian translation of the paper [7] can be found in [8, pp. 370 425]. From (2.1) it follows that the inequality f(t 1 + t 2 ) f(2t 1 )f(2t 2 ) holds for every t 1 R,t 2 R. Thus the alternative takes place: If f is an exponentially convex function, then either f(t) 0, or f(t) > 0 for every t R.

On a special case of the Herbert Stahl theorem 3 Properties of the class of exponentially convex functions. P1. If f(t) is an exponentially convex function and c 0 is a nonnegative constant, then the function cf(t) is exponentially convex. P2. If f 1 (t) and f 2 (t) are exponentially convex functions, then their sum f 1 (t)+f 2 (t) is exponentially convex. P3. If f 1 (t) and f 2 (t) are exponentially convex functions, then their product f 1 (t) f 2 (t) is exponentially convex. P4. Let {f n (t)} 1 n< be a sequence of exponentially convex functions. We assume that for each t R there exists the limit f(t) = lim n f n (t), and that f(t) < t R. Then the limiting function f(t) is exponentially convex. From the functional equation for the exponential function it follows that for each real number µ, for every choice of real numbers t 1,t 2,..., t N and complex numbers ξ 1, ξ 2,...,ξ N, the equality holds N N 2 e (tr+ts)µ ξ r ξ s = e tpµ ξ p 0. (2.2) r,s=1 The relation (2.2) can be formulated as Lemma 2.2. For each real number µ, the function e tµ of the variable t is exponentially convex. The following result is well known. p=1 Theorem 2.3 (The representation theorem). 1. Let σ(dµ) be a nonnegative measure on the real axis, and let the function f(t) be defined as the two-sided Laplace transform of the measure σ(dµ): f(t) = e tµ σ(dµ), (2.3) µ R where the integral in the right hand side of (2.3) is finite for any t R. Then the function f is exponentially convex. 2. Let f(t) be an exponentially convex function. Then this function f can be represented on R as a two-sided Laplace transform (2.3) of a nonnegative measure σ(dµ). (In particular, the integral in the right hand side of (2.3) is finite for any t R.) The representing measure σ(dµ) is unique. The assertion 1 of the representation theorem is an evident consequence of Lemma 2.2, of the properties P1, P2, P4, and of the definition of the integration operation. The proof of the assertion 2 can be found in [9], Theorem 5.5.4, and in [10], Theorem 21. Thus the Herbert Stahl theorem can be reformulated as follows: Let A and B be Hermitian n n matrices. Let the function f A,B (t) is defined by (1.1) for t (, ). Then the function f A,B (t), considered as a function of the variable t, is exponentially convex.

4 V. Katsnelson 3. A special case of the Herbert Stahl theorem Lemma 3.1. Let M be a Hermitian matrix. Assume that all off-diagonal entries of the matrix M are non-negative. Then all entries of the matrix exponential e M are non-negative. Proof. Since the matrix M is Hermitian, its diagonal entries are real. If a def positive number ρ is large enough, all entries of the matrix M ρ = M + ρi, where I is the identity matrix, are non-negative. We choose and fix such ρ. All entries of the matrix e Mρ are nonnegative. Moreover e M = e ρ emρ. Lemma 3.2. Let L and M be Hermitian matrices of the same size, say n n. We assume that 1. The matrix L is diagonal; 2. All off-diagonal entries of the matrix M are non-negative. Then each entry of the matrix function e Lt+M is an exponentially convex function of the variable t. In particular the function f L,M (t) = trace [ e Lt+M] is exponentially convex. Proof. We use the Lie product formula 1 : e X+Y ( X Y ) = lim e p e p p, (3.1) p where X and Y are arbitrary square matrices of the same size. Taking X = Lt, Y = M, we obtain e Lt+M ( Lt M) = lim e p p. (3.2) p According to Lemma 3.1, all entries of the matrixe M/p are non-negative numbers. Since the matrix L is Hermitian, its diagonal entries are real numbers. Therefore, the matrix function e Lt/p is of form e Lt/p = diag ( e l1t/p,...,e lnt/p), where l 1,...,l n are real numbers. The exponentials e ljt/p are exponentially convex functions of t. Each entry of the matrix e Lt M p e p is a linear combination of these exponentials with non-negative coefficients. According to the properties P1 and P2, the entries of the matrix function e Lt M p are exponentially convex functions. Each entry of the matrix function ( e Lt M) p p is a sum of products of some entries of the matrix function e Lt M p. According to the properties P2 and P3, the entries of the matrix function ( e Lt M) p p are exponentially convex functions. From the limiting relation (3.2) and from the property P4 it follows that all entries of the matrix e Lt+M are exponentially convex functions. All the more, the function f L,M (t), which is the sum of the diagonal entries, is exponentially convex. 1 See [11, Theorem 2.10].

On a special case of the Herbert Stahl theorem 5 Theorem 3.3. Let A and B be Hermitian matrices of size n n. Assume moreover that the matrix A is of rank one. Then the function f A,B (t), defined by (1.1), is exponentially convex. Proof. 1. IfU is an unitary matrix, then f A,B (t) = f UAU,UBU (t). Since ranka = 1, we can choose the matrix U such that the matrix UAU are of the form UAU = diag(λ 1,...,λ n 1,λ n ), where λ j = 0, 1 j n 1, λ n 0. We choose and fix such unitary matrix U. 2. Thus from the very beginning we can assume that the matrices A and B are of the form [ ] [ ] 0n 1 0 A = n 1 Bn 1 b 0, B = n 1 n 1 λ n b, n 1 µ n where 0 n 1 is the zero matrix of size (n 1) (n 1), 0 n 1 is the zero column of size n 1, λ n 0 is a real number, B n 1 is the a Hermitian matrix of size (n 1) (n 1), b n 1 is a column of size n 1, µ n is a real number. There exists an unitary matrix V n 1 of size (n 1) (n 1) such that the matrix V n 1 B n 1 V n 1 is diagonal. We choose and fix such matrix V n 1. We define the matrices M n 1 def = V n 1 B n 1 V n 1, g def = V n 1 b n 1. The matrix M n 1 is a diagonal matrix of size (n 1) (n 1): M n 1 = diag(µ 1,...,µ n 1 ). The matrix g is a column of size n 1: g = γ 1. γ n 1 where γ j, 1 j n 1, are complex numbers. Let us define numbers ω j which satisfy the conditions, ω j = 1, ω j γ j = γ j, 1 j n 1. If γ j 0 for some j, then such ω j is unique. If γ j = 0 for some j, we take ω j = 1. Let us define the diagonal matrix Ω n 1 of size (n 1) (n 1) as Ω n 1 = diag(ω 1,...,,ω n 1 ). According to the construction of the matrix Ω n 1, all entries of the column are non-negative numbers: g def = Ω n 1 g g = γ 1.. γ n 1.

6 V. Katsnelson Moreover Ω n 1 M n 1 Ω n 1 = M n 1. We introduce the matrices W n 1 = Ω n 1 V n 1, [ ] Wn 1 0 W = n 1 0 n 1 1 The equalities hold, where L = WAW = L, [ ] 0n 1 0 n 1 0 n 1 λ 3. The matrix W is unitary. Therefore WBW = M, M = f A,B (t) = f L,M (t). [ ] Mn 1 g g µ n The matrix L is diagonal. (Actually L = A.) Off-diagonal entries of the matrix M are non-negative numbers: m j,k = 0, 1 j < k < n; m j,n = γ j, 1 j n 1. According to Lemma 3.2, the function f L,M (t) is exponentially convex. References [1] D. Bessis, P. Moussa, M. Villani. Monotonic converging variational approximations to the functional integrals in quantum statistical mechanics. J. Mat. Phys., 16:11 (1975), 2318-2325. [2] H. Stahl. Proof of the BMV conjecture. arxiv:1107.4875v1, 1-56, 25 Jul2011. [3] H. Stahl. Proof of the BMV conjecture. arxiv:1107.4875v3, 1-25, 17 Aug2012. [4] H.Stahl. Proof of the BMV conjecture. Acta Math., 211 (2013), 255-290. [5] A. Eremenko. Herbert Stahl s proof of the BMV conjecture. arxiv:1312.6003. [6] А.Э.Еременко. Доказательство Герберта Шталя гипотезы БМВ. Матем. Сб., 206:1 (2015), 97-102. (In Russian). English translation: A.Eremenko. Herbert Stahl s proof of the BMV conjecture. Sbornik: Mathematics, 206:1 (2015), 87-92. 211 (2013), 255-290. [7] S.N. Bernstein. Sur les functions absolument monotones. Acta Math. 52 (1928), 1-66. (In French). [8] С.Н. Бернштейн. Собрание Сочинений. Том 1. Издательство АН СССР, 1952. [9] Н.И.Ахиезер. Классическая Проблема Моментов. Физматгиз, Москва, 1965. (In Russian). English Transl.: N.I. Akhiezer. The Clasical Moment Problem. Oliver and Boyd, Edinburgh and London, 1965. [10] D.V. Widder. Laplace Transform. Princeton Univ. Press, Princeton N.J., 1946. [11] B.C. Hall. Lie Groups, Lie Algebras and Representations. Springer-Verlag, New York, 2003.

On a special case of the Herbert Stahl theorem 7 Victor Katsnelson Department of Mathematics The Weizmann Institute 76100, Rehovot Israel e-mail: victor.katsnelson@weizmann.ac.il; victorkatsnelson@gmail.com