计算化学公社

标题: 计算一个UiO-66簇与阴离子之间的bsse相互作用，出现问题 [打印本页]

作者
Author: zyj19831206 时间: 2024-9-5 22:40
标题: 计算一个UiO-66簇与阴离子之间的bsse相互作用，出现问题
计算UiO-66簇与高铬酸根之间的吸附能，按照gjf文件那么设置泛函和基组，SCF不收敛，应该怎么修改呢？

作者
Author: imasen 时间: 2024-9-5 23:24
本帖最后由 imasen 于 2024-9-5 23:27 编辑

SCF不收敛可以参照卢天老师这篇博文解决（解决SCF不收敛问题的方法 http://sobereva.com/61）。
但在此之前你的结构优化了么？UiO66的命名是Hexakis[μ-[1,4-benzenedicarboxylato(2-)-κO1:κO'1]]tetra-μ3-hydroxytetra-μ3-oxohexazirconium，这意味着桥连的8个氧中有四个氧是氢氧根离子，有四个氧是氧负离子，你这个结构是有问题的，好像你为了平衡电荷自己删除了好几个氧。

作者
Author: sobereva 时间: 2024-9-5 23:37
想算准键能，区区6-31G*根本不够。断键附近的原子应该都用def2-TZVP

算化学键作用根本就不该用counterpoise，仔细看下文
计算化学键键能时以counterpoise方式考虑BSSE不仅是多余的甚至是有害的
http://sobereva.com/381（http://bbs.keinsci.com/thread-6047-1-1.html）

guess=save完全多余

用PBE-D3(BJ)毫无好处，改成PBE0-D3(BJ)是更好选择

并且确保结构有意义。其它的，该怎么解决怎么解决
解决SCF不收敛问题的方法
http://sobereva.com/61

作者
Author: zyj19831206 时间: 2024-9-6 09:07

sobereva 发表于 2024-9-5 23:37
想算准键能，区区6-31G*根本不够。断键附近的原子应该都用def2-TZVP

算化学键作用根本就不该用counterpo ...

收到，社长，好久没做高斯了，忘记了一些东西，是要及时复习，

作者
Author: zyj19831206 时间: 2024-9-6 09:09

sobereva 发表于 2024-9-5 23:37
想算准键能，区区6-31G*根本不够。断键附近的原子应该都用def2-TZVP

算化学键作用根本就不该用counterpo ...

主要还是想做这个阴离子与MOF簇之间的吸附能，也不算完全的键能。

作者
Author: sobereva 时间: 2024-9-6 09:30

zyj19831206 发表于 2024-9-6 09:09
主要还是想做这个阴离子与MOF簇之间的吸附能，也不算完全的键能。

当前是化学吸附，牵扯到断键，显然核心就是算键能

作者
Author: zyj19831206 时间: 2024-9-6 09:38

imasen 发表于 2024-9-5 23:24
SCF不收敛可以参照卢天老师这篇博文解决（解决SCF不收敛问题的方法 http://sobereva.com/61）。
但在此之 ...

我没有删除，是从CIF结构直接截取出来的。

作者
Author: zyj19831206 时间: 2024-9-6 09:41

sobereva 发表于 2024-9-6 09:30
当前是化学吸附，牵扯到断键，显然核心就是算键能

收到

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