Density functional theory augmented with an empirical dispersion term. Interaction energies and geometries of 80 noncovalent complexes compared with ab initio quantum mechanics calculations

Standard density functional theory (DFT) is augmented with a damped empirical dispersion term. The damping function is optimized on a small, well balanced set of 22 van der Waals (vdW) complexes and verified on a validation set of 58 vdW complexes. Both sets contain biologically relevant molecules such as nucleic acid bases. Results are in remarkable agreement with reference high-level wave function data based on the CCSD(T) method. The geometries obtained by full gradient optimization are in very good agreement with the best available theoretical reference. In terms of the standard deviation and average errors, results including the empirical dispersion term are clearly superior to all pure density functionals investigated-B-LYP, B3-LYP, PBE, TPSS, TPSSh, and BH-LYP-and even surpass the MP2/cc-pVTZ method. The combination of empirical dispersion with the TPSS functional performs remarkably well. The most critical part of the empirical dispersion approach is the damping function. The damping parameters should be optimized for each density functional/basis set combination separately. To keep the method simple, we optimized mainly a single factor, s(R), scaling globally the vdW radii. For good results, a basis set of at least triple-zeta quality is required and diffuse functions are recommended, since the basis set superposition error seriously deteriorates the results. On average, the dispersion contribution to the interaction energy missing in the DFT functionals examined here is about 15 and 100% for the hydrogen-bonded and stacked complexes considered, respectively.     

Photodissociation dynamics of the Kr-HBr cluster: the effect of the rare gas atom substitution

The ultraviolet photolysis dynamics of Kr-HBr(v=0) is investigated by means of wave packet calculations, focusing on the fragmentation pathway Kr-HBr+ variant Planck's over 2pi omega-->H+Kr-Br. Photolysis is simulated by starting from two different cluster initial states, namely the ground van der Waals (vdW) and an excited vdW bending state, associated with the Kr-H-Br and Kr-Br-H isomers, respectively. The results show that, for the two initial states of the cluster, the Kr-Br product yield is lower than that of Ar-Br radicals found in previous studies on Ar-HBr photolysis. Despite this decrease, the Kr-Br yield is found to be still rather high, in particular for the initial excited vdW state of Kr-HBr(v=0). In addition, the Kr-Br product state distributions exhibit a remarkably higher excitation (mainly rotational) than the corresponding Ar-Br distributions. The lower yield and higher excitation of Kr-Br as compared to Ar-Br, are attributed to a larger share of the energy available for the radical going to internal excitation in the case of the Kr-Br product. The different partition of the energy available for Kr-Br also causes significant deviations in the photolysis behavior of Kr-HBr when compared to that of Ar-HBr, in the case of the initial excited vdW state of both clusters. A common feature of the photodissociation of Kr-HBr and Ar-HBr is the manifestation of quantum interference effects in the Kr-Br and Ar-Br rotational state distributions, in the form of pronounced structures of supernumerary rotational rainbows.     

Assessment of the PW86+PBE+XDM density functional on van der Waals complexes at non-equilibrium geometries

The deficiency of conventional density-functional theory (DFT) in properly describing van der Waals (vdW) (especially dispersion-bound) complexes has been extensively addressed in the past decade. There are now several new methods published in the literature that are capable of accurately capturing weak dispersion interactions in complexes at equilibrium geometries. However, the performance of these new methods at non-equilibrium geometries remains to be assessed. We have previously published [F. O. Kannemann and A. D. Becke, J. Chem. Theory Comput. 6, 1081 (2010); A. D. Becke, A. A. Arabi, and F. O. Kannemann, Can. J. Chem. 88, 1057 (2010)] that the functional PW86+PBE+XDM for exchange + correlation + dispersion, respectively, is a highly accurate functional for general thermochemistry and vdW complexes at equilibrium geometries. Here, we show that this nonempirical, except for two parameters in the dispersion damping part, functional also performs well for vdW complexes at compressed and stretched intermonomer separations. The mean absolute relative error (MARE) is 9.4% overall for vdW complexes in the "S22×5" database incorporating compressed and stretched geometries [J. Rezac, K. E. Riley, and P. Hobza, J. Chem. Theory Comput. 7, 2427 (2011)]. Our largest MARE on the S22×5 database is 13.3% on the compressed geometry set.     

