计算化学公社

标题: 游离态分子加氢的过渡态计算，dist返回值大于5，能垒过低 [打印本页]

---

作者
Author: Feynman_nur    时间: 2024-10-21 12:10
标题: 游离态分子加氢的过渡态计算，dist返回值大于5，能垒过低

各位老师好，我最近计算游离态的苯乙炔分子在掺杂了Cu的Al(001)表面上加氢的过渡态，感觉上这本该是分子从游离态到吸附态然后再进行加氢两个独立的过程，并且dist.pl返回值为6左右，但合作者及老板要求尽可能当做一个过程去处理。结果就得到了一个很小的能垒(≈0.05eV),然后又尝试将该过渡态当做初态，末态保持不变，再进行一步过渡态计算，结果依旧不理想，甚至插点的相对能量均为负值。所以烦请卡卡老师@卡开发发 、社长老师@sobereva及各位大佬提出宝贵建议和批评，初猜的路径文件附上。

目前已经进行的尝试：

1.采用idpp插值
2.先采用IOPT=7跑20步，然后改为IOPT=1，依照neb.dat的切向力决定要不要进行dimer计算。

---

作者
Author: 卡开发发    时间: 2024-10-21 17:52
我认为比较难判断。虽然这里用vdw-optb86b泛函，但这一类的泛函（vdW-DF系列）其实算的比较慢，400eV的动能截断我并不能确定是否充分收敛，这样可能会导致初始结构（也就是能量稍高那个结构）并非真正处于极小点。另外末态的结构我不确定有没有C≡C更贴近金属表面的结构，包括H进攻方向可能也有一些其他可能，这些可能需要再尝试判断。

---

作者
Author: Feynman_nur    时间: 2024-10-21 19:18

卡开发发 发表于 2024-10-21 17:52
我认为比较难判断。虽然这里用vdw-optb86b泛函，但这一类的泛函（vdW-DF系列）其实算的比较慢，400eV的动能 ...

好的，感谢老师，我再试试

---

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

Powered by Discuz! X3.3