Variational method (source code)

= Variational method
{title2=$E_0\leq E_\psi$}
{wiki}

The <variational method> estimates the lowest <energy> of a <Quantum Hamiltonian> by choosing a trial <wavefunction>, evaluating its <energy>, and optimizing the trial parameters. It uses the <quantum variational principle>, the <Rayleigh-Ritz variational principle> expressed for <quantum states>. Let $H$ be a <self-adjoint operator> bounded below, and define $E_0=\inf\sigma(H)$, the bottom of its <spectrum>. Every normalized admissible trial <quantum state> satisfies
$$
E_\psi\geq E_0,\qquad \|\psi\|=1.
$$
When $E_0$ is an <energy eigenvalue>, it is the <ground-state energy>; equality holds precisely for a <quantum state> in the <ground-state subspace>. Otherwise the bottom of the <spectrum> need not be attained by any normalizable <wavefunction>.

Admissible means finite form energy. Choose a real $b\leq E_0$ and use the shifted <quadratic form of a positive quantum Hamiltonian>:
$$
Q(H)=D((H-bI)^{1/2}),\qquad q_H[\psi]=\|(H-bI)^{1/2}\psi\|^2+b\|\psi\|^2,\qquad E_\psi=q_H[\psi].
$$
The form domain $Q(H)$ is independent of this choice of $b$ and can be larger than the <operator domain> $D(H)$. For $\psi\in D(H)$, $E_\psi=\langle\psi,H\psi\rangle$ is the usual <expectation value>; for nonzero unnormalized trials, divide $q_H[\psi]$ by $\|\psi\|^2$ to get the form version of the <Rayleigh quotient>. See the <quantum variational principle> for the general domain justification.

Thus each trial <energy> is an \b[upper bound] on $E_0$, and minimizing within a trial family improves that upper bound. It is not a lower bound on $E_0$. The <Gaussian variational energy for the harmonic oscillator> reaches $E_0=1/2$ because its family contains an exact <ground state>. The <exact variational solution in a two-level system> makes the equality condition visible as vanishing excited-state <probability amplitude>.

The method is useful when solving the <Schrödinger equation> exactly is difficult but integrals of a trial <wavefunction> are tractable. Enlarging a trial family cannot increase its optimal <energy>, although an optimum need not be attained. The beauty is that even an imperfect trial gives a rigorous bound, provided its form energy is computed exactly: numerical integration errors require separate control. The <Gaussian variational estimate for a solvable sextic potential> provides an existing example beyond a <quantum harmonic oscillator>.

A short proof explains the direction. With a complete <orthonormal eigenbasis> of <energy eigenstates>, write $\psi=\sum_n c_n\phi_n$, where $H\phi_n=E_n\phi_n$ and $\sum_n|c_n|^2=1$. Then
$$
E_\psi-E_0=\sum_n|c_n|^2(E_n-E_0)\geq0.
$$
This is a <weighted mean> of nonnegative <energy> differences; equality requires every occupied <energy eigenstate> to have <energy> $E_0$. https://ocw.mit.edu/courses/8-05-quantum-physics-ii-fall-2013/61bc31b8d8bf0680c322733910a71aa0_MIT8_05F13_Chap_01.pdf[Barton Zwiebach, Wave Mechanics, MIT 8.05, 2013, §5, pp. 18–20]{ref}. For <continuous spectrum>, the <spectral theorem> replaces the sum by an integral against the <scalar spectral measure> of $\psi$, with the same nonnegative integrand. This proves the bound without assuming an <orthonormal eigenbasis> exists.

The historical setting predates <quantum mechanics>: Rayleigh's \i[The Theory of Sound], volume I, appeared in 1877, and Ritz's paper on solving variational problems of mathematical physics appeared in volume 135 of Crelle's journal, catalogued by its publisher as 1909, pp. 1–61. These are original works associated with the <Rayleigh-Ritz method>; the dates here identify the publications rather than assert priority. https://archive.org/details/theorysound06raylgoog[Lord Rayleigh, The Theory of Sound, volume I, 1877]{ref}; https://doi.org/10.1515/crll.1909.135.1[Walter Ritz, Über eine neue Methode zur Lösung gewisser Variationsprobleme der mathematischen Physik, original article and publisher metadata]{ref}.