21.2  Convex Functions of a Majorization

Section 1 turned the n inequalities of a majorization into a single identity, x=Dy with D doubly stochastic. This section spends it. Chapter 17 §07 left to this chapter “the general consequences of a majorization x≺y, namely ∑iϕ(xi)≤∑iϕ(yi) for every convex ϕ:R→R, …”; that is the first theorem below, and its converse says that the convex functions see nothing but the majorization. Two more promises are paid from it: Chapter 17 §08’s claim that for nonnegative vectors “x≺y implies x1⋯xn≥y1⋯yn”, and Chapters 18 §10 and §11’s claim that “A↦trf(A)=∑if(λi(A)) is convex on the Hermitian matrices with spectrum in that interval” for every convex f.

Conventions. Vectors are real and ≺, ≺w are Definition 17.8.1. A convex function is one in the sense of Definition 18.10.1, and an interval of R is a convex subset of R, since it contains every number lying between two of its points. In the last subsection F=R or C, matrices in Mn(F) written A,B,C are Hermitian, and f(A) is the spectral functional calculus of Chapter 12 §08, Definition 12.8.2. There Pj, with an integer subscript, is an orthogonal projection of a spectral resolution; Pσ always means a permutation matrix. No analysis is imported here: the only facts used from outside linear algebra reach us through Chapter 18’s already proved results.

21.2.1. Karamata’s inequality

A majorization is a statement about sums of the largest entries. A convex function is one whose graph lies below its chords. The link between them is that an averaging matrix turns one into the other, and §01 says that a majorization is an averaging matrix.

Theorem 21.2.1 (Karamata’s Inequality)

Let I⊆R be an interval, let ϕ:I→R be convex, and let x,y∈Rn satisfy x≺y, with every entry of y lying in I. Then every entry of x lies in I, and ∑i=1nϕ(xi) ≤ ∑i=1nϕ(yi).

Idea. Write x=Dy with D doubly stochastic. Each xi is then a weighted average of the yj, with the weights in row i of D; Jensen’s inequality bounds ϕ(xi) by the same average of the ϕ(yj). Add these n inequalities. The coefficient of ϕ(yj) on the right is the sum of column j of D, which is 1. The rows and the columns do two different jobs, and this is where “doubly” earns its keep.

Proof. By Theorem 21.1.3 there is a doubly stochastic D=(dij) with x=Dy, that is, xi=∑j=1ndijyj(i=1,…,n). Fix i. The numbers di1,…,din are nonnegative and add up to 1, since row i of D is a row of a doubly stochastic matrix (Definition 18.8.4). So xi is a convex combination of y1,…,yn, all of which lie in the convex set I; by Lemma 18.1.2, xi∈I. Now Theorem 18.10.2, applied to the convex ϕ on I with the points y1,…,yn and the weights di1,…,din, gives ϕ(xi)=ϕ(∑jdijyj) ≤ ∑j=1ndijϕ(yj). Summing over i and exchanging the order of the two finite sums, ∑i=1nϕ(xi) ≤ ∑j=1n(∑i=1ndij)ϕ(yj)=∑j=1nϕ(yj), because every column of D adds up to 1. This proves the theorem.

The proof is three lines and uses Theorem 21.1.3 once. That is what §01 was for: without it we would be comparing sorted partial sums with the values of an arbitrary convex function, and there is no route from one to the other.

Example (One majorization, four inequalities)

Let y=(5,3,1,1) and x=(3,3,2,2). Check that x≺y, and write down what Theorem 21.2.1 gives for ϕ(t)=t2, for ϕ(t)=|t−2|, and for ϕ(t)=max{t−2,0}.

Solution. Both vectors are decreasing with total 10, and the running totals are 3,6,8 for x against 5,8,9 for y. So (M1) and (M2) hold, and x≺y.

All three functions are convex on I=R: t2 by Lemma 18.10.5, since its derivative 2t is increasing, and the other two as pointwise maxima of affine functions, by Proposition 18.10.3 (|t−2|=max{t−2, 2−t}). So: ∑ixi2=9+9+4+4=26 ≤ 25+9+1+1=36=∑iyi2,∑i|xi−2|=1+1+0+0=2 ≤ 3+1+1+1=6,∑imax{xi−2,0}=1+1=2 ≤ 3+1=4. Each is an instance of the one theorem. The last one measures how much of each vector sits above the level 2, and it is the test function the converse below is built from.

Quick check. Does Theorem 21.2.1 apply to ϕ(t)=t3 on I=R? Test it on x=(0,0,0)≺y=(1,1,−2).

Solution. No. On R the function t3 is not convex: at a=−2, b=0, and t=12 the midpoint value is ϕ(−1)=−1, while the chord gives 12ϕ(−2)+12ϕ(0)=−4, and −1>−4. The conclusion of the theorem does fail here. Both vectors have total 0, and the running totals 0,0 of x are at most 1,2, so x≺y. But ∑ixi3=0 while ∑iyi3=1+1−8=−6, and 0>−6. On I=(0,∞), where t3 is convex, since its derivative 3t2 is increasing there (Lemma 18.10.5), the theorem does apply, and no such failure is possible for vectors with strictly positive entries.