Water monomer interaction with gold nanoclusters from van der Waals density functional theory

We investigate the interaction between water molecules and gold nanoclusters Au(n) through a systematic density functional theory study within both the generalized gradient approximation and the nonlocal van der Waals (vdW) density functional theory. Both planar (n = 6-12) and three-dimensional (3D) clusters (n = 17-20) are studied. We find that applying vdW density functional theory leads to an increase in the Au-Au bond length and a decrease in the cohesive energy for all clusters studied. We classify water adsorption on nanoclusters according to the corner, edge, and surface adsorption geometries. In both corner and edge adsorptions, water molecule approaches the cluster through the O atom. For planar clusters, surface adsorption occurs in a O-up/H-down geometry with water plane oriented nearly perpendicular to the cluster. For 3D clusters, water instead favors a near-flat surface adsorption geometry with the water O atom sitting nearly atop a surface Au atom, in agreement with previous study on bulk surfaces. Including vdW interaction increases the adsorption energy for the weak surface adsorption but reduces the adsorption energy for the strong corner adsorption due to increased water-cluster bond length. By analyzing the adsorption induced charge rearrangement through Bader's charge partitioning and electron density difference and the orbital interaction through the projected density of states, we conclude that the bonding between water and gold nanocluster is determined by an interplay between electrostatic interaction and covalent interaction involving both the water lone-pair and in-plane orbitals and the gold 5d and 6s orbitals. Including vdW interaction does not change qualitatively the physical picture but does change quantitatively the adsorption structure due to the fluxionality of gold nanoclusters.     

Mixing rules for binary Lennard–Jones chains: theory and Monte Carlo simulation

Theoretically-based van der Waals one-fluid (vdW1) mixing rules are derived for Lennard–Jones (LJ) chain mixtures. The rules provide equivalent one-fluid segment parameters for LJ size (σ) and energy () parameter as well as chain length (m) based on the parameters of the individual mixture components and the component mole fractions. The mixing rules are tested by performing Monte Carlo simulations of eight different binary mixtures and the equivalent vdW1 pure fluid, each at three densities. The simulations test the effects of changing LJ size parameter, LJ energy parameter and chain length individually and together. The effects of mole fraction and density are also examined. The mixing rules are tested for accuracy in predicting compressibility factors and radial distribution functions. It is found that the vdW1 rules provide excellent agreement when size and energy parameter are varied. Good agreement is found for mixtures with different chain lengths. The discrepancy is worst at very high densities when all component parameters are varied simultaneously.     

Naringenin and quercetin – potential anti-HCV agents for NS2 protease targets

Nonstructural proteins of hepatitis C virus had drawn much attention for the scientific fraternity in drug discovery due to its important role in the disease. 3D structure of the protein was predicted using molecular modelling protocol. Docking studies of 10 medicinal plant compounds and three drugs available in the market (control) with NS2 protease were employed by using rigid docking approach of AutoDock 4.2. Among the molecules tested for docking study, naringenin and quercetin revealed minimum binding energy of ? 7.97 and ? 7.95 kcal/mol with NS2 protease. All the ligands were docked deeply within the binding pocket region of the protein. The docking study results showed that these compounds are potential inhibitors of the target; and also all these docked compounds have good inhibition constant, vdW+Hbond+desolv energy with best RMSD value.     

Ab initio characterization of the Ne-I2 van der Waals complex: intermolecular potentials and vibrational bound states

A theoretical study of the potential energy surface and bound states is performed for the ground state of the NeI(2) van der Waals (vdW) complex. The three-dimensional interaction energies are obtained from ab initio coupled-cluster, coupled-cluster single double (triple)/complete basis set, calculations using large basis sets, of quadruple- through quintuple-zeta quality, in conjunction with relativistic effective core potentials for the heavy iodine atoms. For the analytical representation of the surface two different schemes, based on fitting and interpolation surface generation techniques, are employed. The surface shows a double-minimum topology for linear and T-shaped configurations. Full variational quantum mechanical calculations are carried out using the model surfaces, and the vibrationally averaged structures and energetics for the NeI(2) isomers are determined. The accuracy of the potential energy surfaces is validated by a comparison between the present results and the corresponding experimental data available. In lieu of more experimental measurements, we also report our results/predictions on higher bound vibrational vdW levels, and the influence of the employed surface on them is discussed.     

