Skip to content

Eigenvalue Problems

An eigenvalue problem asks for nonzero states or functions that are rescaled by an operator:

Au=λu.A u=\lambda u.

The nonzero solutions uu are eigenvectors or eigenfunctions, and the numbers λ\lambda are eigenvalues. In quantum mechanics, the most important example is the stationary Schrödinger equation

Hψ=Eψ,H\psi=E\psi,

where HH is the Hamiltonian, EE is an allowed energy, and ψ\psi is an energy eigenstate when it is an acceptable state.

This page owns the mathematical problem type. The finite-dimensional algebraic definition is treated in Eigenvalues and Eigenvectors, while the full self-adjoint-operator structure is treated in Spectral Theorem, Practical Version.

The symbolic equation Au=λuA u=\lambda u is incomplete until the following are specified:

  • the space where uu lives;
  • the action of AA;
  • the domain on which AA is allowed to act;
  • any boundary, regularity, decay, or matching conditions.

For a matrix acting on Cn\mathbb C^n, this data is usually implicit. For a differential operator, it is not. The expression

−d2dx2-\frac{d^2}{dx^2}

has different spectra on the real line, on a finite interval with Dirichlet boundary conditions, on a ring with periodic boundary conditions, or on a half-line with a boundary condition at the endpoint.

This is why Boundary Conditions and Domains of Operators are not optional refinements. They are part of the eigenvalue problem.

For an n×nn\times n matrix AA, the equation

Av=λvAv=\lambda v

has a nonzero solution only when

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

This characteristic equation gives the eigenvalues. For each eigenvalue, the eigenspace is the null space of A−λIA-\lambda I:

ker⁡(A−λI)={v:(A−λI)v=0}.\ker(A-\lambda I) = \{v:(A-\lambda I)v=0\}.

If AA is Hermitian, its eigenvalues are real and eigenspaces belonging to distinct eigenvalues are orthogonal. When enough eigenvectors exist to form a basis, the operator can be represented diagonally. The finite-dimensional projector form is Spectral Decomposition.

A differential eigenvalue problem combines a differential equation with admissibility conditions. A simple model is

−u′′(x)=λu(x),0<x<L,-u''(x)=\lambda u(x), \qquad 0<x<L,

with

u(0)=u(L)=0.u(0)=u(L)=0.

The general solution for λ=k2>0\lambda=k^2>0 is

u(x)=Asin⁡(kx)+Bcos⁡(kx).u(x)=A\sin(kx)+B\cos(kx).

The first boundary condition gives B=0B=0. The second gives

sin⁡(kL)=0,\sin(kL)=0,

so

k=nπL,λn=n2π2L2,n=1,2,3,….k=\frac{n\pi}{L}, \qquad \lambda_n=\frac{n^2\pi^2}{L^2}, \qquad n=1,2,3,\ldots.

The boundary conditions have selected a discrete set of allowed values. Without those conditions, the differential equation alone would not determine the same spectrum.

For a particle of mass mm in an infinite square well, this becomes

En=ℏ22mλn=n2π2ℏ22mL2.E_n = \frac{\hbar^2}{2m}\lambda_n = \frac{n^2\pi^2\hbar^2}{2mL^2}.

Some problems have the form

Au=λBu,A u=\lambda B u,

where BB is not the identity. In Sturm–Liouville theory, the standard form is

−ddx[p(x)dudx]+q(x)u=λw(x)u.-\frac{d}{dx} \left[p(x)\frac{du}{dx}\right] +q(x)u = \lambda w(x)u.

The weight w(x)w(x) changes the inner product used for orthogonality:

⟨f,g⟩w=∫abf(x)∗g(x)w(x) dx.\langle f,g\rangle_w = \int_a^b f(x)^*g(x)w(x)\,dx.

The specialized theory, including weight functions and completeness intuition, is Sturm–Liouville Theory.

Many multidimensional eigenvalue problems first become one-dimensional eigenvalue problems by Separation of Variables. Green Functions package the inverse or resolvent of an eigenvalue problem and make poles, boundary conditions, and response to sources explicit.

In finite dimension, eigenvalues form a finite set. In differential-operator problems, several spectral patterns can occur.

A discrete spectrum has isolated eigenvalues with normalizable eigenfunctions. Bound states in an infinite square well or harmonic oscillator are standard examples:

Hψn=Enψn,∥ψn∥=1.H\psi_n=E_n\psi_n, \qquad \lVert\psi_n\rVert=1.

A continuous spectrum is not usually described by normalizable eigenfunctions. For a free particle on the line, the formal plane waves

eikxe^{ikx}

satisfy the eigenvalue equation for kinetic energy, but they are not in L2(R)L^2(\mathbb R). They are best understood through spectral representations, wave packets, or generalized eigenfunctions.

Many Hamiltonians have mixed spectra: bound states at isolated energies and scattering states in a continuum. Continuous Spectra explains why continuum labels require more care than a large sum over ordinary eigenvectors.

