H2和SiCl4双分子体系弱相互作用优化，加弥散和色散矫正后运行一晚上没有结束

healthy · 发表于 Post on 2023-7-5 08:59:14

我想用multiwfn计算两个分子的弱相互作用力和psi4计算能量分解，首先想要得到一个优化后的双分子体系。
我看大家都说计算弱相互作用要加弥散和色散矫正，我看了文章，不知道加的对不对，结果运行了一晚上还卡在一开始的地方。
高斯版本是09W。

输入文件：%chk=D:\Files\Gaussianfiles\7.3\HCl-4\7.41.chk# opt freq B3LYP/6-31g** em=GD3 maxdisk=1500MB

运行后一晚上都卡在输出文件（附件里）里面那个位置。

sobereva · 发表于 Post on 2023-7-5 09:20:23

说明用的CPU
注意下文里说的卡住的问题
Gaussian的安装方法及运行时的相关问题
http://sobereva.com/439（http://bbs.keinsci.com/thread-10814-1-1.html）

你并没有加弥散
做优化用6-311G**足矣，不需要带弥散

healthy · 发表于 Post on 2023-7-5 09:45:32

sobereva 发表于 2023-7-5 09:20
说明用的CPU
注意下文里说的卡住的问题
Gaussian的安装方法及运行时的相关问题

感谢老师的回答。

CPU内核16，逻辑处理器32

电脑里没有自己下载的杀毒软件

%chk=D:\Files\Gaussianfiles\7.3\H2-4.chk
# opt freq B3LYP/6-311g** em=GD3 maxdisk=1400MB
改过之后重新运行也是进行到之前的地方就不动了
我要继续等待吗，还是存在其他什么问题呢

kimariyb · 发表于 Post on 2023-7-5 09:50:42

本帖最后由 kimariyb 于 2023-7-5 10:02 编辑

把 maxdisk=1400MB 删了试试，不知道你为什么要用这个关键词。我把这个删了，随便跑，不过你这个体系几何优化可能会难收敛，我直接拿你的 out 文件在 b3lyp(d3)/6-311g** 下做 opt 和 freq，收敛失败。如果你遇到了，最好把初猜结构摆好一点，或者在《量子化学计算中帮助几何优化收敛的常用方法》（http://sobereva.com/164）里需找解决方法。

sobereva · 发表于 Post on 2023-7-5 10:08:47

healthy 发表于 2023-7-5 09:45
感谢老师的回答。

CPU内核16，逻辑处理器32

用linux版，或者64bit windows版（不建议用windows版，超过10核后效率极低）

healthy · 发表于 Post on 2023-7-5 19:30:52

sobereva 发表于 2023-7-5 10:08
用linux版，或者64bit windows版（不建议用windows版，超过10核后效率极低）

老师，我换了linux  G16版，输入文件为：
%mem=16GB%nprocshared=16
%chk=7.chk
# opt  freq B3LYP/6-311g** em=GD3
可以运行出来了。

然后我换成了H2和SiHCl3体系，用同样的方法却报错了，错误如下：
Error termination request processed by link 9999.
Optimization stopped.
-- Number of steps exceeded,  NStep=28
-- Flag reset to prevent archiving.下面是我用Gview看的输出文件

然后我根据几何优化常见收敛方法那篇文章，因为我这个体系是需要跟H2-SiCl4对比的，所以基组什么的需要保持一致，然后我就加了几个关键词，改成了这样：
%mem=16GB
%nprocshared=16
%chk=H2_3.chk
# opt=(gdiis,maxstep=3,notrust,maxcyc=200) freq B3LYP/6-311g** em=GD3

但还是报错了同样的错误，请问是为什么呢

sobereva · 发表于 Post on 2023-7-6 09:31:37

healthy 发表于 2023-7-5 19:30
老师，我换了linux G16版，输入文件为：
%mem=16GB%nprocshared=16
%chk=7.chk

没输出文件没法说

此文里所有方法都考虑了，还不行再问，别就光试了几个不行就急着问
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

healthy · 发表于 Post on 2023-7-6 09:39:31

sobereva 发表于 2023-7-6 09:31
没输出文件没法说

此文里所有方法都考虑了，还不行再问，别就光试了几个不行就急着问

好的老师，我再试试，感谢。

healthy · 发表于 Post on 2023-7-6 09:40:04

kimariyb 发表于 2023-7-5 09:50
把 maxdisk=1400MB 删了试试，不知道你为什么要用这个关键词。我把这个删了，随便跑，不过你这个体系几何优 ...

好的我尝试一下，感谢。

healthy · 发表于 Post on 2023-7-7 15:42:15

sobereva 发表于 2023-7-5 09:20
说明用的CPU
注意下文里说的卡住的问题
Gaussian的安装方法及运行时的相关问题

请问老师，我用这个优化的结果做弱相互作用分析也不需要加弥散吗

sobereva · 发表于 Post on 2023-7-8 09:23:06

healthy 发表于 2023-7-7 15:42
请问老师，我用这个优化的结果做弱相互作用分析也不需要加弥散吗

说清楚什么方法做什么分析

healthy · 发表于 Post on 2023-7-8 14:50:03

sobereva 发表于 2023-7-8 09:23
说清楚什么方法做什么分析

优化H2和SiCl4、HCl和SiCl4两个双分子体系，用multiwfn的iri方法做分子间弱相互作用分析，用psi4做相互作用能分解，我用Gaussian g16优化的时候加了色散矫正，请问需要加弥散吗

wzkchem5 · 发表于 Post on 2023-7-8 15:02:21

healthy 发表于 2023-7-8 07:50
优化H2和SiCl4、HCl和SiCl4两个双分子体系，用multiwfn的iri方法做分子间弱相互作用分析，用psi4做相互作 ...

算相互作用能应该加弥散
IRI方法不需要加弥散，但加了也不算错

healthy · 发表于 Post on 2023-7-8 20:31:36

wzkchem5 发表于 2023-7-8 15:02
算相互作用能应该加弥散
IRI方法不需要加弥散，但加了也不算错

好的，谢谢

sobereva · 发表于 Post on 2023-7-9 09:28:37

healthy 发表于 2023-7-8 14:50
优化H2和SiCl4、HCl和SiCl4两个双分子体系，用multiwfn的iri方法做分子间弱相互作用分析，用psi4做相互作 ...

这取决于你原本用什么基组
诸如def2-QZVP哪怕算相互作用能都不需要非得加弥散

鉴于你的体系非常小，def2-TZVP优化和做波函数分析，def2-QZVPP算相互作用能就可以（或者ma-def2-TZVP也行）。算相互作用能时考虑counterpoise更好，但不是必须。

PSI4做SAPT适合用的级别在下文有说明，方法和基组搭配合适才能让误差抵消比较理想，推荐的级别都是带弥散函数的
使用PSI4做对称匹配微扰理论(SAPT)能量分解计算
http://sobereva.com/526（http://bbs.keinsci.com/thread-15902-1-1.html）

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

[新手求助] H2和SiCl4双分子体系弱相互作用优化，加弥散和色散矫正后运行一晚上没有结束