17.4  Deleting a Row and a Column

Take a Hermitian matrix, cross out the last row and the last column, and look at what is left. It is again Hermitian, one size smaller, and it has its own eigenvalues. Two lists of numbers are now in play, n of them and n−1 of them, and they cannot be independent of each other: the smaller matrix is the larger one restricted to a hyperplane, and Theorem 17.2.3 describes every eigenvalue by what happens on subspaces. Restricting the supply of subspaces can only make the description worse in one direction and better in the other, and the exact bookkeeping is this section’s subject.

The answer is as tight as it could be. The n−1 eigenvalues of the corner fall into the n−1 gaps of the eigenvalues of the whole, one to a gap. Nothing stronger holds: every interlacing pattern actually occurs (Exercise C2 proves the case n=2).

Throughout this section the field is R or C, and every matrix whose eigenvalues are written down is Hermitian, with eigenvalues real and indexed decreasingly. For an index set I⊆{1,…,n} we write AI,I for the principal submatrix of A on the rows and columns indexed by I, as in Definition 7.7.1. A principal submatrix of a Hermitian matrix is Hermitian, since (AI,I)∗=(A∗)I,I=AI,I.

17.4.1. A matrix restricted to a coordinate hyperplane

Deleting row m and column m is not an arbitrary mutilation of A. It is what A looks like to vectors whose m-th coordinate is zero, and that is the only observation the proof needs.

Lemma 17.4.1 (Deleting a Row and a Column Is Restricting to a Hyperplane)

Let A∈Mn(F) be Hermitian with n≥2, fix m with 1≤m≤n, put I={1,…,n}∖{m} and B=AI,I∈Mn−1(F). Let H={x∈Fn:xm=0}, and let φ:H→Fn−1 delete the m-th coordinate. Then φ is a linear isomorphism, and for every x∈H, ⟨Ax,x⟩=⟨Bφ(x),φ(x)⟩,‖x‖=‖φ(x)‖. In particular RA(x)=RB(φ(x)) for every x∈H with x≠0.

Proof. The map φ is linear, being given by deleting a coordinate, and it is bijective from H to Fn−1 because inserting 0 in position m inverts it. Write I={i1<⋯<in−1}, so that φ(x)r=xir and (B)rs=airis by Definition 7.7.1.

Let x∈H. Expanding the inner product of Fn, which is linear in the first slot and conjugate-linear in the second, ⟨Ax,x⟩=∑p=1n∑q=1napqxqxp―. Every term with p=m or q=m vanishes, because xm=0. The remaining terms are indexed by p,q∈I, say p=ir and q=is, and there the summand is airisφ(x)sφ(x)r―=(B)rsφ(x)sφ(x)r―. Summing over r,s gives ⟨Bφ(x),φ(x)⟩. The same deletion of a zero coordinate gives ‖x‖2=∑p≠m|xp|2=‖φ(x)‖2. Dividing the first identity by the second, which is positive when x≠0, gives the statement about the Rayleigh quotients (Definition 17.1.1). This proves the lemma.

Because φ is an isomorphism, it matches the k-dimensional subspaces of H with the k-dimensional subspaces of Fn−1, for every k. So a statement about RB on subspaces of Fn−1 is a statement about RA on subspaces of H, and the only thing lost in passing from A to B is the subspaces of Fn that are not inside H.

17.4.2. Interlacing

Theorem 17.4.2 (Cauchy’s Interlacing Theorem)

Let A∈Mn(F) be Hermitian with n≥2, let I⊆{1,…,n} with |I|=n−1, and let B=AI,I. Then λk(A) ≥ λk(B) ≥ λk+1(A)(1≤k≤n−1), that is, λ1(A)≥λ1(B)≥λ2(A)≥λ2(B)≥⋯≥λn−1(B)≥λn(A).

Idea. Both halves of Theorem 17.2.3 are used, one for each inequality, and each time the move is the same: a subspace of H is also a subspace of Fn.

① For λk(B)≤λk(A) use the max–min form. The best k-dimensional subspace for B lives inside H; transplanted to Fn it is one competitor among more, so the maximum over all of them is at least as large.

② For λk+1(A)≤λk(B) use the min–max form, and count dimensions. The best subspace for B there has dimension (n−1)−k+1=n−k, and n−k=n−(k+1)+1 is exactly the dimension the min–max form for λk+1(A) asks for. So the same subspace, viewed in Fn, is admissible for A at the index k+1.

The dimension count that makes min–max work at all is Lemma 17.2.1, and it has already been spent inside Theorem 17.2.3; nothing more is needed here than the coincidence of the two numbers in ②.

Proof. Let m be the index not in I, and let H and φ be as in Lemma 17.4.1, so that RA(x)=RB(φ(x)) for every non-zero x∈H. Both A and B are Hermitian, so Theorem 17.2.3 applies to each, to A with n and to B with n−1 in the role of the size. Fix k with 1≤k≤n−1.

(λk(B)≤λk(A)) By the max–min form of Theorem 17.2.3 for B, there is a subspace W′≤Fn−1 with dim⁡W′=k and λk(B)=min0≠y∈W′RB(y). Put W=φ−1(W′)≤H≤Fn. Since φ is an isomorphism, dim⁡W=k, and φ carries the non-zero vectors of W onto the non-zero vectors of W′, so min0≠x∈WRA(x)=min0≠y∈W′RB(y)=λk(B). Now W is one of the k-dimensional subspaces of Fn over which the max–min form for A takes its maximum, so λk(A)≥λk(B).

