World Library  
Flag as Inappropriate
Email this Article

Variational method (quantum mechanics)

Article Id: WHEBN0037487265
Reproduction Date:

Title: Variational method (quantum mechanics)  
Author: World Heritage Encyclopedia
Language: English
Subject: Variational Monte Carlo, Variational perturbation theory, Matrix product state, Theoretical chemistry, Computational chemistry
Publisher: World Heritage Encyclopedia

Variational method (quantum mechanics)

In quantum mechanics, the variational method is one way of finding approximations to the lowest energy eigenstate or ground state, and some excited states. This allows calculating approximate wavefunctions such as molecular orbitals.[1] The basis for this method is the variational principle.[2][3]

The method consists in choosing a "trial wavefunction" depending on one or more parameters, and finding the values of these parameters for which the expectation value of the energy is the lowest possible. The wavefunction obtained by fixing the parameters to such values is then an approximation to the ground state wavefunction, and the expectation value of the energy in that state is an upper bound to the ground state energy. The Hartree–Fock method and the Ritz method both apply the variational method. The Harris functional method is anti-variational (it is a lower bound to the energy).


Suppose we are given a Hilbert space and a Hermitian operator over it called the Hamiltonian H. Ignoring complications about continuous spectra, we look at the discrete spectrum of H and the corresponding eigenspaces of each eigenvalue λ (see spectral theorem for Hermitian operators for the mathematical background):

\sum_{\lambda_1,\,\lambda_2\in\mathrm{Spec}(H)}\lang\psi_{\lambda_1}\mid \psi_{\lambda_2}\rang=\delta_{\lambda_1\lambda_2}

where \delta_{i,j} is the Kronecker delta

\hat{H} \left| \psi_\lambda\right\rangle = \lambda\left|\psi_\lambda \right\rangle.

Physical states are normalized, meaning that their norm is equal to 1. Once again ignoring complications involved with a continuous spectrum of H, suppose it is bounded from below and that its greatest lower bound is E0. Suppose also that we know the corresponding state |ψ>. The expectation value of H is then

\begin{align} \left\langle\psi\mid H\mid \psi\right\rangle & = \sum_{\lambda_1,\lambda_2 \in \mathrm{Spec}(H)} \left\langle\psi|\psi_{\lambda_1}\right\rangle \left\langle\psi_{\lambda_1}|H|\psi_{\lambda_2}\right\rangle \left\langle\psi_{\lambda_2}|\psi\right\rangle \\ & =\sum_{\lambda\in \mathrm{Spec}(H)}\lambda \left|\left\langle\psi_\lambda\mid \psi\right\rangle\right|^2\ge\sum_{\lambda \in \mathrm{Spec}(H)}E_0 \left|\left\langle\psi_\lambda\mid \psi\right\rangle\right|^2=E_0 \end{align}

Obviously, if we were to vary over all possible states with norm 1 trying to minimize the expectation value of H, the lowest value would be E0 and the corresponding state would be an eigenstate of E0. Varying over the entire Hilbert space is usually too complicated for physical calculations, and a subspace of the entire Hilbert space is chosen, parametrized by some (real) differentiable parameters αi (i = 1, 2, ..., N). The choice of the subspace is called the ansatz. Some choices of ansatzes lead to better approximations than others, therefore the choice of ansatz is important.

Let's assume there is some overlap between the ansatz and the ground state (otherwise, it's a bad ansatz). We still wish to normalize the ansatz, so we have the constraints

\left\langle \psi(\alpha_i) \mid \psi(\alpha_i) \right\rangle = 1

and we wish to minimize

\varepsilon(\alpha_i) = \left\langle \psi(\alpha_i)|H|\psi(\alpha_i) \right\rangle.

This, in general, is not an easy task, since we are looking for a global minimum and finding the zeroes of the partial derivatives of ε over αi is not sufficient. If ψ (αi) is expressed as a linear combination of other functions (αi being the coefficients), as in the Ritz method, there is only one minimum and the problem is straightforward. There are other, non-linear methods, however, such as the Hartree–Fock method, that are also not characterized by a multitude of minima and are therefore comfortable in calculations.

There is an additional complication in the calculations described. As ε tends toward E0 in minimization calculations, there is no guarantee that the corresponding trial wavefunctions will tend to the actual wavefunction. This has been demonstrated by calculations using a modified harmonic oscillator as a model system, in which an exactly solvable system is approached using the variational method. A wavefunction different from the exact one is obtained by use of the method described above.

Although usually limited to calculations of the ground state energy, this method can be applied in certain cases to calculations of excited states as well. If the ground state wavefunction is known, either by the method of variation or by direct calculation, a subset of the Hilbert space can be chosen which is orthogonal to the ground state wavefunction.

\left| \psi \right\rangle = \left|\psi_{\text{test}}\right\rangle - \left\langle\psi_{\mathrm{gr}} \mid \psi_{\text{test}}\right\rangle \left|\psi_{\text{gr}}\right\rangle

The resulting minimum is usually not as accurate as for the ground state, as any difference between the true ground state and \psi_{\text{gr}} results in a lower excited energy. This defect is worsened with each higher excited state.

In another formulation:

E_\text{ground} \le \left\langle\phi|H|\phi\right\rangle.

This holds for any trial φ since, by definition, the ground state wavefunction has the lowest energy, and any trial wavefunction will have energy greater than or equal to it.

