Homolytic molecular dissociation in natural orbital functional theory

The dissociation of diatomic molecules of the 14-electron isoelectronic series N(2), O(2)(2+), CO, CN(-) and NO(+) is examined using the Piris natural orbital functional. It is found that the method describes correctly the dissociation limit yielding an integer number of electrons on the dissociated atoms, in contrast to the fractional charges obtained when using the variational two-particle reduced density matrix method under the D, Q and G positivity necessary N-representability conditions. The chemistry of the considered systems is discussed in terms of their dipole moments, natural orbital occupations and bond orders as well as atomic Mulliken populations at the dissociation limit. The values obtained agree well with accurate multiconfigurational wave function based CASSCF results and the available experimental data.