General methodology to optimize damping functions to account for charge penetration effects in electrostatic calculations using multicentered multipolar expansions

We developed a methodology to optimize exponential damping functions to account for charge penetration effects when computing molecular electrostatic properties using the multicentered multipolar expansion method (MME). This methodology is based in the optimization of a damping parameter set using a two-step fast local fitting procedure and the ab initio (Hartree-Fock/6-31G** and 6-31G**+) electrostatic potential calculated in a set of concentric grid of points as reference. The principal aspect of the methodology is a first local fitting step which generates a focused initial guess to improve the performance of a simplex method avoiding the use of multiple runs and the choice of initial guesses. Three different strategies for the determination of optimized damping parameters were tested in the following studies: (1) investigation of the error in the calculation of the electrostatic interaction energy for five hydrogen-bonded dimers at standard and nonstandard hydrogen-bonded geometries and at nonequilibrium geometries; (2) calculation of the electrostatic molecular properties (potential and electric field) for eight small molecular systems (methanol, ammonia, water, formamide, dichloromethane, acetone, dimethyl sulfoxide, and acetonitrile) and for the 20 amino acids. Our results show that the methodology performs well not only for small molecules but also for relatively larger molecular systems. The analysis of the distinct parameter sets associated with different optimization strategies show that (i) a specific parameter set is more suitable and more general for electrostatic interaction energy calculations, with an average absolute error of 0.46 kcal/mol at hydrogen-bond geometries; (ii) a second parameter set is more suitable for electrostatic potential and electric field calculations at and outside the van der Waals (vdW) envelope, with an average error decrease >72% at the vdW surface. A more general amino acid damping parameter set was constructed from the original damping parameters derived for the small fragments and for the amino acids. This damping set is more insensitive to protein backbone and residue side-chain conformational changes and can be very useful for future docking and molecular dynamics protein simulations using ab initio based polarizable classical methods.     

Geometries,vibrational frequencies,and excitation energies of a series of fluorine‐substituted carbenes,FCX (X = H,F, Cl,Br, and I): A high‐level multireference configuration interaction study

High‐level calculations using internally contracted multireference configuration interaction including Davidson correction (icMRCI+Q) method have been carried out for the ground singlet states, the first excited states, and the lowest triplet states of a series of fluorine‐substituted carbenes FCX (X = H, F, Cl, Br, and I). Equilibrium geometries and vibrational frequencies of the three electronic states, adiabatic transition energy of the first excited singlet state, as well as the ground singlet—lowest triplet energy gap (S‐T gap) of each of FCX carbenes have been obtained. Effects of the basis set of icMRCI+Q calculation on the geometries and energies have been investigated. In addition, various corrections, including the scalar relativistic effect, spin‐orbit coupling, and core‐valence correlation, have been studied in calculating the transition energies and the S‐T gaps, especially for heavy‐atom carbenes. This results have been compared with previous calculations using a variety of methods. Our icMRCI+Q results are in very good agreement with the high‐resolution laser‐based spectroscopic results where available. Some structure and spectroscopic constants of the fluorine‐substituted carbenes which are void in the literature have been provided with consistent high‐level calculations. © 2013 Wiley Periodicals, Inc.     

Ab initio calculations, potential representation and vibrational dynamics of He2Br2 van der Waals complex

An intermolecular potential energy surface for He(2)Br(2) complex in the ground state is calculated at the levels of fourth-order (MP4) Moller-Plesset and coupled-cluster [CCSD(T)] approximations, using large-core pseudopotential for Br atoms and the aug-cc-pV5Z basis set for He. The surface is characterized by three minima and the minimum energy pathways through them. The global minimum corresponds to a linear He-Br(2)-He configuration, while the two other ones to "police-nightstick" and tetrahedral structures. The corresponding well depths are -90.39/-89.18, -81.23/-80.78 and -74.40/-74.02 cm(-1), respectively, at MP4/CCSD(T) levels of theory. It is found that results obtained by summing three-body parametrized HeBr(2) interactions and the He-He interaction are in very good accord with the corresponding MP4/CSSD(T) configuration energies of the He(2)Br(2). Variational calculations using a sum of three-body interactions are presented to study the bound states of the vdW He(2)Br(2) complex. The binding energy D(0) and the corresponding vibrationally averaged structure are determined for different isomers of the cluster and their comparison with the available experimental data is discussed.     

