MOF-5连接单元的量化计算收敛问题

Wenchuan · 发表于 Post on 2024-9-30 08:51:23

本帖最后由 Wenchuan 于 2024-9-30 10:37 编辑

1、我想算一个MOF-5连接单元的单点能，需要先对其优化，优化完成后频率分析发现Displacement出现no，尝试了Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)中的多种方法发现只能使用calcall，但是使用calcall后反而不收敛了，输入文件如下：

mix.gjf (5.43 KB, 下载次数 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
没有虚频
这个结果可以用吗？

lemon_electron · 发表于 Post on 2024-9-30 09:28:33

本帖最后由 lemon_electron 于 2024-9-30 09:31 编辑

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

sobereva · 发表于 Post on 2024-10-1 03:31:59

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

Wenchuan · 发表于 Post on 2024-10-1 19:23:00

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

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

opt.rar (750.98 KB, 下载次数 Times of downloads: 4) ，频率分析的输入输出文件为

freq.rar (240.64 KB, 下载次数 Times of downloads: 2) ，频率分析位移结果出现no，但无虚频，这可以用吗？

sobereva · 发表于 Post on 2024-10-1 23:46:43

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

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

Wenchuan · 发表于 Post on 2024-10-2 22:17:43

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] MOF-5连接单元的量化计算收敛问题

浏览过的版块