DFT+D calculations

[bademoor]

Hello,

I'm currently looking into the possibility of using DFT+D calculations for zeolite systems. When using a pure GGA functional such as BLYP or PBE, the keyword $olddisp (Grimme 2004) or $disp (Grimme 2006) has to be added in the control-file to invoke the dispersion correction as mentioned in the TM-manual.

Here's my question:
If I want to use B97-D (implemented in Turbomole-5.10) as proposed by Grimme (2006), do I still need to add $disp in the control file or has this correction been added by default (without $disp in control file)?

Regards,

Bart

[uwe]

Hi,

$disp is always needed if you want to switch on the dispersion correction. B97-D without $disp is just B97-1 (but we do not recommend to use B97-1 for non-DFT+D calculations).

Regards,

Uwe

[cml]

Hello Bart,

 please be careful: in B97-D, the expansion parameters
in the expressions for exchange and correlation energy
have been refitted in the presence of the dispersion
correction (see J.Comput.Chem. 2006, p. 1787).

This means: if you do not turn on dispersion, you get
a modified version of B97-1, that will not yield the same
energies as the original functional ! 
Using B97-D without $disp therefore does not make much
sense and I strongly discourage you from doing so.

Best regards,

Christian Mück-Lichtenfeld

[bademoor]

Hello Uwe and Christian,

Thank you for the very clear comments. I'm happy to use B97-D functional now properly.
I think I was just confused by the name of the functional, B97-D: the "D" in this name wrongly gave me the impression that dispersion correction is already included (which is clearly NOT the case if one doesn't add "$disp" in the control file). Actually, it would be more logic if "$disp" was automatically added in control file if selecting B97-D, since the refitting of the B97-1 functional makes its use without "$disp" not really recommendable (following the comments of Christian).

Best regards,

Bart