No abstract

A new method for calculating the ground state electron density of interacting molecules is presented. The supermolecule electron density is obtained using an iterative procedure. At each step the electron density of one molecule is calculated using previously introduced Kohn-Sham equations with constrained electron density. These equations contain terms representing the coupling between constrained and non-constrained electron densities. The coupling terms also involve a new functional, namely the non-additive kinetic energy functional that is not present in the original Kohn-Sham method. Its first-principles analytical form in not yet known. We examine the analytical form of this functional derived from Thomas-Fermi theory. The electron density obtained is compared with that calculated using the original Kohn-Sham method applied to the supermolecule. Good agreement has been found for a broad range of electron density overlaps.

Although density functional theory (DFT) is more and more commonly used as a very efficient tool for the study of molecules and bulk materials, its applications to weakly bonded systems remain rather sparse in the literature, except studies that consider hydrogen bonding. It is, however, of essential interest to be able to correctly describe weaker van der Waals complexes. This prompted us to investigate more precisely the reliability of several widely-used functionals. The equilibrium geometries and the binding energies of C6H6···X (X = O2, N2, or CO) complexes are determined within the standard Kohn−Sham approach of DFT using different exchange−correlation functionals and at the MP2 level of theory for comparison. It is comprehensively concluded that extreme care must be taken in the choice of the functional since only those that behave properly at large and intermediate values of the reduced density gradient s give relevant results. The PW91 exchange functional, the enhancement factor of which does not diverge at increasing s, appears as the most reliable for the studied systems. It is furthermore demonstrated that the quality of the DFT results is determined by the exchange energy component of the total energy functional.

Ground-state properties of a linear hydrogen-bonded FH...NCH complex are studied by means of the ‘‘freeze-and-thaw’’ cycle of Kohn–Sham Equations with constrained electron density (KSCED) [T. A. Wesolowski and J. Weber, Chem. Phys. Lett. 248, 71, (1996)]. For several geometries of the complex, the electron density and the total energy are compared to the ones obtained by means of the standard Kohn–Sham calculations. The comparisons are made to assess the accuracy of several gradient dependent approximate kinetic energy functionals applied in the KSCED equations. It was found that the closest results to the Kohn–Sham ones were obtained with the functional whose analytical form was proposed by Perdew and Wang for exchange energy [J. P. Perdew and Y. Wang in Electronic Structure of Solids ’91, edited by P. Ziesche and H. Eschrig (Academie Verlag, Berlin, 1991), p. 11] and parametrized by Lembarki and Chermette for kinetic energy [A. Lembarki and H. Chermette, Phys. Rev. A 50, 5328 (1994)]. Around the interaction energy minimum as well as for larger intermolecular distances, the ‘‘freeze-and-thaw’’ cycle of KSCED leads to very similar potential energy surface as the standard supermolecule Kohn–Sham calculations.