<div dir="ltr">Hi,<br><br>Maybe it is right. But I do static simulation, and I need to read pot.cube just for once.<br>You can have a try.<br><br><br>在 2017年12月2日星期六 UTC+8上午3:23:12，Jian Jiang写道：<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>Which means I need to update the "READ_FROM_CUBE" file every step, right? Cuz I am trying to do MD simulation</div><div><br></div><div>Best,</div><div>Jian Jiang<br><br>On Friday, December 1, 2017 at 2:47:33 AM UTC-6, <a>jts2...@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">Hi, <div><br></div><div>Maybe you need to supply a potential file and using "READ_FROM_CUBE" keywords to let CP2K  know the potential at grid points. <br><br>在 2017年11月30日星期四 UTC+8上午5:24:40，Jian Jiang写道：<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>I am not sure if "EXTERNAL_POTENTIAL" works in my case.</div><div>Can I modify EXTERNAL_POTENTIAL to let different atoms feel different potential? Like O and H have different parameters(Well depth,<span style="font-family:sans-serif;font-size:14px">equilibrium bond distance</span>) in term of Morse potential.</div><div>Thank you.</div><div><br>On Saturday, November 25, 2017 at 6:04:42 AM UTC-6, <a>jts2...@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">Hi Jiang,<br>You can search key words "EXTERNAL_POTENTIAL" in the manual page. It will tell you how to add extra potential in CP2K.<br><br>Best wishes!<br><br>在 2017年11月22日星期三 UTC+8上午4:30:59，Jian Jiang写道：<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi guys,<div><br></div><div>I want to introduce smooth walls in z direction, which are going to interact with QM atoms using Morse Potential in NPT DFT MD.</div><div>Which files(force?energy?) are supposed to be modified and recompiled?</div><div><br></div><div>Thank you in advance! Any help would be highly appreciated!</div><div><br></div><div><br></div><div>Best,</div><div>Jian </div></div></blockquote></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div>