21.2.2. The converse, and the weak version

Karamata’s inequality is not merely a consequence of majorization; it is equivalent to it. Two families of test functions do the work: the affine ones pin down the totals, and the “hockey sticks” t↦max{t−c,0} pin down the partial sums.

Proposition 21.2.2 (Convex Functions Detect Majorization)

Let x,y∈Rn. If ∑i=1nϕ(xi) ≤ ∑i=1nϕ(yi) for every convex ϕ:R→R, then x≺y.

Idea. Use ϕ(t)=t and ϕ(t)=−t to get the equal totals. For the k-th partial sum, use the function that is 0 below a level c and rises with slope 1 above it. Summed over a vector it always exceeds “top k sum minus kc”, and it equals that when c is the k-th largest entry. Setting c=yk↓ makes the right-hand side exact and the left-hand side an over-estimate, which is the direction we need.

Proof. Both ϕ(t)=t and ϕ(t)=−t are affine, hence convex by Example (First convex functions) (a). They give ∑ixi≤∑iyi and −∑ixi≤−∑iyi, so the totals are equal. This is (M2).

Fix k with 1≤k≤n−1, put c=yk↓, and let ϕc(t)=max{t−c, 0}, the pointwise maximum of the two affine functions t↦t−c and t↦0, hence convex by Proposition 18.10.3. For any z∈Rn the sum ∑iϕc(zi) does not depend on the order of the entries, so it equals ∑iϕc(zi↓). Every term is ≥0, and each of the first k terms satisfies ϕc(zi↓)≥zi↓−c. Hence (∗)∑i=1nϕc(zi) ≥ ∑i=1kzi↓−kc. For z=y and this particular c, (∗) is an equality: for i≤k we have yi↓≥yk↓=c, so ϕc(yi↓)=yi↓−c; and for i>k we have yi↓≤c, so ϕc(yi↓)=0. Therefore ∑i=1kxi↓−kc ≤(∗) ∑i=1nϕc(xi) ≤ ∑i=1nϕc(yi)=∑i=1kyi↓−kc, the middle step being the hypothesis applied to ϕc. Canceling kc gives (M1) at k. As k was arbitrary, x≺y.

Together with Theorem 21.2.1 this is an equivalence: x≺y if and only if ∑iϕ(xi)≤∑iϕ(yi) for every convex ϕ:R→R. So majorization is exactly the order that all convex functions agree on. The proof also shows how few test functions are needed: the two affine ones and the n−1 hockey sticks at the levels y1↓,…,yn−1↓.

Weak majorization has a matching statement, with the class of test functions cut down.

Proposition 21.2.3 (Karamata for Weak Majorization)

Let x,y∈Rn with x≺wy, and let ϕ:R→R be convex and increasing. Then ∑i=1nϕ(xi) ≤ ∑i=1nϕ(yi).

Idea. The only thing wrong with y is that its total may be too big. Take the excess away from the smallest entry of y. That keeps the vector sorted, leaves all the partial sums before the last one untouched, and produces a genuine majorization, to which Theorem 21.2.1 applies. Lowering an entry can only lower ϕ, which is where “increasing” is used.

Proof. Let s=∑iyi−∑ixi. Taking k=n in (M1) for x≺wy gives s≥0. Define w∈Rn by wi=yi↓  (i≤n−1),wn=yn↓−s. Only the last entry has changed, and it has been lowered, so w is still decreasing and w↓=w. Its total is ∑iyi−s=∑ixi, which is (M2) for x≺w. For k≤n−1, ∑i≤kwi=∑i≤kyi↓≥∑i≤kxi↓ by (M1) for x≺wy. So x≺w.

By Theorem 21.2.1 with I=R, ∑iϕ(xi)≤∑iϕ(wi). Since wi≤yi↓ for every i and ϕ is increasing, ϕ(wi)≤ϕ(yi↓), so ∑iϕ(wi) ≤ ∑iϕ(yi↓)=∑iϕ(yi). Chaining the two inequalities proves the proposition.

Warning. “Increasing” cannot be dropped from Proposition 21.2.3. Take x=(0,0) and y=(1,0). The running totals 0,0 of x are at most 1,1, so x≺wy. The function ϕ(t)=−t is convex but decreasing, and ∑iϕ(xi)=0 while ∑iϕ(yi)=−1. A weak majorization allows y to carry more total than x, and a decreasing ϕ punishes exactly that.

21.2.3. Schur-convex functions

Theorem 21.2.1 is a statement about one particular shape of function, x↦∑iϕ(xi). The property it establishes deserves a name of its own, because many functions have it that are not of that shape.

A function is Schur-convex when spreading a vector out can only increase it.

The domain needs a moment’s thought. If G is to be compared at x and y whenever x≺y, then the domain must contain x as soon as it contains y. Call a set A⊆Rn symmetric if Pσx∈A for every x∈A and every σ∈Sn. A symmetric convex set has the property we need: if y∈A and x≺y, then x∈conv{Pσy} by Corollary 21.1.4, every Pσy lies in A by symmetry, and a convex set contains the convex combinations of its points (Lemma 18.1.2), so x∈A.

