计算化学公社

标题: MOF-5连接单元的量化计算收敛问题 [打印本页]

作者
Author: Wenchuan 时间: 2024-9-30 08:51
标题: MOF-5连接单元的量化计算收敛问题
本帖最后由 Wenchuan 于 2024-9-30 10:37 编辑

1、我想算一个MOF-5连接单元的单点能，需要先对其优化，优化完成后频率分析发现Displacement出现no，尝试了Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)中的多种方法发现只能使用calcall，但是使用calcall后反而不收敛了，输入文件如下： (, 下载次数 Times of downloads: 9) 。尝试opt=(calcall,gdiis)后收敛情况也没改善。应该怎么优化才能得到更精确结果呢？
2、最开始使用近似Hessian矩阵优化后频率计算结果为
      Item             Value    Threshold  Converged?
Maximum Force          0.000229    0.000450    YES
RMS    Force          0.000032    0.000300    YES
Maximum Displacement    0.056907    0.001800    NO
RMS    Displacement    0.018370    0.001200    NO
没有虚频
这个结果可以用吗？

作者
Author: lemon_electron 时间: 2024-9-30 09:28
本帖最后由 lemon_electron 于 2024-9-30 09:31 编辑

之前优化过比这个配体大的，收敛还挺快，不过我选的泛函比较简单。。学习一下

作者
Author: sobereva 时间: 2024-10-1 03:31
“只能使用calcall” 语义不明，怎么“只能”了？博文里该说的都说了，没有任何遗漏
上传振动分析输入输出文件
当前问题没有带D3的必要，否则还得自定义D3参数，增加了使用不当导致问题的风险
确保判断出了实际的点群

作者
Author: Wenchuan 时间: 2024-10-1 19:23

sobereva 发表于 2024-10-1 03:31
“只能使用calcall” 语义不明，怎么“只能”了？博文里该说的都说了，没有任何遗漏
上传振动分析输入输出 ...

不好意思，不是只能使用calcall，是尝试了其它方法比如使用更严收敛限或使用opt=recalc=n等都没有解决问题。
1. 加DFT-D3是因为需要研究气体与MOF的弱相互作用，也计算了气体吸附的结构，这是需要加D3的，而为了计算方法的一致性，用于计算吸附能和对比的纯MOF结构也加了D3，是不是没有与气体的相互作用就应该不加D3呢？
气体吸附的结构虽然增加了气体分子，但使用calcall后结构却更容易收敛。只是纯MOF结构使用calcall后很长时间一直不收敛。
2. 不加D3后优化的输入输出文件为 (, 下载次数 Times of downloads: 4) ，频率分析的输入输出文件为 (, 下载次数 Times of downloads: 2) ，频率分析位移结果出现no，但无虚频，这可以用吗？

作者
Author: sobereva 时间: 2024-10-1 23:46

Wenchuan 发表于 2024-10-1 19:23
不好意思，不是只能使用calcall，是尝试了其它方法比如使用更严收敛限或使用opt=recalc=n等 ...

可以用
然而你这个体系明显是Td点群，显然应当先在gview里做对称化之后再做计算，这是做高斯计算最基本的常识
D3不是不能加，但你若用了IOp自定义了D3参数，必须记得freq的时候也得用
既然都试了calcall，显然试recalc是完全多余的

作者
Author: Wenchuan 时间: 2024-10-2 22:17

sobereva 发表于 2024-10-1 23:46
可以用
然而你这个体系明显是Td点群，显然应当先在gview里做对称化之后再做计算，这是做高斯计算最基本 ...

感谢您的解答，使用Td点群后频率分析中的位移变成了yes，问题已经解决了。

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