<div dir="ltr">Hi Matt<div><br></div><div>Thank you for your help. I had tried to set EPS_SCF as 1.0E-06, but the result is similar to before.  </div><div><br></div><div>And I have a question about POISSON_SOLVER PERIODIC. If I just want to repeat the periodicity in both X and Y directions, it will report an error as "Illegal combination of periodicity and Poisson solver periodic3d". How can I setting parameters to achieve my goal.</div><div><br></div><div>Thanks</div><div>Yulei Wang<br><br>在 2020年6月19日星期五 UTC+8下午6:23:43，Matt W写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><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" rel="nofollow" target="_blank" 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></blockquote></div></div>