24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2865  |  回复: 17
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

身在长沙

新虫 (初入文坛)

[求助] NEB的计算出现负值,测试的其他能量都比初始的00和终止的10低

问题如下:
1.计算NEB的方法:首先利用nebmake.pl产生初始终态结构,用xcrysden查看了路径没问题,然后对初始终态结构进行几何结构优化,接着采用优化后的能量构型进行NEB的计算,按照VASP manual上说明,在INCAR中设置如下:
ISTART = 0; ICHARG = 2
   ENCUT  = 400
   ISMEAR = 2
   LREAL=Auto
NELM=200
NSW=1000
EDIFFG=-0.02
LWAVE = F
LCHARG = F
ICHAIN=0
IMAGES=9
SPRING=-5
LCLIMB=.TRUE.
IBRION=3
POTIM=0
IOPT=7
算完之后用nebbarrier.pl产生的neb.dat如下:
0     0.000000     0.000000     0.000000   0
  1     0.151183    -2.145471    -0.242258   1
  2     0.302705    -2.106921    -0.362263   2
  3     0.452469    -2.051050    -0.403422   3
  4     0.603020    -2.012050    -0.282203   4
  5     0.753502    -1.992645     0.000084   5
  6     0.904120    -2.012146     0.253874   6
  7     1.054939    -2.051164     0.407455   7
  8     1.204729    -2.107090     0.352265   8
  9     1.356227    -2.145477     0.236853   9
10     1.507559     0.000000     0.000000  10
初末位置的能量我多次验算都是最低的,不知道是什么原因,求大神们不吝指教我这个新手!谢谢啦
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

身在长沙

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by emilyoyang at 2013-11-04 16:07:14
参数没什么大的问题
NELM=200 默认的60就行了  不用这么大
IOPT=7 最好用 1  我一般用1 3
IMAGES=9 这么多干嘛 反应路径很长吗  9images的话计算量会很大
方便的话可以把你的initial和final POSCAR贴出来  我看 ...

谢谢你的回复,我的初末位置弛豫用的INCAR都是如下的
SYSTEM =W-bcc
ENCUT = 400
ISTART = 0 ; ICHARG = 2
ISMEAR = 2
NSW = 100 ; IBRION =2
ISIF = 2
EDIFF = 1E-04 ; EDIFFG = -1E-02
LREAL= Auto
PREC = Accurate
虚心求指导!
7楼2013-11-04 16:37:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

liu-biao

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
求楼下解决  我也遇到了 大神们 求支援
2楼2013-11-04 12:16:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
identation: 金币+1, 谢谢交流 2013-11-04 13:53:17
你把03这个结构拿出来 做一下结构优化,看看是否能够回到00结构
不能的话,就是你初始结构没找对,也就是说没有找到能量最低结构
3楼2013-11-04 12:43:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

身在长沙

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by emilyoyang at 2013-11-04 12:43:21
你把03这个结构拿出来 做一下结构优化,看看是否能够回到00结构
不能的话,就是你初始结构没找对,也就是说没有找到能量最低结构

谢谢你的回复,这个我做过了,弛豫之后又回到00结构了。而且我做到这个构型的能量最低位置与文献上的是一致的。我怀疑是参数设置那里出了问题,麻烦你再帮我看看。新手上路,请多多指教啊!
4楼2013-11-04 13:10:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 295分求调剂 +13 ?要上岸? 2026-04-10 13/650 2026-04-12 15:37 by laoshidan
[考研] 求调剂 +16 张番茄不炒蛋 2026-04-10 17/850 2026-04-12 13:58 by 熬夜成!
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 267求调剂 +8 再忙也要吃饭啊 2026-04-09 8/400 2026-04-11 21:42 by cfdbai
[考研] 调剂 +10 只叙离别辞 2026-04-09 12/600 2026-04-11 20:57 by 逆水乘风
[考研] 303求调剂 +14 SereinQ 2026-04-10 15/750 2026-04-11 20:43 by 蓝云思雨
[考研] 生物学调剂 可调剂到生物与医药 +8 李政莹 2026-04-06 9/450 2026-04-11 10:36 by wwj2530616
[考研] 生物学调剂 +8 小冉要努力 2026-04-10 9/450 2026-04-11 10:22 by wwj2530616
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 080500求调剂 +17 黄宇博 2026-04-06 17/850 2026-04-11 08:36 by zhq0425
[考研] 调剂 +12 卷卷卷心菜_ 2026-04-09 13/650 2026-04-10 22:36 by Ftglcn90
[考研] 本科西工大 324求调剂 +4 wysyjs25 2026-04-10 4/200 2026-04-10 20:00 by 来看流星雨10
[考研] 314求调剂 +23 wakeluofu 2026-04-09 24/1200 2026-04-10 15:31 by MOF_Catal
[考研] 285求调剂 +9 AZMK 2026-04-07 11/550 2026-04-10 15:24 by AZMK
[考研] 青岛科技大学材料学院,环境学院调剂补录4月10日以前都可以 +3 1青科大。 2026-04-09 5/250 2026-04-10 09:58 by 翩翩一书生
[考研] 材料专硕初试分332一志愿西北工业大学, +12 故人?? 2026-04-09 12/600 2026-04-09 18:34 by Ccclqqq
[考研] 331求调剂 +5 luoxin0706. 2026-04-08 5/250 2026-04-08 22:15 by zhouyuwinner
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 270求调剂 +3 031127 2026-04-06 4/200 2026-04-08 21:00 by 逆水乘风
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
信息提示
请填处理意见