17.8  The Diagonal and the Spectrum

A Hermitian matrix carries two lists of n real numbers: its diagonal entries, read off at a glance, and its eigenvalues, computed with effort. Since Chapter 9 we have known that the two lists have the same sum, the trace. This section asks what else the eigenvalues force on the diagonal, and answers the question completely. The answer needs a name, majorization. One half of it, Schur’s theorem, is §06’s bound on partial sums of the diagonal in the new language. The other half, Horn’s theorem, says that nothing more is forced; it is proved by an induction driven by a single 2×2 rotation. The section ends with a short second proof of Hadamard’s determinant inequality.

Throughout, the field is R or C, and every matrix whose eigenvalues are indexed is Hermitian, with its eigenvalues real and listed decreasingly with multiplicity, λ1(A)≥⋯≥λn(A). Two vectors in Rn are attached to such an A∈Mn(F): d(A)=(a11,a22,…,ann),λ(A)=(λ1(A),λ2(A),…,λn(A)). The first is the diagonal vector, in the order the entries sit on the diagonal; the second is the eigenvalue vector, always in decreasing order. The diagonal entries are real because A∗=A gives aii=aii―, and the eigenvalues are real by Theorem 12.2.1. The bold d(A) is a vector, not the diagonal matrix diag(a11,…,ann).

Two facts about a Hermitian A are used below. The spectral theorem (Corollary 12.4.3 over C, Corollary 12.5.2 over R) makes A similar to the real diagonal matrix of its eigenvalues, so its characteristic polynomial splits as ∏i(x−λi(A)). Then Theorem 9.3.3 gives trA=∑i=1nλi(A),detA=∏i=1nλi(A).

17.8.1. What the diagonal already knows

Start with the smallest interesting case. The matrix A=(2112) has d(A)=(2,2) and eigenvalues 3 and 1, with eigenvectors (1,1) and (1,−1), so λ(A)=(3,1). The totals agree, 2+2=4=3+1. But the two lists are not alike: the diagonal is perfectly even, while the spectrum is spread out.

A larger example shows the pattern better. The tridiagonal matrix A=(320242025) has eigenvectors (1,2,2), (−2,−1,2) and (2,−2,1), with eigenvalues 7, 4 and 1. To check the first, A(1,2,2)=(3+4, 2+8+4, 4+10)=(7,14,14); the other two are checked the same way. So λ(A)=(7,4,1), while d(A)=(3,4,5). Sort the diagonal decreasingly to (5,4,3) and compare the running totals of the two lists: k123largest k diagonal entries5912largest k eigenvalues71112 Each running total of the eigenvalues is at least the matching running total of the diagonal, and at k=n the two are equal. None of this is new. Corollary 17.6.3 says that any k diagonal entries add up to at most λ1(A)+⋯+λk(A), with equality at k=n, and the k largest diagonal entries are one such choice. What is new is the point of view: to ask which diagonals are possible, the corollary’s list of inequalities has to be read as a single relation between two vectors.

Two features of this comparison deserve notice. First, the diagonal had to be sorted, and this is forced. For a permutation matrix Pσ (Definition 3.4.2), Pσ⊤=Pσ−1 by Lemma 3.4.7 (b), so Pσ⊤APσ is Hermitian with the same eigenvalues as A, and its (j,j) entry is eσ(j)⊤Aeσ(j)=aσ(j)σ(j). Relabeling the basis shuffles the diagonal and fixes the spectrum, so any true relation between d(A) and λ(A) must ignore the order of the diagonal. Second, the comparison is between running totals from the top, and it measures how unevenly a fixed total is spread.

17.8.2. Majorization

One vector is majorized by another when the two have the same total but the second piles more of it onto its largest entries: every “top k” sum of the second is at least that of the first.

The sorting step gets a notation first. For x=(x1,…,xn)∈Rn, the decreasing rearrangement x↓=(x1↓,…,xn↓) is the vector with the same entries, repeated as often as in x, listed so that x1↓≥x2↓≥⋯≥xn↓. For example, (0,3,1,3)↓=(3,3,1,0). In this notation, λ(A)=λ(A)↓ by convention, while d(A) usually needs sorting.

Definition 17.8.1 (Majorization)

Let n≥1 and let x,y∈Rn be real vectors with the same number of entries. We say that x is majorized by y, and write x≺y, if

  • (M1) ∑i=1kxi↓≤∑i=1kyi↓ for every k=1,…,n−1; and
  • (M2) ∑i=1nxi=∑i=1nyi (equal totals).

We say that x is weakly majorized by y, and write x≺wy, if the inequality in (M1) holds for every k=1,…,n, including k=n, and (M2) is dropped.

In words, clause by clause. (M1) sorts both vectors and compares their running totals from the top: the largest entry of x is at most the largest entry of y, the two largest of x add up to at most the two largest of y, and so on up to k=n−1. The word every matters: one failing k refutes x≺y. (M2) says the totals are equal; a sum does not care about order, so no sorting is needed there. Weak majorization keeps (M1), adds its case k=n, which reads ∑ixi≤∑iyi, and asks nothing more. Comparing the definitions: x≺y holds exactly when x≺wy and the totals are equal.

There is an equivalent way to read (M1) from the bottom. Subtracting the k-th running total from the equal totals of (M2) gives ∑i=k+1nxi↓ ≥ ∑i=k+1nyi↓(k=1,…,n−1), so under (M2), the n−k smallest entries of x add up to at least the n−k smallest entries of y. The vector y is more extreme at both ends: its top is heavier and its bottom is lighter.

The definition makes sense because x↓ is determined by x: when two entries are equal, it does not matter which is written first. A consequence worth recording: the relation sees only which numbers occur in each vector and how often, not where they sit. If x′ is a rearrangement of x and y′ of y, then x≺y if and only if x′≺y′, because the rearrangements have the same decreasing rearrangements and the same totals.

Example (A chain, the two extremes, and size one)

Check each claim against Definition 17.8.1.

  • (a) (1,1,1)≺(2,1,0)≺(3,0,0).
  • (b) For every x∈Rn with mean m=1n∑ixi, the flat vector m1=(m,…,m) satisfies m1≺x.
  • (c) For every x∈Rn with non-negative entries and total s=∑ixi, x≺(s,0,…,0).
  • (d) For n=1, x≺y if and only if x=y.

