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