(λk+1(A)≤λk(B)) By the min–max form of Theorem 17.2.3 for B, whose matrices have size n−1, there is a subspace W′≤Fn−1 with dim⁡W′=(n−1)−k+1=n−k and λk(B)=max0≠y∈W′RB(y). Put W=φ−1(W′)≤Fn, of dimension n−k. Since 1≤k+1≤n and n−(k+1)+1=n−k=dim⁡W, the subspace W is admissible in the min–max form for λk+1(A), and that form takes the minimum over all such subspaces. Hence λk+1(A) ≤ max0≠x∈WRA(x)=max0≠y∈W′RB(y)=λk(B). This proves the theorem.

The picture to keep is one of nesting. The interval [λn(A),λ1(A)] contains the interval [λn−1(B),λ1(B)], and inside it the two lists alternate, an eigenvalue of B between each consecutive pair of eigenvalues of A. Deleting a row and a column can only shrink the spectrum inward, never spread it out.

Quick check. Both inequalities in Theorem 17.4.2 can be equalities. Give a matrix A for which λ1(A)=λ1(B), and one for which λ1(B)=λ2(A).

Solution. For the first, A=I2 with I={1}: then B=(1) and λ1(A)=λ1(B)=1. For the second, A=diag(1,0) with I={2}: then B=(0), and λ1(B)=0=λ2(A). Nothing in the theorem forces a gap; when a gap is forced is exactly the question settled later in this section.

Deleting one index at a time and iterating gives the general statement.

Corollary 17.4.3 (Deleting m Rows and the Same m Columns)

Let A∈Mn(F) be Hermitian, let 0≤m≤n−1, let I⊆{1,…,n} with |I|=n−m, and let B=AI,I. Then λi(A) ≥ λi(B) ≥ λi+m(A)(1≤i≤n−m).

Proof. Induction on m. For m=0 we have B=A and both inequalities are equalities.

Let m≥1 and suppose the statement holds for m−1, for every Hermitian matrix. Choose an index j∉I and put I′=I∪{j}, so |I′|=n−m+1, and B′=AI′,I′∈Mn−m+1(F), which is Hermitian. Since I⊆I′, taking the principal submatrix on I of B′ means keeping the rows and columns of A indexed by I, so (B′)I,I=B after relabeling the indices of I inside I′; and |I|=|I′|−1.

Fix i with 1≤i≤n−m. By the inductive hypothesis applied to A and B′, λi(A) ≥ λi(B′) ≥ λi+m−1(A), which is available because i≤n−m<n−m+1. By Theorem 17.4.2 applied to B′ and its principal submatrix B, which is one size smaller, λk(B′) ≥ λk(B) ≥ λk+1(B′)(1≤k≤n−m). Combining, λi(B)≤λi(B′)≤λi(A), which is the first inequality. For the second, take k=i in the display above and then use the inductive hypothesis at the index i+1, legitimate since i+1≤n−m+1: λi(B) ≥ λi+1(B′) ≥ λ(i+1)+m−1(A)=λi+m(A). This proves the corollary.

The extreme case m=n−1 is worth writing down, because it recovers a fact about diagonal entries that also follows from Chapter 13’s description of the extreme eigenvalues.

Example (The diagonal lies inside the spectrum)

Let A∈Mn(F) be Hermitian. Show that λn(A)≤aii≤λ1(A) for every i.

Solution. Take I={i}, so m=n−1 and B=AI,I=(aii), a 1×1 Hermitian matrix whose single eigenvalue is aii — real, as a diagonal entry of a Hermitian matrix must be. Corollary 17.4.3 at i=1 gives λ1(A) ≥ λ1(B)=aii ≥ λ1+(n−1)(A)=λn(A). The same conclusion follows from Lemma 13.5.2 applied to the unit vector ei, since ⟨Aei,ei⟩=aii. What interlacing adds is the version for larger index sets: any k diagonal entries are the diagonal of a k×k principal submatrix, whose i-th eigenvalue is trapped between λi+n−k(A) and λi(A). How much more can be said about the diagonal of a Hermitian matrix, given its spectrum, is the subject of Section 8.

17.4.3. When is the interlacing strict?

Theorem 17.4.2 allows equalities, and the Quick check produced them. Deciding when they occur is the delicate part of the subject, and it is easiest to see when the deleted index is the last one, so that A is presented as B with a border added.

Theorem 17.4.4 (Bordering a Hermitian Matrix)

Let B∈Mn−1(F) be Hermitian with n≥2, let b∈Fn−1, let c∈R, and let A=(Bbb∗c)∈Mn(F), which is Hermitian and has AI,I=B for I={1,…,n−1}. For an eigenvalue ν of B, write Pν for the orthogonal projection of Fn−1 onto the eigenspace Eν(B), and mν=dim⁡Eν(B).

  • (a) If Pνb≠0, then ν is a root of pA of multiplicity exactly mν−1; in particular, if mν=1 then ν is not an eigenvalue of A.
  • (b) If Pνb=0, then ν is a root of pA of multiplicity at least mν.
  • (c) The interlacing λ1(A)≥λ1(B)≥⋯≥λn−1(B)≥λn(A) is strict at every one of its 2(n−1) inequalities if and only if B has n−1 distinct eigenvalues and Pνb≠0 for every eigenvalue ν of B.

