新手初学VASP计算时首先接触到的就是结构优化的计算任务。
在结构优化中,INCAR中的关键参数包括
IBRION ,NSW,ISIF,EDIFF和EDIFFG
各个参数均可在vaspwiki查到可设置的参数以及该参数所具有的设置的含义。
https://www.vasp.at/wiki/index.php/IBRION
https://www.vasp.at/wiki/index.php/NSW
https://www.vasp.at/wiki/index.php/ISIF
https://www.vasp.at/wiki/index.php/EDIFF
https://www.vasp.at/wiki/index.php/EDIFFG
其中,ISIF的可取值范围为;
ISIF = 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8
如在INCAR中没有设置ISIF,则默认设置为等于2,在 IBRION=0 (分子动力学) 或LHFCALC=.TRUE.时则默认设置为0。
不同数值所代表的计算条件为
其中使用较多,也尤为常用的也是ISIF=2的取值,只优化结构中原子位置,不考虑体积变化,适用与二维结构、已明确计算材料的晶格常数和加压或加应变的情况下的结构的计算。
对于体材料,需要弛豫所有方向晶格和原子的话,一般使用ISIF=3。
对于某些需要保留原子间的特殊排列或者有部分原子存在固定(在需要固定的原子坐标最后添加T或F来设置)的情况,可以考虑ISIF=4或6进行优化。
考虑固定晶轴优化的情况下(重新编译vasp的源文件constr_cell_relax.F,并设置OPTCELL,参考:http://blog.wangruixing.cn/2019/05/05/constr),需要确保ISIF的设置是可优化所需要的晶胞的。