Solution.

  • (a) All three vectors are already decreasing, and all have total 3, so (M2) holds for both pairs. Their running totals are (1,1,1): 1, 2, 3;(2,1,0): 2, 3, 3;(3,0,0): 3, 3, 3. At k=1 and k=2 the first row is at most the second (1≤2, 2≤3), and the second is at most the third (2≤3, 3≤3). So (M1) holds for both pairs.

  • (b) The totals are nm=∑ixi, so (M2) holds. For (M1), fix k and let a=1k∑i≤kxi↓ be the average of the k largest entries. Each of these is at least xk↓, so a≥xk↓. Every remaining entry xj↓ with j>k is at most xk↓, hence at most a. Therefore nm=∑i≤kxi↓+∑j>kxj↓≤ka+(n−k)a=na, so m≤a, and multiplying by k gives km≤∑i≤kxi↓. The left side is the k-th running total of m1. This is (M1).

  • (c) The totals are both s. Every running total of (s,0,…,0) equals s, and ∑i≤kxi↓≤s because the omitted entries xk+1↓,…,xn↓ are ≥0. This is (M1). Non-negativity is really needed: (1,−1) has total 0, but (1,−1)≺(0,0) fails at k=1, since 1>0.

  • (d) For n=1 there are no values of k in (M1), since k runs from 1 to 0. So only (M2) remains, and it says x1=y1.

Parts (b) and (c) are the two ends of the order: among vectors with a given total, the flat vector is the most evenly spread, and among non-negative ones, the vector that puts everything in one entry is the least. Part (d) is the degenerate case: with a single entry there is nothing to spread, and the relation collapses to equality.

The chain in (a) might suggest that any two vectors with the same total can be ranked this way. They cannot.

Example (Neither majorizes the other)

Show that x=(4,1,1) and y=(3,3,0) have the same total, that neither x≺y nor y≺x, and that both are majorized by (4,2,0).

Solution. Both totals are 6, and both vectors are already decreasing. Their running totals are x: 4, 5, 6;y: 3, 6, 6. At k=1, 4>3, so (M1) fails for x≺y. At k=2, 6>5, so (M1) fails for y≺x. The vector (4,2,0) has total 6 and running totals 4,6,6, which are at least both 4,5 and 3,6 at k=1,2. So x≺(4,2,0) and y≺(4,2,0).

This is the non-example by minimal change. Start from the comparable pair (4,1,1)≺(4,2,0), move one unit in the second vector from its top entry to its middle entry to get (3,3,0), and the pair becomes incomparable. (M2) still holds; what fails is (M1), at a different k in each direction. Spreading is not a single number.

The picture makes this visible. Plot the running totals of a sorted vector against k, joined by segments, starting from (0,0). Then x≺y says that the curve of y lies on or above the curve of x and ends at the same height. Two curves can also cross.

Diagram. Labels: k, tk, sv, solid: (3,0,0)\\ dashed: (2,1,0)\\ dotted: (1,1,1)\\ nested, same endpoint, solid: (4,1,1)\\ dashed: (3,3,0)\\ the curves cross:\\ incomparable

Why this definition. Each clause earns its place, and changing any one of them changes the subject.

  • Sorting. Relabeling a basis shuffles d(A) and fixes λ(A), so a relation comparing unsorted entries would answer a question the matrix does not pose.
  • Running totals from the top, with ≤. These are the quantities that Theorem 17.6.2 and Corollary 17.6.3 control: the top k eigenvalues sum to a maximum. Totals from the bottom add nothing new, as shown above, because under (M2) they are the same information.
  • Equal totals. This is the trace. Dropping it gives weak majorization, the right notion when there is no conservation law to enforce: (1,0)≺w(2,0), but (1,0)⊀(2,0), because the totals differ.

Warning. x≺y does not mean xi≤yi, not even after sorting both vectors. In (1,1,1)≺(2,1,0), the third entries are 1 and 0, and 1>0. This is no accident: subtracting (M1) at k=n−1 from (M2) gives xn↓≥yn↓ whenever x≺y, so at the bottom the inequality always runs the other way. Majorization compares sums of the largest entries, not entries one at a time.

Quick check. Arrange the five vectors (2,1,1,0), (4,0,0,0), (1,1,1,1), (3,1,0,0) and (2,2,0,0) in a chain x(1)≺x(2)≺⋯≺x(5).

Solution. All five are decreasing and all have total 4. Their running totals at k=1,2,3 are (1,1,1,1): 1, 2, 3;(2,1,1,0): 2, 3, 4;(2,2,0,0): 2, 4, 4;(3,1,0,0): 3, 4, 4;(4,0,0,0): 4, 4, 4. Read in this order, each triple is entrywise at most the next, so (1,1,1,1)≺(2,1,1,0)≺(2,2,0,0)≺(3,1,0,0)≺(4,0,0,0). The ends are the flat vector and the concentrated vector of Example (A chain, the two extremes, and size one) (b) and (c), as they must be. These five happen to be comparable in pairs, but Example (Neither majorizes the other) shows that this is not always so.

17.8.3. Majorization is an order

The first small result is that ≺ behaves like an order, and it says precisely how it falls short of being a partial order on Rn.

Proposition 17.8.1 (Basic Properties of Majorization)

Let x,y,z∈Rn.

  • (a) x≺x.
  • (b) If x≺y and y≺z, then x≺z.
  • (c) If x≺y and y≺x, then x↓=y↓; that is, y is a rearrangement of x.

Statements (a) and (b) hold with ≺w in place of ≺.

