WOLFRAM|DEMONSTRATIONS PROJECT

Simplified Hartree-Fock Computations on Second-Row Atoms

​
atomic number Z
2
3
4
5
6
7
8
9
10
α
0.501
β
0.501
γ
0.501
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
N
occupied spin-orbitals:
Ψ(1⋯N)=
1
N!

ϕ
1
(1)
ϕ
2
(1)
⋮
ϕ
N
(1)
ϕ
1
(2)
ϕ
2
(2)
⋮
ϕ
N
(2)
⋯
⋯
⋱
⋯
ϕ
1
(N)
ϕ
2
(N)
⋮
ϕ
N
(N)

,
where spin-orbitals are products of the form
ϕ(x)=ψ(r)
α
{
β
.
A Slater determinant automatically satisfies the Pauli exclusion principle, requiring that spin-orbitals be, at most, singly occupied. The Hamiltonian for an
N
-electron atom, in hartree atomic units, is given by:
H=
N
Σ
i=1
-
1
2
2
∇
-
Z
r
i
+
N
Σ
i>j=1
1
r
ij
.
The corresponding approximation to the total energy is then given by
E=
N
Σ
i=1
H
i
+
N
Σ
i>j=1

J
ij
-
K
ij

,
where
H
,
J
and
K
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,
Z=2
to 10. Assume a single Slater determinant, thus possibly shortchanging open-shell configurations. The computation involves only
1s
,
2s
and
2p
orbitals, approximated by modified Slater-type orbitals:
ψ
1s
=
3/2
α
π
-αr
e
,
ψ
2s
=
3
5
β
π
2
α
-αβ+
2
β

1-
α+β
3
r
-βr
e
,
ψ
2p{x,y,z}
=
5/2
γ
π
r
-γr
e
{sinθcosϕ,sinθsinϕ,cosθ}
.
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.