<div dir="ltr"><div>Your scf convergence criteria is not good enough, you should have EPS_SCF much smaller, maybe 10-6 to ensure energies are converged to at least micro hartree level.</div><div><br></div><div>I am not at all sure that the electric field and poisson solver work with k-points. I also don't have experience with the implicit psolver to comment on that setup. I'd remove that and check you can optimize a slab with periodic boundaries first then tackle the more complicated problem afterwards.</div><div><br></div><div>Matt</div><div><br></div><div>On Friday, June 19, 2020 at 8:20:26 AM UTC+1, yulei wang wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><span style="color:rgba(0,0,0,0.87);font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:14px">Hello all, </span><div style="color:rgba(0,0,0,0.87);font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:14px"><br></div><div style="color:rgba(0,0,0,0.87);font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:14px">I am fairly new to cp2k so I believe this problem is easily resolvable.  I am doing structural optimization of periodic copper surface under the action of electric field,  and I fixed 3 atomic layers below ( input file is cu.inp and date.coord).  I get the wrong results as the position of atom 7 is too high compared to other atoms in the same layer (the result is in <a href="http://cu-pos-1.xyz" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fcu-pos-1.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHaiph5DlPdxrAvDwkZnpybYObc7g';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fcu-pos-1.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHaiph5DlPdxrAvDwkZnpybYObc7g';return true;">cu-pos-1.xyz</a>). It is undoubtedly unreasonable, and then I tried to do the  geometry optimization without the action of electric field, and the result is just the same. Any ideas?</div><div style="color:rgba(0,0,0,0.87);font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:14px"><br></div><div style="color:rgba(0,0,0,0.87);font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:14px">Thank you for your help or advice. </div><div style="color:rgba(0,0,0,0.87);font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:14px"><br></div><div style="color:rgba(0,0,0,0.87);font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:14px">Best Regards, </div><div style="color:rgba(0,0,0,0.87);font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:14px">Yulei Wang</div></div></blockquote></div></div>