Definition 21.2.1 (Schur-Convex Function)

Let A⊆Rn be symmetric and convex, and let G:A→R. Then G is Schur-convex if x≺y⟹G(x)≤G(y) for all x,y∈A. It is strictly Schur-convex if the inequality is strict whenever x↓≠y↓.

Definition 21.2.2 (Schur-Concave Function)

With A as in Definition 21.2.1, a function G:A→R is Schur-concave if −G is Schur-convex, that is, if x≺y⟹G(x)≥G(y) for all x,y∈A; and strictly Schur-concave if −G is strictly Schur-convex.

In words: a Schur-convex function is one that respects the majorization order, and a Schur-concave one reverses it. The two standing examples of A below are Rn itself and the open orthant {x:xi>0 for every i}, both symmetric and convex.

A Schur-convex function is automatically symmetric, meaning G(Pσx)=G(x) for every permutation. Indeed x and Pσx are rearrangements of each other, so each majorizes the other (Proposition 17.8.1 (a) and the remark after Definition 17.8.1), and two applications of the definition give G(x)≤G(Pσx)≤G(x).

Example (Four functions on Rn)

Decide which of the following are Schur-convex, and which are Schur-concave, on A=Rn.

  • (a) G(x)=∑ixi2.
  • (b) G(x)=maxixi.
  • (c) G(x)=∑ixi.
  • (d) G(x)=∑ixi3, for n=3.

Solution.

  • (a) Schur-convex. ϕ(t)=t2 is convex on R, as in Example (One majorization, four inequalities), so Theorem 21.2.1 gives ∑ixi2≤∑iyi2 whenever x≺y.

  • (b) Schur-convex. maxixi=x1↓, and (M1) at k=1 says exactly x1↓≤y1↓. No convexity is needed here; the definition of majorization already contains the statement.

  • (c) Both Schur-convex and Schur-concave, since (M2) makes it constant on each majorization class. This is the degenerate case, and it shows that the two notions are not exclusive.

  • (d) Neither, and one pair of comparisons settles both. The flat vector is majorized by every vector with the same total (Example (A chain, the two extremes, and size one) (b)), so (0,0,0)≺(2,−1,−1),(0,0,0)≺(1,1,−2), both totals being 0. Now G(0,0,0)=0, while G(2,−1,−1)=8−1−1=6 and G(1,1,−2)=1+1−8=−6. Schur-convexity applied to the second majorization would demand 0≤−6, and Schur-concavity applied to the first would demand 0≥6. Both fail.

Part (d) is the nonexample by minimal change. The function ∑iϕ(xi) with ϕ(t)=t3 differs from the Schur-convex ∑iϕ(xi) with ϕ(t)=t2 only in the exponent, and t3 fails the one clause that matters, convexity on the interval where the entries live.

The next theorem explains why so many natural functions turn out Schur-convex: convexity plus symmetry is enough, and there is no need to be of the form ∑iϕ(xi).

Theorem 21.2.4 (A Symmetric Convex Function Is Schur-Convex)

Let A⊆Rn be symmetric and convex, and let G:A→R be convex and symmetric, that is, G(Pσx)=G(x) for every x∈A and every σ∈Sn. Then G is Schur-convex.

Idea. x≺y puts x inside the permutation polytope of y, whose corners all carry the same value of G, by symmetry. A convex function on a convex combination is at most the corresponding combination of its values, and an average of copies of one number is that number.

Proof. Let x,y∈A with x≺y. By Corollary 21.1.4 and Theorem 18.1.3, x=∑i=1mtiPσiy for some σ1,…,σm∈Sn and weights ti≥0 with ∑iti=1. Each Pσiy lies in A, since A is symmetric. By Theorem 18.10.2 applied to the convex G on A, G(x) ≤ ∑i=1mtiG(Pσiy)=∑i=1mtiG(y)=G(y), the middle equality being the symmetry of G and the last one ∑iti=1. This shows that G is Schur-convex.

Now the promise from Chapter 17 §08. There the second proof of Hadamard’s inequality used only the trace, and a remark recorded that a proof “by majorization in earnest” would need the fact that “for vectors with nonnegative entries, x≺y implies x1⋯xn≥y1⋯yn (the product is Schur-concave)”. Here it is.

Corollary 21.2.5 (The Product Is Schur-Concave)

Let y∈Rn have nonnegative entries and let x≺y. Then every entry of x is nonnegative, and x1x2⋯xn ≥ y1y2⋯yn.

Idea. Take −log, which is convex, and apply Theorem 21.2.1: the product becomes a sum and the inequality flips. The logarithm needs positive entries, so the case of a zero entry is dealt with first — and it is the easy case, because a zero entry makes the right-hand side 0 while the left-hand side is a product of nonnegative numbers.

