Page 53 - Valence Bond Methods. Theory and Applications
P. 53

36
                                                      2H 2 and localized orbitalØ
                                                     2.5.2 Kinetic energy effectp
                             Wheð we gc beyond the simple ccvalent treatment and include an ionic function ið
                             the Heitler–London treatment, we obtaið a further lowering of the calculateł energy.
                             At first glance, perhaps this is surprising, since the ionic function has more electron
                             repulsion than the ccvalent. Although we sŁw ið Section 1.3.à that any additional
                             linear variation term must lower the energy, that does not give any physical picture
                             of the process. We will now give a detaileł analysis of hcw the lowering comes
                             about and its physical origin.
                               Ið the previous section we examineł the variational result of the two-term wave
                             function consisting of the ccvalent and ionic functions. This produces a 2 × 2
                             Hamiltonian, which may be decomposeł intc kinetic energy, nuclear attraction,
                             and electron repulsion terms. Each of these operators produces a 2 × 2 matrix.
                             Along with the overlap matrix these are


                                                   1     S IC            T II  T IC
                                              S =             ;    T =              ;
                                                   S CI   1              T CI  T CC

                                                   V II  V IC            G II  G IC
                                             V n =            ;   G e =              .
                                                   V CI  V CC            G CI  G CC
                                                                                            6
                             As we discusseł abcve, the two functions hŁve the same charge density, and this
                             implies that T II = T CC and V II = V CC , but we expect G II > G CC .
                               The normalization of the wave function requires
                                                                            2
                                                          2
                                                     1 = C + 2S IC C I C C + C .                (2.44)
                                                                            C
                                                          I
                             Two similar expressions give us the expectation values of T amł V n ,
                                                        2
                                                                            2
                                                
T  Ø C T II + 2T IC C I C C + C T CC ,         (2.45)
                                                        I                   C
                                                        2
                                                                             2
                                                
V n  Ø C V II + 2V IC C I C C + C V CC .       (2.46)
                                                        I                    C
                             Multiplying Eq. (2.44) by T CC and V CC ið turð and subtracting the result from the
                             corresponding Eq. (2.45) or Eq. (2.46), we arrive at the equations
                                                
T  à T CC = 2(T IC − T CC S IC )C I C C ,      (2.47)
                                                
V n  à V CC = 2(V IC − V CC S IC )C I C C ,    (2.48)
                             and we see that the differences depend on hcw the off-diagonal matrix elements
                             compare tc the overlap times the diagonal elements. A similar expression for G e is
                             more complicated:
                                                                                          2
                                      
G e  à G CC = 2(G IC − G CC S IC )C I C C + (G II − G CC )C .  (2.49)
                                                                                          I
                             6  They hŁve the same first order density matrices.
   48   49   50   51   52   53   54   55   56   57   58