NEB Energy linearly extrapolates between initial and final
Posted: Tue May 31, 2022 10:04 pm
Hello,
I've been attempting to run an NEB calculation using the attached files. Below are some of the suggestions that I've taken from here to attempt to get relative energies that are higher than the initial and final relative energies (the expected trend according to literature for H2S dissociation on a Cu2O surface):
1) I've attempted the NEB method by both the standard VASP package and the energies trend towards or converge to energies that fall between the initial and final energies...I can see that it's linearly decreasing from the initial to the final for each image.
2) I've attempted the NEB method by using the IOPT 2 and 3 optimizers and it does not converge. I've also looked at the forces that correspond with the energies. The lower the force (i.e. the closer it gets to the EDIFFG setting), the TS energy again falls below the initial.
3) I've kept my ENCUT and KPOINTS grid consistent with the initial and final configuration settings.
4) I've added ISYM to initial, final, and transition state calculations - same issue occurs.
5) It seems as if I'm having a system migration, so I fixed elements that were not involved in the dissociation of the hydrogen element - same issue occurs.
6) I've attempted to add more intermediate images between - same issue occurs.
I know there are different force based optimizers available with the vtst package, but shouldn't the overall trend still give a maximum saddle point regardless of the the force based optimizer used?
I've attached in progress calculations for a standard VASP package NEB run and a vtst package NEB run with IOPT 2. Any and all feedback/help is greatly appreciate.
Also a side note: I terminated the jobs once I saw that the energy would pass the initial energy as I did not want to keep using up resources when I knew the energy would end up converging to a lower energy value.
I've been attempting to run an NEB calculation using the attached files. Below are some of the suggestions that I've taken from here to attempt to get relative energies that are higher than the initial and final relative energies (the expected trend according to literature for H2S dissociation on a Cu2O surface):
1) I've attempted the NEB method by both the standard VASP package and the energies trend towards or converge to energies that fall between the initial and final energies...I can see that it's linearly decreasing from the initial to the final for each image.
2) I've attempted the NEB method by using the IOPT 2 and 3 optimizers and it does not converge. I've also looked at the forces that correspond with the energies. The lower the force (i.e. the closer it gets to the EDIFFG setting), the TS energy again falls below the initial.
3) I've kept my ENCUT and KPOINTS grid consistent with the initial and final configuration settings.
4) I've added ISYM to initial, final, and transition state calculations - same issue occurs.
5) It seems as if I'm having a system migration, so I fixed elements that were not involved in the dissociation of the hydrogen element - same issue occurs.
6) I've attempted to add more intermediate images between - same issue occurs.
I know there are different force based optimizers available with the vtst package, but shouldn't the overall trend still give a maximum saddle point regardless of the the force based optimizer used?
I've attached in progress calculations for a standard VASP package NEB run and a vtst package NEB run with IOPT 2. Any and all feedback/help is greatly appreciate.
Also a side note: I terminated the jobs once I saw that the energy would pass the initial energy as I did not want to keep using up resources when I knew the energy would end up converging to a lower energy value.