<div dir="ltr">Try using the GAPW method (change GPW to GAPW in &QS section ). <div><br></div><div>Na is problematic, the pseudopotential is extremely hard.</div><div><br></div><div>Matt<br><div><br></div><div><br>On Thursday, July 13, 2017 at 4:41:39 PM UTC+2, Hongliu wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Matt,<br><br>it really helps! Now the structure changes indeed. What could be the reason for that?<br><br>But it changes too much :(.  For a 3-layer (111) slab I fixed the first layer and the upper layer becomes a mass after a while. BTW, for the cell optimization of Na bcc crystal structure I got always a disordered state. This was also reported by <span><span><span style="color:rgb(34,34,34)">Aniruddha Dive</span></span></span> at this post: <a href="https://groups.google.com/forum/#!topic/cp2k/aHejgB31V7s" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/forum/#!topic/cp2k/aHejgB31V7s';return true;" onclick="this.href='https://groups.google.com/forum/#!topic/cp2k/aHejgB31V7s';return true;">https://groups.google.com/<wbr>forum/#!topic/cp2k/aHejgB31V7s</a><wbr>. <br><br>My input files for the cell optimization are attached. Any suggestions for that? Thanks!<br><br>Best wishes,<br>Hongliu<br> <br>On Wednesday, July 12, 2017 at 10:21:04 PM UTC+2, Matt W wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Not quite sure. <div><br></div><div>Maybe try deleting all of your GEO_OPT section (i.e. going back to the defaults) and see if that works better?<div><br></div><div>Matt<br><br>On Wednesday, July 12, 2017 at 4:55:12 PM UTC+2, Hongliu wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear CP2Ker:<br><br>I was studying the binding energy of a Na atom to a Na slab. I tried first to optimize the structure. When I used the attached input file the structure of the system changed only very slightly, even if I put the last Na atom to a unrealistic position. I am not sure what goes wrong with my files. Please help me out. Thanks!<br><br>BTW, I am using version 2.6.1.<br><br>Best wishes,<br>Hongliu<br></div></blockquote></div></div></div></blockquote></div></blockquote></div></div></div>