计算化学公社

标题: 阴离子复合物体系不加弥散函数算出来的热力学校正值可以用吗？ [打印本页]

作者
Author: Labber 时间: 2022-2-19 11:54
标题: 阴离子复合物体系不加弥散函数算出来的热力学校正值可以用吗？
本帖最后由 Labber 于 2022-2-19 16:35 编辑

150个原子以上，带一个单位负电荷的阴离子复合物体系，想在小基组下算freq得到热力学校正值，再用高精度的算SP，最后求G。
该体系在B3LYP-D3(BJ)/6-311G**下做opt超算一天以内可以算完，但是加了弥散函数B3LYP-D3(BJ)/6-311+G**之后，即使guess=read，geom=allcheck之后依然收敛的很慢，耗费机时很长，超算24核算了一个晚上opt才走了两步，SCF每步要走25圈左右。
所以请问这样的阴离子体系，不加弥散函数，freq后的热力学校正值能直接使用吗？高精度的SP打算用M062X/ma-TZVp，估计耗时也会很长

请各位大佬指教

作者
Author: gauss98 时间: 2022-2-19 12:00
你用了calcfc
时间接近于算一个频率
相当于20步opt左右
估计你的一晚上大部分时间都用于这个了
再走几步opt
谨慎估计一小时一步左右

作者
Author: Labber 时间: 2022-2-19 12:07

gauss98 发表于 2022-2-19 12:00
你用了calcfc
时间接近于算一个频率
相当于20步opt左右

不是的，一开始算的时候没有加，用的以下这个
#p B3LYP/6-311+G** em=GD3BJ opt=GDIIS freq SCF(novaracc,noincfock) int=acc2e=12
但是还是算了一个晚上

后来我才改回主题所述的条件

所以我还是想知道阴离子体系，不加弥散函数的热力学校正值是否能用呢？

作者
Author: gauss98 时间: 2022-2-19 12:18
按照社长的说法
频率是不需要太大基组的
加上校正因子即可了
150个原子6311g**不低了

作者
Author: sobereva 时间: 2022-2-19 12:33
即便阴离子体系，opt freq的目的用3-zeta基组也够了，顶多是再给负电荷极度富集的原子加弥散。全加弥散又浪费时间又容易不收敛

作者
Author: wzkchem5 时间: 2022-2-19 16:06
阴离子体系并不是必然需要弥散函数，如果负电荷足够弥散，可能不需要加。如果不给具体结构，只告诉我们150个原子、带一个负电荷，条件不足，我们没法判断。

作者
Author: Labber 时间: 2022-2-19 16:36

wzkchem5 发表于 2022-2-19 16:06
阴离子体系并不是必然需要弥散函数，如果负电荷足够弥散，可能不需要加。如果不给具体结构，只告诉我们150 ...

您好，附件log已重新上传，主要是环糊精和对硝基苯甲酸阴离子之间，有氢键作用。

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