计算化学公社

标题: 关于ORR计算OOH、OH吸附构型问题 [打印本页]

作者
Author:
ps悦悦    时间: 2021-6-27 12:19
标题: 关于ORR计算OOH、OH吸附构型问题
最近在用vasp做表面ORR计算,发现吸附OOH和OH的时候,初始OOH摆放位置和角度不同,优化出来的结构也不同,导致每种吸附都需要计算四种可能的情况,请问有什么办法可以减少他的计算量,使他优化到全局最小值点吗。

作者
Author:
jiangning198511    时间: 2021-6-28 15:07
你可以提高收敛精度试试,有些构型就会变成其他的构型了,表面优化不容易收敛到最优结构,需要多尝试初始结构
作者
Author:
ps悦悦    时间: 2021-6-28 19:46
jiangning198511 发表于 2021-6-28 15:07
你可以提高收敛精度试试,有些构型就会变成其他的构型了,表面优化不容易收敛到最优结构,需要多尝试初始结 ...

ediff=1e-6,ediffg=-0.02,精度已经可以了,我现在就是每种吸附都需要每种两三种初始结构,感觉计算量有点大
作者
Author:
wzkchem5    时间: 2021-6-28 22:18
ps悦悦 发表于 2021-6-28 12:46
ediff=1e-6,ediffg=-0.02,精度已经可以了,我现在就是每种吸附都需要每种两三种初始结构,感觉计算量有点 ...

每个吸附物种,最稳定的构型都是同一种构型吗?如果是的话,有可能可以以你算过的物种都是这个构型最稳定为理由,后续的物种都只算这一个构型。
不过不管怎样做,每少算一个构型,严谨性都会打一点折扣,只能谈论如何让计算结果相对没那么不严谨一些,而没有任何办法严格证明某个构型肯定不需要算。
作者
Author:
ps悦悦    时间: 2021-6-28 22:54
wzkchem5 发表于 2021-6-28 22:18
每个吸附物种,最稳定的构型都是同一种构型吗?如果是的话,有可能可以以你算过的物种都是这个构型最稳定 ...

这几天一直在尝试各种OOH和OH吸附结构,大概得出了结论,O-H键通常是朝向得电子较多的原子方向,先计算bader电荷,基本可以推出O-H的方向。
作者
Author:
jiangning198511    时间: 2021-6-29 09:20
ps悦悦 发表于 2021-6-28 19:46
ediff=1e-6,ediffg=-0.02,精度已经可以了,我现在就是每种吸附都需要每种两三种初始结构,感觉计算量有点 ...

我不清楚你的表面是什么情况,以Au111表面而言,就有T H B 等位置,如果是复杂表面 设计的吸附位置就更多了




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3