Proof. For k=1,…,n, write Xk=∑i≤kxi↓, and define Yk and Zk in the same way. Since Xn=∑ixi, Definition 17.8.1 says that x≺y means Xk≤Yk for k<n and Xn=Yn.

  • (a) Xk≤Xk for every k, and Xn=Xn.

  • (b) For k<n, Xk≤Yk≤Zk. At k=n, Xn=Yn=Zn. Hence x≺z.

  • (c) The hypotheses give Xk≤Yk and Yk≤Xk for k<n, and Xn=Yn, so Xk=Yk for every k=1,…,n. Put X0=Y0=0. Then for each k, xk↓=Xk−Xk−1=Yk−Yk−1=yk↓, so x↓=y↓. Two vectors with the same decreasing rearrangement contain the same numbers, each the same number of times, so each is a rearrangement of the other.

For ≺w, the arguments for (a) and (b) are the same with every k≤n treated like k<n. This proves the proposition.

So ≺ is reflexive and transitive. A relation with these two properties is called a preorder. Part (c) is antisymmetry up to rearrangement, and no more: (1,0)≺(0,1)≺(1,0), yet (1,0)≠(0,1). So ≺ is not a partial order on Rn. Restricted to decreasing vectors, where a rearrangement is the vector itself, (c) is exact antisymmetry, and there ≺ is a partial order. By Example (Neither majorizes the other) it is not total, not even among non-negative decreasing vectors with the same total. Like the Loewner order of Chapter 13 (Proposition 13.5.1), it leaves many pairs uncompared.

17.8.4. Schur’s theorem

The running totals in the 3×3 example were not a coincidence. Restated in the language of majorization, §06’s Corollary 17.6.3 says that for every Hermitian matrix the diagonal is more evenly spread than the spectrum:

Theorem 17.8.2 (Schur’s Majorization Theorem)

Let F=R or F=C, and let A∈Mn(F) be Hermitian. Then d(A)≺λ(A). Explicitly, if d(A)↓=(a1↓,…,an↓) lists the diagonal entries of A decreasingly, then a1↓+⋯+ak↓≤λ1(A)+⋯+λk(A)(1≤k≤n−1),a1↓+⋯+an↓=λ1(A)+⋯+λn(A).

Proof. Both vectors lie in Rn, as noted at the start of the section. Fix k with 1≤k≤n. List the positions 1,…,n in an order in which the diagonal entries decrease, breaking ties arbitrarily, and let i1<⋯<ik be the first k positions of that list, written in increasing order. Then ai1i1,…,aikik are a1↓,…,ak↓ in some order, so Corollary 17.6.3 gives a1↓+⋯+ak↓=∑j=1kaijij≤λ1(A)+⋯+λk(A), with equality when k=n. For k≤n−1 this is (M1). For k=n it is (M2), since a sum does not depend on the order of its terms. This proves the theorem.

So Schur’s theorem is Corollary 17.6.3 restated. What the restatement adds is a relation between two vectors, and a relation has a converse that can be asked about. Read from the bottom, (M1) says that the n−k smallest diagonal entries add up to at least the n−k smallest eigenvalues. At the two ends of Example (A chain, the two extremes, and size one), a diagonal matrix has d(A) a rearrangement of λ(A), the most spread a diagonal can be; a constant diagonal t is the flat vector t1, where Schur’s theorem says only that nt=λ1(A)+⋯+λn(A).

Example (Coordinate vectors against eigenvectors)

For the matrix A=(320242025) from the start of the section, whose running totals were tabulated there, exhibit for k=2 the coordinate pair that the proof of Corollary 17.6.3 feeds to Theorem 17.6.2 and an orthonormal pair that attains the maximum in Theorem 17.6.2. Account for the difference between the two sums.

Solution. We found λ(A)=(7,4,1), with unit eigenvectors u1=13(1,2,2), u2=13(−2,−1,2) and u3=13(2,−2,1). They are pairwise orthogonal (for instance ⟨(1,2,2),(−2,−1,2)⟩=−2−2+4=0), and each vector in parentheses has length 3.

For k=2 the proof uses the two largest diagonal entries, a33=5 and a22=4, so the coordinate pair is (e2,e3), and ⟨Ae2,e2⟩+⟨Ae3,e3⟩=9. The pair (u1,u2) gives ⟨Au1,u1⟩+⟨Au2,u2⟩=7+4=11, since ⟨Aui,ui⟩=λi‖ui‖2=λi. The coordinate pair falls short of the eigenvector pair by 2, and that shortfall is exactly the slack in Schur’s theorem at k=2. Ky Fan’s maximum is attained by eigenvectors, and the diagonal only ever sees coordinate vectors.

Quick check. A Hermitian 3×3 matrix A has diagonal entries 4,−1,2. Without knowing anything else about A, what does Theorem 17.8.2 say about λ1(A), about λ1(A)+λ2(A), and about λ3(A)?

Solution. Sorted, d(A)↓=(4,2,−1), with running totals 4,6,5. So λ1(A)≥4, λ1(A)+λ2(A)≥6, and λ1(A)+λ2(A)+λ3(A)=5. Subtracting the middle inequality from the last equality gives λ3(A)≤5−6=−1. In particular A is indefinite, since λ1(A)>0>λ3(A), and that was visible without computing any eigenvalue.

17.8.5. Horn’s converse

Schur’s theorem gives a necessary condition on a diagonal. Is it also sufficient? Given d≺λ, is there a Hermitian matrix with eigenvalues λ1,…,λn and diagonal entries d1,…,dn, in the given order? For n=2 a rotation answers the question. Turning the eigenbasis of diag(λ1,λ2) moves its first diagonal entry through every value between λ2 and λ1, and the trace fixes the second entry.

Lemma 17.8.3 (The Two-by-Two Rotation)

Let λ1≥λ2 and t be real numbers with λ2≤t≤λ1. Put b=(λ1−t)(t−λ2). Then there is an orthogonal R∈O(2) with Rdiag(λ1,λ2)R⊤=(tbbλ1+λ2−t).

Idea. If the first row of R is (c,−s) with c2+s2=1, the (1,1) entry of Rdiag(λ1,λ2)R⊤ is c2λ1+s2λ2, an average of λ1 and λ2 with weights c2 and s2. Every t between them is such an average, and solving for the weights gives c and s.

