The virtual high-throughput screening framework of the Harvard Clean Energy Project allows for the computational assessment of candidate structures for organic electronic materials -in particular photovoltaics -at an unprecedented scale. We report the most promising compounds that have emerged after studying 2.3 million molecular motifs by means of 150 million density functional theory calculations. Our top candidates are analyzed with respect to their structural makeup in order to identify important building blocks and extract design rules for efficient materials. An online database of the results is made available to the community.

Semi-local approximations to the density functional for the exchangecorrelation energy of a many-electron system necessarily fail for lobed one-electron densities, including not only the familiar stretched densities but also the less familiar but closely-related noded ones. The Perdew-Zunger (PZ) self-interaction correction (SIC) to a semi-local approximation makes that approximation exact for all oneelectron ground-or excited-state densities and accurate for stretched bonds. When the minimization of the PZ total energy is made over real localized orbitals, the orbital densities can be noded, leading to energy errors in many-electron systems. Minimization over complex localized orbitals yields nodeless orbital densities, which reduce but typically do not eliminate the SIC errors of atomization energies. Other errors of PZ SIC remain, attributable to the loss of the exact constraints and appropriate norms that the semi-local approximations satisfy, and suggesting the need for a generalized SIC. These conclusions are supported by calculations for oneelectron densities, and for many-electron molecules. While PZ SIC raises and improves the energy barriers of standard generalized gradient approximations (GGA's) and meta-GGA's, it reduces and often worsens the atomization energies of molecules. Thus PZ SIC raises the energy more as the nodality of the valence localized orbitals increases from atoms to molecules to transition states. PZ SIC is applied here in particular to the SCAN meta-GGA, for which the correlation part is already self-interaction-free. That property makes SCAN a natural first candidate for a generalized SIC.

The Fermi-Löwdin orbital self-interaction correction (FLO-SIC) methodology is applied to atoms and molecules from the standard G2-1 test set. For the first time FLO-SIC results for the GGA-type PBE functional are presented. In addition, examples where FLO-SIC like any proper SIC provides qualitative improvements compared to standard DFT functionals are discussed in detail: the dissociation limit for H + 2 , the step-wise linearity behavior for fractional occupation, as well as the significant reduction of the error of static polarizabilities. Further, ionization potentials and enthalpies of formation obtained by means of the FLO-SIC DFT method are compared to other SIC variants and experimental values. The self-interaction correction gives significant improvements if used with the LDA functional but shows worse performance in case of enthalpies of formation if the PBE-GGA functional is used. The errors are analyzed and the importance of the overbinding of hydrogen is discussed.

First principles calculations using density functional theory (DFT) have been performed to investigate the electronic and magnetic properties of DUT-8(Ni) (DUT - Dresden University of Technology). This flexible metal-organic framework (MOF) exists in two crystalline forms: DUT-8(Ni)open and DUT-8(Ni)closed. To identify the energetically favoured magnetic ordering, the density of states (DOS) and the energy difference between a low-spin (LS) and a high-spin (HS) coupling ΔELS-HS for those crystalline structures have been computed. Calculations on supercells have been carried out to include a variety of different magnetic couplings beyond a single unit cell. Several molecular model systems have been employed to further investigate the magnetic behaviour by introducing a diversity of chemical environments to the magnetic centers. The magnetic ground state of both crystalline structures has been found to be the low-spin state (S = 0). This low-spin ordering can be seen in the DOS as well as from ΔELS-HS calculations. Additionally, the calculations on the supercells confirm that the local character of the ordering (i.e. within the Ni dimers) is the most favoured one. However, the model systems indicate a change from the low-spin (S = 0) to a high-spin (S ≠ 0) ordering by introducing certain alterations into the chemical environment. Such alterations could be incorporated into the crystalline systems which should lead to similar results.