1Exact exchange and Wilson–Levy correlation: a pragmatic device for studying complex weakly-bonded systems
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj1
1
Exact exchange and Wilson–Levy correlation: a pragmatic device for studying complex weakly-bonded systems
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj1
2
The Wilson–Levy (WL) correlation functional is used together with Hartree–Fock (HF) theory to evaluate interaction energies at intermediate separations (i.e. around equilibrium separation) for several weakly-bonded systems.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj2
3
The HF+WL approach reproduces binding trends for all complexes studied: selected rare-gas dimers, isomers of the methane dimer, benzene dimer and naphthalene dimer, and base-pair stacking structures for pyrimidine, cytosine, uracil and guanine dimers.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res1
4
These HF+WL data are contrasted against results obtained from some popular functionals (including B3LYP and PBE), as well as two newly-developed functionals, X3LYP and xPBE.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj3
5
The utility of HF+WL, with reference to exact-exchange (EXX) density-functional theory, is discussed in terms of a suggested EXXWL exchange–correlation functional.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con1
Introduction
6
The interaction energy between two weakly-interacting fragments comprises many distinct and well-understood (although perhaps not easy to calculate) contributions.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac1
7
Following Stone,1 the interaction energy can be decomposed into first-order, second-order and higher-order contributions.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac1
8
Hartree–Fock (HF) theory adequately captures (at small and large separations) first-order contributions such as the electrostatic interaction and the exchange–repulsion interaction.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met1
9
The performance of HF theory for second-order contributions is less comprehensive.
Type: Method |
Advantage: No |
Novelty: Old |
ConceptID: Met1
10
Although HF theory can yield a reasonable estimate of the induction contribution, HF theory cannot recover the dispersion energy.
Type: Method |
Advantage: No |
Novelty: Old |
ConceptID: Met1
11
Higher-order terms are recovered (with varying degrees of accuracy) via use of post-HF methods such as Møller–Plesset second-order perturbation theory, MP2, and coupled-cluster singles and doubles with triples corrections, CCSD(T).
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met2
12
Both induction and dispersion are present in the regime of large fragment separations (the ‘long-range’ form) where the charge distributions on each fragment do not overlap, and are also present in the small to intermediate separation regime (the ‘short-range’ form) where the charge distributions do overlap.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met2
13
Where charge distributions overlap, second-order contributions are modified by general overlap effects as well as electron exchange, leading to damped induction, exchange-induction, damped dispersion and exchange–dispersion contributions.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met2
14
These short-range terms should behave as a function of overlap between fragments.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met2
15
Therefore, in this article, the general phrase ‘dispersion energy’ does not solely refer to the London dispersion (with r−6 dependence) that characterises the long-range contribution, rather, it refers to the dispersion interaction at all separations.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod1
16
The reader is directed to the body of work concerned with symmetry-adapted perturbation theory (SAPT), where such contributions to the interaction energy are rigorously defined (e.g., the review article of Jeziorski et al2)..
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac2
17
Accurate calculation of these second-order contributions, at any separation, for large and/or complex weakly-interacting systems remains a challenging prospect in theoretical chemistry.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac3
18
In the worst case, the use of post-HF methods such as MP2 or CCSD(T) can be limited at best to rather modest basis set sizes.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met3
19
These methods may yield unreliable results for such small basis sets.
Type: Method |
Advantage: No |
Novelty: Old |
ConceptID: Met3
20
For example, even with use of large basis sets, MP2 can substantially overestimate the attraction in weakly-bound aromatic systems.3
Type: Method |
Advantage: No |
Novelty: Old |
ConceptID: Met3
21
However, recent progress in the development of algorithms for the rapid evaluation of MP2 energies and gradients (e.g. see Saebø et al4)., should in future enable practical geometry optimisations using larger basis sets.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met4
22
Second-order contributions to the interaction energy arise from electronic correlation effects (beyond electron exchange alone).
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac4
23
Density-functional theory (DFT) is often proposed as a means of obtaining the correlation contribution to the total electronic energy, in a way that is computationally economical.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met5
24
However, it is now well-known that no mainstream, practical correlation functional can recover the r−6 dependence of London dispersion at large fragment separations (see for example refs. 7 and 25, and references therein).
Type: Method |
Advantage: No |
Novelty: Old |
ConceptID: Met5
25
This fact will not be disputed here; lack of London dispersion remains a major limitation of all approaches that involve Kohn–Sham DFT in some way.
Type: Method |
Advantage: No |
Novelty: Old |
ConceptID: Met5
26
However, many studies have amply demonstrated that application of DFT to weakly-bound systems can not yield globally reliable estimates of the interaction energy even at intermediate intermonomer separations5–9 where there is overlap between fragments.
Type: Method |
Advantage: No |
Novelty: Old |
ConceptID: Met5
27
Several groups have made progress in the development of ‘seamless DFT’ approaches that can describe dispersion contributions at all separations10–12 (including in the overlapping regime).
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met6
28
However, all of these approaches are at present too computationally intensive to be applied to large and/or complex weakly-bonded systems.
Type: Method |
Advantage: No |
Novelty: Old |
ConceptID: Met6
29
On the basis of Heitler–London calculations, Cybulski and Seversen13 correctly identified that the correlation functionals used in their study only showed binding when it was due to exchange effects from the correlation functional, such that the total intermolecular exchange was negative (for their decomposition scheme).
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac5
30
Furthermore, they noted that, as expected, this intermolecular exchange contribution vanishes faster with respect to separation than r−6.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac5
31
This overlap effect has been observed (but not quantified) by others, e.g. see .ref. 7
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac5
32
The decomposition analysis used by Cybulski and Seversen, while helpfully demonstrating the overlap effect of binding for some correlation functionals, does not explain why these correlation functionals yield negative contributions to the intermolecular exchange.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac5
33
While it is often the correlation functional which is identified as deficient in recovering second-order contributions in the overlapping regime, several studies have instead identified the exchange functional as problematic,6,8,13–16 especially under the conditions of low-density/high-density-gradient.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac5
34
An accurate description of exchange under such conditions is vital for describing weakly-bonded systems.
Type: Motivation |
Advantage: None |
Novelty: None |
ConceptID: Mot1
35
Empirical scaling of DFT exchange has also been examined by several groups,14,16,17 with a view to correcting this problem.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met7
36
However, the results of these studies do not show transferability of the scaling factor across different systems.
Type: Method |
Advantage: No |
Novelty: Old |
ConceptID: Met7
37
Therefore, the fact that mainstream DFT cannot recover short-range second-order contributions must be due, in part, to problems with the exchange functional.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac6
38
Despite these well-known deficiencies of DFT, a number of schemes for empirically correcting DFT calculations by addition of damped dispersion (DFT+D) have been proposed.6,18–20
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac7
39
As an alternative to post-HF theories and also DFT+D, another route to gaining economical (and reasonable) estimates of the interaction energy for large weakly-bound systems is to use a Hartree–Fock (HF) plus dispersion approach (HF+D).
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met8
40
The use of HF implies that, for a given density, the exchange contribution will be exact (by definition), including under low-density/high-density-gradient conditions.
Type: Method |
Advantage: Yes |
Novelty: Old |
ConceptID: Met8
Type: Method |
Advantage: Yes |
Novelty: Old |
ConceptID: Met8
42
One issue to note with both HF+D and DFT+D is that these methods require good quality dispersion coefficients.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met8
43
Also, the form of the damping function, including anisotropic damping, can be a very important factor in gaining good results from these DFT+D models.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met8
44
While some studies have reported progress in this area,18,19 further studies into damping are probably warranted.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac8
45
Both of these ‘+dispersion’ approaches may succeed for the long-range part of the dispersion interaction, but it is not entirely clear that the DFT+D approach is correct in the overlapping regime.25
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac8
46
The problem of calculating the short-range interaction energy can also be approached via the addition of a correlation contribution derived from DFT onto the total HF energy of the system, bearing in mind the deficiencies pointed out by Cybulski and Seversen.13
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met9
47
As in the HF+D approach, the exchange contribution should behave correctly at long range.
Type: Method |
Advantage: Yes |
Novelty: Old |
ConceptID: Met9
Type: Method |
Advantage: Yes |
Novelty: Old |
ConceptID: Met9
49
Some other correlation functionals also showed binding, but did not recover qualitatively correct equilibrium binding energies.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac9
50
Moreover, this study also showed that in all cases where binding was found, the interaction energy curve decayed to zero faster than r−6, clearly not capturing any long-range dispersion contribution (as expected).
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac9
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac9
52
This HF+DFT approach is not the same as Hartree–Fock Kohn–Sham (HFKS) theory, but Pérez-Jordá et al26. reported that the HF+DFT results for rare-gas dimers were almost identical to the HFKS results.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac10
53
In particular, these authors demonstrated that the Wilson–Levy correlation functional (WL)27 worked extremely well in partnership with HF theory, in reproducing the experimental well-depths of the rare-gas dimer systems studied “with remarkable agreement”.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac10
54
Their data, shown in Fig. 1, emphatically reinforces their statement.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac10
55
This approach will herein be referred to as HF+WL.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac10
56
The WL correlation functional is based on the Wigner correlation functional30 (W38), contains gradient corrections, and includes four empirically-determined parameters.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac10
57
The parameters were determined via comparison with the exact (calculated) correlation energy (in the traditional quantum chemistry definition) for several atoms (neutral and cationic), plus the experimental correlation energy of the helium atom.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac10
58
Hartree–Fock densities of these atoms were used in the parametrisation of this functional.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac10
59
As noted by Wilson and Levy,27 since the traditional quantum chemistry definition of correlation energy is not always equivalent to the DFT definition, the WL functional may not be successful if applied to systems where the HF density is substantially different to the true density of the system.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac10
60
Since the WL correlation functional was designed for use with HF densities, it may be an appropriate approximate correlation functional to use with exact exchange methods, at least for studying pure van der Waals interactions.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac10
61
However, the utility of this functional has principally been demonstrated only as a post-HF correction8,31,32 for atoms and molecules, not as a dynamic correlation contribution that is self-consistently evaluated.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac10
62
The performance of WL correlation for covalent bonds has not been as extensively tested.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac11
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac11
64
Their data indicate that, compared with other correlation functionals used with B88, WL underestimates covalent bond lengths.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac11
65
The HF+WL method should essentially be regarded as an empirical correction to the HF interaction energy, and presumably has inherent deficiencies with the description of intermolecular exchange.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac11
66
It is the applicability of this approach as a pragmatic device that will be further explored in this work.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj4
67
Here, the HF+WL approach is applied to all homoatomic and heteroatomic rare-gas dimers RgKr {Rg = He, Ne, Ar, Kr}, the methane dimer, benzene dimer, naphthalene dimer and also DNA base-pair stacking interactions.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj4
68
Some comments on the choice of example systems studied here are warranted.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met10
69
It is noted that, in general, any perturbative correction to the HF interaction energy should only hold under conditions such that the HF wavefunction will not substantially distort if (some) second-order effects are dynamically included.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod2
70
Therefore, the method has been tested on systems where dynamic inclusion of such effects would not substantially distort the wavefunction, compared with the HF wavefunction.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod2
71
Systems exhibiting pure van der Waals interactions (dispersion and repulsion) are ideal in this case, and the systems chosen here reflect this.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod2
72
In general, all HF + dispersion methods are subject to this condition, which does limit the applicability of this type of approach.
Type: Method |
Advantage: No |
Novelty: New |
ConceptID: Met10
73
However, it appears possible that the HF+WL method can be extended such that all types of non-bonded system could be treated in future.
Type: Hypothesis |
Advantage: None |
Novelty: None |
ConceptID: Hyp1
74
Specifically, it is conjectured that use of the WL correlation functional, applied self-consistently in conjunction with exact-exchange DFT (e.g. the optimised effective potential, OEP35), denoted here as the EXXWL functional, should have a performance similar to the HF+WL approach.
Type: Hypothesis |
Advantage: None |
Novelty: None |
ConceptID: Hyp1
75
The EXXWL functional would open up possibilities for qualitatively accurate calculations on complex systems not yet accessible by other methodologies, e.g. physisorption onto nanotubes, gas–surface interactions, etc.
Type: Method |
Advantage: Yes |
Novelty: New |
ConceptID: Met10
76
The example systems selected in this work examine aromatic ring interactions.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj5
77
Such interactions are ubiquitous in Nature, and should be described properly by any theory that claims to model ‘biological systems’.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac12
78
Finally, it must be emphasised that, as demonstrated by Pérez-Jordá et al.,26 the HF+WL method has no pretensions to capturing London dispersion (the long-range contribution).
Type: Method |
Advantage: No |
Novelty: New |
ConceptID: Met10
79
However, this caveat does not necessarily exclude the HF+WL method yielding binding when there is overlap between fragments, crucially, around the equilibrium separation for many weakly-interacting systems.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met10
80
Having accepted the caveats of HF+WL, the utility of this approach as a pragmatic device for locating low-energy structures of large van der Waals systems (some of which have low symmetry) in a simple and quick way will be demonstrated.
Type: Goal |
Advantage: None |
Novelty: None |
ConceptID: Goa1
Type: Goal |
Advantage: None |
Novelty: None |
ConceptID: Goa2
82
It has been recently suggested that both of these functionals should adequately describe non-bonded interactions for systems involving the first ten elements in the periodic table36,37.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac13
Computational methods
83
All calculations in this work were done using GAUSSIAN.9838
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp1
84
The monomers were treated as rigid bodies throughout.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod3
85
Unless otherwise stated, the monomer structures were fully optimised at the MP2/6-31G** level of theory.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod3
86
For a given arrangement of rigid monomers, the HF wavefunction was converged and the total electronic energy (EHF) obtained, for both the complex, and also for each monomer in the complex.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp2
87
The monomer wavefunctions and energies were obtained in the full set of dimer basis functions, such that the interaction energies were corrected for basis-set superposition error (BSSE) via the counterpoise correction39 (CP).
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp2
88
The correlation energy was then evaluated using the corresponding HF density.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp2
89
The integration was done over a pruned grid with 99 radial shells and 590 angular points per shell.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp2
90
The quadrature used here ensured that the correlation energy was converged to within 1–2 μEh per atom (depending on atom type).
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp2
91
The WL correlation energy (Ecorr) was then extracted from the calculation for both the complex and individual monomers (again using the CP correction).
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp2
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp3
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp3
94
The basis set dependence of the HF+WL energies was examined.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj6
95
This was done to confirm the dependence of the interaction energy on overlap, for short- to medium-range separations.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj6
96
The interaction energy as calculated here by the supermolecule method is of course not variational, so we do not necessarily expect to see convergence of this energy from above.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj6
97
Furthermore, since the HF+WL method used here is an ad hoc post-Hartree–Fock correction (and not variationally obtained), this will cause further problems when seeking monotonically decreasing behaviour of interaction energies as a function of basis-set size.
Type: Object |
Advantage: No |
Novelty: New |
ConceptID: Obj6
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod4
99
The HF+WL interaction energy was calculated as a function of separation in all cases, and the minimum interaction energy was reported in each case.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp4
100
For the methane dimer, three different relative orientations were examined, as shown in Fig. 2.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod5
101
For the benzene dimer, three geometries were considered: a T-shaped structure (C2v symmetry), a slipped-parallel (SP) structure (C2h symmetry) and an undisplaced parallel (F) structure (D6h symmetry), also shown in Fig. 2.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod6
102
For the naphthalene dimer, four low-energy structures and one high-energy structure were examined: two crossed structures (D2d and C2 symmetries) and two slipped-parallel structures (C2h and Ci symmetries), and one undisplaced parallel structure (D2h).
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod7
103
The four low-energy structures are shown in Fig. 3.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod7
104
Optimised naphthalene dimer geometries were taken from previous work.44
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp5
105
In this case, the HF+WL interaction energies were calculated as a function of separation (along the vector that connects centre of mass points), while the relative orientations of each dimer structure were kept fixed.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp5
106
Finally, for the base-pair stacking dimers, the geometries were taken from Hobza and Šponer.45
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod8
107
That is, antiparallel undisplaced structures were used for the pyrimidine and cytosine dimers, and the vertical displacement was optimised using HF+WL.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod8
108
Compared with pyrimidine and cytosine, the procedure for creating the input structures for uracil and guanine was a little different.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod8
109
Following Hobza and Šponer, the uracil dimer geometry was optimised in this work at the MP2/6-31G** level, and agrees closely with the structural data of Hobza and Šponer.46
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod8
110
Using this MP2 structure, the separation (along the vector connecting the centre of mass of each molecule) was optimised using HF+WL, while keeping the angle of the monomer normal vectors (normal to the ring plane) constant.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod8
111
The guanine dimer structure was an undisplaced stacked geometry, where one monomer was rotated clockwise about its centre of mass (in the plane of the molecule).
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod8
112
The optimum twist angle was identified at the HF+WL/aug-cc-pVDZ level, for a fixed vertical separation of 3.4 Å.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res2
113
The vertical separation (for the fixed optimum twist angle) was then optimised for this structure using HF+WL.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp6
114
The geometries used for the base-pair calculations can be found in Fig. 4.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod8
Results and discussion
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs1
116
The interaction energies are also plotted in Fig. 1.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs2
117
In agreement with the previous findings of Pérez-Jordá et al.,26 the trend of the well depths presented here are in very good qualitative agreement with the experimental results, while not being quantitatively accurate.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res3
118
The minimum energy separation is consistently underestimated by not more than 0.5 a0, as is also demonstrated in the results of Pérez-Jordá et al. The HF+WL approach has not captured the slight contraction in separation in going from KrHe to KrNe.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res4
119
Examination of these data plus the data of Pérez-Jordá et al26. show that, with the exception of dimers containing helium, the HF+WL underestimates the interaction energy.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res5
120
The interaction energy for each system was also calculated for separations beyond the minimum-energy point.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp7
121
These calculations confirmed what has been reported previously by Pérez-Jordá et al.,26 namely that the interaction energy decays to zero too rapidly compared with the r−6 decay observed with long-range London dispersion interactions.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con2
122
These results re-emphasise that weak interactions for non-overlapping fragments cannot be captured by the HF+WL approach.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con2
123
However, these data also suggest that failure to recover long-range London dispersion does not preclude a DFT-based methodology from recovering binding at distances where fragments do overlap.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con3
124
As shown in Fig. 1, the qualitative agreement found using HF+WL clearly exceeds the performance of any other mainstream DFT approach (demonstrated, for example, in the results of Patton and Pederson.9)
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con4
125
Fig. 1 also re-emphasises the point made by others that the PW91 and PBE functionals do not show great variation in energy between each rare-gas dimer system.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con4
126
The PW91 data fortuitously coincides with the HeKr and NeKr points, while the PBE data fortuitously coincides with the ArKr and Kr2 points.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res6
127
Neither can be said to describe rare-gas equilibrium interaction energies satisfactorily, while the HF+WL equilibrium interaction energies follow the trend of the Ogilvie and Wang data relatively closely.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con5
128
The interaction energies of these rare-gas systems, calculated using HF+WL for a range of basis sets aug-cc-pVXZ {X = D, T, Q, 5} is shown in Fig. 5.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs3
129
For each basis the separation between atoms was optimised, but in most cases did not change appreciably for the different basis sets.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res7
130
It seems clear that the energy has converged by the aug-cc-pVTZ basis, and that the aug-cc-pVDZ basis still follows a reasonable qualitative trend (especially when compared with the PBE and PW91 data), although the trend for HeKr and NeKr is not captured.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con6
131
The basis set dependence was also checked for the smaller 6-31G* basis.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp8
132
For all systems but KrNe, the smaller basis set performed no worse than the aug-cc-pV5Z basis, in terms of capturing the correct trend.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res8
133
However, the disagreement for HeNe is striking – the 6-31G* basis recovering an interaction energy of only −5.5 K. The optimum seperation was increased (to 7.0 a0).
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res8
134
As will be confirmed with some of the other systems, the overlap obtained using the 6-31G* basis appears insufficient to ensure confident use of HF+WL here.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con7
135
Clearly the HF+WL approach is not perfect here – but it is outperforming other DFT-based approaches.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con7
136
Three structures of the methane dimer were studied, as shown in Fig. 2; one D3d structure and two D3h structures.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj7
137
The D3d structure has been identified by Tsuzuki and coworkers49 as the global minimum structure, at least at the MP3/aug(df,pd)-6-311G** level of theory.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj7
138
The D3d structure corresponds with isomer H, and the two D3h structures as isomers A and G, in the nomenclature of Tsuzuki et al.49
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj7
139
Isomers A and G are other low-energy structures, and therefore these two geometries probe the interaction energy in regions of configuration space away from the global minimum.
Type: Object |
Advantage: Yes |
Novelty: New |
ConceptID: Obj7
140
Separations and interaction energies for these three structures, calculated using the HF+WL method, as well as the X3LYP and xPBE functionals, are shown in Table 2, along with the ab initio data of Tsuzuki et al.49
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs4
141
The aug-cc-pVTZ basis was used here for HF+WL, X3LYP and xPBE.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp9
142
These data demonstrate that the HF+WL approach captures the binding trends for these three structures, for all basis sets considered, while clearly not displaying quantitative accuracy.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con8
143
Again, the HF+WL approach underestimates the minimum-energy separation, found to be 4.6, 3.5 and 3.4 Å for the A, G and H structures, respectively.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res9
144
For comparison, the MP3 separations from Tsuzuki et. al. are 4.8, 3.8 and 3.8 Å for the A, G and H structures, respectively.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac14
145
In contrast, the X3LYP and xPBE data seem less reliable here.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con9
146
Although the xPBE data show binding for each structure, and also that the trend of the interaction energies seems correct (A > G > H), the binding is reduced compared with the MP3 results (except in the case of isomer A).
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res10
147
The range of xPBE interaction energies is significantly compressedwith respect to the MP3 data.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res10
148
The X3LYP data show binding only for structure A. Out to a separation of 6 Å, the other two structures were not bound using the X3LYP functional.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res11
149
Also in contrast with the rare-gas systems, the isomers studied here showed little sensitivity to changes in basis set using the HF+WL approach.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res12
150
For each basis set used (6-31G*, 6-31G*(0.25), aug-cc-pVXZ {X = D, T, Q}) the trend of interaction energies (A > G > H) was always followed.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res13
151
The minimum energy separations and interaction energies for the three benzene dimer structures, calculated at the HF+WL/aug-cc-pVDZ level, are shown in Tables 3 and 4, respectively.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs5
152
Again, the D6h (F) structure provides a probe of the HF+WL method away from the two low-energy structures.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res14
153
For comparison, the data of Tsuzuki et al50. are shown.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac15
154
These data were obtained using their extrapolation method where the interaction energy is calculated at the MP2 basis set limit plus a CCSD(T) correction evaluated with a large basis set.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac15
155
Again, it can be seen that the HF+WL separations are underestimated compared with the best estimate data.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res15
156
The interaction energies, although overestimated for the T and SP structures, still show a reasonable trend compared with the F structure, although the T and SP structures are not isoenergetic as suggested in the results of the ab initio calculations.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res16
157
The F structure is bound in the HF+WL model, but underestimates the attraction, again emphasising that the HF+WL method cannot compete with post-HF approaches.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res17
158
Also shown in Table 4 for comparison are results from full DFT calculations, as reported by Tsuzuki and Lüthi51 and also from current work, that confirm the unsuitability of these functionals as applied to weakly-bound systems.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res18
159
The HF+WL approach provides a reasonable estimate of the energetic and structural trends of this system (while not being quantitatively accurate), where full DFT can not.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con10
160
The similarity between the PBE and PW91 results is also noted.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res19
161
Table 4 shows the interaction energy of the T, F and SP structures calculated at the X3LYP/aug-cc-pVDZ and xPBE/aug-cc-pVDZ levels of theory, at the geometries taken from Tsuzuki et al.,50 revealing a similar performance of these two functionals compared with existing functionals.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res20
162
The interaction energy as a function of centre-of-mass separation was also investigated for these two functionals.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj8
163
The F structure was unbound at all separations up to and including 6.0 Å at both the xPBE/aug-cc-pVDZ and X3LYP/aug-cc-pVDZ levels.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs6
164
The T structure was bound at a separation of 5.7 Å at the xPBE/aug-cc-pVDZ level, with an optimum interaction energy of −0.48 kcal mol−1.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res21
165
Similarly, at the X3LYP/aug-cc-pVDZ level, the T structure was bound at an optimum distance of 5.6 Å with an interaction energy of −0.37 kcal mol−1.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res22
166
The SP structure was investigated in the neighbourhood of the structure of Tsuzuki et al., such that the vertical and horizontal separations, as well as the distance along the vector connecting the two centre of mass points, was varied by displacements of ±0.2 Å.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj9
167
None of these locally-perturbed structures showed binding at the xPBE/aug-cc-pVDZ level.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res23
168
The performance was very similar at the X3LYP/aug-cc-pVDZ level.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res24
169
The only structure that yielded an interaction energy of less than +1 (but greater than zero) kcal mol−1, for either xPBE or X3LYP, was for the case when the distance along the vector connecting the two centre of mass points was extended by 0.2 Å.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs7
170
These results, as well as those summarised in Table 4, demonstrate that the recovery of binding for molecular complexes, even at separations where the fragments overlap, remains an elusive goal for mainstream DFT exchange–correlation functionals.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con11
171
The performance of all of the functionals examined here is poor in comparison with the HF+WL approach.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con12
172
The basis set dependence of the the HF+WL approach was also gauged for the benzene dimer structures.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj10
173
At the minimum-energy separations found at the HF+WL/aug-cc-pVDZ level, the interaction energies were recalculated using the 6-31G* basis set.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp10
174
While the T and SP structures remained bound (interaction energies of −2.83 kcal mol−1 and −1.85 kcal mol−1, respectively), the F structure is unbound at this level of theory, at +0.78 kcal mol−1.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs8
175
Exploration of the interaction energy of the F structure at a range of separations confirmed that this structure remained unbound.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs9
176
Increasing flexibility by using a triple split valence basis, in this case using 6-311G*, leads to binding of the F structure (−0.53 kcal mol−1).
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs10
177
At the HF+WL/6-311G* level, the T and SP structures also remained with reasonable interaction energies (−3.16 and −2.84 kcal mol−1, respectively).
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs11
178
Similarly, diffuse functions were added to the 6-31G* basis, such that the 6-31+G* basis set was considered.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp11
179
At this level, all structures are bound, with F, T and SP structures having interaction energies −0.20, −3.08 and −2.58 kcal mol−1, respectively.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs12
180
These data suggest that a basis with at least split valence plus diffuse and polarization functions, or triple split valence plus polarization should be used to obtain qualitatively correct HF+WL interaction energies.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con13
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac16
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac16
183
The naphthalene dimer presents a significant challenge in the calculation of competitive low-energy structures.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj11
184
Calculations at the MP2 level tend to overestimate the attraction in the naphthalene dimer, most significantly for the structures lowest in energy.52
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac17
185
Application of DFT to this system is also particularly unsuitable: the high-energy structures have attractive electrostatic energy contributions and relatively small dispersion contributions, and the low-energy structures have electrostatic contributions that are less favourable, while having relatively large dispersion contributions.44,52
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac17
186
This is opposed to the benzene dimer where the low-energy structures have favourable electrostatic interactions.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac17
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac17
188
In other words, for the naphthalene dimer, fortuitous location of the low-energy structures, by virtue of attractive electrostatic interactions, is not possible, making this system an excellent test case for new functionals.
Type: Motivation |
Advantage: None |
Novelty: None |
ConceptID: Mot2
189
Most published theoretical studies of the naphthalene dimer have made use of MP2 calculations.44,52,54–56
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac18
190
Full CCSD(T) calculations with any reasonable basis set are currently out of reach, at least for the low symmetry structures.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac18
191
Several studies have reported low-energy structures of the naphthalene dimer using MP2 calculations with a range of basis sets,44,52,54,55 including correlation corrections to MP2 interaction energies, using CCSD(T) calculations.52
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac19
192
It appears that, at present, the lowest-energy structure corresponds to a crossed C2 geometry, with a slipped Ci structure also energetically competitive.44
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac19
193
This result has also been reported by Gonzalez and Lim, using a HF+D approach.24
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac19
194
Here, four low-energy structures (Fig. 3) and one high-energy structure are considered.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj12
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs13
196
For each low-energy structure, the separations were varied with displacements of ±0.2 Å away (along the vector connecting the centre ofmass points) from the optimum separation (as was determined by MP2/6-31G* geometry optimizations.44)
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp12
197
For the high-energy D2h structure, the optimum separation was taken from Tsuzuki et al.52
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp12
198
This separation was varied along the vector connecting the centre of mass points in the same way.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp12
199
It is seen in Table 5 that the HF+WL separations all seem reasonable when compared with the ab initio data, although are consistently too short.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res25
200
The trend in energetics for the D2d, Ci and C2h structures agrees with the best estimates, although there is not quantitative agreement.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res25
201
In particular, the best estimates indicate that the D2d and C2h structures are almost isoenergetic, while the HF+WL data support a difference of over 0.5 kcal mol−1.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res25
202
Tsuzuki et al. did not examine the C2 structure, so we have made comparisons with the trend from the HF+D results of Gonzalez and Lim.24
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac20
203
In agreement with their findings, and with previous MP2 optimisations,44 the HF+WL approach also supports the C2 structure as the putative global minimum for the naphthalene dimer system.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con14
204
The absolute HF+WL interaction energy for the D2h structure deviates significantly from the ab initio value, although it is bound, and exhibits the smallest interaction energy of the set, in accord with Tsuzuki et al. Again, these data emphasise that HF+WL cannot outperform high-quality post-HF calculations.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con14
205
The X3LYP and xPBE functionals were also tested for all five naphthalene dimer structures.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj13
206
For both of these functionals, none of the five structures were bound, out to a centre-of-mass separation of 6.0 Å.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs14
207
For the xPBE functional, at a 6.0 Å separation, interaction energies were 0.19, 0.12, 0.72, 0.06 and 0.23 kcal mol−1 for the D2d, C2, C2h, Ci and D2h structures, respectively.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs15
208
The X3LYP values showed a similar trend, with absolute interaction energies slightly higher in all cases.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res26
209
This behaviour is similar to that noted for other popular functionals.53
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res26
210
It is clear that neither of these functionals can reliably yield binding for aromatic ring interactions at a qualitative level.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con15
211
The basis set dependence of the HF+WL approach was also tested for the naphthalene dimer.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj14
212
At first the 6-31G* basis was used.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp13
213
This appeared to give satisfactory results for the four low-energy structures–yielding −4.77, −3.10, −4.70 and −3.93 kcal mol−1 for the C2, D2d, Ci and C2h structures, respectively.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs16
214
At this level of theory, the D2h was not bound for any separation.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res27
215
This result confirms the behaviour of the HF+WL/6-31G* level as noted for the benzene dimer calculations.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res27
216
The D2h was checked instead using both the 6-31+G* and 6-311G* basis sets, with both yielding binding, optimised at −1.07 kcal mol−1 and −1.80 kcal mol−1, respectively.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res28
217
The HF+WL/6-31+G* also yields reasonable results for the low-energy structures, giving interaction energies of −5.83, −4.14, −5.75 and −4.86 kcal mol−1 for the C2, D2d, Ci and C2h structures respectively.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res29
218
This confirms the necessity for using at least split-valence plus diffuse functions in the basis set when using HF+WL.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con16
219
The final test system comprises the base-stacking interaction for dimers of pyrimidine, cytosine, uracil and guanine.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj15
220
The stabilisation of stacked base-pairs is thought to arise from van der Waals interactions.57
Type: Hypothesis |
Advantage: None |
Novelty: None |
ConceptID: Hyp2
221
The geometries of these systems are shown in Fig. 4.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod9
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac21
223
The pyrimidine system has been best characterised, at the CCSD(T)/aug-cc-pVDZ(+f functions and bond functions) level of theory,45 but only at one geometry (see Fig. 4) with the vertical separation fixed to 3.3 Å.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac21
224
For the remaining systems (fixed at the geometries in Fig. 4) examined by Hobza and Šponer, uracil was characterised up to the MP2/aug-cc-pVTZ(+ bond functions) level, cytosine up to the MP2/aug-cc-pVTZ leveland guanine up to the MP2/aug-cc-pVDZ level.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac21
225
In this by work of Hobza and Šponer,45 the vertical separation of cytosine and guanine dimers were fixed at 3.3 Å and 3.4 Å, respectively, and the uracil dimer geometry came from geometry optimisation at the MP2/6-31G** level.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac21
226
The minimum energy separations and interaction energies for these four systems, calculated at the HF+WL/aug-cc-pVDZ level, are reported in Table 6, as are ab initio data obtained by Hobza and Šponer.45
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs17
227
The ab initio data of Hobza and Šponer are the best estimates (obtained using an extrapolation procedure involving MP2 complete basis set data with CCSD(T) corrections) available at present.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac22
228
The agreement with the best energetic estimates from Hobza and Šponer is good for all cases.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res30
229
The trend in the interaction energies is reproduced, although the absolute values appear fortuitously close to the best estimates.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res30
230
It is interesting to note that the lowest-energy separation for pyrimidine is larger than the separation estimated by Hobza and Šponer.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res31
231
All other separations obtained using the HF+WL approach are, as expected, underestimated.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res32
232
These systems were also investigated using the xPBE and X3LYP functionals, both with the aug-cc-pVDZ basis set.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj16
233
The results are summarised in Table 6.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs18
234
Every structure was bound for both functionals, although the attraction is decreased in comparison with the data of Hobza and Šponer.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res33
235
Furthermore, the separation corresponding to the optimum interaction energy is increased relative to the results from Hobza and Šponer. xPBE manages to capture the trend of the interaction energies, although as already observed for the methane dimer, the range of interaction energies is considerably smaller than is seen in the ab initio results.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res34
236
X3LYP is seen to perform similarly to xPBE for these systems.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res35
237
The basis set dependency was again explored by repeating these calculations at the HF+WL/6-31G* level.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj17
238
Again, the vertical separation between monomers was optimised to yield the optimum interaction energy in each case.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp14
239
The results with the 6-31G* basis do not show much difference with the data obtained with the aug-cc-pVDZ basis.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res36
240
The interaction energies of pyrimidine, uracil, cytosine and guanine were −3.1, −11.3, −10.3 and −14.0 kcal mol−1, respectively.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs19
241
The corresponding optimised separations at the HF+WL/6-31G* level were 3.4, 2.97, 3.1 and 3.0 Å, respectively.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs20
242
This behaviour in stark contrast with the behaviour noted for the D2h isomer of the naphthalene dimer, and the F structure of the benzene dimer, both of which were unbound at the HF+WL/6-31G* level.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res37
243
It must be re-emphasised that, as a general rule, it appears that HF+WL calculations must be used with at least split-valence plus diffuse functions in order to capture qualitative behaviour.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con17
244
It is clear from all of the test cases examined that the HF+WL approach is capable of recovering qualitative trends in interaction energies for pure van der Waals (dispersion and repulsion) systems, albeit with shorter equilibrium distances (most of the time).
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con17
245
Relative orientations seemed to be captured well for the HF+WL approach.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con17
246
However, no claim is made to suggest that the HF+WL approach can describe a pure van der Waals system with quantitative accuracy, and therefore is not competitive with post-Hartree–Fock methodologies.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con17
247
To quantify the conditions under which this approach appears to be of value it is noted that, for example, in the case of the rare gas dimers, the HF+WL approach yields reasonable results for systems with overlaps no smaller than around 10−5 (with the exception of the helium dimer).
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con17
248
A HF+WL-type of approach may be extremely useful for optimising geometries of low-symmetry structures.
Type: Hypothesis |
Advantage: None |
Novelty: None |
ConceptID: Hyp3
249
Low-symmetry structures can be very time-consuming to optimise at any post-HF level (with a sufficiently large basis set), although advances are being made in this area (e.g. see .ref. 4)
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac23
250
The results of HF+WL optimisations would at least indicate which of the structures should be considered for further refinement by more expensive post-HF methods.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con18
251
For example, the location of the C2 structure of the naphthalene dimer was missed by several MP2-based geometry optimisation studies, presumably because of the low symmetry and the computational demands this creates.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con18
252
Another example is the geometry optimisation of small molecules adsorbed onto periodic substrates.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con18
253
Periodic EXX calculations are now possible via some mainstream plane-wave codes, making EXXWL presumably more efficient than currently available periodic MP2 implementations (as found in GAUSSIAN03).
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac24
254
It is also noted that the HF+WL approach appeared to work reliably only when a basis set featuring at least split-valence plus diffuse functions was used.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con19
255
The HF+WL/6-31G* level was demonstrated to be insufficient in several instances.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con19
256
Furthermore, these conclusions only hold at separations where the charge distribution on each fragment overlaps.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con19
257
The HF+WL approach cannot be successfully applied in the large separation regime.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con19
258
It is also clear that the current group of exchange–correlation functionals, B3LYP, PBE, X3LYP and xPBE, cannot be reliably used even for qualitative work.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con20
259
For these functionals, of those structures that showed any binding at all, the binding was substantially underestimated compared with best estimates from ab initio calculations.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con20
260
Two clear limitations of these current calculations are acknowledged.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac25
261
The first is that for all the systems studied here, the monomers were treated as rigid body fragments.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con21
262
However, the intramolecular degrees of freedom could in principle be optimised using the HF+WL approach.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con21
263
The second limitation is on the type of system that can be studied with this method (the same restrictions apply to HF+D methods), and is a consequence of how the correlation contribution forms a perturbative correction to the total Hartree–Fock energy.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con22
264
This limitation could be lifted if the WL correlation could be included dynamically, within Kohn–Sham DFT, via the implementation of exact-exchange methods (denoted here as the proposed EXXWL functional).
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con22
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac26
266
The fully self-consistent EXXWL functional would allow treatment of systems where a general HF+D approach currently fails (such as the NaAr complex, molecular physisorption on metal surfaces, etc)., where the HF wavefunction cannot distort if dynamic correlation effects are not accounted for.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con23
267
Work is currently underway to enable calculations with the EXXWL functional in future.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac27
268
Compared with the HF+WL approach, a much wider range of systems could be treated with the EXXWL functional.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac27
269
Furthermore, promising developments of inclusion of ‘left–right’ correlation with exact exchange within the KS-DFT framework62 suggest that an approach that combines non-dynamical correlation with the EXXWL functional may in future permit reasonable modelling of bond-breaking in weakly-bonded systems.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac28
270
Such an EXXWL functional is not intended to act as a replacement for more accurate post-HF methods, but rather as a device for optimising the structure and energetics of weakly-bonded systems that are not amenable even to MP2 optimisations.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac28
271
Compelling reasons for why the WL functional works so well in the cases studied here remain unclear.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con24
272
Of course, there is no way that the WL correlation functional can provide a realistic description of London dispersion interactions at long range.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con24
273
WL was parametrised from high-quality quantum chemical data for a small number of atoms.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met11
274
It also, by design, should work well with Hartree–Fock densities.
Type: Hypothesis |
Advantage: None |
Novelty: None |
ConceptID: Hyp4
275
Wilson and Levy point out that, due to the way in which the functional was parametrised, it may not perform well for systems where the Hartree–Fock density differs substantially from the true density.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac29
276
Because of the constraints imposed by the HF+D-like approach used here, the systems chosen for study should have HF densities that are reasonably close to the true density.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac29
277
Therefore, an investigation of systems where there is less similarity between HF and real densities is, at present, desirable.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac29
278
However, these systems should only be studied within a fully self-consistent implementation of the EXXWL functional.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac29
279
The WL correlation potential has been examined, and appears to have unremarkable properties in comparison with other correlation potentials.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con25
280
It is similar to other correlation potentials, e.g. the Perdew(86)63 correlation functional, in the sense that the potential decays to the asymptotic limit exponentially.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con25
281
This behaviour is reasonable in the sense that the correlation potential should decay faster than the exchange potential (in the case of the HF potential, the asymptotic decay goes like −1/r).
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con25
282
However, it has been remarked by others64 that perhaps a more realistic asymptotic behaviour of the correlation potential would be −α/r4, where α is the polarizability.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac30
283
The WL correlation functional is based on the local Wigner correlation functional30 (W38).
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met12
284
The LYP correlation functional,65 due to its basis on the Colle–Salvetti formula, is also based on the form of the W38 functional.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met13
285
The LYP functional comprises a (local) W38 component multiplied by another functional that contains gradient correction terms.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met13
286
This is in common with many exchange functionals that analogously comprise the Dirac (local) exchange, multiplied by another functional that contains gradient correction terms.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met13
287
Functionals of a form such as this lend themselves easily to an analysis in the low-density/high-density-gradient regime, since one can construct an ‘enhancement factor’ (the function that is multiplied by the local component of the functional).
Type: Method |
Advantage: Yes |
Novelty: Old |
ConceptID: Met13
288
Unfortunately, the WL functional is not multiplicatively separable in this sense from the W38 functional.
Type: Method |
Advantage: No |
Novelty: Old |
ConceptID: Met13
289
Since it is difficult to disentangle the W38 functional from the WL functional, the WL functional is therefore not amenable to an analogous analysis of the ‘enhancement factor’.
Type: Method |
Advantage: No |
Novelty: Old |
ConceptID: Met13
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac31
291
Finally, it is plausible that the WL functional yields unphysical attractive contributions for some of the short-ranged contributions to the interaction energy.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac32
292
A decomposition scheme based on Heitler–London type calculations, as reported by Cybulski and Seversen,13 may be helpful in this instance.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac33
293
The fact that the HF+WL almost always underestimates monomer separations seems consistent with this idea.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con26
294
Even if the origins of this unphysical behaviour can be revealed, further work is required to understand why these (potentially) unphysical terms combine such that very good interaction energies are ultimately recovered, albeit fortuitously.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con26
295
Thus, while admittedly not having a satisfying and physically rigorous theoretical justification, the HF+WL approach (or in future, the EXXWL functional) may become a useful pragmatic device for optimising the geometries of large and/or complicated van der Waals complexes that are currently beyond the reach of optimisation by more sophisticated ab initio methods.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con26
296
Other functionals, including those that are claimed to describe weak interactions such as xPBE and X3LYP, while also not capable of recovering London dispersion, cannot reliably perform the role of a pragmatic device, even when monomer charge distributions overlap.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con26
Conclusions
297
The Wilson–Levy (WL) correlation functional was used as a post-Hartree–Fock (HF) perturbative correction in the evaluation of the interaction energy for a range of weakly-bound dimers, in the regime of small to intermediate fragment separations.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj18
298
It was found that this HF+WL approach works well in a pragmatic and qualitative sense for the energetics of closed-shell van der Waals systems, but without a satisfying physical justification, and without quantitative accuracy.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con27
299
Two newly-developed exchange–correlation functionals, designed for applications to biological systems (X3LYP and xPBE), were also investigated and found to perform not as well as the HF+WL approach.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res38
300
The results of this work suggest an apparently promising device in the form of the EXXWL functional, for determining the energetic trends of molecules physisorbed on surfaces, nanotubes and other periodic structures.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con28