Get mag.sh

go back to Main Page, Computational Resources, Scripts, Scripts for VASP or VASP

File:Get mag.tgz 

This script will do the following:

1) Reads CONTCAR and calculate the total atom numbers;

2) Extracts the magnetization information from the OUTCAR to a file named mag-temp ;

3) Prints the magnetization of the atoms you want to check.

Note 1:

To get the magnetization information, we need set:

ISPIN = 2 MAGMOM = xxx (It depends on your system) LORBIT = 11

Note 2:

The total magnetization in OUTCAR is less than the value from OSZICAR. This is because:

The difference in the values roots in the fact that, when you project the magnetization on the atoms,

you use a default value (RWIGS) to define the radius of each atom.

For this reason, even if you compute the magnetization for all the atoms,

the sum of their "volume" does not fill up the total volume of the unit cell,

and so the sum of the projections of the magnetization on each atom is different from the overall magnetization, written in the OSZICAR.

PS: I copied the explanations above from http://cms.mpi.univie.ac.at/vasp-forum/viewtopic.php?t=15642

To use this script:

1) Get one atom's magnetization: for example, atom 15:

get_mag.sh 15

2) Get many atoms' magnetization, for example, atoms: 1 3 5 7 9

for i in 1 3 5 7 9  ; do get_mag.sh $i ; done

3) You can also save the atom list in a file (atom-list): and do the following:

for i in $(cat atom-list) ; do get_mag.sh $i ; done

Enjoy!