24小时热门版块排行榜

返回列表

vigaryang

木虫 (正式写手)

应助: 16 (小学生)
金币: 3466.3
散金: 677
红花: 6
帖子: 471
在线: 495.7小时
虫号: 227097
注册: 2006-03-23
性别: GG
专业: 应用化学

[求助] 关于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

回复此楼

» 猜你喜欢

1楼 2016-08-24 00:43:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 vigaryang 的主题更新

返回列表