Proof. First, every entry of x is nonnegative. For n=1, (M2) gives x=y. For n≥2, subtracting (M1) at k=n−1 from (M2) gives xn↓ ≥ yn↓ ≥ 0, and xn↓ is the smallest entry of x.

Case 1: some yi=0. Then y1⋯yn=0, while x1⋯xn≥0 as a product of nonnegative numbers. The inequality holds.

Case 2: every yi>0. Then yn↓>0, so xn↓>0 by the display above, and every entry of both vectors lies in I=(0,∞). By Lemma 18.11.1 (c) the function log is concave on I, so ϕ=−log is convex there. By Theorem 21.2.1, ∑i=1n(−log⁡xi) ≤ ∑i=1n(−log⁡yi),that is,∑i=1nlog⁡yi ≤ ∑i=1nlog⁡xi. By log⁡(ab)=log⁡a+log⁡b (Lemma 18.11.1 (c)) and induction on n, ∑ilog⁡xi=log⁡(x1⋯xn) and likewise for y. So log⁡(y1⋯yn)≤log⁡(x1⋯xn). Since log is strictly increasing on I (Lemma 18.11.1 (c)), two positive numbers with log⁡u≤log⁡v satisfy u≤v; otherwise u>v would give log⁡u>log⁡v. Hence y1⋯yn≤x1⋯xn. This proves the corollary.

So the flat vector, the bottom of the majorization order, has the largest product among nonnegative vectors with a fixed total — which is Lemma 17.8.6 again, now as a special case rather than a separate induction.

Warning. Schur-convex does not mean convex. By Corollary 21.2.5 the function G(x)=−x1x2 is Schur-convex on the symmetric convex set A={x∈R2:x1,x2≥0}. It is not convex: at x=(1,1), y=(2,2), and t=12, G(32,32)=−94=−2.25,12G(1,1)+12G(2,2)=12(−1)+12(−4)=−2.5, and −2.25>−2.5. Theorem 21.2.4 is a one-way street.

Quick check. Is the function G(x)=x1↓+x2↓ on R3 Schur-convex? Is it convex?

Solution. Schur-convex: the inequality x1↓+x2↓≤y1↓+y2↓ is (M1) at k=2. Convex: yes, and this gives a second proof. The two largest entries have the largest sum among all pairs, so G(x)=maxi<j(xi+xj), a maximum of three linear functions, convex by Proposition 18.10.3. It is also symmetric, so Theorem 21.2.4 applies. (This G is the Ky Fan sum s2 of Chapter 17 §06, read on vectors instead of eigenvalues.)

21.2.4. Convex functions of a Hermitian matrix

The last payoff is about matrices. For a Hermitian A and a convex f on an interval containing spec(A), the number trf(A) turns out to be ∑if(λi(A)). Chapter 18 §10 stated, and deferred, that this is a convex function of A. The proof needs two lemmas: a bookkeeping one about traces, and one saying that the quadratic form of f(A) dominates f of the quadratic form of A.

Lemma 21.2.6 (Reading a Trace on an Orthonormal Basis)

Let F=R or C.

  • (a) Let (u1,…,un) be an orthonormal basis of Fn and M∈Mn(F). Then trM=∑i=1n⟨Mui,ui⟩.
  • (b) Let I⊆R, let A∈Mn(F) be Hermitian with spec(A)⊆I, and let g:I→R. Then g(A)u=g(λ)u for every eigenvector u of A with eigenvalue λ, and trg(A)=∑i=1ng(λi(A)).

Proof.

  • (a) Let U∈Mn(F) be the matrix whose i-th column is ui. Then (U∗U)ij=ui∗uj=⟨uj,ui⟩=δij, so U is unitary, and (U∗MU)ii=ui∗Mui=⟨Mui,ui⟩. By Theorem 0.9.6 (3), used with the factors U∗ and MU, tr(U∗MU)=tr(MUU∗)=trM. Summing the diagonal of U∗MU gives (a).

  • (b) Let A=∑jμjPj be the spectral resolution of Definition 12.8.1, with μ1,…,μk the distinct eigenvalues; it exists by Theorem 12.8.1, since a Hermitian matrix has an orthonormal basis of eigenvectors (Corollary 12.4.3 over C, Corollary 12.5.2 over R). Let Au=λu with u≠0. Then λ=μj0 for exactly one j0, and u lies in the eigenspace of μj0, which is the image of Pj0; so Pj0u=u, and Pju=PjPj0u=0 for j≠j0 by Theorem 12.8.1 (b). Hence g(A)u=∑jg(μj)Pju=g(λ)u, by Definition 12.8.2.

Now take an orthonormal basis (u1,…,un) of Fn with Aui=λi(A)ui, which exists by the spectral theorem just cited. By the previous paragraph, ⟨g(A)ui,ui⟩=g(λi(A))⟨ui,ui⟩=g(λi(A)), and (a) adds these up to trg(A).

Lemma 21.2.7 (Peierls’ Inequality)

