Efficient exact exchange approximations in density-functional theory.

Two approaches to approximate the Slater potential component of local exact exchange of density-functional theory are investigated. The first approach employs density fitting of the electrostatic potential integrals over two occupied orbitals and the other approach approximates the "exact" Slater potential with the potential derived from the Becke-Roussel [Phys. Rev. A. 39, 3761 (1989)] model of the exchange hole. In both cases significant time savings can be achieved for larger systems compared to the calculation of the numerical Slater potential. It is then analyzed how well the orbitals obtained from the various total exchange potentials reproduce Hartree-Fock energies and molecular properties. A large range of atoms and small molecules has been utilized, including the three DNA bases adenine, thymine, and cytosine.

Local correlation potentials from Brueckner coupled-cluster theory.

Local correlation potentials have been obtained from the nonlocal Brueckner coupled-cluster correlation potentials for the rare-gas atoms He, Ne, and Ar and the CO molecule. It is shown that the local correlation potential can mainly be expressed as a sum of two components: a "pure" correlation part and a relaxation contribution. While the total correlation potentials show an oscillating behavior near the nuclei, indicating the atomic shell structure, their components decrease rather monotonously, with a step structure in case of Ne and Ar. By looking at the determinantal overlap and one-electron properties it has been found that the orbitals obtained from these local potentials form a determinant which very well corresponds with the Brueckner determinant. Thus the previously found closeness between the Hartree-Fock determinant and the exchange-only Kohn-Sham determinant [Della Sala and Gorling, J. Chem. Phys. 115, 5718 (2001)] is confirmed also for the correlated case.

New CO-CO interaction potential tested by rovibrational calculations.

A four-dimensional potential energy surface (PES) for the CO dimer consisting of rigid molecules has been calculated, using a scheme that combines density functional theory to describe the monomers and symmetry adapted perturbation theory for the interaction energy (DFT-SAPT). The potential is fitted in terms of analytic functions, and the fitted potential is used to compute the lowest rovibrational states of the dimer. The quality of the PES is comparable to that of a previously published surface, which was calculated with the coupled cluster single double and perturbative triples [CCSD(T)] method. It is shown that a weighted average of the DFT-SAPT and the CCSD(T) potential gives results that are in very good agreement with experimental data, for both ((12)CO)(2) and ((13)CO)(2). The relative weight was determined by adjusting the energy gap between the origins of the lowest two stacks of rotational levels of ((12)CO)(2) to the measured value.

Density-functional theory-symmetry-adapted intermolecular perturbation theory with density fitting: a new efficient method to study intermolecular interaction energies.

The previously developed DFT-SAPT approach, which combines symmetry-adapted intermolecular perturbation theory (SAPT) with a density-functional theory (DFT) representation of the monomers, has been implemented by using density fitting of two-electron objects. This approach, termed DF-DFT-SAPT, scales with the fifth power of the molecular size and with the third power upon increase of the basis set size for a given dimer, thus drastically reducing the cost of the conventional DFT-SAPT method. The accuracy of the density fitting approximation has been tested for the ethyne dimer. It has been found that the errors in the interaction energies due to density fitting are below 10(-3) kcal/mol with suitable auxiliary basis sets and thus one or two orders of magnitude smaller than the errors due to the use of a limited atomic orbital basis set. An investigation of three prominent structures of the benzene dimer, namely, the T shaped, parallel displaced, and sandwich geometries, employing basis sets of up to augmented quadruple-zeta quality shows that DF-DFT-SAPT outperforms second-order Moller-Plesset theory (MP2) and gives total interaction energies which are close to the best estimates inferred from combining the results of MP2 and coupled-cluster theory with single, double, and perturbative triple excitations.