NEB calculations

Question on input files/tags, interpreting output, etc.

Please check whether the answer to your question is given in the VASP online manual or has been discussed in this forum previously!

Moderators: Global Moderator, Moderator

Post Reply
Message
Author
mhyman
Newbie
Newbie
Posts: 19
Joined: Wed Jul 20, 2005 10:37 pm

NEB calculations

#1 Post by mhyman » Wed Jul 20, 2005 10:56 pm

Hi,

I've been trying to perform NEB calculations for my system. I've gotten convergence using 1 image and am now trying to expand to 3. I'm at the point where the energies of the image files oscillate with a difference between high and low energies of 0.006 eV, so it doesn't seem like the calculations will ever converge. The only thing that I can think that might be a problem is that images 01 and 02 are less than 0.01 eV less stable that image 00. That is, the energy barrier from 00 is basically nonexistent. Also, there are numerous local minima with similar geometries, but I'm fairly certain that my two endpoints are true minima. So my question is...should this converge (is there any reason why it shouldn't)? If not, how do I handle this system to calculate a transition state energy?

Thanks,
Matt

The pertinent settings are

SPRING = -5.0
IBRION = 1
POTIM = .1
NFREE = 2
NELMIN = 4
EDIFF = 1e-5
Last edited by mhyman on Wed Jul 20, 2005 10:56 pm, edited 1 time in total.

graeme
Jr. Member
Jr. Member
Posts: 54
Joined: Sat Jul 02, 2005 6:11 am
Location: Austin, TX
Contact:

NEB calculations

#2 Post by graeme » Thu Jul 21, 2005 7:24 am

Try IBRION=3. This routine does not rely on estimating second derivaties, so it has a larger radius of convergence and tends to be more stable (and less efficient). Also, there is no Hessian consistent with the NEB forces, and I don't think it's well understood how the quasi-Newton methods deal with non-conservative forces.
Last edited by graeme on Thu Jul 21, 2005 7:24 am, edited 1 time in total.

Post Reply