X

Search Results

Searching....

10.2 Wave Function Analysis

10.2.6 Oxidation State Localized Orbitals

(July 14, 2022)

Oxidation State Localized Orbitals (OSLO) is a new localization scheme focused on molecular fragments for the purposes of oxidation state assignment. 393 Gimferrer M. et al.
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

di=[A(NAi)2]-1 (10.9)

where the NAi is the population of the ith orbital on center A and the summation runs over all centers. When an orbital is entirely localized on a given center, its delocalization measure is di=1. If the orbital is perfectly delocalized among two centers A and B then NA=NB=1/2 and the delocalization measure is di=2, 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 Pipek J., Mezey P. G.
J. Chem. Phys.
(1989), 90, pp. 4916.
Link
However, for our purposes, out of the localized orbitals generated from fragment A with low delocalization measure, we are interested in those that are also highly localized on fragment A. Defining the fragment orbital localization index (FOLI),

diA=(diNA)1/2, (10.10)

it is easy to see diA tends to unity for fragment A when orbital i is perfectly localized on A (di=NA=1), tends to 2 when the localized orbital is perfectly delocalized over two fragments (di=2,NA=0.5), and gradually increases as the orbital becomes more delocalized as well as less centered on fragment A. Thus, among the redundant set of Mnocc 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 c Use threshold of c/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 abcdf, 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:
       abcd, 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

Example 10.2  An OSLO job to analyze the oxidation state (or formal charges) of the water molecule, H2O

$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