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.

7.6 PRACTICAL ISSUES IN APPLICATION 231<br />

this term is infinite! Note that this is not a case <strong>of</strong> the actual electronic state being degenerate,<br />

but is entirely an artifact <strong>of</strong> using a single-reference wave function.<br />

As a general rule, whenever the frontier orbital separation becomes small, the magnitude<br />

<strong>of</strong> the MPn energy terms will become large because <strong>of</strong> their inverse dependence on orbital<br />

energy separation, and perturbation theory will be very slowly convergent in such instances.<br />

Such decreased separations also increase the degree <strong>of</strong> multireference character in the wave<br />

function, so the tight coupling between these two phenomena is rationalized. Figure 7.5<br />

provides an example <strong>of</strong> this phenomenon for the case <strong>of</strong> the energy <strong>of</strong> carbonyl oxide,<br />

which has a moderate degree <strong>of</strong> multireference character, relative to its isomer dioxirane.<br />

The comparatively small size <strong>of</strong> these systems permits the extension <strong>of</strong> perturbation theory<br />

through fifth order with the 6-31G(d) basis set, but even the difference between the MP4<br />

and MP5 results remains a fairly large 1.3 kcal mol −1 .<br />

CCSD is similarly sensitive to multireference character, although it is less obvious that this<br />

should be so based on the formalism presented above. However, inclusion <strong>of</strong> triples in the<br />

CCSD wave function is usually very effective in correcting for a single-reference treatment<br />

<strong>of</strong> a weakly to moderately multireference problem. Of course, the most common way to<br />

include the triples is by perturbation theory, i.e., CCSD(T), and as noted above, this level<br />

too can be unstable if singles amplitudes are large. In such an instance, BD(T) calculations,<br />

which eliminate the singles amplitudes, can be efficacious.<br />

To illustrate some <strong>of</strong> these points in greater detail, consider the aryne diradicals p-benzyne<br />

and its isoelectronic but charged congener, N-protonated 2,5-pyridyne (Figure 7.6). In these<br />

systems, the frontier orbitals <strong>of</strong> interest are the bonding and antibonding combinations <strong>of</strong> the<br />

‘singly occupied’ σ orbitals left after abstraction <strong>of</strong> two hydrogen atoms from the aromatic<br />

ring. Because the orbitals are para-related to one another, the interaction between them<br />

is weak, and thus the frontier orbital energy separation is small. As such, the closed-shell<br />

singlet and triplet states lie relatively near one another; in the case <strong>of</strong> p-benzyne, negative ion<br />

Energy difference (kcal mol −1 )<br />

−28<br />

−30<br />

−32<br />

−34<br />

−36<br />

−38<br />

−40<br />

−42<br />

0<br />

1 2 3 4 5<br />

Order <strong>of</strong> perturbation theory<br />

6<br />

H<br />

O<br />

O<br />

H H<br />

O O<br />

Figure 7.5 Slowly oscillatory behavior <strong>of</strong> MPn/6-31G(d)//HF/6-31G(d) theory for the energy separation<br />

between carbonyl oxide and dioxirane. Accurate extrapolation from this perturbation series is<br />

an unlikely prospect<br />

?<br />

H

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

Saved successfully!

Ooh no, something went wrong!