Trace inequality

From testwiki
Revision as of 02:26, 27 December 2024 by imported>LR.127 (Adding local short description: "Concept in Hlibert spaces mathematics", overriding Wikidata description "inequalities involving linear operators on Hilbert spaces")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Template:Short description In mathematics, there are many kinds of inequalities involving matrices and linear operators on Hilbert spaces. This article covers some important operator inequalities connected with traces of matrices.[1][2][3][4]

Basic definitions

Let 𝐇n denote the space of Hermitian n×n matrices, 𝐇n+ denote the set consisting of positive semi-definite n×n Hermitian matrices and 𝐇n++ denote the set of positive definite Hermitian matrices. For operators on an infinite dimensional Hilbert space we require that they be trace class and self-adjoint, in which case similar definitions apply, but we discuss only matrices, for simplicity.

For any real-valued function f on an interval I, one may define a matrix function f(A) for any operator A𝐇n with eigenvalues λ in I by defining it on the eigenvalues and corresponding projectors P as f(A)jf(λj)Pj, given the spectral decomposition A=jλjPj.

Operator monotone

Template:Main

A function f:I defined on an interval I is said to be operator monotone if for all n, and all A,B𝐇n with eigenvalues in I, the following holds, ABf(A)f(B), where the inequality AB means that the operator AB0 is positive semi-definite. One may check that f(A)=A2 is, in fact, not operator monotone!

Operator convex

A function f:I is said to be operator convex if for all n and all A,B𝐇n with eigenvalues in I, and 0<λ<1, the following holds f(λA+(1λ)B)λf(A)+(1λ)f(B). Note that the operator λA+(1λ)B has eigenvalues in I, since A and B have eigenvalues in I.

A function f is Template:Visible anchor if f is operator convex;=, that is, the inequality above for f is reversed.

Template:AnchorJoint convexity

A function g:I×J, defined on intervals I,J is said to be Template:Visible anchor if for all n and all A1,A2𝐇n with eigenvalues in I and all B1,B2𝐇n with eigenvalues in J, and any 0λ1 the following holds g(λA1+(1λ)A2,λB1+(1λ)B2)λg(A1,B1)+(1λ)g(A2,B2).

A function g is Template:Visible anchor if −g is jointly convex, i.e. the inequality above for g is reversed.

Trace function

Given a function f:, the associated trace function on 𝐇n is given by ATrf(A)=jf(λj), where A has eigenvalues λ and Tr stands for a trace of the operator.

Convexity and monotonicity of the trace function

Let f: be continuous, and let Template:Mvar be any integer. Then, if tf(t) is monotone increasing, so is ATrf(A) on Hn.

Likewise, if tf(t) is convex, so is ATrf(A) on Hn, and it is strictly convex if Template:Mvar is strictly convex.

See proof and discussion in,[1] for example.

Löwner–Heinz theorem

For 1p0, the function f(t)=tp is operator monotone and operator concave.

For 0p1, the function f(t)=tp is operator monotone and operator concave.

For 1p2, the function f(t)=tp is operator convex. Furthermore,

f(t)=log(t) is operator concave and operator monotone, while
f(t)=tlog(t) is operator convex.

The original proof of this theorem is due to K. Löwner who gave a necessary and sufficient condition for Template:Mvar to be operator monotone.[5] An elementary proof of the theorem is discussed in [1] and a more general version of it in.[6]

Template:AnchorKlein's inequality

For all Hermitian Template:Mvar×Template:Mvar matrices Template:Mvar and Template:Mvar and all differentiable convex functions f: with derivative Template:Math, or for all positive-definite Hermitian Template:Mvar×Template:Mvar matrices Template:Mvar and Template:Mvar, and all differentiable convex functions Template:Mvar:(0,∞) → , the following inequality holds,

Tr[f(A)f(B)(AB)f(B)]0.

In either case, if Template:Mvar is strictly convex, equality holds if and only if Template:Mvar = Template:Mvar. A popular choice in applications is Template:Math, see below.

Proof

Let C=AB so that, for t(0,1),

B+tC=(1t)B+tA,

varies from B to A.

Define

F(t)=Tr[f(B+tC)].

By convexity and monotonicity of trace functions, F(t) is convex, and so for all t(0,1),

F(0)+t(F(1)F(0))F(t),

which is,

F(1)F(0)F(t)F(0)t,

and, in fact, the right hand side is monotone decreasing in t.

Taking the limit t0 yields,

F(1)F(0)F(0),

which with rearrangement and substitution is Klein's inequality:

tr[f(A)f(B)(AB)f(B)]0

Note that if f(t) is strictly convex and C0, then F(t) is strictly convex. The final assertion follows from this and the fact that F(t)F(0)t is monotone decreasing in t.

Golden–Thompson inequality

Template:Main

In 1965, S. Golden [7] and C.J. Thompson [8] independently discovered that