Idea. Write the characteristic polynomial of A down. With B diagonalized, a Schur complement turns det(xI−A) into (∏k(x−μk))((x−c)−∑k|b~k|2x−μk), and clearing the denominators leaves a polynomial identity in which each μk appears with one factor fewer than one might expect — provided Pνb≠0, that is, provided the |b~k|2 over the group belonging to that eigenvalue do not all vanish (for a repeated eigenvalue the individual b~k depend on the choice of basis, but their sum does not). Part (c) then costs almost nothing: strict interlacing means the combined list of 2n−1 numbers is strictly decreasing, hence has no repetitions at all, and (a) and (b) say exactly which repetitions are possible.

Proof. Step 1: diagonalize the corner. By the spectral theorem (Corollary 12.4.3 over C, Corollary 12.5.2 over R) write B=UDU∗ with U∈Mn−1(F) unitary and D=diag(μ1,…,μn−1) carrying the eigenvalues of B in decreasing order, the columns u1,…,un−1 of U being an orthonormal basis of eigenvectors. Let V=(U00∗1)∈Mn(F), which is unitary because V∗V=In. Multiplying in blocks (Theorem 8.1.1) gives V∗AV=(Db~b~∗c),b~=U∗b. Similar matrices have the same characteristic polynomial (Theorem 7.10.3), so we may compute with the right-hand side. Moreover b~k=(U∗b)k=⟨b,uk⟩, and the uk with μk=ν form an orthonormal basis of Eν(B): they lie in it, and a vector x=∑kckuk with Bx=νx has (μk−ν)ck=0 for every k, so ck=0 unless μk=ν. Hence exactly mν of the μk equal ν, and Pνb=∑k:μk=ν⟨b,uk⟩ukand‖Pνb‖2=∑k:μk=ν|b~k|2. In particular Pνb=0 exactly when b~k=0 for every k with μk=ν.

Step 2: the characteristic polynomial. We claim that, as polynomials in x,

pA(x)=(x−c)∏k=1n−1(x−μk)−∑k=1n−1|b~k|2∏l≠k(x−μl).
(17.4.1)

Claim. Eq. 17.4.1 holds.

Proof. Fix a scalar x0 that is not an eigenvalue of B. Then x0In−1−D is invertible, being diagonal with non-zero diagonal entries, so Theorem 8.3.2 (a), applied to M=x0In−V∗AV=(x0In−1−D−b~−b~∗x0−c) with the invertible block in the top left, gives detM=det(x0I−D)⋅det(M/(x0I−D)), where the Schur complement is the 1×1 matrix (x0−c)−(−b~∗)(x0I−D)−1(−b~)=(x0−c)−∑k|b~k|2x0−μk. Multiplying out det(x0I−D)=∏k(x0−μk) against this scalar, and canceling one factor (x0−μk) in the k-th summand, gives exactly the right-hand side of Eq. 17.4.1 evaluated at x0. So the two sides of Eq. 17.4.1, which are polynomials of degree n, agree at every scalar outside the finite set spec(B). Their difference is therefore a polynomial of degree at most n with infinitely many roots in F, and a non-zero polynomial of degree at most n has at most n roots (Corollary 6.1.4); so the difference is zero.

Step 3: multiplicities. Let ν1>⋯>νp be the distinct eigenvalues of B, with multiplicities m1,…,mp summing to n−1, and put βt=‖Pνtb‖2. Grouping equal μk in Eq. 17.4.1 — the product ∏l≠k(x−μl) is the same for all k in one group, and the coefficients |b~k|2 add up to βt by Step 1 — turns it into pA(x)=(x−c)∏t=1p(x−νt)mt−∑t=1pβt(x−νt)mt−1∏s≠t(x−νs)ms. Fix t. Every term on the right is divisible by (x−νt)mt−1, so pA(x)=(x−νt)mt−1g(x) with g(x)=(x−c)(x−νt)∏s≠t(x−νs)ms−βt∏s≠t(x−νs)ms−∑u≠tβu(x−νt)(x−νu)mu−1∏s≠u,t(x−νs)ms. Evaluate at x=νt. The first term vanishes because of its factor (x−νt), and so does every summand of the third, for the same reason. What is left is g(νt)=−βt∏s≠t(νt−νs)ms, and the product is non-zero because the νs are distinct. So g(νt)≠0 exactly when βt≠0, that is, exactly when Pνtb≠0. In that case νt is a root of pA of multiplicity exactly mt−1, which proves (a); in the other case (x−νt) divides g, so the multiplicity is at least mt, which proves (b).

Step 4: strictness. (⇐) Suppose B has n−1 distinct eigenvalues, so every mt=1, and Pνb≠0 for each of them. By (a), no eigenvalue of B is an eigenvalue of A. Theorem 17.4.2 gives λk(A)≥λk(B)≥λk+1(A) for each k; since λk(B) is an eigenvalue of B and λk(A), λk+1(A) are eigenvalues of A, neither inequality can be an equality. So all 2(n−1) of them are strict.

(⇒) Suppose all of them are strict, that is, λ1(A)>λ1(B)>λ2(A)>⋯>λn−1(B)>λn(A). Reading two steps at a time gives λk(B)>λk+1(A)>λk+1(B) for 1≤k≤n−2, so the eigenvalues of B are distinct. And the displayed chain is strictly decreasing, so its 2n−1 entries are pairwise distinct; in particular no eigenvalue of B is an eigenvalue of A. By (b), read as a contrapositive, Pνb≠0 for every eigenvalue ν of B. This proves the theorem.

