24小时热门版块排行榜

返回列表

saikun

木虫 (正式写手)

应助: 4 (幼儿园)
金币: 2483.2
散金: 891
红花: 5
帖子: 424
在线: 89.1小时
虫号: 197920
注册: 2006-02-25
性别: GG
专业: 古陶瓷与传统陶瓷

[求助] 关于优化结构中的 total drift

想优化一个结晶结构（IBRION=2,EDIFFG=-0.02,ISIF=3)，力都收敛了，但total drift很大。
试从ISIF＝３改为ISIF＝２，ENCUT＝850，PREC=Accurate，　NGX，NGXF等看过是都够了，KPOINTS肯定是收敛了。再优化，仍是力收敛了，但total drift很大。

奇怪是total drift应该是力的总和，但我把力加起来，明明是零的，但total drift却不是零!

那位大哥知道是怎回事吗？

POSITION                                     TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
   3.03037    6.79848    -0.95638       -0.001622    0.001648    0.016269
   1.11615    5.23970    2.04034       0.003515    -0.005411    -0.012153
   3.89655    3.68530    -2.99248       -0.006450    0.000185    -0.012153
   4.22533    6.12013    4.07643       0.002251    -0.000527    0.016269
   4.18524    8.55078    -3.66112       0.012518    -0.011454    0.000029
   0.44635    2.02962    1.90239       -0.003328    0.002883    0.007199
   1.49779    1.43188    6.93520       0.004194    -0.001341    0.007199
   5.12068    8.01820    1.37170       -0.016295    0.004729    0.000029
   5.01700    9.30325    -2.68292       0.007169    0.004409    -0.017688
   3.24417    7.65227    -2.94239       -0.015314    0.002170    0.004015
   3.37162    2.95685    -1.07530       0.001326    0.003261    -0.001576
   1.53269    1.63532    2.84314       -0.003387    -0.008844    0.005736
   0.20145    0.93312    0.92251       0.001939    -0.006475    -0.003204
   0.83215    3.26323    1.17449       -0.011749    0.009088    0.003209
   4.82261    5.69824    -0.68234       0.001664    -0.000203    0.004746
   5.03770    7.76765    -4.60115       0.012907    -0.003121    0.004268
   0.68891    0.65215    5.95532       -0.006538    -0.001713    -0.003204
   2.35043    2.40329    6.20730       0.013875    -0.005304    0.003209
   2.35314    7.07801    4.35048       -0.001039    0.001315    0.004746
   4.00833    8.33860    0.43166       -0.009382    0.009398    0.004268
   5.33039    9.12004    2.34990       0.000034    0.008416    -0.017688
   4.84314    6.74701    2.09043       0.009823    -0.011947    0.004015
   0.76728    4.41237    3.95751       0.002095    0.002829    -0.001576
   0.59576    2.15434    7.87595       -0.005789    -0.007494    0.005736
   2.35591    5.41769    0.48397       -0.000887    -0.008388    -0.006146
   2.32523    5.21864    -2.15697       -0.001825    0.005419    -0.001853
   3.24321    4.69585    2.87585       0.005577    0.001261    -0.001853
   3.39723    4.82564    5.51679       -0.006701    -0.005122    -0.006146
   1.76370    6.24487    6.30817       0.010600    -0.000561    -0.006153
   1.48523    3.92494    -0.86091       0.005357    0.004973    0.003301
   2.57562    3.30532    4.17191       0.001459    0.007162    0.003301
   4.41791    4.76134    1.27535       -0.005995    0.008760    -0.006153
-----------------------------------------------------------------------------------
total drift:                            -0.002100    -0.003738    0.118488

回复此楼