Skip to content

Eigenvalues and Eigenvectors

An eigenvector of a linear operator is a nonzero vector whose one-dimensional span is preserved by that operator. The corresponding eigenvalue is the scalar by which the vector is multiplied. Eigenspaces collect all vectors with the same eigenvalue and expose invariant directions or subspaces of the map.

Eigenvalues are properties of the operator, not of one matrix representation. Coordinates and eigenvector components change with basis; the eigenvalues and eigenspace dimensions do not.

Eigenproblems identify modes on which a linear transformation acts without mixing directions. They appear whenever one asks for:

  • stationary modes of a Hamiltonian or differential equation;
  • principal axes or normal modes;
  • invariant subspaces of a symmetry;
  • decay, growth, or oscillation rates in linear dynamics;
  • a basis that simplifies powers, exponentials, or functions of an operator;
  • possible sharp values of a quantum observable.

The final item requires additional physical structure. An arbitrary linear operator can have eigenvalues without representing an observable.

Let A:V→VA:V\to V be a linear operator on a vector space over a field F\mathbb F. A scalar λ∈F\lambda\in\mathbb F is an eigenvalue when there exists a nonzero vector v∈Vv\in V such that

Av=λv.Av=\lambda v.

The vector vv is an eigenvector belonging to λ\lambda. Equivalently,

(A−λI)v=0.(A-\lambda I)v=0.

Thus λ\lambda is an eigenvalue exactly when A−λIA-\lambda I has a nontrivial kernel.

The zero vector is excluded because

A0=λ0A0=\lambda0

holds for every scalar λ\lambda and therefore distinguishes nothing. By contrast, zero is a valid eigenvalue when a nonzero vector lies in ker⁡A\ker A.

For a fixed eigenvalue λ\lambda, the eigenspace is

Eλ=ker⁡(A−λI).\mathcal E_\lambda = \ker(A-\lambda I).

It is a vector subspace containing the zero vector and all eigenvectors with eigenvalue λ\lambda. If u,v∈Eλu,v\in\mathcal E_\lambda, then

A(au+bv)=λ(au+bv).A(au+bv) = \lambda(au+bv).

Every nonzero linear combination within one eigenspace remains an eigenvector with the same eigenvalue. A sum of eigenvectors from distinct eigenspaces is generally not an eigenvector.

Eigenvectors are directions, not normalized objects. If vv is an eigenvector, then every cvcv with c≠0c\ne0 is another eigenvector for the same eigenvalue. Normalization becomes convenient after an inner product is supplied, but it is not part of the eigenvalue definition.

In a dd-dimensional space, choose a basis and represent AA by a d×dd\times d matrix. A nonzero solution of

(A−λI)v=0(A-\lambda I)v=0

exists exactly when the matrix is singular. Therefore the eigenvalues are the roots of the characteristic polynomial

pA(λ)=det⁡(λI−A).p_A(\lambda) = \det(\lambda I-A).

This is a monic polynomial of degree dd. Over C\mathbb C, it factors as

pA(λ)=∏j=1r(λ−λj)mj,p_A(\lambda) = \prod_{j=1}^{r} (\lambda-\lambda_j)^{m_j},

where λ1,…,λr\lambda_1,\ldots,\lambda_r are the distinct eigenvalues and mjm_j are their algebraic multiplicities. Counting multiplicity,

∑j=1rmj=d.\sum_{j=1}^{r}m_j=d.

For small symbolic matrices, the determinant equation is useful. For numerical matrices, explicitly forming a characteristic polynomial is usually a poor algorithm because its coefficients can be sensitive and root finding can magnify error. Stable eigensolvers work with the matrix more directly.

If a change of basis replaces AA by the similar matrix

A′=S−1AS,A'=S^{-1}AS,

then

