OPEN(n1,n2)

RHF description: The M.O. occupancy during the SCF calculation can be defined in terms of doubly occupied, empty, and fractionally occupied M.O.s. By default, RHF SCF calculations are run using doubly occupied M.O.s, with, at most, one singly occupied M.O.  For some systems, the symmetry of the M.O.s can be preserved only if two or more M.O.s have fractional occupancies. Such fractional occupancies can be defined using OPEN(n1,n2), where n1 = number of electrons in the open-shell manifold, and n2 = number of open-shell M.O.s;  The ratio of n1 to n2 is restricted to  2 >  n1/n2  > 0.  Example: if OPEN(3,4)were used then the occupancy near the HOMO-LUMO gap would be ...2, 2, 2, 0.75, 0.75, 0.75, 0.75, 0, 0, 0, ...

The two π M.O.s in  molecular oxygen would not be degenerate if a closed-shell SCF were run.  In order to maintain the degeneracy, OPEN(2,2) would need to be used.

Do not use OPEN(n1,n2) for ground-state systems except for high symmetry systems with open shells, such as twisted (D2d) ethylene or molecular oxygen O2, or if there is a very small band-gap (such as in metal clusters).

Examples of OPEN(n1,n2) are given in the Table. OPEN(1,1) will be assumed for odd-electron systems unless an OPEN keyword is used. Errors introduced by use of fractional occupancy are automatically corrected [136] in a MECI calculation when OPEN(n1,n2) is used. See also C.I.=n.

Notes:

UHF description: The M.O. occupancy during the SCF calculation can be defined in terms of singly occupied, empty, and fractionally occupied M.O.s. By default, UHF SCF calculations are run using singly occupied M.O.s. When OPEN(n1,n2) is present, the highest n2 α molecular orbitals are each given a population of n1/n2 electrons.  To use OPEN(n1,n2) in a UHF calculation, UHF must also be present.  Because fractional electrons are used, small errors in energy, on the order of 1-2 kcal/mol, are introduced when OPEN(n1,n2) is used.  These errors are not corrected by C.I., however they are small and should not affect the results significantly.

Because OPEN only applies to the α electrons, care must be exercised in defining the number of α and β electrons.  For example, if the cation of tetrahedral methane were to be calculated, the keywords needed would be "OPEN(2,3)MS=-0.5 UHF CHARGE=1" If MS=-0.5 was not present, then there would be four α and three β electrons, instead of the required three α and four β electrons.

The main use of OPEN(n1,n2) in UHF calculations is to mimic the dynamic Jahn-Teller effect in transition metal complexes.  Without OPEN(n1,n2)an octahedral complex  system that could undergo Jahn-Teller distortion would be distorted.  However, in many X-ray structures, systems that should Jahn-Teller distort do not in fact distort and instead retain high symmetry, typically Oh.  This high symmetry can be regarded as a dynamic J-T effect. Transition metal complexes are, as their name suggests, complicated, and correctly defining the various J-T states is also complicated.  The following set of examples is provided to illustrate how to define keywords for various octahedral systems.  These systems are assumed to have a three-fold degenerate set of "d" M.O.s, with dominant contributions from the d-orbitals dxy, dyz, dxz, that can be related to the t2g set of point group Oh, and a two-fold degenerate set, composed mainly of the other two d-orbitals; these are equivalent to the  eg set of point group Oh.

d0: No keyword needed - closed shell, i.e., 1A1g.
d
1: keywords: "UHF OPEN(1,3)"  One α d electron in a  t2g set of M.O.s.  The resulting state is 2T2g.
d2: keywords: "UHF OPEN(2,3) MS=1"  Two α d electrons in a  t2g set of M.O.s.  The resulting state is 3T1g.
d3: keywords: "UHF MS=1.5"  The state is 4A2g.
d4: keywords for low spin: "UHF OPEN(1,3)MS=-1"  Because OPEN(n1,n2), applies only to α electrons, and because there are four electrons in a  t2g set of M.O.s, the only way to have a fractional population of the α set of M.O.s is to have the β set filled, i.e., α1β3 . By using "MS=-1" the number of  β electrons is defined as being two more than the number of α electrons  The resulting state is 3T1g.
d
4: keywords for high spin: "UHF OPEN(1,2) MS=2"  In high-spin d4 complexes, there are three α electrons in the t2g  set of M.O.s, and a single α electron in the two-fold degenerate eg M.O.  "MS=2" results in there being four more α than β electrons.  The resulting state is 5Eg.
d5: keywords for low spin: "UHF OPEN(2,3) MS=-0.5"  See d4 low spin. This gives rise to a  2T2g state.
d
5: keywords for high spin: "UHF MS=2.5"  The state is 6A1g.
d6: No keyword needed - closed shell, i.e., 1A1g.
d
6: keywords for high spin: "UHF OPEN(1,3) MS=-2"  The resulting state is 5T2g = A2g (from the β2 of eg symmetry)*A2g (from the β3 of t2g symmetry)*T2g (from the α1 of t2g symmetry).
d7: keywords for low spin: "UHF OPEN(1,2)"  As with d1, the extra electron has, by default, α spin, and so the keyword "MS=0.5" is unnecessary.  The resulting state is 2Eg.
d
7: keywords for high spin: "UHF OPEN(1,2) MS=-1.5"  With seven d-electrons and  MS=-1.5 there are three more α than β electrons, i.e., α2β5.  The resulting state is  4T1g.
d8: keywords: "UHF MS=1"   The resulting state is  3A2g.
d9: keywords: "UHF OPEN(1,2)"   The resulting state is  2Eg.
d10: No keyword needed - closed shell, i.e., 1A1g.

A similar exercise could be done with the Td complexes. As this is simple, albeit tedious, it is left as an exercise.

Table: Use of OPEN(n,m)
UHF Keywords

RHF Keywords

Example

Number of M.O.s

No. of Electrons

State
MS=1

OPEN(2,2)

Twisted Ethylene

2

2

3A2
OPEN(1,2)

OPEN(1,2)

O2+

2

1

2πg

MS=-0.5 OPEN(2,3)

OPEN(5,3)

CH4+

3

5

2T2
MS=1.5 OPEN(3,3) [CrIIIF6]3- 3 3 4A2g
MS=2.5 OPEN(5,5) [MnII(H2O)6]2+ 5 5 6A1
MS=-0.5 OPEN(2,3) OPEN(5,3) [FeIII(CN)6]3- 3 5 2T2g