计算化学公社

标题: 求助：Ni111表面化学吸附cooh结构优化 [打印本页]

作者
Author: lhx-xwzj 时间: 2024-1-25 17:29
标题: 求助：Ni111表面化学吸附cooh结构优化
本帖最后由 lhx-xwzj 于 2024-1-25 17:41 编辑

各位老师好，学生最近用QE7.0计算Ni111表面化学吸附cooh，输入文件是用QEtooltik生成的，无论初始结构怎么摆（见图优化前1和2），最后都会优化成CO和OH两个片段（见优化后），但文献中是能存在Ni111和cooh稳定吸附的情况，下面我的是输入文件 (, 下载次数 Times of downloads: 12) ，输出文件过大无法上传，我贴了优化前后的图如下，不知道什么原因，请大家赐教，谢谢。

(, 下载次数 Times of downloads: 16) (, 下载次数 Times of downloads: 19) (, 下载次数 Times of downloads: 20)

作者
Author: fanzhichao001 时间: 2024-1-25 18:02
会不会是初始结构离表面过近的原因呢？

作者
Author: 123wo 时间: 2024-1-25 18:27
楼上说的有道理，另外可以检查一下吸附能有多少，有时候它确实解离了，可是吸附能是很正的。另外COOH应该按照顺式和反式两种样子摆放。试试吧

作者
Author: pal 时间: 2024-1-25 19:01
COOH是个什么状态呢，有没有单个电子之类的呢

作者
Author: Exodus 时间: 2024-1-25 20:50
波函数和电荷密度截断能太低了，至少50和500，这个体系的k点至少取到441，电子步收敛标准太低，一般设置1e-8，degauss取的太大，取0.005足够了

作者
Author: lhx-xwzj 时间: 2024-1-25 21:34

fanzhichao001 发表于 2024-1-25 18:02
会不会是初始结构离表面过近的原因呢？

谢谢，我之前摆的远一点就会优化到物理吸附构型，物理吸附是没有问题的，但我需要化学吸附构型，我再探索一下。

作者
Author: lhx-xwzj 时间: 2024-1-25 21:37

pal 发表于 2024-1-25 19:01
COOH是个什么状态呢，有没有单个电子之类的呢

没有电荷，就单纯多了个H，自旋多重度为2

作者
Author: lhx-xwzj 时间: 2024-1-25 21:39

Exodus 发表于 2024-1-25 20:50
波函数和电荷密度截断能太低了，至少50和500，这个体系的k点至少取到441，电子步收敛标准太低，一般设置1e- ...

谢谢老师，因为算力有限，我按以上设置后算的太慢了，准备先把K点提高试一下，您觉得哪几个参数比较优先呢？

作者
Author: lhx-xwzj 时间: 2024-1-25 21:41

123wo 发表于 2024-1-25 18:27
楼上说的有道理，另外可以检查一下吸附能有多少，有时候它确实解离了，可是吸附能是很正的。另外COOH应该按 ...

感觉和摆的初始结构关系不大，试过好几个初始结构了，谢谢您，我再试试

作者
Author: 123wo 时间: 2024-1-26 00:08

lhx-xwzj 发表于 2024-1-25 21:41
感觉和摆的初始结构关系不大，试过好几个初始结构了，谢谢您，我再试试

我以前可能做过Cu（111）上的，发现COOH吸附是可以的，但是不知道为什么Ni会出现这样的问题

作者
Author: lhx-xwzj 时间: 2024-1-29 17:50
本帖最后由 lhx-xwzj 于 2024-2-28 12:39 编辑

顶一下

作者
Author: lhx-xwzj 时间: 2024-2-28 12:40
目前还未解决，还有解决办法嘛，求指点

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