|q⟩ Bad Qubits

intermediate · Physics · Operators, Eigenvalues & Spectral Theorem

Eigenvalue Equations

An eigenvalue equation asks a sharp question of an operator: is there a direction the operator leaves unchanged except for stretching? Those special directions — eigenvectors — and their stretch factors — eigenvalues — are where the physics lives. The possible outcomes of a measurement are the eigenvalues of the observable, and the states that give a definite outcome are its eigenvectors.

The equation

A nonzero vector v|v\rangle is an eigenvector of the operator AA with eigenvalue λ\lambda if

Av=λv.A|v\rangle = \lambda |v\rangle.

Applying AA does not rotate v|v\rangle into a new direction; it merely scales it by the number λ\lambda (possibly complex). We exclude v=0|v\rangle = 0, which would satisfy the equation trivially for any λ\lambda and tell us nothing. Note that if v|v\rangle is an eigenvector, so is any nonzero multiple cvc|v\rangle — eigenvectors come as whole directions (one-dimensional subspaces, or larger if the eigenvalue is degenerate).

Finding eigenvalues: the characteristic equation

Rewrite the equation as (Aλ1)v=0(A - \lambda\mathbb{1})|v\rangle = 0. A nonzero solution v|v\rangle exists precisely when the operator Aλ1A - \lambda\mathbb{1} is not invertible, which happens when its determinant vanishes:

det(Aλ1)=0.\det(A - \lambda\mathbb{1}) = 0.

This is the characteristic equation, a polynomial of degree nn in λ\lambda for an n×nn\times n matrix. Its roots are the eigenvalues.

For a 2×22\times 2 matrix A=(abcd)A = \begin{pmatrix} a & b \\ c & d\end{pmatrix}, expanding the determinant gives the compact form

λ2(a+dtrA)λ+(adbcdetA)=0,\lambda^2 - (\underbrace{a+d}_{\operatorname{tr} A})\,\lambda + (\underbrace{ad - bc}_{\det A}) = 0,

so the two eigenvalues are

λ±=trA±(trA)24detA2.\lambda_\pm = \frac{\operatorname{tr}A \pm \sqrt{(\operatorname{tr}A)^2 - 4\det A}}{2}.

A useful built-in check: the eigenvalues always satisfy λ++λ=trA\lambda_+ + \lambda_- = \operatorname{tr}A and λ+λ=detA\lambda_+ \lambda_- = \det A.

Worked example

Take the real symmetric matrix

A=(3113).A = \begin{pmatrix} 3 & 1 \\ 1 & 3 \end{pmatrix}.

Here trA=6\operatorname{tr}A = 6 and detA=91=8\det A = 9 - 1 = 8, so the characteristic equation is λ26λ+8=0=(λ4)(λ2)\lambda^2 - 6\lambda + 8 = 0 = (\lambda - 4)(\lambda - 2). The eigenvalues are λ+=4\lambda_+ = 4 and λ=2\lambda_- = 2, both real — as guaranteed for a Hermitian (here symmetric) operator. Their sum is 6=trA6 = \operatorname{tr}A and product is 8=detA8 = \det A, confirming the answer.

Finding the eigenvectors

Once an eigenvalue is known, solve the linear system (Aλ1)v=0(A - \lambda\mathbb{1})|v\rangle = 0. For λ+=4\lambda_+ = 4,

A41=(1111),A - 4\mathbb{1} = \begin{pmatrix} -1 & 1 \\ 1 & -1 \end{pmatrix},

and the equation v1+v2=0-v_1 + v_2 = 0 forces v1=v2v_1 = v_2. A normalized eigenvector is v+=12(11)|v_+\rangle = \tfrac{1}{\sqrt2}\begin{pmatrix} 1 \\ 1 \end{pmatrix}. For λ=2\lambda_- = 2 the same method gives v=12(11)|v_-\rangle = \tfrac{1}{\sqrt2}\begin{pmatrix} 1 \\ -1 \end{pmatrix}. The two eigenvectors are orthogonal, as they must be for distinct eigenvalues of a Hermitian operator.

Try it

Solve Av=λvA|v\rangle = \lambda|v\rangle for A=(3113)A = \begin{pmatrix} 3 & 1 \\ 1 & 3\end{pmatrix} using the trace–determinant formula, and return the larger of the two eigenvalues.

Run your code to see the quantum state.

Sign in on the full site to ask questions and join the discussion.