24小时热门版块排行榜    

查看: 3979  |  回复: 9

艾小疯

银虫 (小有名气)

[求助] NEB和CINEB的INCAR设置有什么区别? 已有2人参与

各位前辈,我刚接触用VASP算过渡态,也看过一些资料,但还是有些不明白的地方,希望得到大家的帮助。
1. VASP自带的NEB和编译后的CINEB在INCAR设置上有什么不同?只是LCLIMB= .FALSE. 或.TURE.的区别吗?
2. IOPT应该怎么使用?NEB和CINEB都可以用吗?
3. 初态和终态的文件里除了POSCAR 还要放入OUTCAR吗?有什么作用呢?
4.算了一个过渡态2000步还没有收敛,有什么方法可以让收敛快一点或者更好的找到过渡态吗?
希望能得到前辈的解答,感激不尽~祝大家科研顺利。
回复此楼

» 收录本帖的淘帖专辑推荐

VASP and MS

» 猜你喜欢

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

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

porominix

银虫 (小有名气)


ljw4010: 金币+1, 谢谢交流! 2014-11-15 14:36:52
同是沦落人…
我现在设IOPT=7(才发现VASP手册里居然提都没提IOPT,不过这个帖子里说IOPT=1更合理些:http://muchong.com/bbs/viewthread.php?tid=2139534),EDIFFG=-0.05(这好像还有点粗了,他们都设成0.02甚至0.01)。
初态和终态的文件夹不用放OUTCAR,里面就只有POSCAR,需要自己事先优化好之后再放进来,NEB不对这两个态进行优化。
你的2000步是电子步(NELM)还是原子步(NSW)?是原子步的话也太凶残了。。。我现在设NSW=60,已经跑到第9轮了。。。
2楼2014-11-14 22:21:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

艾小疯

银虫 (小有名气)

引用回帖:
2楼: Originally posted by porominix at 2014-11-14 22:21:51
同是沦落人…
我现在设IOPT=7(才发现VASP手册里居然提都没提IOPT,不过这个帖子里说IOPT=1更合理些:http://muchong.com/bbs/viewthread.php?tid=2139534),EDIFFG=-0.05(这好像还有点粗了,他们都设成0.02甚至0. ...

最近真的被这个过渡态折磨的不行。算了俩月了,一个有用的数据都没出来,调了各种参数都不见效。那2000是NSW,很凶残吧?可是我会告诉你现在跑到3000步了!还是不收敛啊。。。累觉不爱。希望有大神来解救一下啊~~
3楼2014-11-14 23:30:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
艾小疯: 金币+5, ★★★很有帮助, 谢谢。 2014-11-14 23:36:22
一言难尽。

编译了vtstcode代码的vasp才带有LCIMB标签和定义IOPT. 请搜索neb+henkelman找具体说明和用法。

另外,你在学习计算的时候,同课题组没有师兄师姐?
4楼2014-11-14 23:31:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

艾小疯

银虫 (小有名气)

引用回帖:
4楼: Originally posted by cenwanglai at 2014-11-14 23:31:28
一言难尽。

编译了vtstcode代码的vasp才带有LCIMB标签和定义IOPT. 请搜索neb+henkelman找具体说明和用法。

另外,你在学习计算的时候,同课题组没有师兄师姐?

有师兄师姐,但他们都不算过渡态的,说是很不好算,但是老板又要求我做这个,没办法...自己硬啃...区长,能再指点一二吗?感谢~
5楼2014-11-14 23:35:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

porominix

银虫 (小有名气)

引用回帖:
3楼: Originally posted by 艾小疯 at 2014-11-14 23:30:36
最近真的被这个过渡态折磨的不行。算了俩月了,一个有用的数据都没出来,调了各种参数都不见效。那2000是NSW,很凶残吧?可是我会告诉你现在跑到3000步了!还是不收敛啊。。。累觉不爱。希望有大神来解救一下啊~~...

你的IMAGE是多少,还有体系里有多少原子?
我设IMAGE=5,体系80多个原子
6楼2014-11-15 11:11:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

porominix

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
艾小疯(ljw4010代发): 金币+2, 谢谢交流! 2014-11-15 14:37:30
艾小疯: 金币+5, ★★★很有帮助 2014-11-16 11:34:03
引用回帖:
3楼: Originally posted by 艾小疯 at 2014-11-14 23:30:36
最近真的被这个过渡态折磨的不行。算了俩月了,一个有用的数据都没出来,调了各种参数都不见效。那2000是NSW,很凶残吧?可是我会告诉你现在跑到3000步了!还是不收敛啊。。。累觉不爱。希望有大神来解救一下啊~~...

参考:http://theory.cm.utexas.edu/vasp/optimizers.html?highlight=potim

你的IBRION和POTIM分别是3和0吗?
先设成这样吧,然后调IOPT。
如果你的两三千步虽然不收敛、但能量的“震荡”也已经缩小到一个很小的阈值(比如0.5或者0.1ev),那试着把IOPT改成1或者2吧。
我现在的IOPT=7,快速,但是不够精确,现在我的能量“震荡”已经缩得很小了。IOPT=1或者2时,精确的优点就凸显出来了,我打算等这一轮NSW算完了就修改IOPT。
7楼2014-11-15 13:17:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

艾小疯

银虫 (小有名气)

引用回帖:
7楼: Originally posted by porominix at 2014-11-15 13:17:30
参考:http://theory.cm.utexas.edu/vasp/optimizers.html?highlight=potim

你的IBRION和POTIM分别是3和0吗?
先设成这样吧,然后调IOPT。
如果你的两三千步虽然不收敛、但能量的“震荡”也已经缩小到一个很小 ...

多谢多谢,我现在外地开会,等回去一定改改看。

[ 发自手机版 http://muchong.com/3g ]
8楼2014-11-16 11:33:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxh0402

铁虫 (小有名气)

引用回帖:
8楼: Originally posted by 艾小疯 at 2014-11-16 11:33:42
多谢多谢,我现在外地开会,等回去一定改改看。
...

你好,你的问题解决了吗
9楼2015-06-01 21:39:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿克飞顿

金虫 (正式写手)

引用回帖:
3楼: Originally posted by 艾小疯 at 2014-11-14 23:30:36
最近真的被这个过渡态折磨的不行。算了俩月了,一个有用的数据都没出来,调了各种参数都不见效。那2000是NSW,很凶残吧?可是我会告诉你现在跑到3000步了!还是不收敛啊。。。累觉不爱。希望有大神来解救一下啊~~...

80多个原子的体系,插入4个构象,NEB算了3天,才600原子步,还没算完。感受到了这个世界深深地恶意。
10楼2015-07-13 18:58:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 艾小疯 的主题更新
信息提示
请填处理意见