For the first time, a complete implementation of coupled perturbed density functional theory (CPDFT) for the calculation of NMR spin–spin coupling constants (SSCCs) with pure and hybrid DFT is presented. By applying this method to several hydrides, hydrocarbons, and molecules with multiple bonds, the performance of DFT for the calculation of SSCCs is analyzed in dependence of the XC functional used. The importance of electron correlation effects is demonstrated and it is shown that the hybrid functional B3LYP leads to the best accuracy of calculated SSCCs. Also, CPDFT is compared with sum-over-states (SOS) DFT where it turns out that the former method is superior to the latter because it explicitly considers the dependence of the Kohn–Sham operator on the perturbed orbitals in DFT when calculating SSCCs. The four different coupling mechanisms contributing to the SSCC are discussed in connection with the electronic structure of the molecule.

Unrestricted density functional theory (UDFT) can be used for the description of open-shell singlet (OSS) biradicals provided a number of precautions are considered. Biradicals that require a two-determinantal wave function (e.g. OSS state of carbenes) cannot be described by UDFT for principal reasons. However, if the overlap between the open-shell orbitals is small (the single electrons are located at different atomic centers) errors become small and, then, the principal failure of UDFT in these cases is not apparent and may even be disguised by the fact that UDFT has the advantage of describing spin polarization better than any restricted open shell DFT method. In the case of OSS biradicals with two- or multiconfigurational character (but a onedeterminantal form of the leading configuration), reasonable results can be obtained by broken-symmetry (BS)-UDFT, however in each case this has to be checked. In no case is it reasonable to lower the symmetry of a molecule to get a suitable UDFT description. Hybrid functionals such as B3LYP perform better than pure DFT functionals in BS-UDFT calculations because the former reduce the self-interaction error of DFT exchange functionals, which mimics unspecified static electron correlation effects, so that the inclusion of specific static electron correlation effects via the form of the wavefunction becomes more effective.

Self-interaction corrected density functional theory was used to determine the self-interaction error for dissociating one-electron bonds. The self-interaction error of the unpaired electron mimics nondynamic correlation effects that have no physical basis where these effects increase for increasing separation distance. For short distances the magnitude of the self-interaction error takes a minimum and increases then again for decreasing R. The position of the minimum of the magnitude of the self-interaction error influences the equilibrium properties of the one-electron bond in the radical cations H2+ (1), B2H4+ (2), and C2H6+ (3), which differ significantly. These differences are explained by hyperconjugative interactions in 2 and 3 that are directly reflected by the self-interaction error and its orbital contributions. The density functional theory description of the dissociating radical cations suffers not only from the self-interaction error but also from the simplified description of interelectronic exchange. The calculated differences between ionic and covalent dissociation for 1, 2, and 3 provide an excellent criterion for determining the basic failures of density functional theory, self-interaction corrected density functional theory, and other methods. Pure electronic, orbital relaxation, and geometric relaxation contributions to the self-interaction error are discussed. The relevance of these effects for the description of transition states and charge transfer complexes is shown. Suggestions for the construction of new exchange-correlation functionals are given. In this connection, the disadvantages of recently suggested self-interaction error-free density functional theory methods are emphasized.