Units of magnetic moment on CHG file

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
salva
Newbie
Newbie
Posts: 11
Joined: Thu Oct 12, 2006 6:57 pm
License Nr.: 592
Contact:

Units of magnetic moment on CHG file

#1 Post by salva » Tue Dec 26, 2006 6:46 pm

Dear users/developers.
I am working with a magnetic system, having Mn atoms on it. I wrote a script to get the magnetic moment within a sphere centered on the Mn atoms. Specifics follow:

The calculation is spin-polarized. I take the CHG file and read the density, dividing each entry by the volume of the cell.

The radii around each Mn atom is about 1.2 Ang. The calculation is done with PAW pseudos.

After adding up the entries from the second dataset on the CHG file (related to the magnetization) around the radius mentioned about the selected ions, I get for the magnetization some big values, say, ~14000 for the magnetization about this ion.

I am confident about the script as also the electronic charges around the ions are overestimated by a factor of thousands. So even though the values are high, they seem to be consistent within each other (consistently bigger I should say).

The value should be close to three, as this is the value assigned through the MAGMOM flag. So, even though I can just "normalize" the value I get so that it becomes the value I am expecting to get, in this case 3, I'd rather learn about the units on the file so that I can make the necessary change of units that (hopefully) will give me a consistent answer.

Thank you for your time.

Salvador
Last edited by salva on Tue Dec 26, 2006 6:46 pm, edited 1 time in total.

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

Units of magnetic moment on CHG file

#2 Post by admin » Wed Dec 27, 2006 1:42 pm

why do you not make use of the LORBIT input tag? It will give the local charges and magnetisations (spin-up - spin-down charge) in OUTCAR.
please mind that
if you set LORBIT >= 10, the sphere radius will be chosen to be RWIGS as read from POTCAR
to set RWIGS explicitely, please define it in INCAR (for each atom) and set LORBIT=1
Last edited by admin on Wed Dec 27, 2006 1:42 pm, edited 1 time in total.

Post Reply