Mehmet Emre Tiryaki
Back to Notes

Majorization

Majorization

Let x,y∈Rnx, y \in \mathbb{R}^n and order their entries increasingly, xℓ1≤xℓ2≤⋯≤xℓnx_{\ell_1} \le x_{\ell_2} \le \cdots \le x_{\ell_n} and ym1≤ym2≤⋯≤ymny_{m_1} \le y_{m_2} \le \cdots \le y_{m_n}. xx majorizes yy when for all k=1,…,nk = 1, \ldots, n

∑i=1kxℓi≥∑i=1kymi\sum_{i=1}^{k} x_{\ell_i} \ge \sum_{i=1}^{k} y_{m_i}

with equality at k=nk = n.

Note this is stated with increasing order. Many references order decreasingly instead, which reverses the direction of every statement below (e.g. they say λ(A)\lambda(A) majorizes diag(A)\mathrm{diag}(A)).

Diagonal and Eigenvalues

For A∈MnA \in M_n Hermitian, diag(A)=(a11,…,ann)\mathrm{diag}(A) = (a_{11}, \ldots, a_{nn}) majorizes λ(A)=(λ1(A),…,λn(A))\lambda(A) = (\lambda_1(A), \ldots, \lambda_n(A)).

Proof By induction, n=1n = 1 is trivial. Assume it holds for n−1n - 1. Order the diagonal aℓ1ℓ1≤⋯≤aℓnℓna_{\ell_1 \ell_1} \le \cdots \le a_{\ell_n \ell_n} and let BB be the submatrix from deleting row and column ℓn\ell_n. Then diag(B)\mathrm{diag}(B) is the n−1n - 1 smallest diagonal entries of AA, and for all k=1,…,n−1k = 1, \ldots, n - 1, by interlacing and then induction,

∑i=1kλi(A)≤∑i=1kλi(B)≤∑i=1kaℓiℓi\sum_{i=1}^{k} \lambda_i(A) \le \sum_{i=1}^{k} \lambda_i(B) \le \sum_{i=1}^{k} a_{\ell_i \ell_i}

For k=nk = n, ∑i=1nλi(A)=tr(A)=∑i=1naℓiℓi\sum_{i=1}^{n} \lambda_i(A) = \mathrm{tr}(A) = \sum_{i=1}^{n} a_{\ell_i \ell_i}.

Eigenvalue Sums as Traces

For A∈MnA \in M_n Hermitian and all r=1,…,nr = 1, \ldots, n,

∑k=1rλk(A)=min⁡U∈Mn,rorthonormal colstr(U∗AU),∑k=0r−1λn−k(A)=max⁡U∈Mn,rorthonormal colstr(U∗AU)\sum_{k=1}^{r} \lambda_k(A) = \min_{\substack{U \in M_{n,r} \\ \text{orthonormal cols}}} \mathrm{tr}(U^*AU), \qquad \sum_{k=0}^{r-1} \lambda_{n-k}(A) = \max_{\substack{U \in M_{n,r} \\ \text{orthonormal cols}}} \mathrm{tr}(U^*AU)

Proof

  1. For any U∈Mn,rU \in M_{n,r} with orthonormal columns, there is a unitary V=[ U∣∗ ]∈MnV = [\, U \mid * \,] \in M_n by Gram-Schmidt.
  2. U∗AUU^*AU is the leading r×rr \times r principal submatrix of V∗AVV^*AV.
  3. So λk(A)=λk(V∗AV)≤λk(U∗AU)\lambda_k(A) = \lambda_k(V^*AV) \le \lambda_k(U^*AU) by interlacing.
  4. Summing over kk, ∑k=1rλk(A)≤∑k=1rλk(U∗AU)=tr(U∗AU)\sum_{k=1}^{r} \lambda_k(A) \le \sum_{k=1}^{r} \lambda_k(U^*AU) = \mathrm{tr}(U^*AU), so we just need a UU attaining the min.
  5. Let W=[w1∣⋯∣wr]W = [w_1 | \cdots | w_r] with wiw_i orthonormal eigenvectors of λ1,…,λr\lambda_1, \ldots, \lambda_r. Then W∗AW=diag(λ1,…,λr)W^*AW = \mathrm{diag}(\lambda_1, \ldots, \lambda_r) and tr(W∗AW)=∑k=1rλk(A)\mathrm{tr}(W^*AW) = \sum_{k=1}^{r} \lambda_k(A).

Eigenvalues of a Sum

For A,BA, B Hermitian, λ(A+B)\lambda(A + B) majorizes λ(A)+λ(B)\lambda(A) + \lambda(B).

Proof For all k=1,…,nk = 1, \ldots, n,

∑i=1kλi(A+B)=min⁡U∈Mn,korthonormal colstr(U∗(A+B)U)=min⁡(tr(U∗AU)+tr(U∗BU))≥min⁡tr(U∗AU)+min⁡tr(U∗BU)=∑i=1kλi(A)+∑i=1kλi(B)\begin{aligned} \sum_{i=1}^{k} \lambda_i(A + B) &= \min_{\substack{U \in M_{n,k} \\ \text{orthonormal cols}}} \mathrm{tr}(U^*(A + B)U) = \min \big(\mathrm{tr}(U^*AU) + \mathrm{tr}(U^*BU)\big) \\ &\ge \min \mathrm{tr}(U^*AU) + \min \mathrm{tr}(U^*BU) = \sum_{i=1}^{k} \lambda_i(A) + \sum_{i=1}^{k} \lambda_i(B) \end{aligned}

with equality at k=nk = n since both sides are tr(A+B)\mathrm{tr}(A + B).

Hadamard's Inequality

Let A∈MnA \in M_n be PSD. Then

det⁡A≤∏i=1naii\det A \le \prod_{i=1}^{n} a_{ii}

Proof For all ii, 0≤λ1(A)≤aii0 \le \lambda_1(A) \le a_{ii} by PSD and interlacing, so ∏iaii≥0\prod_i a_{ii} \ge 0. If AA is singular, det⁡A=0\det A = 0 and we are done.

Else AA is PD, so 0<λ1(A)≤aii0 < \lambda_1(A) \le a_{ii}. Let D=diag(1/a11,…,1/ann)D = \mathrm{diag}(1/\sqrt{a_{11}}, \ldots, 1/\sqrt{a_{nn}}). DADDAD is Hermitian PD with ones on the diagonal, so by AM-GM

det⁡A∏i=1naii=det⁡D2det⁡A=det⁡(DAD)=∏i=1nλi(DAD)≤(1n∑i=1nλi(DAD))n=(1ntr(DAD))n=(nn)n=1\frac{\det A}{\prod_{i=1}^{n} a_{ii}} = \det D^2 \det A = \det(DAD) = \prod_{i=1}^{n} \lambda_i(DAD) \le \Big(\frac{1}{n} \sum_{i=1}^{n} \lambda_i(DAD)\Big)^n = \Big(\frac{1}{n} \mathrm{tr}(DAD)\Big)^n = \Big(\frac{n}{n}\Big)^n = 1