Proof. If λ1=λ2, then t=λ1=λ2 and b=0, and R=I works. Otherwise λ1−λ2>0, and c=t−λ2λ1−λ2,s=λ1−tλ1−λ2 are real numbers, because both radicands are ≥0 by the hypothesis on t. They satisfy c2+s2=1. Let R=(c−ssc). Then R⊤R=(c2+s2)I=I, so R∈O(2). Multiplying out, Rdiag(λ1,λ2)R⊤=(c2λ1+s2λ2cs(λ1−λ2)cs(λ1−λ2)s2λ1+c2λ2). Since s2=1−c2, the (1,1) entry is λ2+c2(λ1−λ2)=λ2+(t−λ2)=t. The two diagonal entries add up to (c2+s2)(λ1+λ2)=λ1+λ2, so the (2,2) entry is λ1+λ2−t. Finally cs(λ1−λ2)=(t−λ2)(λ1−t)=b. This proves the lemma.

This is the converse for n=2: if λ1≥λ2 and (d1,d2)≺(λ1,λ2), then d1,d2≤d1↓≤λ1 by (M1), so by (M2) d1=λ1+λ2−d2≥λ2, so t=d1 is allowed. For larger n, the lemma is used once, on a well-chosen pair of eigenvalues, to put d1 in the corner, and induction does the rest.

Theorem 17.8.4 (Horn’s Theorem)

Let n≥1, and let d=(d1,…,dn) and λ=(λ1,…,λn) be vectors in Rn with d≺λ. Then there is an orthogonal Q∈O(n) such that the matrix A=Qdiag(λ1,…,λn)Q⊤ has diagonal entries d1,d2,…,dn, in this order. So A is a real symmetric matrix with d(A)=d and λ(A)=λ↓.

Idea. Induction on n, one diagonal entry at a time.

① Permutation matrices are orthogonal and shuffle diagonal entries, so it is enough to treat d and λ both decreasing. ② The top and bottom clauses of majorization put d1 between λn and λ1, hence between two consecutive eigenvalues, λk≥d1≥λk+1. The lemma rotates that pair so that d1 sits in the corner, with η=λk+λk+1−d1 next to it. ③ The Claim: what is left over, the diagonal d2,…,dn and the spectrum with λk,λk+1 replaced by η, satisfies the hypothesis again, in size n−1. ④ Induction rotates the lower (n−1)×(n−1) block into place without touching the corner.

Step ③ is where the choice of k matters.

Proof. For every Q∈O(n), the matrix A=Qdiag(λ1,…,λn)Q⊤ is real and symmetric, and by Corollary 12.5.2 its eigenvalues, with multiplicity, are λ1,…,λn; listed decreasingly they form λ↓. So only the diagonal needs work. We use two facts repeatedly. A product of orthogonal matrices is orthogonal, since (XY)⊤XY=Y⊤X⊤XY=I. And for a permutation π∈Sn and any M∈Mn(R), Definition 3.4.2 gives (†)(Pπ⊤MPπ)ij=eπ(i)⊤Meπ(j)=mπ(i)π(j), where Pπ is orthogonal by Lemma 3.4.7 (b).

Step 1: reduction to decreasing vectors. Suppose the theorem holds whenever d and λ are both decreasing, and let d≺λ be arbitrary. Majorization sees only rearrangements (as noted after Definition 17.8.1), so d↓≺λ↓, and there is Q0∈O(n) such that B=Q0diag(λ↓)Q0⊤ has diagonal d↓. Here diag(λ↓) is the diagonal matrix with diagonal λ1↓,…,λn↓. Since λ↓ is a rearrangement of λ, and d↓ of d, there are permutations σ,τ∈Sn with λσ(i)=λi↓ and dτ(j)↓=dj for all i,j. By (†), Pσ⊤diag(λ1,…,λn)Pσ=diag(λ↓), and the (j,j) entry of Pτ⊤BPτ is bτ(j)τ(j)=dτ(j)↓=dj. So Q=Pτ⊤Q0Pσ⊤ is orthogonal and Qdiag(λ1,…,λn)Q⊤=Pτ⊤Q0diag(λ↓)Q0⊤Pτ=Pτ⊤BPτ has diagonal d. It therefore suffices to prove the theorem for decreasing d and λ, and we do so by induction on n.

Step 2: the corner. For n=1, d≺λ means d1=λ1 (Example (A chain, the two extremes, and size one) (d)), and Q=(1) works. Let n≥2, assume the theorem for n−1 (for vectors in any order, by Step 1), and let d and λ be decreasing with d≺λ. By (M1) at k=1, d1≤λ1. Subtracting (M1) at k=n−1 from (M2) gives dn≥λn, so d1≥dn≥λn. Let k be the smallest index in {1,…,n−1} with λk+1≤d1; one exists, since λn≤d1. Then λk≥d1: for k=1 this is d1≤λ1, and for k≥2 the minimality of k gives λk>d1. So λk≥d1≥λk+1. By Lemma 17.8.3 with t=d1, there are R∈O(2) and a real b with Rdiag(λk,λk+1)R⊤=(d1bbη),η=λk+λk+1−d1. Moreover λk+1≤η≤λk, since η−λk+1=λk−d1≥0 and λk−η=d1−λk+1≥0.

Step 3: the leftover data. Put d′=(d2,…,dn),λ′=(η,λ1,…,λk−1,λk+2,…,λn), two vectors in Rn−1. So λ′ is λ with the pair λk,λk+1 replaced by the single number η, written first.

Claim. d′≺λ′.

Proof. By Step 2, λk−1≥λk≥η≥λk+1≥λk+2, for those of these indices that exist. So (μ1,…,μn−1)=(λ1,…,λk−1,η,λk+2,…,λn) is (λ′)↓, and d′ is already decreasing. For (M2), ∑idi′=∑idi−d1=∑iλi−d1=∑iλi′, by (M2) for d≺λ and the definition of η. For (M1), fix m with 1≤m≤n−2; when n=2 there is nothing to check.

Case m≤k−1. Then μi=λi for i≤m. For each such i, di+1≤d1≤λk≤λi, because d and λ are decreasing and i<k. Summing over i≤m gives d2+⋯+dm+1≤μ1+⋯+μm.

