| 查看: 1428 | 回复: 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 |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
100