应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 用VASP计算NEB都十几天了还不收敛，急啊！！！
91/1返回列表
查看: 4803  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lijuanshu

铜虫 (正式写手)

[交流] 用VASP计算NEB都十几天了还不收敛，急啊！！！

我用VASP计算H2分子在合金表面的催化裂解，由于想知道详细的催化裂解过程，取了14个点，NEB从这个月的13号就开始计算了，到目前为止还没收敛，感觉力振动的很厉害，有时候看到都快收敛了，但是过会又变得很大，下面是我的INCAR，还请大家看看设置的是否合理，不胜感激

SYSTEM = *****
ENCUT = 280
NBANDS = 280
ISTART = 0
ICHARG = 2
NPAR = 2
NELMIN = 5
ISMEAR = 0
SIGMA = 0.2
NSW = 1000
IBRION = 3
POTIM = 0
IOPT = 1
ISPIN = 2
NELM = 2000
EDIFF = 1E-4
EDIFFG = -0.04
LREAL = Auto
GGA = 91
VOSKOWN = 1
ALGO = Fast
LPLANE = .TRUE.
NSIM = 4
LCLIMB = .TRUE.
ICHAIN = 0
LGOLOBAL = .TRUE.
MAXMOVE = 0.2
INVCURV = 0.01
SPRING = -5
IMAGES = 14
LWAVE = .FALSE.
LCHARG = .FALSE.
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2012-12-24 10:00:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)
★ ★ ★
lijuanshu(金币+1): 谢谢参与
franch: 金币+2, 鼓励交流。。 2012-12-24 21:35:09
你首先查看一下输出 看看已经跑了多少步了， 14个image有点多 先取少点算
赞一下 回复此楼
2楼2012-12-24 12:00:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阚敏

新虫 (小有名气)
★ ★ ★
lijuanshu(金币+1): 谢谢参与
franch: 金币+2, 鼓励交流。。 2012-12-24 21:35:20
表示我测试VASP的NEB里面，除了跑单个原子的扩散比较快，别的都很慢，而且经常得不到想要的结果呀。另外，你的IMAGES确实取得太多了，你可以先取少一点，完了再在此基础上插入新的IMAGE。
赞一下 回复此楼
3楼2012-12-24 14:10:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
2楼: Originally posted by emilyoyang at 2012-12-24 12:00:23
你首先查看一下输出 看看已经跑了多少步了， 14个image有点多 先取少点算

都算了121个离子步了，由于要知道详细过程，所以点取的很多
赞一下 回复此楼
4楼2012-12-24 22:32:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)
★ ★
fzx2008: 金币+2, 谢谢交流！ 2012-12-25 12:15:26
引用回帖:
4楼: Originally posted by lijuanshu at 2012-12-24 22:32:10
都算了121个离子步了，由于要知道详细过程，所以点取的很多...

看一下 每个离子步都多少步收敛  看看是有收敛的趋势
或者  查看一下 能量和力的收敛趋势  如果变化太剧烈 那估计就是路径跑错了
如果很平稳的话，那就快要收敛了
赞一下 回复此楼
5楼2012-12-24 22:55:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
5楼: Originally posted by emilyoyang at 2012-12-24 22:55:34
看一下 每个离子步都多少步收敛  看看是有收敛的趋势
或者  查看一下 能量和力的收敛趋势  如果变化太剧烈 那估计就是路径跑错了
如果很平稳的话，那就快要收敛了...

中间有一个IMAGE点基本上每个离子步有20步才收敛，其他的差不多是10步收敛，力的变化很大，来回振荡，我是用线性插入产生的点，但是它自己跑的路径是折线，位置也是来回移动的，比如03文件的位置在1, 04文件的位置在2, 05文件的位置应该是在3的，但是却跑回到位置1附近，但是如果路径不对，那究竟该怎么选取路径呢？
赞一下 回复此楼
6楼2012-12-25 09:20:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)
★
fzx2008: 金币+1, 谢谢指教！ 2012-12-25 12:16:06
引用回帖:
6楼: Originally posted by lijuanshu at 2012-12-25 09:20:26
中间有一个IMAGE点基本上每个离子步有20步才收敛，其他的差不多是10步收敛，力的变化很大，来回振荡，我是用线性插入产生的点，但是它自己跑的路径是折线，位置也是来回移动的，比如03文件的位置在1, 04文件的位置 ...

有时候线性插值 生成初始images不太合理  你需要人为调整一下 然后再跑 可能收敛会很快收敛到MEP上
赞一下 回复此楼
7楼2012-12-25 12:13:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
7楼: Originally posted by emilyoyang at 2012-12-25 12:13:52
有时候线性插值 生成初始images不太合理  你需要人为调整一下 然后再跑 可能收敛会很快收敛到MEP上...

但是该怎么取路径我也不是很清楚，所以每次都是线性插入的
赞一下 回复此楼
8楼2012-12-25 15:36:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

04nylxb

木虫 (正式写手)
★ ★
lijuanshu(金币+1): 谢谢参与
sunyang1988: 金币+1, 谢谢交流 2012-12-27 19:53:18
先不跑CNEB，先跑NEB，感觉CNEB不好收敛，尤其是在初始结构不太好的时候
赞一下 回复此楼
9楼2012-12-26 11:46:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lijuanshu 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