Hello Prof. Henkelman and group
I am trying to find the TS of a hydride shift step on Al2O3. Initially, I started a NEB+climb calculation using a single kpoint, iopt=3, default time step of 0.1. This calculation showed issues with force convergence; the forces were oscillating between 1 and 0.5 eV/A. To avoid this problem, I lowered the time step to 0.01. Even after trying this, the forces do not converge, and then eventually I see the atoms showing some meaningless movements (CH bonds dissociating randomly).
I am stuck at this point. It would be of great help if you can provide some suggestions regarding how to troubleshoot this NEB.
Also, I have a question, does the TIMESTEP parameter of an NEB calculation serve the same purpose of POTIM ? Is are the units of TIMESTEP in angstrom?
Attached below is the zipped folder of the NEB with all necessary files. Kindly have a look.
NEB calculation of hydride shift step not converging
Moderator: moderators
NEB calculation of hydride shift step not converging
 Attachments

 NEB_climb_new_idpp_4.zip
 (85.91 MiB) Downloaded 2233 times