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.

6.2 BASIS SETS 171<br />

Table 6.1 STO-3G 2sp basis set for oxygen<br />

α2sp c2s c2p<br />

5.0331527 −0.099967 0.155916<br />

1.1695944 0.399513 0.607684<br />

0.3803892 0.700115 0.391957<br />

defined for each type <strong>of</strong> orbital core through valence. Thus for H and He, there is only a 1s<br />

function. For Li to Ne, there are five functions, 1s, 2s, 2px, 2py, and 2pz. For Na to Ar, 3s,<br />

3px, 3py, and 3pz are added to the second-row set, making a total <strong>of</strong> nine functions, etc. This<br />

number is the absolute minimum required, and it is certainly nowhere near the infinite basis<br />

set limit. Other minimal basis sets include the MINI sets <strong>of</strong> Huzinaga and co-workers, which<br />

are named MINI-1, MINI-2, etc., and vary in the number <strong>of</strong> primitives used for different<br />

kinds <strong>of</strong> functions.<br />

One way to increase the flexibility <strong>of</strong> a basis set is to ‘decontract’ it. That is, we might<br />

imagine taking the STO-3G basis set, and instead <strong>of</strong> constructing each basis function as a<br />

sum <strong>of</strong> three Gaussians, we could construct two basis functions for each AO, the first being<br />

a contraction <strong>of</strong> the first two primitive Gaussians, while the second would simply be the<br />

normalized third primitive. This prescription would not double the size <strong>of</strong> our basis set, since<br />

we would have all the same individual integrals to evaluate as previously, but the size <strong>of</strong> our<br />

secular equation would be increased. A basis set with two functions for each AO is called a<br />

‘double-ζ ’ basis. Of course, we could decontract further, and treat each primitive as a fullfledged<br />

basis function, in which case we would have a ‘triple-ζ ’ basis, and we could then<br />

decide to add more functions indefinitely creating higher and higher multiple-ζ basis sets.<br />

Modern examples <strong>of</strong> such basis sets are the cc-pCVDZ, cc-pCVTZ, etc. sets <strong>of</strong> Dunning<br />

and co-workers, where the acronym stands for ‘correlation-consistent polarized Core and<br />

Valence (Double/Triple/etc.) Zeta’ (Woon and Dunning 1995); correlation consistency and<br />

polarization are described in more detail below.<br />

The advantage <strong>of</strong> such a scheme is, naturally, that these increasingly large basis sets must<br />

come closer and closer to the HF limit. Let us step back for a moment, however, and consider<br />

the chemical consequences <strong>of</strong> providing extra basis functions for a given AO. Recall that a<br />

final MO from an HF calculation is a linear combination <strong>of</strong> all <strong>of</strong> the basis functions. Indeed,<br />

if we were to examine the 1s core orbital resulting from an HF calculation on atomic oxygen<br />

using the fully uncontracted set <strong>of</strong> STO-3G Gaussian primitives as a basis (i.e., a triple-ζ<br />

basis), we might well find it to be a linear combination <strong>of</strong> the 1s functions very similar to<br />

that defining the STO-3G contracted oxygen 1s function. And, if we were to look at the<br />

MOs resulting from an equivalent calculation on, say, formaldehyde (H2C=O), we would<br />

probably find another orbital which we would assign as the oxygen 1s orbital having very<br />

similar AO coefficients. Indeed, we would find this same orbital little changed in almost any<br />

molecule incorporating oxygen we might choose to examine. The reason for this is that core<br />

orbitals are only weakly affected by chemical bonding.<br />

Valence orbitals, on the other hand, can vary widely as a function <strong>of</strong> chemical bonding.<br />

Atoms bonded to significantly more electronegative elements take on partial positive charge

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

Saved successfully!

Ooh no, something went wrong!