sobereva 发表于 2023-7-11 10:37 好的老师,抱歉老师,发第一条的时候跨页了,我没看见,我以为没发出来,就重新发了一个,后来发现之后重复的那个删不掉了,以后注意 |
healthy 发表于 2023-7-10 14:17 如果之前帖子内容有问题,直接重新编辑,不要再发一个。 QZ级别的基组自然非常贵,做优化opt freq极其不划算 仔细看 浅谈为什么优化和振动分析不需要用大基组 http://sobereva.com/387(http://bbs.keinsci.com/thread-6600-1-1.html) 谈谈量子化学中基组的选择 http://sobereva.com/336(http://bbs.keinsci.com/thread-3545-1-1.html) 别给氢加弥散,纯粹白浪费时间,此文说了 谈谈弥散函数和“月份”基组 http://sobereva.com/119 中性分子优化目的6-311G**足矣,重原子都不需要带弥散。此基组下的波函数做IRI分析完全够了,优化出的结构做SAPT计算也够了。 若无特殊情况,B3LYP一律结合DFT-D3(BJ) |
sobereva 发表于 2023-7-9 09:28 老师我尝试了# opt freq B3LYP/def2QZVPP em=GD3 运行时间很长,然后一直不收敛,我通过之前的文章尝试了很多收敛方法到现在也没有收敛 但我之前用# opt freq B3LYP/6-311+G** em=GD3或者# opt freq B3LYP/6-311++G** em=GD3计算收敛就很快,但就是不知道这个基组是否适合我用psi4 SAPT2+(3)δMP2/aug-cc-pVTZ级别做相互作用能分解,和用multiwfn做iri弱相互作用分析? |
sobereva 发表于 2023-7-9 09:28 老师,我用opt freq B3LYP/def2QZVPP em=GD3计算,最后不收敛,我又加了opt=recalc=3,结果算了三四个小时了还没算完,不知道怎么回事。 请问我能不能用opt freq B3LYP/6-311+G** em=GD3或者opt freq B3LYP/6-311++G** em=GD3来算呢,我之前用这个很快出结果了,就是不知道这个基组对我用psi4 SAPT2+(3)δMP2/aug-cc-pVTZ级别算相互作用能和用multiwfn算弱相互作用力来说合不合适? |
healthy 发表于 2023-7-8 14:50 这取决于你原本用什么基组 诸如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) |
wzkchem5 发表于 2023-7-8 15:02 好的,谢谢 |
healthy 发表于 2023-7-8 07:50 算相互作用能应该加弥散 IRI方法不需要加弥散,但加了也不算错 |
sobereva 发表于 2023-7-8 09:23 优化H2和SiCl4、HCl和SiCl4两个双分子体系,用multiwfn的iri方法做分子间弱相互作用分析,用psi4做相互作用能分解,我用Gaussian g16优化的时候加了色散矫正,请问需要加弥散吗 |
healthy 发表于 2023-7-7 15:42 说清楚什么方法做什么分析 |
sobereva 发表于 2023-7-5 09:20 请问老师,我用这个优化的结果做弱相互作用分析也不需要加弥散吗 |
kimariyb 发表于 2023-7-5 09:50 好的我尝试一下,感谢。 |
sobereva 发表于 2023-7-6 09:31 好的老师,我再试试,感谢。 |
healthy 发表于 2023-7-5 19:30 没输出文件没法说 此文里所有方法都考虑了,还不行再问,别就光试了几个不行就急着问 量子化学计算中帮助几何优化收敛的常用方法 http://sobereva.com/164 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-27 02:19 , Processed in 0.248570 second(s), 26 queries , Gzip On.