Intermolecular potentials of the silane dimer calculated with Hartree-Fock theory, Møller-Plesset perturbation theory, and density functional theory

We have calculated the intermolecular interaction potentials of the silane dimer at the D3d conformation using the Hartree-Fock (HF) self-consistent theory, the correlation-corrected second-order M?ller-Plesset (MP2) perturbation theory, and the density functional theory (DFT) with 108 functionals chosen from the combinations of 9 exchange and 12 correlation functionals. Single-point coupled cluster [CCSD(T)] calculations have also been carried out to calibrate the correlation effect. The HF calculations yield unbound potentials largely because of the exchange-repulsion interaction. In the MP2 calculations, the basis set effects on the repulsion exponent, the equilibrium bond length, the binding energy, and the asymptotic behavior of the calculated intermolecular potentials have been thoroughly studied. We have employed basis sets from the Slater type orbitals fitted with Gaussian functions (STO-nG, n = 3 approximately 6), Pople's medium size basis sets [up to 6-311++G(3df,3pd)], to Dunning's correlation consistent basis sets (cc-pVXZ and aug-cc-pVXZ, X = D, T, Q). With increasing basis size, the repulsion exponent and the equilibrium bond length converge at the 6-31G** basis set and the 6-311++G(3d,3p) basis set, respectively, while a large basis set (aug-cc-pVTZ) is required to converge the binding energy at a chemical accuracy ( approximately 0.05 kcal/mol). Up to the largest basis set used, the asymptotic dispersion coefficient has not converged to the expected C6 value from molecular polarizability calculations. We attribute the slow convergence partly to the inefficacy of using the MP2 calculations with Gaussian type functions to model the asymptotic behavior. Both the basis set superposition error (BSSE) corrected and uncorrected results are presented to emphasize the importance of including such corrections. Only the BSSE corrected results systematically converge to the expected potential curve with increasing basis size. The DFT calculations generate a wide range of interaction patterns, from purely unbound to strongly bound, underestimating or overestimating the binding energy. The binding energies calculated using the OPTXHCTH147, PBEVP86, PBEP86, PW91TPSS, PW91PBE, and PW91PW91 functionals and the equilibrium bond lengths calculated using the MPWHCTH93, TPSSHCTH, PBEVP86, PBEP86, PW91TPSS, PW91PBE, and PW91PW91 functionals are close to the MP2 results using the 6-311++G(3df,3pd) basis set. A correlation between the calculated DFT potentials and the exchange and correlation enhancement factors at the low-density region has been elucidated. The asymptotic behaviors of the DFT potentials are also analyzed.     

甲烷晶体的晶格能和弹性性质: 不同方法及泛函的评估

通过对甲烷晶体进行结构、晶格能和弹性特性的研究, 评估了不包含和包含色散能量修正的密度泛函理论的性能. 我们分别利用不包含色散能量修正的密度泛函理论(DFT) (包含不同的标准泛函和杂化泛函)和包含色散能量修正的密度泛函理论(DFT-D)计算了甲烷晶体特性, 并与实验作对比. 尽管DFT-D 与传统密度泛函理论及杂化密度泛函理论相比, 修正了甲烷晶体中的范德华(vdW)相互作用, 但是一些修正方案过分修正了这种相互作用. 因此, 人们在使用DFT-D方法时务必谨慎.     

Critical examination of the supermolecule density functional theory calculations of intermolecular interactions

The results of calculations employing twelve different combinations of exchange and correlation functionals are compared with results of ab initio calculations for two different configurations of the water dimer and three different configurations of the thymine-adenine complex. None of the density functional theory (DFT) treatments could properly reproduce the results of coupled-cluster calculations for all configurations examined. The DFT approaches perform well when the interaction energy is dominated by the electrostatic component and the dispersion energy is less important. Two mechanisms that compensate for the missing dispersion component were identified. The first one is the decrease of the magnitude of the intermolecular exchange-repulsion and the second one is the increase of the magnitude of the attractive deformation energy. For some functionals both effects are observed together, but for some other ones only the second effect occurs. The three correlation functionals that were examined were found to make only very small contributions to the deformation energy. The examination of angular and distance dependence of the interactions shows that the currently available DFT approaches are not suitable for developing intermolecular potential energy surfaces. They could however be used to find global minima on potential energy surfaces governed by intermolecular electrostatic interactions. Additional single point ab initio calculations are recommended as the means of validating optimized structures.     