On geometries of stacked and H-bonded nucleic acid base pairs determined at various DFT, MP2, and CCSD(T) levels up to the CCSD(T)/complete basis set limit level

The geometries and interaction energies of stacked and hydrogen-bonded uracil dimers and a stacked adeninecdots, three dots, centeredthymine pair were studied by means of high-level quantum chemical calculations. Specifically, standard as well as counterpoise-corrected optimizations were performed at second-order Moller-Plesset (MP2) and coupled cluster level of theory with single, double, and perturbative triple excitations [CCSD(T)] levels with various basis sets up to the complete basis set limit. The results can be summarized as follows: (i) standard geometry optimization with small basis set (e.g., 6-31G(*)) provides fairly reasonable intermolecular separation; (ii) geometry optimization with extended basis sets at the MP2 level underestimates the intermolecular distances compared to the reference CCSD(T) results, whereas the MP2/cc-pVTZ counterpoise-corrected optimization agrees well with the reference geometries and, therefore, is recommended as a next step for improving MP2/cc-pVTZ geometries; (iii) the stabilization energy of stacked nucleic acids base pairs depends considerably on the method used for geometry optimization, so the use of reliable geometries, such as counterpoise-corrected MP2/cc-pVTZ ones, is recommended; (iv) the density functional theory methods fail completely in locating the energy minima for stacked structures and when the geometries from MP2 calculations are used, the resulting stabilization energies are strongly underestimated; (v) the self-consistent charges-density functional tight binding method, with inclusion of the empirical dispersion energy, accurately reproduces interaction energies and geometries of dispersion-bonded (stacked) complexes; this method can thus be recommended for prescanning the potential energy surfaces of van der Waals complexes.     

Effect of cation size and charge on the interaction between silica surfaces in 1:1, 2:1, and 3:1 aqueous electrolytes

Application of two complementary AFM measurements, force vs separation and adhesion force, reveals the combined effects of cation size and charge (valency) on the interaction between silica surfaces in three 1:1, three 2:1, and three 3:1 metal chloride aqueous solutions of different concentrations. The interaction between the silica surfaces in 1:1 and 2:1 salt solutions is fully accounted for by ion-independent van der Waals (vdW) attraction and electric double-layer repulsion modified by cation specific adsorption to the silica surfaces. The deduced ranking of mono- and divalent cation adsorption capacity (adsorbability) to silica, Mg(2+) < Ca(2+) < Na(+) < Sr(2+) < K(+) < Cs(+), follows cation bare size as well as cation solvation energy but does not correlate with hydrated ionic radius or with volume or surface ionic charge density. In the presence of 3:1 salts, the coarse phenomenology of the force between the silica surfaces as a function of salt concentration resembles that in 1:1 and 2:1 electrolytes. Nevertheless, two fundamental differences should be noticed. First, the attraction between the silica surfaces is too large to be attributed solely to vdW force, hence implying an additional attraction mechanism or gross modification of the conventional vdW attraction. Second, neutralization of the silica surfaces occurs at trivalent cation concentrations that are 3 orders of magnitude smaller than those characterizing surface neutralization by mono- and divalent cations. Consequently, when trivalent cations are added to our cation adsorbability series the correlation with bare ion size breaks down abruptly. The strong adsorbability of trivalent cations to silica contrasts straightforward expectations based on ranking of the cationic solvation energies, thus suggesting a different adsorption mechanism which is inoperative or weak for mono- and divalent cations.     

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

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

Direct ab initio dynamics studies on the hydrogen-abstraction reactions of OH radicals with HOX (X = F, Cl, and Br)

