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" data-atf="true" src="cid:fbca3de0-0a46-4be1-a588-74b2b2971a58"><img alt="图片2.png" width="207px" height="70px" data-atf="true" src="cid:f963e023-a211-4ac7-81c9-032345b7b305"></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>

<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 on the web visit <a href="https://groups.google.com/d/msgid/cp2k/6e23b0a7-018d-430a-83ed-f8e31398963dn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/6e23b0a7-018d-430a-83ed-f8e31398963dn%40googlegroups.com</a>.<br />