pA′(λ)=det⁡(λI−S−1AS)=det⁡(S−1(λI−A)S)=pA(λ).\begin{aligned} p_{A'}(\lambda) &= \det(\lambda I-S^{-1}AS) \\ &= \det\left( S^{-1}(\lambda I-A)S \right) \\ &= p_A(\lambda). \end{aligned}

Similar matrices therefore have the same eigenvalues with the same algebraic multiplicities. Their eigenvector coordinate columns are related by v′=S−1vv'=S^{-1}v.

The characteristic polynomial also connects eigenvalues to familiar invariants. Over C\mathbb C,

tr⁡A=∑j=1dλj,det⁡A=∏j=1dλj,\operatorname{tr}A = \sum_{j=1}^{d}\lambda_j, \qquad \det A = \prod_{j=1}^{d}\lambda_j,

with eigenvalues repeated according to algebraic multiplicity.

For an eigenvalue λ\lambda:

  • its algebraic multiplicity mλm_\lambda is its multiplicity as a root of pAp_A;
  • its geometric multiplicity is
gλ=dim⁡ker⁡(A−λI).g_\lambda = \dim\ker(A-\lambda I).

They satisfy

1≤gλ≤mλ.1\leq g_\lambda\leq m_\lambda.

An eigenvalue with mλ>1m_\lambda>1 is a repeated root. In quantum mechanics, “degenerate eigenvalue” usually means that its eigenspace has dimension gλ>1g_\lambda>1. For Hermitian matrices, algebraic and geometric multiplicities agree, so this distinction causes no conflict. For a general matrix, a repeated eigenvalue can have too few eigenvectors.

The total degeneracy relevant to a sharp quantum outcome is the dimension of its eigenspace. A basis inside that eigenspace is not unique; the subspace is the invariant object.

Distinct Eigenvalues Give Independent Eigenvectors

Section titled “Distinct Eigenvalues Give Independent Eigenvectors”

Eigenvectors belonging to distinct eigenvalues are linearly independent. A short induction proof shows why.

Suppose v1,…,vnv_1,\ldots,v_n have distinct eigenvalues λ1,…,λn\lambda_1,\ldots,\lambda_n and

∑j=1ncjvj=0.\sum_{j=1}^{n}c_jv_j=0.

Apply A−λnIA-\lambda_n I:

∑j=1n−1cj(λj−λn)vj=0.\sum_{j=1}^{n-1} c_j(\lambda_j-\lambda_n)v_j =0.

By the induction hypothesis, v1,…,vn−1v_1,\ldots,v_{n-1} are independent. Because λj−λn≠0\lambda_j-\lambda_n\ne0, one obtains c1=⋯=cn−1=0c_1=\cdots=c_{n-1}=0, and the original relation then gives cn=0c_n=0.

This result guarantees enough eigenvectors when a d×dd\times d matrix has dd distinct eigenvalues. Repeated eigenvalues require inspection of their eigenspaces.

An operator is diagonalizable when its eigenspaces together span the whole space:

V=⨁λEλ.V = \bigoplus_{\lambda} \mathcal E_\lambda.

Equivalently,

∑λgλ=d.\sum_{\lambda}g_\lambda=d.

A characteristic polynomial that splits over C\mathbb C guarantees eigenvalues, but not enough independent eigenvectors. A defective matrix has at least one gλ<mλg_\lambda<m_\lambda and cannot be diagonalized.

The complete criterion, Jordan-block counterexamples, and change-of-basis construction belong to Diagonalization. Normal operators form the especially important class whose eigenvectors can be chosen orthonormal.

Worked Example: A Coupled Two-Level Matrix

Section titled “Worked Example: A Coupled Two-Level Matrix”

Consider

A=(2112).A = \begin{pmatrix} 2&1\\ 1&2 \end{pmatrix}.

The characteristic polynomial is

pA(λ)=det⁡(λ−2−1−1λ−2)=(λ−2)2−1=(λ−1)(λ−3).\begin{aligned} p_A(\lambda) &= \det \begin{pmatrix} \lambda-2&-1\\ -1&\lambda-2 \end{pmatrix} \\ &= (\lambda-2)^2-1 \\ &= (\lambda-1)(\lambda-3). \end{aligned}

For λ=3\lambda=3,

(A−3I)v=0(A-3I)v=0

gives v1=v2v_1=v_2. For λ=1\lambda=1, it gives v1=−v2v_1=-v_2. Normalized eigenvectors are therefore

v+=12(11),v−=12(1−1),v_+ = \frac{1}{\sqrt2} \begin{pmatrix} 1\\ 1 \end{pmatrix}, \qquad v_- = \frac{1}{\sqrt2} \begin{pmatrix} 1\\ -1 \end{pmatrix},

with

Av+=3v+,Av−=v−.Av_+=3v_+, \qquad Av_-=v_-.

The matrix is Hermitian, so the eigenvalues are real and the two eigenspaces are orthogonal. Those conclusions are not consequences of the determinant calculation alone; they follow from Hermiticity.

Let

B=(200020005).B = \begin{pmatrix} 2&0&0\\ 0&2&0\\ 0&0&5 \end{pmatrix}.

Its characteristic polynomial is

pB(λ)=(λ−2)2(λ−5).p_B(\lambda) = (\lambda-2)^2(\lambda-5).

The eigenspaces are

E2=span⁡{(100),(010)}\mathcal E_2 = \operatorname{span} \left\lbrace \begin{pmatrix}1\\0\\0\end{pmatrix}, \begin{pmatrix}0\\1\\0\end{pmatrix} \right\rbrace

and

E5=span⁡{(001)}.\mathcal E_5 = \operatorname{span} \left\lbrace \begin{pmatrix}0\\0\\1\end{pmatrix} \right\rbrace.

The eigenvalue 22 has algebraic and geometric multiplicity two. Every nonzero vector in the plane E2\mathcal E_2 is an eigenvector with eigenvalue 22. Rotating the basis inside that plane changes the eigenvectors used in a coordinate description but does not change the eigenspace.

Whether eigenvalues exist can depend on the field. The real rotation matrix

R=(0−110)R = \begin{pmatrix} 0&-1\\ 1&0 \end{pmatrix}

has characteristic polynomial

pR(λ)=λ2+1.p_R(\lambda)=\lambda^2+1.

It has no real eigenvalues. On the complexified space, it has eigenvalues +i+i and −i-i, with eigenvectors proportional to

(1−i),(1i),\begin{pmatrix}1\\-i\end{pmatrix}, \qquad \begin{pmatrix}1\\i\end{pmatrix},

respectively. Quantum state spaces are complex, so finite-dimensional characteristic polynomials always split into linear factors, although a matrix may still be defective.

The eigenvalue equation alone does not guarantee real eigenvalues, orthogonality, or a complete eigenbasis.

Operator classEigenvalue facts
General complex matrixeigenvalues may be complex; eigenvectors may be incomplete or nonorthogonal
Diagonalizable matrixhas a basis of eigenvectors, not necessarily orthogonal
Normal matrixhas an orthonormal eigenbasis
Hermitian matrixhas an orthonormal eigenbasis and real eigenvalues
Unitary matrixhas an orthonormal eigenbasis and eigenvalues of unit modulus

Hermitian Operators and Normal Operators own the proofs and spectral consequences for those classes.

For a non-Hermitian matrix, one may need both right eigenvectors,

Avj=λjvj,Av_j=\lambda_jv_j,

and left eigenvectors, represented by functionals satisfying

wj†A=λjwj†.w_j^\dagger A = \lambda_jw_j^\dagger.

They need not be related by ordinary conjugate transpose, and their sensitivity can be much worse than in the Hermitian case. Do not silently import orthogonality or projector formulas from Hermitian spectral theory.

In the finite-dimensional projective measurement model, a Hermitian observable AA has real eigenvalues that label possible sharp outcomes. A nondegenerate eigenvector determines a one-dimensional outcome ray; a degenerate eigenvalue corresponds to its entire eigenspace and orthogonal projector.

This interpretation has conditions:

  • the operator must represent the observable under the quantum postulates;
  • eigenvalues carry the physical units of that operator;
  • states must be normalized before probabilities are assigned;
  • probabilities come from spectral projectors and the Born rule;
  • an arbitrary operator’s eigenvalues are not automatically measurable quantities.

The physical treatment is Eigenvalues and Eigenstates. This page supplies the finite-dimensional linear algebra that treatment uses.

Differential and Infinite-Dimensional Eigenproblems

Section titled “Differential and Infinite-Dimensional Eigenproblems”

For differential operators, the equation

Aψ=λψA\psi=\lambda\psi

still defines an eigenproblem, but admissible functions, boundary conditions, and the operator domain are part of the problem. A formal solution of the differential equation need not be an allowed eigenvector.

Infinite-dimensional operators can also have continuous spectrum, where a spectral value has no normalizable eigenvector. Characteristic determinants are not the general tool. See Eigenvalue Problems for boundary-value formulations and Discrete and Continuous Spectra for the physical distinction.

For a computed pair (λ~,v~)(\widetilde\lambda,\widetilde v), check the residual

r=Av~−λ~v~.r = A\widetilde v -\widetilde\lambda\widetilde v.

A small relative residual confirms that the pair nearly satisfies the matrix equation. It does not by itself guarantee that an eigenvector is accurately determined when eigenvalues are clustered or the matrix is highly nonnormal.

Within a degenerate eigenspace, individual numerical eigenvectors are not unique: a solver may return any orthonormal basis of that subspace. Compare the projector or invariant subspace rather than matching vectors component by component. Numerical algorithms, conditioning, residual scaling, and truncation checks belong to Matrix Diagonalization.

  • Allowing the zero vector as an eigenvector. It satisfies every formal eigenvalue equation and therefore carries no eigenvalue information.
  • Treating zero as an invalid eigenvalue. It is valid when the kernel is nontrivial.
  • Finding roots but not eigenspaces. The characteristic polynomial alone does not provide eigenvectors or geometric multiplicities.
  • Confusing repeated roots with complete degeneracy data. Distinguish algebraic from geometric multiplicity.
  • Assuming every matrix is diagonalizable. Defective matrices lack enough independent eigenvectors.
  • Assuming eigenvectors are automatically orthogonal. This is guaranteed for normal operators, not general matrices.
  • Treating eigenvector coordinates as unique. Scaling, phase, basis choice, and rotations within a degenerate eigenspace all change them.
  • Assigning measurement meaning to every eigenvalue. The observable postulate and spectral projectors supply the physical interpretation.
  • Using characteristic polynomials as numerical eigensolvers. Use stable matrix algorithms and residual checks.
  1. Find normalized eigenvectors as well as eigenvalues of
σx=(0110).\sigma_x = \begin{pmatrix} 0&1\\ 1&0 \end{pmatrix}.
Solution

The characteristic equation is

det⁡(λI−σx)=λ2−1=0.\det(\lambda I-\sigma_x) = \lambda^2-1 =0.

Thus λ±=±1\lambda_\pm=\pm1. For λ=+1\lambda=+1, the equation gives v1=v2v_1=v_2; for λ=−1\lambda=-1, it gives v1=−v2v_1=-v_2. Normalized choices are

v+=12(11),v−=12(1−1).v_+ = \frac{1}{\sqrt2} \begin{pmatrix}1\\1\end{pmatrix}, \qquad v_- = \frac{1}{\sqrt2} \begin{pmatrix}1\\-1\end{pmatrix}.

They are orthogonal, as expected because σx\sigma_x is Hermitian.

  1. For
A=(410040002),A = \begin{pmatrix} 4&1&0\\ 0&4&0\\ 0&0&2 \end{pmatrix},

find the algebraic and geometric multiplicities and decide whether AA is diagonalizable.

Solution

The matrix is triangular, so

pA(λ)=(λ−4)2(λ−2).p_A(\lambda) = (\lambda-4)^2(\lambda-2).

Thus m4=2m_4=2 and m2=1m_2=1. For λ=4\lambda=4,

A−4I=(01000000−2),A-4I = \begin{pmatrix} 0&1&0\\ 0&0&0\\ 0&0&-2 \end{pmatrix},

so the kernel is spanned by (1,0,0)T(1,0,0)^{\mathsf T} and g4=1g_4=1. For λ=2\lambda=2, the eigenspace is spanned by (0,0,1)T(0,0,1)^{\mathsf T}, so g2=1g_2=1.

Because

g4+g2=2<3,g_4+g_2=2<3,

the matrix does not have an eigenbasis and is not diagonalizable.

  1. Let v1v_1 and v2v_2 be eigenvectors with distinct eigenvalues λ1≠λ2\lambda_1\ne\lambda_2. Prove directly that they are linearly independent.
Solution

Assume

c1v1+c2v2=0.c_1v_1+c_2v_2=0.

Apply A−λ2IA-\lambda_2I:

c1(λ1−λ2)v1=0.c_1(\lambda_1-\lambda_2)v_1=0.

Because v1≠0v_1\ne0 and λ1−λ2≠0\lambda_1-\lambda_2\ne0, one has c1=0c_1=0. The original relation then gives c2=0c_2=0. Hence the vectors are linearly independent.

  1. A numerical diagonalization of a Hermitian matrix with a two-fold degenerate eigenvalue returns different normalized eigenvectors on two machines, but both pairs span the same two-dimensional subspace. Is this a disagreement?
Solution

No. An eigenspace of dimension two has infinitely many orthonormal bases. Floating-point details may select different bases within that same invariant subspace. The basis-independent comparison is the orthogonal projector

P=∣v1⟩⟨v1∣+∣v2⟩⟨v2∣,P = \lvert v_1\rangle\langle v_1\rvert +\lvert v_2\rangle\langle v_2\rvert,

or an equivalent subspace-distance diagnostic. If the projectors agree within numerical tolerance and the residuals are small, the eigenspace results agree.

  • S. Axler, Linear Algebra Done Right, 3rd ed., Springer, 2015.
  • R. A. Horn and C. R. Johnson, Matrix Analysis, 2nd ed., Cambridge University Press, 2013.
  • B. C. Hall, Quantum Theory for Mathematicians, Springer, 2013.
  • L. N. Trefethen and D. Bau III, Numerical Linear Algebra, SIAM, 1997.