Choice of State to be Optimized

MECI can calculate a large number of states of various total spin. Two schemes are provided to allow a given state to be selected. First, ROOT=n will, when used on its own, select the n'th state, irrespective of its total spin. By default, n=1. If ROOT=n is used in conjunction with a keyword from the set SINGLET, DOUBLET, TRIPLET, QUARTET, QUINTET, SEXTET, SEPTET, OCTET, or NONET. , then the n'th root of that spin-state will be used. For example, ROOT=4 and SINGLET will select the 4th singlet state. If there are two triplet states below the fourth singlet state then this will mean that the sixth state will be selected.

Sometimes the energy required to form an excited state is wanted. By this we mean the energy of the excited state relative to the energy of the ground state, and not the heat of formation of the excited state. To calculate this quantity, the keywords PRECISE, GNORM=0.01, MECI and C.I.=2 should be used. For formaldehyde, these keywords would produce the output shown in the Figure.

 Figure: Energies of Excited States

State

Energy (eV)

Q.N.

Spin

Symmetry

 

Absolute

Relative

 

 

 

1

-0.0049

0.0000

1

Singlet

A1

2

2.7109

2.7158

1

Triplet

A2

3

3.1029

3.1078

2

Singlet

A2

4

7.8630

7.8679

2

Singlet

A1


This output can be read as follows: The first state (the one at -0.004891eV) is the new ground state. C.I. will lower the energy of the ground state, relative to the SCF ground state, and for formaldehyde this extra stabilization amounts to 0.0049 eV. The ground state is a singlet, and has A1 symmetry. The second state is a triplet, with energy 2.7109eV above the SCF energy, or 2.7158eV above the ground state, and has A2 symmetry. The third and fourth states are both singlets.

Using the two keywords given, the system would optimize on the ground singlet state, and the bond orders and density matrix would reflect this. If the first excited singlet state were wanted, then the extra keywords ROOT=2 and SINGLET would also be used. Alternatively, the single extra keyword ROOT=3 could be used. If the first triplet state were wanted, then TRIPLET or ROOT=2 (but not both!) could be used.