<div dir="ltr">Dear Juerg,<div><br></div><div>Thank you for taking care of this. I just checked out the updated wB97X-V.inp from cp2k trunk, and it seems the parameters for SCALE_LONGRANGE and SCALE_COULOMB are reversed in this version. For this functional, SCALE_COULOMB should be 0.167 and SCALE_LONGRANGE should be 0.833. Would you mind fixing it again?</div><div><br></div><div>Thanks a lot,</div><div>Yuezhi</div><div><br>On Saturday, July 14, 2018 at 9:39:19 AM UTC-7, jgh wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Thank you. I have updated the test.
<br>
<br>best regards
<br>
<br>Juerg Hutter
<br>------------------------------<wbr>------------------------------<wbr>--
<br>Juerg Hutter                         Phone : ++41 44 635 4491
<br>Institut für Chemie C                FAX   : ++41 44 635 6838
<br>Universität Zürich                   E-mail: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="Q3PxS7E9BAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">hut...@chem.uzh.ch</a>
<br>Winterthurerstrasse 190
<br>CH-8057 Zürich, Switzerland
<br>------------------------------<wbr>------------------------------<wbr>---
<br>
<br>-----<a href="javascript:" target="_blank" gdf-obfuscated-mailto="Q3PxS7E9BAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a> wrote: -----
<br>To: "cp2k" <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="Q3PxS7E9BAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>>
<br>From: "Yuezhi Mao" 
<br>Sent by: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="Q3PxS7E9BAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>
<br>Date: 07/14/2018 04:17AM
<br>Subject: [CP2K:10530] wB97X-V in CP2K
<br>
<br>Dear CP2K developers,
<br>
<br>I was trying to figure out how to run wB97X-V calculations with CP2K, as AIMD with this functional using CP2K was reported earlier this year. By testing the binding energy for the water dimer in S22 set (the results are summarized in the attached .log file), I found the follows:
<br>(i) The setup of this functional in tests/QS/regtest-hybrid-4/<wbr>wB97X-V.inp seems incorrect. The line "FRACTION 0.167" should be deleted. As a mixed Coulomb operator is used and the 1/r and erf(wr)/r operators have been scaled separately, there is definitely no need to apply an additional scaling factor to scale down the exact exchange. From the output generated with the same setup as this regtest ("water_dimer_gapw_wB97X-V_<wbr>regtest_def2-qzvp.out"), we can find that HF exchange energy (as well as the total energy) is particularly small.
<br>
<br>(ii) Without the global scaling factor, the total energy given by CP2K (using the built-in parameterization) looks too low for this system (the corresponding output is "water_dimer_gapw_wB97X-V_<wbr>regtest_def2-qzvp.out") compared to the q-chem reference. It was reported before that wB97X-V, as a parameterization of BECKE97, was not implemented correctly because it was not attenuated as an SR functional (see a previous discussion on this issue: <a href="https://groups.google.com/forum/#!msg/cp2k/JEvnnfAEgPY/CFwlbsroSOcJ" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/forum/#!msg/cp2k/JEvnnfAEgPY/CFwlbsroSOcJ';return true;" onclick="this.href='https://groups.google.com/forum/#!msg/cp2k/JEvnnfAEgPY/CFwlbsroSOcJ';return true;">https://groups.google.com/<wbr>forum/#!msg/cp2k/JEvnnfAEgPY/<wbr>CFwlbsroSOcJ</a>). We suspected that this problem has not been fixed since then.
<br>
<br>(iii) Replacing the built-in implementation of the XC functional with the counterpart in libxc, the results produced by CP2K (especially the GAPW/def2-QZVP calculation) are in great agreement with the Q-Chem numbers. Therefore, I believe that the setup in "water_dimer_gapw_wB97X-V_<wbr>def2-qzvp.inp" is the correct way of running this functional.
<br>
<br>I suggest that the regtest job should be modified, as it is often the first thing a user would refer to when one wants to try out some new features. The current one turns out to be misleading. Also, I think the built-in implementation of this functional should be either fixed or disabled, as it produces incorrect results.
<br>
<br>Best regards,
<br>
<br>Yuezhi Mao
<br>Postdoctoral scholar
<br>Department of Chemistry, 
<br>Stanford University  
<br>  -- 
<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="javascript:" target="_blank" gdf-obfuscated-mailto="Q3PxS7E9BAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+uns...@googlegroups.<wbr>com</a>.
<br> To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="Q3PxS7E9BAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.
<br> Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>.
<br> For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.
<br> 
<br>
<br>[attachment "wB97X-V_water_dimer.log" removed by Jürg Hutter/at/UZH]
<br>[attachment "tests.zip" removed by Jürg Hutter/at/UZH]
<br></blockquote></div></div>