Helium Atom, Approximate Methods: 24th April 2008
Helium Atom, Approximate Methods: 24th April 2008
Helium Atom, Approximate Methods: 24th April 2008
The hydrogen atom wavefunctions and energies, we have seen, are deter-
mined as a combination of the various quantum ”dynamical” analogues of
classical motions (translation, vibration, rotation) and a central-force inter-
action (i.e, the Coulomb interaction between an electron and a nucleus).
Now, we consider the Helium atom and will see that due to the attendant
3-body problem for which we cannot determine a closed-form, first-principles
analytic solution, we will have to find recourse in approximate methods.
1
h̄2 2e2 1 1 e2
−∇21 − ∇22 ψ(r1 , r2 )− + + ψ(r1 , r2 ) = Eψ(r1 , r2 )
2me 4πo r1 r2 4πo |r2 − r1 |
The ∇2 terms represent the kinetic energy of the two electrons. The r11
and r12 terms represent the nucleus-electron Coulomb interaction. The last
term on the left hand side of the equation represents the electron-electron
repulsion taken as a Coulomb interaction based on the absolute value of the
electron-electron separation.
1 ∂ ∂ 1 ∂ ∂ 1 ∂2
r12 + sinθ1 +
r12 ∂r1 ∂r1 2
r1 sinθ1 ∂θ1 ∂θ1 r12 sin2 θ1 ∂φ21
• Solve the relative motion problem (separate out the center of mass
motion as we have seen earlier)
• Center of mass is assumed to be the nucleus; good approximation for
heavier nuclei
Ĥ = (KE)e1 + (KE)e2 + VN e1 + VN e2
h̄2 2e2 1 1
2 2
= −∇1 − ∇2 − +
2me 4πo r1 r2
! !
−h̄2 2 2e2 1 −h̄2 2 2e2 1
= ∇ − + ∇ −
2me 1 4πo r1 2me 2 4πo r2
= Ĥ 1 + Ĥ 2
2
Under the independent electron approximation, if we take the total He atom
wavefunction as a product of the individual electron wavefunctions (here ap-
proxmiated as hydrogen-like wavefunctions):
• Ψnlm (r) = Rnl (Zr/ao )Ylm (θ, φ) (Z=1 for hydrogen, Z=2 for helium)
2
• Energies are: En = − Z 2Eh n12
He → He+ + e−
∆E = EHe,1s − EHe,(1s)2
= −54.4eV − (−108.8eV ) = 54.4eV
3
The ground state energy for Helium can also be contrasted as:
22 Eh 1 1
E 0 = E1s,1s = − 2
+ 2
2 1 1
= −4Eh
= −108.8eV signif icantly too negative
The experimental value for the ionization potential (from mass spectroscopic
measurements) is 24.6 eV; so the independent electron approach entails sig-
nifican error. Because the electrons are not allowed Coulombic repulsion,
the energy required to remove a particular electron is higher than the ex-
perimental value. This is also validated by the much lower (and thus more
attractive/favorable) ground state energy predicted by the independent elec-
tron model compared to experiment.
h̄2 2e2 1 1 e2
2 2
Ĥ = −∇1 − ∇2 − + +
2me 4πo r1 r2 4πo |r2 − r1 |
= Ĥ 0 + Ĥ 1 Ĥ 0 ψ 0 = E 0 ψ 0
4
subscript ”n” refers to any energy level in general (perturbation theory is
capable of giving us energy levels and wavefunctions in addition to those
for the ground state; as we will see below, variational approaches give us
ground state information only):
Z
ψn = ψn0 + ψn1 ψn0∗ ψn1 dτ =0
all space
En = En0 + En1
Ĥ 0 + Ĥ 1 ψn0 + ψn1 = En0 + En1 ψn0 + ψn1
Ĥ 0 ψn0 + Ĥ 1 ψn0 + Ĥ 0 ψn1 + Ĥ 1 ψn1 = En0 ψn0 + En1 ψn0 + En0 ψn1 + En1 ψn1
The first terms on each side are equivalent. Last terms are approxi-
mated/assumed to be neglibly small–this is a perturbation energy
and wavefunction. This leaves:
Solve for ψn1 and En1 . Left multiply by ψn0∗ and integrate:
Z Z Z Z
dτ ψn0∗ Ĥ 1 ψn0 + dτ ψn0∗ Ĥ 0 ψn1 = dτ ψn0∗ En1 ψn0 + dτ ψn0∗ En0 ψn1
∗
Since Ĥ 0 is Hermitian, dτ ψn0∗ Ĥ 0 ψn1 = dτ Ĥ 0 ψn0 ψn1 .
R R
Z Z Z
dτ ψn0∗ Ĥ 1 ψn0 + dτ ψn1 En0 ψn0∗ = En0 dτ ψn0∗ ψn1 + En1
Solving for En1 , the fist order perturbative correction to the inde-
pendent electron energy level for helium gives:
Z
En1 = dτ ψn0∗ Ĥ 1 ψn0
5
Expilcitly, this means, En = En0 + En1 ; we have added a ”small” perturbative
correction to the reference independent electron energy for a given energy
level, n!
Now, what about the correction to the wavefunction? Let’s recall:
j6=n
Z Z Z Z
dτ ψk0∗ Ĥ 1 ψn0 + dτ ψk0∗ Ĥ 0 ψn1 = dτ ψk0∗ En1 ψn0 + dτ ψk0∗ En0 ψn1
Z Z Z Z
dτ ψk0∗ Ĥ 1 ψn0 + dτ ψk0∗ Ĥ 0 anj ψj0 = dτ ψk0∗ En1 ψn0 + dτ ψk0∗ En0 anj ψj0
X X
j6=n j6=n
Z Z Z Z
dτ ψk0∗ Ĥ 1 ψn0 anj Ej0 dτ ψk0∗ ψj0 En1 dτ ψk0∗ ψn0 En0 dτ ψk0∗ ψj0
X X
+ = + anj
j6=n j6=n
Z Z
ank Ek0 − En0 = En1 dτ ψk0∗ ψn0 − dτ ψk0∗ Ĥ 1 ψn0
Z
k=n → En1 = dτ ψn0∗ Ĥ 1 ψn0
dτ ψk0∗ Ĥ 1 ψn0 1
R
Ĥkn
k 6= n → ank = 0 ≡
En0 − Ek En0 − Ek0
6
Thus, the first order correction to the wavefuntion is:
j6=n
1 dτ ψj0∗ Ĥ 1 ψn0
! R !
Ĥjn
ψn0 ψj0 ψn0 ψj0
X X
= + = +
j6=n
En0 − Ej0 j6=n
En0 − Ej0
Z Z 1
Ĥjn
En2 dτ ψn0∗ Ĥ 1 ψn1 dτ ψn0∗ Ĥ 1 ψ0
X
= = 0 − E0 j
E
j6=n n j
1 Ĥ 1
X Ĥnj jn
=
j6=n
En0 − Ej0
Z
En1 = dτ ψn0∗ Ĥ 1 ψn0
7
!1/2 2
Z3
ψ 0 (r1 , r2 ) = φ1s (r1 )phi1s (r2 ) = e−Zr1 e−Zr2
π
1
Ĥ 1 =
r12
Z6 1
Z Z
En1 = dτ ψn0∗ Ĥ 1 ψn0
= 2 dτ e−Zr1 e−Zr2 e−Zr1 e−Zr2
π r12
Z 6 Z
e −2Zr 1 e −2Zr 2 5
= 2
dτ = ZEh
π r12 8
5 −11
E0 = E00 + E01 = −4Eh + ZEh = Eh = −74.8eV
8 4
• Experiment: E = -79.0 eV
• Perturbation theory results may be greater or less than the true energy,
unlike variational results
8
the coordinates of a single electron. This is the orbital approxima-
tion. The wavefunctions are written:
NOTE:
• Due to the electron-electron repulsion term in the He atom Hamilto-
nian not being spherically symmetric, the Schrodinger equation cannot
be solved analytically. Numerical methods are applied along with ap-
proximations.
• One such approxmiation is the neglect of full electron-electron cor-
relation (motion of electrons is correlated; electrons stay out of each
other’s way; correlation opposes repulsion, and thus contributes an
energy-lowering influence).
• Advanced electron correlation methods are well-developed and rou-
tinely applied.
−2e2 e2
Z
V1ef f (r1 ) = + φ∗2 (r2 ) φ2 (r2 )dr2
4πo |r1 | rπo |r1 − r2 |
9
This model is the Hartree model and serves as a starting point for treating
many-electron atoms.
• Ag atoms have a closed inner shell (no net angular momentum for
closed shells; symmetry); outer valence electron, 5s must contain all
magnetic moment.
• Note that this idea and picture of intrinsic particle spin is a convenient
classical way of accounting for experiment.
10
Introduce spin angular momentum.
h̄2 1
2 2
ŝ α = h̄ s(s + 1)α(σ) = + 1 α(σ)
2 2
h̄2 1
ŝ2 β = h̄2 s(s + 1)β(σ) = + 1 β(σ)
2 2
h̄
ŝz α = ms h̄α(σ) = α(σ)
2
h̄
ŝz β = ms h̄β(σ) = − β(σ)
2
Z Z
∗
α(σ) β(σ)dσ = β(σ)∗ α(σ)dσ = 0
Z Z
α(σ)∗ α(σ)dσ = β(σ)∗ β(σ)dσ = 1
11
If:
ψ 2 (r1 , σ1 , r2 , σ2 ) = ψ 2 (r2 , σ2 , r1 , σ1 )
ψ(r1 , σ1 , r2 , σ2 ) = ±ψ(r2 , σ2 , r1 , σ1 )
For fermions, the negative sign applies, and the wavefunction is said
to be antisymmetric. For bosons, the positive sign applies and the
wavefunctions are said to be symmetric. The nature of these par-
ticles is described by different statistics, from which the names arise
(Fermi-Dirac for the former, and Bose-Einstein for the latter).
Thus:
ψ(r1 , r2 ) = −ψ(r2 , r1 )
Now let’s consider a few actual systems, beginning with the Helium
atom. If we write an approximate He atom wavefunction as a product
of 1-electron wavefunctions with adjustable value of nuclear charge, Z̃,
we have:
Including the spin orbitals and invoking Pauli Exclusion, we can write:
1
Ψ(r1 , r2 ) = √ |1s(r1 )α(1)1s(r2 )β(2) − 1s(r2 )α(2)1s(r1 )β(1)|
2
12
Pauli exclusion tells us that for two electrons in the same spatial
orbital, the spins must be opposite (no two electrons share all four
quantum numbers). This avoids distinguishing between the 2 elec-
trons, and is antisymmetric with respect to interchange of identical
particles.
1
Ψ(r1 , r2 ) = √ |1s(r1 )α(1)1s(r2 )β(2) − 1s(r2 )α(2)1s(r1 )β(1)|
2
1 1s(r1 )α(1) 1s(r1 )β(1)
= √
2 1s(r2 )α(2) 1s(r2 )β(2)
For the Lithium atom, the Slater determinant for the ground state is:
or
13
1s(r1 )α(1) 1s(r1 )β(1) 2s(r1 )β(1)
1
Ψ(r1 , σ1 , r2 , σ2 , r3 , σ3 ) = √ 1s(r2 )α(2) 1s(r2 )β(2) 2s(r2 )β(2)
3! 1s(r3 )α(3) 1s(r3 )β(3) 2s(r3 )β(3)
Finally, consider the caes of ground state Helium. Here the two elec-
trons are in the same orbital and have different spin. Thus, the Slater
Determinant for this case can be separated into its spatial and spin
parts. The spin part is then the antisymmetric component of the wave-
function:
1
Ψ(r1 , σ1 , r2 , σ2 ) = √ 1s(r1 )1s(r2 ) [α(1)β(2) − α(2)β(1)]
2
1 α(1) β(1)
= √ 1s(r1 )1s(r2 )
2 α(2) β(2)
To summarize:
Ĥψo = Eo ψo
Z Z
ψo∗ Ĥψo dr = ψo∗ Eo ψo dr = Eo
14
ψo not known, so assume a trial funtion, χ. Determine the expectation
of the Hamiltonian, Ĥ and define it to be the Variational energy,
Evar .
χ∗ Ĥχdr
R Z
Evar = R ∗ = χ∗ Ĥχdr
χ χdr
X
χ= an ψn
n
Z
χ∗ χdr = |an |2 = 1
X
normalization
n
Z
χ∗ Ĥχdr = |an |2 En
X
Evar =
n
Now, subtract the ground state energy, whatever that may be, from
the approximate variational energy:
since all the individual energies must be equal to or greater than the
ground state energy! Thus, as an approximation, one can posit sev-
eral trial wavefunctions, determine the variational energy, and then
15
optimize in the space of wavefunctions. No doubt, this is not the
most efficient approach to determine a valid functional form, but such
information is generally obtained from chemical intuition.
The variational method also requires that variational parameters ap-
pear in the formulation so as to allow optimization, as the following
example shows.
h̄2 d kx2
Ĥ = − +
2µ dx2 2
h̄2 d kx2
Ĥ = − 2
+
2µ dx 2
Trial wavefunction:
1/4
γ 2 /2
χ= e−γx ; γ variable (variational parameter)
π
h̄2 d
1/2 Z !
γ kx2
Z
−γx2 /2 2 /2
Evar = χ∗ Ĥχdr = e − + e−γx dx
π 2µ dx2 2
h̄2 h̄2 γ k
1/2 Z ∞ ∞ 1/2 Z ∞
γ k γ
Z
−γx2 2 −γx2 2
=− γe dx + (γx) e dx + x2 e−γx dx = + ≥ Eo
2µ π −∞ −∞ 2 π −∞ 4µ 4γ
dEvar h̄2 k
= − =0
dγ 4µ 4γ 2
(kµ)1/2
γ= ≡α
h̄
This leads, with no surprise, to the exact solutions to the 1-D H.O.
16
1/4
α 2 /2 h̄ω
χ = ψo (x) = e−αx → Evar = = Eo
π 2
Now, returning to the Helium atom, we can attempt to use the varia-
tional principle to determine appropriate wavefunctions and energies.
Recall from earlier the Helium atom Hamiltonian operator:
!
h̄2 2e2 1 1 e2
(−∇21 − ∇22 ψ(r1 , r2 )− + + ψ(r1 , r2 ) = Eψ(r1 , r2 )
2me 4πo r1 r2 4πo |r2 − r1 |
1 1 1 1
ĤHe = − ∇21 + ∇22 − 2 + +
2 r1 r2 |r2 − r1 |
Z
Evar = χ∗ Ĥχdr
!3/2
1 Z̃
φZ̃
1s =√ e−Z̃r/ao (1 − electron H atom orbitals)
π ao
!1/2 !1/2
Z̃ 3 −Z̃r1 Z̃ 3 −Z̃r2
χ(r1 , r2 ) = e ao e ao
π π
17
Z̃ 6 −1 2 1 1 1
Z −Z̃r1 −Z̃r2
−Z̃r1 −Z̃r2
Evar = e ao e ao ∇1 + ∇22 − 2 + + e ao e ao dτ
π 2 r1 r2 r12
27
= Evar (Z̃) = Z̃ 2 − Z̃
8
dEvar (Z̃) 27 27
= 2Z̃ − =0 Z̃ = <2
dZ̃ 8 16
= −77.5eV
18
Since we don’t know yet what these orbitals may look like (though
we can guess that they will be much like the hydrogenic orbital wave-
functions) we can label these individual orbitals as φ HF (r1 ), φHF (r2 ).
These are the Hartree-Fock orbitals. We can also solve for the or-
bital energies in a variational manner; this suggests that the orbitals
we use will incorporate a variational parameters. In practice, it
is the radial portion of the Hartree-Fock orbitals that differ from the
standard hydrogenic orbitals. The radial components incorporate an
”effective nuclear charge” as opposed to the actual nuclear charge of
a many-electron atom. This effective charge attempts to capture the
shielding effect of other electrons on the particular electron in a given
orbital. The effective charge is usually referred to as ζ.
1
Z
U1ef f (r1 ) =e 2
φ∗ (r2 ) φ(r2 )dr2
r12
1 Z
Ĥ1ef f (r1 ) = − ∇21 − + U1ef f (r1 )
2 r1
Ĥ1ef f (r1 )φHF (r1 ) = 1 φHF (r1 )
19
from knowing anything about electron 2, since we need the distrubition
function (probability density) of electron 2 to construct the effective
potential for electron 1. Likewise, for a many-electron problem, we
would need information on all other electrons to solve this one-electron
problem. The approach thus requires a self-consistent approach as
outlined below.
−1 2 1 1 1
Z
HF HF ∗ HF ∗
E = φ (r1 )φ (r2 ) ∇1 + ∇22 − Z + + φHF (r1 )φHF (r2 )dr1 r. 2
2 r1 r2 r12
= I1 + I2 + J12
20
Several terms are evident:
" #
∇2 Z
Z
I1 = φ HF ∗
(r1 ) − 1 − φHF (r1 )dr1
2 r1
" #
∇2 Z
Z
I2 = φHF ∗ (r2 ) − 2 − φHF (r2 )dr2 = I1
2 r2
1
Z
J12 = φHF ∗ (r1 )φHF ∗ (r2 ) φHF (r1 )φHF (r2 )dr1 dr2 Coulomb Integral
r12
For the Helium Atom, the ground state energies based on HF approach
compare well to experiment:
EHF = −77.87eV
Eexpt = −79.0eV
= I1 + J12 = E HF − I2
21
To a good approximation, the association of orbital energies to ion-
ization energies provides a reasonable estimate for most systems. In
certain cases, there are systematic errors which can be dealt with via
correlation to experimental data. One must also bear in mind that
the orbital energies will vary based on the form of the wavefunctions
employed for the calculations. Correcting for electron correlation will
also lead to differences.
Regarding the relationship between the orbital energies and the to-
tal energy of a particular electronic configuration of an atom, we can
consider the sum of the two orbital energies, 1 and 2 for the He atom.
1 + 2 = I1 + I2 + 2J12 = E HF + J12
E HF = 1 + 2 − J12
22
Shell 1 1s
Shell 2 2s 2p
Shell 3 3s 3p 3d
Shell 4 4s 4p 4d 4f
Shell 5 5s 5p 5d 5f
Shell 6 6s 6p 6d 6f
Thus, we see that based on total energies, the Aufbau Principle pro-
vides a framework for ordering the available many-electron atom or-
bitals based on energy. We see that in the many-electron atom, the
s,p,d,f,g,... orbitals are not degenerate as they were for the Hydrogen
atom. This arises due to differing shielding of electrons in different or-
bitals based on the probabilities of finding electrons closer or further
away from the nucleus.
23