Remark. Part (b) is deliberately one-sided. When Pνb=0 the multiplicity of ν in pA may be mν or mν+1 — never more, since by Cauchy interlacing q equal eigenvalues ν of A squeeze q−1 eigenvalues of B to ν — and both occur: B=diag(3,1), b=(0,1), c=0 gives pA(x)=(x−3)(x2−x−1), where 3 keeps multiplicity 1=m3; while B=(μ), b=0, c=μ gives A=μI2, where μ has multiplicity 2=mμ+1. Deciding between the two needs the value of c as well, and Theorem 17.4.4 does not attempt it; what it does settle completely is the strictness question (c), which is what interlacing asks.

Warning. Strictness needs both clauses, and a non-zero border is not enough. It is tempting to say that the interlacing is strict as soon as b≠0, or as soon as b meets every eigenspace of B. The first is wrong even when the eigenvalues of B are distinct (Exercise B3(a)); the second fails exactly when B has a repeated eigenvalue. Take B=2I2, b=(1,1) and c=2. Here b has a non-zero component in the one eigenspace E2(B)=F2, and yet λ1(B)=λ2(B)=2, so the middle of the interlacing chain reads 2≥λ2(A)≥2 and λ2(A)=2 with no freedom at all. A repeated eigenvalue of B is passed on to A, by Theorem 17.4.4 (a) with mν=2.

To apply Theorem 17.4.4 when the deleted index is not the last one, permute. Let I={i1<⋯<in−1}={1,…,n}∖{m}, and let P be the matrix with columns ei1,…,ein−1,em, a permutation matrix. Its columns are orthonormal, so P⊤P=I and P is unitary (Lemma 3.4.7 (b)); hence P⊤AP is Hermitian and similar to A, with the same eigenvalues. Its (r,s) entry is eir⊤Aeis=airis for r,s≤n−1, so its leading (n−1)×(n−1) corner is AI,I, and its border is the column (ai1m,…,ain−1m) with corner entry amm. So P⊤AP is a bordered matrix with corner AI,I, and Theorem 17.4.4 applies to it verbatim.

The strictness criterion pays a debt from Chapter 11. The solution of Exercise (C1: Interlacing roots) proved by hand that the roots of consecutive orthogonal polynomials interlace for degrees 1,2,3, and then said: “The general case is not proved this way. The coefficients ak and bk assemble into a symmetric tridiagonal matrix whose characteristic polynomial is pk, and the interlacing then follows from the interlacing theorem for the eigenvalues of a symmetric matrix and its leading submatrices, which is Chapter 17’s business.” Two things in that sentence were asserted, not proved: that the characteristic polynomial is pk, and that the interlacing is strict, which is what “exactly one root strictly between” demands. Plain Theorem 17.4.2 gives only ≥. Strictness needs Theorem 17.4.4 (c), and (c) needs one fact about tridiagonal matrices.

Lemma 17.4.5 (Eigenvectors of an Unreduced Tridiagonal Matrix)

Let T∈Mk(R) be symmetric and tridiagonal, with every entry ti,i+1=ti+1,i (1≤i≤k−1) non-zero. Then every eigenvector u of T has uk≠0. Consequently every eigenspace of T has dimension 1, so T has k distinct eigenvalues.

Proof. Let Tu=νu and suppose uk=0; we show u=0. If k=1 this is immediate. Otherwise row k of Tu=νu reads tk,k−1uk−1+tkkuk=νuk, so tk,k−1uk−1=0 and uk−1=0, since tk,k−1≠0. Now suppose ui=ui+1=0 for some 2≤i≤k−1. Row i reads ti,i−1ui−1+tiiui+ti,i+1ui+1=νui, so ti,i−1ui−1=0 and ui−1=0. Walking up from i=k−1 to i=2 gives uk−2=⋯=u1=0, so u=0, and u is not an eigenvector.

If some eigenspace Eν(T) contained two independent vectors u,u′, then uk′u−uku′ would be a non-zero vector of Eν(T) — non-zero by independence, since uk≠0 — with last coordinate 0, which we have just excluded. So each eigenspace is a line. Since T is real symmetric, its eigenspaces have dimensions adding up to k (Corollary 12.5.2), so there are k of them, that is, k distinct eigenvalues. This proves the lemma.

Corollary 17.4.6 (Zeros of Orthogonal Polynomials Interlace) Optional: nothing later depends on this.

Let (pk) be the monic orthogonal polynomial sequence of an inner product on R[x] to which Theorem 11.10.3 applies, with coefficients ak and bk>0; for the node inner product on R[x]≤n read everything below for k+1≤n. For k≥1 let Tk (the Jacobi matrix of the sequence) be the k×k real symmetric tridiagonal matrix with diagonal a0,a1,…,ak−1 and entries b1,…,bk−1 just above and just below it.

  • (a) pTk=pk; so the roots of pk are the eigenvalues of Tk, and they are real and distinct.
  • (b) If r1>⋯>rk+1 are the roots of pk+1 and s1>⋯>sk those of pk, then r1>s1>r2>s2>⋯>sk>rk+1.

Idea. The polynomials pk are the characteristic polynomials of a nested family of symmetric tridiagonal matrices, each the leading corner of the next. So their zeros are eigenvalues of a matrix and of its corner, and Cauchy interlacing applies at once — non-strictly. Strictness is the bordered-matrix criterion, and both of its clauses come from one fact about tridiagonal matrices with non-zero off-diagonal entries: an eigenvector can never end in a zero.