Density, structure, and dynamics of water: the effect of van der Waals interactions

It is known that ab initio molecular dynamics (AIMD) simulations of liquid water at ambient conditions, based on the generalized gradient approximation (GGA) to density functional theory (DFT), with commonly used functionals fail to produce structural and diffusive properties in reasonable agreement with experiment. This is true for canonical, constant temperature simulations where the density of the liquid is fixed to the experimental density. The equilibrium density, at ambient conditions, of DFT water has recently been shown by Schmidt et al. [J. Phys. Chem. B, 113, 11959 (2009)] to be underestimated by different GGA functionals for exchange and correlation, and corrected by the addition of interatomic pair potentials to describe van der Waals (vdW) interactions. In this contribution we present a DFT-AIMD study of liquid water using several GGA functionals as well as the van der Waals density functional (vdW-DF) of Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)]. As expected, we find that the density of water is grossly underestimated by GGA functionals. When a vdW-DF is used, the density improves drastically and the experimental diffusivity is reproduced without the need of thermal corrections. We analyze the origin of the density differences between all the functionals. We show that the vdW-DF increases the population of non-H-bonded interstitial sites, at distances between the first and second coordination shells. However, it excessively weakens the H-bond network, collapsing the second coordination shell. This structural problem is partially associated to the choice of GGA exchange in the vdW-DF. We show that a different choice for the exchange functional is enough to achieve an overall improvement both in structure and diffusivity.     

Exchange and correlation energy in density functional theory

Several different versions of density functional theory (DFT) that satisfy Hohenberg–Kohn theorems are characterized by different definitions of a reference or model state determined by an N‐electron ground state. A common formalism is developed in which exact Kohn–Sham equations are derived for standard Kohn–Sham theory, for reference‐state density functional theory, and for unrestricted Hartree–Fock (UHF) theory considered as an exactly soluble model Hohenberg–Kohn theory. A natural definition of exchange and correlation energy functionals is shown to be valid for all such theories. An easily computed necessary condition for the locality of exchange and correlation potentials is derived. While it is shown that in the UHF model of DFT the optimized effective potential (OEP) exchange satisfies this condition by construction, the derivation shows that this condition is not, in general, sufficient to define an exact local exchange potential. It serves as a test to eliminate proposed local potentials that are not exact for ground states. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 521–525, 2000     

Comparison of DFT methods for molecular orbital eigenvalue calculations

We report how closely the Kohn-Sham highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) eigenvalues of 11 density functional theory (DFT) functionals, respectively, correspond to the negative ionization potentials (-IPs) and electron affinities (EAs) of a test set of molecules. We also report how accurately the HOMO-LUMO gaps of these methods predict the lowest excitation energies using both time-independent and time-dependent DFT (TD-DFT). The 11 DFT functionals include the local spin density approximation (LSDA), five generalized gradient approximation (GGA) functionals, three hybrid GGA functionals, one hybrid functional, and one hybrid meta GGA functional. We find that the HOMO eigenvalues predicted by KMLYP, BH&HLYP, B3LYP, PW91, PBE, and BLYP predict the -IPs with average absolute errors of 0.73, 1.48, 3.10, 4.27, 4.33, and 4.41 eV, respectively. The LUMOs of all functionals fail to accurately predict the EAs. Although the GGA functionals inaccurately predict both the HOMO and LUMO eigenvalues, they predict the HOMO-LUMO gap relatively accurately (approximately 0.73 eV). On the other hand, the LUMO eigenvalues of the hybrid functionals fail to predict the EA to the extent that they include HF exchange, although increasing HF exchange improves the correspondence between the HOMO eigenvalue and -IP so that the HOMO-LUMO gaps are inaccurately predicted by hybrid DFT functionals. We find that TD-DFT with all functionals accurately predicts the HOMO-LUMO gaps. A linear correlation between the calculated HOMO eigenvalue and the experimental -IP and calculated HOMO-LUMO gap and experimental lowest excitation energy enables us to derive a simple correction formula.     

Evaluation of DFT methods for computing the interaction energies of homomolecular and heteromolecular dimers of monosubstituted benzene

