vasp优化4*4石墨烯报错

小书童

vasp优化4*4石墨烯报错，INCAR和K-POINTs以及报错如图，请各位老师帮帮忙看看，哪里的问题。

卡开发发

看看ISYM=0是否还报错，如果不报错的话检查一下结构对称性或者修改一下对称性识别的阈值。

另外：
1、石墨烯是二维材料，ISIF=3很可能stress是收敛不了的，除非你的stress计算方式是修改过的。建议你可以采用固定体积的计算或者ISIF=2索性不要优化xy方向了。
2、POTIM=0.1太小了，如果结构优化过程能量和力的振荡不是很严重，其实并不需要这么小，默认值就很好。
3、EDIFF = 1e-4是不是还是太大了些，想想Gaussian的SCF收敛标准要严格得多啊。
4、如果只是石墨烯的话，不需要自旋极化计算，关掉自旋极化能省不少的时间。

小书童

多谢老师回复！
ISYM=0不再报错，POTIM修改为默认的0.5.
EDIFF=1e-4是vasp的默认设置，提升为10-5，对应的EDIFF提升为10-4

对于您说了“检查一下结构对称性或者修改一下对称性识别的阈值”，老师，可以再展开说下吗？

卡开发发

小书童 发表于 2017-2-19 15:55
多谢老师回复！
ISYM=0不再报错，POTIM修改为默认的0.5.
EDIFF=1e-4是vasp的默认设置，提升为10-5，对应 ...

4x4的石墨烯肯定有对称性，本身可以利用之进行加速计算，但是现在不知道对称性识别过程出了什么样的bug，你可以先从自己的结构检查，比如导入会MS找一下对称性等，或者索性修改一下VASP对 对称性识别的标准。

EDIFFG其实1E-4也不大，这是能量收敛的标准，我个人还是觉得力收敛的标准让人放心一些。

小书童

卡开发发 发表于 2017-2-19 16:22
4x4的石墨烯肯定有对称性，本身可以利用之进行加速计算，但是现在不知道对称性识别过程出了什么样的bug， ...

谢谢老师的回复，我利用MS找了下对称性，并导入到结构，用VESTA.exe导入vasp计算文件的时候，发现变成了单胞了，我想要4*4的超胞，导入过程如图

卡开发发

小书童 发表于 2017-2-19 16:46
谢谢老师的回复，我利用MS找了下对称性，并导入到结构，用VESTA.exe导入vasp计算文件的时候，发现变成了 ...

没有，只是让你导入到MS看看能不能识别出对称性，如果能识别出对称性我就不知道是什么因素导致的。目前我猜测是密度截断不够导致的（换句话说，就是密度的网格不符合对称性），要不然看看PREC=A看看怎样，再不行就将就着ISYM=0吧。

helpme

感谢 @卡开发发 对新手的耐心指导。
跟着一起学习。

小书童

卡开发发 发表于 2017-2-19 12:54
看看ISYM=0是否还报错，如果不报错的话检查一下结构对称性或者修改一下对称性识别的阈值。

另外：

今天读manual才明白老师您的意思，更改为EDIFF  = 0.1E-04，EDIFFG = -0.01这样可以了吧，还有一个疑问，就是能量和力这两个参数都用EDIFFG 来设定，就是说能量和力只能选择一个作为收敛不准，不能想Gaussian两个都设定，对吗？