计算化学公社

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

[新手求助] 频率分析失败( Inaccurate quadrature in CalDSu.)

[复制链接 Copy URL]

76

帖子

0

威望

863

eV
积分
939

Level 4 (黑子)

请教各位老师:

对优化收敛的结构进行频率分析,最后失败结束,原因是:
Inaccurate quadrature in CalDSu.
Error termination via Lnk1e in /home/gauss/Gaussian16_AVX/tar/g16/l502.exe at Wed Nov  9 18:39:35 2022.

按照论坛上的帖子尝试了加入关键词:int=superfine scf=(novaracc,noincfock) guess=INDO 但是还是出现这个错误。输入文件关键词:#p int=superfine b3lyp/6-31g(d) em=gd3bj iop(7/127=-99,8/117=-99) scf=(novaracc,noincfock) guess=read freq CPHF=grid=fine guess=INDO

已上传结构优化后的.log文件和freq计算的输入输出文件,还请老师指点一二,万分感谢!(不好意思刚发现freq输入文件多输了一遍scf=(novaracc,noincfock) ,还请忽略这个失误)

附件.rar (441.99 KB, 下载次数 Times of downloads: 4)

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

2#
发表于 Post on 2022-11-9 20:29:35 | 只看该作者 Only view this author
为什么一定要加guess=INDO?
既然结构优化正常结束了,而频率计算又是在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
本团队长期招收研究生,有意者可私信联系

63

帖子

0

威望

867

eV
积分
930

Level 4 (黑子)

3#
发表于 Post on 2022-11-9 20:31:44 | 只看该作者 Only view this author
你咋输入了俩guess;我之前也是用了INDO没用,就用了单点的初猜,你可以都试试

76

帖子

0

威望

863

eV
积分
939

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2022-11-9 20:43:16 | 只看该作者 Only view this author
本帖最后由 lr8616 于 2022-11-9 20:47 编辑
星纹c 发表于 2022-11-9 20:31
你咋输入了俩guess;我之前也是用了INDO没用,就用了单点的初猜,你可以都试试

我用的是结构优化的初猜,你是用更高级别算单点,然后用单点能的初猜?

76

帖子

0

威望

863

eV
积分
939

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-11-9 20:44:34 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-9 20:29
为什么一定要加guess=INDO?
既然结构优化正常结束了,而频率计算又是在SCF里报错的,那么可以尝试一下读 ...

感谢回复,我就是读取的结构优化的波函数作为初猜进行的频率计算

76

帖子

0

威望

863

eV
积分
939

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2022-11-9 20:46:24 | 只看该作者 Only view this author
星纹c 发表于 2022-11-9 20:31
你咋输入了俩guess;我之前也是用了INDO没用,就用了单点的初猜,你可以都试试

不好意思,想请问是guess=read 和 guess=INDO只能用一个吗?我开始只有guess=read读取结构优化的初猜波函数,也是出现这个错误,才加上的guess=INDO

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

7#
发表于 Post on 2022-11-9 20:59:05 | 只看该作者 Only view this author
lr8616 发表于 2022-11-9 13:46
不好意思,想请问是guess=read 和 guess=INDO只能用一个吗?我开始只有guess=read读取结构优化的初猜波函 ...

对,只能用一个。只要了解guess=read和guess=INDO的原理就知道为什么。
guess=INDO其实是把初猜变得比默认初猜更差了,即使偶尔能解决问题也是碰巧,不是一个特别可靠的解决方法。
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
本团队长期招收研究生,有意者可私信联系

76

帖子

0

威望

863

eV
积分
939

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2022-11-9 21:05:28 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-9 20:59
对,只能用一个。只要了解guess=read和guess=INDO的原理就知道为什么。
guess=INDO其实是把初猜变得比默 ...

我本来确实不理解,只是跟着帖子加。。看到您的回复忽然才意识到这个问题。我用guess=read出现了这个错误,如果想尝试guess=INDO应该把guess=read删掉才对。还想请教一下不知道有什么其它可以尝试的方法吗?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

9#
发表于 Post on 2022-11-9 21:17:09 | 只看该作者 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

63

帖子

0

威望

867

eV
积分
930

Level 4 (黑子)

10#
发表于 Post on 2022-11-9 21:24:16 | 只看该作者 Only view this author
lr8616 发表于 2022-11-9 20:46
不好意思,想请问是guess=read 和 guess=INDO只能用一个吗?我开始只有guess=read读取结构优化的初猜波函 ...

我也是新手,只不过之前遇到过,我是之前算超极化率,scf不收敛,三个都加上了没用,最后用了单点的波函数,基组都不一样但就是收敛了。你这种情况实在没办法就换个基组算出来试试,用他的波函数,感觉算是野路子

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

11#
发表于 Post on 2022-11-9 22:14:18 | 只看该作者 Only view this author
lr8616 发表于 2022-11-9 14:05
我本来确实不理解,只是跟着帖子加。。看到您的回复忽然才意识到这个问题。我用guess=read出现了这个错误 ...

还有一种可以尝试的方法是,先加int=noxctest收敛以后,再读这个波函数,做一个不加int=noxctest的计算。
可能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
本团队长期招收研究生,有意者可私信联系

76

帖子

0

威望

863

eV
积分
939

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2022-11-10 09:19:58 | 只看该作者 Only view this author

谢谢卢老师,现在能算下去了,感谢!

76

帖子

0

威望

863

eV
积分
939

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2022-11-10 09:20:23 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-9 22:14
还有一种可以尝试的方法是,先加int=noxctest收敛以后,再读这个波函数,做一个不加int=noxctest的计算。 ...

谢谢您!我试试看

76

帖子

0

威望

863

eV
积分
939

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2022-11-10 09:20:59 | 只看该作者 Only view this author
星纹c 发表于 2022-11-9 21:24
我也是新手,只不过之前遇到过,我是之前算超极化率,scf不收敛,三个都加上了没用,最后用了单点的波函 ...

多谢啦 也是个思路

3

帖子

0

威望

131

eV
积分
134

Level 2 能力者

15#
发表于 Post on 2022-11-17 10:39:53 | 只看该作者 Only view this author
星纹c 发表于 2022-11-9 21:24
我也是新手,只不过之前遇到过,我是之前算超极化率,scf不收敛,三个都加上了没用,最后用了单点的波函 ...

想问下您,算polar不收敛时加了什么关键词还是换了基组和方法,还有出现Inaccurate quadrature in CalDSu这种错误过么这种是怎么解决的呢

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 23:09 , Processed in 0.195960 second(s), 29 queries , Gzip On.

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