计算化学公社

标题: 沸石吸附Zn离子的吸附能过大 [打印本页]

作者
Author: DILIGENT 时间: 2023-7-10 00:29
标题: 沸石吸附Zn离子的吸附能过大
我计算了沸石吸附Zn2+的吸附能，结果如下，吸附能为-2000kj/mol。。。。
我上传了我的输入文件和输出文件。

作者
Author: wzkchem5 时间: 2023-7-10 04:45
Zn2+必须做显式溶剂化

作者
Author: sobereva 时间: 2023-7-10 08:55
-2000kj/mol不一定就是错的，毕竟形成那么多Zn-O键
算吸附能的时候建议改为TZVP-MOLOPT-GTH（对于有定义的元素来说）以提升精度

作者
Author: DILIGENT 时间: 2023-7-11 12:07

sobereva 发表于 2023-7-10 08:55
-2000kj/mol不一定就是错的，毕竟形成那么多Zn-O键
算吸附能的时候建议改为TZVP-MOLOPT-GTH（对于有定义的 ...

谢谢卢老师，还想请教下您计算吸附能的几个问题：1、将优化完复合物结构中的沸石、Zn结构提取出来用TZVP-MOLOPT-GTH计算单点能，然后用复合物能量 - 沸石能量 - zn能量的方式算结合能。2 、分别用DZVP-MOLOPT-SR-GTH优化完结构后再用TZVP-MOLOPT-GTH计算单点能用于计算吸附能。
以上两种方式计算结合强弱都可以吗。
3、用PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH计算MOF吸附CO2体系（OT），用PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH计算氧化铜吸附CO2体系（对角化），然后用TZVP-MOLOPT-GTH计算单点能用于吸附能的计算；这样得到两个体系的吸附能可以用于比较吗

作者
Author: DILIGENT 时间: 2023-7-11 12:12

wzkchem5 发表于 2023-7-10 04:45
Zn2+必须做显式溶剂化

谢谢老师。

作者
Author: sobereva 时间: 2023-7-11 12:43

DILIGENT 发表于 2023-7-11 12:07
谢谢卢老师，还想请教下您计算吸附能的几个问题：1、将优化完复合物结构中的沸石、Zn结构提取出来用TZVP- ...

没看懂
优化一般用DZVP-MOLOPT-SR-GTH就够了。基于优化完的结构，对整体、Zn、沸石再都用TZVP-MOLOPT-GTH算单点，按你的公式求差算结合能。

OT和对角化得到的能量可以直接对比（假设都收敛到了稳定波函数）

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