24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者lingfaling将赠送您 95 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lingfaling

金虫 (著名写手)

应助: 16 (小学生)
金币: 6643.8
散金: 560
红花: 4
帖子: 1249
在线: 95.4小时
虫号: 3901838
注册: 2015-06-01
专业: 凝聚态物性 II ：电子结构

[求助] 求助大佬CL-NEB计算问题已有1人参与

计算一步收敛，受力确为0。但这也只是针对某些体系，有的体系能够正常计算，并得出结果。检查过中间的结构，没得问题。

主要参数设置：
NSW = 200    # number of steps for ionic upd
IBRION =    3    # ionic relaxation 0=MD 1=quasi-New 2=CG
POTIM  = 0
ICHAIN = 0
IMAGES = 5
SPRING = -5
LCLIMB = T
LTANGENT = T
LTANGENTOLD = F
LDNEB = F
NEBCELL = F

IOPT = 1
MAXMOVE = 0.2
ILBFGSMEM = 20
LGLOBAL = .TRUE.
INVCURV = 0.01
LLINEOPT = .FALSE

结果1:
1 F= -.13251778E+03 E0= -.13248957E+03  d E =-.132518E+03
2 F= -.13251783E+03 E0= -.13248962E+03  d E =-.433186E-04
3 F= -.13251784E+03 E0= -.13248963E+03  d E =-.123478E-04
4 F= -.13251785E+03 E0= -.13248964E+03  d E =-.836182E-05
5 F= -.13251785E+03 E0= -.13248964E+03  d E =-.579447E-05

结果2：
0       0.000000    -128.224300       0.000000
1       0.000000    -132.489600       -4.265300
2       0.000000    -132.340900       -4.116600
3       0.000000    -132.208900       -3.984600
4       0.000000    -132.154000       -3.929700
5       0.000000    -132.165600       -3.941300
6       0.000000    -127.784600       0.439700 ：

回复此楼