24小时热门版块排行榜

返回列表

xiaoluoj

金虫 (正式写手)

应助: 12 (小学生)
金币: 690.1
散金: 5
红花: 17
帖子: 487
在线: 504.1小时
虫号: 1526088
注册: 2011-12-06
专业: 通信理论与系统

[求助] 有关 NEB 的问题

最近想用NEB算一些问题，由于不确定是否编译了ci-nbe，所以采用vasp自带的NEB方法，INCAR 设置如下
---------------------------------------------------------------------------------------------------
ISTART=0
LWAVE=.False.
PREC=high
#LVTOT=.True.
#IDIPOL=4
ALGO=Fast
LREAL=Auto
#LELF=.TRUE.
LCHARG=.FALSE.
################electronic relaxation############
NELM=80
ENCUT=500
EDIFF=1E-05
#ISPIN=2
###################ion relaxtion################
NSW=400
IBRION=1
#POTIM=0
ISIF=2
EDIFFG=-0.015
####################################################
ISMEAR=0
SIGMA=0.05
#LORBIT=11
#NEDOS=
#NBANDS =
ISYM=0
#NELECT
###################NEB########################################
ICHAIN=0
IMAGES=5
SPRING=-5
#LCLIMB=.TRUE.
#IOPT=3
------------------------------------------------------------------------------------------------------
  运算过程中，用脚本nebef查看受力情况
[dyh2000@ln01 4_NEB_IBRION=1_EDIFFG=-0.015]$ ./nebef.pl
0       0.000000    -2247.462489       0.000000
1       0.000000    -2247.442673       0.019816
2       0.000000    -2247.385011       0.077478
3       0.000000    -2247.315143       0.147346
4       0.000000    -2247.278576       0.183913
5       0.000000    -2247.286221       0.176268
6       0.000000    -2247.302143       0.160346
受力情况怎么为0，是不是这个脚本只能处理ci-nbe而不能处理NEB，请高手指教

回复此楼