<div dir="ltr"><span style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;">Dear Carlo,</span><br style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;"><br style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;"><span style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;">since you used the smearing of the occupation numbers, the multiplicity is not constrained. </span><br style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;"><span style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;">The occupation of the eigenstates is assigned according to their energy, irrespective of alpha and beta.</span><br style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;"><span style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;">The Fermi energy is then determined consistently for the full electronic system, and it is -0.14381945698278 a.u. in your case.</span><br style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;"><br style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;"><span style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;">If this is really meaningful for your system, I do not know. </span><br style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;"><span style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;">Is this really a metallic system, as it apparently results from this calculation?</span><br style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;"><span style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;">Is PBE0 appropriate for a metallic system?</span><br style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;"><br style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;"><span style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;">ciao</span><br style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;"><span style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px;">Marcella</span><br><br>Il giorno martedì 14 gennaio 2014 12:08:05 UTC+1, c.pignedoli ha scritto:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear all,<div> I am doing a PBE0 calculation for a periodic quasi 1D system (a graphene nanoribbon)</div><div><br></div><div>I did restart the calculation from a converged antiferromagnetic PBE calculation</div>

<div><br></div><div>The PBE0 calculation terminated with a unequal number of alpha and beta orbitals thus violating (I guess) the MULTIPLICITY 1 of the input and of the restart.</div><div>(</div><div>355 alpha and 349 beta while it should be  352 alpha and 352 beta</div>
<div><div>this can be seen checking   the added mos</div><div>declared at the end of the file: should have been 100 and 100 instead it is 97 102</div><div>and form the file with occupation numbers of the states (which I do not include)</div>
</div><div>) </div><div><br></div><div>I enclose the input file and the output, what is scaring is a sudden jump of the energy at </div>
<div>the beginning of the SCF (from -980 to -1200 at step 24 of the scf) that could be related to problems with the coulomb cutoff.</div><div><br></div><div><br></div><div><br></div><div>Does anybody have suggestions about what could be wrong?</div>

<div>Could it be that the cell size is too small for this particular system (however it is big enough</div><div>to have a fair PBE description)?</div><div>(Basically I have 12 unit cells of a zig-zag graphene nanoribbon with at least 20A of vacuum in the non periodic directions, the PBE UKS calculationn produces the correct edge states but of course these states are nearly metallic and have, in PBE a small gap)</div>
<div><br></div><div>Kind regards</div><div><br></div><div><br></div><div>Carlo</div></div>
</blockquote></div>