<div dir="ltr">Hi Natalie,<div><br></div><div>if your system is metallic, and you are using smearing, there is not necessarily any spin polarization, even if there is an odd number of electrons in the system. Imagine running a system with 27 Na atoms in your supercell, for instance. Then RKS and UKS give the same results.</div><div><br></div><div>Do you expect any spin polarization in your case?</div><div><br></div><div>Matt<br><br>On Monday, October 17, 2016 at 9:29:35 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, <div><br></div><div>I was running a calculation (calculation 1) with an odd number of electrons. The default multiplicity was set to doublet. This calculation ran without the UKS tag. </div><div><br></div><div>I compared this calculation to one in which I included the UKS keyword in the input (calculation 2). The electronic energies of each calc1 and cal 2 were comparable to 10E-7. </div><div><br></div><div>I found it strange that CP2K would run a system with an odd number of electrons as RKS, is this because of the smearing?</div><div><br></div><div>Input:</div><div><div>&FORCE_EVAL</div><div>  METHOD QS</div><div>  &DFT</div><div>  Â  BASIS_SET_FILE_NAME BASIS_MOLOPT</div><div>  Â  POTENTIAL_FILE_NAME GTH_POTENTIALS</div><div>  Â  &MGRID</div><div>  Â  Â  CUTOFF 500</div><div>  Â  Â  REL_CUTOFF 60</div><div>  Â  Â  NGRIDS 4</div><div>  Â  &END MGRID</div><div>  Â  &QS</div><div>  Â  Â  Â EPS_DEFAULT 1.0E-14</div><div>  Â  Â  Â MAP_CONSISTENT</div><div>  Â  &END QS</div><div>  &SCF</div><div>  Â  Â SCF_GUESS ATOMIC</div><div>  Â  Â EPS_SCF 1.0E-8</div><div>  Â  Â MAX_SCF 500</div><div>  Â  Â ADDED_MOS Â 200</div><div>  Â  Â CHOLESKY INVERSE</div><div>  Â  Â &SMEAR Â ON</div><div>  Â  Â  Â  Â METHOD FERMI_DIRAC</div><div>  Â  Â  Â  Â ELECTRONIC_TEMPERATURE [K] 300</div><div>  Â  Â &END SMEAR</div><div>  Â  Â &DIAGONALIZATION</div><div>  Â  Â  Â  Â  ALGORITHM STANDARD</div><div>  Â  Â &END DIAGONALIZATION</div><div>  Â  Â &MIXING</div><div>  Â  Â  Â  Â  METHOD BROYDEN_MIXING</div><div>  Â  Â  Â  Â  ALPHA Â  0.1</div><div>  Â  Â  Â  Â  BETA Â  Â 1.5</div><div>  Â  Â  Â  Â  NBROYDEN Â 8</div><div>  Â  Â &END MIXING</div><div>  &END SCF</div><div>  Â  &XC</div><div>  Â  Â  &XC_FUNCTIONAL PBE</div><div>  Â  Â  &END XC_FUNCTIONAL</div><div>  Â  &END XC</div><div>CHARGE -1</div><div>  Â  &POISSON</div><div>  Â  Â  POISSON_SOLVER WAVELET</div><div>  Â  Â  PERIODIC NONE</div><div>  Â  &END POISSON</div><div># UKS .TRUE. Â  Â  Â  Â => this is uncommented in the UKS calculation</div><div>  &END DFT</div></div><div><br></div><div><br></div><div>Calculation 1:</div><div><div> DFT| Spin restricted Kohn-Sham (RKS) calculation Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â RKS</div><div> DFT| Multiplicity Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  2</div><div> DFT| Number of spin states Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 1</div><div> DFT| Charge Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â -1</div></div><div><br></div><div><br></div><div>Calculation 2</div><div><div> DFT| Spin unrestricted (spin-polarized) Kohn-Sham calculation Â  Â  Â  Â  Â  Â  Â  UKS</div><div> DFT| Multiplicity Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  2</div><div> DFT| Number of spin states Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 2</div><div> DFT| Charge Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â -1</div></div><div><br></div></div></blockquote></div></div>