Oxidation State Localized Orbitals (OSLO) is a new localization scheme focused
on molecular fragments for the purposes of oxidation state assignment.
393
J. Chem. Theory Comput.
(2022),
18,
pp. 309.
Link
The method has been developed to avoid some pitfalls
encountered in the LOBA method where the PM orbitals can spread across many
fragments without reaching the 60% threshold. OSLO starts by looking for
fragments’ centers of charges, and localizes on radial spread (in the real
space) to those centers. Then, it admits orbitals to its list of orbitals if
they are above a slowly increasing threshold of a criterion called the
fragment orbital localization index (FOLI). To understand this criterion, one needs to
know Pipek’s delocalization measure. This is defined as
(10.9) |
where the is the population of the th orbital on center and the
summation runs over all centers. When an orbital is entirely localized on a
given center, its delocalization measure is . If the orbital is perfectly
delocalized among two centers and then and the
delocalization measure is , etc. The square helps make it less
sensitive to the ratio of the population to each other compared to how many are
related. In fact, minimizing the sum over all occupied orbitals is precisely
what leads to Pipek-Mezey localization procedure.
954
J. Chem. Phys.
(1989),
90,
pp. 4916.
Link
However, for our purposes,
out of the localized orbitals generated from fragment with low delocalization
measure, we are interested in those that are also highly localized on fragment .
Defining the fragment orbital localization index (FOLI),
(10.10) |
it is easy to see tends to unity for fragment when orbital is perfectly localized on (), tends to 2 when the localized orbital is perfectly delocalized over two fragments (), and gradually increases as the orbital becomes more delocalized as well as less centered on fragment . Thus, among the redundant set of localized orbitals, one selects all orbitals above the smallest FOLI by a threshold (0.01 by default) and assigns each orbital to the originator fragment. The orbitals which are then symmetrically orthogonalized, and projected out from the space of remaining unassigned occupied space. In case a set of orbitals is redundant due to a symmetry in the system or simple covalency, then the orbitals are split over all contributing fragments (originators). After the iterative process, each molecular fragment has associated a set of localized orbitals derived from the simplest orbital spread criterion, which in turns determines the fragment’s formal charge or oxidation state in a natural manner.
The new method expects fragments, otherwise it localizes on atomic centers instead. Although this method was developed for oxidation state, it produces a set of localized orbitals on fragments or atoms that can be used like any other localization method.
These orbitals do not suffer from the multiple minima problem like most localization procedure that do iterative rotations discussed in the ER section. This is mainly because the iterative process essentially takes the low-hanging fruit orbitals first before going to higher up ones that are more ambiguous. In very few examples, it was found that the orbitals above the threshold could demonstrate another solution, i.e., another Lewis picture. Therefore, the algorithm was augmented to look into the next set of orbitals, print them out, and see if they are similar enough to the current ones, by looking at the singular values of the overlap of the current set and the next set.
OSLO
OSLO
Triggers OSLO procedure after a converged SCF
TYPE:
INTEGER
DEFAULT:
0
OPTIONS:
0
Don’t perform OSLO
1
Perform the OSLO procedure
RECOMMENDATION:
None
More precise control of OSLO goes under the $loco input section, with keywords that are introduced below. See the example for reference.
OSLO_PopAnalType
Specifies which population analysis is used for OSLOs.
INPUT SECTION: $loco
TYPE:
INTEGER
DEFAULT:
2
OPTIONS:
0
Mulliken
1
Lowdin
2
IAOs
RECOMMENDATION:
Use IAOs since they are insenstive to the basis set. When using IAOs, use
AUTOSAD = TRUE in the $loco input section to use the new AUTOSAD-IAO algorithm, which
uses AUTOSAD atomic densities instead of some minimal basis.
OSLO_THRESH
Specifies the threshold for FOLI when adding new orbitals for each iteration.
INPUT SECTION: $loco
TYPE:
INTEGER
DEFAULT:
0
OPTIONS:
0
Use the default, 0.01
Use threshold of /1000
RECOMMENDATION:
Use the default
OSLO_BRANCHING
Specifies which iteration with too much overlap between its orbitals and next
sets orbitals to choose the next set. The input here can be , which is a
binary string for which iteration (only consider those with warning prints) to
branch in, e.g., 010011, means to branch in iterations 2, 5, and 6.
INPUT SECTION: $loco
TYPE:
INTEGER
DEFAULT:
0
OPTIONS:
, for each digit:
0
Follow the default behavior
1
Favor the next set of orbitals over the current set
RECOMMENDATION:
Use the default, and only worry about this if the system involves multiple Lewis structures
$molecule 0 1 -- 0 2 H 0.7493679 0.000000 0.44243272 -- 0 3 O 0.0000000 0.000000 -0.16535063 -- 0 2 H -0.7493679 0.000000 0.44243272 $end $rem jobtype sp method wB97X-V basis def2-TZVP ecp def2-ecp mem_total 64000 mem_static 4000 sym_ignore true symmetry false scf_algorithm DIIS scf_convergence 9 thresh 14 max_scf_cycles 1000 oslo 6 $end $loco autosad 1 ibo_mem 500 print_level 0 oslo_popanaltype 2 $end