Proof.

  • (a) Put Dk(x)=det(xIk−Tk) and D0=1. For each real x0, the matrix x0Ik−Tk is tridiagonal in the sense of Chapter 7, with diagonal entries x0−ai−1 and off-diagonal entries −bi in both positions (i,i+1) and (i+1,i), so Theorem 7.8.3 gives D1(x0)=x0−a0 and Dk+1(x0)=(x0−ak)Dk(x0)−(−bk)2Dk−1(x0)=(x0−ak)Dk(x0)−bkDk−1(x0) for k≥1. Here Dk(x0)=det(x0I−Tk) is the characteristic polynomial evaluated at x0 (Lemma 7.10.1), and two real polynomials that agree at every real number are equal (Theorem 6.1.5 (b)), so these identities hold as polynomials. They are the recurrence of Theorem 11.10.3 with the same starting values p0=1, p1=x−a0, so Dk=pk for every k by induction. The roots of pTk are the eigenvalues of Tk (Theorem 9.2.1), which are real because Tk is real symmetric (Theorem 12.2.1) and distinct by Lemma 17.4.5, whose hypothesis holds because every bi>0.

  • (b) Write Tk+1 in bordered form: Tk+1=(Tkbkekbkek⊤ak),ek∈Rk. The border of Tk+1 is w=bkek, whose k-th entry is bk. We call it w rather than b because bk already names a recurrence coefficient. We check the two clauses of Theorem 17.4.4 (c) for B=Tk, with w in place of b. By Lemma 17.4.5, Tk has k distinct eigenvalues. For each eigenvalue ν, the eigenspace Eν(Tk) is spanned by one unit vector u, so Pνw=⟨w,u⟩u=bkuku, which is non-zero because bk>0 and uk≠0 by the same lemma. So Theorem 17.4.4 (c) applies, and the eigenvalues of Tk+1 and of Tk interlace strictly. By (a) these are the roots of pk+1 and of pk. This proves the corollary.

Part (a) gives a second proof, with no integrals, that the roots of pk are real and simple — half of Theorem 11.10.5. The other half, that they lie inside the interval of the weight, is not visible in the matrix and still needs the argument of Chapter 11. For the monic Legendre polynomials of Example (The monic Legendre polynomials, continued), ak=0, b1=13 and b2=415, so T3=(01/301/302/1502/150), and expanding along the first row, pT3(x)=x3−(13+415)x=x3−35x=p3. Its leading corner T2 has characteristic polynomial x2−13=p2 and unit eigenvectors 12(1,±1), both with non-zero last coordinate, so the border 215e2 meets both eigenspaces. The prediction is strict interlacing, and indeed 3/5>1/3>0>−1/3>−3/5, that is 0.7746>0.5774>0>−0.5774>−0.7746, where the strict comparisons 3/5>1/3 and 1/3>0 need no decimals.

17.4.4. A worked example

Example (One matrix and all three of its corners)

Let A=(210121012). Compute the eigenvalues of A and of each of its three 2×2 principal submatrices, exhibit the interlacing in each case, and explain the pattern of strict and non-strict inequalities using Theorem 17.4.4.

Solution. The eigenvalues of A. Expanding det(xI−A) along the first row, pA(x)=(x−2)[(x−2)2−1]+1⋅[−(x−2)]=(x−2)[(x−2)2−2], so λ(A)=(2+2, 2, 2−2)≈(3.41421, 2, 0.58579).

The three corners. With I the index set kept, A{1,2},{1,2}=(2112),A{1,3},{1,3}=(2002),A{2,3},{2,3}=(2112). A matrix (2112) has trace 4 and determinant 3, so its eigenvalues are 3 and 1; and 2I2 has the eigenvalue 2 twice.

The interlacing. In all three cases Theorem 17.4.2 is confirmed: I={1,2}:3.41421>3>2>1>0.58579,I={1,3}:3.41421>2=2=2>0.58579,I={2,3}:3.41421>3>2>1>0.58579. The first and third are strict everywhere; the second has two equalities in the middle. The exact comparisons need no decimals: 2+2>3 is 2>1, and 1>2−2 is 2>1 again.

Why. Read A as a bordered matrix in each case, using the permutation remark above.

For I={1,2} the deleted index is already the last one, so B=(2112), b=(0,1) and c=2. The eigenvalues 3 and 1 of B are distinct, with unit eigenvectors 12(1,1) and 12(1,−1). Both inner products with b are non-zero: ⟨b,12(1,1)⟩=12,⟨b,12(1,−1)⟩=−12. So both clauses of Theorem 17.4.4 (c) hold and the interlacing is strict, as computed. The case I={2,3} is the same: moving index 1 to the end gives the corner (2112) again, with border (a21,a31)=(1,0), whose inner products with the two unit eigenvectors 12(1,±1) are both 12, non-zero.

For I={1,3} the deleted index is 2. Conjugating by the permutation matrix that sends the coordinate order (1,2,3) to (1,3,2) turns A into (201021112),B=2I2,b=(1,1),c=2. Here b≠0 and P2b=b≠0, so the first clause of Theorem 17.4.4 (c) is the one that fails: B has the repeated eigenvalue 2. Part (a) with m2=2 then says that 2 is an eigenvalue of A of multiplicity exactly 1, which is visible in pA, and the interlacing collapses at the middle to 2=2=2. This is exactly the situation of the warning above, occurring inside a perfectly ordinary integer matrix.

17.4.5. Inertia under compression

Interlacing controls not only where the eigenvalues are but how many of them have each sign, and that is a statement about inertia.

Corollary 17.4.7 (A Principal Submatrix Keeps Most of the Inertia)