We present density functional theory (DFT) interaction energies for the sandwich and T‐shaped conformers of substituted benzene dimers. The DFT functionals studied include TPSS, HCTH407, B3LYP, and X3LYP. We also include Hartree–Fock (HF) and second‐order Møller–Plesset perturbation theory calculations (MP2), as well as calculations using a new functional, P3LYP, which includes PBE and HF exchange and LYP correlation. Although DFT methods do not explicitly account for the dispersion interactions important in the benzene–dimer interactions, we find that our new method, P3LYP, as well as HCTH407 and TPSS, match MP2 and CCSD(T) calculations much better than the hybrid methods B3LYP and X3LYP methods do. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Some thoughts about the stability and reliability of commonly used exchange–correlation functionals – coverage of dynamic and nondynamic correlation effects

 The self-interaction error (SIE) of commonly used density functional theory (DFT) exchange functionals mimics long-range (nondynamic) pair correlation effects in an unspecified way. Slater exchange suffers from a larger SIE and, therefore, covers more nondynamic correlation effects than Becke exchange, which is the reason why exchange–correlation (XC) functionals based on Slater exchange lead to stabler restricted DFT solutions than those based on Becke exchange. However, the stability of an XC functional does not guarantee higher accuracy. On the contrary, if system-specific nondynamic correlation effects have to be introduced via the form of the wave function, these will be suppressed by nondynamic correlation effects already covered by the exchange functional. Hybrid functionals suffer less from the SIE and, therefore, cover a smaller number of nondynamic electron correlation effects. Accordingly, they are better suited when nondynamic electron correlation has to be introduced by the form of the wave function. It is shown that, for example, broken-symmetry unrestricted DFT calculations are more accurate when carried out with B3LYP than BLYP contrary to claims made in the literature. Received: 8 November 2001 / Accepted: 30 January 2002 / Published online: 8 April 2002     

Coupled-cluster, Möller Plesset (MP2), density fitted local MP2, and density functional theory examination of the energetic and structural features of hydrophobic solvation: water and pentane

The interaction potentials between immiscible polar and non-polar solvents are a major driving force behind the formation of liquid:liquid interfaces. In this work, the interaction energy of water-pentane dimer has been determined using coupled-cluster theory with single double (triple) excitations [CCSD(T)], 2nd order M?ller Plesset perturbation theory (MP2), density fitted local MP2 (DF-LMP2), as well as density functional theory using a wide variety of density functionals and several different basis sets. The M05-2X exchange correlation functionals exhibit excellent agreement with CCSD(T) and DF-LMP2 after taking into account basis set superposition error. The gas phase water-pentane interaction energy is found to be quite sensitive to the specific pentane isomer (2,2-dimethylpropane vs. n-pentane) and relative orientation of the monomeric constituents. Subsequent solution phase cluster calculations of 2,2-dimethylpropane and n-pentane solvated by water indicate a positive free energy of solvation that is in good agreement with available experimental data. Structural parameters are quite sensitive to the density functional employed and reflect differences in the two-body interaction energy calculated by each method. In contrast, cluster calculations of pentane solvation of H(2)O solute are found to be inadequate for describing the organic solvent, likely due to limitations associated with the functionals employed (B3LYP, BHandH, and M05-2X).     

Density functional theory and molecular clusters

Density functional theory (DFT) methods, including nonlocal density gradient terms in the exchange and correlation energy functionals, were applied to various types of molecular clusters: H-bonded, ionic, electrostatic, and London. Reliable results on the structure and stabilization energy were obtained for the first two types of cluster as long as Becke3LYP and Becke3P86 functionals and basis sets of at least DZ + P quality were used. DFT methods with currently available functionals failed completely, however, for London-type clusters, for which no minimum was found on the potential energy surfaces. DFT interaction energy exhibits the same basis set extension dependence as the Hartree-Fock (HF) interaction energy. Therefore, the Boys-Bernardi function counterpoise procedure should be employed for elimination of the DFT basis set superposition error. © 1995 John Wiley & Sons, Inc.     

Density Functional Theory as a Data Science

