Dear PW-Users,<br><br>While performing a test calculation for a non-collinear Fe2 dimer without spin-orbit, I get the following results (just printing the relevant part):<br><br><br>==============================================================================<br>
atom number 1 relative position : 0.1234 0.0000 0.0000<br> charge : 5.500331<br> magnetization : 0.000000 0.000000 2.549039<br> magnetization/charge: 0.000000 0.000000 0.463434<br>
polar coord.: r, theta, phi [deg] : 2.549039 0.000000 360.000000<br><br> ==============================================================================<br><br> ==============================================================================<br>
atom number 2 relative position : 0.0001 0.0000 0.0000<br> charge : 5.500187<br> magnetization : 0.000000 0.000000 2.548980<br> magnetization/charge: 0.000000 0.000000 0.463435<br>
polar coord.: r, theta, phi [deg] : 2.548980 0.000000 360.000000<br><br> ==============================================================================<br><br>(little more output here, then)<br><br> total magnetization = 0.00 0.00 6.00 Bohr mag/cell<br>
absolute magnetization = 6.36 Bohr mag/cell<br><br><br>The cell magnetization prints 6.0 Bohr/cell, which is the expected result. However, the atomic magnetization prints 2.549039 + 2.548980 ~ 5.10 mB.<br><br>
I'm wondering: why is not the atomic magnetization sum equal (or very close to) the cell magnetization? If needed, the input file is attached.<br><br>Thanks in advance.<br><br>Giovani Faccin<br>UFMS / Brazil<br><br><br>
<br><br><br><br><br><br clear="all">Giovani<br>