Equilibrium Geometries of Noncovalently Bound Intermolecular Complexes Derived from Subsystem Formulation of Density Functional Theory

Abstract

Abstract:The subsystem formulation of density functional theory is used to obtain equilibrium geometries and interaction energies for a representative set of noncovalently bound intermolecular complexes. The results are compared with literature benchmark data. The range of applicability of two considered approximations to the exchange-correlation-and nonadditive kinetic energy components of the total energy is determined. Local density approximation, which does not involve any empirical parameters, leads to excellent intermolecular equilibrium distances for hydrogen-bonded complexes (maximal error 0.13 Å for NH 3 -NH 3 ). It is a method of choice for a wide class of weak intermolecular complexes including also dipole-bound and the ones formed by rare gas atoms or saturated hydrocarbons. The range of applicability of the chosen generalized gradient approximation, which was shown in our previous works to lead to good interaction energies in such complexes, where π-electrons are involved in the interaction, remains limited to this group because it improves neither binding energies nor equilibrium geometries in the wide class of complexes for which local density approximation is adequate. An efficient energy minimization procedure, in which optimization of the geometry and the electron density of each subsystem is made simultaneously, is proposed and tested.

Embedded density functional theory for covalently bonded and strongly interacting subsystems

“…V, we discuss a technique to robustly implement the ZMP step for NAKP calculations in large systems. We note that the calculation of gradients using the e-DFT framework is well established, 21,37 although demonstration of gradient calculations with the EE method is left for future work.…”

Chromophore-specific theoretical spectroscopy: From subsystem density functional theory to mode-specific vibrational spectroscopy

“…[192]): FDE employing LDA functionals for both E xc and T nadd s leads to very good interaction energies for weakly overlapping densities, such as in rare gas dimers or the CH 4 · · · CH 4 complex. For interactions involving hydrogen-bonded systems, GGA-type approximations result in closer agreement with KS-DFT calculations [119], while the study in Ref.…”

The energy-differences based exact criterion for testing approximations to the functional for the kinetic energy of non-interacting electrons

“…As far as interaction energies and equilibrium geometries are concerned, local density approximation applied in Kohn-Sham calculations is known to lead to unsatisfactory results whereas the same approximation applied in the subsystem based calculations simultaneously for exchange-correlation-and non-additive kinetic components of the total energy leads to good results in many cases [22,23]. This is probably the result of mutual cancellation of errors in the exchange-and kinetic energy contributions to the interaction energy in the weakly overlapping case.…”