For any matrices A,B𝐇n,

TreA+BTreAeB.

This inequality can be generalized for three operators:[9] for non-negative operators A,B,C𝐇n+,

TrelnAlnB+lnC0TrA(B+t)1C(B+t)1dt.

Peierls–Bogoliubov inequality

Let R,F𝐇n be such that Tr eR = 1. Defining Template:Math, we have

TreFeRTreF+Reg.

The proof of this inequality follows from the above combined with Klein's inequality. Take Template:Math.[10]

Gibbs variational principle

Let H be a self-adjoint operator such that eH is trace class. Then for any γ0 with Trγ=1,

TrγH+TrγlnγlnTreH,

with equality if and only if γ=exp(H)/Trexp(H).

Lieb's concavity theorem

The following theorem was proved by E. H. Lieb in.[9] It proves and generalizes a conjecture of E. P. Wigner, M. M. Yanase, and Freeman Dyson.[11] Six years later other proofs were given by T. Ando [12] and B. Simon,[3] and several more have been given since then.

For all m×n matrices K, and all q and r such that 0q1 and 0r1, with q+r1 the real valued map on 𝐇m+×𝐇n+ given by

F(A,B,K)=Tr(K*AqKBr)
  • is jointly concave in (A,B)
  • is convex in K.

Here K* stands for the adjoint operator of K.

Lieb's theorem

For a fixed Hermitian matrix L𝐇n, the function

f(A)=Trexp{L+lnA}

is concave on 𝐇n++.

The theorem and proof are due to E. H. Lieb,[9] Thm 6, where he obtains this theorem as a corollary of Lieb's concavity Theorem. The most direct proof is due to H. Epstein;[13] see M.B. Ruskai papers,[14][15] for a review of this argument.

Ando's convexity theorem

T. Ando's proof [12] of Lieb's concavity theorem led to the following significant complement to it:

For all m×n matrices K, and all 1q2 and 0r1 with qr1, the real valued map on 𝐇m++×𝐇n++ given by

(A,B)Tr(K*AqKBr)

is convex.

Template:AnchorJoint convexity of relative entropy

For two operators A,B𝐇n++ define the following map

R(AB):=Tr(AlogA)Tr(AlogB).

For density matrices ρ and σ, the map R(ρσ)=S(ρσ) is the Umegaki's quantum relative entropy.

Note that the non-negativity of R(AB) follows from Klein's inequality with f(t)=tlogt.

Statement

The map R(AB):𝐇n++×𝐇n++𝐑 is jointly convex.

Proof

For all 0<p<1, (A,B)Tr(B1pAp) is jointly concave, by Lieb's concavity theorem, and thus

(A,B)1p1(Tr(B1pAp)TrA)

is convex. But

limp11p1(Tr(B1pAp)TrA)=R(AB),

and convexity is preserved in the limit.

The proof is due to G. Lindblad.[16]

Jensen's operator and trace inequalities

The operator version of Jensen's inequality is due to C. Davis.[17]

A continuous, real function f on an interval I satisfies Jensen's Operator Inequality if the following holds

f(kAk*XkAk)kAk*f(Xk)Ak,

for operators {Ak}k with kAk*Ak=1 and for self-adjoint operators {Xk}k with spectrum on I.

See,[17][18] for the proof of the following two theorems.

Jensen's trace inequality

Let Template:Mvar be a continuous function defined on an interval Template:Mvar and let Template:Mvar and Template:Mvar be natural numbers. If Template:Mvar is convex, we then have the inequality

Tr(f(k=1nAk*XkAk))Tr(k=1nAk*f(Xk)Ak),

for all (Template:Mvar1, ... , Template:Mvarn) self-adjoint Template:Mvar × Template:Mvar matrices with spectra contained in Template:Mvar and all (Template:Mvar1, ... , Template:Mvarn) of Template:Mvar × Template:Mvar matrices with

k=1nAk*Ak=1.

Conversely, if the above inequality is satisfied for some Template:Mvar and Template:Mvar, where Template:Mvar > 1, then Template:Mvar is convex.

Jensen's operator inequality

For a continuous function f defined on an interval I the following conditions are equivalent:

  • f is operator convex.
  • For each natural number n we have the inequality
f(k=1nAk*XkAk)k=1nAk*f(Xk)Ak,

for all (X1,,Xn) bounded, self-adjoint operators on an arbitrary Hilbert space with spectra contained in I and all (A1,,An) on with k=1nAk*Ak=1.

  • f(V*XV)V*f(X)V for each isometry V on an infinite-dimensional Hilbert space and

every self-adjoint operator X with spectrum in I.

  • Pf(PXP+λ(1P))PPf(X)P for each projection P on an infinite-dimensional Hilbert space , every self-adjoint operator X with spectrum in I and every λ in I.

Araki–Lieb–Thirring inequality

