Exact solutions to the Schrödinger equation are rare treasures in physics, restricted to simple systems like the Hydrogen atom or the Infinite Square Well. For the vast majority of real-world quantum systems—from multi-electron atoms to solid-state materials—exact analytical solutions are impossible. The Variational Method is a powerful approximation technique that allows us to estimate the ground state energy of these complex systems using a rigorous inequality rather than guesswork.
Figure 1: Conceptual illustration of the Variational Principle. The calculated energy \( E(\alpha) \) (blue curve) is always greater than the true ground state energy \( E_0 \) (red dashed line). The minimum point represents the best approximation.
The Variational Method: Approximation by Inequality
1. Mathematical Setup
Consider a quantum system described by a Hamiltonian \( \hat{H} \). We assume this system has a set of discrete bound states with energies ordered as:
\[ E_0 < E_1 < E_2 < \dots \]
Here, \( E_0 \) is the true ground state energy, which is often unknown. We introduce a normalized trial wavefunction \( |\psi\rangle \) that satisfies the boundary conditions of the system.
We define the Energy Functional as the expectation value of the Hamiltonian in this trial state:
2. The Variational Theorem
The fundamental theorem of this method states that the expectation value of the Hamiltonian is an upper bound to the true ground state energy.
\[ \langle \psi | \hat{H} | \psi \rangle \ge E_0 \]
Equality holds if and only if the trial state is exactly the ground state (\( |\psi\rangle = |0\rangle \)). This theorem ensures that any approximate energy we calculate will never be lower than the true physical energy.
3. Proof of the Theorem
The proof relies on the completeness of the eigenstates of the Hamiltonian.
Step 1: Expansion
Since the eigenstates \( |n\rangle \) form a complete basis, we can expand our arbitrary trial state \( |\psi\rangle \) as:
\[ |\psi\rangle = \sum_n c_n |n\rangle \]
where \( c_n = \langle n | \psi \rangle \).
Step 2: Normalization
The normalization condition implies:
\[ \langle \psi | \psi \rangle = \sum_n |c_n|^2 = 1 \]
Step 3: Expectation Value
Now we calculate the energy expectation value:
\[ \langle \psi | \hat{H} | \psi \rangle = \left( \sum_m c_m^* \langle m| \right) \hat{H} \left( \sum_n c_n |n\rangle \right) \]
Using \( \hat{H}|n\rangle = E_n|n\rangle \) and the orthogonality \( \langle m|n \rangle = \delta_{mn} \):
\[ \langle \psi | \hat{H} | \psi \rangle = \sum_n |c_n|^2 E_n \]
Step 4: The Inequality
Since \( E_0 \) is the lowest energy, we know that \( E_n \ge E_0 \) for all \( n \). Therefore:
\[ \sum_n |c_n|^2 E_n \ge \sum_n |c_n|^2 E_0 \]
Factoring out \( E_0 \):
\[ \langle \psi | \hat{H} | \psi \rangle \ge E_0 \left( \sum_n |c_n|^2 \right) \]
Since the sum is 1, we prove:
\[ \langle \psi | \hat{H} | \psi \rangle \ge E_0 \]
4. Method of Optimization
To make this theorem useful, we choose a family of trial functions \( \psi(\mathbf{r}; \alpha_1, \alpha_2, \dots) \) that depend on adjustable parameters \( \alpha_i \).
- Compute the energy expectation value as a function of the parameters: \( E(\alpha) \).
- Minimize \( E(\alpha) \) with respect to these parameters: \[ \frac{\partial E}{\partial \alpha_i} = 0 \]
- The resulting \( E_{min} \) is the best possible approximation for the ground state energy using that specific trial function.
5. Example: The Helium Atom
The variational method is famously used to solve the Helium atom, where the Hamiltonian includes a repulsion term between the two electrons: \[ \hat{H} = -\frac{\hbar^2}{2m}(\nabla_1^2 + \nabla_2^2) - \frac{2e^2}{r_1} - \frac{2e^2}{r_2} + \frac{e^2}{|\mathbf{r}_1 - \mathbf{r}_2|} \]
The Trial Function: Instead of assuming the full nuclear charge \( Z=2 \), we treat the charge as a variational parameter \( Z_{eff} < 2 \). This accounts for screening, where one electron "shields" the other from the full pull of the nucleus.
Figure 2: Representation of the Helium atom. The electron-electron repulsion makes exact solution impossible, but variational parameters allow us to model the effective nuclear charge.
Result: Minimizing the energy with respect to \( Z_{eff} \) yields \( Z_{eff} \approx 1.69 \). This simple variational calculation brings the theoretical energy within roughly 2% of the experimental value, a massive improvement over standard perturbation theory.
6. Conclusion
The Variational Method connects quantum mechanics to optimization theory. It tells us that nature "prefers" the state of lowest energy, and by mathematically searching for that minimum, we can approximate the behavior of complex systems like molecules and solids. Modern techniques like Density Functional Theory (DFT) and Hartree-Fock are essentially sophisticated extensions of this variational principle.