<div dir="ltr">Dear Katarina,<br><br>first of all, most results are not that different. <br>1) Regarding the literature. Surely, different basis sets have been used. Also please note that CP2K employs pseudopotentials. Truhlar has performed all-electron calculations. Thus, I would say that the results rather agree, than not.<br>2) As for the tests, there will be some tolerance margin differing on you compilation, system etc. So, only the last number looks suspicious, although not really.<br><br>Yours,<br><br>Vladimir<br><br>вторник, 17 января 2017 г., 14:07:00 UTC+1 пользователь katarína stančiaková написал:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear CP2K users,<br>I did some basic testing on bond dissociation energies on a small set of molecules using PBE functional.<br>I compared my results with literature (<a href="http://dx.doi.org/10.1007/s00214-007-0401-8" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fdx.doi.org%2F10.1007%2Fs00214-007-0401-8\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGvnQruapicy7fydx9968R7v5iKLw';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fdx.doi.org%2F10.1007%2Fs00214-007-0401-8\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGvnQruapicy7fydx9968R7v5iKLw';return true;">http://dx.doi.org/10.1007/<wbr>s00214-007-0401-8</a> ) and with results obtained using ADF. <br><br>However CP2K results are quite different from my reference and ADF. I tried to play with settings (e.g. turn off the periodicity), but this did not lead to any improvement.<br>Afterwards I decided to do some basic testing, using test files. Especially I tested following input files:<br><br>/tests/QS/regtest-hybrid-4/<wbr>CH3-PBE0_TC.inp<br>/tests/QS/regtest-hybrid-4/<wbr>CH3-PBE0_TC_LRC.inp<br>/tests/QS/regtest-hybrid-4/<wbr>CH4-PBE0_TC.inp<br>/tests/QS/regtest-hybrid-4/<wbr>CH4-PBE0_TC_LRC.inp<br><br>However, for CH4 case I obtained quite a difference between the numbers given in txt TEST_FILES, where I assume numbers which I should obtain are stores. I tested two versions of CP2K (4.0 and 2.6), and numbers are still the same. Now I am struggling what might be wrong. Please, do you have any suggestions?<br><br>Below you can find example of input files and my data<br><br>(for all cases we assume homolytic dissociation)<br><br>CH3-CH3 -> CH3 + CH3        92.59 (96.79) kcal/mol <br>CH3-OCH3 -> CH3 + OCH3       81.48 (87.24) kcal/mol<br>CH3-isopropyl -> CH3 + isopropyl         86.01 (89.65) kcal/mol<br>OCH3-isopropyl -> OCH3 + isopropyl                   77.69 (84.08) kcal/mol<br><br>Numbers in parentheses are from literature.<br><br>The results from testing files:<br>CH3-PBE0_TC.inp         -7.98384  (-7.94392) a.u.<br>CH3-PBE0_TC_LRC.inp    -7.34353  (-7.34351) a.u.<br>cCH4-PBE0_TC.inp        -7.35850  (-7.35848) a.u.<br>CH4-PBE0_TC_LRC.inp    -8.00130  (-7.96099) a.u.<br>Numbers in parentheses are from corresponding TXT file TEST_FILES (/tests/QS/regtest-hybrid-4/)<br><br><br>Thank you for any help<br>Best,<br>Katarina<br><br></div></blockquote></div>