<div dir="ltr">Dear CP2K developers, <br>I want to report a serious problem with the preconditioner when computing the energy of Fe-MOF-74 with UKS settings.<br>Long story short the energy goes down, it is close to convergence but then, after recomputing the preconditioner in the outer step it starts to diverge: see the image.<br><div><br></div><div><a imageanchor="1" href="about:invalid#zClosurez" style="clear: left; margin-bottom: 1em; float: left; margin-right: 1em;"><img src="cid:9f22886b-88cc-4ec0-9c51-1e023ac018da" alt="image.png" width="397" height="400" style="" border="0"></a><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>Now, I attached the input and the output, but I already tried very different settings:<br>-smearing<br>-diagonalization without smearing<br>-different BS<br>-UKS false (converged, but to a higher energy, forcing the closed shell)<br>- starting magnetization on metals<br>- CG minimizer<br>- lower multiplicity<br>- 2x and 3x replicated unit cell in the shortest dimension<br>- lowering ENERGY_GAP to 0.001<br>These calculation ALL failed (I made 16 attempts in total, changing settings!).<br><br>And at the end the only one working (the 17th attempt, ironically the unlucky number for Italians!) was to use CG minimizer with <u>no outer SCF</u>, hence not recomputing the preconditioner. It converged after ~600 inner steps, to:<br></div><div><br></div><div><div> HOMO - LUMO gap [eV] :    1.917671</div><div> HOMO - LUMO gap [eV] :    0.235602</div><div><br></div><div> ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -1170.045288008604075</div></div><div><br></div><div>Please let me know what was the problem, if it is a known issue of the preconditioner and how to circumvent it: for the moment I'm using no outer steps to have a robust convergence, but I know that this may cause the convergence to a local miniimum!<br><br>Thanks a lot!<br><br>Daniele <br>PhD candidate, LSMO, EPFL Sion</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><p class="separator" style="text-align: center; clear: both;"></p></div>