Chemistry - Do modern dispersion-corrected DFT methods give more accurate molecular geometries?

Solution 1:

I am afraid I do not know of a paper specifically discussing the effects of DFT dispersion corrections on intramolecular properties. However, it should be noted that for short-range descriptors (such as bond lengths and bending angles) the dispersion corrections will not influence those as they correspond to distances below the typical distance at which the corrections take place. The various schemes differ in how they work, but typically the corrections (in Grimme's schemes) are damped below a distance of ~6 Å. It is also my experience that such local structural characteristics are not affected by dispersion corrections.

However, you are right that dispersion corrections will influence intramolecular structure at larger scale, and in particular for “soft” or “flexible” long molecules with many conformations. There, it is worth noting that the so-called “dispersion” corrections actually correspond to the lack of accounting in local DFT functionals of middle-to-long range correlation effects. Thus, it is actually expected that the corrections — or nonlocal functionals including dispersion terms — actually perform better for those too.

Unless, like Grimme's 2006 two-body dispersion correction scheme (often called “D2”), the dispersion corrections happen to be often too strong and make intramolecular distance somewhat too small (instead of much too large)!


Edit: Now that I come to think of it, one extreme example of intramolecular system that has been studied in detail w.r.t. dispersion corrections is that of helicenes. It was shown (for example there and there) that dispersion corrections improve the accuracy of computational structures and energies compared to experimental (or high-level computational) data.

Solution 2:

Disclaimer: These articles are from my old group and myself. But they may serve as a starting point for someone's own investigation, even if one disagrees with our conclusions.

The Grimme group, me included, looked at experimental, back-corrected rotational constants of smallish to medium-sized organic molecules a while back. Our conclusion was that a modern dispersion-correction, i.e. D3(BJ) is beneficial for geometry optimization at virtually non-existent computational cost. The references are: T. Risthaus, M. Steinmetz, S. Grimme, J. Comput Chem., 35, 1509 (2014), M. Steinmetz, S. Grimme, Phys. Chem. Chem. Phys., 15, 16031 (2013).


Solution 3:

Our group just published a benchmark, considering ~6500 conformer geometries across ~650 molecules. Using high-level DLPNO-CCSD(T) / def2-TZVP energies, we found a huge increase in accuracy (from several metrics) when using dispersion correction at minuscule cost in time.

'Assessing conformer energies using electronic structure and machine learning methods' Int. J. Quantum Chem. 2020

For example, the mean $R^2$ correlation between B3LYP and DLPNO-CCSD(T) energies was 0.706 without dispersion and 0.920 with -D3BJ dispersion correction. Similar effects were seen with other functionals.

Our conclusion is that intramolecular non-bonded interactions are handled better, and as mentioned by F'x reflect medium-to-long range behavior of standard chemical functionals.

tl;dr - use dispersion correction for any density functional calculation. Even on molecular geometries it's generally better.