Exponential decay of spin-spin correlation between distant defect states produced by contour hydrogenation of polycyclic aromatic hydrocarbon molecules

The first few low-lying spin states of alternant polycyclic aromatic hydrocarbon (PAH) molecules of several shapes showing defect states induced by contour hydrogenation have been studied both by ab initio methods and by a precise numerical solution of Pariser-Parr-Pople (PPP) interacting model. In accordance with Lieb's theorem, the ground state shows a spin multiplicity equal to one for balanced molecules, and it gets larger values for imbalanced molecules (that is, when the number of π electrons on both subsets is not equal). Furthermore, we find a systematic decrease of the singlet-triplet splitting as a function of the distance between defects, regardless of whether the ground state is singlet or triplet. For example, a splitting smaller than 0.001 eV is obtained for a medium size C46H28 PAH molecule (di-hydrogenated [11]phenacene) showing a singlet ground state. We conclude that π electrons unbound by lattice defects tend to remain localized and unpaired even when long-range Coulomb interaction is taken into account. Therefore they show a biradical character (polyradical character for more than two defects) and should be studied as two or more local doublets. The implications for electron transport are potentially important since these unpaired electrons can trap traveling electrons or simply flip their spin at a very small energy cost.