Case m≥k. The first m entries of (λ′)↓ are λ1,…,λk−1, then η, then λk+2,…,λm+1. So μ1+⋯+μm=(λ1+⋯+λm+1)−λk−λk+1+η=(λ1+⋯+λm+1)−d1. Since m+1≤n−1, (M1) for d≺λ applies at m+1 and gives d2+⋯+dm+1=(d1+⋯+dm+1)−d1≤(λ1+⋯+λm+1)−d1=μ1+⋯+μm. In both cases (M1) holds, which proves the claim.

Step 4: assembly. Let σ∈Sn list the indices k,k+1 first and the others after them in increasing order: σ(1)=k, σ(2)=k+1, and σ(3),…,σ(n) are 1,…,k−1,k+2,…,n. By (†), Pσ⊤diag(λ1,…,λn)Pσ=(diag(λk,λk+1)00D1), where D1=diag(λ1,…,λk−1,λk+2,…,λn)∈Mn−2(R); when n=2 the blocks involving D1 are absent. Let R1=(R00In−2). Multiplying blocks and using Step 2, M=R1Pσ⊤diag(λ1,…,λn)PσR1⊤=(d1b0bη000D1)=(d1be1⊤be1E), where e1 is the first standard basis vector of Rn−1 and E=diag(η,λ1,…,λk−1,λk+2,…,λn) is the diagonal matrix of λ′. By the Claim and the inductive hypothesis, there is W∈O(n−1) such that WEW⊤ has diagonal entries d2,…,dn. Let W1=(10⊤0W). Then W1MW1⊤=(d1b(We1)⊤bWe1WEW⊤), whose diagonal entries are d1,d2,…,dn. The block matrices R1 and W1 are orthogonal, since R1⊤R1 and W1⊤W1 are block diagonal with blocks R⊤R=I, I, and 1, W⊤W=I. So Q=W1R1Pσ⊤ is orthogonal, and Qdiag(λ1,…,λn)Q⊤=W1MW1⊤. This completes the induction and proves the theorem.

Over C nothing more is needed, because an orthogonal matrix is unitary and a real symmetric matrix is Hermitian. Together with Schur’s theorem, this settles the question the section began with.

Theorem 17.8.5 (Schur–Horn Theorem)

Let F=R or F=C, and let d,λ∈Rn. There is a Hermitian A∈Mn(F) with λ(A)=λ↓ and d(A)=d if and only if d≺λ. In words: the diagonals of the Hermitian matrices with spectrum λ are exactly the vectors majorized by λ, and real symmetric matrices already give all of them.

Proof. (⇒) By Theorem 17.8.2, d=d(A)≺λ(A)=λ↓. Since majorization sees only rearrangements, d≺λ.

(⇐) Theorem 17.8.4 gives a real symmetric A∈Mn(R)⊆Mn(F) with d(A)=d and λ(A)=λ↓, and a real symmetric matrix is Hermitian.

Example (Building a matrix from its diagonal and spectrum)

Find a real symmetric matrix with eigenvalues 6,1,0 and diagonal (4,2,1).

Solution. Both vectors are decreasing, and (4,2,1)≺(6,1,0): the running totals are 4,6 against 6,7, and both totals are 7. Follow the proof of Theorem 17.8.4.

The corner. d1=4 lies between λ1=6 and λ2=1, so k=1, η=6+1−4=3, and Lemma 17.8.3 gives b=(6−4)(4−1)=6. After the first rotation the matrix is M=(460630000).

The leftover. d′=(2,1) and λ′=(3,0), with 2≤3 and 2+1=3+0, so d′≺λ′, as the Claim promises. The lemma with t=2 gives c=2/3, s=1/3, and Wdiag(3,0)W⊤=(2221),We1=(c,s).

Assembly. The first row becomes (4, 6c, 6s)=(4,2,2), so A=(422222221). Check directly: A(3,2,2)=(18,12,62)=6(3,2,2), A(−2,2,1)=(−2,2,1), and A(0,−1,2)=(0,0,0). Three eigenvectors with distinct eigenvalues 6,1,0 give the whole spectrum, and the diagonal is (4,2,1).

17.8.6. Hadamard’s inequality from the diagonal

Chapter 13 §06 proved Theorem 13.6.4: for a positive semidefinite A, detA≤a11a22⋯ann, with equality if and only if A is diagonal or some aii=0. That proof inducts on the size through Fischer’s inequality (Corollary 13.6.3). A much shorter route scales A so that its diagonal is flat, and then compares the product of the eigenvalues with their sum. That comparison is the inequality between the arithmetic and geometric means, which the book has not yet proved; it is pure algebra, and it comes first.

Lemma 17.8.6 (The Arithmetic–Geometric Mean Inequality)

Let n≥1 and let t1,…,tn be non-negative real numbers. Then t1t2⋯tn≤(t1+t2+⋯+tnn)n, with equality if and only if t1=t2=⋯=tn.

Idea. Scale so that the mean is 1. If the numbers are not all 1, one of them lies above 1 and another below. Merge that pair into the single number tj+tk−1. This leaves n−1 non-negative numbers with mean 1, so induction applies to them. The merge can only increase the product, because (tj+tk−1)−tjtk=(tj−1)(1−tk), and both factors are positive. So the original product is at most, and when it is non-zero strictly less than, a product that induction bounds by 1.

Proof. Step 1: mean 1. We prove by induction on n≥1: if t1,…,tn≥0 and t1+⋯+tn=n, then t1⋯tn≤1, with equality only when every ti=1. For n=1 the hypothesis says t1=1, and the product is 1.

Let n≥2 and assume the claim for n−1 numbers. If every ti=1, the product is 1. Otherwise some tj>1: if every ti≤1 with at least one strictly smaller, the sum would be less than n. Similarly some tk<1, since if every ti≥1 with at least one strictly larger, the sum would exceed n. Then j≠k. Neither the sum nor the product depends on the order of the numbers, so we may relabel so that j=n−1 and k=n. Put s=tn−1+tn−1. Then s>0, since tn−1>1 and tn≥0. The n−1 numbers t1,…,tn−2,s are non-negative and add up to n−1. So by the inductive hypothesis ps≤1, where p=t1⋯tn−2≥0 (with p=1 when n=2). Moreover s−tn−1tn=(tn−1−1)(1−tn)>0. If p=0, then t1⋯tn=0<1. If p>0, then t1⋯tn=ptn−1tn<ps≤1. In either case the product is strictly less than 1 when the ti are not all 1. This completes the induction.

