A method is presented for localizing molecular orbitals, based on diagonalizing subunits of the density matrix. First, nonbonding orbitals are found by diagonalizing the monatomic subunits; then, diatomic a or -r bonding and antibonding orbitals are obtained from the diatomic subunits for all bonded pairs of atoms; finally, the delocalized ir-orbitals for particular chromophores are found by projecting the first set out of the self-consistent field (SCF) Hamiltonian. The results show good general agreement with other localization methods, with advantages in the ability to display group orbitals in complex molecules which most closely resemble the SCF orbitals for simple prototypes. The impetus for obtaining localized orbitals is well recognized and much has been written about the applications and uses of localized orbitals (1-3).A set of localized orbitals may be obtained by a direct calculation of self-consistent field (SCF) localized molecular orbitals (4) or, given an arbitrary set of delocalized orbitals, one may construct a unitary matrix that will transform the delocalized orbitals into localized orbitals. The two most commonly used procedures for doing the latter are those of Boys (5) and Edmiston and Ruedenberg (6). Closely related methods (7, 8), improvements in efficiency (9), and general discussions of the various methods have been presented elsewhere (10-14).In this work a method is presented that is based on the density matrix formalism and can be used with either ab initio or semi-empirical programs. In the latter case the required input consists of the overlap matrix, the molecular Hamiltonian, and the first-order density matrix.
METHODThe localization process is based on the diagonalization of appropriate subunits of the density matrix (15). First, the nonbonding and vacant atomic orbitals are obtained by diagonalizing the individual monatomic components, which are individually delineated by the basis orbitals for each atom. Ideally, a nonbonding orbital would be associated with the eigenvalue 2 and a vacant orbital would have the value 0. In practice it has been found convenient to take all orbitals with ni > 1.6 as nonbonding and those for which nj < 0.4 as vacant. The remaining eigenvectors are hybrids or linear combinations thereof associated with bonding to nearest neighbors.Although it is tempting to use the remaining eigenvectors as hybrids for localized bonding and antibonding orbitals, this process is defeated by the near degeneracy (ideally, 1) for such orbitals, and extensive mixing occurs. The problem is avoided by the second step.The diatomic portions for all pairs of atoms are isolated and diagonalized. Again