Towards extending the applicability of density functional theory to weakly bound systems

Abstract

While the attempts currently in progress in several groups for the rigorous inclusion of dispersion interactions in density functional theory (DFT) calculations mature and evolve into practical methodology, we contribute to the debate on the applicability of current functionals to the calculation of weak interaction with a systematic investigation of a few, typical, weakly bound systems. We have used both pure DFT and a hybrid approach in which the total interaction energy is partitioned into two parts: (a) the dispersion energy which, in a first approximation is the contribution due to intermonomer correlations and (b) all other interactions. The first component is accurately obtained at all distances of interest by means of a well-known damped multipolar expansion of the dispersion energy while for the second component different approximations will be evaluated. The need to avoid double counting a fraction of the correlation energy when using the hybrid approach and the choice of the appropriate functional are also discussed. We consider four systems of increasing binding strength, namely the Ar2 and Kr2 dimers, the benzene dimer, the water dimer, and a few metal carbonyls. For pure DFT calculations we confirm the conclusion reached by others concerning (a) the strong dependence of the results on the choice of the GGA functional for dispersion-dominated interaction (noble gases and benzene) with the overall tendency to yield underbinding and (b) the relatively accurate, functional-independent, description for that DFT gives of water, which we attribute to the fact that this system is dominated by electrostatic interactions. For the carbonyls we find that DFT yields results which area again strongly dependent on the choice of the functional and show a tendency to give overbinding. Our hybrid method shows instead shortcomings only for the noble gases. The problem in this case is traceable to the well-known difficulties that all current functionals experience at medium–large intermonomer separations. The quality of the hybrid results improves markedly for benzene due to the large value of both dispersion and repulsive interactions at the equilibrium distance for this dimer, which makes the balance between the two, less delicate. Excellent results are also obtained for water (for the same reason as indicated above) and more significantly for the carbonyls where we find that dispersion contributes to the binding more than it could be guessed a priori. We do not claim to have found a general solution to this difficult problem, but we aim at providing a quantitative assessment to where the problems are pointing at directions from which a general solution may, eventually, emerge.

Comment on “Noncovalent functionalization of carbon nanotubes by aromatic organic molecules” [Appl. Phys. Lett. 82, 3746 (2003)]

“…5 The results in Table I of Ref. 1 for both benzene and cyclohexane adsorbed on nanotubes are perfectly consistent with physisorption, as described by GGA: 6 large equilibrium molecule-nanotube distance, very small binding energy, no sizable charge transfer. In spite of such evidence, Fig.…”

Long-range corrected hybrid density functionals with damped atom–atom dispersion corrections

“…Following the general form of the DFT-D scheme [30][31][32][33][34][35][36][37][38][39][40][41][42][43], our total energy…”

“…Such an approach was used long ago to correct HF calculations [29], and over the past 8 years has been incorporated into density functional theory [30][31][32][33][34][35][36][37][38][39][40][41][42][43]. These DFT-D (density functional theory with empirical dispersion corrections) schemes have shown generally satisfactory performance on a large set of non-covalent systems [35,40].…”

Adenine versus guanine quartets in aqueous solution: dispersion-corrected DFT study on the differences in π-stacking and hydrogen-bonding behavior

“…Calculations were done using dispersion-corrected DFT-D as developed by Grimme [41][42][43][44][45][46][47] for a correct treatment of the stacking interactions between the DNA bases. The density functionals are augmented with an empirical correction for long-range dispersion effects, described by a sum of damped interatomic potentials of the form C 6 R -6 added to the usual DFT energy [41][42][43][44][45][46][47].…”