Definition of some C.I. Keywords

It has been my policy, ever since the first release of MOPAC in 1983, to resist changing the definition of keywords. This policy has allowed users to confidently use a new MOPAC in the belief that old keywords will have their old, familiar, meaning. However, an ambiguity was found in certain keywords, an ambiguity which has, at times, resulted in severe frustration.

Consider the word TRIPLET. In MOPAC 6, this meant (but no longer means) "Do an SCF calculation in which the M.O. populations are [...,2,2,2,1,1,0,0,...], then do a C.I. on the two half-occupied M.O.s, and select the triplet state." This definition meant that twisted ethylene would have the correct symmetry, as would triplet oxygen. However, if a user wanted to examine triplet formaldehyde, and compare it with the ground state, problems arose. The keywords C.I.=2 ROOT=2 would generate the correct energy, but a user might expect that TRIPLET should achieve the same result. Because of the definition of TRIPLET, the SCF starting configuration was different, and as a result, the ΔHf was also different. Under earlier MOPACs, there was no way to set up a calculation using the keyword TRIPLET and go SCF on a closed-shell configuration as the precursor to a C.I. calculation.

Because of the limitations of the earlier definitions of spin-states (TRIPLET, QUARTET, QUINTET, SEXTET, etc.), these words were all redefined in 1993, in MOPAC 93. In order to reproduce the earlier keywords, pairs of keywords, such as TRIPLET OPEN(2,2) or SEXTET OPEN(5,5) must now be used. Spin-states which result from SCF calculations on ground-state configurations can now be specified by the following pairs of keywords: TRIPLET C.I.=2; QUARTET C.I.=3; QUINTET C.I.=4; SEXTET C.I.=5.

Using these new definitions, spin-states of a system can now be more easily related. Consider the various states of formaldehyde (Table), in which all calculations use the ground-state geometry and 1SCF.

 

Table :

Examples of Use of C.I. Keywords

Keywords Used

$Delta H_f$

(No keywords)

-32.9040

C.I.=1

-32.9040

C.I.=2

-33.0166

C.I.=3

-39.7234

C.I.=4

-39.9665

C.I.=5

-40.1743

C.I.=2 TRIPLET

29.6348

C.I.=3 ROOT=2

28.2840

C.I.=3 TRIPLET

28.2840

C.I.=3 TRIPLET MS=0

28.2840

OPEN(2,2) TRIPLET

27.9318

Now we see that C.I.=3 ROOT=2 and C.I.=3 TRIPLET do, in fact, give the same result. The "old" MOPAC (pre-1993) result of using TRIPLET can still be generated by OPEN(2,2) TRIPLET. Note that C.I.=1 generates the normal ΔHf of CH2O, and that increasing the C.I. lowers the energy steadily.