<div dir="ltr">Just to be a bit clearer - your ground state density is converged - which means that your occupied orbitals are converged.<div><br></div><div>But, the OT method only works with the occupied orbitals. To get unoccupied orbitals an extra computation needs to be done. </div><div><br></div><div>You ask for extra unoccupied states in the PDOS section (maybe elsewhere too, I've only looked very quickly).</div><div><br></div><div>CP2K then uses some algorithm, maybe Gram Schmitt or similar to build unoccupied orbitals that are orthogonal to the occupied ones.<br></div><div><br></div><div>This is the bit of the calculation that is not converging and where Samuel suggests you allow CP2K to carry out extra iterations.</div><div><br></div><div>HTH,</div><div><br></div><div>Matt<br><br>On Monday, October 26, 2015 at 9:28:46 AM UTC, Samuel Andermatt wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Well ot might not be converged, but your scf is converged.<br>*** SCF run converged in     1 steps ***<br>and<br>  outer SCF iter =    1 RMS gradient =   0.80E-07 energy =       -784.0104898082<br>  outer SCF loop converged in   1 iterations or    1 steps<br><br>Therefore restarting with the same density doesn't help.<br><br>You need to change max_:iter_lumos in the scf section.<br></div></blockquote></div></div>