Numerical Methods in Quantum Mechanics - Dipartimento di Fisica
Numerical Methods in Quantum Mechanics - Dipartimento di Fisica
Numerical Methods in Quantum Mechanics - Dipartimento di Fisica
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
nucleus. Even <strong>in</strong> open-shell atoms, this can be imposed as an approximation,<br />
by spherically averag<strong>in</strong>g ρ k . The simplification is considerable, s<strong>in</strong>ce we know<br />
a priori that the orbitals will be factorized as <strong>in</strong> Eq.(2.9). The angular part is<br />
given by spherical harmonics, labelled with quantum numbers l and m, while<br />
the ra<strong>di</strong>al part is characterized by quantum numbers n and l. Of course the<br />
accidental degeneracy for <strong>di</strong>fferent l is no longer present. It turns out that even<br />
<strong>in</strong> open-shell atoms, this is an excellent approximation.<br />
Let us consider the case of two-electron atoms. The Hartree equations,<br />
Eq.(6.14), for orbital k = 1 reduces to<br />
− ¯h2<br />
2m e<br />
∇ 2 1φ 1 (1) − Zq2 e<br />
r 1<br />
φ 1 (1) +<br />
[ ∫<br />
φ ∗ 2(2) q2 e<br />
r 12<br />
φ 2 (2) dv 2<br />
]<br />
φ 1 (1) = ɛ 1 φ 1 (1) (6.19)<br />
For the ground state of He, we can assume that φ 1 and φ 2 have the same spherically<br />
symmetric coord<strong>in</strong>ate part, φ(r), and opposite sp<strong>in</strong>s: φ 1 = φ(r)v + (σ),<br />
φ 2 = φ(r)v − (σ). Eq.(6.19) further simplifies to:<br />
− ¯h2<br />
2m e<br />
∇ 2 1φ(r 1 ) − Zq2 e<br />
r 1<br />
φ(r 1 ) +<br />
6.4 Code: helium hf ra<strong>di</strong>al<br />
[ ∫ q<br />
2<br />
e<br />
r 12<br />
|φ(r 2 )| 2 d 3 r 2<br />
]<br />
φ(r 1 ) = ɛφ(r 1 ) (6.20)<br />
Code helium hf ra<strong>di</strong>al.f90 1 (or helium hf ra<strong>di</strong>al.c 2 ) solves Hartree equations<br />
for the ground state of He atom. helium hf ra<strong>di</strong>al is based on code<br />
hydrogen ra<strong>di</strong>al and uses the same <strong>in</strong>tegration algorithm based on Numerov’s<br />
method. The new part is the implementation of the method of self-consistent<br />
field for f<strong>in</strong>d<strong>in</strong>g the orbitals.<br />
The calculation consists <strong>in</strong> solv<strong>in</strong>g the ra<strong>di</strong>al part of the Hartree equation<br />
(6.20). The effective potential V scf is the sum of the Coulomb potential of the<br />
nucleus, plus the (spherically symmetric) Hartree potential<br />
V scf (r) = − Zq2 e<br />
r<br />
+ V H (r), V H (r 1 ) = q 2 e<br />
∫ ρ(r2 )<br />
r 12<br />
d 3 r 2 . (6.21)<br />
We start from an <strong>in</strong>itial estimate for V H (r), calculated <strong>in</strong> rout<strong>in</strong>e <strong>in</strong>it pot (<br />
V (0)<br />
H<br />
(r) = 0, simply). With the ground state R(r) obta<strong>in</strong>ed from such potential,<br />
we calculate <strong>in</strong> rout<strong>in</strong>e rho of r the charge density ρ(r) = |R(r)| 2 /4π (note that<br />
ρ is here only the contribution of the other orbital, so half the total charge, and<br />
remember the presence of the angular part!). Rout<strong>in</strong>e v of rho re-calculates<br />
the new Hartree potential VH<br />
out (r) by <strong>in</strong>tegration, us<strong>in</strong>g the Gauss theorem:<br />
∫ r<br />
ṼH<br />
out (r) = V 0 + qe<br />
2 Q(s)<br />
r max s<br />
∫r