Testing POTIM for structure relaxation

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


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
Tomás
Newbie
Newbie
Posts: 17
Joined: Fri Jan 20, 2006 11:49 am

Testing POTIM for structure relaxation

#1 Post by Tomás » Mon Dec 18, 2006 2:30 pm

Hi
I have a monoclinic molecular crystal with C,H,O and Br. Which range of values for parameter POTIM could be more appropiate? I know that each system has its optimum value, but in your experience what values should be tested? I am optimizing the structure of the crystal with ISIF = 3, ISYM=2, EDIFFG = -0.01 and IBRION = 2. I have tested POTIM = 0.15, 0.3, 0.5 AND 0.8 but I do not find out which the best. I only run 10 cycles per value of POTIM, maybe is a very short number of cycles in order to test that parameter. What do you think? Any ideas, hints, clues ...?
Last edited by Tomás on Mon Dec 18, 2006 2:30 pm, edited 1 time in total.
Tomás

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

Testing POTIM for structure relaxation

#2 Post by admin » Mon Dec 18, 2006 3:46 pm

there is no general rule of thumb for how to choose POTIM: it depends on the shape of the potential hypersurface of your system how large it can be chosen. Please mind that for the CG method, POTIM only influences the step width of the first step.
If you are not sure that the system converges to its global minimum, decrease POTIM (even below 0.15) and/or choose an alternative optimization algorithm.
Also, page 25 of optionic.pdf of the vasp-workshop tutorials may help if you are not sure about how to best optimize your system.
Last edited by admin on Mon Dec 18, 2006 3:46 pm, edited 1 time in total.

Post Reply