Three-Parameter Variational Wavefunctions for the Helium Isoelectronic Series

Initializing live version
Download to Desktop

Requires a Wolfram Notebook System

Interact on desktop, mobile and cloud with the free Wolfram Player or other Wolfram Language products.

An optimized three-parameter variational wavefunction for the helium isoelectronic series , , , …, can give quite accurate values for the ground-state energies. This computation is based on a simplified version of Hylleraas' classic work on two-electron atoms [1]. The Schrödinger equation for a two-electron atom with nuclear charge can be written as



with the Hamiltonian

expressed in atomic units . For S-states, including the ground state, the wavefunction depends on just three variables: , , and . Hylleraas used these alternative variables: , , and . Consider simple trial functions of the form . We need the combination in order that is an even function of , thus a symmetric function with respect to the interchange of and .

According to the variational principle, the ground-state energy can be approximated by , where


The parameters , , and are then selected such that takes a minimum value. This is always an upper bound on the exact ground-state energy . To run this Demonstration, choose a value of corresponding to one of the two-electron atoms or ions. Then, vary the sliders determining various combinations of α, β, and γ. The ionization potential () is the energy to remove one electron, which leaves the atom in the one-electron energy state . The is expressed in (). As the blue and red areas become larger, exact values of the energy and ionization potential are being approached more closely. Also shown is the ratio of expectation values of potential and kinetic energies, . By the virial theorem for Coulombic systems, this ratio should equal . This value would be obtained by the exact solution of the Schrödinger equation, but also for an approximate solution with optimally chosen variational parameters.

If you want to look at the answers, the optimal values of , , and are tabulated in the Details section.


Contributed by: S. M. Blinder (April 2014)
Open content licensed under CC BY-NC-SA



The calculated optimal values are given in the following table.


[1] H. A. Bethe and E. E. Salpeter, Quantum Mechanics of One- and Two-Electron Atoms, Berlin: Springer–Verlag, 1957 pp. 146–154.

Feedback (field required)
Email (field required) Name
Occupation Organization
Note: Your message & contact information may be shared with the author of any specific Demonstration for which you give feedback.