Dear cp2k,<div><br /></div><div>I am no using CP2K-v2025.1 to perform the GEO_OPT for my system (FeN4). I am now using DIAG algorithm with 3*3*1 MONKHORST PACK grid, along with FERMI-DIRAC smearing method. But I find that the SCF is very hard to converge each ionic step of 130, even with DZVP-MOLOPT-SR-GTH basis set. </div><div><br /></div><div>Previous tests I have found that the DIAG is always harder to converge than OT, but OT cannot support k-points, so I have to switch to DIAG with k-points. And, after the GEO_OPT with DZVP-MOLOPT-SR-GTH, if I want to calculate with TZVP-MOLOPT-GTH for highly accurate energy, hard convergence again. </div><div><br /></div><div>Attached is my .inp file and the unfinished .out file. I want to know if some other good solutions. Many thanks~</div><div><br /></div><div>I learn that if the MOLOPT basis set is not very suitable for slab of surface calculations? Because I find that the same system is very easy to converge in VASP. <br /><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">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/9ca211f5-c29f-4883-988d-da5e9b592e80n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/9ca211f5-c29f-4883-988d-da5e9b592e80n%40googlegroups.com</a>.<br />