Force acting on one of the images in CL-NEB is huge
Posted: Sun Aug 21, 2011 10:32 am
Hi dear all,
My nebef.pl results show this after 200 electronic cycle:
0 0.00000000 -251.49423800 0.00000000
1 0.08026200 -251.13039000 0.36384800
2 0.08919400 -250.65757900 0.83665900
3 0.16627700 -251.13418700 0.36005100
4 5.31365100 -249.68733300 1.80690500
5 0.00000000 -251.26829100 0.22594700
I dont know why the force acting on image 4, is so huge.
This is my INCAR:
ISTART = 0
ICHARG = 2
ISPIN = 2
MAGMOM = ...
ENCUT = 450
ISMEAR = 1
SIGMA = 0.1
NELM = 100
EDIFF = 1E-4
ALGO = F
ICHAIN = 0
IMAGES = 4
SPRING = -5
LSCALAPACK = .FALSE.
EFIRST = -251.48
ELAST = -251.26
GGA = PE
VOSKOWN = 1
NBANDS = 230
LREAL = .TRUE.
ADDGRID = .TRUE.
PREC = Medium
NSW = 550
IBRION = 3
POTIM = 0.15
EDIFFG = -0.05
LMAXMIX = 4
ISYM = 0
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = ...
LDAUU = ...
LDAUJ = ...
LDAUPRINT = 2
I appreciate any help,
Regards
My nebef.pl results show this after 200 electronic cycle:
0 0.00000000 -251.49423800 0.00000000
1 0.08026200 -251.13039000 0.36384800
2 0.08919400 -250.65757900 0.83665900
3 0.16627700 -251.13418700 0.36005100
4 5.31365100 -249.68733300 1.80690500
5 0.00000000 -251.26829100 0.22594700
I dont know why the force acting on image 4, is so huge.
This is my INCAR:
ISTART = 0
ICHARG = 2
ISPIN = 2
MAGMOM = ...
ENCUT = 450
ISMEAR = 1
SIGMA = 0.1
NELM = 100
EDIFF = 1E-4
ALGO = F
ICHAIN = 0
IMAGES = 4
SPRING = -5
LSCALAPACK = .FALSE.
EFIRST = -251.48
ELAST = -251.26
GGA = PE
VOSKOWN = 1
NBANDS = 230
LREAL = .TRUE.
ADDGRID = .TRUE.
PREC = Medium
NSW = 550
IBRION = 3
POTIM = 0.15
EDIFFG = -0.05
LMAXMIX = 4
ISYM = 0
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = ...
LDAUU = ...
LDAUJ = ...
LDAUPRINT = 2
I appreciate any help,
Regards