<p>CP2K version 2025.2</p>
<p> Dear all, </p>
<p> I am running some test calculations using the D4 dispersion
together with the m06 functional. </p>
<p>I am confused about some of the parameters in the input and
how these correspond to information in the output. </p>
<p>My VDW_POTENTIAL section looks like the following:</p>
<p> </p>
<p> &VDW_POTENTIAL<br />
POTENTIAL_TYPE PAIR_POTENTIAL<br />
&PAIR_POTENTIAL<br />
TYPE DFTD4<br />
D4_CN_CUTOFF 5.29177209<br />
D4_CUTOFF 1.05835442E+001<br />
D4_DEBUG T<br />
D4_REFERENCE_CODE F<br />
FACTOR_S9_TERM 1.00000000E+000<br />
REFERENCE_FUNCTIONAL m06<br />
R_CUTOFF 1.05835442E+001<br />
VERBOSE_OUTPUT T<br />
&PRINT_DFTD DEBUG<br />
&END PRINT_DFTD<br />
&END PAIR_POTENTIAL<br />
&END VDW_POTENTIAL</p>
<p> </p>
<p>The value of D4_CN_CUTOFF is the suggested value in the
manual. </p>
<p> </p>
<p>In the output, the following information is printed:</p>
<p> </p>
<p> !-----------------------------------------------------------------------------!<br />
DEBUG D4 DISPERSION<br />
!-----------------------------------------------------------------------------!<br />
DEBUG D4| Reference functional
m06<br />
DEBUG D4| Scaling parameter (s6)
1.0000<br />
DEBUG D4| Scaling parameter (s8)
0.1637<br />
DEBUG D4| BJ Damping parameter (a1)
0.5346<br />
DEBUG D4| BJ Damping parameter (a2)
6.0619<br />
DEBUG D4| Cutoff value coordination numbers
0.1000E-05<br />
DEBUG D4| Cutoff radius coordination numbers
11.1000<br />
DEBUG D4| Coordination number function type
3<br />
DEBUG D4| Cutoff radius 2-body terms [bohr]
40.0000<br />
DEBUG D4| Cutoff radius 3-body terms [bohr]
40.0000<br />
DEBUG D4| CN differences (max)
0.000001<br />
DEBUG D4| CN differences (ave)
0.000000</p>
<p> </p>
<p>I have the following questions regarding this output
section:</p><p style="margin-left: 18pt;"></p><ol><li>I don't understand what "Cutoff
radius coordination numbers" means. I could guess that it is the
value of the input parameter D4_CN_CUTOFF, but the value does not match. If D4_CN_CUTOFF
is converted from Ångstrøm to Bohr in the output, the value should be 10 and
not 11.1. Is the value of D4_CN_CUTOFF in the input ignored?</li><li>What does “Coordination number
function type” refer to?</li><li>By comparing the values of the
scaling and BJ damping parameters to the values of the dftd4 Github web page, I
see that the ATM values are used. Is there an input option for choosing the
many-body dispersion parameters instead?</li></ol><p></p><p>
</p>
<p style="text-indent: -18pt;"><br /></p>
<p>Thanks in advance for your help. </p>
<p>Best regards,</p>
<p>Torstein</p>
<p> </p>
<p></p>
-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/aad7a74c-a98b-40fb-a228-d155b83b890dn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/aad7a74c-a98b-40fb-a228-d155b83b890dn%40googlegroups.com</a>.<br />