计算化学公社

标题: 求助：电催化CO2还原计算CO*偶联失败 [打印本页]

作者
Author: 愚不可及 时间: 2021-8-29 21:14
标题: 求助：电催化CO2还原计算CO*偶联失败
本帖最后由愚不可及于 2021-8-29 21:20 编辑

各位大佬，小弟现在在做电催化二氧化碳的计算，在CO偶联到OCCO这一步总是偶联失败
参照的文献是Noskov的文献J. Phys. Chem. Lett. 2015, 6, 11, 2032–2037。然后按照他的说法在表面加入了一层六元环状的水分子 (, 下载次数 Times of downloads: 36) (, 下载次数 Times of downloads: 30)
左边是我隐藏了除去水之外其他原子的模型，排布和文献是基本一致的。
然后这是我的INCAR信息
ISTART = 0; ICHARG = 2

general:
SYSTEM = demo surface
ENCUT = 500
ISMEAR = 2; SIGMA = 0.2
EDIFF = 1E-3
IVDW=11

dynamic:
NSW = 1500
POTIM = 0.5
IBRION = 2
EDIFFG = -0.2
ISIF = 3

special:
LREAL = A
NSIM = 4
NCORE = 4
LCHARG =.FALSE.
LWAVE = .FALSE.
ISOL = .FALSE.

优化完结构全乱了，如图
(, 下载次数 Times of downloads: 33)
请问这种情况应该咋办呢，求指教

作者
Author: wzkchem5 时间: 2021-8-29 22:54
找到优化过程中最后一个没有乱掉的结构，从那个结构开始重新优化，如果还不解决问题，把最后一个没有乱掉的结构贴上来

作者
Author: 愚不可及 时间: 2021-8-31 12:09
本帖最后由愚不可及于 2021-8-31 12:15 编辑

wzkchem5 发表于 2021-8-29 22:54
找到优化过程中最后一个没有乱掉的结构，从那个结构开始重新优化，如果还不解决问题，把最后一个没有乱掉的 ...

大佬，感谢帮助。
我再尝试一下，多谢多谢

作者
Author: 愚不可及 时间: 2021-8-31 14:50

wzkchem5 发表于 2021-8-29 22:54
找到优化过程中最后一个没有乱掉的结构，从那个结构开始重新优化，如果还不解决问题，把最后一个没有乱掉的 ...

大佬你好，
我按照你说的准备优化，然后发现之前优化完第一步是这样的
(, 下载次数 Times of downloads: 33)
然后第二步是这样 (, 下载次数 Times of downloads: 38)
然后我接着刚才的结果重跑，然后第二步又就直接飞了，是这样的
(, 下载次数 Times of downloads: 37)
看样子和重跑之前是基本一致的。
求大佬指教

作者
Author: wzkchem5 时间: 2021-8-31 17:06

愚不可及发表于 2021-8-31 07:50
大佬你好，
我按照你说的准备优化，然后发现之前优化完第一步是这样的

首先确定一下初始结构的键长键角是否合理。从图的这个角度看不出来。比如金属和碳的距离如果太近了，结构优化时可能就会把碳原子崩飞，因为碳原子向外的受力太大。

作者
Author: 愚不可及 时间: 2021-9-1 16:29

wzkchem5 发表于 2021-8-31 17:06
首先确定一下初始结构的键长键角是否合理。从图的这个角度看不出来。比如金属和碳的距离如果太近了，结构 ...

您好，按照您说的尝试了一下，将CO原子和水层都往远离Cu的方向拖了一下，情况稍微有点好转，不是马上就会结构乱掉。但是在几步之后仍然会乱掉。
我的模型信息已经上传了。能否拜托您帮忙看一下
十分感谢！ (, 下载次数 Times of downloads: 7)

欢迎光临计算化学公社 (http://bbs.keinsci.com/)