DFT新手教程：VASP中ISIF取值设置

    新手初学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的设置是可优化所需要的晶胞的。