<div dir="ltr"><div>Hi Marcella,</div><div><br></div><div>One of my students has encountered this problem. He tried changing the OT and increasing the number of MOs (we don't know what you mean by "smearing the occupation numbers around the Fermi energy") but it didn't correct the problem. I’m posting this on his behalf because he was unable to- he got a message that he had “limited access to public groups in his domain.”</div><div><br></div><div>Do you (or anyone else) have another suggestion? In his pdos files all the levels look fine except for the HOMO and LUMO.</div><div><br></div><div>From his output:</div><div> Fermi Energy [eV] :   -5.966796</div><div>  Lowest Eigenvalues of the unoccupied subspace spin            1</div><div> -----------------------------------------------------</div><div>  Reached convergence in          149  iterations </div><div>      -0.23489989     -0.21167052     -0.18538586     -0.17943768</div><div>      -0.13796253     -0.13007839</div><div> HOMO - LUMO gap [eV] :   -0.425155</div><div> ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -1076.972424168399129</div><div><br></div><div>Best regards,</div><div>Ian Haminton</div><br>On Sunday, November 22, 2009 at 6:13:06 AM UTC-5, Claudio wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Dear all,
<br>
<br>sometimes after I converge a calculation to an accuracy of say 1.0E-5
<br>I get things like
<br>
<br>HOMO - LUMO gap [eV] :   -0.980594
<br>
<br>This is because the Homo Energy is 0.22417642 and the Lumo Energy is
<br>0.18814027 in AU.
<br>
<br>I have requested the same accuracy for homo and lumo  and have added
<br>enough iterations to the lumo calculation loop to make sure the
<br>program claims it is converged. I am using the OT method. How do I
<br>prevent this from happening? I mean how do I ensure that the
<br>calculation yields Lumo energies that are higher than the homo energy?
<br>
<br>
<br>Sorry to bother you and thanks for your help!
<br>
<br>Claudio
<br>
<br></blockquote></div>