Here’s one more attempt to discern the failure of DFT to handle simple alkanes (see this earlier post for a previous attempt to answer this question). Tsuneda and co-workers1 have employed long-range corrected (LC) DFT to the problem of the energy associated with “protobranching”, i.e., from the reaction
CH3(CH2)nCH3 + n CH4 → (n+1) CH3CH3
They computed the energy of this reaction for the normal alkanes propane through decane using a variety of functionals, and compared these computed values with experimentally-derived energies. Table 1 gives the mean unsigned error for a few of the functionals. The prefix “LC” indicated inclusion of long-range corrections, “LCgau” indicates the LC scheme with a gaussian attenuation, and “LRD” indicates inclusion of long-range dispersion.
Table 1. Mean unsigned errors of the “protobranching”
reaction energy of various functional compared to experiment.
Functional |
MUE |
|
|
LC2gau-BLYP+LRD |
0.09 |
LC-PBE+LRD |
0.17 |
SVWN5 |
0.27 |
LCgau-PBE |
1.56 |
M06 |
1.98 |
LC-PBE |
2.24 |
M06-2x |
3.40 |
B3LYP |
5.97 |
HF |
6.96 |
|
A number of important conclusions can be drawn. First, with both LC and LRD very nice agreement with experiment can be had. If only LC is included, the error increases on average by over 1 kcal mol-1. The MO6-2x functional, touted as a fix of the problem, does not provide complete correction, though it is vastly superior to B3LYP and other hybrid functionals. The authors conclude that the need for LC incorporation points out that the exchange functional lacks the ability to account for this effect. Medium-range correlation is not the main source of the problem as large discrepancies in the reaction energy error occur when different functionals are used that are corrected for LC and LRD. Choice of functional still matters, but LC correction appears to be a main culprit and further studies of its addition to standard functionals would be most helpful.
References
(1) Song, J.-W.; Tsuneda, T.; Sato, T.; Hirao, K., "Calculations of Alkane Energies Using Long-Range Corrected DFT Combined with Intramolecular van der Waals Correlation," Org. Lett. 2010, 12, 1440–1443, DOI: 10.1021/ol100082z
Henry Rzepa responded on 28 May 2010 at 6:39 am #
Its not just proto-branching where long range corrections are needed. In an entirely different context, that of predicting CD (circular dichroism) spectra as a reference for assigning absolute configurations of molecules, LC corrections can improve the match with experiment significantly. One method of choice is CAM-B3LYP, which includes the Handy long range corrections (this will be reported in an article currently in press).
Steven Bachrach responded on 28 May 2010 at 7:54 am #
Be sure to let me know when and where the article on LC effects on CD comes out. I’ll want to blog it!
Henry Rzepa responded on 30 May 2010 at 2:58 am #
An addendum to the above observation. Another chiro-optical property where long range corrections seem to improve matters is optical rotation. M-hexahelicene has an observed rotation of -3640° at 589 nm in chloroform. A solvation corrected, frequency dependent, B3LYP/6-31G(d,p) calculation predicts -4456°. The Cam-B3LYP value is -3225°. It is known that improving the basis set quality will increase the absolute value of the rotation (by ~200° in the limit for this system) which improves the fit further. For details of such a calculation, see 10042/to-4945.
Henry Rzepa responded on 01 Jun 2010 at 2:18 am #
I recently came across this article by Hans-Joachim Werner (DOI: 10.1063/1.3160675) in which linear scaling of the computational cost with molecular size is claimed, with an RMS error of 0.6 kJ/mol for 54 different reaction energies, and a largest deviation of 2.5 kJ/mol. It is based on a local approximation to the coupled-cluster CCSD method, and as such offers an interesting alternative to the more empirical DFT approach, with its plethora of short and long range corrections for correlation/dispersion effects. Is a strong competitor to the DFT methods for larger molecules on the horizon?
Computational Organic Chemistry » Origin of DFT failures – part III responded on 08 Nov 2010 at 1:25 pm #
[…] I have blogged about many times) has recently been attributed to long-range exchange1 (see this post) or simply just DFT exchange2 (see this post). Grimme now responds by emphatically claiming that it […]
Computational Organic Chemistry » AIM analysis and molecular structure responded on 03 Dec 2010 at 8:29 am #
[…] the ortho hydrogens of biphenyl or between hydrogens across the bay region of phenanthrene (see this post). He argues that the barrier for biphenyl rotation arises from the stretching of the C-C bond when […]