For molecules which are believed to have biradicaloid character the option exists to optimize the lowest singlet energy state which results from the mixing of three states. These states are, in order, (1) the (micro)state arising from a one electron excitation from the HOMO to the LUMO, which is combined with the microstate resulting from the time-reversal operator acting on the parent microstate, the result being a full singlet state; (2) the state resulting from de-excitation from the formal LUMO to the HOMO; and (3) the state resulting from the single electron in the formal HOMO being excited into the LUMO.
Microstate 1 Microstate 2 Microstate 3
Alpha Beta Alpha Beta Alpha Beta Alpha Beta
LUMO * * * *
--- --- --- --- --- --- --- ---
+
HOMO * * * *
--- --- --- --- --- --- --- ---
A configuration interaction calculation is involved here. A biradical
calculation done without C.I. at the RHF level would be meaningless.
Either rotational invariance would be lost, as in the D2d form of
ethylene, or very artificial barriers to rotations would be found, such
as in a methane molecule ``orbiting'' a D2d ethylene. In both cases the
inclusion of limited configuration interaction corrects the error.
BIRADICAL should not be used if either the HOMO or LUMO is degenerate; in
this case, the full manifold of HOMO
LUMO should be included in the
C.I., using MECI options. The user should be aware of this situation.
When the biradical calculation is performed correctly, the result is
normally a net stabilization. However, if the first singlet excited
state is much higher in energy than the closed-shell ground state,
BIRADICAL can lead to a destabilization. Abbreviation: BIRAD. See also
MECI, C.I., OPEN, SINGLET.