<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>