NEB calculation between previously converged images.
Posted: Thu Sep 21, 2023 8:32 pm
Dear all,
I have been using the NEB method for a while now, and I am wondering what are the pro's and con's to run NEB calculations between images of previous runs (provided these images are converged). My question arise from the need of a more fine sampling of the reaction path only around a certain region. I am aware of the existence of other methods that could solve this problem in a different manner, yet I would like to find out how valid it is to do so using the NEB/cNEB method. As an example I made an attempt to sample the reaction path of a radical attack on a slab (~200 atoms), that led to the energy profile (in attachment).
A first 5-images NEB run was performed between initial and final states as usual, then two 5-images additional ones between images 1-2 and 2-3, then three additional one between images 1-2, 2-3 and 3-4 of the "sub" NEB 2-3 previously mentioned. Eventually I found what I was looking for (i.e. the "bump" you see on the energy profile) and ran a cNEB to accurately determine the position and energy of this local maximum. Every NEB calculation was performed with EDIFF=5E-09 and EDIFFG=-1E-03, but endpoints of "sub" NEB were not relaxed and taken as they were computed on the previous NEB calculations (but of course endpoints of the initial NEB calculation were converged using the same criteria).
An analysis of the coordinates of the hydrogen atoms involved in this reaction shows that the reaction is well defined, with a smooth displacement of the atoms along the reaction path as one would expect for this particular reaction.
Are there concerns that should be raised with respect to the methodology I used ?
Best regards,
Emerick
PS: I am sorry if this is a duplicate topic, if it is could you refer me to the original question ?
I have been using the NEB method for a while now, and I am wondering what are the pro's and con's to run NEB calculations between images of previous runs (provided these images are converged). My question arise from the need of a more fine sampling of the reaction path only around a certain region. I am aware of the existence of other methods that could solve this problem in a different manner, yet I would like to find out how valid it is to do so using the NEB/cNEB method. As an example I made an attempt to sample the reaction path of a radical attack on a slab (~200 atoms), that led to the energy profile (in attachment).
A first 5-images NEB run was performed between initial and final states as usual, then two 5-images additional ones between images 1-2 and 2-3, then three additional one between images 1-2, 2-3 and 3-4 of the "sub" NEB 2-3 previously mentioned. Eventually I found what I was looking for (i.e. the "bump" you see on the energy profile) and ran a cNEB to accurately determine the position and energy of this local maximum. Every NEB calculation was performed with EDIFF=5E-09 and EDIFFG=-1E-03, but endpoints of "sub" NEB were not relaxed and taken as they were computed on the previous NEB calculations (but of course endpoints of the initial NEB calculation were converged using the same criteria).
An analysis of the coordinates of the hydrogen atoms involved in this reaction shows that the reaction is well defined, with a smooth displacement of the atoms along the reaction path as one would expect for this particular reaction.
Are there concerns that should be raised with respect to the methodology I used ?
Best regards,
Emerick
PS: I am sorry if this is a duplicate topic, if it is could you refer me to the original question ?