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.

9.1 PROPERTIES RELATED TO CHARGE DISTRIBUTION 313<br />

Mulliken partial atomic charge is then defined as<br />

qk = Zk − Nk<br />

(9.15)<br />

where Z is the nuclear charge and Nk is computed according to Eq. (9.14).<br />

With minimal or small split-valence basis sets, Mulliken charges tend to be reasonably<br />

intuitive, certainly in sign if not necessarily in magnitude. Analysis <strong>of</strong> changes in charge as<br />

a function <strong>of</strong> substitution or geometric change tends to be the best use to which Mulliken<br />

charges may be put, and this can <strong>of</strong>ten provide chemically meaningful insight, as illustrated<br />

in Figure 9.3.<br />

The use <strong>of</strong> a non-orthogonal basis set in the Mulliken analysis, however, can lead to<br />

some undesirable results. For instance, if one divides up the total number <strong>of</strong> electrons over<br />

AO basis functions (in a fashion exactly analogous to that used for atoms), it is possible<br />

for individual basis functions to have occupation numbers greater than 1 (which would be<br />

greater than 2 in a restricted theory) or less than 0, and such a situation obviously can<br />

have no physical meaning. In addition, the rule that all shared electrons should be divided<br />

up equally between the atoms on which the sharing basis functions reside would seem to<br />

ignore the possibly very different electronegativities <strong>of</strong> these atoms. Finally, Mulliken partial<br />

charges prove to be very sensitive to basis-set size, so that comparisons <strong>of</strong> partial charges<br />

Mulliken charge (a.u.)<br />

0.6<br />

0.5<br />

0.4<br />

0.3<br />

0.2<br />

0.1<br />

0<br />

−0.1<br />

H<br />

O<br />

H<br />

H<br />

H<br />

O<br />

H<br />

C O<br />

H<br />

H<br />

C O<br />

H<br />

H<br />

0 20 40 60 80 100<br />

Dihedral angle (°)<br />

H<br />

H<br />

H<br />

O<br />

120 140 160 180<br />

H<br />

C O<br />

Figure 9.3 AM1 Mulliken charges <strong>of</strong> the hydroxyl (circles) and oxonium (squares) oxygen atoms in<br />

protonated dihydroxymethane as a function <strong>of</strong> HOCO + dihedral angle. Standard precepts <strong>of</strong> conformational<br />

analysis suggest that hyperconjugation <strong>of</strong> hydroxyl oxygen lone-pair density (acting as a<br />

donor) into the C–O + σ ∗ orbital (acting as an acceptor) may occur, and the effect is expected to be<br />

maximal at a dihedral angle <strong>of</strong> 90 ◦ , and minimal at 0 ◦ and 180 ◦ . The computed Mulliken charges<br />

on the oxygen atoms support hyperconjugation being operative, with about one-tenth <strong>of</strong> a positive<br />

charge being transferred from the oxonium oxygen to the hydroxyl oxygen at a dihedral angle <strong>of</strong> 90 ◦<br />

compared to 180 ◦ (an interpretation also consistent with geometric and energetic analysis, see Cramer<br />

1992)<br />

H<br />

H

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

Saved successfully!

Ooh no, something went wrong!