Let I⊆R be an interval, let A∈Mn(F) be Hermitian with spec(A)⊆I, let f:I→R be convex, and let u∈Fn with ‖u‖=1. Then ⟨Au,u⟩∈I and f(⟨Au,u⟩) ≤ ⟨f(A)u,u⟩, with equality when u is an eigenvector of A.

Idea. Write A=∑jμjPj, its spectral resolution. Then ⟨Au,u⟩=∑jμj‖Pju‖2 is a weighted average of the eigenvalues, with weights ‖Pju‖2 that add up to ‖u‖2=1. The same weights appear in ⟨f(A)u,u⟩=∑jf(μj)‖Pju‖2. So the two sides of the inequality are f of an average against the average of f, which is Jensen.

Proof. Let A=μ1P1+⋯+μkPk be the spectral resolution of Definition 12.8.1, with μ1,…,μk the distinct eigenvalues of A, which exists by Theorem 12.8.1 since a Hermitian matrix has an orthonormal basis of eigenvectors (Corollary 12.4.3 over C, Corollary 12.5.2 over R). Put wj=⟨Pju,u⟩. By Theorem 12.8.1 (b), Pj=Pj2 and Pj∗=Pj, so wj=⟨PjPju,u⟩=⟨Pju,Pj∗u⟩=⟨Pju,Pju⟩=‖Pju‖2 ≥ 0, and by Theorem 12.8.1 (c), ∑jwj=⟨u,u⟩=1. Hence, by Theorem 12.8.1 (d) and Definition 12.8.2, ⟨Au,u⟩=∑j=1kμjwj,⟨f(A)u,u⟩=∑j=1kf(μj)wj. The first is a convex combination of μ1,…,μk, all of which lie in the convex set I, so ⟨Au,u⟩∈I by Lemma 18.1.2. Applying Theorem 18.10.2 to the convex f on I, with points μj and weights wj, f(⟨Au,u⟩)=f(∑jwjμj) ≤ ∑jwjf(μj)=⟨f(A)u,u⟩.

For the equality clause, let Au=λu with ‖u‖=1. By Lemma 21.2.6 (b), f(A)u=f(λ)u, so the right side is f(λ)⟨u,u⟩=f(λ); and the left side is f(⟨λu,u⟩)=f(λ) as well. This proves the lemma.

Quick check. Take A=(2112), f(t)=t2, and u=e1. Compute both sides of Lemma 21.2.7, and then repeat with u=12(1,1).

Solution. ⟨Ae1,e1⟩=a11=2, so the left side is f(2)=4. Here f(A)=A2=(5445), so the right side is (A2)11=5. Indeed 4≤5, and the inequality is strict because e1 is not an eigenvector of A.

For u=12(1,1), an eigenvector with Au=3u, the left side is f(3)=9 and the right side is ⟨A2u,u⟩=9⟨u,u⟩=9. Equality, as the lemma says.

Theorem 21.2.8 (The Trace of a Convex Function Is Convex)

Let I⊆R be an interval and let f:I→R be convex. Let A,B∈Mn(F) be Hermitian with spec(A)⊆I and spec(B)⊆I, and let 0≤θ≤1. Then C=θA+(1−θ)B is Hermitian with spec(C)⊆I, and trf(C) ≤ θtrf(A)+(1−θ)trf(B). In words: the Hermitian matrices with spectrum in I form a convex set, and A↦trf(A) is a convex function on it.

Idea. Fix an orthonormal eigenbasis of the midpoint matrix C=θA+(1−θ)B. Along each of those n directions the value of f at the eigenvalue of C is exactly the left side of Peierls’ inequality, because an eigenvector gives equality; and the eigenvalue itself splits as θ times a quadratic form of A plus (1−θ) times one of B. Ordinary convexity of f handles that split, and Peierls’ inequality, now as an inequality, converts the two quadratic forms into quadratic forms of f(A) and f(B). Summing over the basis turns quadratic forms into traces.

Proof. Step 1: the spectrum of C. The matrix C=θA+(1−θ)B is Hermitian, being a real combination of Hermitian matrices. Let μ be an eigenvalue of C, with unit eigenvector x. Then μ=⟨Cx,x⟩=θ⟨Ax,x⟩+(1−θ)⟨Bx,x⟩. By Lemma 13.5.2, ⟨Ax,x⟩ lies in [λn(A),λ1(A)], an interval whose endpoints are eigenvalues of A, hence contained in I; likewise ⟨Bx,x⟩∈I. So μ is a convex combination of two points of I and lies in I. Hence spec(C)⊆I, and f(C) is defined.

