<div dir="ltr"><div>I am attempting to carry out a calculation for ferromagnetic bcc Fe using CP2K 6.1. Attached are the input and output files.</div><div><br></div><div>The calculation runs to completion, but I am confused by the output. The multiplicity is not 1</div><div><br></div><div> <span style="font-family: courier new, monospace;">Orbital energies  State     Spin  L     Occupation   Energy[a.u.]    Energy[eV]<br><br>                       1    alpha  0          1.000      -3.533760    -96.158487<br>                       2    alpha  0          1.000      -0.171394     -4.663878<br>                       1     beta  0          1.000      -3.422092    -93.119858<br>                       2     beta  0          1.000      -0.151638     -4.126285<br><br>                       1    alpha  1          3.000      -2.201116    -59.895424<br>                       1     beta  1          3.000      -2.092954    -56.952177<br><br>                       1    alpha  2          4.250      -0.282066     -7.675405<br>                       1     beta  2          1.750      -0.191730     -5.217225<br></span><br></div><div>yet the system appears to have zero spin moment:<br></div><div><br></div><div><span style="font-family: courier new, monospace;"> !-----------------------------------------------------------------------------!<br>                     Mulliken Population Analysis<br><br> #  Atom  Element  Kind  Atomic population (alpha,beta) Net charge  Spin moment<br>       1     Fe       1         8.000000     8.000000     0.000000     0.000000<br>       2     Fe       1         8.000000     8.000000     0.000000    -0.000000<br> # Total charge and spin       16.000000    16.000000     0.000000    -0.000000<br><br> !-----------------------------------------------------------------------------!</span><br></div><div><br></div><div>Between these blocks of output seems to be a rescaling step which might explain the discrepancy:<br></div><div><br></div><div><span style="font-family: courier new, monospace;">Spin 1<br> Re-scaling the density matrix to get the right number of electrons for spin 1<br>                  # Electrons              Trace(P)               Scaling factor<br>                           16                18.500                        0.865<br><br> Spin 2<br> Re-scaling the density matrix to get the right number of electrons for spin 2<br>                  # Electrons              Trace(P)               Scaling factor<br>                           16                13.500                        1.185</span><br><br></div><div>Can anyone explain what is happening here and how I allow the system to be magnetic? Is there a keyword I need to add?<br></div><div><br></div></div>