NEB criterion, forces

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


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
neogien

NEB criterion, forces

#1 Post by neogien » Tue Apr 15, 2008 4:39 am

Hello VASP users and admin,

I did several NEB calculations and found some TSs and reaction paths I wonder however my TSs are resonable.

In my knowledge for normal geometry optimization calculations EDIFFG is the criterion. And for the NEB I should judge whether my TS is resonable with forces in the OUTCAR.
Is there any resonable guide to force values?
1eV/Ang is enough? Or should it be smaller than 0.x eV/Ang?


One more thing, sometimes my NEB jobs finish just after first ionic step even though the forces on the atoms are still high (2~3eV/Ang). I'm using IBRION=1 and default POTIM (0.5).
Do anybody can explain what the possible reason for this abnormal situation is? Shoud I change optimizer IBRION=3 with small POTIM (0.01)?

Many Thanks

HY KIM
Last edited by neogien on Tue Apr 15, 2008 4:39 am, edited 1 time in total.

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

NEB criterion, forces

#2 Post by admin » Thu Apr 24, 2008 10:37 am

please note that in NEB (as in the standard implementation of VASP) the atoms only move along the hypertangents normal to the reaction path. If all relaxations were allowed, the system would relax back to the ground state. I would rather choose the change in total energies (EDIFFG>0) as convergence criterium.
Last edited by admin on Thu Apr 24, 2008 10:37 am, edited 1 time in total.

Post Reply