计算化学公社

标题: 请教溶剂化计算出错问题 [打印本页]

作者
Author: keshsong 时间: 2015-11-18 14:31
标题: 请教溶剂化计算出错问题
输入文件：
#p b3lyp/6-31++g(d,p) scrf(pcm,read,dovacuum,externaliteration)

输出文件（最后部分）：

            Baa -0.0062    0.451    0.161    0.150  0.7141  0.1210 -0.6896
   9 O(17)  Bbb -0.0041    0.295    0.105    0.098  0.6816 -0.3449  0.6453
            Bcc    0.0103 -0.746 -0.266 -0.249  0.1598  0.9308  0.3287

            Baa -0.0087    0.628    0.224    0.209  0.7135 -0.1794  0.6773
10 O(17)  Bbb -0.0032    0.234    0.083    0.078  0.6628  0.4861 -0.5695
            Bcc    0.0119 -0.861 -0.307 -0.287 -0.2271  0.8553  0.4657

---------------------------------------------------------------------------------

No NMR shielding tensors so no spin-rotation constants.
Leave Link  601 at Thu Jul 24 17:50:40 2014, MaxMem= 1073741824 cpu:    3.6
(Enter /opt//g09/l124.exe)
SC-PCM: Using the SCF density.
SC-PCM: The equilibrium contribution to the free energy is -0.068199183 a.u.
SC-PCM: The total energy including the PCM contribution is    -342.461159696 a.u.
SC-PCM: Erms=  0.52528D-05 Emax=  0.17248D-04 Qconv=  0.10000D-06
SC-PCM: Iterating to achieve self-consistent PCM reaction field, IExtIt=  50.
SC-PCM: Convergence not achieved within  50 iterations ... aborting job.
Error termination via Lnk1e in /opt//g09/l124.exe at Thu Jul 24 17:50:41 2014.
Job cpu time:  0 days  0 hours  5 minutes 36.1 seconds.
File lengths (MBytes):  RWF=    35 Int=    0 D2E=    0 Chk=    3 Scr=    1

请问这个问题怎么解决啊？谢谢

作者
Author: keshsong 时间: 2015-11-18 14:55
请大神给点建议或者相关文献。

作者
Author: sobereva 时间: 2015-11-18 23:23
去掉弥散，弥撒函数总会阻碍收敛。
不要写externaliteration，很多人不懂外迭代是什么意思却总是乱用。这仅仅仅仅对于后HF计算才有意义！对于DFT基态计算毫无意义！所以去掉externaliteration

你要计算溶剂化能，最简单最好的办法就是用M052X/6-31G* scrf=smd的单点能减去M052X/6-31G*的单点能。不要用dovacuum。

作者
Author: keshsong 时间: 2015-11-19 08:36
本帖最后由 keshsong 于 2015-11-19 09:34 编辑

sobereva 发表于 2015-11-18 23:23
去掉弥散，弥撒函数总会阻碍收敛。
不要写externaliteration，很多人不懂外迭代是什么意思却总是乱用。这 ...

谢谢你的详细解答，还有如下问题请指教。
1. 如果优化用了弥散函数，溶剂化计算不用，对结果是否有影响？
抑或优化也就不用弥散函数？因为我的体系是负离子，所以才使用弥散函数进行优化。
2. 因为想给出G03的结果模式，才根据手册选用externaliteration和dovacuum.
你的建议是否认为上述2个参数在非post-HF的scrf中没有使用的必要。

谢谢！

作者
Author: sobereva 时间: 2015-11-19 11:25

keshsong 发表于 2015-11-19 08:36
谢谢你的详细解答，还有如下问题请指教。
1. 如果优化用了弥散函数，溶剂化计算不用，对结果是否有影 ...

1 无影响。阴离子看具体负电荷怎么分布的，明显分布在一处就给那部分用弥散就完了，如果是均摊开，远不如用3-zeta基组。氢的弥散几乎总是多余的，除了诸如LiH这样H带明显负电荷的情况。
2 不知道你为什么非得得到03的结果。直接用我说的SMD的方法直接就得到溶解自由能了，比起PCM更准。

作者
Author: keshsong 时间: 2015-11-19 15:24

sobereva 发表于 2015-11-19 11:25
1 无影响。阴离子看具体负电荷怎么分布的，明显分布在一处就给那部分用弥散就完了，如果是均摊开，远不如 ...

非常感谢，明白。

作者
Author: lei234 时间: 2016-4-2 08:08

sobereva 发表于 2015-11-18 23:23
去掉弥散，弥撒函数总会阻碍收敛。
不要写externaliteration，很多人不懂外迭代是什么意思却总是乱用。这 ...

我想请教您算溶剂化能用的基组和方法可以和之前优化或者单点计算的方法和基组不一样吗？就是不管当前的体系中用的是什么方法和基组，只要是算溶剂化能就用M052X/6-31G* scrf=smd的单点能减去M052X/6-31G*的单点能吗？@sobereva

作者
Author: sobereva 时间: 2016-4-2 12:46

lei234 发表于 2016-4-2 08:08
我想请教您算溶剂化能用的基组和方法可以和之前优化或者单点计算的方法和基组不一样吗？就是不管当前的体 ...

可以，本来也应当这样

作者
Author: lei234 时间: 2016-4-3 06:34

sobereva 发表于 2016-4-2 12:46
可以，本来也应当这样

知道了，谢谢~

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