Template:Distinguish

E. H. Lieb and W. E. Thirring proved the following inequality in [19] 1976: For any A0, B0 and r1, Tr((BAB)r)Tr(BrArBr).

In 1990 [20] H. Araki generalized the above inequality to the following one: For any A0, B0 and q0, Tr((BAB)rq)Tr((BrArBr)q), for r1, and Tr((BrArBr)q)Tr((BAB)rq), for 0r1.

There are several other inequalities close to the Lieb–Thirring inequality, such as the following:[21] for any A0, B0 and α[0,1], Tr(BAαBBA1αB)Tr(B2AB2), and even more generally:[22] for any A0, B0, r1/2 and c0, Tr((BAB2cAB)r)Tr((Bc+1A2Bc+1)r). The above inequality generalizes the previous one, as can be seen by exchanging A by B2 and B by A(1α)/2 with α=2c/(2c+2) and using the cyclicity of the trace, leading to Tr((BAαBBA1αB)r)Tr((B2AB2)r).

Additionally, building upon the Lieb-Thirring inequality the following inequality was derived: [23] For any A,B𝐇n,Tn×n and all 1p,q with 1/p+1/q=1, it holds that |Tr(TAT*B)|Tr(T*T|A|p)1pTr(TT*|B|q)1q.

Effros's theorem and its extension

E. Effros in [24] proved the following theorem.

If f(x) is an operator convex function, and L and R are commuting bounded linear operators, i.e. the commutator [L,R]=LRRL=0, the perspective

g(L,R):=f(LR1)R

is jointly convex, i.e. if L=λL1+(1λ)L2 and R=λR1+(1λ)R2 with [Li,Ri]=0 (i=1,2), 0λ1,

g(L,R)λg(L1,R1)+(1λ)g(L2,R2).

Ebadian et al. later extended the inequality to the case where L and R do not commute .[25]

Template:Visible anchor, named after its originator John von Neumann, states that for any n×n complex matrices A and B with singular values α1α2αn and β1β2βn respectively,[26] |Tr(AB)|i=1nαiβi, with equality if and only if A and B share singular vectors.[27]

A simple corollary to this is the following result:[28] For Hermitian n×n positive semi-definite complex matrices A and B where now the eigenvalues are sorted decreasingly (a1a2an and b1b2bn, respectively), i=1naibni+1Tr(AB)i=1naibi.

See also

References

Template:Reflist

  1. 1.0 1.1 1.2 E. Carlen, Trace Inequalities and Quantum Entropy: An Introductory Course, Contemp. Math. 529 (2010) 73–140 Template:Doi
  2. R. Bhatia, Matrix Analysis, Springer, (1997).
  3. 3.0 3.1 B. Simon, Trace Ideals and their Applications, Cambridge Univ. Press, (1979); Second edition. Amer. Math. Soc., Providence, RI, (2005).
  4. M. Ohya, D. Petz, Quantum Entropy and Its Use, Springer, (1993).
  5. Template:Cite journal
  6. W.F. Donoghue, Jr., Monotone Matrix Functions and Analytic Continuation, Springer, (1974).
  7. Template:Cite journal
  8. Template:Cite journal
  9. 9.0 9.1 9.2 Template:Cite journal
  10. D. Ruelle, Statistical Mechanics: Rigorous Results, World Scient. (1969).
  11. Template:Cite journal
  12. 12.0 12.1 Template:Cite journal
  13. Template:Cite journal
  14. Template:Cite journal
  15. Template:Cite journal
  16. Template:Cite journal
  17. 17.0 17.1 C. Davis, A Schwarz inequality for convex operator functions, Proc. Amer. Math. Soc. 8, 42–44, (1957).
  18. Template:Cite journal
  19. E. H. Lieb, W. E. Thirring, Inequalities for the Moments of the Eigenvalues of the Schrödinger Hamiltonian and Their Relation to Sobolev Inequalities, in Studies in Mathematical Physics, edited E. Lieb, B. Simon, and A. Wightman, Princeton University Press, 269–303 (1976).
  20. Template:Cite journal
  21. Z. Allen-Zhu, Y. Lee, L. Orecchia, Using Optimization to Obtain a Width-Independent, Parallel, Simpler, and Faster Positive SDP Solver, in ACM-SIAM Symposium on Discrete Algorithms, 1824–1831 (2016).
  22. L. Lafleche, C. Saffirio, Strong Semiclassical Limit from Hartree and Hartree-Fock to Vlasov-Poisson Equation, arXiv:2003.02926 [math-ph].
  23. V. Bosboom, M. Schlottbom, F. L. Schwenninger, On the unique solvability of radiative transfer equations with polarization, in Journal of Differential Equations, (2024).
  24. Template:Cite journal
  25. Template:Cite journal
  26. Template:Cite journal
  27. Template:Cite journal
  28. Template:Cite book