Step 2: the general case. Let m=1n(t1+⋯+tn)≥0. If m=0, then every ti=0, since non-negative numbers with sum 0 are all 0. Both sides are then 0, and the ti are all equal. If m>0, the numbers ti/m are non-negative with sum n. By Step 1, t1m⋯tnm≤1,that is,t1⋯tn≤mn, with equality if and only if every ti/m=1, that is, every ti=m. Numbers that are all equal to some c have mean m=c, so “every ti=m” is the same as “all the ti are equal”. This proves the lemma.

In the language of this section, the lemma says that among non-negative vectors with a fixed total, the flat vector, which is the bottom of the majorization order (Example (A chain, the two extremes, and size one) (b)), has the largest product.

Corollary 17.8.7 (Hadamard’s Inequality, Second Proof)

Let F=R or F=C, and let A∈Mn(F) with A⪰0. Then detA≤a11a22⋯ann, with equality if and only if A is diagonal or some aii=0.

This is Theorem 13.6.4 (a) again, with the same hypothesis and the same equality condition. What is new is the proof.

Idea. Suppose first that every diagonal entry equals 1. Then trA=n, so the eigenvalues, which are non-negative, add up to n, and Lemma 17.8.6 gives detA=∏iλi(A)≤1=∏iaii. A general diagonal is made flat by a diagonal congruence A↦SAS with S=diag(a11−1/2,…,ann−1/2). This congruence preserves positivity and divides both sides of the inequality by the same number a11⋯ann. Zero diagonal entries, where S does not exist, are handled first and separately.

Proof. Each aii=⟨Aei,ei⟩≥0 by Definition 13.1.1, and each λi(A)≥0 by Theorem 13.1.4 (b).

Case 1: some aii=0. The right-hand side is 0. Suppose detA≠0. By Theorem 9.3.3, detA=∏jλj(A), so every λj(A) is non-zero, hence positive. By Theorem 13.1.5 (b), A≻0, so aii=⟨Aei,ei⟩>0, a contradiction. Hence detA=0, and the inequality holds with equality, as the statement says it should when some aii=0.

Case 2: every aii>0. Let S=diag(a11−1/2,…,ann−1/2), a real diagonal matrix with non-zero diagonal, so S∗=S and S is invertible. Put B=S∗AS=SAS. By Proposition 13.1.2 (a), B⪰0. Its entries are bij=aij/aiiajj, so bii=1 for every i, and trB=n. By Theorem 9.3.3, λ1(B)+⋯+λn(B)=n. The λi(B) are non-negative by Theorem 13.1.4 (b). So Lemma 17.8.6 and Theorem 9.3.3 give detB=∏i=1nλi(B)≤(nn)n=1, with equality if and only if all λi(B) are equal, hence all equal to 1. By Theorem 7.5.4, and since detS=(a11⋯ann)−1/2, detB=(detS)2detA=detAa11a22⋯ann. Multiplying detB≤1 by the positive number a11⋯ann gives detA≤a11⋯ann.

For equality in Case 2: equality holds exactly when every λi(B)=1. If so, the spectral theorem (Corollary 12.4.3, or Corollary 12.5.2 over R) gives B=UIU∗=I for some unitary (or orthogonal) U. Then A=S−1BS−1=S−2=diag(a11,…,ann) is diagonal. Conversely, if A is diagonal, then B=SAS=I, and every λi(B)=1. So in Case 2, equality holds if and only if A is diagonal. Together with Case 1, this is the stated equality condition. This proves the corollary.

For the matrix of Example (Coordinate vectors against eigenvectors), which is positive definite since its eigenvalues are 7,4,1, the inequality reads detA=7⋅4⋅1=28≤60=3⋅4⋅5.

The two proofs show different things. The first proves more along the way: Fischer’s inequality for blocks of any sizes. The second shows that once the diagonal is scaled flat, the trace and the arithmetic–geometric mean inequality suffice. Of d(B)≺λ(B) it uses only clause (M2), which for a matrix is just the trace.

Remark. A proof by majorization in earnest would use all of d(A)≺λ(A), without scaling, together with the fact that for vectors with non-negative entries, x≺y implies x1⋯xn≥y1⋯yn (the product is Schur-concave). That fact is true but not proved in this book yet; Chapter 21 proves it, for the whole class of functions that reverse majorization.

Majorization is the language this chapter’s eigenvalue inequalities were already using, and the Schur–Horn theorem shows that it fits the diagonal problem exactly. The same language compresses the inequality of Section 7 for arbitrary index sets, Theorem 17.7.2, into one line. The k largest entries of the vector λ(A)−λ(B) add up to the largest value of ∑j(λij(A)−λij(B)) over index sets of size k, and the totals are trA−trB=tr(A−B). So a bound of every such sum by ∑j≤kλj(A−B) says exactly that λ(A)−λ(B)≺λ(A−B). The diagonal of a Hermitian matrix and the eigenvalues of a difference are both controlled by one ordering of vectors in Rn.

17.8.7. Exercises

A. Check your understanding

Exercise (A1)
  • (a) State the definitions of x≺y and x≺wy for x,y∈Rn.
  • (b) Decide whether (2,1,1)≺(2,2,0), and whether (2,2,0)≺(2,1,1).
  • (c) True or false: if x≺y, then x1≤y1. Justify your answer.
  • (d) State Theorem 17.8.2. Which result of §06 is it a restatement of?
  • (e) State Theorem 17.8.5. Which half is proved by induction on n, and which lemma does each step of that induction use?

