24小时热门版块排行榜    

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

julyrush5530

铁虫 (初入文坛)

[求助] NEB计算问题求助 已有1人参与

模型:溴苯在Cu111上断裂C-Br键, NEB过渡态搜索。
问题:不能收敛,计算差不多2个小时,导出发现各个Image苯环结构明显变形,(因为初末态原子不匹配,我把它们的溴苯原子先匹配好,然后都采用末态的基底,初态结构几乎没有变化,这样产生的Image苯环都是完整的,反应路径也是和我预想大致一样)。
我尝试过,固定断键的原子,先让每个Image独自优化,然后再替换POSCAR,删除每个Image不必要的文件,再做NEB。这个的话跑了很久,苯环结构没有坏,但和以前一样,01Image能量总是不能和好的收敛,波动很大,力更是变动得离谱。

NEB计算问题求助
INCAR01.png


NEB计算问题求助-1
INCAR02.png


NEB计算问题求助-2
nebef.png


NEB计算问题求助-3
nebef01.png
回复此楼

» 猜你喜欢

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

看清楚看明白
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

julyrush5530

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by 490105179 at 2015-05-27 08:58:36
“因为初末态原子不匹配,我把它们的溴苯原子先匹配好,然后都采用末态的基底,初态结构几乎没有变化,这样产生的Image苯环都是完整的,反应路径也是和我预想大致一样”1.初末态原子都不匹配如何做NEB?2.你处理后的 ...

谢谢大侠的回复。
1、2:我因为要把基地的最下面两层固定,使用wxDragon产生的POSCAR(它把原子按Z轴高度排列,打乱了原来的排序)方便我固定。请问您是用什么产生POSCAR的,VESTA可以,但好像得自己一个个挑出来固定,我的体系有140多个原子,这样子太不方便了。
我觉得您说得是对的,我也考虑这个稳态问题,所以替换基底后把初态跑了一次发现整个过程能量几乎没变。吸附态能量在表面上几乎一样的,所以我也就默认了,主要是我没有很好的方法产生POSCAR。
3,4: 这里是因为我跑了一个把小时苯环结构就散掉了,当初我认为是差值产生的IMAGE结构不合理(其实没多大问题,键长什么的都还在合理范围内),我也没办法,想着先固定断键的两个原子然他们单独优化(不是NEB),优化完了我发现能量基本是一个合理的抛物线,然后我再把,每个IMAGE的CONTCAR拷成POSCAR(并且把POSCAR的元素符去掉,松开固定原子),删掉INCAR,KPOINTS,脚本文件。让NEB接着这个跑,这样的话跑了接近一天,其中有一个态的能量和文献的过渡态能量很接近了,一直在那个值附近波动。但第一个IMAGE总是不能很好收敛其受的力大的离谱,能量时高时低。

希望你给我指点下,我觉得产生POSCAR很关键,因为以前做H2分解时,NEB就可以跑,只是没加VDW能量不准,但其他都很接近。
看清楚看明白
3楼2015-05-27 10:10:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

490105179

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
julyrush5530: 金币+2, ★★★很有帮助 2015-05-27 10:13:29
“因为初末态原子不匹配,我把它们的溴苯原子先匹配好,然后都采用末态的基底,初态结构几乎没有变化,这样产生的Image苯环都是完整的,反应路径也是和我预想大致一样”1.初末态原子都不匹配如何做NEB?2.你处理后的初态没有再次优化,说明结构是不稳定的,也不能进行下一步计算。
“我尝试过,固定断键的原子,先让每个Image独自优化,然后再替换POSCAR,删除每个Image不必要的文件,再做NEB。”3.你把NEB中运动的原子都固定了再计算还有什么意义?4.VASP是等间距插点,你删除其中的某些对计算肯定有影响。
我不做NEB对你的INCAR不评论。关于nebef中第二列是Image间距,你那明显不合理。
所以你还是老老实实的一步步来,先优化好初态,然后由初态得到末态,再根据初末态做NEB吧
2楼2015-05-27 08:58:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

julyrush5530

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by 490105179 at 2015-05-27 08:58:36
“因为初末态原子不匹配,我把它们的溴苯原子先匹配好,然后都采用末态的基底,初态结构几乎没有变化,这样产生的Image苯环都是完整的,反应路径也是和我预想大致一样”1.初末态原子都不匹配如何做NEB?2.你处理后的 ...

您说的nebef中第二列IMAGE间距不合理是指那个数字大小分布不均?
看清楚看明白
4楼2015-05-27 10:21:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

490105179

木虫 (正式写手)

【答案】应助回帖


liliangfang: 金币+1, 谢谢交流 2015-05-29 07:41:27
VASP可以自己插点,你装上nebmake.pl插件,分别把初末态设为POSCAR1 POSCAR2,执行
nebmake.pl POSCAR1 POSCAR2 X  X为插点数,可以自己产生Images
5楼2015-05-27 14:45:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 274求调剂 +9 S.H1 2026-03-18 9/450 2026-03-20 20:54 by zhukairuo
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 271材料工程求调剂 +7 .6lL 2026-03-18 7/350 2026-03-20 09:10 by xingguangj
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
信息提示
请填处理意见