| 查看: 1444 | 回复: 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一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
30