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.

14.2 SINGLY EXCITED STATES 495<br />

Table 14.1 Energies (kcal mol −1 ), where available, for lowest<br />

singlet excited states <strong>of</strong> phenylnitrene relative to the 3 A2 ground<br />

state a<br />

Source 1 A2 1 1 A1 2 1 A1<br />

HF/6-31G(d) 64.6 80.1<br />

BLYP/6-31G(d) 14.5 b , 22.8 c 31.4 48.7<br />

BPW91/cc-pVDZ (SCF) d 14.3 b 33.9 43.0<br />

BLYP/cc-pVTZ (SCF) e 29.5 41.0<br />

CCSD(T)/cc-pVDZ e 35.2 47.2<br />

CAS(8,8)/cc-pVDZ d 17.8 42.1 76.2<br />

CASPT2/cc-pVDZ d 19.3 37.4 57.8<br />

CASPT2/cc-pVTZ e 19.3 34.8 54.5<br />

MRCISD/DZP f 21.0 39.8 (52.0)<br />

Experiment g 18. 30. n.a.<br />

a Zero-point vibrational energies are not included in the theoretical energies,<br />

but ZPVE differences between alternative electronic states are predicted to<br />

be small at the few levels where they have been evaluated.<br />

b Determined from sum method.<br />

c Determined from spin projection.<br />

d Johnson and Cramer (2001).<br />

e Smith and Cramer (1996).<br />

f Hrovat, Waali, and Borden (1992); Kim, Hamilton, and Schaefer (1992).<br />

g Travers et al. (1992); Ellison, G. B., unpublished results.<br />

levels <strong>of</strong> theory that would be expected to be reasonably accurate suggests that the DFT<br />

predictions are substantially too low. The DFT wave function for the 2 1 A1 state is that <strong>of</strong><br />

Eq. (14.8) and, as discussed above, is found by fortuitous convergence <strong>of</strong> the SCF equations<br />

for this occupation scheme where variational collapse to the 1 1 A1 state <strong>of</strong> Eq. (14.7) would<br />

otherwise be expected. The apparently rather poor accuracy <strong>of</strong> the energy for the higher state<br />

suggests that this orthogonality issue cannot be ignored here, and the SCF procedure must<br />

be regarded as unreliable.<br />

As for the HF level, the SCF approach for the closed-shell singlet states is identical to<br />

that in the DFT case (in this instance, the two-determinantal nature <strong>of</strong> the lower energy openshell<br />

singlet requires an MCSCF description, so HF values are not reported for this state).<br />

However, both <strong>of</strong> the closed-shell singlets are subject to large non-dynamical correlation<br />

effects (as a consequence, in part, <strong>of</strong> being so close in energy to one another). Since HF<br />

theory is much more sensitive to such correlation than DFT, the energies <strong>of</strong> these two states<br />

are predicted to be much too high. This error is in some sense even worse than it appears,<br />

because severe spin contamination in the triplet, which exhibits an expectation value for S 2<br />

in excess <strong>of</strong> 2.7, probably causes it too to be poorly represented at the HF level.<br />

Of course, with HF wave functions in hand, it is possible to carry out post-HF calculations<br />

to partially correct for electron correlation effects. The poor quality <strong>of</strong> the HF wave functions,<br />

however, militate against any treatment much less sophisticated than coupled-cluster. At the<br />

CCSD(T)/cc-pVDZ level, the predicted energy <strong>of</strong> the lowest closed-shell singlet is in fair<br />

agreement with experiment (other data in the table suggest that use <strong>of</strong> a triple-ζ basis set<br />

would improve the CCSD(T) estimate). The energy <strong>of</strong> the second closed-shell singlet state

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

Saved successfully!

Ooh no, something went wrong!