The development of density functional theory (DFT) functionals and physical corrections are reviewed focusing on the physical meanings and the semiempirical parameters from the viewpoint of data science. This review shows that DFT exchange‐correlation functionals have been developed under many strict physical conditions with minimizing the number of the semiempirical parameters, except for some recent functionals. Major physical corrections for exchange‐correlation function‐ als are also shown to have clear physical meanings independent of the functionals, though they inevitably require minimum semiempirical parameters dependent on the functionals combined. We, therefore, interpret that DFT functionals with physical corrections are the most sophisticated target functions that are physically legitimated, even from the viewpoint of data science.     

An examination of density functionals on aldol, Mannich and ��-aminoxylation reaction enthalpy calculations

The reaction enthalpies of aldol, Mannich and ??-aminoxylation reactions were calculated by density functional theory (DFT) using long-range-corrected (LC), hybrid B3LYP and other up-to-date functionals to show why conventional DFT including B3LYP has given poor enthalpies for these reactions. As a result, we found that long-range exchange interactions significantly affect the reaction enthalpies. We therefore proposed that the poor enthalpies of B3LYP are due to its insufficient long-range exchange effect. On the other hand, LC functionals accurately reproduce reaction enthalpies for these reactions. However, we noticed that even LC functionals present poor reaction enthalpies for specific reactions, in which many branches are produced or very small molecules such as methane molecule participate.     

Electronic structure and excitations in oligoacenes from ab initio calculations

Oligoacenes C(4n+2)H(2n+4) (n=2,...,6) are studied using a variety of ab initio methods. Density functional theory (DFT) optimized geometries were in good agreement with experiment. Vertical and adiabatic ionization potentials and electron affinities were computed with DFT and it was found that standard exchange-correlation (xc) functionals underestimate ionization potentials in oligoacenes. Possible reasons for this underestimation are discussed. Low lying electronic excitations were computed using time-dependent density functional theory, configuration interaction singles, and configuration interaction singles with approximate treatment of doubles. In agreement with earlier work, time-dependent DFT in conjunction with standard xc-energy functionals substantially underestimates the lowest (p) singlet-singlet electronic transition.     

Efficient hybrid density functional calculations in solids: assessment of the Heyd-Scuseria-Ernzerhof screened Coulomb hybrid functional

The present work introduces an efficient screening technique to take advantage of the fast spatial decay of the short range Hartree-Fock (HF) exchange used in the Heyd-Scuseria-Ernzerhof (HSE) screened Coulomb hybrid density functional. The screened HF exchange decay properties and screening efficiency are compared with traditional hybrid functional calculations on solids. The HSE functional is then assessed using 21 metallic, semiconducting, and insulating solids. The examined properties include lattice constants, bulk moduli, and band gaps. The results obtained with HSE exhibit significantly smaller errors than pure density functional theory (DFT) calculations. For structural properties, the errors produced by HSE are up to 50% smaller than the errors of the local density approximation, PBE, and TPSS functionals used for comparison. When predicting band gaps of semiconductors, we found smaller errors with HSE, resulting in a mean absolute error of 0.2 eV (1.3 eV error for all pure DFT functionals). In addition, we present timing results which show the computational time requirements of HSE to be only a factor of 2-4 higher than pure DFT functionals. These results make HSE an attractive choice for calculations of all types of solids.     

PW86–PW91 density functional calculation of vertical ionization potentials: Some implications for present‐day functionals

