<div dir="ltr">Thank you very much.<br><br>在 2018年3月12日星期一 UTC+8下午9:10:18，Matt W写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi,<div><br></div><div>can't comment on what is / will be available in CP2K, but that water model was shown to be quite flawed, freezing at rather low temperatures.</div><div><br></div><div>There are some refitted parameters in the literature somewhere, I think. But I am not sure it is a very good idea for production calculations....</div><div><br></div><div>Matt<br><br>On Monday, March 12, 2018 at 12:27:22 PM UTC, <a>zuoh...@gmail.com</a> wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">If there is no function in current version, could developer add a morse potential with coulombic subtraction and add a parameter to control the contribution of the term.<div><br>在 2018年3月12日星期一 UTC+8下午8:04:24，<a>zuohen...@<wbr>gmail.com</a>写道：<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>I am trying to use FIST to do a box of water molecules with a force field (de Leeuw N. H., and Parker S. C. (1998a) Molecular-dynamics simulations of MgO surfaces in liquid water using a shell-model potential for water. Phys. Rev. B 58, 13901–13908.). </div><div><br></div><div>There is the coulombic subtraction in Morse potential for bond.</div><div><br></div><div>In CP2K the morse potential is in the form</div><div><br></div><div><em style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:medium">K1*[(1-exp(-K2*(R-R0)))^2-1])</em><span style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:small"> </span><br></div><div><span style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:small"><br></span></div><div><span style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:small">But I want to know how to add the the term of coulombic interaction.</span></div><div><br></div><div>I am not sure if it is possible in CP2K. </div><div><br></div><div>Please give me some suggestions. Thank you very much.</div></div></blockquote></div></div></blockquote></div></div></blockquote></div>