| 查看: 1400 | 回复: 0 | ||
vigaryang木虫 (正式写手)
|
[求助]
关于total drift的一点疑惑
|
|
计算某3x3超胞表面，固定底部的两层原子，对顶部的两层原子进行弛豫，计算结束后OUTCAR显示”reached required accuracy - stopping structural energy minimisation“，看起来应该是正常结束了。然而当我在OUTCAR里检索total drift时（最末尾的total drift），发现三个方向上的total drift都很大，感到很困惑。特别地，从表格里数据来看，每个原子在x、y方向上的TOTAL-FORCE都为0，为何这两个方向上的total drift也很大？ POSITION TOTAL-FORCE (eV/Angst) ----------------------------------------------------------------------------------- 1.28478 2.22570 6.26981 0.000000 0.000000 0.000293 3.85503 2.22570 6.26981 0.000000 0.000000 0.000293 -0.00031 4.45154 6.26981 0.000000 0.000000 0.000293 1.28482 6.67744 6.26981 0.000000 0.000000 0.000293 6.42520 2.22570 6.26981 0.000000 0.000000 0.000293 5.14012 4.45154 6.26981 0.000000 0.000000 0.000293 2.56994 4.45154 6.26981 0.000000 0.000000 0.000293 -1.28543 6.67744 6.26981 0.000000 0.000000 0.000293 3.85499 6.67744 6.26981 0.000000 0.000000 0.000293 1.28478 0.74181 4.18448 0.000000 0.000000 0.000635 3.85503 0.74181 4.18448 0.000000 0.000000 0.000635 6.42520 0.74181 4.18448 0.000000 0.000000 0.000635 2.56994 2.96765 4.18448 0.000000 0.000000 0.000635 -0.00031 2.96765 4.18448 0.000000 0.000000 0.000635 1.28486 5.19348 4.18448 0.000000 0.000000 0.000635 3.85503 5.19348 4.18448 0.000000 0.000000 0.000635 -1.28540 5.19348 4.18448 0.000000 0.000000 0.000635 5.14012 2.96765 4.18448 0.000000 0.000000 0.000635 -0.00031 1.48376 2.09852 0.000000 0.000000 -0.046038 2.56994 1.48376 2.09852 0.000000 0.000000 -0.046038 5.14012 1.48376 2.09852 0.000000 0.000000 -0.046038 1.28486 3.70959 2.09852 0.000000 0.000000 -0.046038 -1.28540 3.70959 2.09852 0.000000 0.000000 -0.046038 -0.00023 5.93543 2.09852 0.000000 0.000000 -0.046038 2.56994 5.93543 2.09852 0.000000 0.000000 -0.046038 -2.57048 5.93543 2.09852 0.000000 0.000000 -0.046038 3.85503 3.70959 2.09852 0.000000 0.000000 -0.046038 1.28478 2.22570 0.00000 0.000000 0.000000 0.045109 3.85503 2.22570 0.00000 0.000000 0.000000 0.045109 6.42520 2.22570 0.00000 0.000000 0.000000 0.045109 2.56994 4.45154 0.00000 0.000000 0.000000 0.045109 -0.00031 4.45154 0.00000 0.000000 0.000000 0.045109 1.28482 6.67744 0.00000 0.000000 0.000000 0.045109 3.85499 6.67744 0.00000 0.000000 0.000000 0.045109 -1.28543 6.67744 0.00000 0.000000 0.000000 0.045109 5.14012 4.45154 0.00000 0.000000 0.000000 0.045109 ----------------------------------------------------------------------------------- total drift: 1.286626 0.724770 0.031983 （上表中前面18行原子是free的，后面18行原子是fix的） 查阅资料时，找到了一份国外的ppt，里面有这么一页：  我在计算时所使用的LREAL、PREC、ENCUT与上图最左边这三个菱形框图是一致的（ENCUT比ENMAX的1.3倍还要更大），那为何仍然会得到这么大的total drift。。上述total drift是表面结构优化中的正常现象么？在这种情况下，计算的结果是正确的么？ 诚心求教，谢谢！ 我的INCAR如下： SYSTEM = Cu(111)-P3x3 ISTART = 0 ICHARG = 2 PREC=Accurate LREAL = .F. IBRION = 2 ISIF=2 NSW = 100 POTIM = 0.5 EDIFFG = -1.0E-3 ENCUT = 550 eV ！对Cu bulk进行encut测试时获得 NELM = 60 EDIFF = 1E-5 LCHARG = .F. LWAVE = .F. ISMEAR = 1 SIGMA = 0.07 ！对Cu bulk进行sigma测试时获得 SYMPREC = 1E-3 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
上海大学实验技术岗位非升即走
已经有3人回复
-
有带发论文的吗
已经有4人回复
-
风把牡丹吹跑了
已经有9人回复
-
青C资助名额大幅增加！
已经有15人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
评审感受-评审感受-评审感受
已经有18人回复
-
这年头没有找到涵评专家，还有中面上的可能吗
已经有10人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
精华III:
20