In quantum mechanics, observables are represented by self-adjoint operators. For eigenvalue problems this matters because self-adjointness gives:

  • real spectral values;
  • orthogonality of eigenfunctions associated with distinct discrete eigenvalues;
  • unitary time evolution for Hamiltonians;
  • the projection-valued spectral structure used for probabilities.

Formal symmetry is not enough in infinite-dimensional problems. Boundary conditions can make a differential expression self-adjoint, fail to make it self-adjoint, or define different self-adjoint operators with different spectra. The distinction is developed in Symmetric versus Self-Adjoint Operators.

An eigenvalue is degenerate when the eigenspace has dimension greater than one. In finite dimensions, degeneracy means

dim⁡ker⁡(A−λI)>1.\dim\ker(A-\lambda I)>1.

In quantum mechanics, degeneracy often reflects symmetry. For example, rotational symmetry can make several angular-momentum states share one energy. Degenerate eigenspaces should be treated as subspaces, not as a single preferred eigenvector. The correct measurement object is usually the projector onto the whole eigenspace.

For a practical eigenvalue problem:

  1. Specify the space, operator, domain, and boundary conditions.
  2. Solve the differential or algebraic equation for a general spectral parameter.
  3. Impose boundary, matching, regularity, and normalizability conditions.
  4. Identify the allowed spectral values and eigenspaces.
  5. Normalize normalizable eigenfunctions using the correct inner product.
  6. Decide whether the spectrum is discrete, continuous, or mixed.
  7. Use the spectral decomposition or spectral theorem appropriate to the problem.

For numerical work, a discretized differential operator turns the problem into a matrix eigenvalue problem, but the discretization should preserve the essential boundary conditions and self-adjoint structure. See Matrix Diagonalization.

  • Solving only the differential equation and treating every value of λ\lambda as allowed.
  • Forgetting that the zero function is not an eigenfunction.
  • Confusing a formal plane wave with a normalizable eigenstate.
  • Applying finite-dimensional diagonalization intuition to continuous spectra without the spectral theorem.
  • Changing boundary conditions without recognizing that the operator and spectrum have changed.
  • Ignoring degeneracy and summing over arbitrary eigenvectors when a projector onto the eigenspace is the invariant object.
  • Treating self-adjointness as a technicality rather than part of the physical definition of an observable.
  • G. Teschl, Mathematical Methods in Quantum Mechanics, 2nd ed., American Mathematical Society, 2014.
  • B. C. Hall, Quantum Theory for Mathematicians, Springer, 2013.
  • G. B. Arfken, H. J. Weber, and F. E. Harris, Mathematical Methods for Physicists, 7th ed., Academic Press, 2013.
  • S. Axler, Linear Algebra Done Right, 3rd ed., Springer, 2015.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  1. For
A=(1004),A= \begin{pmatrix} 1 & 0\\ 0 & 4 \end{pmatrix},

find the eigenvalues and eigenspaces.

Solution

The eigenvalues are 11 and 44. The eigenspace for λ=1\lambda=1 is the span of (1,0)T(1,0)^T, and the eigenspace for λ=4\lambda=4 is the span of (0,1)T(0,1)^T.

  1. Solve the eigenvalue problem
−u′′=λu,u(0)=u(L)=0.-u''=\lambda u, \qquad u(0)=u(L)=0.

Why is λ=0\lambda=0 not an allowed eigenvalue?

Solution

For λ=k2>0\lambda=k^2>0, the solution satisfying u(0)=0u(0)=0 is u(x)=Asin⁡(kx)u(x)=A\sin(kx). The condition u(L)=0u(L)=0 requires kL=nπkL=n\pi, so

λn=n2π2L2,n=1,2,3,….\lambda_n=\frac{n^2\pi^2}{L^2}, \qquad n=1,2,3,\ldots.

For λ=0\lambda=0, the equation gives u(x)=Ax+Bu(x)=Ax+B. The two boundary conditions force A=B=0A=B=0, leaving only the zero function, which is not an eigenfunction.

  1. Explain why changing from Dirichlet to periodic boundary conditions changes the eigenvalue problem, even if the differential expression is still −d2/dx2-d^2/dx^2.
Solution

The operator is not just the differential expression. Its domain includes the boundary conditions. Dirichlet conditions require u(0)=u(L)=0u(0)=u(L)=0, while periodic conditions identify endpoint values and derivatives. These domains define different self-adjoint operators and generally produce different spectra.

  1. Why are free-particle plane waves not ordinary normalized eigenvectors in L2(R)L^2(\mathbb R)?
Solution

The function eikxe^{ikx} has constant modulus, so

∫−∞∞∣eikx∣2 dx=∞.\int_{-\infty}^{\infty} \lvert e^{ikx}\rvert^2\,dx = \infty.

It is not square-integrable. Plane waves are useful as generalized eigenfunctions or as components in wave-packet and spectral representations, not as ordinary normalized Hilbert-space vectors.