- Search
- Download PDF

(June 30, 2021)

The maximum overlap method (Sec. 4.5.11) is successful in many cases, but when optimizing excited state orbitals it can be prone to variational collapse (falling from the target configuration to the ground state). This behavior was improved with the “initial” MOM (IMOM) procedure, which uses the initial guess molecular orbital coefficients as the reference for the overlap criterion (see Sec. 7.6 for details). While IMOM lends a significant improvement to the MOM, it occasionally also falls victim to variational collapse, incentivizing the pursuit of alternatives to these methods.

One such alternative, proposed by Kevin Carter-Fenk and John Herbert, is the
state-targeted energy projection (STEP) procedure, which applies a simple level-shift
formalism to converge to the target state.
^{
158
}
J. Chem. Theory Comput.

(2020),
16,
pp. 5067.
Link
The STEP algorithm retains the cost-effectiveness of the MOM procedures
(about the same cost per cycle as a normal SCF), while simultaneously being
far more robust in converging to the target state.
Further details on the STEP procedure can be found in Sec. 7.8.3, and
for a complete account of the STEP algorithm the reader is referred to Ref. 158.

The STEP-related *$rem* variables are the following:

STEP

Activates the STEP procedure.

TYPE:

LOGICAL

DEFAULT:

FALSE

OPTIONS:

FALSE
Do not apply the STEP level-shift algorithm.
TRUE
Apply the STEP level-shift algorithm.

RECOMMENDATION:

None

STEP_EPSILON

Scales the size of the occupied/virtual gap imposed by the level-shift by $N$/100 Hartree.

TYPE:

INTEGER

DEFAULT:

10

OPTIONS:

$N$

RECOMMENDATION:

Use the default unless convergence issues arise, in which case a larger value
can be used until the desired state is found. Be aware that increasing the occupied/virtual
gap in level-shift algorithms slows convergence so it may be advisable to
increase SCF_MAX_CYCLES if large shifts are required.

STEP_PRINT

Controls the print level for STEP algorithm information.

TYPE:

INTEGER

DEFAULT:

1

OPTIONS:

0
Do not print any information about STEP between SCF cycles.
1
Print the level-shift applied at each SCF cycle (R- and U-STEP).
2
Print the level-shift for both mixed and triplet states at each SCF cycle (RO-STEP).

RECOMMENDATION:

Use the default. Level shifts of 0 indicate that an aufbau criterion
is sufficient to determine orbital occupation, and shifts $>0$ imply
non-aufbau selection of the occupied space.