The hydrogen abstract reactions of OH radicals with HOF (R1), HOCl (R2), and HOBr (R3) have been studied systematically by a dual-level direct-dynamics method. The geometries and frequencies of all the stationary points are optimized at the MP2/6-311+G(2d, 2p) level of theory. A hydrogen-bonded complex is located at the product channel for the OH + HOBr reaction. To improve the energetics information along the minimum energy path (MEP), single-point energy calculations are carried out at the CCSD(T)/6-311++G(3df, 3pd) level of theory. Interpolated single-point energy (ISPE) method is employed to correct the energy profiles for the three reactions. It is found that neither the barrier heights (DeltaE) nor the H-O bond dissociation energies [D(H-O)] exhibit any clear-cut linear correlations with the halogen electronegative. The decrease of DeltaE and D(H-O) for the three reactions are in order of HOF > HOBr > HOCl. Rate constants for each reaction are calculated by canonical variational transition-state theory (CVT) with a small-curvature tunneling correction (SCT) within 200-2000 K. The agreement of the rate constants with available experimental values for reactions R2 and R3 at 298 K is good. Our results show that the variational effect is small while the tunneling correction has an important contribution in the calculation of rate constants in the low-temperature range. Due to the lack of the kinetic data of these reactions, the present theoretical results are expected to be useful and reasonable to estimate the dynamical properties of these reactions over a wide temperature range where no experimental value is available.     

Comparison of different mixing rules for prediction of density and residual internal energy of binary and ternary Lennard–Jones mixtures

Mixing rules are necessary when equations of state for pure fluids are used to calculate various thermodynamic properties of fluid mixtures. The well-known van der Waals one-fluid (vdW1) mixing rules are proved to be good ones and widely used in different equations of state. But vdW1 mixing rules are valid only when molecular size differences of components in a mixture are not very large. The vdW1 type density-dependent mixing rule proposed by Chen et al. [1] is superior for the prediction of pressure and vapor–liquid equilibria when components in the mixture have very different sizes. The extension of the mixing rule to chain-like molecules and heterosegment molecules was also made with good results. In this paper, the comparison of different mixing rules are carried out further for the prediction of the density and the residual internal energy for binary and ternary Lennard–Jones (LJ) mixtures with different molecular sizes and different molecular interaction energy parameters. The results show that the significant improvement for the prediction of densities is achieved with the new mixing rule [1], and that the modification of the mixing rule for the interaction energy parameter is also necessary for better prediction of the residual internal energy.     

Non-hard sphere thermodynamic perturbation theory

A non-hard sphere (HS) perturbation scheme, recently advanced by the present author, is elaborated for several technical matters, which are key mathematical details for implementation of the non-HS perturbation scheme in a coupling parameter expansion (CPE) thermodynamic perturbation framework. NVT-Monte Carlo simulation is carried out for a generalized Lennard-Jones (LJ) 2n-n potential to obtain routine thermodynamic quantities such as excess internal energy, pressure, excess chemical potential, excess Helmholtz free energy, and excess constant volume heat capacity. Then, these new simulation data, and available simulation data in literatures about a hard core attractive Yukawa fluid and a Sutherland fluid, are used to test the non-HS CPE 3rd-order thermodynamic perturbation theory (TPT) and give a comparison between the non-HS CPE 3rd-order TPT and other theoretical approaches. It is indicated that the non-HS CPE 3rd-order TPT is superior to other traditional TPT such as van der Waals/HS (vdW/HS), perturbation theory 2 (PT2)/HS, and vdW/Yukawa (vdW/Y) theory or analytical equation of state such as mean spherical approximation (MSA)-equation of state and is at least comparable to several currently the most accurate Ornstein-Zernike integral equation theories. It is discovered that three technical issues, i.e., opening up new bridge function approximation for the reference potential, choosing proper reference potential, and/or using proper thermodynamic route for calculation of f(ex-ref), chiefly decide the quality of the non-HS CPE TPT. Considering that the non-HS perturbation scheme applies for a wide variety of model fluids, and its implementation in the CPE thermodynamic perturbation framework is amenable to high-order truncation, the non-HS CPE 3rd-order or higher order TPT will be more promising once the above-mentioned three technological advances are established.     

The Cohesive Energy and Vibration Characteristics of Parallel Single-Walled Carbon Nanotubes

