SSNEB for phase transition in Oxide
Posted: Fri Mar 15, 2013 12:28 am
I have an issue with SSNB to study a martensitic phase transition in a metal oxide and any tips would be very appreciated.
In all cases I did linear interpolation between the first and final states. First I tried to do energy cutoff convergence on 5 intermediate images. The lower energy cuttoffs (<450eV) converge to a saddle point, while the higher ones keep oscillating (every now and then the max force on an atom in an intermediate image shoots up to 50+ ev/A and then go down after few ionic steps).
So I tried these high energy cuttoff on larger number of images (11), and again I observe these oscillations in the forces.
Any suggestions?
The VTST version I use is VTST: version 3.0b, (11/11/12), and here is a typical INCAR:
SSNEB Metal Oxide
ISTART=0 ; ICHARG = 2
ISMEAR = 0 ; SIGMA =0.05
ENCUT= 600
IBRION=3
POTIM=0
IOPT=3 ; ISIF=3; NSW=5000
EDIFF = 1E-06
EDIFFG = -0.01
LWAVE = .FALSE.
LCHARG= .FALSE.
PREC=Accurate
IMAGES=5
SPRING=-5
LCLIMB=.TRUE.
NPAR=1
LNEBCELL=.TRUE.
TIMESTEP=0.04
In all cases I did linear interpolation between the first and final states. First I tried to do energy cutoff convergence on 5 intermediate images. The lower energy cuttoffs (<450eV) converge to a saddle point, while the higher ones keep oscillating (every now and then the max force on an atom in an intermediate image shoots up to 50+ ev/A and then go down after few ionic steps).
So I tried these high energy cuttoff on larger number of images (11), and again I observe these oscillations in the forces.
Any suggestions?
The VTST version I use is VTST: version 3.0b, (11/11/12), and here is a typical INCAR:
SSNEB Metal Oxide
ISTART=0 ; ICHARG = 2
ISMEAR = 0 ; SIGMA =0.05
ENCUT= 600
IBRION=3
POTIM=0
IOPT=3 ; ISIF=3; NSW=5000
EDIFF = 1E-06
EDIFFG = -0.01
LWAVE = .FALSE.
LCHARG= .FALSE.
PREC=Accurate
IMAGES=5
SPRING=-5
LCLIMB=.TRUE.
NPAR=1
LNEBCELL=.TRUE.
TIMESTEP=0.04