physics2u
Depth
PU-104 / L22 · lecture

The Min-Max Principle and Eigenvalue Bounds

~7 min read · Tier toggle changes depth (T1 plain · T2 standard · T3 rigour)

BEFORE YOU START

You should be comfortable with the definition of a self-adjoint operator, the concept of an inner product in a Hilbert space, and the basic properties of eigenvalues and eigenvectors for finite-dimensional matrices.

Why this matters

The min-max principle provides a rigorous framework for bounding eigenvalues without solving the full eigenvalue problem.

In quantum mechanics and classical normal-mode analysis, we often cannot find exact analytical solutions for complex systems. However, we frequently need to know the ground state energy or the frequency of a specific mode. The min-max principle allows us to construct upper and lower bounds for these values using trial functions or subspaces, providing a powerful tool for approximation and stability analysis that is far more robust than simple perturbation theory.

The Core Idea

Imagine you are trying to find the lowest point in a landscape, but you are not allowed to look at the whole map. You can only walk on a specific path. The lowest point you find on that path is an upper bound for the true lowest point of the entire landscape. If you restrict your search to a larger area, your best guess can only get better (lower). Now, imagine you want to find the second lowest point. You must ignore the absolute lowest point. The min-max principle formalizes this intuition: to find the \(n\)-th eigenvalue, you look at all possible \(n\)-dimensional subspaces, and within each, you find the maximum possible Rayleigh quotient. The minimum of these maxima gives you the exact \(n\)-th eigenvalue. It is a game of "best worst-case scenario" that yields a precise mathematical result.

Formal Statement of the Principle

Let \(\hat{A}\) be a self-adjoint operator on a Hilbert space \(\mathcal{H}\) with a discrete, bounded-below spectrum. Let the eigenvalues be ordered such that \(\lambda_1 \le \lambda_2 \le \lambda_3 \le \cdots\). The Rayleigh quotient for a non-zero vector \(\psi\) is defined as \(R(\psi) = \dfrac{\langle \psi, \hat{A}\psi\rangle}{\langle \psi,\psi\rangle}\). The Courant–Fischer min–max principle states that the \(n\)-th eigenvalue can be characterized in two equivalent ways: as the minimum over all \(n\)-dimensional subspaces of the maximum Rayleigh quotient within that subspace, or as the maximum over all \((n-1)\)-dimensional subspaces of the minimum Rayleigh quotient for vectors orthogonal to that subspace.

\[ \lambda_n = \min_{\dim S = n}\ \max_{\psi \in S \setminus \{0\}} R(\psi) = \max_{\dim T = n-1}\ \min_{\psi\perp T} R(\psi). \]

This duality is crucial. The first expression tells us that if we pick any \(n\)-dimensional subspace \(S\), the largest eigenvalue of the projection of \(\hat{A}\) onto \(S\) is at least \(\lambda_n\). The second expression tells us that if we remove any \((n-1)\)-dimensional subspace \(T\), the smallest eigenvalue of \(\hat{A}\) in the remaining orthogonal complement is at most \(\lambda_n\).

Deriving Bounds for Approximations

In practical applications, we often approximate the Hilbert space with a finite-dimensional subspace \(V\) of dimension \(N\). If we compute the eigenvalues \(\mu_1 \le \mu_2 \le \cdots \le \mu_N\) of the operator restricted to \(V\), the min-max principle guarantees that \(\mu_k \ge \lambda_k\) for all \(k \le N\). This is the basis of the Rayleigh-Ritz method. The approximation provides an upper bound for the true eigenvalues. To improve the bound for \(\lambda_1\), we simply need to choose a subspace \(V\) that contains a vector close to the true ground state eigenvector. The quality of the bound depends entirely on how well the subspace captures the relevant physics.

For the \(n\)-th eigenvalue, the condition is stricter. We need a subspace of dimension \(n\) that is "aligned" with the first \(n\) eigenvectors. If the trial subspace is exactly the span of the first \(n\) true eigenvectors, the maximum Rayleigh quotient in it is exactly \(\lambda_n\). Any other \(n\)-dimensional subspace gives a maximum greater than or equal to \(\lambda_n\), still a valid upper bound.

Connection to Variational Methods

The min-max principle generalizes the simple variational principle for the ground state, which states that \(\lambda_1 = \min_{\psi \neq 0} R(\psi)\). While the ground state energy is easy to bound from above, excited states are notoriously difficult because the variational principle alone does not distinguish between the ground state and excited states. The min-max principle resolves this by introducing the dimension of the subspace as a parameter. It ensures that we are not accidentally finding the ground state again when looking for the first excited state. This is why in quantum chemistry, finding the first excited state often requires imposing orthogonality constraints to the ground state wavefunction, effectively applying the second part of the min-max formula.

The derivation

◆ DERIVATION EMBED · D-086 · not re-derived here

Variational Characterization of Eigenvalues

