<div dir="ltr">Hi Natalie,<div><br></div><div>the warning is just for the calculation of the S**2 operator - i.e the expectation value of the total spin for the system. The value of it for DFT is a bit uncertain anyway...</div><div><br></div><div>The spin density is a totally different thing, and should be fine to use in analysis. Same goes for all properties that can be output as cube files I can think of.</div><div><br></div><div>Matt<br><br>On Wednesday, September 30, 2015 at 4:13:43 PM UTC+1, Natalie Austin wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hello,<br><br>I have a system consisting of Cu55<br>In this section I've included smearing and added mos: <br><br>  &SCF<br>     SCF_GUESS  ATOMIC<br>     EPS_SCF 1.0E-6<br>     MAX_SCF 500<br>     ADDED_MOS  200<br>     CHOLESKY INVERSE<br>     &SMEAR  ON<br>         METHOD FERMI_DIRAC<br>         ELECTRONIC_TEMPERATURE [K] 300<br>     &END SMEAR<br>     &DIAGONALIZATION<br>          ALGORITHM STANDARD<br>     &END DIAGONALIZATION<br>     &MIXING<br>          METHOD BROYDEN_MIXING<br>          ALPHA   0.1<br>          BETA    1.5<br>          NBROYDEN  8<br>     &END MIXING<br> &END SCF<br>UKS .TRUE.<br><br><br>At the end of my output file I get the warning: <br><br>  Integrated absolute spin density  :                             <wbr>  0.0003914226<br>  WARNING: S**2 computation does not yet treat fractional occupied orbitals<br><br>Is this warning relevant to just the spin density of does it have an effect on the total energy? <br>If this means the spin density is unrealible, how would I go about in my analysis of total density,homo-lumo gap. Basically, if I were to generate cube files based on this calculation would they be accurate?<br><br><br>Thanks,<br><br>Natalie<br><br></div></blockquote></div></div>