计算化学公社

标题: 溶剂化后出现Error on total polarization charges问题，请问大神如何处理？ [打印本页]

作者
Author: lastzealot 时间: 2015-4-16 17:34
标题: 溶剂化后出现Error on total polarization charges问题，请问大神如何处理？
Error on total polarization charges =  0.05638
SCF Done:  E(UB3LYP) =  -555.504532264    A.U. after 2001 cycles
         NFock=***  Conv=0.22D-03    -V/T= 2.1273
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 1.0000 <S**2>= 2.0239 S= 1.0079
<L.S>= 0.000000000000E+00
KE= 4.927864503262D+02 PE=-2.267216731101D+03 EE= 7.163955074496D+02
SMD-CDS (non-electrostatic) energy    (kcal/mol) =    1.20
(included in total energy above)
Annihilation of the first spin contaminant:
S**2 before annihilation    2.0239, after    2.0003
Convergence failure -- run terminated.
Error termination via Lnk1e in /db/home/renym/g09/l502.exe at Thu Apr 16 18:45:52 2015.
Job cpu time:    0 days 12 hours 29 minutes 16.5 seconds.
File lengths (MBytes):  RWF= 230 Int=    0 D2E=    0 Chk=    10 Scr=    2

这是输出结果。气相计算没有问题，但是加了水溶剂就出现问题。结果文件中的能量和RMS都是单调下降。
请问怎么处理？
谢谢。

作者
Author: sobereva 时间: 2015-4-18 04:32
还是普通的SCF不收敛问题，虽然是溶剂下，但处理方式和气相下一样，见
解决SCF不收敛问题的方法
http://sobereva.com/61

步数上限设到2000不会解决任何问题

也可以尝试改为scrf=CPCM看是否恰好能收敛。

作者
Author: lastzealot 时间: 2015-4-20 10:33

sobereva 发表于 2015-4-18 04:32
还是普通的SCF不收敛问题，虽然是溶剂下，但处理方式和气相下一样，见
解决SCF不收敛问题的方法
http://h ...

大神，那个解决SCF不收敛的方法有很多，请问是按照优先的顺序写的吗？我如何来选择用哪种方法呢？
比如，先尝试哪一种，然后...

作者
Author: sobereva 时间: 2015-4-20 12:58
没有万能的方法，应在理解原理的基础上根据实际情况进行尝试

作者
Author: laoman 时间: 2017-8-29 16:19

sobereva 发表于 2015-4-20 12:58
没有万能的方法，应在理解原理的基础上根据实际情况进行尝试

请问Sob老师，我的几个加溶剂的体系也都出现了这样的提示，只有两个体系是出现了虚频，其他体系没有虚频，且freq之后也是4个YES，这种情况是不是重算出现虚频的两个体系就可以？

作者
Author: sobereva 时间: 2017-8-29 20:57

laoman 发表于 2017-8-29 16:19
请问Sob老师，我的几个加溶剂的体系也都出现了这样的提示，只有两个体系是出现了虚频，其他体系没有虚频 ...

不知你说的是哪个提示。如果是指Error on total polarization charges，这根本就不是报错
有虚频的结构肯定没法用来算极小点结构下的性质，必须把虚频搞掉。有虚频时即便都四个YES也没意义，此时有N个虚频说明这个结构是N阶鞍点。

作者
Author: laoman 时间: 2017-8-31 15:17

sobereva 发表于 2017-8-29 20:57
不知你说的是哪个提示。如果是指Error on total polarization charges，这根本就不是报错
有虚频的结 ...

对的，就是所有体系的log文件在SCF的时候都有Error on total polarization charges的提示，原来不算是报错的意思啊。
两个出现虚频的体系还在尝试用这个帖子提到的方法继续优化。。。
http://sobereva.com/61

作者
Author: avatar_man 时间: 2021-8-18 20:20
你好，请问你的这个问题最终解决了吗？怎么解决的？我的也碰到了类似的情况，谢谢

作者
Author: zhaojie1986 时间: 2023-5-27 09:53

同问，请问你的这个问题最终解决了吗？怎么解决的？我的也碰到了类似的情况，谢谢

作者
Author: wzkchem5 时间: 2023-5-27 15:04

zhaojie1986 发表于 2023-5-27 02:53
同问，请问你的这个问题最终解决了吗？怎么解决的？我的也碰到了类似的情况，谢谢

Error on total polarization charges的error是误差的意思，不是错误的意思。看输出文件的时候必须对最基本的英语一词多义现象有概念，不要看到error就以为是错误，即便程序确实报错退出了，那也可能是因为另一个原因报错退出的，而描述那个原因的话里未必有error这个词。
真正报错信息是Convergence failure -- run terminated. 按SCF不收敛的解决方法解决即可。

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