Accurate Dissociation of Chemical Bonds Using DFT-in-DFT Embedding Theory with External Orbital Orthogonality

Our recent density functional theory (DFT)-in-DFT embedding protocol, which enforces intersubsystem (or external orbital) orthogonality, is used for the first time to investigate covalent bond dissociation and is shown to do so accurately. Full potential energy curves for the dissociation of a H-O bond in H₂O and the C-C bond in H₃C-CH₃ have been constructed using the new embedding method, as have the challenging ionic bonds in LiH and LiF, and were found to match the reference Kohn-Sham (KS)-DFT curves to at least one part in 10⁶. The added constraint of external orbital orthogonality allows for the formulation of an embedding protocol that does not rely on approximate kinetic energy functionals for the evaluation of the so-called nonadditive kinetic potential, does not introduce compensatory potentials, and does not require a total system calculation at any stage. The present work extends the demonstrated applicability of the external orthogonality variant of embedding theory by more than a factor of 2 to the interaction strength range of strong single bonds. In particular, it is demonstrated that homolytic cleavage of both covalent and ionic bonds into radicals can be accomplished.