07.04.2013 Views

Essentials of Computational Chemistry

Essentials of Computational Chemistry

Essentials of Computational Chemistry

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

or<br />

7.4 PERTURBATION THEORY 219<br />

cj = 〈(0) j |V| (0)<br />

0 〉<br />

a (0)<br />

0<br />

− a(0)<br />

j<br />

(7.39)<br />

With the first-order eigenvalue and wave function corrections in hand, we can carry out<br />

analogous operations to determine the second-order corrections, then the third-order, etc.<br />

The algebra is tedious, and we simply note the results for the eigenvalue corrections, namely<br />

and<br />

a (3)<br />

0<br />

= <br />

j>0,k>0<br />

a (2)<br />

0<br />

〈 (0)<br />

0 |V|(0) j 〉[〈 (0)<br />

= <br />

j>0<br />

|〈 (0)<br />

j |V| (0)<br />

0 〉|2<br />

a (0)<br />

0<br />

− a(0)<br />

j<br />

j |V| (0)<br />

(0)<br />

k 〉−δjk〈 0 |V|(0) 0 〉]〈(0) k |V| (0)<br />

0 〉<br />

(a (0)<br />

0 − a(0)<br />

j )(a (0)<br />

0 − a(0)<br />

k )<br />

(7.40)<br />

(7.41)<br />

Let us now examine the application <strong>of</strong> perturbation theory to the particular case <strong>of</strong> the<br />

Hamiltonian operator and the energy.<br />

7.4.2 Single-reference<br />

We now consider the use <strong>of</strong> perturbation theory for the case where the complete operator<br />

A is the Hamiltonian, H. Møller and Plesset (1934) proposed choices for A (0) and V with<br />

this goal in mind, and the application <strong>of</strong> their prescription is now typically referred to by<br />

the acronym MPn where n is the order at which the perturbation theory is truncated, e.g.,<br />

MP2, MP3, etc. Some workers in the field prefer the acronym MBPTn, to emphasize the<br />

more general nature <strong>of</strong> many-body perturbation theory (Bartlett 1981).<br />

The MP approach takes H (0) to be the sum <strong>of</strong> the one-electron Fock operators, i.e., the<br />

non-interacting Hamiltonian (see Section 4.5.2)<br />

H (0) =<br />

n<br />

fi<br />

i=1<br />

(7.42)<br />

where n is the number <strong>of</strong> basis functions and fi is defined in the usual way according to<br />

Eq. (4.52). In addition, (0) is taken to be the HF wave function, which is a Slater determinant<br />

formed from the occupied orbitals. By analogy to Eq. (4.36), it is straightforward to show<br />

that the eigenvalue <strong>of</strong> H (0) when applied to the HF wave function is the sum <strong>of</strong> the occupied<br />

orbital energies, i.e.,<br />

H (0) (0) occ.<br />

= εi (0)<br />

(7.43)<br />

where the orbital energies are the usual eigenvalues <strong>of</strong> the specific one-electron Fock operators.<br />

The sum on the r.h.s. thus defines the eigenvalue a (0) .<br />

i

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

Saved successfully!

Ooh no, something went wrong!