07.04.2013 Views

Essentials of Computational Chemistry

Essentials of Computational Chemistry

Essentials of Computational Chemistry

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

252 8 DENSITY FUNCTIONAL THEORY<br />

for the former. In particular, he suggested that the exchange hole about any position could<br />

be approximated as a sphere <strong>of</strong> constant potential with a radius depending on the magnitude<br />

<strong>of</strong> the density at that position. Within this approximation, the exchange energy Ex is<br />

determined as<br />

Ex[ρ(r)] =− 9α<br />

8<br />

1/3 <br />

3<br />

π<br />

ρ 4/3 (r)dr (8.7)<br />

Within Slater’s derivation, the value for the constant α is 1, and Eq. (8.7) defines so-called<br />

‘Slater exchange’.<br />

Starting from the uniform electron gas, Bloch and Dirac had derived a similar expression<br />

several years previously, except that in that case α = 2<br />

3 (Bloch, F. 1929 and Dirac,<br />

P. A. M. 1930). The combination <strong>of</strong> this expression with Eqs. (8.3)–(8.5) defines the<br />

Thomas–Fermi–Dirac model, although it too remains sufficiently inaccurate that it fails<br />

to see any modern use.<br />

Given the differing values <strong>of</strong> α in Eq. (8.7) as a function <strong>of</strong> different derivations, many<br />

early workers saw fit to treat it as an empirical value, and computations employing Eq. (8.7)<br />

along these lines are termed Xα calculations (or sometimes Hartree–Fock–Slater calculations<br />

in the older literature). Empirical analysis in a variety <strong>of</strong> different systems suggests that α = 3<br />

4<br />

provides more accurate results than either α = 1orα = 2<br />

3 . This particular DFT methodology<br />

has largely fallen out <strong>of</strong> favor in the face <strong>of</strong> more modern functionals, but still sees occasional<br />

use, particularly within the inorganic community.<br />

8.2 Rigorous Foundation<br />

The work described in the previous section was provocative in its simplicity compared to<br />

wave-function-based approaches. As a result, early DFT models found widespread use in<br />

the solid-state physics community (where the enormous system size required to mimic the<br />

properties <strong>of</strong> a solid puts a premium on simplicity). However, fairly large errors in molecular<br />

calculations, and the failure <strong>of</strong> the theories to be rigorously founded (no variational principle<br />

had been established), led to their having little impact on chemistry. This state <strong>of</strong> affairs was<br />

set to change when Hohenberg and Kohn (1964) proved two theorems critical to establishing<br />

DFT as a legitimate quantum chemical methodology. Each <strong>of</strong> the two theorems will be<br />

presented here in somewhat abbreviated form.<br />

8.2.1 The Hohenberg–Kohn Existence Theorem<br />

In the language <strong>of</strong> DFT, electrons interact with one another and with an ‘external potential’.<br />

Thus, in the uniform electron gas, the external potential is the uniformly distributed positive<br />

charge, and in a molecule, the external potential is the attraction to the nuclei given by<br />

the usual expression. As noted previously, to establish a dependence <strong>of</strong> the energy on the<br />

density, and in the Hohenberg–Kohn theorem it is the ground-state density that is employed,<br />

it is sufficient to show that this density determines the Hamiltonian operator. Also as noted<br />

previously, integration <strong>of</strong> the density gives the number <strong>of</strong> electrons, so all that remains to

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!