计算化学公社

 找回密码 Forget password
 注册 Register
Views: 5376|回复 Reply: 5
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] 求助:电催化CO2还原计算CO*偶联失败

[复制链接 Copy URL]

33

帖子

0

威望

193

eV
积分
226

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 愚不可及 于 2021-8-29 21:20 编辑

各位大佬,小弟现在在做电催化二氧化碳的计算,在CO偶联到OCCO这一步总是偶联失败
参照的文献是Noskov的文献J. Phys. Chem. Lett. 2015, 6, 11, 2032–2037。然后按照他的说法在表面加入了一层六元环状的水分子
左边是我隐藏了除去水之外其他原子的模型,排布和文献是基本一致的。
然后这是我的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.

优化完结构全乱了,如图

请问这种情况应该咋办呢,求指教



1万

帖子

0

威望

9854

eV
积分
22088

Level 6 (一方通行)

2#
发表于 Post on 2021-8-29 22:54:28 | 只看该作者 Only view this author
找到优化过程中最后一个没有乱掉的结构,从那个结构开始重新优化,如果还不解决问题,把最后一个没有乱掉的结构贴上来
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

33

帖子

0

威望

193

eV
积分
226

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2021-8-31 12:09:32 | 只看该作者 Only view this author
本帖最后由 愚不可及 于 2021-8-31 12:15 编辑
wzkchem5 发表于 2021-8-29 22:54
找到优化过程中最后一个没有乱掉的结构,从那个结构开始重新优化,如果还不解决问题,把最后一个没有乱掉的 ...

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

33

帖子

0

威望

193

eV
积分
226

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2021-8-31 14:50:34 | 只看该作者 Only view this author
wzkchem5 发表于 2021-8-29 22:54
找到优化过程中最后一个没有乱掉的结构,从那个结构开始重新优化,如果还不解决问题,把最后一个没有乱掉的 ...

大佬你好,
我按照你说的准备优化,然后发现之前优化完第一步是这样的

然后第二步是这样
然后我接着刚才的结果重跑,然后第二步又就直接飞了,是这样的

看样子和重跑之前是基本一致的。
求大佬指教



1万

帖子

0

威望

9854

eV
积分
22088

Level 6 (一方通行)

5#
发表于 Post on 2021-8-31 17:06:55 | 只看该作者 Only view this author
愚不可及 发表于 2021-8-31 07:50
大佬你好,
我按照你说的准备优化,然后发现之前优化完第一步是这样的

首先确定一下初始结构的键长键角是否合理。从图的这个角度看不出来。比如金属和碳的距离如果太近了,结构优化时可能就会把碳原子崩飞,因为碳原子向外的受力太大。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

33

帖子

0

威望

193

eV
积分
226

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2021-9-1 16:29:52 | 只看该作者 Only view this author
wzkchem5 发表于 2021-8-31 17:06
首先确定一下初始结构的键长键角是否合理。从图的这个角度看不出来。比如金属和碳的距离如果太近了,结构 ...

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




本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-18 16:38 , Processed in 0.186718 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list