Let \(\hat{A}\) be a self-adjoint operator on a Hilbert space \(\mathcal{H}\) with a discrete, bounded-below spectrum and eigenvalues ordered \(\lambda_1 \le \lambda_2 \le \lambda_3 \le \cdots\). Define the Rayleigh quotient \(R(\psi) = \dfrac{\langle \psi, \hat{A}\psi\rangle}{\langle \psi,\psi\rangle}\) for \(\psi \neq 0\). Then the stationary points of \(R\) are exactly the eigenvectors of \(\hat{A}\), the stationary values are the eigenvalues, and the ordered eigenvalues obey the Courant–Fischer min–max principle \(\lambda_n = \min_{\dim S = n}\ \max_{\psi \in S \setminus \{0\}} R(\psi)\).

\[ \lambda_1 = \min_{\psi\neq 0} R(\psi), \qquad \lambda_n = \min_{\dim S = n}\ \max_{\psi\in S\setminus\{0\}} \frac{\langle\psi,\hat{A}\psi\rangle}{\langle\psi,\psi\rangle} = \max_{\dim T = n-1}\ \min_{\psi\perp T} R(\psi). \]
Open full derivation →
Tier 3 · watch out

A common pitfall is assuming that if a trial function is orthogonal to the first \(n-1\) approximate eigenvectors, it will automatically yield a good bound for \(\lambda_n\). This is only true if the approximate eigenvectors are themselves good approximations to the true eigenvectors. If the lower-lying states are poorly approximated, the orthogonality constraint may force the trial function into a region of high energy, leading to a loose upper bound that is far from \(\lambda_n\). The accuracy of the bound for \(\lambda_n\) is limited by the accuracy of the bounds for \(\lambda_1, \dots, \lambda_{n-1}\). C

Worked example

Problem. Consider a 3x3 symmetric matrix \(A\) representing a self-adjoint operator in a 3-dimensional space. The eigenvalues of \(A\) are known to be \(\lambda_1 = 1\), \(\lambda_2 = 4\), and \(\lambda_3 = 9\). We construct a 2-dimensional subspace \(S\) spanned by the vectors \(\mathbf{v}_1 = (1, 0, 0)\) and \(\mathbf{v}_2 = (0, 1, 1)\). Determine the upper bound for \(\lambda_2\) provided by this subspace using the min-max principle.

  1. Identify the subspace \(S\). The dimension is \(n=2\). According to the min-max principle, \(\lambda_2 \le \max_{\psi \in S \setminus \{0\}} R(\psi)\).
  2. Find the eigenvalues of the projection of \(A\) onto \(S\). First, we need the matrix representation of \(A\) in the standard basis. Since the eigenvalues are 1, 4, 9, and assuming the standard basis aligns with the eigenvectors for simplicity (or we can just use the Rayleigh quotient directly on the subspace), let's assume the eigenvectors are the standard basis vectors \(\mathbf{e}_1, \mathbf{e}_2, \mathbf{e}_3\). Then \(A = \text{diag}(1, 4, 9)\).
  3. Any vector in \(S\) is of the form \(\psi = a \mathbf{v}_1 + b \mathbf{v}_2 = (a, b, b)\). Calculate \(R(\psi) = \frac{\langle \psi, A \psi \rangle}{\langle \psi, \psi \rangle}\).
  4. \(\langle \psi, A \psi \rangle = a^2(1) + b^2(4) + b^2(9) = a^2 + 13b^2\).
  5. \(\langle \psi, \psi \rangle = a^2 + b^2 + b^2 = a^2 + 2b^2\).
  6. Maximize \(f(a,b) = \frac{a^2 + 13b^2}{a^2 + 2b^2}\). Let \(x = a^2\) and \(y = b^2\). We maximize \(g(x,y) = \frac{x + 13y}{x + 2y}\). Taking the derivative with respect to \(x\) and setting to zero, or observing that the function is monotonic in \(x/y\), we see that as \(x \to 0\) (i.e., \(a \to 0\)), the value approaches \(\frac{13y}{2y} = 6.5\). As \(x \to \infty\), it approaches 1. The maximum is 6.5.
  7. Thus, the upper bound for \(\lambda_2\) is 6.5. Since the true \(\lambda_2 = 4\), the bound is valid but loose.
\[ \lambda_2 \le 6.5 \]

Units check. The Rayleigh quotient is dimensionless if \(A\) is dimensionless, or has the same units as \(A\). Here, the numbers are pure scalars, so the result is consistent.

CHECK YOURSELF

Three quick tests

  1. Does the min-max principle provide lower or upper bounds for eigenvalues when using a finite-dimensional subspace?
    AnswerUpper bounds. The maximum Rayleigh quotient in a subspace is always greater than or equal to the true eigenvalue.
  2. If you want to bound \(\lambda_3\), what is the minimum dimension of the subspace \(S\) you must consider?
    AnswerDimension 3. The principle requires \(\dim S = n\) to bound \(\lambda_n\).
  3. Why is finding excited states harder than finding the ground state using variational methods?
    AnswerBecause the variational principle minimizes the energy, it naturally finds the ground state. To find excited states, you must impose orthogonality constraints to exclude lower states, which is less robust and more sensitive to errors in the lower states.
PROBLEMS

Problem set for PU-104 →

Each derivation in this unit ends with five graded problems with full solutions.

NOTES

Lecture notes drafted with a local model and checked line by line by the editor; every number in the worked example was recomputed.

Derivations cited: D-086 · Threads: symmetry, waves, matter