← Exercise library← All exercises for this theme
Theme 2 — Hilbert Spaces and Dirac Notation

Diagonalisation of a two-by-two Hermitian matrix

Keywords: diagonalisation · Hermitian matrix · two-level system · mixing angle · avoided crossing · level repulsion

Exercise 1 : Diagonalisation of a two-by-two Hermitian matrix

Consider the most general 2×22\times2 Hermitian matrix,

H=(abb∗d),a,d∈R,b=∣b∣eiφ∈C.\begin{aligned} H=\begin{pmatrix}a&b\\ b^*&d\end{pmatrix},\qquad a,d\in\R,\quad b=|b|e^{i\varphi}\in\C . \end{aligned}

For example, it represents the Hamiltonian of a two-level system in a basis in which aa and dd are the energies of the two states and bb is their coupling. Let

m=a+d2,Δ=(a−d2)2+∣b∣2.m=\frac{a+d}{2},\qquad\Delta=\sqrt{\Bigl(\frac{a-d}{2}\Bigr)^2+|b|^2}.
Question 1
Determine the eigenvalues of HH.

Solution
The characteristic polynomial is

det⁡(H−λ1)=(a−λ)(d−λ)−∣b∣2=λ2−(a+d)λ+ad−∣b∣2.\det(H-\lambda\mathbf{1})=(a-\lambda)(d-\lambda)-|b|^2=\lambda^2-(a+d)\lambda+ad-|b|^2 .

Its discriminant is (a+d)2−4ad+4∣b∣2=(a−d)2+4∣b∣2=4Δ2≥0(a+d)^2-4ad+4|b|^2=(a-d)^2+4|b|^2=4\Delta^2\geq0, and the eigenvalues are

λ±=a+d2±12(a−d)2+4∣b∣2=m±Δ.\lambda_\pm=\frac{a+d}{2}\pm\frac12\sqrt{(a-d)^2+4|b|^2}=m\pm\Delta .

They are real, as they must be for a Hermitian matrix. They are equal if and only if Δ=0\Delta=0, that is, if a=da=d and b=0b=0: the matrix is then proportional to the identity.

Question 2
Suppose that Δ>0\Delta>0, and define the angle θ∈[0,π]\theta\in[0,\pi] by cos⁡θ=a−d2Δ\cos\theta=\frac{a-d}{2\Delta} and sin⁡θ=∣b∣Δ\sin\theta=\frac{|b|}{\Delta}. Show that H=m 1+Δ KH=m\,\mathbf{1}+\Delta\,K with

K=(cos⁡θeiφsin⁡θe−iφsin⁡θ−cos⁡θ),\begin{aligned} K=\begin{pmatrix}\cos\theta&e^{i\varphi}\sin\theta\\ e^{-i\varphi}\sin\theta&-\cos\theta\end{pmatrix}, \end{aligned}

then show that the vectors

∣+⟩=(cos⁡θ2e−iφsin⁡θ2),∣−⟩=(−eiφsin⁡θ2cos⁡θ2)\ket{+}=\begin{pmatrix}\cos\frac\theta2\\ e^{-i\varphi}\sin\frac\theta2\end{pmatrix},\qquad \ket{-}=\begin{pmatrix}-e^{i\varphi}\sin\frac\theta2\\ \cos\frac\theta2\end{pmatrix}

form an orthonormal basis of eigenvectors of HH.

Solution
Direct calculation gives m+Δcos⁡θ=a+d2+a−d2=am+\Delta\cos\theta=\frac{a+d}{2}+\frac{a-d}{2}=a, m−Δcos⁡θ=dm-\Delta\cos\theta=d, and Δeiφsin⁡θ=∣b∣eiφ=b\Delta e^{i\varphi}\sin\theta=|b|e^{i\varphi}=b. The eigenvectors of HH are those of KK, with eigenvalues m+Δkm+\Delta k if kk is an eigenvalue of KK. Let us calculate K∣+⟩K\ket+, using the identities cos⁡θ=cos⁡2θ2−sin⁡2θ2\cos\theta=\cos^2\frac\theta2-\sin^2\frac\theta2 and sin⁡θ=2sin⁡θ2cos⁡θ2\sin\theta=2\sin\frac\theta2\cos\frac\theta2. The first component is

cos⁡θcos⁡θ2+eiφsin⁡θ e−iφsin⁡θ2=cos⁡θcos⁡θ2+sin⁡θsin⁡θ2=cos⁡(θ−θ2)=cos⁡θ2,\cos\theta\cos\tfrac\theta2+e^{i\varphi}\sin\theta\,e^{-i\varphi}\sin\tfrac\theta2=\cos\theta\cos\tfrac\theta2+\sin\theta\sin\tfrac\theta2=\cos\bigl(\theta-\tfrac\theta2\bigr)=\cos\tfrac\theta2,

and the second is

