求助：Ni111表面化学吸附cooh结构优化

lhx-xwzj · 发表于 Post on 2024-1-25 17:29:09

本帖最后由 lhx-xwzj 于 2024-1-25 17:41 编辑

各位老师好，学生最近用QE7.0计算Ni111表面化学吸附cooh，输入文件是用QEtooltik生成的，无论初始结构怎么摆（见图优化前1和2），最后都会优化成CO和OH两个片段（见优化后），但文献中是能存在Ni111和cooh稳定吸附的情况，下面我的是输入文件

Ni111-cooh.relax.in (4.89 KB, 下载次数 Times of downloads: 12) ，输出文件过大无法上传，我贴了优化前后的图如下，不知道什么原因，请大家赐教，谢谢。

fanzhichao001 · 发表于 Post on 2024-1-25 18:02:49

会不会是初始结构离表面过近的原因呢？

123wo · 发表于 Post on 2024-1-25 18:27:05

楼上说的有道理，另外可以检查一下吸附能有多少，有时候它确实解离了，可是吸附能是很正的。另外COOH应该按照顺式和反式两种样子摆放。试试吧

pal · 发表于 Post on 2024-1-25 19:01:29

COOH是个什么状态呢，有没有单个电子之类的呢

Exodus · 发表于 Post on 2024-1-25 20:50:35

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

lhx-xwzj · 发表于 Post on 2024-1-25 21:34:50

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

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

lhx-xwzj · 发表于 Post on 2024-1-25 21:37:34

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

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

lhx-xwzj · 发表于 Post on 2024-1-25 21:39:41

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

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

lhx-xwzj · 发表于 Post on 2024-1-25 21:41:27

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

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

123wo · 发表于 Post on 2024-1-26 00:08:54

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

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

lhx-xwzj · 发表于 Post on 2024-1-29 17:50:06

本帖最后由 lhx-xwzj 于 2024-2-28 12:39 编辑

顶一下

lhx-xwzj · 发表于 Post on 2024-2-28 12:40:39

目前还未解决，还有解决办法嘛，求指点

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Quantum ESPRESSO] 求助：Ni111表面化学吸附cooh结构优化