Proof: φ can be expanded as a linear combination of the actual eigenfunctions of the Hamiltonian (which we assume to be normalized and orthogonal):

\phi = \sum_n c_n \psi_n. \,

Then, to find the expectation value of the Hamiltonian:

\begin{align} & \left\langle\phi|H|\phi\right\rangle \\ & = \left\langle\sum_n c_n \psi_n |H|\sum_m c_m\psi_m\right\rangle \\ & = \sum_n\sum_m \left\langle c_n \psi_{n}|E_m|c_m\psi_m\right\rangle \\ & = \sum_n\sum_m c_n^*c_m E_m\left\langle\psi_n\mid\psi_m\right\rangle \\ & = \sum_{n} |c_n|^2 E_n. \end{align}

Now, the ground state energy is the lowest energy possible, i.e. E_{n} \ge E_{g}. Therefore, if the guessed wave function φ is normalized:

\left\langle\phi|H|\phi\right\rangle \ge E_g \sum_n |c_n|^2 = E_g. \,

In general

For a hamiltonian H that describes the studied system and any normalizable function Ψ with arguments appropriate for the unknown wave function of the system, we define the functional

\varepsilon\left[\Psi\right] = \frac{\left\langle\Psi|\hat{H}|\Psi\right\rangle}{\left\langle\Psi \mid \Psi\right\rangle}.

The variational principle states that

  • \varepsilon \geq E_0, where E_0 is the lowest energy eigenstate (ground state) of the hamiltonian
  • \varepsilon = E_0 if and only if \Psi is exactly equal to the wave function of the ground state of the studied system.

The variational principle formulated above is the basis of the variational method used in quantum mechanics and quantum chemistry to find approximations to the ground state.

The Harris functional method is anti-variational (it is a lower bound on the energy).

Another facet in variational principles in quantum mechanics is that since \Psi and \Psi^\dagger can be varied separately (a fact arising due to the complex nature of the wave function), the quantities can be varied in principle just one at a time.[4]

Helium atom ground state

The helium atom consists of two electrons with mass m and electric charge −e, around an essentially fixed nucleus of mass Mm and charge +2e. The Hamiltonian for it, neglecting the fine structure, is:

H = -\frac{\hbar^2}{2m} (\nabla_1^2 + \nabla_2^2) - \frac{e^2}{4\pi\epsilon_0} \left(\frac{2}{r_1} + \frac{2}{r_2} - \frac{1}{|\mathbf{r}_1 - \mathbf{r}_2|}\right)

where ħ is the reduced Planck constant, ε0 is the vacuum permittivity, ri (for i = 1, 2) is the distance of the i-th electron from the nucleus, and |r1 − r2| is the distance between the two electrons.

If the term Vee = e2/(4πε0|r1 − r2|), representing the repulsion between the two electrons, were excluded, the Hamiltonian would become the sum of two hydrogen-like atom Hamiltonians with nuclear charge +2e. The ground state energy would then be 8E1 = −109 eV, where E1 is the Rydberg constant, and its ground state wavefunction would be the product of two wavefunctions for the ground state of hydrogen-like atoms:[5]

\psi(\mathbf{r}_1,\mathbf{r}_2) = \frac{Z^3}{\pi a_0^3} e^{-Z(r_1+r_2)/a_0}.

where a0 is the Bohr radius and Z = 2, helium's nuclear charge. The expectation value of the total Hamiltonian H (including the term Vee) in the state described by ψ0 will be an upper bound for its ground state energy. <Vee> is −5E1/2 = 34 eV, so is 8E1 − 5E1/2 = −75 eV.

A tighter upper bound can be found by using a better trial wavefunction with 'tunable' parameters. Each electron can be thought to see the nuclear charge partially "shielded" by the other electron, so we can use a trial wavefunction equal with an "effective" nuclear charge Z < 2: The expectation value of H in this state is:

\langle H \rangle = \left[-2Z^2 + \frac{27}{4}Z\right] E_1

This is minimal for Z = 27/16; Shielding reduces the effective charge to ~1.69. Substituting this value of Z into the expression for H yields 729E1/128 = −77.5 eV, within 2% of the experimental value, −78.975 eV.[6]

Even closer estimations of this energy have been found using more complicated trial wave functions with more parameters. This is done in physical chemistry via Variational Monte Carlo.


  1. ^ Lorentz Trial Function for the Hydrogen Atom: A Simple, Elegant Exercise Thomas Sommerfeld Journal of Chemical Education 2011 88 (11), 1521–1524 doi:10.1021/ed200040e
  2. ^  
  3. ^  
  4. ^ see Landau, Quantum Mechanics , pg. 58 for some elaboration.
  5. ^ Griffiths (1995), p. 262.
  6. ^ G.W.F. Drake and Zong-Chao Van (1994). "Variational eigenvalues for the S states of helium", Chem. Phys. Lett. 229 486–490. [1]
This article was sourced from Creative Commons Attribution-ShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for and content contributors is made possible from the U.S. Congress, E-Government Act of 2002.
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a non-profit organization.

Copyright © World Library Foundation. All rights reserved. eBooks from World eBook Library are sponsored by the World Library Foundation,
a 501c(4) Member's Support Non-Profit Organization, and is NOT affiliated with any governmental agency or department.