Step 2: the inequality. Choose an orthonormal basis (u1,…,un) of Fn consisting of eigenvectors of C, which exists by Corollary 12.4.3 over C and Corollary 12.5.2 over R. Fix i. By Lemma 21.2.7 applied to C and the eigenvector ui, with equality, ⟨f(C)ui,ui⟩=f(⟨Cui,ui⟩). By Step 1, ⟨Cui,ui⟩=θαi+(1−θ)βi with αi=⟨Aui,ui⟩∈I and βi=⟨Bui,ui⟩∈I. So Definition 18.10.1 and then Lemma 21.2.7, now applied to A and to B with the unit vector ui, give ⟨f(C)ui,ui⟩=f(θαi+(1−θ)βi)≤θf(αi)+(1−θ)f(βi)≤θ⟨f(A)ui,ui⟩+(1−θ)⟨f(B)ui,ui⟩, where the last step uses θ≥0 and 1−θ≥0. Summing over i=1,…,n and using Lemma 21.2.6 (a) three times, once for each of f(C), f(A), and f(B), trf(C) ≤ θtrf(A)+(1−θ)trf(B). Since A, B, and θ were arbitrary, this is Definition 18.10.1 for A↦trf(A) on the set of Hermitian matrices with spectrum in I, which Step 1 shows is convex. This proves the theorem.

This pays Chapter 18 §10’s “it is proved in Chapter 21, and nothing in this chapter uses it”, and Chapter 18 §11’s closing note that “the convexity of A↦trf(A) for convex f, and the consequences of a majorization for every convex function, belong to Chapter 21”.

Example (The logarithm of the determinant)

Let f(t)=−log⁡t on I=(0,∞), so that the matrices in play are the positive definite Hermitian ones. Identify trf(A), and check the resulting inequality for A=(1004),B=(2112),θ=12.

Solution. A Hermitian A has spec(A)⊆(0,∞) exactly when A≻O, by Theorem 13.1.5 (b). For such an A, Lemma 21.2.6 (b), then log⁡(ab)=log⁡a+log⁡b (Lemma 18.11.1 (c)) with induction on n, and then Theorem 9.3.3 give trf(A)=−∑ilog⁡λi(A)=−log⁡∏iλi(A)=−log⁡detA. By Theorem 21.2.8, A↦−log⁡detA is convex on the positive definite matrices; equivalently, det is log-concave there.

For the data given, detA=4 and detB=4−1=3, and C=12(A+B)=(3212123),detC=92−14=174. Both A and B are positive definite (A is diagonal with positive entries; B has eigenvalues 3 and 1), so the theorem says −log⁡174≤−12log⁡4−12log⁡3, that is, after multiplying by −1 and exponentiating, detC=174=4.25 ≥ 4⋅3=12≈3.464, which holds, since 4.252=18.0625≥12.

Remark. Theorem 21.2.8 holds for every convex f, which is a very large class. It is nonetheless much weaker than the operator convexity of Definition 17.11.1 (b), which asks for f(θA+(1−θ)B)⪯θf(A)+(1−θ)f(B) in the Loewner order, a statement about matrices rather than about one number extracted from them. Taking traces turns the Loewner inequality into the one proved here, since the trace of a positive semidefinite matrix is the sum of its nonnegative eigenvalues (Theorem 13.1.4 (b), Theorem 9.3.3); so operator convexity implies trace convexity, but the converse fails badly, and a later section of this chapter exhibits convex functions that are not operator convex.

21.2.5. Exercises

A. Check your understanding

Exercise (A1)
  • (a) State Karamata’s inequality, with all its hypotheses.
  • (b) Which test functions prove Proposition 21.2.2, and which clause of Definition 17.8.1 does each one give?
  • (c) Define Schur-convex and Schur-concave.
  • (d) True or false: every Schur-convex function on Rn is convex. Justify your answer.
  • (e) True or false: every convex function on Rn is Schur-convex. Justify your answer.
  • (f) State what Theorem 21.2.8 says for f(t)=t2, and identify the resulting function of A.

Solution.

  • (a) If I is an interval, ϕ:I→R is convex, x≺y in Rn, and every entry of y lies in I, then every entry of x lies in I and ∑iϕ(xi)≤∑iϕ(yi) (Theorem 21.2.1).

  • (b) The affine functions ϕ(t)=t and ϕ(t)=−t give the equal totals, (M2). The hockey stick ϕc(t)=max{t−c,0} with c=yk↓ gives (M1) at k.

  • (c) For a symmetric convex A⊆Rn, G:A→R is Schur-convex if x≺y implies G(x)≤G(y) (Definition 21.2.1), and Schur-concave if x≺y implies G(x)≥G(y) (Definition 21.2.2).

  • (d) False. G(x)=−x1x2 on the nonnegative quadrant is Schur-convex by Corollary 21.2.5, but the warning after that corollary shows that G(32,32)=−2.25 exceeds the average −2.5 of G(1,1) and G(2,2).

  • (e) False. A Schur-convex function must be symmetric, as shown after Definition 21.2.2, and G(x)=x1 is linear, hence convex, but not symmetric: G(1,0)=1≠0=G(0,1), while (1,0) and (0,1) majorize each other. Symmetry is exactly the missing hypothesis of Theorem 21.2.4.

  • (f) For I=R and f(t)=t2, it says that A↦tr(A2) is convex on the Hermitian matrices. Since A is Hermitian, tr(A2)=tr(A∗A)=‖A‖F2, the squared Frobenius norm.

B. Practice

Exercise (B1: Four consequences of one majorization)

