Dear CP2K developers,<div><br /></div><div>Hi, We are running a calculations for Cd33Se33 with Ferrocine ligands (containing Fe). The input and the structure are attached. However, when running a single-point calculations we face multiple problems. The first thing is that the &DIAGONALIZATIO approach does not converge and it gives large convergence values (like 100, 200, etc). Changing the basis set to DZVP or even TZVP did not have any effect on this procedure. We also tried adding smearing, increasing the cutoff, lowering the EPS_DEFAULT etc but no effects on the convergence. <br /><br />So, we started using the OT procedure and it converged the results in a large number of steps (about 500). However, when printing out the energy levels, we see that the LUMO energy is lower than the occupied orbitals including HOMO, HOMO-1, HOMO-2, and HOMO-3. The results seem to be no physical at all. So, even trying with unrestricted case, we get a similar thing. Is this a problem of the OT procedure? Because the diagonalization should print the eigenvalues in ascending order but in this case they are not. We appreciate any suggestions to helps us understand the problem. <br /><br />Thank you in advance.<br /><br />The attached files are the input, the coordinates, and one sample pdos file showing the energies of the molecular orbitals.</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/c0490612-eb84-4139-958f-dea542150ef4n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/c0490612-eb84-4139-958f-dea542150ef4n%40googlegroups.com</a>.<br />