We present the procedure for transforming delocalized molecular orbitals into the localized property-optimized orbitals (LPOs) designed for building the most accurate, in the Frobenius norm sense, approximation to the first-order reduced density matrix in form of the sum of localized monoatomic and diatomic terms. In this way, a decomposition of molecular properties into contributions associated with individual atoms and the pairs of atoms is obtained with the a priori known upper bound for the decomposition accuracy. Additional algorithm is proposed for obtaining the set of "the Chemist's LPOs" (CLPOs) containing a single localized orbital, with nearly double occupancy, per a pair of electrons. CLPOs form an idealized Lewis structure optimized for the closest possible reproduction of one-electron properties derived from the original many-electron wavefunction. The computational algorithms for constructing LPOs and CLPOs from a general wavefunction are presented and their implementation within the open-source freeware program JANPA (http://janpa.sourceforge.net/) is discussed. The performance of the proposed procedures is assessed using the test set of density matrices of 33 432 small molecules obtained at both Hartree-Fock and second-order Moller-Plesset theory levels and excellent agreement with the chemist's Lewis-structure picture is found. K E Y W O R D S algorithms and software, chemist's Lewis-structure picture, localized orbitals, properties decomposition, quantum-chemical methods
| INTRODUCTIONLocalized molecular orbitals resulting from a unitary transformation of occupied canonical molecular orbitals (MOs) [1] play essential role in physical chemistry as the "building blocks" or "descriptors" through which the complicated electronic structure of atoms and molecules can be modeled or interpreted in a comprehensible way. In addition, the localized orbitals concentrated in a limited spatial region of a molecule proved useful in making the high-level correlated quantum-chemical methods more computationally tractable. [2][3][4][5][6][7][8] Although the concept of an orbital itself has been much methodologically debatable and is too "fuzzy" [34] to be defined more precisely than just as the function of coordinates of a single electron somehow related to the system's wavefunction or electron density, this concept still remains virtually the best one proposed so far for moving the ideas of valence electrons and electron pairs (including bonding and antibonding orbitals, lone pairs etc.), which are the key elements of the chemist's Lewis-structure picture, [35] from qualitative to a quantum-mechanical ground. From this perspective, the localized orbitals are used to decompose (typically in an approximate manner) true many-electron wavefunction of the molecule into the components allowing a chemically meaningful interpretation.To transform delocalized orbitals (either the canonical MOs obtained as a solutions of self-consistent field Hartree-Fock or Kohn-Sham equations, [36] or the Lowding natural orbitals [37] ...