<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi August,<div><br></div><div>The correlation part seems to be missing from that input, try the following: FUNCTIONAL XC_MGGA_X_M06_L XC_MGGA_C_M06_L</div><div><br></div><div>With the M06_L PP parameters below and a couple quick calculations (500Ry, TZV2P, 10 Angstrom^3), I get </div><div><br></div><div>BE = 2*E_Monomer - E_Dimer = 5.48 kcal/mol.</div><div><br></div><div><br></div><div>chris</div><div><br></div><div><br></div><div>=============</div><div><br></div><div><div>#</div><div>H GTH-M06_L-q1</div><div> 1</div><div> 0.200000000000 2 -4.30444 0.799586</div><div> 0</div></div><div><div>#</div><div>O GTH-M06_L-q6</div><div> 2 4</div><div> 0.24084771 2 -17.00218284 2.49799489</div><div> 1</div><div> 0.22187778 1 18.01262971</div></div><div><br></div><div><br></div><div><br></div><div><br><div><div>On Jun 26, 2014, at 4:57 PM, August Melcher wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div dir="ltr">Hello again,<div><br></div><div>I'm running into a problem computing water dimer potential energy curves with GPW. I tried three different pseudopotentials (PBE, PBE0, and PBE "re-optimized" for M06-L) and two different basis sets (TZV2P and MOLOPT-TZV2PX) and all of them significantly underbind the dimer. </div><div><br></div><div>Attached are the curves as well as a sample input file. The reference data was taken from "Further analysis and comparative study of intermolecular interactions using dimers from the S22 database" Laszlo Fusti Molnar, Xiao He, Bing Wang and Kenneth M. Merz Jr. J. Chem. Phys. 131, 065102 (2009); <a href="http://dx.doi.org/10.1063/1.3173809">http://dx.doi.org/10.1063/1.3173809</a></div><div><br></div><div>I noticed the comment below in another thread. Could this be the cause of the underbinding I'm noticing?</div><div><br></div><div>"I looked at some of the source code files related to the calling of the XC functionalities. In particular, the xc_libxc.f file states the following as a comment at the beginning of the file.</div><br>!> WARNING: In the subroutine libxc_lsd_calc, it could be that the<br>!> ordering for the 1st index of v2lapltau, v2rholapl, v2rhotau,<br>!> v2sigmalapl and v2sigmatau is not correct. For the moment it does not<br>!> matter since the calculation of the 2nd derivatives for meta-GGA<br>!> functionals is not implemented in CP2K.<br><br>If this is the case, then it is possible that a meta-GGA functional (which generally is a function of the laplacian of the density) such as M06-2X may not work properly in CP2K. Sorry I could not be of more help."<div><br></div><div>Thanks,</div><div>August</div></div><div><br class="webkit-block-placeholder"></div>
-- <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+uns...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k">http://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout">https://groups.google.com/d/optout</a>.<br>
<span><m06-L.inp></span><span><pseudobench.xlsx></span></blockquote></div><br></div></body></html>