| 查看: 4698 | 回复: 34 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
VASP静态计算中体系不收敛，恳请各位前辈指教 已有1人参与
|
|||
|
体系含32个原子，从木虫上搜索过各种处理方法，经修改测试，结果还是不收敛。 INCAR文件如下： SYSTEM = ******** ISTART = 0 ICHARG = 2 PREC = Accurate LREAL = .F. ENCUT = 220 ISPIN = 2 ISYM = 0 NELM = 200 EDIFF = 1E-04 ISMEAR = 1 SIGMA = 0.5 IBRION = -1 NSW = 300 POTIM = 0.5 EDIFFG = -0.05 ALGO = Fast NELMDL = -5 在上面参数中，尝试修改过NELM,ISMEAR,SIGMA,ALGO,NELMDL等参数设置，均未能达到收敛。ENCUT设置为220，是因为现在在做前期的测试。K点为3x3x3。POSCAR由驰豫收敛后的CONTCAR得到，但驰豫是使用的ENCUT为370（小虫想这里应该没有问题）。 在以上述INCAR参数进行计算时，使用grep 'energy without entropy' OUTCAR命令，抓到的结果如下图。恳请前辈指教。如需要其他文件，小虫继续上传。  grep energy 截图.png |
回复此楼» 猜你喜欢
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
-
A期刊撤稿
已经有5人回复
-
临港实验室与上科大联培博士招生1名
已经有8人回复
-
26申博自荐
已经有7人回复
-
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有4人回复
-
带资进组求博导收留
已经有9人回复
-
求助大佬们，伤口沾上了乙腈
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有9人回复
|
看OUTCAR得情况，应该就不是能量没有收敛了；剩下的就是力收敛问题。今天做过测试，分别设置如下参数：EDIFF=1E-03；EDIFFG=-0.1。结果还是未能收敛。然后我请师兄给我发了一个他做的可以收敛的算例，先正在测试，其中参数主要变化有：NSW采用默认(0)，IBRION采用默认（-1），ISIF采用默认（2），NELM=100；力收敛采用的是EDIFFG=-0.05，EDIFF=1E-04。 很不好的是，在编辑给您的这个回帖的时候，结果出来了 ，虽然OUTCAR已给出了‘Total cpu time used (sec): ’，但是并没有提示‘reach the accuracy’。 之所以在后面的计算中采用了ENCUT=650的设置，是因为昨天晚上经教导，说用大些的ENCUT来算，说ENCUT太小的话可能一辈子都不会收敛。（虽然觉得个人ENCUT不会影响是否收敛，但是为了解决问题，还是采用了尽可能大的ENCUT） 前辈，下一步，我想再刚才的设置的基础上再测试一下EDIFF和EDIFFG，如果还不能收敛，我真的找不到怎样再继续测试以达到收敛了。请前辈不吝指教。 |
16楼2016-09-06 21:00:29
2楼2016-09-05 07:50:12
★
蛋蛋小童鞋: 金币+1, 谢谢参与 2016-09-05 11:49:24
蛋蛋小童鞋: 金币+1, 谢谢参与 2016-09-05 11:49:24
3楼2016-09-05 08:22:41
4楼2016-09-05 10:26:15