e−iφsin⁡θcos⁡θ2−cos⁡θ e−iφsin⁡θ2=e−iφsin⁡(θ−θ2)=e−iφsin⁡θ2.e^{-i\varphi}\sin\theta\cos\tfrac\theta2-\cos\theta\,e^{-i\varphi}\sin\tfrac\theta2=e^{-i\varphi}\sin\bigl(\theta-\tfrac\theta2\bigr)=e^{-i\varphi}\sin\tfrac\theta2 .

Thus K∣+⟩=∣+⟩K\ket+=\ket+, and one similarly shows that K∣−⟩=−∣−⟩K\ket-=-\ket-. The two vectors are normalised, and ⟨+|−⟩=−cos⁡θ2eiφsin⁡θ2+eiφsin⁡θ2cos⁡θ2=0\braket+-=-\cos\frac\theta2e^{i\varphi}\sin\frac\theta2+e^{i\varphi}\sin\frac\theta2\cos\frac\theta2=0. They are therefore eigenvectors of HH, with eigenvalues m±Δm\pm\Delta.

Question 3
Apply these results to H=(1i−i−1)H=\begin{pmatrix}1&i\\-i&-1\end{pmatrix}. Verify the first eigenvector by direct calculation.

Solution
Here a=1a=1, d=−1d=-1, b=ib=i, so m=0m=0, ∣b∣=1|b|=1, φ=π/2\varphi=\pi/2, and Δ=1+1=2\Delta=\sqrt{1+1}=\sqrt2. The eigenvalues are ±2\pm\sqrt2. The angle satisfies cos⁡θ=sin⁡θ=12\cos\theta=\sin\theta=\frac{1}{\sqrt2}, hence θ=π/4\theta=\pi/4, and e−iφ=−ie^{-i\varphi}=-i:

∣+⟩=(cos⁡π8−isin⁡π8),∣−⟩=(−isin⁡π8cos⁡π8).\ket+=\begin{pmatrix}\cos\frac\pi8\\-i\sin\frac\pi8\end{pmatrix},\qquad \ket-=\begin{pmatrix}-i\sin\frac\pi8\\ \cos\frac\pi8\end{pmatrix}.

Let us check the first component of H∣+⟩H\ket+: 1⋅cos⁡π8+i⋅(−i)sin⁡π8=cos⁡π8+sin⁡π81\cdot\cos\frac\pi8+i\cdot(-i)\sin\frac\pi8=\cos\frac\pi8+\sin\frac\pi8. It must equal 2cos⁡π8\sqrt2\cos\frac\pi8, which it does because cos⁡π8+sin⁡π8=2cos⁡(π8−π4)=2cos⁡π8\cos\frac\pi8+\sin\frac\pi8=\sqrt2\cos\bigl(\frac\pi8-\frac\pi4\bigr)=\sqrt2\cos\frac\pi8. Numerically, 0,924+0,383=1,307=2×0,9240{,}924+0{,}383=1{,}307=\sqrt2\times0{,}924.

Question 4
Consider the case a≠da\neq d with weak coupling, ∣b∣≪∣a−d∣|b|\ll|a-d|. Give an expansion of the eigenvalues to the first non-zero order in ∣b∣|b|, and describe the eigenvectors. What happens instead when a=da=d? Sketch qualitatively the two eigenvalues as functions of a−da-d, with mm and bb fixed.

Solution
Suppose a>da>d for definiteness. Expanding the square root gives

Δ=a−d21+4∣b∣2(a−d)2≃a−d2+∣b∣2a−d,λ+≃a+∣b∣2a−d,λ−≃d−∣b∣2a−d.\Delta=\frac{a-d}{2}\sqrt{1+\frac{4|b|^2}{(a-d)^2}}\simeq\frac{a-d}{2}+\frac{|b|^2}{a-d}, \qquad \lambda_+\simeq a+\frac{|b|^2}{a-d},\quad\lambda_-\simeq d-\frac{|b|^2}{a-d}.

The coupling pushes the two levels apart: the upper level rises and the lower level falls, by an amount proportional to the square of the coupling and inversely proportional to the initial separation. The angle θ\theta is small, θ≃2∣b∣/(a−d)\theta\simeq2|b|/(a-d), and the eigenvectors are close to the basis vectors: each is only slightly mixed with the other. When a=da=d, by contrast, θ=π/2\theta=\pi/2, and the eigenvectors 12(1,±e−iφ)\frac{1}{\sqrt2}(1,\pm e^{-i\varphi}) are equal-weight superpositions of the two basis states: the mixing is maximal, even for weak coupling, and the separation between the levels is 2∣b∣2|b|. As a−da-d varies, the two eigenvalues m±(a−d)2/4+∣b∣2m\pm\sqrt{(a-d)^2/4+|b|^2} trace the two branches of a hyperbola: they never cross, unlike the uncoupled energies aa and dd, which cross at a=da=d. This is the phenomenon of an avoided crossing, or level repulsion, which occurs throughout quantum physics, from molecules to qubits.