Based on the van der Waals (vdW) interaction between carbon atoms, the interface cohesive energy between parallel single-walled carbon nanotubes was studied using continuous mechanics theory, and the influence of the diameter of carbon nanotubes and the distance between them on the cohesive energy was analyzed. The results show that the size has little effect on the cohesive energy between carbon nanotubes when the length of carbon nanotubes is over 10 nm. At the same time, we analyzed the cohesive energy between parallel carbon nanotubes with the molecular dynamics simulation method. The results of the two methods were compared and found to be very consistent. Based on the vdW interaction between parallel carbon nanotubes, the vibration characteristics of the two parallel carbon nanotube system were analyzed based on the continuous mechanical Euler-beam model. The effects of the vdW force between carbon nanotubes, the diameter and length of carbon nanotubes on the vibration frequency of carbon nanotubes was studied. The obtained results are helpful in improving the understanding of the vibration characteristics of carbon nanotubes and provide an important theoretical basis for their application.     

Conformations and energetics of sulfur and selenium diimides

The geometries and energetics of different conformations of sulfur and selenium diimides E(NR)(2) (E = S, Se; R = H, Me, (t)Bu, C(6)H(3)Me(2)-2,6, SiMe(3)) have been studied by using various ab initio and DFT molecular orbital techniques. The syn,syn conformation is found to be most stable for parent E(NH)(2), but in general, the preferred molecular conformation for substituted chalcogen diimides is syn,anti. In the case of E(NH)(2) the present calculations further confirm that syn,syn and syn,anti conformations lie energetically close to each other. From the three different theoretical methods used, B3PW91/6-31G proved to be the most suitable method for predicting the geometries of chalcogen diimides. The optimized geometrical parameters are in a good agreement with all available experimental data. While qualitative energy ordering of the different conformations is independent of the level of theory, the quantitative energy differences are dependent on the method used. The performance and reliability of higher level ab initio calculations and DFT methods using large basis sets were tested and compared with experimental information where available. All of the higher level ab inito methods give very similar results, but the use of large basis sets with the B3PW91 method does not increase the reliability of the results. The combination of CCSD(T)/cc-pVDZ with the B3PW91/6-31G-optimized geometries is found to be the method of choice to study energetic properties of chalcogen diimides.     

Assessment of density functional theory to calculate the phase transition pressure of ice

To assess the accuracy of density functional theory (DFT) methods in describing hydrogen bonding in condensed phases, we benchmarked their performance in describing phase transitions among different phases of ice. We performed DFT calculations of ice for phases Ih, II, III, VI and VII using BLYP, PW91, PBE, PBE-D, PBEsol, B3LYP, PBE0, and PBE0-D, and compared the calculated phase transition pressures between Ih-II, Ih-III, II-VI, and VI-VII with the 0 K experimental values of Whalley [J. Chem. Phys., 1984, 81, 4087]. From the geometry optimization of many different candidates, we found that the most stable proton orientation as well as the phase transition pressure does not show much functional dependence for the generalized gradient approximation and hybrid functionals. Although all these methods overestimated the phase transition pressure, the addition of van der Waals (vdW) correction using PBE-D and PBE0-D reduced the transition pressure and improved the agreement for Ih-II. On the other hand, energy ordering between VI and VII reversed and gave an unphysical negative transition pressure. Binding energy profiles of a few conformations of water dimers were calculated to understand the improvement for certain transitions and failures for others with the vdW correction. We conclude that vdW dispersion forces must be considered to accurately describe the hydrogen bond in many different phases of ice, but the simple addition of the R(-6) term with a small basis set tends to over stabilize certain geometries giving unphysical ordering in the high density phases.     

Power series expansion of the random phase approximation correlation energy: The role of the third- and higher-order contributions

We derive a power expansion of the correlation energy of weakly bound systems within the random phase approximation (RPA), in terms of the Coulomb interaction operator, and we show that the asymptotic limit of the second- and third-order terms yields the van der Waals (vdW) dispersion energy terms derived by Zaremba-Kohn and Axilrod-Teller within perturbation theory. We then show that the use of the second-order expansion of the RPA correlation energy results in rather inaccurate binding energy curves for weakly bonded systems, and discuss the implications of our findings for the development of approximate vdW density functionals. We also assess the accuracy of different exchange energy functionals used in the derivation of vdW density functionals.