Let A∈Mn(F) be Hermitian, let 0≤m≤n−1, let |I|=n−m and B=AI,I. Write n+(M) and n−(M) for the number of positive and of negative eigenvalues of a Hermitian M, counted with multiplicity. Then n+(B)≤n+(A)≤n+(B)+m,n−(B)≤n−(A)≤n−(B)+m. In particular A⪰0 implies B⪰0, and A≻0 implies B≻0.

Proof. Write r=n+(B) and s=n+(A), and abbreviate λi=λi(A), μi=λi(B). By Corollary 17.4.3, λi≥μi≥λi+m for 1≤i≤n−m.

(n+(B)≤n+(A)) For i≤r we have μi>0, hence λi≥μi>0. So A has at least r positive eigenvalues.

(n+(A)≤n+(B)+m) If r+m≥n there is nothing to prove, since s≤n. Otherwise r+1≤n−m, so the index i=r+1 is admissible and μr+1≤0 by the definition of r; hence λr+1+m≤μr+1≤0, and the eigenvalues λr+m+1,…,λn are all ≤0 by the decreasing indexing. So s≤r+m.

For the two statements about n−, apply what has just been proved to −A, whose principal submatrix on I is −B. By Lemma 17.3.2 the eigenvalue lists are negated and reversed, so n+(−A)=n−(A) and n+(−B)=n−(B), and the two displayed chains are the same statement.

Finally, A⪰0 means every eigenvalue of A is ≥0 (Theorem 13.1.4 (b)), that is, n−(A)=0; then n−(B)≤n−(A)=0, so B⪰0. And A≻0 means n+(A)=n (Theorem 13.1.5 (b)); then n+(B)≥n+(A)−m=n−m, which is the size of B, so every eigenvalue of B is positive and B≻0. This proves the corollary.

The last two sentences are a second proof of something Chapter 13 already knew. That a principal submatrix of a positive semidefinite matrix is positive semidefinite is the step (a) ⇒ (e) in the proof of Theorem 13.1.4, where it comes from the observation that a principal submatrix is what the quadratic form does to vectors supported on the corresponding coordinates — one line, needing no interlacing. What is new here is the quantitative form. Sylvester’s law of inertia (Theorem 14.5.1) says that the triple of counts is an invariant of a real symmetric matrix up to congruence, and Theorem 14.5.3 identifies Chapter 14’s counts with the eigenvalue counts n+ and n− used here; the corollary says that deleting m rows and the matching m columns can destroy at most m positive and at most m negative eigenvalues, and create none of either. Section 10 gives a second proof of Sylvester’s law from Courant–Fischer, and bounds how far a low-rank perturbation can move the inertia.

17.4.6. Interlacing is a Hermitian phenomenon

Everything above rests on Theorem 17.2.3, which rests on the quadratic form being real. Without the Hermitian hypothesis, nothing survives — not a weakened interlacing, not even a bound.

Example (A corner outside the spectrum)

Let A=(51−8−1). Show that A has real eigenvalues, so that the conclusion of Theorem 17.4.2 can be stated for it, and that the conclusion is false.

Solution. trA=4 and detA=5(−1)−1(−8)=3, so pA(x)=x2−4x+3=(x−3)(x−1) and the eigenvalues are 3 and 1, both real. The two 1×1 principal submatrices are A{1},{1}=(5) and A{2},{2}=(−1), with eigenvalues 5 and −1.

Interlacing would require 3≥5≥1 in the first case and 3≥−1≥1 in the second. Both fail, and they fail on opposite sides: the diagonal entries 5 and −1 lie outside the interval [1,3], which is the convex hull of the spectrum. For a Hermitian matrix that is impossible, by Example (The diagonal lies inside the spectrum).

Remark. The failure is total, not marginal. For any a∈F, the matrix Aa=(a1−(a−1)(a−3)4−a) has trace 4 and determinant a(4−a)+(a−1)(a−3)=4a−a2+a2−4a+3=3, hence the same characteristic polynomial (x−3)(x−1) as above. Its (1,1) principal submatrix is (a). So a general matrix with spectrum {1,3} can have any prescribed number whatsoever as its (1,1) entry: the eigenvalue of a single principal submatrix is subject to no constraint at all. (The two diagonal entries together still add up to the trace, 4; that constraint holds for every matrix and has nothing to do with interlacing.) The choice a=5 recovers A, and a=0 gives (01−34).

Warning. Principal means the same index set on both sides. Interlacing says nothing about a submatrix obtained by deleting row i and column j with i≠j, even for a Hermitian A. Take A=(−1011−10), whose eigenvalues are −9 and −11. Deleting row 1 and column 2 leaves the 1×1 matrix (1), and 1 is nowhere near [−11,−9]. The submatrix in Theorem 17.4.2 is AI,I, and the two I’s are not decoration.

17.4.7. Exercises

A. Check your understanding

Exercise (A1)
  • (a) State Theorem 17.4.2, and say which hypothesis on A and which hypothesis on the submatrix are used.
  • (b) Let A∈M6(C) be Hermitian and let B=AI,I with |I|=4. Which eigenvalues of A bound λ3(B) above and below?
  • (c) Determine whether the following is correct, and justify your answer: if b≠0 then the interlacing between (Bbb∗c) and B is strict.
  • (d) Explain why a principal submatrix of a positive semidefinite matrix is positive semidefinite, in one sentence, without using interlacing.
  • (e) Give a 2×2 matrix with real eigenvalues whose (1,1) entry is larger than both of them.

