Simplified Hartree-Fock Computations on Second-Row Atoms

Modern computational quantum chemistry has developed largely from applications of the Hartree–Fock method to atoms and molecules [1–3]. A simple representation of a many-electron atom is given by a Slater determinant constructed from occupied spin-orbitals:
where spin-orbitals are products of the form
A Slater determinant automatically satisfies the Pauli exclusion principle, requiring that spin-orbitals be, at most, singly occupied. The Hamiltonian for an -electron atom, in hartree atomic units, is given by:
The corresponding approximation to the total energy is then given by
where , and are, respectively, the core, coulomb and exchange integrals.
This Demonstration carries out a simplified Hartree–Fock computation on the ground states of the atoms He to Ne, to 10. Assume a single Slater determinant, thus possibly shortchanging open-shell configurations. The computation involves only , and orbitals, approximated by modified Slater-type orbitals:
These are normalized and mutually orthogonal, which simplifies the computation. The orbital parameters are chosen so as to minimize the total energy, in accordance with the variational principle. You can vary these using the sliders. Actually, since the possible values of vary over inconveniently large ranges, the sliders vary internal parameters that determine them. The values of are shown in the graphic.
The blue bar shows the exact energy of the atomic ground state (actually, the exact nonrelativistic value) in hartrees. The red bar shows the calculated energy, given the displayed values of . You should try to adjust the three parameters to get the lowest possible energy. You can never reach the exact energy, however, since a single Hartree–Fock determinant fails to account for correlation energy, which involves instantaneous electron-electron interactions.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


Following are the optimized values of the energies and orbital parameters, found using the built-in Mathematica function FindMinimum:
Z = 2, {-2.84766, {α 1.6875}}
Z = 3, {-7.42617, {α
2.69716, β 0.721027}}
Z = 4, {-14.5413, {α
3.71867, β 1.10264}}
Z = 5, {-24.4442, {α
4.72562, β 1.50386, γ 1.2137}}
Z = 6, {-37.4687, {α
5.72987, β 1.90153, γ 1.53991}}
Z = 7, {-53.9194, {α
6.73255, β 2.29803, γ 1.86317}}
Z = 8, {-74.1007, {α
7.73422, β 2.69406, γ 2.18487}}
Z = 9, {-98.3167, {α
8.73521, β 3.08988, γ 2.50566}}
Z = 10, {-126.872, {α
9.73572, β 3.48558, γ 2.82585}}
These can be slightly improved by hand.
[1] S. M. Blinder, Introduction to Quantum Mechanics, Amsterdam: Elsevier, 2004 section 9.7.
[2] A. Szabo and N. S. Ostlund, Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory, Mineola, NY: Dover, 1996.
[3] S. M. Blinder, "Introduction to the Hartree–Fock Method," in Mathematical Physics in Theoretical Chemistry (S. M. Blinder and J. E. House, eds.), Elsevier, 2018, in press.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2017 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+