neb of spontaneous reactions

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

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

neb of spontaneous reactions

#1 Post by mhyman » Sun Aug 14, 2005 2:26 am

Will a neb calculation of a spontaneous reaction pathway converge in VASP?

Thanks,
Matt
Last edited by mhyman on Sun Aug 14, 2005 2:26 am, edited 1 time in total.

tjf
Full Member
Full Member
Posts: 107
Joined: Wed Aug 10, 2005 1:30 pm
Location: Leiden, Netherlands

neb of spontaneous reactions

#2 Post by tjf » Thu Aug 18, 2005 2:39 pm

What exactly do you mean by "spontaneous reaction pathway"?
Last edited by tjf on Thu Aug 18, 2005 2:39 pm, edited 1 time in total.

mhyman
Newbie
Newbie
Posts: 19
Joined: Wed Jul 20, 2005 10:37 pm

neb of spontaneous reactions

#3 Post by mhyman » Thu Aug 18, 2005 7:34 pm

One in which there is no activation barrier. The starting and end points of the pathway were extremely close to each other. Since the first post I ran single point calculations of three geometries linearly interpolated between image 00 and 01, which reveal an energy barrier of 0.02 eV.
Last edited by mhyman on Thu Aug 18, 2005 7:34 pm, edited 1 time in total.

tjf
Full Member
Full Member
Posts: 107
Joined: Wed Aug 10, 2005 1:30 pm
Location: Leiden, Netherlands

neb of spontaneous reactions

#4 Post by tjf » Fri Aug 26, 2005 8:59 am

Well, I was going to say that if there's no barrier how is the reactant defined? The only option is for a bimolecular reaction. If there's no reaction, there's no chemical species.

Have you simply tried a NEB calc on your system?
Last edited by tjf on Fri Aug 26, 2005 8:59 am, edited 1 time in total.

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

neb of spontaneous reactions

#5 Post by admin » Fri Sep 16, 2005 12:25 pm

If there is no reaction barrier, you will probably not be able to find
any minimum for the 'initial state', because the system relaxes to
the 'final state' immediately, as it is the ground state structure.
So you have to choose the starting geometry according to some other criterium.

The NEB method should work though, because the relaxation of the
geometries only is done on the hyper-planes normal to the reaction
path. There only has to be taken care that the entire hyper-plane
does not drift away, which can be done eg. by choosing a rather small the number of ionic steps (NSW)
Last edited by admin on Fri Sep 16, 2005 12:25 pm, edited 1 time in total.

Post Reply