Solution.

  • (a) See Theorem 17.4.2. The matrix A must be Hermitian, so that Theorem 17.2.3 applies and the eigenvalues are real and can be put in order; and B must be a principal submatrix, B=AI,I with the same index set on rows and columns, so that Lemma 17.4.1 identifies RB with RA restricted to a coordinate hyperplane.

  • (b) Here n=6 and m=2, so Corollary 17.4.3 gives λ3(A)≥λ3(B)≥λ5(A).

  • (c) Incorrect. The warning after Theorem 17.4.4 gives B=2I2, b=(1,1), c=2, where b≠0 and the interlacing is not strict, because B has a repeated eigenvalue. The correct criterion is Theorem 17.4.4 (c), which needs distinct eigenvalues of B as well as Pνb≠0 for each ν.

  • (d) Because ⟨By,y⟩=⟨Ax,x⟩≥0, where x∈Fn is y placed in the coordinates of I and 0 elsewhere (Lemma 17.4.1); this is how Theorem 13.1.4 obtained its clause (e).

  • (e) (51−8−1), with eigenvalues 3 and 1 and (1,1) entry 5 (Example (A corner outside the spectrum)). It is necessarily non-Hermitian, by Example (The diagonal lies inside the spectrum).

B. Practice

Exercise (B1: Interlacing in integers)

Let A=(012102223). Compute the eigenvalues of A and of all three of its 2×2 principal submatrices, and exhibit the interlacing in each case. Hence say, for each corner, whether the interlacing is strict.

Solution. The eigenvalues of A. The vector (1,−1,0) satisfies A(1,−1,0)=(−1,1,0), so −1 is an eigenvalue. Dividing, pA(x)=x3−3x2−9x−5=(x−5)(x+1)2, so λ(A)=(5, −1, −1). (As a check, trA=3=5−1−1, and expanding along the first row detA=0−1(3−4)+2(2−0)=5=5⋅(−1)⋅(−1).)

The corners. A{1,2},{1,2}=(0110),A{1,3},{1,3}=A{2,3},{2,3}=(0223). The first has trace 0 and determinant −1, so its eigenvalues are 1 and −1. The second has trace 3 and determinant −4, so its characteristic polynomial is x2−3x−4=(x−4)(x+1) and its eigenvalues are 4 and −1.

Interlacing. I={1,2}:5>1>−1=−1=−1,I={1,3}:5>4>−1=−1=−1,I={2,3}:5>4>−1=−1=−1. No corner gives strict interlacing, and the reason is on the side of A: whichever index is deleted, λ2(A)=λ3(A)=−1, so the chain λ2(A)≥λ2(B)≥λ3(A) squeezes λ2(B) to −1. Strict interlacing would make the combined chain strictly decreasing, and so would force A to have three distinct eigenvalues. Theorem 17.4.4 sees the same thing from the side of the corner. For I={1,2} the border is b=(2,2), orthogonal to (1,−1), which spans E−1(B); so P−1b=0, and part (b) predicts that −1 is an eigenvalue of A. For I={1,3} the border is (a12,a32)=(1,2), orthogonal to (2,−1), which spans E−1 of (0223), with the same conclusion.

Exercise (B2: Bounds without a characteristic polynomial)

Let A∈M4(R) be symmetric with A=(7∗∗∗∗2∗∗∗∗−3∗∗∗∗1), the entries marked ∗ being unknown.

Solution.

  • (a) Each diagonal entry is a 1×1 principal submatrix, so λ1(A)≥max{7,2,−3,1}=7 and λ4(A)≤min{7,2,−3,1}=−3. In particular A is indefinite whatever the unknown entries are.

  • (b) The submatrix B=(7442) has trace 9 and determinant 14−16=−2, so its characteristic polynomial is x2−9x−2 and its eigenvalues are λ1(B)=9+892≈9.2170,λ2(B)=9−892≈−0.2170. Here n=4 and m=2, and the left-hand inequality λi(A)≥λi(B) of Corollary 17.4.3 at i=1,2 gives λ1(A) ≥ 9+892≈9.2170,λ2(A) ≥ 9−892≈−0.2170. The first improves the bound 7 from (a). The second is new: the diagonal entries alone say nothing about λ2(A). The right-hand inequalities λi(B)≥λi+2(A) give λ3(A)≤9.2170, which is true but of little use, and λ4(A)≤−0.2170, which is weaker than the bound λ4(A)≤−3 of (a). Each principal submatrix contributes what it can, and one keeps the best of the bounds.

Exercise (B3: Deciding strictness)

For each of the following bordered matrices A=(Bbb⊤c), decide without computing pA whether the interlacing between A and B is strict, and justify with Theorem 17.4.4.

  • (a) B=diag(3,1), b=(0,1), c=0.
  • (b) B=diag(3,1), b=(2,−1), c=5.
  • (c) B=(4004), b=(1,2), c=0.

Solution. In each case B is diagonal, so the standard basis vectors are an orthonormal eigenbasis and Pνb is just the part of b in the coordinates where the diagonal entry equals ν.

  • (a) Not strict. The eigenvalues 3,1 of B are distinct, but P3b=(0,0)=0. By Theorem 17.4.4 (b), 3 is an eigenvalue of A, so one of the inequalities λ1(A)≥3≥λ2(A) is an equality. (Indeed pA(x)=(x−3)(x2−x−1).)

  • (b) Strict. The eigenvalues 3,1 are distinct and P3b=(2,0)≠0, P1b=(0,−1)≠0. Both clauses of Theorem 17.4.4 (c) hold. (Indeed pA(x)=(x−2)(x2−7x+4), with roots 7+332≈6.372, 2 and 7−332≈0.628, and 6.372>3>2>1>0.628.)

  • (c) Not strict. P4b=b=(1,2)≠0, so the second clause holds, but B has the repeated eigenvalue 4 and the first clause fails. By Theorem 17.4.4 (a) with m4=2, the number 4 is an eigenvalue of A of multiplicity exactly 1, and the middle of the chain reads 4≥λ2(A)≥4. (Indeed pA(x)=(x−5)(x−4)(x+1).)

