应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于total drift的一点疑惑
11/1返回列表
查看: 1390  |  回复: 0
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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，里面有这么一页：
关于total drift的一点疑惑

我在计算时所使用的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
回复此楼

» 猜你喜欢
1楼 2016-08-24 00:43:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 vigaryang 的主题更新
11/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工调剂 +7 A-哆啦Z梦 2026-03-23 12/600 2026-03-23 23:16 by 星空星月
[考研] 材料专业求调剂 +11 hanamiko 2026-03-18 11/550 2026-03-23 23:12 by peike
[考研] 一志愿南航材料专317分求调剂 +3 炸呀炸呀炸薯条 2026-03-23 3/150 2026-03-23 20:47 by pswait
[考研] 303求调剂 +4 元夕元 2026-03-20 4/200 2026-03-23 19:00 by macy2011
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 333求调剂 +6 87639 2026-03-21 10/500 2026-03-23 10:41 by Iveryant
[考研] 一志愿东华大学化学070300，求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 317求调剂 +12 申子申申 2026-03-19 18/900 2026-03-22 22:23 by luoyongfeng
[考研] 311求调剂 +6 冬十三 2026-03-18 6/300 2026-03-22 20:18 by edmund7
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 考研调剂 +4 来好运来来来 2026-03-21 4/200 2026-03-22 12:15 by 星空星月
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[基金申请] 学校已经提交到NSFC，还能修改吗？ 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
信息提示
关闭
请填处理意见
关闭