Solution.

  • (a) x≺y means that ∑i≤kxi↓≤∑i≤kyi↓ for every k=1,…,n−1, and ∑i=1nxi=∑i=1nyi (Definition 17.8.1). x≺wy means that the same inequality holds for every k=1,…,n, with no condition on the totals.

  • (b) Both vectors are decreasing with total 4. Their running totals at k=1,2 are 2,3 and 2,4. Since 2≤2 and 3≤4, (2,1,1)≺(2,2,0). The reverse fails at k=2, since 4>3.

  • (c) False. The definition compares sorted vectors, so the unsorted first entries can go either way. Take x=(1,0) and y=(0,1). Both have decreasing rearrangement (1,0), so (M1) at k=1 reads 1≤1, and both totals are 1; hence x≺y. But x1=1>0=y1. What is true is x1↓≤y1↓, which is (M1) at k=1.

  • (d) For a Hermitian A∈Mn(F), F=R or C, d(A)≺λ(A). It restates Corollary 17.6.3: applied to the positions of the k largest diagonal entries, the corollary gives (M1), and its equality at k=n gives (M2).

  • (e) For d,λ∈Rn, there is a Hermitian matrix over R or C with diagonal d, in the given order, and eigenvalues λ1,…,λn if and only if d≺λ. The half “d≺λ implies such a matrix exists” is Horn’s theorem, Theorem 17.8.4, proved by induction on n. Each step uses Lemma 17.8.3, applied to a pair of consecutive eigenvalues λk≥d1≥λk+1, to put d1 in the corner. The other half is Schur’s theorem.

B. Practice

Exercise (B1: Determine which pairs are comparable)

For each pair, determine whether x≺y and whether x≺wy. Where a relation fails, name the exact clause and the value of k at which it fails. Justify your answer.

  • (a) x=(2,2,1,1), y=(3,2,1,0).
  • (b) x=(0,3,1), y=(4,0,0).
  • (c) x=(1,1,0), y=(2,1,0).
  • (d) x=(−1,2,2), y=(3,0,0).
  • (e) x=(3,1,1,1), y=(2,2,2,0); decide also whether y≺x.

Solution. In each case, sort first, then compare running totals.

  • (a) Both are decreasing with total 6. Running totals: 2,4,5,6 and 3,5,6,6. Every entry of the first list is at most the matching entry of the second, so x≺y and hence x≺wy.

  • (b) x↓=(3,1,0) with running totals 3,4,4, and y has running totals 4,4,4. The totals agree and 3≤4, 4≤4, so x≺y and x≺wy. This is Example (A chain, the two extremes, and size one) (c), applied after sorting.

  • (c) Running totals 1,2,2 and 2,3,3. Every inequality holds, including at k=3, so x≺wy. But the totals are 2≠3, so (M2) fails and x⊀y.

  • (d) x↓=(2,2,−1) with running totals 2,4,3, and y has running totals 3,3,3. The totals agree, but at k=2, 4>3, so (M1) fails at k=2. Neither x≺y nor x≺wy holds.

  • (e) Both are decreasing with total 6. Running totals: 3,4,5,6 and 2,4,6,6. At k=1, 3>2, so (M1) fails at k=1 and x⊀y, and also x⊀wy. For y≺x, the inequalities 2≤3 and 4≤4 hold, but 6>5 at k=3, so (M1) fails at k=3. The pair is incomparable, like Example (Neither majorizes the other).

Exercise (B2: Schur’s theorem for one matrix)

Let A=(322240207).

  • (a) Verify that (2,1,4), (−1,−2,1) and (−3,2,1) are eigenvectors of A, and find λ(A).
  • (b) Write down d(A) and verify every clause of d(A)≺λ(A).
  • (c) Hence verify Hadamard’s inequality (Theorem 13.6.4) for A.

Solution.

  • (a) Multiply out: A(2,1,4)=(6+2+8, 4+4, 4+28)=(16,8,32)=8(2,1,4),A(−1,−2,1)=(−3−4+2, −2−8, −2+7)=(−5,−10,5)=5(−1,−2,1),A(−3,2,1)=(−9+4+2, −6+8, −6+7)=(−3,2,1)=1(−3,2,1). Three distinct eigenvalues of a 3×3 matrix account for the whole spectrum, so λ(A)=(8,5,1).

  • (b) d(A)=(3,4,7), and d(A)↓=(7,4,3). (M1): at k=1, 7≤8; at k=2, 7+4=11≤13=8+5. (M2): 3+4+7=14=8+5+1. Hence d(A)≺λ(A), as Theorem 17.8.2 requires.

  • (c) All eigenvalues are positive, so A≻0 by Theorem 13.1.5 (b), and Corollary 17.8.7 applies. By Theorem 9.3.3, detA=8⋅5⋅1=40, while a11a22a33=3⋅4⋅7=84. Hence 40≤84. The inequality is strict, as it must be, since A is not diagonal and no diagonal entry is 0.

Exercise (B3: Which diagonals are ruled out?)

A Hermitian 3×3 matrix has eigenvalues 5,1,0. For each proposed list of diagonal entries, determine whether Theorem 17.8.2 rules it out. Justify your answer.

  • (a) 4,2,0.
  • (b) 6,0,0.
  • (c) 3,3,1.
  • (d) 2,2,2.
  • (e) −1,4,3.

For the lists that are not ruled out, what does Theorem 17.8.5 add?

Solution. Here λ=(5,1,0), with running totals 5,6 and total 6.

  • (a) Running totals 4,6, total 6. Since 4≤5 and 6≤6, the list is not ruled out.

  • (b) At k=1, 6>5. Ruled out. No diagonal entry can exceed λ1.

  • (c) The total is 7≠6. Ruled out by (M2): the trace would be wrong.

  • (d) This is the flat vector with total 6, which is majorized by every vector with total 6 (Example (A chain, the two extremes, and size one) (b)). Not ruled out.

  • (e) Sorted, it is (4,3,−1), with running totals 4,7 and total 6. At k=2, 7>6. Ruled out. Independently, −1<0=λ3 contradicts the bottom-end inequality of the warning.

Schur’s theorem is only a necessary condition, so for (a) and (d) it shows only that nothing forbids these diagonals. Theorem 17.8.5 adds that they occur: for each of them there is a real symmetric matrix with eigenvalues 5,1,0 and that diagonal.

C. Going deeper

