计算化学公社

标题: 频率分析失败（ Inaccurate quadrature in CalDSu.） [打印本页]

作者
Author: lr8616 时间: 2022-11-9 20:14
标题: 频率分析失败（ Inaccurate quadrature in CalDSu.）
请教各位老师：

对优化收敛的结构进行频率分析，最后失败结束，原因是：
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) ，还请忽略这个失误）

(, 下载次数 Times of downloads: 4)

作者
Author: wzkchem5 时间: 2022-11-9 20:29
为什么一定要加guess=INDO？
既然结构优化正常结束了，而频率计算又是在SCF里报错的，那么可以尝试一下读结构优化的收敛的波函数作为初猜

作者
Author: 星纹c 时间: 2022-11-9 20:31
你咋输入了俩guess；我之前也是用了INDO没用，就用了单点的初猜，你可以都试试

作者
Author: lr8616 时间: 2022-11-9 20:43
本帖最后由 lr8616 于 2022-11-9 20:47 编辑

星纹c 发表于 2022-11-9 20:31
你咋输入了俩guess；我之前也是用了INDO没用，就用了单点的初猜，你可以都试试

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

作者
Author: lr8616 时间: 2022-11-9 20:44

wzkchem5 发表于 2022-11-9 20:29
为什么一定要加guess=INDO？
既然结构优化正常结束了，而频率计算又是在SCF里报错的，那么可以尝试一下读 ...

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

作者
Author: lr8616 时间: 2022-11-9 20:46

星纹c 发表于 2022-11-9 20:31
你咋输入了俩guess；我之前也是用了INDO没用，就用了单点的初猜，你可以都试试

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

作者
Author: wzkchem5 时间: 2022-11-9 20:59

lr8616 发表于 2022-11-9 13:46
不好意思，想请问是guess=read 和 guess=INDO只能用一个吗？我开始只有guess=read读取结构优化的初猜波函 ...

对，只能用一个。只要了解guess=read和guess=INDO的原理就知道为什么。
guess=INDO其实是把初猜变得比默认初猜更差了，即使偶尔能解决问题也是碰巧，不是一个特别可靠的解决方法。

作者
Author: lr8616 时间: 2022-11-9 21:05

wzkchem5 发表于 2022-11-9 20:59
对，只能用一个。只要了解guess=read和guess=INDO的原理就知道为什么。
guess=INDO其实是把初猜变得比默 ...

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

作者
Author: sobereva 时间: 2022-11-9 21:17
没有可补充的

(, 下载次数 Times of downloads: 3)

作者
Author: 星纹c 时间: 2022-11-9 21:24

lr8616 发表于 2022-11-9 20:46
不好意思，想请问是guess=read 和 guess=INDO只能用一个吗？我开始只有guess=read读取结构优化的初猜波函 ...

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

作者
Author: wzkchem5 时间: 2022-11-9 22:14

lr8616 发表于 2022-11-9 14:05
我本来确实不理解，只是跟着帖子加。。看到您的回复忽然才意识到这个问题。我用guess=read出现了这个错误 ...

还有一种可以尝试的方法是，先加int=noxctest收敛以后，再读这个波函数，做一个不加int=noxctest的计算。
可能SCF迭代路径上经过的某些中间的波函数有数值精度不够的问题，但是收敛的波函数恰好没有这个问题，所以可以先关掉检查收敛以后，再打开检查。

作者
Author: lr8616 时间: 2022-11-10 09:19

sobereva 发表于 2022-11-9 21:17
没有可补充的

谢谢卢老师，现在能算下去了，感谢！

作者
Author: lr8616 时间: 2022-11-10 09:20

wzkchem5 发表于 2022-11-9 22:14
还有一种可以尝试的方法是，先加int=noxctest收敛以后，再读这个波函数，做一个不加int=noxctest的计算。 ...

谢谢您！我试试看

作者
Author: lr8616 时间: 2022-11-10 09:20

星纹c 发表于 2022-11-9 21:24
我也是新手，只不过之前遇到过，我是之前算超极化率，scf不收敛，三个都加上了没用，最后用了单点的波函 ...

多谢啦也是个思路

作者
Author: 刘文博 时间: 2022-11-17 10:39

星纹c 发表于 2022-11-9 21:24
我也是新手，只不过之前遇到过，我是之前算超极化率，scf不收敛，三个都加上了没用，最后用了单点的波函 ...

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

作者
Author: 星纹c 时间: 2022-11-17 12:59
本帖最后由星纹c 于 2022-11-17 13:01 编辑

刘文博发表于 2022-11-17 10:39
想问下您，算polar不收敛时加了什么关键词还是换了基组和方法，还有出现Inaccurate quadrature in CalDSu ...

我没换基组和方法。我关键词是CAM-B3LYP/Aug-cc-pvdz，出了Inaccurate quadrature in CalDSu，就是积分精度不够的意思。这个帖子前面社长发的图里面的前两条我加上了，第三条没用，但把用B3LYP/def2tzvp下做单点的波函数做初猜，总算是收敛了。
不收敛的话，你可以去看看社长发的一篇关于SCF收敛的文章。我就用了前几条就收敛了。

作者
Author: 刘文博 时间: 2022-11-17 14:35
你好，我是用CAM-B3LYP，BHandHYLP和M062X分别计算了我的分子的第一超极化率（polar），结构中含有过渡金属Zn就用了赝势基组，然后在计算Zn离子和原子的时候分别用SDD和LAN2DZ计算一遍，总是会有出现Inaccurate quadrature in CalDSu这种情况，要么就是用SDD算的时候离子状态下出现这种情况，要么就是用LAN2DZ时Zn离子的出现这种情况请问用卢老师图中的那几条会对结果产生影响么，请问这两种赝势不是一样的么，为什么会出现这种情况，新手小白

作者
Author: wzkchem5 时间: 2022-11-17 15:04

刘文博发表于 2022-11-17 07:35
你好，我是用CAM-B3LYP，BHandHYLP和M062X分别计算了我的分子的第一超极化率（polar），结构中含有过渡金属 ...

当然不一样，仔细阅读http://sobereva.com/373

作者
Author: 刘文博 时间: 2022-11-17 15:12
好的谢谢老师

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