A total of 181 vertical ionization potentials (VIPs) of 41 molecules were calculated by density functional theory (DFT) employing the Perdew–Wang 1986 (PW86) exchange and Perdew–Wang 1991 (PW91) correlation functionals and using the aug‐cc‐pV5Z basis and experimental ground‐state geometries. The overall average absolute deviation (AAD) from experiment was found to be 0.55 eV and only 0.31 eV for linear molecules but 0.86 eV for nonplanar molecules. A number of VIPs were in error by over 2 eV. In particular, DFT performed most poorly when ionization was from an orbital with highly varying density gradients (which arise from the orbital's shape or compactness or through its density being distributed over a number of atoms). Indications are that many or all present‐day functionals suffer from the same failings. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 34–52, 2001     

A CSOV study of the difference between HF and DFT intermolecular interaction energy values: the importance of the charge transfer contribution

Intermolecular interaction energy decompositions using the Constrained Space Orbital Variation (CSOV) method are carried out at the Hartree-Fock level on the one hand and using DFT with usual GGA functionals on the other for a number of model complexes to analyze the role of electron correlation in the intermolecular stabilization energy. In addition to the overall stabilization, the results provide information on the variation, with respect to the computational level, of the different contributions to the interaction energy. The complexes studied are the water linear dimer, the N-methylformamide dimer, the nucleic acid base pairs, the benzene-methane and benzene-N2 van der Waals complexes, [Cu+ -(ImH)3]2, where "ImH" stands for the Imidazole ligand, and ImH-Zn++. The variation of the frozen core energy (the sum of the intermolecular electrostatic energy and the Pauli repulsion energy) calculated from the unperturbed orbitals of the interacting entities indicates that the intramolecular correlation contributions can be stabilizing as well as destabilizing, and that general trends can be derived from the results obtained using usual density functionals. The most important difference between the values obtained from HF and DFT computations concerns the charge transfer contribution, which, in most cases, undergoes the largest increase. The physical meaning of these results is discussed. The present work gives reference calculations that might be used to parametrize new correlated molecular mechanics potentials.     

Increasing the applicability of density functional theory. II. Correlation potentials from the random phase approximation and beyond

Density functional theory (DFT) results are mistrusted at times due to the presence of an unknown exchange correlation functional, with no practical way to guarantee convergence to the right answer. The use of a known exchange correlation functional based on wave-function theory helps to alleviate such mistrust. The exchange correlation functionals can be written exactly in terms of the density-density response function using the adiabatic-connection and fluctuation-dissipation framework. The random phase approximation (RPA) is the simplest approximation for the density-density response function. Since the correlation functional obtained from RPA is equivalent to the direct ring coupled cluster doubles (ring-CCD) correlation functional, meaning only Coulomb interactions are included, one can bracket RPA between many body perturbation theory (MBPT)-2 and CCD with the latter having all ring, ladder, and exchange contributions. Using an optimized effective potential strategy, we obtain correlation potentials corresponding to MBPT-2, RPA (ring-CCD), linear-CCD, and CCD. Using the suitable choice of the unperturbed Hamiltonian, Kohn-Sham self-consistent calculations are performed. The spatial behavior of the resulting potentials, total energies, and the HOMO eigenvalues are compared with the exact values for spherical atoms. Further, we demonstrate that the self-consistent eigenvalues obtained from these consistent potentials used in ab initio dft approximate all principal ionization potentials as demanded by ionization potential theorem.     

Comparing ab initio density-functional and wave function theories: the impact of correlation on the electronic density and the role of the correlation potential

The framework of ab initio density-functional theory (DFT) has been introduced as a way to provide a seamless connection between the Kohn-Sham (KS) formulation of DFT and wave-function based ab initio approaches [R. J. Bartlett, I. Grabowski, S. Hirata, and S. Ivanov, J. Chem. Phys. 122, 034104 (2005)]. Recently, an analysis of the impact of dynamical correlation effects on the density of the neon atom was presented [K. Jankowski, K. Nowakowski, I. Grabowski, and J. Wasilewski, J. Chem. Phys. 130, 164102 (2009)], contrasting the behaviour for a variety of standard density functionals with that of ab initio approaches based on second-order M?ller-Plesset (MP2) and coupled cluster theories at the singles-doubles (CCSD) and singles-doubles perturbative triples [CCSD(T)] levels. In the present work, we consider ab initio density functionals based on second-order many-body perturbation theory and coupled cluster perturbation theory in a similar manner, for a range of small atomic and molecular systems. For comparison, we also consider results obtained from MP2, CCSD, and CCSD(T) calculations. In addition to this density based analysis, we determine the KS correlation potentials corresponding to these densities and compare them with those obtained for a range of ab initio density functionals via the optimized effective potential method. The correlation energies, densities, and potentials calculated using ab initio DFT display a similar systematic behaviour to those derived from electronic densities calculated using ab initio wave function theories. In contrast, typical explicit density functionals for the correlation energy, such as VWN5 and LYP, do not show behaviour consistent with this picture of dynamical correlation, although they may provide some degree of correction for already erroneous explicitly density-dependent exchange-only functionals. The results presented here using orbital dependent ab initio density functionals show that they provide a treatment of exchange and correlation contributions within the KS framework that is more consistent with traditional ab initio wave function based methods.