<div dir="ltr">Some suggestions for what to try:<div><br></div><div>(1) Precision</div><div> Your value of EPS_PGF_ORB is too large and should be substantially decreased. I think the default value is too large for a lot of systems. Most numerical issues stem from the overlap matrix precision in my opinion. Set it to something like 1e-16 and see what happens. That might be too fine, but its not very much extra computational cost and it can make a huge difference.<br></div><div><br></div><div>(2) Try these settings</div><div><br></div><div>&OT</div><div> algorithm irac # do not enforce strict orthogonality</div><div> minimizer cg # safer than DIIS, though often its too safe and makes the calculation slow.</div><div> linesearch 2pnt # 3pnt is not compatible with irac+rotation</div><div> rotation true # allow for subspace rotations (fractional occupations)</div><div> preconditioner full_single_inverse # full_all not compatible with rotations</div><div>&end</div><div><br></div><div>If that doesn't work, decrease STEPSIZE in &OT to be 0.05 </div><div><br></div><div>Hope that helps,</div><div>Nick</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jun 21, 2022 at 8:45 AM 'Börries von Seggern' via cp2k <<a href="mailto:cp2k@googlegroups.com">cp2k@googlegroups.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div>Dear Cp2k Forum,</div><div>For a while now, I have been trying to calculate XAS Spectra of transition metal oxides using the LR-TDDFT method. Unfortunately, I wasn`t able to make the preceeding Hybrid-DFT SCF calculation converge for open-shell systems (CuO, CoO, NiO), while other closed shell systems (MgO, ZnO) show rapid convergence quite easily. So far, I've tried using the OT methods, which do seem to converge until reaching a sort of plateau, where the energy starts increasing in miniscule amounts. Further, Iǜe been trying to use diaganolization, however these calculations seem to oscillate even when including very little amounts (1e-2) of the new density.</div><div>Concerning the OT methods, I ve tried both DIIS and CG (with 2pnt and 3pnt linesearch), however, all calculations show the same behaviour.</div><div>Could anyone please help me with these calculations?</div><div><br></div><div>Best regards,</div><div><br></div><div>Börries</div><div><br></div><div>PS: Please find an example input, output and structure files attached, I am using cp2k - 8.2. As even this calculation using the CG optimizer with linesearch 3pnt reached the plateau and was very expensive, I unfortunately cancelled the calculation. (But could restart it if necessary) <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/149a0cc5-5ea3-4a94-b3ea-b893b725b9ddn%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/149a0cc5-5ea3-4a94-b3ea-b893b725b9ddn%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/CAE8AJWGYVOg73%3DEyF5-AOacrNruq4Rwv9rD%2BqMvkKiZModm7vw%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAE8AJWGYVOg73%3DEyF5-AOacrNruq4Rwv9rD%2BqMvkKiZModm7vw%40mail.gmail.com</a>.<br />