<div>Hi Xiao,</div><div><br></div><div>As Lenard pointed out your SCF does not converge, consequently the forces are not converged and BFGS cannot optimize the structure.</div><div>You should adjust a few parameters:<br>MAX_SCF 200 ! change if needed to make sure SCF is really converged<br>ELECTRONIC_TEMPERATURE     2000 ! large smearing helps convergence. If needed you can reduce this once the optimization is converged and restart</div><div>ALPHA 0.01 !  you can play with this but it probably needs to be << 0.1</div><div><br></div><div>To speed up the first few steps of geometry optimization you could consider reducing the k-point grid to something like</div><div><br></div><div>    SCHEME  MONKHORST-PACK  2  3  1<br>    !FULL_GRID  .TRUE.</div><div><br></div><div>Once the the geometry is more reasonable you can switch back the the more accurate k-point sampling.<br></div><div><br></div><div>Cheers,</div><div>Fabian<br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Sunday, 28 February 2021 at 20:36:45 UTC+1 ASSIDUO Network wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">If I'm not mistaken, I think your structure is too far away from the optimized geometry. Your best bet would probably be to get the best structure possible before starting the geometry optimization.<br><br>This is true because I see the following in your output file:<br>At the 50th cycle your convergence criterion is around 0.9, nowhere near 1.0E-07, and you also get a message of "SCF run NOT converged".<br><br>I don't really see any issues with the flags you chose, so I think your issue is probably with the structure. The same thing happened to me the other day.<div class="gmail_quote"><div dir="auto" class="gmail_attr">On Saturday, February 27, 2021 at 7:15:49 PM UTC+2 <a href data-email-masked rel="nofollow">ti...@hotmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">




<div dir="ltr">
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Dear CP2K Users,</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I'm trying to run a geometry optimization for CoF2 (110) slabs, but the energy does not seem to decrease. Instead, it fluctuates a lot and sometimes increases significantly. I already used BFGS as the optimizer, and the eps_scf was set to 10^(-7). I'm wondering
 if anyone knows how to deal with such a problem. My input file, the generated xyz file, and the output file are attached below.</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Thanks in advance!</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Best regards,</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Xiao</div>
</div>

</blockquote></div></blockquote></div>