C. Going deeper

Exercise (C1: Rank under compression)

Let A∈Mn(F) be Hermitian of rank r, and let B=AI,I with |I|=n−m.

  • (a) Prove that rankB≥r−2m.
  • (b) Show that the bound is attained, by exhibiting for each n and each m≤n/2 a Hermitian A of rank 2m with a principal submatrix of size n−m and rank 0.

Solution.

  • (a) For a Hermitian matrix M, diagonalizing by the spectral theorem gives rankM=n+(M)+n−(M), the total number of non-zero eigenvalues. By Corollary 17.4.7, n+(B)≥n+(A)−m and n−(B)≥n−(A)−m. Adding, rankB=n+(B)+n−(B) ≥ (n+(A)+n−(A))−2m=r−2m.

  • (b) Let m≤n/2 and let I={m+1,m+2,…,n}, of size n−m. Define A by ai,m+i=am+i,i=1 for 1≤i≤m, and all other entries 0. Then A is real symmetric. It is a direct sum of m copies of (0110), acting on the coordinate pairs {i,m+i}, together with zeros on the remaining coordinates, so its eigenvalues are ±1 each m times and 0 otherwise; hence rankA=2m and n+(A)=n−(A)=m. Every non-zero entry of A has one index in {1,…,m}, which I omits, so AI,I=0 and rankAI,I=0=rankA−2m. The bound in (a) is therefore sharp for every m.

Exercise (C2: The converse, in the smallest case)

Theorem 17.4.2 says that the spectrum of a corner interlaces the spectrum of the whole. This exercise asks whether every interlacing pattern occurs.

  • (a) Let λ1>μ1>λ2 be real numbers. Find b∈R and c∈R such that A=(μ1bbc) has eigenvalues λ1 and λ2.
  • (b) Deduce that for n=2 the interlacing condition of Theorem 17.4.2 is not merely necessary but also sufficient.
  • (c) Explain why the strict hypothesis λ1>μ1>λ2 may be relaxed to λ1≥μ1≥λ2.

Hint: match the trace and the determinant.

Solution.

  • (a) A real symmetric 2×2 matrix with eigenvalues λ1,λ2 must have trA=λ1+λ2 and detA=λ1λ2. The first forces c=λ1+λ2−μ1, and the second then forces μ1c−b2=λ1λ2, that is, b2=μ1(λ1+λ2−μ1)−λ1λ2=−(μ1−λ1)(μ1−λ2)=(λ1−μ1)(μ1−λ2), where the middle equality is the expansion −(μ12−(λ1+λ2)μ1+λ1λ2). The hypothesis λ1>μ1>λ2 makes both factors positive, so b=(λ1−μ1)(μ1−λ2) is a real number. With this b and c, the matrix A is real symmetric with the prescribed trace and determinant, hence with characteristic polynomial x2−(λ1+λ2)x+λ1λ2=(x−λ1)(x−λ2), hence with eigenvalues λ1 and λ2.

    For example λ1=4, μ1=2, λ2=1 gives c=3 and b2=2, that is A=(2223), whose trace is 5 and determinant 6−2=4, as required.

  • (b) Necessity is Theorem 17.4.2 with n=2 and I={1}. Sufficiency is part (a): given any λ1>μ1>λ2, a Hermitian A with spectrum {λ1,λ2} and A{1},{1}=(μ1) exists.

  • (c) If μ1=λ1 or μ1=λ2 then the formula gives b=0, and A=diag(μ1,λ1+λ2−μ1) is diagonal with the two entries λ1 and λ2 in one order or the other. So the construction still works, with a zero border; consistently with Theorem 17.4.4, a zero border is exactly what makes the interlacing non-strict.

Exercise (C3: How far apart can the extremes be)

Let A∈Mn(F) be Hermitian.

  • (a) Prove that λ1(A)−λn(A)≥maxi,j |aii−ajj|.
  • (b) Prove that λ1(A)−λn(A)≥2maxi≠j|aij|.

Hint for (b): look at a 2×2 principal submatrix.

Solution.

  • (a) By Example (The diagonal lies inside the spectrum), λn(A)≤akk≤λ1(A) for every k. So for any i,j, both aii−ajj and ajj−aii are at most λ1(A)−λn(A), which gives the bound after taking the maximum over i,j.

  • (b) Fix i≠j and let I={i,j}, so that B=AI,I=(aiiaijaij―ajj) is Hermitian of size 2. Its trace is aii+ajj and its determinant is aiiajj−|aij|2, so its eigenvalues are aii+ajj2±(aii−ajj2)2+|aij|2, and their difference is 2((aii−ajj)/2)2+|aij|2≥2|aij|. Now apply Corollary 17.4.3 with m=n−2: it gives λ1(A)≥λ1(B) and λ2(B)≥λ2+(n−2)(A)=λn(A), so λ1(A)−λn(A) ≥ λ1(B)−λ2(B) ≥ 2|aij|. Taking the maximum over the pairs i≠j finishes the proof. (For n=2 the step with m=0 is the identity B=A, and the conclusion is the computation itself.)

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.