Let y=(6,2,1,1) and x=(4,3,2,1).

  • (a) Verify that x≺y.
  • (b) Compute ∑ixi2 and ∑iyi2, and ∑imax{xi−3,0} and ∑imax{yi−3,0}. Name the theorem each comparison illustrates.
  • (c) Compute ∏ixi and ∏iyi. Hence check Corollary 21.2.5.
  • (d) Give a convex ϕ:R→R for which the two sides are equal.

Solution.

  • (a) Both vectors are decreasing with total 10; the running totals are 4,7,9 for x against 6,8,9 for y. So (M1) and (M2) hold.

  • (b) ∑ixi2=16+9+4+1=30 and ∑iyi2=36+4+1+1=42, so 30≤42. Also ∑imax{xi−3,0}=1+0+0+0=1 and ∑imax{yi−3,0}=3, so 1≤3. Both illustrate Theorem 21.2.1, the first with ϕ(t)=t2 and the second with the hockey stick at level 3.

  • (c) ∏ixi=4⋅3⋅2⋅1=24 and ∏iyi=6⋅2⋅1⋅1=12. Indeed 24≥12, as Corollary 21.2.5 requires for nonnegative vectors.

  • (d) Any affine ϕ, for instance ϕ(t)=t: both sides are the common total 10, by (M2). Affine functions are the ones on which majorization gives no information.

Exercise (B2: Which are Schur-convex?)

Determine which of the following functions are Schur-convex on the stated symmetric convex set. Justify your answer, naming the theorem used or giving a counterexample.

  • (a) G(x)=maxixi−minixi on Rn.
  • (b) G(x)=∑i|xi| on Rn.
  • (c) G(x)=x1+2x2 on R2.
  • (d) G(x)=∑ixi on A={x∈Rn:xi>0 for every i}.

Solution.

  • (a) Schur-convex. G(x)=x1↓−xn↓, and x≺y gives x1↓≤y1↓ by (M1) at k=1 and xn↓≥yn↓ by subtracting (M1) at k=n−1 from (M2). Subtracting the second from the first gives G(x)≤G(y).

  • (b) Schur-convex, by Theorem 21.2.1 with ϕ(t)=|t|, convex as a maximum of two affine functions (Proposition 18.10.3).

  • (c) Not Schur-convex. It is not symmetric, and a Schur-convex function must be. Concretely (1,0) and (0,1) are rearrangements of each other, so (1,0)≺(0,1) and (0,1)≺(1,0). Schur-convexity would give G(1,0)≤G(0,1) from the first and G(0,1)≤G(1,0) from the second, hence G(1,0)=G(0,1). But G(1,0)=1 and G(0,1)=2.

  • (d) Not Schur-convex; it is Schur-concave. First, t↦t is concave on (0,∞): for a,b>0 and 0≤θ≤1, Lemma 17.8.6 gives 2ab≤a+b, so (θa+(1−θ)b)2=θ2a+(1−θ)2b+2θ(1−θ)ab≤θ2a+(1−θ)2b+θ(1−θ)(a+b)=θa+(1−θ)b, using θ2+θ(1−θ)=θ and (1−θ)2+θ(1−θ)=1−θ; taking square roots of two nonnegative numbers preserves the inequality. So ϕ(t)=−t is convex on (0,∞), and Theorem 21.2.1 turns x≺y into ∑i(−xi)≤∑i(−yi), that is, G(x)≥G(y). For a concrete instance, (2,2)≺(3,1) and G(2,2)=22≈2.83, while G(3,1)=3+1≈2.73, so G went down as the vector spread out.

Exercise (B3: Trace convexity on two-by-two matrices)
  • (a) Let f(t)=t2, A=diag(2,0), B=(1111), and θ=12. Compute trf(A), trf(B), and trf(C) for C=12(A+B), and verify Theorem 21.2.8.
  • (b) Let f(t)=−log⁡t, A=I2, B=diag(4,14), and θ=12. Do the same, and state the inequality in terms of determinants.

Solution.

  • (a) Both matrices are real symmetric, hence Hermitian. A2=diag(4,0), so trf(A)=4; and B2=(2222), so trf(B)=4. Next C=12(3111)=(32121212),C2=(521112), so trf(C)=52+12=3. The theorem asserts 3≤12⋅4+12⋅4=4, which holds.

  • (b) Both are positive definite, so their spectra lie in (0,∞). By Example (The logarithm of the determinant), trf(M)=−log⁡detM. Here detA=1 and detB=4⋅14=1, so both traces are −log⁡1=0. The midpoint is C=diag(52,58), with detC=2516, so trf(C)=−log⁡2516. The theorem asserts −log⁡2516 ≤ 12⋅0+12⋅0=0, which holds because 2516>1 and log is increasing with log⁡1=0 (Lemma 18.11.1 (c)). In determinant form: detC=2516≥1=detAdetB.

C. Going deeper

Exercise (C1: Entropy)

