Hello, it is ok for seting a potential to prevent molecules motion during AIMD simulations? I get the same problem. <br /><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2022年4月21日星期四 UTC+8 18:42:13<D Y> 写道：<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear cp2k users and experts,<div>  I want to set up a function with spherical symmetry in CP2K_INPUT
/ FORCE_EVAL / EXTERNAL_POTENTIAL section during AIMD simulation, to prevent particles escape form slab. So i looked up the xTB manual for reference, and found there were two ways to set up such potential. one is logfermi potential, the other is polynomial potential. </div><div><br><img alt="图片1.png" width="416px" height="64px" src="https://groups.google.com/group/cp2k/attach/8e00eae770221/%E5%9B%BE%E7%89%871.png?part=0.2&view=1"><img alt="图片2.png" width="207px" height="70px" src="https://groups.google.com/group/cp2k/attach/8e00eae770221/%E5%9B%BE%E7%89%872.png?part=0.1&view=1"></div><div>Hence, i set up something like these in cp2k input file:</div><div>1. &EXTERNAL_POTENTIAL</div><div>        ATOMS_LIST   129..138<br>        FUNCTION   A*(sqrt((X-a)^2+(Y-b)^2+(Z-c)^2)/R)^8<br>        PARAMETERS A R a b c<br>        VALUES     1.0E-2 5 9.2 10.8 5.6<br>        UNITS  eV angstrom angstrom angstrom angstrom<br>  &END EXTERNAL_POTENTIAL<br></div><div><br></div><div>2.   &EXTERNAL_POTENTIAL</div>        ATOMS_LIST   129..138<br>        FUNCTION   A*log10(1+exp(8*sqrt((X-a)^2+(Y-b)^2+(Z-c)^2)-R))<br>        PARAMETERS A R a b c<br>        VALUES     1.0E-2 5 9.2 10.8 5.6<br>        UNITS  eV angstrom angstrom angstrom angstrom<br>  &END EXTERNAL_POTENTIAL<div><br></div><div>I really want to know whether the above mentioned functions are correctly setted? if it is wrong, could you please tell me  how to correctly set up this function?<br><div><br></div></div></blockquote></div>

<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/0e3d1f35-08bb-484e-a868-0fe2d6792c5an%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/0e3d1f35-08bb-484e-a868-0fe2d6792c5an%40googlegroups.com</a>.<br />