Oxidation states represent the ionic distribution of charge in a molecule and are significant in tracking redox reactions and understanding chemical bonding. While effective algorithms already exist based on formal Lewis structures as well as using localized orbitals, they exhibit differences in challenging cases where effects such as redox noninnocence are at play. Given a density functional theory (DFT) calculation with chosen total charge and spin multiplicity, this work reports a new approach to obtaining fragment-localized orbitals that is termed oxidation state localized orbitals (OSLO), together with an algorithm for assigning the oxidation state using the OSLOs and an associated fragment orbital localization index (FOLI). Evaluating the FOLI requires fragment populations, and for this purpose a new version of the intrinsic atomic orbital (IAO) scheme is introduced in which the IAOs are evaluated using a reference minimal basis formed from on-the-fly superposition of atomic density (IAO-AutoSAD) calculations in the target basis set and at the target level of theory. The OSLO algorithm is applied to a range of challenging cases including high valent metal oxide complexes, redox noninnocent NO and dithiolate transition metal complexes, a range of carbene-containing TM complexes, and other examples including the potentially inverted ligand field in [Cu(CF₃)₄]^-. Across this range of cases, OSLO produces generally satisfactory results. Furthermore, in borderline cases, the OSLOs and associated FOLI values provide direct evidence of the emergence of covalent interactions between fragments that nicely complements existing approaches.

The oxidation state (OS) of metals and ligands in inorganic complexes may be defined by carefully curated rules, such as from IUPAC, or by computational procedures such as the effective oxidation state (EOS) or localized orbital bonding analysis (LOBA). Such definitions typically agree for systems with simple ionic bonding and innocent ligands but may disagree as the boundary between ionic and covalent bonds is approached, or as the role of ligand noninnocence becomes nontrivial, or high oxidation states of metals are supported by heavy dative bonding, and so on. This work systematically compares IUPAC, EOS, and LOBA across a series of complexes where OS assignment is challenging. These systems include high-valent transition metal oxides, transition metal complexes with noninnocent ligands such as dithiolate and nitrosyl, metal sulfur dioxide adducts, and two transition metal carbene complexes. The differences in OS assignment by the three methods are carefully discussed, demonstrating the synergy between EOS and LOBA. In addition, a clarity index for LOBA OS assignments is introduced that provides an indication of whether or not its predictions are close to the ionic-covalent boundary.

Computational quantum chemistry can be more than just numerical experiments when methods are specifically adapted to investigate chemical concepts. One important example is the development of energy decomposition analysis (EDA) to reveal the physical driving forces behind intermolecular interactions. In EDA, typically the interaction energy from a good-quality density functional theory (DFT) calculation is decomposed into multiple additive components that unveil permanent and induced electrostatics, Pauli repulsion, dispersion, and charge-transfer contributions to noncovalent interactions. Herein, we formulate, implement, and investigate decomposing the forces associated with intermolecular interactions into the same components. The resulting force decomposition analysis (FDA) is potentially useful as a complement to the EDA to understand chemistry, while also providing far more information than an EDA for data analysis purposes such as training physics-based force fields. We apply the FDA based on absolutely localized molecular orbitals (ALMOs) to analyze interactions of water with sodium and chloride ions as well as in the water dimer. We also analyze the forces responsible for geometric changes in carbon dioxide upon adsorption onto (and activation by) gold and silver anions. We also investigate how the force components of an EDA-based force field for water clusters, namely MB-UCB, compare to those from force decomposition analysis.