For x in A={x∈Rn:xi>0 for every i}, let H(x)=−∑ixilog⁡xi.

  • (a) Prove that g(t)=tlog⁡t is convex on (0,∞). Hint: apply Theorem 18.10.2 to the concave log, with the points 1/a, 1/b and cleverly chosen weights.
  • (b) Deduce that H is Schur-concave on A.
  • (c) Deduce that H(x)≤log⁡n for every x∈A with ∑ixi=1, with equality when x=(1n,…,1n).

Solution.

  • (a) Let a,b>0 and 0≤θ≤1, and put c=θa+(1−θ)b>0. Set w1=θac,w2=(1−θ)bc, which are nonnegative and satisfy w1+w2=(θa+(1−θ)b)/c=1. The function −log is convex on (0,∞) by Lemma 18.11.1 (c), so Theorem 18.10.2 with the points 1/a and 1/b gives −log⁡(w11a+w21b) ≤ −w1log⁡1a−w2log⁡1b. Now w1/a+w2/b=θ/c+(1−θ)/c=1/c, and log⁡(1/a)=−log⁡a since log⁡a+log⁡(1/a)=log⁡1=0 (Lemma 18.11.1 (c)). So the display reads log⁡c ≤ w1log⁡a+w2log⁡b=θalog⁡a+(1−θ)blog⁡bc. Multiplying by c>0 gives clog⁡c≤θalog⁡a+(1−θ)blog⁡b, which is Definition 18.10.1 for g.

  • (b) By (a) and Theorem 21.2.1 with I=(0,∞), x≺y in A implies ∑ixilog⁡xi≤∑iyilog⁡yi, that is, H(x)≥H(y). By Definition 21.2.2, H is Schur-concave. (The set A is symmetric and convex, so Definition 21.2.2 applies.)

  • (c) Let x∈A with ∑ixi=1. Its mean is 1n, so the flat vector f=(1n,…,1n) satisfies f≺x by Example (A chain, the two extremes, and size one) (b), and f∈A. By (b), H(f)≥H(x). Finally H(f)=−∑i=1n1nlog⁡1n=−log⁡1n=log⁡n, using log⁡1n=−log⁡n. So H(x)≤log⁡n, with equality at x=f.

Exercise (C2: Three norms under majorization)

Let x≺y in Rn. Prove that ‖x‖1≤‖y‖1,‖x‖2≤‖y‖2,‖x‖∞≤‖y‖∞. Which of the three needs Theorem 21.2.1, and which follows from Definition 17.8.1 alone?

Solution. For ‖⋅‖1, apply Theorem 21.2.1 with ϕ(t)=|t|, convex on R by Proposition 18.10.3: ∑i|xi|≤∑i|yi|.

For ‖⋅‖2, apply Theorem 21.2.1 with ϕ(t)=t2, convex on R by Lemma 18.10.5, since its derivative 2t is increasing: ∑ixi2≤∑iyi2. Both sides are nonnegative and t↦t is increasing on [0,∞), so ‖x‖2≤‖y‖2.

For ‖⋅‖∞, no convexity is needed. Since ‖z‖∞=maxi|zi|=max{z1↓, −zn↓}, and (M1) at k=1 gives x1↓≤y1↓ while subtracting (M1) at k=n−1 from (M2) gives xn↓≥yn↓, hence −xn↓≤−yn↓. The maximum of the two smaller numbers is at most the maximum of the two larger ones.

So the ∞-norm follows from Definition 17.8.1 directly, and the other two from Theorem 21.2.1. A later section of this chapter proves the same for every norm that is invariant under permuting and changing the signs of the coordinates.

Exercise (C3: One convex function is not enough)

Let x=(2,−1,−1) and y=(1,1,−2).

  • (a) Show that ∑ixi=∑iyi and ∑ixi2=∑iyi2, yet neither x≺y nor y≺x.
  • (b) Exhibit a convex ϕ with ∑iϕ(xi)>∑iϕ(yi), and another with the inequality reversed. Explain how this is consistent with Proposition 21.2.2.

Solution.

  • (a) Both totals are 0, and ∑ixi2=4+1+1=6=1+1+4=∑iyi2. Both vectors are already decreasing, with running totals 2,1 for x and 1,2 for y. (M1) at k=1 fails for x≺y, since 2>1, and (M1) at k=2 fails for y≺x, since 2>1.

  • (b) Take ϕ(t)=max{t−1,0}, convex by Proposition 18.10.3. Then ∑iϕ(xi)=1+0+0=1 and ∑iϕ(yi)=0+0+0=0, so 1>0. This is the test function of Proposition 21.2.2 at c=y1↓=1, and it certifies x⊀y.

Take ψ(t)=max{t+1,0}. Then ∑iψ(xi)=3+0+0=3 and ∑iψ(yi)=2+2+0=4, so 3<4, certifying y⊀x by the same proposition applied with the roles exchanged, at c=x2↓=−1.

There is no contradiction: Proposition 21.2.2 asks the inequality to hold for every convex ϕ. Agreeing on the affine functions and on t2 is far from agreeing on all of them, and here two hockey sticks pull in opposite directions.

Found a typo or a mistake? Report it, or read this page’s source

The text of this book is licensed under CC BY-NC-SA 4.0.