Page 10 - 71 the abc of dft_opt
P. 10

1.3. REINTERPRETING MOLECULAR ORBITALS                              19     20                                                CHAPTER 1. INTRODUCTION
             4
                                                                                                   !
             2     density
             0
             -2    Kohn-Sham potential
                                                                                                           1s
                                                                                           1s
                                                                                          1s
             -4             -2/r
                                                                                                   !
             -6
                             He atom
             -8
              0   0.2  0.4  r  0.6  0.8  1
        Figure 1.2: The external and Kohn-Sham potentials for the He atom, thanks to Cyrus’ Umrigar’s very accurate density.  Figure 1.3: Orbital diagram of H 2 bond formation.

       electrons. If we can figure out some way to approximate this potential accurately, we have  And a highly accurate approximate density functional calculation produces the full electronic
       a much less demanding set of equations to solve than those of the true system. Thus we  energy from these orbitals, resolving the paradox.
       are always trying to improve a non-interacting calculation of a non-interacting wavefunction,
       rather than that of the full physical system. In Fig. 1.1 there are also plotted the local
                                                                                  1.4  Sampler of applications
       density approximation (LDA) and generalized gradient approximation (GGA) curves. LDA
       is the simplest possible density functional approximation, and it already greatly improves on
       HF, although it typically overbinds by about 1/20 of a Hartree (or 1 eV or 30 kcal/mol),  In this section, we highlight a few recent applications of modern DFT, to give the reader a
                                                                                  feeling for what kinds of systems can be tackled.
       which is too inaccurate for most quantum chemical purposes, but sufficiently reliable for many
       solid-state calculations. More sophisticated GGA’s (and hybrids) reduce the typical error in  Example from biochemistry, using ONIOM.
                                                                                    Example from solid-state, eg ferromagnetism.
       LDA by about a factor of 5 (or more), making DFT a very useful tool in quantum chemistry.
                                                                                    Example using CP molecular dynamics, eg phase-diagram of C.
       In section 1.5, we show how density functionals work with a simple example from elementary
       quantum mechanics.                                                           Example from photochemistry, using TDDFT.

       1.3  Reinterpreting molecular orbitals                                     1.5  Particle in a box

       The process of bonding between molecules is shown in introductory chemistry textbooks as  In this section, we take the simplest example from elementary quantum mechanics, and apply
       linear combinations of atomic orbitals forming molecular orbitals of lower energy, as in Fig.  density functional theory to it. This provides the underlying concepts behind what is going
       1.3.                                                                       on in the more previous sections of this chapter. Much of the rest of the book is spent
         But later, in studying computational chemistry, we discover this is only the Hartree-Fock  connecting these two.
       picture, which, as stated above, is rarely accurate enough for quantum chemical calculations.  In general, we write our Hamiltonian as
       In this picture, we need a more accurate wavefunction, but then lose this simple picture of
                                                                                                                 ˆ
                                                                                                                     ˆ
                                                                                                                H = T + V ˆ                         (1.7)
       chemical bonding. This is a paradox, as chemical reactivity is usually thought of in terms of
       frontier orbitals.                                                               ˆ                                     ˆ
                                                                                  where T denotes the operator for the kinetic energy and V the potential energy. We begin
         In the Kohn-Sham approach, the orbitals are exact and unique, i.e., there exists (at most)  with the simplest possible case. The Hamiltonian for a 1-D 1-electron system can be written
       one external potential that, when doubly occupied by two non-interacting electrons, yields  as
       the exact density of the H 2 molecule. So in this view, molecular orbital pictures retain their        ˆ    1 d 2
       significance, if they are the exact Kohn-Sham orbitals, rather than those of Hartree-Fock.             H = − 2 dx 2  + V (x)                  (1.8)
   5   6   7   8   9   10   11   12   13   14   15