关于ORR计算OOH、OH吸附构型问题

ps悦悦 · 发表于 Post on 2021-6-27 12:19:18

最近在用vasp做表面ORR计算，发现吸附OOH和OH的时候，初始OOH摆放位置和角度不同，优化出来的结构也不同，导致每种吸附都需要计算四种可能的情况，请问有什么办法可以减少他的计算量，使他优化到全局最小值点吗。

jiangning198511 · 发表于 Post on 2021-6-28 15:07:19

你可以提高收敛精度试试，有些构型就会变成其他的构型了，表面优化不容易收敛到最优结构，需要多尝试初始结构

ps悦悦 · 发表于 Post on 2021-6-28 19:46:10

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

ediff=1e-6,ediffg=-0.02,精度已经可以了，我现在就是每种吸附都需要每种两三种初始结构，感觉计算量有点大

wzkchem5 · 发表于 Post on 2021-6-28 22:18:32

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

每个吸附物种，最稳定的构型都是同一种构型吗？如果是的话，有可能可以以你算过的物种都是这个构型最稳定为理由，后续的物种都只算这一个构型。
不过不管怎样做，每少算一个构型，严谨性都会打一点折扣，只能谈论如何让计算结果相对没那么不严谨一些，而没有任何办法严格证明某个构型肯定不需要算。

ps悦悦 · 发表于 Post on 2021-6-28 22:54:42

wzkchem5 发表于 2021-6-28 22:18
每个吸附物种，最稳定的构型都是同一种构型吗？如果是的话，有可能可以以你算过的物种都是这个构型最稳定 ...

这几天一直在尝试各种OOH和OH吸附结构，大概得出了结论，O-H键通常是朝向得电子较多的原子方向，先计算bader电荷，基本可以推出O-H的方向。

jiangning198511 · 发表于 Post on 2021-6-29 09:20:41

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

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

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

[VASP] 关于ORR计算OOH、OH吸附构型问题