计算化学公社

 找回密码 Forget password
 注册 Register
Views: 26474|回复 Reply: 9
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 溶剂化后出现Error on total polarization charges问题,请问大神如何处理?

[复制链接 Copy URL]

720

帖子

0

威望

1252

eV
积分
1972

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
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都是单调下降。
请问怎么处理?
谢谢。

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

2#
发表于 Post on 2015-4-18 04:32:45 | 只看该作者 Only view this author
还是普通的SCF不收敛问题,虽然是溶剂下,但处理方式和气相下一样,见
解决SCF不收敛问题的方法
http://sobereva.com/61

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

也可以尝试改为scrf=CPCM看是否恰好能收敛。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

720

帖子

0

威望

1252

eV
积分
1972

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2015-4-20 10:33:42 | 只看该作者 Only view this author
sobereva 发表于 2015-4-18 04:32
还是普通的SCF不收敛问题,虽然是溶剂下,但处理方式和气相下一样,见
解决SCF不收敛问题的方法
http://h ...

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

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

4#
发表于 Post on 2015-4-20 12:58:59 | 只看该作者 Only view this author
没有万能的方法,应在理解原理的基础上根据实际情况进行尝试
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

143

帖子

3

威望

4430

eV
积分
4633

Level 6 (一方通行)

5#
发表于 Post on 2017-8-29 16:19:08 | 只看该作者 Only view this author
sobereva 发表于 2015-4-20 12:58
没有万能的方法,应在理解原理的基础上根据实际情况进行尝试

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

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

6#
发表于 Post on 2017-8-29 20:57:42 | 只看该作者 Only view this author
laoman 发表于 2017-8-29 16:19
请问Sob老师,我的几个加溶剂的体系也都出现了这样的提示,只有两个体系是出现了虚频,其他体系没有虚频 ...


不知你说的是哪个提示。如果是指Error on total polarization charges,这根本就不是报错
有虚频的结构肯定没法用来算极小点结构下的性质,必须把虚频搞掉。有虚频时即便都四个YES也没意义,此时有N个虚频说明这个结构是N阶鞍点。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

143

帖子

3

威望

4430

eV
积分
4633

Level 6 (一方通行)

7#
发表于 Post on 2017-8-31 15:17:07 | 只看该作者 Only view this author
sobereva 发表于 2017-8-29 20:57
不知你说的是哪个提示。如果是指Error on total polarization charges,这根本就不是报错
有虚频的结 ...

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

36

帖子

0

威望

401

eV
积分
437

Level 3 能力者

8#
发表于 Post on 2021-8-18 20:20:33 | 只看该作者 Only view this author
你好,请问你的这个问题最终解决了吗?怎么解决的?我的也碰到了类似的情况,谢谢

8

帖子

0

威望

31

eV
积分
39

Level 2 能力者

9#
发表于 Post on 2023-5-27 09:53:04 | 只看该作者 Only view this author

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

1万

帖子

0

威望

9855

eV
积分
22090

Level 6 (一方通行)

10#
发表于 Post on 2023-5-27 15:04:29 | 只看该作者 Only view this author
zhaojie1986 发表于 2023-5-27 02:53
同问,请问你的这个问题最终解决了吗?怎么解决的?我的也碰到了类似的情况,谢谢

Error on total polarization charges的error是误差的意思,不是错误的意思。看输出文件的时候必须对最基本的英语一词多义现象有概念,不要看到error就以为是错误,即便程序确实报错退出了,那也可能是因为另一个原因报错退出的,而描述那个原因的话里未必有error这个词。
真正报错信息是Convergence failure -- run terminated. 按SCF不收敛的解决方法解决即可。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-19 05:35 , Processed in 0.161153 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list