<div dir="ltr">As a first glance, your EPS_DEFAULT is 10⁻⁷ which is too big, I usually use 10⁻¹⁰ at least, or 10⁻¹² for more precise AIMD. Also a good relation between EPS_DEFAULT and EPS_SCF is to set the last one as the square root of the former.<div>All this affects the quality of the forces which dominates the MD... maybe if you check the geometry of the long SCF frames you could find some weird bond lengths or unphysical structures which can explain the behaviour.</div><div><br></div><div>Regards - Lucas</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El lun, 5 jun 2023 a las 14:09, Dune André (<<a href="mailto:dune.andre@gmail.com">dune.andre@gmail.com</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Dear all,<br></div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br></div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">I
am facing an issue while running AIMD simulations of a box of ~1000
atoms. The system is an amorphous MOF made of Fe linked by benzene
tricarboxylate (named Fe-BTC).</div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br></div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">I
use PBE functional with DZVP basis set and GTH pseudopotential ; DIIS
method and FULL_KINETIC preconditioner. I perform the calculation from a
restart file. The first steps of the MD are quite fast with ~10 SCF
steps per MD step. But along the AIMD run, the SCF cycles become longer
and longer (~100 SCF steps per MD step). <br></div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br></div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">I
tried to use other OT method such as CG or other preconditioner such as
FULL_SINGLE_INVERSE but in both cases the SCF cycles are taking too
much time and sometimes fail to converge. Besides, FULL_ALL
preconditioner is too computationally expensive in terms of memory.<br></div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br></div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Do you have some ideas of the origin of this increase of the SCF cycles? I attached the restart file to the message. <br></div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br></div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Thanks in advance for your help,</div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br></div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Kind regards,</div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br></div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Dune André<br></div>
<p></p>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/5d10d338-5a72-475a-adb4-d4ea7ecbd2d2n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/5d10d338-5a72-475a-adb4-d4ea7ecbd2d2n%40googlegroups.com</a>.<br>
</blockquote></div>
<p></p>
-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAOFT4PJAgDG-cAm4QS7HZV330bKfpF8D85RV_WR-2zK%2BpT2Y%2BA%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAOFT4PJAgDG-cAm4QS7HZV330bKfpF8D85RV_WR-2zK%2BpT2Y%2BA%40mail.gmail.com</a>.<br />