计算化学公社

标题: 如何计算显示溶剂模型下的吸附能 [打印本页]

作者
Author: calamity 时间: 2021-1-24 22:21
标题: 如何计算显示溶剂模型下的吸附能
如题，我在Ag（111）面上计算CO的吸附能时，加入了4个H2O分子，首先slab加了4个H2O然后结构优化，随后再CO吸附体系中加了4个H2O然后优化，计算得到的吸附能达到了0.23eV，感觉算的结果很差，请问这样的方法对吗。
是否可以直接优化一个CO+4个H2O分子的结构，然后把CO去掉后算一个单点，由此得到的能量来计算其吸附能？

作者
Author: Penson 时间: 2021-1-25 09:27
我认为原来的计算方法才合理，有无CO都应该做结构优化才能判断两个对比的体系是否都是在给定精度下结构稳定的状态

作者
Author: 卡开发发 时间: 2021-1-25 10:16
不太确定CO初始状态是气态还是溶液。我觉得可能处理的手段是，应该开着VASPSol，这样即便是clean的表面也把溶剂效应考虑了。对于有CO的体系，显式H2O加上连续介质都当做溶剂处理，类似于10.1021/ct050164b，算出溶剂化能加到CO的自由能上就行了。

作者
Author: aaronzjw 时间: 2021-1-25 14:28
4个水可能不一定够

作者
Author: calamity 时间: 2021-1-25 14:40

aaronzjw 发表于 2021-1-25 14:28
4个水可能不一定够

当前计算能力只能算4个水了，请问方法对吗？第一个还是第二个方法

作者
Author: calamity 时间: 2021-1-25 14:41

Penson 发表于 2021-1-25 09:27
我认为原来的计算方法才合理，有无CO都应该做结构优化才能判断两个对比的体系是否都是在给定精度下结构稳定 ...

但是原来的方法对slab上水的构型要求太大了，然后跑AIMD感觉不太来得及。后来我想了一下，可以用第二种方法然后去掉CO以后不算单点，做结构优化来算吸附能吗

作者
Author: calamity 时间: 2021-1-25 14:45

卡开发发发表于 2021-1-25 10:16
不太确定CO初始状态是气态还是溶液。我觉得可能处理的手段是，应该开着VASPSol，这样即便是clean的表面也把 ...

谢谢，我去看一下文献

作者
Author: Penson 时间: 2021-1-25 15:07

calamity 发表于 2021-1-25 14:41
但是原来的方法对slab上水的构型要求太大了，然后跑AIMD感觉不太来得及。后来我想了一下，可以用第二种方 ...

没太理解，又看了一遍你说的第二种方法，好像意思是不考虑slab？只算CO和H2O的体系？

作者
Author: 卡开发发 时间: 2021-1-25 17:00

calamity 发表于 2021-1-25 14:40
当前计算能力只能算4个水了，请问方法对吗？第一个还是第二个方法

4个H2O应该配合隐式溶剂使用（甚至用不了4个），作为隐式溶剂不好描述的部分进行讨论就行（比如氢键或其他H2O显式参加的相互作用），如果没有这种情况甚至直接考虑使用隐式溶剂就行。在处理各个结构时应该各自优化一下。

作者
Author: calamity 时间: 2021-1-26 17:08

Penson 发表于 2021-1-25 15:07
没太理解，又看了一遍你说的第二种方法，好像意思是不考虑slab？只算CO和H2O的体系？

不是，考虑了slab，在slab上吸附4个H2O分子以及一个CO+4个H2O分子

作者
Author: calamity 时间: 2021-1-26 17:09

卡开发发发表于 2021-1-25 17:00
4个H2O应该配合隐式溶剂使用（甚至用不了4个），作为隐式溶剂不好描述的部分进行讨论就行（比如氢键或其 ...

加了隐式溶剂模型，再加了4个H2O分子

作者
Author: 奥利给 时间: 2023-10-17 09:02
请问你的问题解决了没？可否给俺说一下

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