Dear Professor Henkelman,
I'm studying the diffusion of a Li atom on the GeP monolayer using CINEB method, but I have found a strange problem. I'm using the vtstcode and I can see the related parts in OUTCAR file for that but I can not see any difference between the results obtained from NEB and CINEB calculations. In other words, it seems that CINEB is not working in my calculations. I have no idea if the problems is coming form vtstcode installation or I'm doing something wrong. Below you can see the part for vtstcode in my OUTCAR file.
VTST: version 3.1, (03/28/14)
CHAIN: initializing optimizer
OPT: Using VASP QUASInewton optimizer
CHAIN: Read ICHAIN 0
CHAIN: Running the NEB
NEB: SPRING 5.000000
NEB: LCLIMB T
NEB: LTANGENTOLD F
NEB: LDNEB F
NEB: LDNEBORG F
NEB: LNEBCELL F
NEB: EFIRST 0.000000
NEB: ELAST 0.000000
I have also attached the MEP profile created by nebspline.pl, for a CINEB calculation that shows the Image 2 is not located at the expected position.
I run CINEB calculation ans follows, Please correct me if I am doing something wrong. First I run a regular NEB calculation. After convergence, I move CONTCAR file of each 0N directories to POSCAR file. Finally I add LCLIMB = .TRUE. to the INCAR and rerun the job.
I will be grateful if you help me solve this problem.
CINEB problem
Moderator: moderators
CINEB problem
 Attachments

 mep.jpg (588.45 KiB) Viewed 8411 times
Re: CINEB problem
The energy path looks reasonable.
One comment: you can run the climbing image right from the start. This will be more efficient than converging a regular NEB and then turning on the climbing image.
One thing that I can not tell from what you have posted is how well converged the band is and what your convergence criterion is. If you use EDIFFG=0.002, the highest energy image should converge very close to the saddle (top of the MEP). Check the residual force, to see if the climbing image is still converging or if you need to tighten your convergence criterion.
One comment: you can run the climbing image right from the start. This will be more efficient than converging a regular NEB and then turning on the climbing image.
One thing that I can not tell from what you have posted is how well converged the band is and what your convergence criterion is. If you use EDIFFG=0.002, the highest energy image should converge very close to the saddle (top of the MEP). Check the residual force, to see if the climbing image is still converging or if you need to tighten your convergence criterion.
Re: CINEB problem
So many thanks for your prompt answer. I used EDIFFG=0.02 in my calculations. Just one more question: Does CINEB work well with PBED2 or PBED3?
Re: CINEB problem
it should be fine.