<html>
  <head>
    <meta content="text/html; charset=utf-8" http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    Dear Martin,<br>
    <br>
    you may have a look at this keyword in the &SMEAR section<br>
    <br>
    <a
href="http://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/SCF/SMEAR.html#list_FIXED_MAGNETIC_MOMENT">http://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/SCF/SMEAR.html#list_FIXED_MAGNETIC_MOMENT
    </a><br>
    <div class="moz-cite-prefix"><br>
      I guess it imposes the constraint you are looking for.<br>
      <br>
      Best regards,<br>
      <br>
      Matthias<br>
      <br>
      On 17.02.2015 20:50, M. Brehm wrote:<br>
    </div>
    <blockquote
      cite="mid:3ff69495-0b0b-44...@googlegroups.com"
      type="cite">
      <meta http-equiv="Content-Type" content="text/html; charset=utf-8">
      <div dir="ltr">Dear Community,<br>
        <br>
        I am performing DFT calculations of a solid-state lattice with
        CP2k, using the classical mixing scheme and
        Fermi-Dirac-Smearing, like usual for lattices. As the total
        number of electrons in my system is odd, I perform unrestricted
        DFT calculations, setting the spin multiplicity to 2. In the
        beginning, CP2k displays 288 alpha electrons and 287 beta
        electrons, like it should be. The SCF cycle converges, and in
        the end, I print out the MO energies and occupation numbers.
        Suddenly, there are 287.5 alpha electrons and 287.5 beta
        electrons! For my understanding, the total number of alpha/beta
        electrons is not a degree of freedom of a unrestricted DFT
        calculation, and should remain fixed during the calculations.
        How can it be that in the end I get these fractional electron
        numbers? With slightly different systems or geometries, this is
        not the case, and the total electron count for alpha/beta
        remains constant.<br>
        <br>
        I assumed that it could have something to do with the initial
        SCF guess (poor guess might lead to convergence into a local
        minimum which is not the ground state), but I re-ran the
        calculation with SCF guess "random", and still get exactly the
        same results.<br>
        <br>
        My input file and output file are attached to this posting.<br>
        <br>
        Any hints/comments are appreciated.<br>
        <br>
        Best regards,<br>
        Martin<br>
      </div>
      -- <br>
      You received this message because you are subscribed to the Google
      Groups "cp2k" group.<br>
      To unsubscribe from this group and stop receiving emails from it,
      send an email to <a moz-do-not-send="true"
        href="mailto:cp2k+uns...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
      To post to this group, send email to <a moz-do-not-send="true"
        href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
      Visit this group at <a moz-do-not-send="true"
        href="http://groups.google.com/group/cp2k">http://groups.google.com/group/cp2k</a>.<br>
      For more options, visit <a moz-do-not-send="true"
        href="https://groups.google.com/d/optout">https://groups.google.com/d/optout</a>.<br>
    </blockquote>
    <br>
  </body>
</html>