计算化学公社

标题: 求助ORCA计算过渡金属吸附小分子机组推荐 [打印本页]

作者
Author: XxxXOO 时间: 2025-3-16 11:11
标题: 求助ORCA计算过渡金属吸附小分子机组推荐
体系一共有60个原子，58个银，吸附气体小分子，结构优化用PBE0 def2-SV(P),单点能计算用PBE0 def2-TZVP
算出的吸附能与文献有一定差距，有什么好的推荐吗？

作者
Author: cokie 时间: 2025-3-16 11:27
本帖最后由 cokie 于 2025-3-16 11:29 编辑

基组

给的信息有些少，需要补充大家才好知道问题在哪里。计算结果的不合理/“错误”等等问题不一定只是因为泛函/基组选择不合适，还有其他可能没注意到的问题：比如模拟的是固体表面吸附，还是银团簇在溶液中的吸附；读取的自由能还是电子能；银团簇的结构是否构建合理；是否考虑了BSSE等等

都用orca了，60个原子也可以考虑 wB97M-V/def2-TZVP/gCP 来算吸附能，也不贵。

最后吐槽一下，既然要还原文献计算结果了，为何不用和文献里一样的计算方法呢

作者
Author: XxxXOO 时间: 2025-3-16 12:51

cokie 发表于 2025-3-16 11:27
基组

给的信息有些少，需要补充大家才好知道问题在哪里。计算结果的不合理/“错误”等等问题不一定只是 ...

我是研究固体表面吸附，读取电子能，文献都是用的vsap，我想我ORCA算一下，对比吸附能差0.1eV

作者
Author: Stardust0831 时间: 2025-3-16 16:50
既然是研究固体表面吸附，那更好的策略是用第一性程序（比如ABACUS、CP2K、VASP等）来当成周期性体系计算。非要在ORCA中计算的话，应当确认簇模型建模合理，见：《使用量子化学程序基于簇模型计算金属表面吸附问题》（http://sobereva.com/540），结构优化可以用PBE0-D4/def2-SVP/gCP，单点可以按楼上说的wB97M-V/def2-TZVP/gCP。

作者
Author: XxxXOO 时间: 2025-3-16 20:06

Stardust0831 发表于 2025-3-16 16:50
既然是研究固体表面吸附，那更好的策略是用第一性程序（比如ABACUS、CP2K、VASP等）来当成周期性体系计算。 ...

我看sob老师之前的文章不是说在计算过程中尽量不要更换泛函吗？

作者
Author: cokie 时间: 2025-3-16 20:47

XxxXOO 发表于 2025-3-16 20:06
我看sob老师之前的文章不是说在计算过程中尽量不要更换泛函吗？

http://bbs.keinsci.com/thread-9878-1-1.html（量子化学研究中切换泛函应当注意的问题）

不要切换泛函，是指没有什么特别的原因时不要切换，如果在计算单点的时候选择明显精度更高且更加适合的泛函当然没有问题。简单来说，由于几何优化以及频率计算相比算单点，耗时耗力程度几何倍数的增加，但几何优化和频率计算却又不是特别吃泛函和基组的精度级别（只要合理选择即可），因此几何优化和频率计算可以选择精度并不那么高的计算方法同时最小化对结果精度的影响。

都用ORCA了，没有任何理由在算力足够的前提下，不选择wB97M-V/def2-TZVP/gCP，6.0.1算着飞速

作者
Author: Stardust0831 时间: 2025-3-16 21:42

XxxXOO 发表于 2025-3-16 20:06
我看sob老师之前的文章不是说在计算过程中尽量不要更换泛函吗？

当前有足够理由切换泛函。wB97M-V只有解析梯度，没有解析hessian，所以如果直接用wB97M-V来算振动分析的话会相当耗时；而wB97M-V又是目前杂化泛函里算能量名列前茅的，用wB97M-V来算单点能在大多数情况下可以得到比PBE0更精确的结果。就和社长在博文量子化学研究中切换泛函应当注意的问题中说的那样：

A泛函的耗时明显比B泛函低，而精度也明显比B泛函差。因此用A泛函计算那些对计算级别要求低而且计算量较大的任务，而用B泛函较准确地计算那些对计算级别要求高的问题，即好钢用在刀刃上。

此外，关于wB97M-V的一些介绍可以看：《求简单科普wb97家族》

作者
Author: Stardust0831 时间: 2025-3-16 21:48

cokie 发表于 2025-3-16 20:47
http://bbs.keinsci.com/thread-9878-1-1.html（量子化学研究中切换泛函应当注意的问题）

不要切换泛 ...

其实是有理由不上gCP的。原因很简单，gCP的参数只支持前四周期，第五周期及以后的元素会用同族的第四周期的元素的参数来近似替代。Ag在第五周期，所以上gCP时，Ag其实在用Cu的参数，在研究Ag表面的吸附时，这种近似替代的可靠性就值得怀疑。

作者
Author: XxxXOO 时间: 2025-3-18 09:38

Stardust0831 发表于 2025-3-16 21:48
其实是有理由不上gCP的。原因很简单，gCP的参数只支持前四周期，第五周期及以后的元素会用同族的第四周期 ...

谢谢谢谢

作者
Author: XxxXOO 时间: 2025-3-18 09:39

cokie 发表于 2025-3-16 20:47
http://bbs.keinsci.com/thread-9878-1-1.html（量子化学研究中切换泛函应当注意的问题）

不要切换泛 ...

好的，谢谢谢谢

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