Exercise (C1: Majorization squeezes the range)
  • (a) Prove that if x≺y, then yn↓≤xn↓≤x1↓≤y1↓.
  • (b) Deduce that every diagonal entry of a Hermitian A lies in [λn(A),λ1(A)].
  • (c) Show that the conclusion of (b), together with equal totals, does not imply majorization. Prove that no Hermitian matrix has eigenvalues 4,1,1,0 and diagonal entries 3,3,0,0, although every one of these entries lies in [0,4] and the totals agree.

Solution.

  • (a) The inequality x1↓≤y1↓ is (M1) at k=1 (for n=1, x=y by Example (A chain, the two extremes, and size one) (d), and there is nothing to prove). The middle inequality holds because x↓ is decreasing. For the left inequality, subtract (M1) at k=n−1 from (M2): xn↓=∑i=1nxi−∑i=1n−1xi↓≥∑i=1nyi−∑i=1n−1yi↓=yn↓.

  • (b) By Theorem 17.8.2, d(A)≺λ(A). Apply (a) with x=d(A) and y=λ(A), using λ(A)↓=λ(A). Every diagonal entry lies between the smallest and the largest diagonal entry, which (a) places between λn(A) and λ1(A). (This also follows from Lemma 13.5.2 at x=ei. The point of (a) is that it is a property of majorization itself.)

  • (c) The vector (3,3,0,0) has running totals 3,6,6, and (4,1,1,0) has 4,5,6. At k=2, 6>5, so (3,3,0,0)⊀(4,1,1,0), and by Theorem 17.8.2 no Hermitian matrix has this diagonal and this spectrum. Yet 0≤3,3,0,0≤4, and both totals are 6. The range condition controls only the running totals at k=1 and k=n−1, while majorization constrains the middle ones too, which is why the example needs n≥4.

Exercise (C2: A weaker determinant bound)

Let A∈Mn(F) with A⪰0.

  • (a) Prove that detA≤(1ntrA)n, with equality if and only if A=cI for some real c≥0.
  • (b) Prove that a11⋯ann≤(1ntrA)n. Deduce that Hadamard’s inequality (Theorem 13.6.4) always implies the bound in (a).
  • (c) For A=(202052226), whose eigenvectors are (1,2,3), (1,−2,1) and (−4,−1,2), compute the three numbers detA, a11a22a33 and (13trA)3.

Solution.

  • (a) The eigenvalues λi(A) are non-negative by Theorem 13.1.4 (b). By Theorem 9.3.3 and Lemma 17.8.6, detA=∏iλi(A)≤(1n∑iλi(A))n=(1ntrA)n. Equality holds if and only if all λi(A) equal a common value c≥0. By the spectral theorem this happens if and only if A=U(cI)U∗=cI. Conversely, A=cI gives equality directly.

  • (b) The aii are non-negative (as in the proof of Corollary 17.8.7) and add up to trA. So Lemma 17.8.6 gives ∏iaii≤(1ntrA)n. Combined with Corollary 17.8.7, this gives detA≤∏iaii≤(1ntrA)n, so the Hadamard bound always sits between detA and the bound of (a). The bound in (a) sees only the total of the diagonal, while Hadamard’s bound sees each entry, and so it is never worse.

  • (c) Check the eigenvectors: A(1,2,3)=(2+6, 10+6, 2+4+18)=(8,16,24); A(1,−2,1)=(2+2, −10+2, 2−4+6)=(4,−8,4); A(−4,−1,2)=(−8+4, −5+4, −8−2+12)=(−4,−1,2). So λ(A)=(8,4,1), and A≻0. Then detA=8⋅4⋅1=32,a11a22a33=2⋅5⋅6=60,(trA3)3=(133)3=219727≈81.4, and 32≤60≤81.4, with both inequalities strict.

Exercise (C3: The possible diagonals form a convex set)

For x∈Rn and 1≤k≤n, write sk(x)=x1↓+⋯+xk↓.

  • (a) Prove that sk(x)=max{∑i∈Ixi:I⊆{1,…,n}, |I|=k}.
  • (b) Deduce that sk(x+y)≤sk(x)+sk(y) for all x,y∈Rn.
  • (c) Let λ∈Rn. Prove that if d≺λ and d′≺λ, then td+(1−t)d′≺λ for every 0≤t≤1.
  • (d) Hence show that if A and B are Hermitian n×n matrices with the same eigenvalues, then for every 0≤t≤1 there is a real symmetric matrix with those eigenvalues and diagonal td(A)+(1−t)d(B).

Hint for (a): if v1≥⋯≥vk are the entries xi, i∈I, in decreasing order, count how many entries of x are at least vj.

Solution.

  • (a) Let |I|=k, and list the entries xi, i∈I, decreasingly as v1≥⋯≥vk. For each j≤k, the entries v1,…,vj sit at j different positions of x and are all ≥vj. If xj↓<vj, then, since x↓ is decreasing, only x1↓,…,xj−1↓ could be ≥vj, which is fewer than j entries. Hence vj≤xj↓ for every j, and summing gives ∑i∈Ixi≤sk(x). The bound is attained by the positions of the k largest entries, chosen as in the proof of Theorem 17.8.2. This proves (a).

  • (b) Let I be a set of k positions at which x+y attains the maximum in (a). Applying (a) three times, sk(x+y)=∑i∈Ixi+∑i∈Iyi≤sk(x)+sk(y).

  • (c) Put z=td+(1−t)d′. Its total is t∑iλi+(1−t)∑iλi=∑iλi, by (M2) for both hypotheses. For a scalar c≥0, multiplying by c does not change the order of the entries, so sk(cx)=csk(x). For 1≤k≤n−1, (b) and (M1) for both hypotheses give sk(z)≤tsk(d)+(1−t)sk(d′)≤tsk(λ)+(1−t)sk(λ)=sk(λ). So z≺λ.

  • (d) Let λ=λ(A)=λ(B). By Theorem 17.8.2, d(A)≺λ and d(B)≺λ. By (c), td(A)+(1−t)d(B)≺λ. By Theorem 17.8.4 there is a real symmetric matrix with eigenvalues λ and exactly this diagonal. So among Hermitian matrices with a fixed spectrum, the set of diagonals contains the segment between any two of its points: it is convex.

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.