MAGMOM-tag

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


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
fabella
Newbie
Newbie
Posts: 19
Joined: Fri Nov 11, 2005 4:26 pm
Location: University of South Florida
Contact:

MAGMOM-tag

#1 Post by fabella » Tue Feb 21, 2006 4:00 am

Hello, I have read the manual for MAGMOM-tag. I have 4 atoms; C H S and Ni. I assume that this tag requests the initial magnetic moment for each atom. In this case I put MAGMOM = 0 0 0 28*0.6*1.3 (1.3 for scaling of expt. val. Ni mag.mom.). So, Is this correct? When trying to run I receive error reading 'MAGMOM' from INCAR. Error code was IERR=0 ... . Found N= 4 data.

Question: What is the proper format of this tag?

A clearer explanation in VASP UG may be needed. (?)

Thank You

W
Last edited by fabella on Tue Feb 21, 2006 4:00 am, edited 1 time in total.

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

MAGMOM-tag

#2 Post by admin » Tue Feb 21, 2006 10:45 am

28*0.6*1.3 is not a valid format to read a magnetic moment.
for 3 non-magnetic and 28 magnetic atoms,
MAGMOM can be given like
3*0 28*0.9
0 0 0 28*0.9
or you can give each moment explicitely

Please mind that the MAGMOM specified in INCAR is just a starting guess, the actual self-consistent magnetic moments will be calculated as the electronic structure converges.
Last edited by admin on Tue Feb 21, 2006 10:45 am, edited 1 time in total.

fabella
Newbie
Newbie
Posts: 19
Joined: Fri Nov 11, 2005 4:26 pm
Location: University of South Florida
Contact:

MAGMOM-tag

#3 Post by fabella » Wed Feb 22, 2006 3:29 am

TY for the info.

Now, it seems I have a problem with convergence (Continuation job ISTART=1, ICHARG=1, i have WAVECAR and CHGCAR. I can guess that the problem started when adding the ISPIN=2 (spin-polarization). I ran two calculations one with MAGMOM set appropriate, the other without MAGMOM.

-The first case does not even finish the first sc loop (-c gives 35).

-The second case gets through the first sc loop then mid-way through the second 'stops' without notification ... just "hangs" there.

One thing in common, they both are sitting 'about' the same energy (within 1eV of each other). "They are screaming to converge"

How do I successfully troubleshoot "solve" this issue ... I was thinking of a couple things. 1) Maybe I should rerun previous relaxation of structure with strict force criteria ( <0.05 eV/A instead of <0.08 eV/A ), then run with spin polarization w/o relaxation. 2) I think my initial structure is close to equilibrium so I use IBRION=1, maybe I should use IBRION=2 (much slower) to relax and find equilibrium positions ... then turn on spin-polarization. (?) I think this with some uncertainty and time constraints. Could you be of any assistance here, maybe this situation is quite common if the structure is not 'well' relaxed before switching on ISPIN=2.

THANKX

W
Last edited by fabella on Wed Feb 22, 2006 3:29 am, edited 1 time in total.

Post Reply